%%?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