空间平滑music算法

N=2;%信源数目
C=3e8;
f0=10e6;
lamda=C/f0;
d_X=0.5;
d=d_X*lamda;
fs=1000;
ts=1/fs;
M_zhenyuan=7;
kelm = 6;   %分成两个子阵，每个子阵6个阵元
SNAP_kuaipai=200;
t=（0:SNAP_kuaipai-1）*ts;
a=[0:M_zhenyuan-1]‘;%阵列矢量
b=[0:kelm-1]‘;
thea_0=10;
thea_1=40;
SNR=10;
j=sqrt（-1）;
s0=exp（j*2*pi*f0*t）;
s1=exp（j*2*pi*f0*t）;
%阵列流行矢量
a_theta0=exp（j*2*pi*d/lamda*a*sin（thea_0/180*pi））;
a_theta1=exp（j*2*pi*d/lamda*a*sin（thea_1/180*pi））;
A=[a_theta0 a_theta1];
S=[s0;s1];
X0=A*S;
X=awgn（X0SNR）;
Rxxm=X*X‘/length（t）;

[MMM]=size（Rxxm）;
N_zizhen=M-ke