clc
clear?all
format?long?
%雙向平滑算法
%FBSS適用于最一般的情況
%%
L???????????????????????=100;?????????????????????%快拍數(shù)
M???????????????????????=9;??????????????????????%陣元數(shù)
j=sqrt（-1）;
doa?????????????????????=[-20?0?40]/180*pi;????%入射角度
doa1????????????????????=[-20?0?40];???????
w???????????????????????=[pi/2pi/2pi/3]‘;??%信號頻率
N???????????????????????=length（doa）;?????????????%信源數(shù)
coefficient?????????????=[1-0.64-j*0.650.75-j*0.35];%衰減系數(shù)
lamda???????????????????=150;?????????????????????%波長
d???????????????????????=lamda/2;?????????????????%間距
snr?????????????????????=10;???????????????????????%信噪比
noise_power?????????????=2;???????????????????????%噪聲功率
m???????????????????????=6;???????????????????????%每個子陣陣元數(shù)量
p???????????????????????=M-m+1;???????????????????%子陣數(shù)量
%%
B=zeros（NM）;
for?k=1:N
????B（k:）=exp（-j*2*pi*d*sin（doa（k））/lamda*[0:M-1]）;
end
B=B.‘;
w_noise?=?randn（ML）+j*randn（ML）;
Signal_Amplitude?=?sqrt（2*noise_power*10^（snr/10））;
s=Signal_Amplitude?*exp（j*（w*[1:L]））;
xx=B*diag（coefficient）*s;
X=xx+w_noise;
tic;
%%
%?截取子陣，作協(xié)方差矩陣再取平均
H=zeros（mm）;
for?k=1:p?
????X1=X（k:k+m-1:）;%直接對總的加噪信號模型進(jìn)行截取即可
????X2=X1*X1‘;%每個子陣的協(xié)方差矩陣
????H=H+X2;
end
H=H/p;%前向平滑矩陣
XX=flipud（X）