MUSIC算法程序

MUSIC算法 应用于DOA估计的超分辨算法 多重信号旋转子空间不变算法 波达方向估计最经典的算法

% DOA estimation by MUSIC 
% Developed by xiaofei zhang （南京航空航天大学 电子工程系 张小飞）
% EMAIL:zhangxiaofei@nuaa.edu.cn
clear all
close all
derad = pi/180;        % deg -> rad
radeg = 180/pi;
twpi = 2*pi;
kelm = 12;               % 阵列数量
dd = 0.5;               % space 
d=0:dd:（kelm-1）*dd;     % 
iwave = 1;              % number of DOA
theta =25;     % 角度
snr = -10;               % input SNR （dB）
n = 100;                 % 
A=exp（-j*twpi*d.‘*sin（theta*derad））;%%%% direction matrix
S=randn（iwaven）;
X=A*S;

X1=awgn（Xsnr‘measured‘）;
rank（X1）
Rxx=X1*X1‘/n;
rank（Rxx）
InvS=inv（Rxx）; %%%%
[EVD]=eig（Rxx）;%%%% 
EVA=diag（D）‘;
[EVAI]=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
   
% 
SP=abs（SP）;
SPmax=max（SP）;
SP=20*log10（SP/SPmax）;
h=plot（angleSP）;
set（h‘Linewidth‘2）
xlabel（‘angle （degree）‘）
ylabel（‘magnitude （dB）‘）
axis（[-90 90 -60 0]）
set（gca ‘XTick‘[-90:30:90]）
grid on  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-06-13 11:15  1.MUSIC算法MATLAB程序\
     文件        1135  2018-04-10 19:40  1.MUSIC算法MATLAB程序\doa_music.m
     文件        8653  2017-05-10 19:49  1.MUSIC算法MATLAB程序\RMSEofMUSICandSVM.m