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

68 阅读 0 评论 45 点赞

我是靠谱客的博主俭朴纸飞机，最近开发中收集的这篇文章主要介绍AOA估计中的MUSIC算法（Matlab代码实现），觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

????‍????个人主页：研学社的博客

????????????????欢迎来到本博客❤️❤️????????

????博主优势：????????????博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

????????????本文目录如下：????????????

目录

????1 概述

????2 运行结果

????3 参考文献

????4 Matlab代码实现

????1 概述

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

????2 运行结果

部分代码：%%
% 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