%%產生Stripmap?SAR的回波
clear?all
c=3e8;%光速
%%測繪帶區域
X0=204.8;%方位向[-X0X0]，條帶長度的“一半”
Rc=10000;%條帶中心參考距離
R0=480;%距離向[Rc-R0Rc+R0]條帶寬度的“一半”
lambda=0.02;%波長
fc=c/lambda;%載頻

%%距離向（Range）r/t?domain
Tr=40e-6;%LFM信號脈寬
Br=10e6;?%LFM信號帶寬（=Kr*Tr）Hz
Kr=Br/Tr;?%調頻斜率
Fsr=40e6;%徑向快時間采樣率，應>Br
?????????%還應盡量使下面的Nr＝2若干次冪，以利于FFT。
Nr=ceil（Fsr*4*R0/c）;
r=Rc+linspace（-R0R0Nr）;
t=2*r/c;%t域序列
dt=R0*4/c/Nr;%采樣周期
f=linspace（-1/2/dt1/2/dtNr）;%f域序列

%%方位向（AzimuthCross-Range）x/u?domain
Va=200;%SAR?平臺速度
Lsar=150;%合成孔徑長度
PRF=1000;%Hz應>多普勒帶寬Ba（k*Tsar）還應盡量使下面的Na＝2若干次冪，以利于FFT
Ta=1/PRF;
D=Va*Ta;%兩孔徑間距
Na=ceil（PRF*2*X0/Va）;
x=linspace（-X0X0Na）;%u域序列
u=x/Va;
du=2*X0/Va/Na;
fu=linspace（-1/2/du1/2/duNa）;%fu域序列
fdc=0;%Doppler調頻中心頻率
fdr=-2*（Va）^2/lambda/Rc;%Doppler調頻斜率
Ba=abs（fdr*Lsa