%?頻域波束形成

clear?all?
clc?
f0?=?2.2*10^3;??????????????????%?信號中心頻率?
fr=2.2*10^3;??
B?=?3.6*10^3;?
fl?=?f0-B/2;????????????????????%?信號起始頻率?
fh?=?f0+B/2;?
Tr?=?50*10^（-6）;????????????????%?工作周期?
T?=?10*10^（-6）;?????????????????%?信號持續(xù)時間，要求大于孔徑渡越時間?
%?Tr?=?5;????????????????%?工作周期?
%?T?=?1;?????????????????%?信號持續(xù)時間，要求大于孔徑渡越時間?
fs?=?100*10^6;??????????????????%?采樣頻率?
snr?=?[40;0];???????????????????%?信噪比?
w0?=?40/180*pi;??????????????????%?指向角?
M?=8;??????????????????????????%?陣元數(shù)為M?
Kr=1;?
Ks?=?2;?????????????????????????%?信號數(shù)目?
seta?=[0/180*pi40/180*pi];?????%?干擾信號方向?
Nm?=?6;????????????????????????%?駐留周期個數(shù)?
Nr?=?Tr*fs;?????????????????????%?采樣點數(shù)?
N?=?T*fs;?
J?=?1000;???????????????????????%?FFT的點數(shù)?
NN?=?Nm*Nr/J;???????????????????%?頻域快拍數(shù)??為么？?
c=?340;?????????????????????%?光速?
d?=?0.5*c/f0;???????????????????%?陣元間距?
x?=?zeros（M1）;?????????????????%?數(shù)據(jù)矢量?
R?=?zeros（MM）;?????????????????%?接收數(shù)據(jù)協(xié)方差矩陣?
?
tic?
%?產(chǎn)生陣列接收的寬帶數(shù)據(jù)?
x2?=?LFMsource（Mf0flBfsTsetasnr2）;?%?期望信號的脈沖寬度為T，出現(xiàn)時間不定?
?
%?產(chǎn)生干擾信號?
x1?=?LFMsource（Mf0flBfsTrsetasnr1）;?
t1?=?1000;?????????????????????????????????%?任意選擇的期望信號出現(xiàn)位置?
x1（:（t1+1）:（t1+N））?=?x1（:（t1+1）:（t1+N））?+?x2;?
for?nm=1:Nm?
????x（:（（nm-1）*Nr+1）:（nm*Nr））?=?x1;?
end?
noise?=?randn（MNm*Nr）+i*randn（MNm*Nr）;???%?噪聲?
x?=?x?+?noise;?
?
%?時域信號變換為頻域信號?
fft_8_1;?
toc?
tic?
%?構(gòu)造方向矢量?
w?=?-90*pi/180:0.01:90*pi/180;?
WW?=?length（w）;?
P?=?zeros（1WW）;?
Wav=zeros（M1）;?
%?計算不同子頻點的加權(quán)向量???????????????????????????????????????????????????????????????
F?=?fl:B/（J/2-1）:fh;?
for?k?=?1:250??????????????????????%為么是250？?
????for?m?=?1:M