clear?all;clc;
%=========參數(shù)初始化=========%
m=8;%陣元數(shù)
f=10e5;%載波頻率
theta1=30;%信號1方向
theta2=35;%信號2方向
th=[theta1;theta2];
nn=1024;%每個陣元采樣點數(shù)
Ts=[0:1/nn:1-1/nn];
s1=exp（-1i*2*pi*f*Ts）;%信號1
s2=s1;%信號2
S=[s1;s2];
degrad=pi/180;
n=0;
for?SN=-10:5:30
%信號1的信噪比
%SN=20;%信號2的信噪比
SN=[500
????0SN];

%--------------------------------------------

%生成陣元接收信號及噪聲
%--------------------------------------------
S2=（10.^（SN/10））*S;
wsum=0;
for?k=1:50
nr=randn（mnn）;%noise
ni=randn（mnn）;
U=sqrt（2）*（nr+1i*ni）;?%每路噪聲能量為1.
or=[];
%--------------------------------------------

%協(xié)方差矩陣及前后向平滑
%--------------------------------------------
tmp2=[0:m-1]‘;?%陣元位置
A=exp（-1i*pi*tmp2*sin（th‘*degrad））;?%方向矩陣
X=A*S2+U;
X1=X（1:7:）;X2=X（2:8:）;
Rxx1=X1*X1‘/nn;%第一個子陣
Rxx2=X2*X2‘/nn;%第二個子陣
RXX=（Rxx1+Rxx2）/2;
%--------------------------------------------

%奇異值分解
%--------------------------------------------
[QSSW]=svd（RX