%%%%%%可能是程序?qū)懙牟缓茫惴ǖ姆抡娼Y(jié)果不行，分辨力較差@@@@@@@@@
clc;
clear?;
close?all;
tic
M=20;%陣元數(shù)目
N=3;%信源數(shù)目
b=5;%形成的波束個(gè)數(shù)
snap=1000;%快拍數(shù)目
C=3e8;
f0=10e6;
lamda=C/f0;
d=0.5*lamda;
%?k=d/lamda;
theta0=50;
theta1=30;

fs=1000;
ts=1/fs;
t=（0:snap-1）*ts;
a=[0:M-1]‘;%陣列矢量
u0=5;
u1=10;


s0=exp（j*2*pi*（f0*t+0.5*u0*t.^2））;
s1=exp（j*2*pi*（f0*t+0.5*u1*t.^2））;


%陣列流行矢量

a_theta0=exp（j*2*pi*d/lamda*a*sin（theta0/180*pi））;
a_theta1=exp（j*2*pi*d/lamda*a*sin（theta1/180*pi））;


A=[a_theta0?a_theta1];
S=[s0;s1];
X0=A*S;
SNR=10;
randn（‘state‘0）;
real_noise=randn（size（X0））;
randn（‘state‘3）;
imag_noise=randn（size（X0））;
noise0=（real_noise+j*imag_noise）/2^0.5;
noise=10^（-SNR/20）*noise0;
X=X0+noise;
Rx=X*X‘/length（t）;

jj=1;