tu32.m

9单元线阵MUSIC算法 9单元线阵MUSIC算法

                           %九单元线阵DOA估计 单信号加噪时Capon法%
 %信号处理部分                          
d=1; %天线阵元的间距
a=2;  %信号中心波长
b=20;%信号幅度
q1=30*(pi/180);%输入信号的方向
A1=[exp(-2*pi*j*d*[0:8]*cos(q1)/a)]';%求阵因子
A=[A1];%得出A 矩阵
n=1:2000;
v1=.015;%信号的频率

d=[b*cos(v1*n)];%输入信号矢量
randn('state',1);
en=randn(9,2000);%高斯白噪声
U=A*d;%总的输入信号


%算法部分
U1=(U)';
c=cov(U*U1);%总输入信号的协方差矩阵
ci=inv(c);%求协方差矩阵的逆矩阵
q1b=[pi/180:pi/180:pi];
for n2=1:length(q1b)
    q1a(n2)=q1b(n2);
    A1a=[exp(-2*pi*j*1*[0:8]*cos(q1a(n2))/a)]';
    Pcap(n2)=inv((A1a)'*ci*(A1a));%应用Capon 法估计输出
    T(n2)=q1a(n2);
    P2=abs(Pcap);
end

%绘出应用 法估Capon 计的波达方向极坐标图
subplot(1,2,1)
polar(T,P2)
xlabel('Capon 法估计的波达方向极坐标图')
%绘出应用Capon 法估计的波达方向图
subplot(1,2,2)
T1=T*180/pi;
semilogy(T1,P2);grid
xlabel('Capon 法估计的波达方向图')