SAR动目标回波模拟与成像

SAR运动目标回波模拟以及对回波数据成像的仿真

%%========================================================
clear;clc;close all;
%%========================================================
%%Parameter--constant
C=3e8;                           %propagation speed
%%Parameter--radar characteristics
Fc=1e9;                          %carrier frequency 1GHz
lambda=C/Fc;                     %wavelength 
%%Parameter--target area
Xmin=0;                          %target area in azimuth is within[XminXmax]
Xmax=50;                  
Yc=10000;                        %center of imaged area
Y0=500;                          %target area in range is within[Yc-Y0Yc+Y0]
                                 %imaged width 2*Y0
%%Parameter--orbital information
V=100;                           %SAR velosity 100 m/s
H=5000;                          %height 5000 m
R0=sqrt（Yc^2+H^2）;
%%Parameter--antenna
D=4;                            %antenna length in azimuth direction
Lsar=lambda*R0/D;               %SAR integration length
Tsar=Lsar/V;                     %SAR integration time
%%Parameter--slow-time domain
Ka=-2*V^2/lambda/R0;            %doppler frequency modulation rate
Ba=abs（Ka*Tsar）;                %doppler frequency modulation bandwidth
PRF=2*Ba;                        %pulse repitition frequency
PRT=1/PRF;                     %pulse repitition time
ds=PRT;                         %sample spacing in slow-time domain
Nslow=ceil（（Xmax-Xmin+Lsar）/V/ds）; %sample number in slow-time domain
% Nslow=2^nextpow2（Nslow）;          %for fft
sn=linspace（（Xmin-Lsar/2）/V（Xmax+Lsar/2）/VNslow）;%discrete time array in slow-time domain
PRT=（Xmax-Xmin+Lsar）/V/Nslow;    %refresh
PRF=1/PRT;
ds=PRT;
%%Parameter--fast-time domain
Tr=5e-6;                         %pulse duration 10us
Br=30e6;                        %chirp frequency modulation bandwidth 30MHz
Kr=Br/Tr;                        %chirp slope
Fsr=3*Br;                        %sampling frequency in fast-time domain
dt=1/Fsr;                         %sample spacing in fast-time domain
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）;%sample number in fast-time domain
% Nfast=2^nextpow2（Nfast）;                   %for fft
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）; %discrete time array in fast-time domain
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;    %refresh
Fsr=1/dt;
%%Parameter--resolution
DY=C/2/Br;                          %range resolution
DX=D/2;                            %cross-range resolution
%%Parameter--point targets
Ntarget=1;                           %number of targets
%format [x y reflectivity]
Ptarget=[XminYc1
                 XminYc+10*DY1
              Xmin+20*DXYc+50*DY1];  
disp（‘参数:‘）
disp（‘Sampling Rate in fast-time domain‘）;disp（Fsr/Br）
disp（‘快时间域采样点数‘）;disp（Nfast）
disp（‘Sampling Rate in slow-time domain‘）;disp（PRF/Ba）
disp（‘慢时间域采样点数‘）;disp（Nslow）
disp（‘距离分辨率‘）;disp（DY）
disp（‘方位分辨率‘）;disp（DX）     
disp（‘合成孔径长度‘）;disp（Lsar）     
disp（‘目标位置【（0-5

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

     文件       4856  2010-03-13 07:53  movingtarget\dian.m

     文件        696  2010-03-14 19:35  movingtarget\dianmian.m

     文件        203  2010-03-13 09:51  movingtarget\dmstap.m

     文件         78  2010-03-08 09:44  movingtarget\ftx.m

     文件         79  2010-03-08 09:44  movingtarget\fty.m

     文件         81  2010-03-08 09:44  movingtarget\iftx.m

     文件         84  2010-03-08 09:44  movingtarget\ifty.m

     文件       5685  2010-03-14 19:30  movingtarget\jdian.m

     文件       5558  2010-03-13 08:00  movingtarget\mdian.m

     文件       5643  2010-03-14 19:35  movingtarget\mdian1.m

     文件       3686  2010-03-14 19:03  movingtarget\mmian.m

     文件       5224  2010-06-14 21:37  movingtarget\tdian.m

     目录          0  2012-10-03 18:49  movingtarget

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

                31873                    13