clear?all;
close?all;
derad=pi/180;
radeg=180/pi;
twpi=2*pi;
kelm=8;?%陣元數(shù)
dd=0.5;%陣元間距
d=0:dd:（kelm-1）*dd;
iwave=3;%信號源數(shù)
theta=[30?60?120];?%波達方向
snr=10;

A=exp（j*twpi*d.‘*cos（theta*derad））;?%方向矢量

n=[10?50?100?200?400];;%采樣數(shù)
l_n=length（n）;

SP_y=zeros（180l_n）;

for?q=1:l_n
????S=randn（iwaven（q））;
X=A*S;%接收信號
X1=awgn（Xsnr‘measured‘）;%添加噪聲
Rxx=X1*X1‘/n（q）;%計算協(xié)方差
InvS=inv（Rxx）;
[EVD]=eig（Rxx）;?%特征分解
EVA=diag（D）‘;
[EVAI]=sort（EVA）;?%特征值降序排序
EVA=fliplr（EVA）;?%左右翻轉(zhuǎn)?從大到小順序
EV=fliplr（EV（:I））;?%?對應(yīng)特征向量排序

%構(gòu)造MUSIC譜函數(shù)
for?iang=1:180
????%angle（iang）=（iang-181）/2;
????phim=derad*iang;
?