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

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