SAR成像算法RDCSwk

经典的SAR成像算法合集，包含RD（距离多普勒）算法，CS（chirp scaling）算法，wk算法，点目标仿真。

 %% 正侧视点目标成像--CSA
 clc;
 clear;
 close all;
 %% 参数设置
 c=3e8;                         
 fc=3e9;                        
 lambda=c/fc;                  %载波频率（Hz）、波长0.1（m）
 Ralpha=1.2;  Aalpha=1.3;      %距离向、方位向过采样率
 %----------观察场景相关参数----------%
 H=8000;                       %平台运行高度
 Yc=6000;                      %场景中心线长度
 R0=sqrt（H^2+Yc^2）;
 Xmin=-200;                    %[Xmin Xmax]
 Xmax=200;                     %在合成孔径长度Ls基础上，额外方位向范围
 Yw=300;                      
  %----------天线孔径设置----------%
 La=1.5;
 theta=0.886*lambda/La;
 Ls=R0*theta;
 %----------慢时间参数设置----------%
 v=150;                        %平台运行速度
 Ts=Ls/v;
 Xwid=Ls+Xmax-Xmin;            %慢时间域时间窗长度
 Twid=Xwid/v;
 Ka=2*v^2/lambda/R0;           %方位向调频率
 Ba=abs（Ka*Ts）;                %方位向带宽
 PRF=Aalpha*Ba;
 PRT=1/PRF;
 dx=PRT;                       %方位向采样间隔
 N=ceil（Twid/dx）;
 N=2^nextpow2（N）;              %提高fft效率
 x=linspace（（Xmin-Ls/2）/v（Xmax+Ls/2）/vN）;
                               %慢时间域时间序列
 X=v*x;                        %慢时间域时间序列对应方位向距离
 PRT=Twid/N;                   %更新
 PRF=1/PRT;
 dx=PRT;
 %----------快时间参数设置----------%
 Tr=5e-6;                      %脉冲宽度
 Br=100e6;                     %带宽100MHz
 Kr=Br/Tr;                     %调频率
 Fs=Ralpha*Br;                 %距离向采样率
 dt=1/Fs;
 Rmin=sqrt（（Yc-Yw）^2+H^2）;
 Rmax=sqrt（（Yc+Yw）^2+H^2）;
 Rm=Rmax-Rmin+c*Tr/2;          %斜距测绘带宽
 M=ceil（2*Rm/c/dt）;            %采样点数 
 M=2^nextpow2（M）;       
 t=linspace（2*Rmin/c-Tr/22*Rmax/c+Tr/2M）;
                               %快时间域时间序列
 r=c*t/2;                      %快时间域时间序列对应斜距
 dt=（2*Rmax/c+Tr-2*Rmin/c）/M;  %更新
 Fs=1/dt;
 %----------目标参数设置----------%
 Ptarget=[0R01;50R0+1000.8;700R0+1000.8;];
 Ntarget=length（Ptarget）;
 %% 生成SAR回波
 s0=zeros（NM）;
 h=waitbar（0‘SAR回波生成‘）;
 for k=1:Ntarget
    R=sqrt（Ptarget（k2）^2+（X-Ptarget（k1））.^2）; 
    delay=2*R/c;               %距离徙动带来的距离时间延时
    Delay=ones（N1）*t-delay‘*ones（1M）;   %生成N*M的矩阵
    Phase=1j*pi*Kr*Delay.^2-1j*4*pi*fc*（R‘*ones（1M））/c;
    s0=s0+Ptarget（k3）*rectpuls（Delay/Tr）.*...
      rectpuls（（X-Ptarget（k1））‘*ones（1M）/Ls）.*exp（Phase）;
   %s0=s0+Ptarget（k3）*exp（Phase）.*（abs（Delay）<=Tr/2）.*（abs（（X-Ptarget（k1））‘*ones（1M））<=Lsar/2）;
    waitbar（k/Ntarget）;
 end
 close（h）;
 %----------绘制原始信号----------%
 h=waitbar（0‘绘制原始信号‘）;
 figure（1）;
 imagesc（r*1e-3Xabs（s0））;
 axis tight;
 xlabel（‘距离向（km）\rightarrow‘）;
 ylabel（‘\leftarrow方位向（m）‘）;
 title（‘原始信号‘）;
 close（h）;
 clear R delay Delay Phase;
 %% 方位向fft、补余RCMC（在距离多普勒域中进行补余RCMC）
 h=waitbar（0‘方位向fft、补余RCMC‘）;
 fr=（-M/2:M/2-1）*Fs/M;         %距离向频率
 fa=（-N/2:N/2-1）*PRF/N;        %方位向频率
 Srd=ftx（s0）;                  %距离多普勒域
 Km=Kr;                        %这里我们近似等于
 D_fa_v=sqrt（1-（c^2*fa.^2/4/v^2/fc^2））.‘*ones（1M）;
 Ssc=exp（1j*pi*Km*（1./D_fa_v-1）.*（（ones（N1）*t-2*R0/c./D_fa_v）.^2））;
 S1=Srd.*Ssc;
 figure（2）;
 imagesc（r*1e-3faabs（iftx（S1）））;
 axis tight;
 xlabel（‘距离向（km）\rightarrow‘）;
 ylabel

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5801  2017-12-28 20:45  CS\CSA.m

     文件       5871  2014-03-28 16:30  CS\CSA_basic.m

     文件        966  2014-03-25 14:22  CS\CSA_tu7_3.m

     文件        956  2014-04-06 15:16  CS\CSA_tu7_5.m

     文件       2812  2014-11-18 20:19  CS\CS_scsj.m

     文件       2808  2014-11-18 20:24  CS\cs_scsj2.m

     文件       4643  2014-11-14 09:25  CS\CS_theta.m

     文件       4718  2014-11-16 10:44  CS\CS_theta1.m

     文件        115  2014-03-23 14:42  CS\ftx.m

     文件        116  2014-03-23 14:42  CS\fty.m

     文件        119  2014-03-23 14:42  CS\iftx.m

     文件        120  2014-03-23 14:43  CS\ifty.m

     文件        481  2014-03-23 20:14  CS\islrfunc.m

     文件        270  2006-06-30 12:06  CS\plot_img.m

     文件        277  2014-03-23 20:14  CS\pslrfunc.m

     文件        115  2014-03-23 14:42  RD\ftx.m

     文件        116  2014-03-23 14:42  RD\fty.m

     文件        205  2014-11-18 14:57  RD\Hist.m

     文件        713  2014-11-18 14:58  RD\Histeq_f.m

     文件        119  2014-03-23 14:42  RD\iftx.m

     文件        120  2014-03-23 14:43  RD\ifty.m

     文件        481  2014-03-23 20:14  RD\islrfunc.m

     文件        666  2014-11-18 14:59  RD\medfilt_f.m

     文件       1045  2014-11-18 14:59  RD\Pinghua_f.m

     文件        335  2014-11-12 20:12  RD\plot_img.m

     文件       1651  2014-11-20 10:52  RD\poumian.m

     文件        277  2014-03-23 20:14  RD\pslrfunc.m

     文件       8287  2017-05-11 10:06  RD\RDA.m

     文件       7515  2014-10-28 10:35  RD\RDA_basic.m

     文件       2584  2014-11-20 09:14  RD\RDA_scsj.m

............此处省略31个文件信息