我模拟了一个相控阵天线。我还模拟了一个点源。我不知道如何在不使用源位置的情况下计算伪谱矢量。Capon 和 Bartlett 等算法使用源的位置来计算自相关矩阵,该矩阵进一步用于计算伪谱。如果首先需要源位置来运行算法,如何估计到达角?
这是代码:
M=6;
sig2=.1;
th1=-3*pi/180;
th2=3*pi/180;
a1=[1];
a2=[1];
一个=[1];
对于 i=2:M
a1=[a1 exp(-1j*i*pi*sin(th1))];
a2=[a2 exp(-1j*i*pi*sin(th2))];
结尾
A=[a1'a2'];
Rss=[1 0;0 1];
Rrr=A Rss A'+sig2*eye(6);
对于 k=1:180;
th(k)=-pi/6+pi*k/(3*180);
清除一个
一个=[1];
对于 jj=2:M
a = [a exp(-1j*jj*pi*sin(th(k)))];
结尾
P(k)=实数(1./(conj(a)*inv(Rrr)*a.'));
结尾
数字;
绘图(th*180/pi,10*log10(P/max(P)),'k')
网格上
xlabel('角度')
ylabel('|P(\theta)| (dB)')