sar_sig.m

利用matlab编程语言

%%SAR回波
%Stripmap 模式下SAR的回波生成
%NUDT,Hezhihua
%E-mail:[e-mail]
%last update:6/7/2005

clear all
%%range:x domain
%x=c*t/2,kx=2*f/c
Tx=200;%时宽200m(1.33us)
Bx=1;%带宽1(1/m)(150MHz)
a=Bx/Tx;%调频斜率
kxc=4;%载频4(1/m)
Nx=512;
Xc=1200;X0=150;
x=Xc+linspace(-X0,X0,Nx);%x域序列:Xc-X0~Xc+X0
dx=2*X0/Nx;
kx=linspace(-1/dx/2,1/dx/2,Nx);%kx域序列

%%cross-range:y domain
Ty=300;%时宽300m,合成孔径长度
By=1;%带宽1(1/m)
b=kxc/Xc;%调频斜率 b=By/Ty=kxc/Xc
Ny=1024;
Y0=200;
y=linspace(-Y0,Y0,Ny);%y域序列:-Y0~Y0
dy=2*Y0/Ny;
ky=linspace(-1/dy/2,1/dy/2,Ny);%ky域序列

%%target geometry
Ntar=1;%目标个数
Ptar=[Xc,0,1+0j  %x坐标,y坐标,复后向散射系数             
          Xc+50,-50,1+0j
          Xc+50,50,1+0j];
      
%%SAR raw data
s_xy=zeros(Nx,Ny);
for i=1:1:Ntar
    xn=Ptar(i,1);yn=Ptar(i,2);sigma=Ptar(i,3);
    X=x'*ones(1,Ny);%扩充为矩阵
    Y=ones(Nx,1)*y;
    DX=X-sqrt(xn^2+(Y-yn).^2);
    phase=pi*a*DX.^2-2*pi*kxc*sqrt(xn^2+(Y-yn).^2);
    s_xy=s_xy+sigma*exp(j*phase).*(abs(DX)<Tx/2).*(abs(Y-yn)<Ty/2);%SAR回波信号
end