music.m

ULA1一个天线阵的方向图,可以用来估计信号的来波方向,face_direction_pic面阵的方向图 music,MN-music,capon,ESPRINT为四种进测量来波方向的算法 该程序包主

clear;close all;clc;
%=========参数初始化=========%
num=8;%阵元数
f=10e9;%载波频率
f1=20;%信号1载波频率
f2=10;%信号2载波频率
a1=20;%信号1方向
a2=30;%信号2方向
c=3e8;%光速
lmda=c/f;%波长
d=lmda/2;%阵元间距
thita=0:1:180;
SNR1=20;%信噪比
SNR2=10;%信噪比
fs=500;
Ts=[0:1/fs:1];
snapshot=length(Ts);
s1=10^(SNR1/20)*sin(2*pi*f1*Ts);%信号1包络
s2=10^(SNR2/20)*sin(2*pi*f2*Ts);%信号2包络
athita1=exp(j*2*pi*d/lmda*cos(a1*pi/180)*[0:num-1]).';%信号导向矢量
athita2=exp(j*2*pi*d/lmda*cos(a2*pi/180)*[0:num-1]).';%信号导向矢量
x=athita1*s1+athita2*s2+sqrt(0.5)*(randn(num,snapshot)+j*randn(num,snapshot));%输入
%==============MUSIC===============%
R=x*x'/snapshot;
[V,D]=eig(R);
P=V(:,1:6)*V(:,1:6)';
for i=1:length(thita)
    A=exp(j*2*pi*d/lmda*cos(thita(i)*pi/180)*[0:num-1]).';
    S(i)=1/(A'*P*P'*A);   %P=P',P*P'=P*P=P
    %S(i)=1/(A'*P*A);%A信号空间,P*A噪声空间，
end
figure(1)
plot(thita,10*log10(abs(S)));
grid on;hold on;
title('MUSIC测角')
xlabel('方位角（度）')
ylabel('输出功率（dB）')