AOA估计中的MUSIC算法（Matlab代码实现）

           目录

💥1 概述

📚2 运行结果

🎉3 参考文献

👨‍💻4 Matlab代码

---

💥1 概述

随着阵列信号处理技术的不断发展,到达角估计(Angle Of Arrival)的研究在移动通信系统中具有重要意义。通过分析经典MUSIC算法,针对其搜索空间较大,受噪声影响较大的因素,提出改进的根值MUSIC算法。对影响MUSIC算法性能的信噪比、采样数、阵元数、入射角度等因素以及根值MUSIC算法进行仿真,仿真结果表明改进的根值MUSIC算法是有效的。

📚2 运行结果

 

 

🎉3 参考文献

[1]韩卫杰. 改进MUSIC算法在AOA估计中的研究[D].西南交通大学,2006.

👨‍💻4 Matlab代码

主函数部分代码：%%
% Code name: 2D MUSIC algorithm 
clc
clear 
close all
format long
N=200;fs=2e11;
doa=[40 60]/180*pi;
w=[pi/4 pi/4]'*95e9;
M=10;
Msub=3;
P=length(w);
c=3e8;
lambda=c*2*pi/w(1);
deltad=lambda/2;
% deltad=lambda/1.5;
snr=10;
D=zeros(P,M);
for k=1:P
D(k,:)=exp(-1i*2*[0:M-1]*pi*deltad*sin(doa(k))/lambda);
end
s=2*exp(1i*(w*[1:N]));
x=D'*s;
x=x+awgn(x,snr);
figure,
%% Without spatial smoothing
R_old=x*x';
J=fliplr(eye(M));
R_old=R_old+J*conj(R_old)*J;
[N,~]=eig(R_old);
NN=N(:,1:M-P);
theta=-90:0.5:90;
for ii=1:length(theta)
SS=zeros(1,length(M));
for jj=0:M-1
SS(1+jj)=exp(-1i*2*jj*pi*deltad*sin(theta(ii)/180*pi)/lambda);
end
PP= SS*NN*NN'*SS';
Pmusic_im(ii)=abs(1/PP);
end
Pmusic_im=10*log10(Pmusic_im/max(Pmusic_im));
plot(theta,Pmusic_im,'c');
hold on
%% When we use only one sub sub-array
R_sub_ma=[]; 
for t=1:M-Msub+1
x_sub=x(t:t+Msub-1,:);
R_sub=x_sub*x_sub';
R_sub_ma(t,:,:)=R_sub;
end