MUSIC算法DOA估计 MATLAB程序

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（