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=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=2;                           %number of targets
%format [x y reflectivity]
Ptarget=[XminYc1                    %position of targets
              XminYc+10*DY1
              Xmin+20*DXYc+50*DY1];  
disp（‘Parameters:‘）
disp（‘Sampling Rate in fast-time domain‘）;disp（Fsr/Br）
disp（‘Sampling Number in fast-time domain‘）;disp（Nfast）
disp（‘Sampling Rate in slow-time domain‘）;disp（PRF/Ba）
disp（‘Sampling Number in slow-time domain‘）;disp（Nslow）
disp（‘R

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

     文件       5748  2010-01-27 09:16  rawdata\dian.m

     文件         78  2010-01-31 20:05  rawdata\ftx.m

     文件         79  2010-01-31 20:05  rawdata\fty.m

     文件       3593  2010-01-28 11:03  rawdata\fwxz.m

     文件       3477  2010-01-31 20:42  rawdata\hh.m

     文件         81  2010-01-31 20:05  rawdata\iftx.m

     文件         84  2010-01-31 20:05  rawdata\ifty.m

     文件       3915  2010-01-28 10:44  rawdata\jlxz.m

     文件      14194  2010-01-27 15:20  rawdata\mian\1.jpg

     文件      24462  2010-01-27 15:27  rawdata\mian\2.jpg

     文件     112202  2010-01-27 15:27  rawdata\mian\3.jpg

     文件     122494  2010-01-27 15:28  rawdata\mian\4.jpg

     文件      45993  2010-01-27 16:30  rawdata\mian\5.jpg

     文件         78  2010-01-13 15:25  rawdata\mian\ftx.m

     文件         79  2010-01-13 15:25  rawdata\mian\fty.m

     文件         81  2010-01-13 15:25  rawdata\mian\iftx.m

     文件         84  2010-01-13 15:25  rawdata\mian\ifty.m

     文件       6012  2010-01-27 10:21  rawdata\mian\mian.m

     文件        452  2010-01-30 19:10  rawdata\mian\RCSmat.m

     文件     276084  2010-01-28 12:03  rawdata\mian\新建 Microsoft Office Word 文档.docx

     文件      14146  2010-01-27 15:23  rawdata\mian\旋转 1.jpg

     文件      45772  2010-01-27 16:30  rawdata\mian\旋转 5.jpg

     文件         78  2010-01-13 15:25  rawdata\mian1\ftx.m

     文件         79  2010-01-13 15:25  rawdata\mian1\fty.m

     文件         81  2010-01-13 15:25  rawdata\mian1\iftx.m

     文件         84  2010-01-13 15:25  rawdata\mian1\ifty.m

     文件       4511  2010-01-30 19:10  rawdata\mian1\mian1.m

     文件        462  2010-01-30 19:09  rawdata\mian1\RCSmat.m

     文件      12853  2010-01-28 16:38  rawdata\mian1\untitled.jpg

     文件      12739  2010-01-28 16:40  rawdata\mian1\旋转 untitled.jpg

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