经典的music算法的程序仿真

%%%%%%可能是程序写的不好，算法的仿真结果不行，分辨力较差@@@@@@@@@
clc;
clear ;
close all;
tic
M=20;%阵元数目
N=3;%信源数目
b=5;%形成的波束个数
snap=1000;%快拍数目
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;