%?寬帶信號頻域處理?

?
clear?all?
clc?
%?f0?=?1.3*10^9;??????????????????%?信號中心頻率?
%?fr=1.3*10^9;??
%?B?=?50*10^6;????????????????????%?信號帶寬?
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）;?????????????????%?信號持續時間，要求大于孔徑渡越時間?
%?Tr?=?5;????????????????%?工作周期?
%?T?=?1;?????????????????%?信號持續時間，要求大于孔徑渡越時間?
fs?=?100*10^6;??????????????????%?采樣頻率?
snr?=?[40;0];???????????????????%?信噪比?
w0?=?40/180*pi;??????????????????%?指向角?
M?=8;??????????????????????????%?陣元數為M?
Kr=1;?
Ks?=?2;?????????????????????????%?信號數目?
seta?=[0/180*pi40/180*pi];?????%?干擾信號方向?
Nm?=?6;????????????????????????%?駐留周期個數?
Nr?=?Tr*fs;?????????????????????%?采樣點數?
N?=?T*fs;?
J?=?1000;???????????????????????%?FFT的點數?
NN?=?Nm*Nr/J;???????????????????%?頻域快拍數??為么？?
c=?340;?????????????????????%?光速?
d?=?0.5*c/f0;???????????????????%?陣元間距?
x?=?zeros（M1）;?????????????????%?數據矢量?
R?=?zeros（MM）;?????????????????%?接收數據協方差矩陣?
?
tic?
%?產生陣列接收的寬帶數據?
x2?=?LFMsource（Mf0flBfsTsetasnr2）;?%?期望信號的脈沖寬度為T，出現時間不定?
?
%?產生干擾信號?
x1?=?LFMsource（Mf0flBfsTrsetasnr1）;?
%?u?=?1:Nr;?
%?for?k=1:Kr?
%?for?n=1:Nr?
%?????fai（kn）?=?rand;?
%??s（kn）?=?exp（i*2*pi*（fr*u（n）/fs+fai（kn）））;???????%?源信號?
%?end?
%?end?
%?for?k=1:Kr?
%?for?m=1:M?
%?????????A（mk）=exp（-i*2*pi*d*fr*（m-1）*sin（seta（k））/c）;?
%?????end?
%?end?
%?%?陣列接收信號?
%?x1?=?10.^（snr（1）/20）.*A*s;??
%?陣列接收的和信號?
t1?=?1000;?????????????????????????????????%?任意選擇的期望信號出現位置?
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?
%?構造方向矢量?
w?=?-90*pi/180:0.01:90*pi/180;?
WW?=?length（w）;?
P?=?zeros（1WW）;?
Wav=zeros（M1）;?
%?計算不同子頻點的加權向量???????????????????????????????????????????????????????????????
F?=?fl:B/（J/2-1）:fh;?
for?k?=?1:250??????????????????????%為么是250？?
????for?m?=?1:M?????????????????????%?估計子頻點對應的協方差矩陣?
????????for?n?=?1:NN?
????????????xf（nm）?=?X（knm）;?
????????end?
????end?
????R?=?xf.‘*conj（xf）./NN;%???xf.‘???
????Ri?=?inv（R）;?
????for?m=1:M?
????????a0（m1）?=?exp（-i*2*pi*d*F（k+250）*（m-1）*sin（w0）/c）;%?指向向量??
????end??
????W?=?Ri*a0/（a0‘*Ri*a0）;????