clear?all;
close?all;
derad=pi/180;
radeg=180/pi;
twpi=2*pi;
kelm=8;?%陣元數
dd=0.5;%陣元間距
d=0:dd:（kelm-1）*dd;
iwave=3;%信號源數
theta=[30?60?120];?%波達方向
snr=[-5051015];
n=200;%采樣數
A=exp（j*twpi*d.‘*cos（theta*derad））;?%方向矢量
S=randn（iwaven）;
X=A*S;%接收信號
l_snr=length（snr）;
SP_y=zeros（180l_snr）;
for?q=1:l_snr
X1=awgn（Xsnr（q）‘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:180
????%angle（iang）=（iang-181）/2;
????phim=derad*iang;
????a=exp