经典MUSIC算法MATLAB仿真，带详细注释

经典MUSIC算法MATLAB仿真，带详细注释，仿真通过。可以直接使用。MATLAB阵列信号处理，DOA估计，MATLAB信号处理，信息与信号处理 专业，通信专业等学科学习的例程，经典仿真案例，学习入门必备，程序正确，可运行。

% 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 = 8;               % 阵列数量
dd = 0.5;               % space 
d=0:dd:（kelm-1）*dd;     % 
iwave = 3;              % number of DOA
theta = [10 30 60];     % 角度
snr = 10;               % input SNR （dB）
n = 500;                 % 
A=exp（-j*twpi*d.‘*sin（theta*derad））;%%%% direction matrix
S=randn（iwaven）;
X=A*S;
X1=awgn（Xsnr‘measured‘）;
Rxx=X1*X1‘/n;
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=10*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  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1122  2010-08-16 09:41  张小飞-MUSIC算法MATLAB程序\doa_music.m

     目录          0  2017-11-02 16:19  张小飞-MUSIC算法MATLAB程序

----------- ---------  ---------- -----  ----

                 1122                    2