%SAR.m
%產生SAR的同波
clear?all
thetaT=0;?%波束寬度角
thetaT=thetaT*pi/180;%轉換成弧度
thetaR=0;?%波束高度角
thetaR=thetaR*pi/180;
c=3e8;?%光速
fc=10e9;?%載頻
lambda=c/fc;?%波長
%測繪帶區域
X0=200;%方位向【-X0，X0】
Rtc=3000;
Rrc=3000;
Rc=（Rtc+Rrc）/2;
R0=150;?%距離向[Rc-R0，Rc+R0】
%距離向（Range），快時間域
Tr=1.5e-6;??%LFM信號脈寬1.5us（200m）
Br=1.5e8;?%LFM信號帶寬150MHz
Kr=Br/Tr;???%調頻斜率
Nr=512;?????%快時間采樣點數
r=Rc+linspace（-R0R0Nr）;?%瞬時斜距
t=2*r/c;?%時域序列
dt=R0*4/c/Nr;?%采樣周期
f=linspace（-1/2/dt1/2/dtNr）;??%頻域序列
%方位向（Azimuth，Cross—Range），慢時間域
v=100;?%SAR平臺速度
D=1;%天線方位尺寸
Lsar=lambda*Rc/D;?%合成孔徑長度
Na=1024;?%快時間采樣點數
x=linspace（-X0X0Na）;?%慢時域序列
u=x/v;?%慢時間
du=2*X0/v/Na;%采樣周期
fu=linspace（-1/2/du1/2/duNa）;?%慢頻域序列
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;?%距離走動增量
fdr=ftdr+frdr;?%多普勒調頻率
%目標位置
Ntar=3;?%目標個數
Ptar=[Rrc01
??????Rrc+50-501
??????Rrc+50501];?%距離向坐標，方位向坐標，目標散射系數（sigma）
%產生同波
s_ut=zeros（NrNa）;
U=ones（Nr1）*u;?%擴充為矩陣
T=t‘*ones（1Na）;
for?i=1:1:Ntar
??????rn=Ptar（i1）;xn=Ptar（i2）;sigma=Ptar（i3）;%獲取點目標散射系數
??????rtn=rn+Rtc-Rrc;
??????RT=sqrt（rtn^2+（rtn*tan（thetaT）+xn-v*U）.^2）;
??????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（1i*phase）.*（abs（DT）?end;
%距離壓縮
p0_t=exp（1i*pi*Kr*（t-2*Rc/c）.^2）.*（abs（t-2*Rc/c）p0_f=fftshift（fft（fftshift（p0_t）））;
s_uf=fftshift（fft（fftshift（s_ut）））;?%距離向FFT
src_uf=s_u