bodafangxiang.m

用Capon法和MUSIC法对七单元线形天线阵在四信号输入情况下的编程仿真

d=1;
lma=2;
q1=1*pi/4;
q2=pi/3;q3=pi/6;q4=3*pi/4;
A1=[exp(-2*pi*j*d*[0:6]*cos(q1)/lma)]';
A2=[exp(-2*pi*j*d*[0:6]*cos(q2)/lma)]';
A3=[exp(-2*pi*j*d*[0:6]*cos(q3)/lma)]';
A4=[exp(-2*pi*j*d*[0:6]*cos(q4)/lma)]';
A=[A1,A2,A3,A4];
n=1:1900;
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);
A1a=[exp(-2*pi*j*1*[0:6]*cos(q1a(n))/lma)]';
Pmusic(n)=(A1a)'*A1a*(inv((A1a)'*Vn*(Vn)'*A1a));
Pcap(n)=inv((A1a)'*ci*(A1a));
   T(n)=q1a(n);
P1=abs(Pmusic);
P2=abs(Pcap);
end
figure(1)
polar(T,P1)
figure(2)
T1=T*180/pi;
semilogy(T1,P1);grid
figure(3)
polar(T,P2)
figure(4)
T1=T*180/pi;
semilogy(T1,P2);grid