RTPC sar面目标成像

基于SAR RD算法的面目标成像，有关于面目标的具体布置方式，成像算法，很不错，望大家多多指教。

%%========================================================
%%RTPC方法：SAR原理的直观解释
%%========================================================
clear;clc;close all;
%%========================================================
%%参量--常数
C=3e8;                           %传播速度（光速）
%%参量--雷达特性
Fc=1e9;                          %载频
lambda=C/Fc;                     %波长
%%参数--目标区域
Xmin=0;                          %方位向坐标范围[XminXmax]
Xmax=50;                  
Yc=10000;                        %测绘带中心
Y0=500;                          %距离向坐标范围[Yc-Y0Yc+Y0]
                                 %测绘带宽2*Y0
%%参数--平台信息
V=100;                           %SAR 飞行速度 100 m/s
H=5000;                          %高度 5000 m
R0=sqrt（Yc^2+H^2）;               %SAR距测绘带中心距离
%%参数--天线
D=4;                             %天线孔径最大尺寸
Lsar=lambda*R0/D;                %合成孔径长度：Lsar决定于天线运动过程中所能接收到的来自同一目标单元的回波信号的最大作用范围，
                                 %             它等于真实天线所能覆盖的最大范围                                             
Tsar=Lsar/V;                     %合成孔径时间
%%参数--慢时间域
Ka=-2*V^2/lambda/R0;            %多普勒调频斜率
Ba=abs（Ka*Tsar）;                %多普勒调频带宽
PRF=Ba;                        %PRF：（PRF的取值范围：1.1BdPRT=1/PRF;                     %PRT
ds=PRT;                         %慢时间域抽样间隔
Nslow=ceil（（Xmax-Xmin+Lsar）/V/ds）; %慢时间域抽样数
Nslow=2^nextpow2（Nslow）;          %取2的阶乘数便于FFT计算
sn=linspace（（Xmin-Lsar/2）/V（Xmax+Lsar/2）/VNslow）;%慢时间域离散时间阵（SAR对所有点有作用的时间）
PRT=（Xmax-Xmin+Lsar）/V/Nslow;    %更新
PRF=1/PRT;
ds=PRT;
%%参数--快时间域
Tr=5e-6;                         %脉冲持续 10us
Br=30e6;                        %LFM信号带宽 30MHz
Kr=Br/Tr;                        %LFM斜率
Fsr=3*Br;                        %快时间抽样频率：理论上1.5倍以上的抽样频率即可，但是在matlab中，3倍Br时波形仍成锯齿状
dt=1/Fsr;                         %快时间抽样间隔
Rmin=sqrt（（Yc-Y0）^2+H^2）;
Rmax=sqrt（（Yc+Y0）^2+H^2+（Lsar/2）^2）;
Nfast=ceil（2*（Rmax-Rmin）/C/dt+Tr/dt）;%快时间域抽样个数
Nfast=2^nextpow2（Nfast）;                   %取2的阶乘便于FFT计算
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）; %快时间域离散时间阵
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;    %更新
Fsr=1/dt;
%%参数--分辨率
DY=C/2/Br;                          %距离分辨率
DX=D/2;                            %方向分辨率
%%参数--点目标个数
Ntarget=1;                           %点目标个数
%目标参数格式 [x y 散射系数]
Ptarget=[XminYc1                    %目标的坐标
             XminYc+DY1
              XminYc+2*DY1
               XminYc+3*DY1
              XminYc+4*DY1
               XminYc+5*DY1
                XminYc+6*DY1
                 XminYc+7*DY1
                  XminYc+8*DY1
                   XminYc+9*DY1
                   XminYc+10*DY1
                    XminYc+11*DY1
                     XminYc+12*DY1
                      XminYc+13*DY1
                       XminYc+14*DY1
                        XminYc+15*DY1
           Xmin+DXYc1                    %目标的坐标
             Xmin+DXYc+DY1
              Xmin+DXYc+2*DY1
               Xmin+DXYc+3*DY1
              Xmin+DXYc+4*DY1
   

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

     文件       9738  2010-01-14 14:08  sar_p.m

----------- ---------  ---------- -----  ----

                 9738                    1