水声方位估计 - MUSIC算法

最新推荐文章于 2024-03-25 02:44:34 发布

温逗死

最新推荐文章于 2024-03-25 02:44:34 发布

阅读量893

点赞数 1

文章标签：算法矢量水听器 DOA MUSIC MATLAB

本文链接：https://blog.csdn.net/qq_43142509/article/details/121176014

版权

算法基本思想

算法的基本思想是对任意阵列输出数据的协方差矩阵进行特征值分解，从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间，然后利用这两个子空间的正交性构造空间谱函数，通过峰值搜素，估测信号的角度。

基于 MUSIC 算法的 MATLAB 仿真

% DOA estimation by MUSIC 
clear all
close all
derad = pi/180;             % deg -> rad
radeg = 180/pi;
twpi = 2*pi;
kelm = 8;                   % 阵列数量
dd = 0.5;                   % space 
d=0:dd:(kelm-1)*dd;         % 
iwave = 3;                  % number of DOA
theta = [10 30 90];         % 角度
snr = 10;                   % input SNR (dB)
n = 500;                    % 
A=exp(-j*twpi*d.'*sin(theta*derad));%%%% direction matrix
S=randn(iwave,n);
X=A*S;
X1=awgn(X,snr,'measured');
Rxx=X1*X1'/n;
InvS=inv(Rxx); 
[EV,D]=eig(Rxx);
EVA=diag(D)';
[EVA,I]=sort(EVA);
EVA=fliplr(EVA);
EV=fliplr(EV(:,I));

% MUSIC
for iang = 1:361
        angle(iang)=(iang-181)/2;
        phim=derad*angle(iang);
        a=exp(-j*twpi*d*sin(phim)).';
        L=iwave;    
        En=EV(:,L+1:kelm);
        SP(iang)=(a'*a)/(a'*En*En'*a);
end
   
figure()
SP=abs(SP);
SPmax=max(SP);
SP=10*log10(SP/SPmax);
h=plot(angle,SP);
set(h,'Linewidth',2)
xlabel('angle (degree)')
ylabel('magnitude (dB)')
axis([-90 90 -60 0])
set(gca, 'XTick',[-90:30:90])
grid on