clc;
clear?all;
%?---------------------------參數初始化------------------------------------
M=32;?????????????????????????%?陣元數目
L=1024;??????????????????????%?信號長度/快拍數?
lam=0.15;????????????????????%?信號波長
d=1/2*lam;???????????????????%?陣元間距
angle=[-pi/3?pi/6?35*pi/180];??????%?信號入射角度
snr=-10;?????????????????????????????%信噪比
N=3;?????????????????????????%信源數
%?---------------------------信號模型--------------------------------------
A=zeros（NM）;
for?k=1:N
A（k:）=exp（-j*2*pi/lam*（[0:M-1]*d）*sin（angle（k）））;???
end
A=A‘;?????????????????????????????????????????????%?導向向量
SS=2*（randn（1L）+j*randn（1L））;
Wn=[0.36?0.38];
[ba]?=?fir1（512Wn‘bandpass‘）;??????%?設計帶通濾波器
SS1=filter（baSS）;?
Wn=[0.38?0.4];
[ba]?=?fir1（512Wn‘bandpass‘）;??????%?設計帶通濾波器
SS2=filter（baSS）;?
Wn=[0.4?0.42];
[ba]?=?fir1（512Wn‘bandpass‘）;??????%?設計帶通濾波器
SS3=filter（baSS）;?
SS=zeros（NL）;???SS=[SS1;SS2;SS3];
RS=SS*SS‘/L;?
gls=trace（RS）/N;
R=A*RS*A‘+（gls*10^（-snr/10））*eye（M）;?????????????%白噪聲協方差矩陣為對角陣
%?---------------------------角度估計--------------------------------------
theta=-pi/2:pi/3600:pi/2;???????%?ULA估計角度變化的范圍和頻率選擇?
for?k2=1:length（theta）??????%?角度估計
????AA=exp（-j*2*pi/lam*（[0:M-1]*d）*si