clear?all;
close?all;
clc;
%%設(shè)置信號為窄帶信號，幅度為1
S=[111]
imag1=sqrt（-1）;
element_num=8;%陣元數(shù)為8
d_lamda=1/2;%陣元間距d與波長lamda的關(guān)系
theta=linspace（-pi/2pi/22000）;
theta0=0/180*pi;%來波方向?（我覺得應(yīng)該是天線陣的指向）
%theta0=0;%來波方向
theta1=12/180*pi;
theta2=30/180*pi;
L=20;
kk_sum=0;
SNR=8;
SNR1=16;
SNR2=24;
w_l0=exp（-imag1*2*pi*d_lamda*sin（theta0）*[0:element_num-1]‘）;
w_l1=exp（-imag1*2*pi*d_lamda*sin（theta1）*[0:element_num-1]‘）;
w_l2=exp（-imag1*2*pi*d_lamda*sin（theta2）*[0:element_num-1]‘）;
A=[w_l0w_l1w_l2];
S=randn（length（S）L）
x=A*S;
X=awgn（xSNR‘measured‘）;
X1=awgn（xSNR1‘measured‘）;
X2=awgn（xSNR2‘measured‘）;
R=X*X‘/L;
R1=X1*X1‘/L;
R2=X2*X2‘/L;
[VD]=eig（R）;
[V1D1]=eig（R1）;
[V2D2]=eig（R2）;

Sp=diag（D）‘;
Sp1=diag（D1）‘;
Sp2=diag（D2）‘;
[Sp_randI]=sort（Sp‘descend‘）;
[Sp_rand1I1]=sort（Sp1‘descend‘）;
[Sp_rand2I2]=sort（Sp2‘d