%?寬帶信號頻域處理用一致聚焦算法進行波束形成
%?陣元數為8，間距為半波長，信號中心頻率為3.1G，帶寬50M
%?干擾為窄帶頻率為3.1G
%RSS算法
%?2008.3.10



%%%添加注釋研究
clear?all;
clc;
f0?=?3.1*10^9;??????????????????%?信號中心頻率
fr?=?3.1*10^9;?????????????????%?干擾信號頻率
B?=?300*10^6;????????????????????%?信號帶寬
fl?=?f0-B/2;????????????????????%?信號起始頻率
fh?=?f0+B/2;????????????????????%信號截止頻率
Tr?=?50*10^（-6）;????????????????%?工作周期，脈沖重復周期？？
T?=?10*10^（-6）;?????????????????%?信號持續時間，要求大于孔徑渡越時間
????????????????????????????????%T和Tr的關系
fs=400*10^6;????????????????????%采樣頻率？
snr?=?[40;0;40;40];???????????????????%?信噪比？？？
theta0?=?0/180*pi;??????????????%?指向俯仰角，方位角為0
%?phi0?=?0/180*pi;????????????????%?指向方位角
M?=?16;??????????????????????????%?行陣元數為M?
%?Q?=?8;??????????????????????????%?列陣元數為M?
Ks?=?2;?????????????????????????%?信號數目
Theta?=[35/180*pi0/180*pi-40/180*pi60/180*pi];?????%?干擾信號方向
%?Phi?=[0/180*pi0/180*pi];?????%?干擾信號方向
c=3*10^8;
d=0.5*c/f0;
N=T*fs;???????????????????????%一個周期采樣點數，4000
Nr=Tr*fs;?????????????????????%%五個周期采樣點數
Nm?=?5;????????????????????????%?駐留周期個數
J?=?200;???????????????????????%?FFT的點數？？？FFT點數在頻域寬帶波束形成中起到什么作用？？
???????????????????????????????%當信號帶寬較小時，相對于固定采樣率，則信號所占的頻率點數越少，則無法較好的左頻帶劃分。
???????????????????????????????%將fft點數增多，則對同樣帶寬信號具有更詳細的表示，同時更容易做頻帶劃分。
fw=linspace（0fsJ+1）;??????????????%前nfft點是fft后的頻域頻率
????????????????????????????????????%%%即將信號對應的頻率點數，實際中這些點值有少部分點有值，對應帶通信號
????????????????????????????????????%%%還是變換到基帶進行處理，所以其頻帶
????????????????????????????????%%應當好奇的是這些點所對應的頻帶范圍始終在一個頻率周期以內即0—fs之內
????????????????????????????????????
fwn=[];????????????????????????????????%fft后落在帶寬內頻率
kn=[];
m=fix（f0/fs）;??????????????????%%此處m=7;
for?k=1:J???????????????????????????????%%%這里到底是用的帶通采樣，根據帶通采樣定理
????if?fw（k）>=f0-m*fs-B/2?&&?fw（k）<=f0-m*fs+B/2??%%判斷帶通信號在一個頻率周期內的數據是那哪些點，信號所占頻帶點數與信號帶寬跟采樣頻率之比相等。
????????fwn=[fwn?fw（k）];??????????????%fft后落在帶寬內頻點，保存這些數據，fwn對應頻帶的橫軸
????????kn=[kn?k];????????????????????%%kn為索引
????end
end
%%%%%對信號頻帶做子帶劃分和處理應僅對這部分完全的信號頻帶點做處理

G=length（kn）;?????????????????????????%?信號頻率帶寬所對應的點數，頻帶劃分，劃分了多少個子帶
NN?=?Nm*Nr/J;???????????????????%?頻域快拍數500個，
tic?????????????????????????????%%%%倒數計時器
%?產生陣列接收的寬帶數據
x=[];
for?nm=1:Nm
%?????phase_rand=randn（Ks1）;
????x2?=?LFMsource（Mf0flBfsTThetasnr2）;?%?期望信號的脈沖寬度為T，出現時間不定，跟噪聲功率相同，方向為0
????%?產生干擾信號
????x1?=?LFMsource（Mf0flBfsTrThetasnr1）;%干擾信號1，整個重復周期的數據，功率大，方向為40
%?????ft?=?0:1/fs:Tr-1/fs;????
%?????N1?=?length（ft）;???????????????????%?快拍數
%?????????????????????????%?調頻速率
%?????for?m=1:M
%?????for?n=1:N1
%?????????x3（mn）?=?10^（snr（3）/20）*exp（i*（2*pi*fl*（ft（n）-（m-1）*d/c*sin（Theta（3）））））;
%?????????????????????????????%%%%fl是否該改成f0
%?????end
%????end
%??????????????????????????%?調頻速率

????x3?=?LFMsource（Mf0flBfsTrThetasnr3）;
????x4?=?LFMsource（Mf0flBfsTrThetasnr4）;
%?????nn=1:Nr;
%?????for?n=1:Nr
%?????????fai（1n）?=?rand;
%?????end
%?????s（1:）?=?exp（i*2*pi*（fr*nn/fs+k*fai））;
%?????for?k=1:Kr
%?????????for?m=1:M
%?????????