% 均匀平面阵Music算法
clear all;clc;
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
daoxiang1(m,i)=exp(-j*2*pi*dw*(m-1)*cos(fangwei(i)*pi/180)*cos(yangjiao(i)*pi/180));
end
for mm=1:my
daoxiang2(mm,i)=exp(-j*2*pi*dw*mm*sin(fangwei(i)*pi/180)*cos(yangjiao(i)*pi/180));
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
daoxiang=[daoxiang1;daoxiang2];
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
daoxiang3(m,1)=exp(-j*2*pi*dw*(m-1)*cos(azi*pi/180)*cos(ele*pi/180));
end
for mm=1:my
daoxiang4(mm,1)=exp(-j*2*pi*dw*mm*sin(azi*pi/180)*cos(ele*pi/180));
end
AQ1=[daoxiang3;daoxiang4];
Power=AQ1'*Nspace*Nspace'*AQ1; %在1-180度范围内进行计算
P(ele,azi)=-10*log10(abs(Power));
end
end
figure;
mesh(P);
title('九元L阵;信噪比:[50,50,50,50];距离波长比:0.5');
xlabel('方位角');ylabel('仰角');
zlabel('空间谱/db');
grid;
评论0