水声方位估计 - MUSIC算法

算法基本思想

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

基于 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  

结果