%%產生Stripmap?SAR的回波
clear?all

thetaT=0;%T平臺波束斜視角
thetaT=thetaT*pi/180;%rad
thetaR=0;%R平臺波束斜視角
thetaR=thetaR*pi/180;
c=3e8;%光速
fc=1.5e9;%載頻
lambda=c/fc;%波長

%%測繪帶區域
X0=200;%方位向[-X0X0]
Rtc=3000;?
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;?%調頻斜率
Nr=512;
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
v=100;%SAR?平臺速度
Lsar=300;%合成孔徑長度
Na=1024;
x=linspace（-X0X0Na）;%u域序列
u=x/v;
du=2*X0/v/Na;
fu=linspace（-1/2/du1/2/duNa）;%fu域序列
ftdc=v*sin（thetaT）;
ftdr=-（v*cos（thetaT））^2/lambda/Rtc;
frdc=v*sin（thetaR）;
frdr=-（v*cos（thetaR））^2/lambda/Rrc;
fdc=ftdc+frdc;%Doppler調頻中心頻率
fdr=ftdr+frdr;%Doppler調頻斜率
%％目標位置
Ntar=3;%目標個數
Ntar=1;
Ptar=[Rrc01]%距離向坐標方位向坐標sigma?????????????
??????????%Rrc+50-501
?????????%Rrc+50501];


%%產生回波
s_ut=zeros（NrNa）;
U=ones（Nr1）*