clc
clear?all?
close?all
ima=sqrt（-1）;
element_num=20;???????%陣元數
d_lamda=0.5;?????????%陣元間距與波長的關系
theta=-90:0.5:90;??????%搜索范圍
theta0=0;????????????%來波方向
theta1=20;???????????%干擾方向
L=32;??????????????%采樣單元數
for?i=1:L;
????amp0=10*randn（1）;%信號的幅度隨機產生，保證信號之間是不相關的
????amp1=200*randn（1）;
ampn=1;
????s（:i）=amp0*exp（ima*2*pi*d_lamda*sin（theta0*pi/180）*[0:element_num-1]‘）;%來波信號
????j（:i）=amp1*exp（ima*2*pi*d_lamda*sin（theta1*pi/180）*[0:element_num-1]‘）;%干擾信號
????n（:i）=ampn*（randn（element_num1）+ima*randn（element_num1））;????????%噪聲信號
????x（:i）=?s（:i）+j（:i）+n（:i）;???????%接受信號
%表示接收信號，且噪聲的幅度和相位隨機產生，保證噪聲與信號之間是不相關的
end
Rs=1/L*s*?s‘;
Rnj=1/L*（j*j‘+n*n‘）;??
[VD]=eig（RsRnj）;????%廣義特征值分解，V是特征向量，D是特征值
[DI]=sort（diag（D））;??%取出D的對角元素即所有特征值，對其排序I存放序號
Wopt=V（:I（element_num））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-04-16?16:10??MSNR\
?????文件????????1989??2017-04-16?16:10??MSNR\MSNR.m