capon.m

capon的程序

clear all 
d=1;                                 %天线阵元间距 
lma=2;                               %信号中心波长 
q1=1*pi/4; 
q2=pi/3;q3=pi/6;q4=3*pi/4;           %四个输入信号的方向 
A1=[exp(-2*pi*i*d*[0:6]*cos(q1)/lma)]';      %求阵因子 
A2=[exp(-2*pi*i*d*[0:6]*cos(q2)/lma)]'; 
A3=[exp(-2*pi*i*d*[0:6]*cos(q3)/lma)]'; 
A4=[exp(-2*pi*i*d*[0:6]*cos(q4)/lma)]'; 
A=[A1,A2,A3,A4];                     %得出A矩阵 
n=1:2000; 
v1=.015;                             %四个信号的频率 
v2=.05; 
v3=.02; 
v4=.035; 
d=[1.3*cos(v1*n);1*sin(v2*n);1*sin(v3*n);1*sin(v4*n)];  %输入信号矢量 
U=A*d;                               %总的输入信号 
U1=(U)'; 
c=cov(U*U1);                         %总输入信号的协方差矩阵 
[s,h]=eig(c);                        %求协方差矩阵的特征矢量及特征值 
Vn=s(:,[5:7]);                       %取出与零特征值对应的特征矢量 
ci=inv(c);                           %求协方差矩阵的逆矩阵 
q1b=[pi/180:pi/180:pi]; 
for n=1:length(q1b) 
    q1a(n)=q1b(n); 
    Ala=[exp(-2*pi*i*1*[0:6]*cos(q1a(n))/lma)]'; 
Pmusic(n)=(Ala)'*Ala*(inv((Ala)'*Vn*(Vn)'*Ala));  %应用MUSIC法估计输出 
Pcap(n)=inv((Ala)'*ci*(Ala)); 
T(n)=qla(n); 
P1=abs(Pmusic); 
P2=Abs(Pcap); 
end 
figure(1);                            %绘出应用MUSIC法估计的波达方向图 
polar(T,P1); 
figure(2)                           %绘出应用MUSIC法估计的波达方向图 
T1=T*180/pi 
semilogy(T1,P1);grid 
figure(3)                            %绘出应用Capon法估计的波达方向图 
polar(T,P2) 
figure(4)                            %绘出应用Capon法估计的波达方向图 
T1=T*180/pi; 
semilogy(T1,P2);grid