music.txt

《空间谱估计原理与算法》多重信号分类法估计matlab源代码

clear;
q1=pi*2/3;q2=pi*1/3;q3=pi*3/2;q4=pi*5.8/6;                  %四个输入信号的方向
E=1;                                                        
lam=1;                                                      %信号波长
r=1;                                                        %天线阵的半径
m1=1;m2=1;m3=1;m4=1;m5=1;m6=1;m7=1;m8=1;                    %信号幅度
a=0:7;%a的取值
for n=1:length(a)
    A1(n)=[exp(-j*pi*cos(q1-(pi*a(n)/4))/lam)];%第n个阵元对第1个信号的影响
    A2(n)=[exp(-j*pi*cos(q2-(pi*a(n)/4))/lam)];%第n个阵元对第2个信号的影响
    A3(n)=[exp(-j*pi*cos(q3-(pi*a(n)/4))/lam)];%第n个阵元对第3个信号的影响
    A4(n)=[exp(-j*pi*cos(q4-(pi*a(n)/4))/lam)];%第n个阵元对第4个信号的影响
end
A=[A1',A2',A3',A4'];%得出A矩阵——8x4                                         
n=1:1900;%n的取值（抽样）
v1=.06;%v1-v4代表四个信号的频率
v2=.02;
v3=.03;
v4=.073;
D=[1*cos(v1*n);1*sin(v2*n);1*sin(v3*n);1*square(v4*n)];      %四个输入信号
U=A*D;   %8x1900                                             %总的输入信号
U1=(U)'; %1900x8
c=cov(U*U1);   %8阶方阵                                      %总输入信号的协方差矩阵
[s,z]=eig(c)                                                %求协方差矩阵的特征矢量及特征值
Vn=s(:,[1:4]);                                               %取出与零特征值对应的特征矢量
ci=inv(c); %求逆矩阵，得到8阶方阵
%bb=[1 0 0 0]';
%Wopte=A'\bb;                                                 %求解线性方程组（求wopte)
q1b=[2*pi:-2*pi/180:2*pi/180];  %0-2pi内均匀分割——1x180
b=0:7; %b的取值——1x8
for t=1:length(q1b)%开始第一重循环
    h(t)=q1b(t);%0-2*pi内均匀取值
    for m=1:length(b)
        Aa(m)=[exp(-j*pi*cos(h(t)-(pi*a(m)/4))/lam)];%1x8矩阵——阵列流型
    end
    A1a=Aa';%求转置矩阵
    Pmusic(t)=(A1a)'*A1a*(inv((A1a)'*Vn*(Vn)'*A1a));
    Pcap(t)=inv((A1a)'*ci*(A1a));
    T(t)=q1b(t); %将q1b取值反应到T坐标轴上                     
    P1=abs(Pmusic);                                           %MUSIC算法估计 
    P2=abs(Pcap);                                             %CAPON算法估计
    %Ye(t)=Wopte'*A1a;
    %p=abs(Ye);
end
figure name 'MUSIC'
T1=T*180/pi;
semilogy(T1,P1);grid                                          %MUSIC算法波达方向估计
figure name 'Polar plans-MUSIC'%MUSIC极坐标图
polar(T1,P1)
figure name 'CAPON'%CAPON极坐标图
T1=T*180/pi;
semilogy(T1,P2);grid                                          %CAPON算法波达方向估计
figure name 'Polar plans-CAPON'
polar(T1,P2)
%figure(3)
%polar(T,p)                                                    %绘出应用矩阵运算求解加权系数后的波束