clear?all;??
clc;??

%---------------------------------start-----------------------------------%??
%-------------------------------基本參數(shù)設(shè)置-------------------------------%??
thetaT=0;??????????????????????????????%T平臺波束斜視角
thetaT=thetaT*pi/180;??????????????????%rad弧度
thetaR=0;??????????????????????????????%R平臺波束斜視角
thetaR=thetaR*pi/180;??????????????????%rad弧度
c=3e8;?????????????????????????????????%光速
fc=1.5e9;????????????????????%載波頻率1.5G
lambda=c/fc;?????????????????%載波的波長
%%-----------測繪帶區(qū)域----------------
X0=50;??????????%方位向[-X0X0]，范圍有自己根據(jù)實際情況確定
Rtc=3000;???????%載波發(fā)射距離
Rrc=3000;???????%載波接收距離
Rc=（Rtc+Rrc）/2;?%載波距離向
R0=150;?????????%距離向[Rc-R0Rc+R0]
%%-----------距離向（Range）r/t?domain---
Tr=1.5e-6;??%LFM信號脈寬?1.5us?（200m）

Br=150e6;?????%LFM信號帶寬?150MHz
Kr=Br/Tr;?????%調(diào)頻斜率
Nr=512;???????%快時間采樣點數(shù)
?
%+++++距離域序列+++++
r=Rc+linspace（-R0R0Nr）;????%距離域序列

t=2*r/c;???????????????????????%距離域t值對應(yīng)
dt=R0*4/c/Nr;??????????????????%快時間采樣周期
%+++++頻率域序列+++++
f=linspace（-1/2/dt1/2/dtNr）;?%f域序列

%%----------方位向（AzimuthCross-Range）x/u?domain----
v=100;???????????%SAR?平臺移動速度
%{
??合成孔徑觀點
??通過PRF進行方位向全息圖的采樣;
??采樣間隔必須固定不變，因此要求PRF與地速成正比.
%}
Lsar=300;?????????????%合成孔徑長度
Na=1024;??????????????%?慢時間采樣點數(shù)
%+++++u域序列+++++
x=linspace（-X0X0Na）;?????????%u域序列
u=x/v;?????????????????????????%u域序列t值對應(yīng)
du=2*X0/v/Na;??????????????????%慢時間采樣間隔
%+++++fu域序列+++++
fu=linspace（-1/2/du1/2/duNa）;%fu域序列
ftdc=v*sin（thetaT）;???????????%SAR-T平臺波束速度
ftdr=-（v*cos（thetaT））^2/lambda/Rtc;
frdc=v*sin（thetaR）;???????????%SAR-R平臺波束速度?
frdr=-（v*cos（thetaR））^2/lambda/Rrc;
fdc=ftdc+frdc;????????????????%Doppler調(diào)頻中心頻率
fdr=ftdr+frdr;????????????????%Doppler調(diào)頻斜率
%-------目標位置-----------
Ntar=3;%目標個數(shù)
Ptar=[??Rrc?????0???1?%參數(shù)對應(yīng)：距離向坐標方位向坐標sigma?????????????
??????Rrc+50???-50??1
??????Rrc+50????50??1];
??%%-----------LFM產(chǎn)生回波---------
s_ut=zeros（NrNa）;??%生成m×n的double類零矩陣，設(shè)定數(shù)值內(nèi)存
%{
??概念普及：
??zeros:
??zeros（mn）產(chǎn)生m×n的double類零矩陣，zeros（n）產(chǎn)生n×n的全0方陣。
%}
U=ones（Nr1）*u;????????????%擴充為矩陣
T=t‘*ones（1Na）;???????????%?將快時間拓展為?Na?列
%{
??概念普及：
??ones:
??????ones函數(shù)——生成全1陣
??zeros（mn）產(chǎn)生m×n的double類零矩陣，zeros（n）產(chǎn)生n×n的全0方陣。
%}
for?i=1:1:Ntar
????rn=Ptar（i1）;??????????%?目標距離向坐標
????xn=Ptar（i2）;??????????%?目標方位向坐標
????sigma=Ptar（i3）;???????%?目標RCS
????rtn=rn+Rtc-Rrc;????????
????RT=sqrt（rtn^2+（rtn*tan（thetaT）+xn-v*U）.^2）;?????%?發(fā)射?目標斜距
????RR=sqrt（rn^2+（rn*tan（thetaT）+xn-v*U）.^2）;???????%?接收?目標斜距
????R=RT+RR;
????DT=T-R/c;
????phase=pi*Kr*DT.^2-2*pi/lambda*R;
????s_ut=s_ut+sigma*exp（1j*phase）.*（abs（DT）end;