clear;?clc;?close?all;
%%?Parameter?defination
j?=?sqrt（-1）;
C?=?3e8;
Fc?=?300e6;???????????%載波頻率為100MHz
lamda?=?C/Fc;?????????%波長為1
Nsensor?=?8;?????????%陣元數目
Nsample?=?200;??????%信號采樣長度
Nsnapshot?=?Nsample;%快拍數
Nsource?=?2;??????????%假設信源數已知
d?=?lamda/2;??????????%陣元間距為半波長

source_doa=[10?60];%兩個信號的入射角度相差10度
Rxx?=?zeros（NsensorNsnapshot）;
w?=?[0.10.1];????????????????????%來波信號的頻率?
SNR?=[10?10];
for?i=1:Nsnapshot;?????%信號的幅度隨機產生，保證信號之間是不相關的
????signal0=sqrt（10.^（SNR（1）/10））*randn（1）*exp（j*w（1）*（i-1））;??%信噪比是15dB
????signal1=sqrt（10.^（SNR（2）/10））*randn（1）*exp（j*w（2）*（i-1））;??%信噪比是30dB
????noise=1/sqrt（2）;???????????????????????%噪聲功率為0dB
????x（:i）=signal0*exp（j*2*pi*1/2*sin（source_doa（1）*pi/180）*[0:Nsensor-1]‘）+...
???????????signal1*exp（j*2*pi*1/2*sin（source_doa（2）*pi/180）*[0:Nsensor-1]‘）+...
???????????noise*（randn（Nsensor1）+j*randn（Nsensor1））;
??%表示接收信號，且噪聲的幅度和相位隨機產生，保證噪聲與信號之間是不相關的
end
%%------用L1-SVD準則恢復原始信號?xx為恢復得到的信號-------
%%稀疏基矩陣為：角度估計精度為0.1度
Searching_doa?=?-90:1:90;
theta_len?=?length（Searching_doa）;?????????????????%角度參量搜索長度
A?=?zeros（Nsensortheta_len）;??%初始化稀疏基矩陣
for?m