SAR CS成像算法matlab

%% Chirp Scaling 算法CSA
%NUDT ?
%2006/7/8

clear allclc
close all
%%距离向参数range:x domain
%x=c*t/2kx=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（-X0X0Nx）;%x域序列:Xc-X0~Xc+X0
dx=2*X0/Nx;
kx=linspace（-1/dx/21/dx/2Nx）;%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（-Y0Y0Ny）;%y域序列:-Y0~Y0
dy=2*Y0/Ny;
ky=linspace（-1/dy/21/dy/2Ny）;%ky域序列

%%目标几何关系target geometry
%x坐标y坐标复后向散射系数 
Ptar=[Xc01+0j              
          Xc+50-501+0j
          Xc+50501+0j
          Xc-50-501+0j
          Xc-50501+0j];

%%生成SAR正交解调后的回波数据
s_xy=zeros（NxNy）;
N=size（Ptar1）;%目标个数
h = waitbar（0‘SAR回波生成‘）;
for i=1:1:N
    xn=Ptar（i1）;yn=Ptar（i2）;sigma=Ptar（i3）;%提取每个目标的信息
    X=x.‘*ones（1Ny）;%扩充为矩阵
    Y=ones（Nx1）*y;%扩充为矩阵
    DX=sqrt（xn^2+（Y-yn）.^2）;%中间变量
    phase=pi*a*（X-DX）.^2-2*pi*kxc*DX;%回波相位
    s_xy=s_xy+sigma*exp（j*phase）.*（abs（X-DX）    waitbar（i/N）
end
close（h）

%%数据准备
phi0=-x‘*sqrt（kxc^2-ky.^2）;
phi1=-kxc*x‘*（1./sqrt（kxc^2-ky.^2））;
phi2=1/2*x‘*（ky.^2./（kxc^2-ky.^2）.^1.5）;
Cs=ones（Nx1）*（kxc./sqrt（kxc^2-ky.^2）-1）;
Ks=1./（1/a-2*phi2）;

%%CSA:7 steps
s_xky=ff