MUSIC 算法MATLAB仿真源代码
clc
clear all
format long %将数据显示为长整型科学计数
N=200;%快拍数
doa=[20 60]/180*pi; %信号到达角
w=[pi/4 pi/3]';%信号频率
M=10;%阵元数
P=length(w); %信号个数
lambda=150;%波长
d=lambda/2;%阵元间距
snr=20;%信噪比
B=zeros(P,M); %创建一个P行M列的0矩阵
for k=1:P
B(k,:)=exp(-j*2*pi*d*sin(doa(k))/lambda*[0:M-1]); %矩阵赋值
end
B=B';
xx=2*exp(j*(w*[1:N])); %仿真信号
x=B*xx;　
x=x+awgn(x,snr);%加入高斯白噪声
R=x*x'; %数据协方差矩阵
[U,V]=eig(R); %求R的特征值和特征向量
UU=U(:,1:M-P); %估计噪声子空间
theta=-90:0.5:90;
%%谱峰搜索
for ii=1:length(theta)
AA=zeros(1,length(M));
for jj=0:M-1
　 AA(1+jj)=exp(-j*2*jj*pi*d*sin(theta(ii)/180*pi)/lambda);
end
WW=AA*UU*UU'*AA';
Pmusic(ii)=abs(1/ WW);
end
Pmusic=10*log10(Pmusic/max(Pmusic)); %空间谱函数
plot(theta,Pmusic,'-k')
xlabel('角度 \theta/degree')
ylabel('谱函数P(\theta) /dB')
title('MUSIC算法的DOA估计谱')
grid on(MUSIC algorithm MATLAB simulation source code)