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=[10?30?60];%入射角度
snr=6;%信噪比
n=500;%采樣數(shù)
A=exp（-j*twpi*d.‘*sin（theta*derad））;%方向矢量
S=randn（iwaven）;%信源信號
X=A*S;%接收信號
X1=awgn（Xsnr‘measured‘）;%添加噪聲
Rxx=X1*X1‘/n;%計算協(xié)方差矩陣
[EVD]=eig（Rxx）;%特征值分解
EVA=diag（D）‘;
[EVAI]=sort（EVA）;%特征值從小到大排序
EVA=fliplr（EVA）;%左右翻轉，從大到小排序
EV=fliplr（EV（:I））;%對應特征矢量