Music算法实现均匀平面阵、直线阵圆阵附matlab代码

% 均匀平面阵Music算法
clear all;clc;
close all
mx=5;my=4;%x轴和y轴阵元个数; 
sn=4;%信号个数
dw=0.2;%半径波长比   
snr1=[50,50,50,50];
snr=10.^(snr1/20);%信噪比（幅度表示）
If=[20 30 30 50];%信号中频率(单位：MHz);
f=[1.5,4.0,2.0,5.0];%信号调制频率(单位：MHz);
fs=120;%采样频率(满足采样定理)(单位：MHz)；
N=4096;n=(1:N);%采样点数；
fangwei=[10,25,135,170];%信号方位角
yangjiao=[60 80 20 10]
Ss=zeros(sn,N);
  for i=1:sn
      for m=1:mx
          for mm=1:my
            daoxiang((m-1)*my mm,i)=exp(-j*2*pi*dw*(cos(2*pi*(m-1)/mx-fangwei(i)*pi/180) sin(2*pi*(m-1)/mx-fangwei(i)*pi/180))*cos(yangjiao(i)*pi/180));%导向矢量
          end
      end
      ss(i,:)=snr(i)*(1 0.3*sin(2*pi*f(i)*n/fs)).*exp(j*2*pi*n*If(i)/fs);%AM调制信号（S（t））
  end
 Signal=daoxiang*ss;
 noise=randn(mx*my,N);
 noise_h=(hilbert(noise.')).'/sqrt(2);%对噪声进行希尔伯特变化映射到复数空间
 x=Signal noise_h;%接收信号（y（t））
 R=x*x'/N;
 [tzxiangliang,tzzhi]=eig(R);
 Nspace=tzxiangliang(:,1:mx*my-sn);%噪声子空间对应小的特征值（从小到大排列）
 for azi=1:1:180
     for ele=1:1:90
       for m=1:mx
           for mm=1:my
           AQ1((m-1)*my mm,1)=exp(-j*2*pi*dw*(cos(2*pi*(m-1)/mx-azi*pi/180) sin(2*pi*(m-1)/mx-azi*pi/180))*cos(ele*pi/180));;
           end
       end
     
     Power=AQ1'*Nspace*Nspace'*AQ1; %在1-180度范围内进行计算
     P(ele,azi)=-10*log10(abs(Power));
     end
 end
 figure;
 mesh(P);
 title('5*4均匀平面阵;信噪比：[50,50,50,50]；距离波长比：0.5');
 xlabel('方位角');ylabel('仰角');
 zlabel('空间谱/db');
 grid;