sarecho.asv

条带式合成孔径雷达的cs。RDA。RMA算法 作者何志华

%%SAR回波
%Stripmap 模式下SAR的回波生成
%NUDT,Hezhihua
%E-mail:[e-mail]
%last update:23/6/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=10;%载频10(1/m)(1.5GHz)
Nx=512;
Xc=3000;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=3;%目标个数
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

%%产生相关的数据并存储
p0_x=exp(j*pi*a*(x-Xc).^2).*(abs(x-Xc)<Tx/2);%距离向LFM信号
p0_kx=fftshift(fft(fftshift(p0_x)));
p0_y=exp(-j*pi*b*y.^2).*(abs(y)<Ty/2);%方位向LFM信号
p0_ky=fftshift(fft(fftshift(p0_y)));
 
s_kxy=fftshift(fft(fftshift(s_xy)));%对raw data距离向FFT
sxc_kxy=s_kxy.*(conj(p0_kx).'*ones(1,Ny));
sxc_kxky=fftshift(fft(fftshift(sxc_kxy).')).';%距离压缩后的2D频域信号
sxc_xy=fftshift(ifft(fftshift(sxc_kxy)));%距离压缩后的信号

sxc_xky=fftshift(fft(fftshift(sxc_xy).')).';%距离压缩后，距离-多普勒域
sxcyc_xky=sxc_xky.*(ones(Nx,1)*conj(p0_ky));
sxcyc_xy=fftshift(ifft(fftshift(sxcyc_xky).')).';%距离，方位压缩后的信号

save sarecho.mat Tx Bx a Ty By b Xc...
                             dx dy ...
                             Nx Ny kxc...
                             x y kx ky...
                             p0_x p0_kx...
                             p0_y p0_ky...
                             sxc_xy...
                             sxc_xky...
                             sxc_kxky...
                             sxcyc_xy
clear all