sar 点目标分析

sar点目标质量分析，峰值旁瓣比，分辨率等。

%%========================================================
clear;clc;close all;
%%========================================================
%%参数——常量
C=3e8;                           %电磁波速度

%%参数——雷达规格
Fc=3e9;                          %载波频率
lambda=C/Fc;                     %波长 

%%参数——目标区域
X0=500;                          %方位向范围：[-X0X0]
Y=10000;                         %测绘中心Y坐标
Y0=500;                          %成像宽度：2*Y0

%%参数
theta=0;                         %平台波束斜视角
theta=theta*pi/180;              %rad
V=100;                           %平台速度：100 m/s
H=4000;                          %平台高度：4 km
Rb=sqrt（H^2+Y^2）;
Rc=Rb/cos（theta）;                %目标区域距离向范围：[Rc-Y0Rc+Y0]

%%参数——雷达天线
D=4;                             %天线方位向尺寸
Lsar=lambda*Rc/D;                %SAR合成孔径长度
Tsar=Lsar/V;                     %SAR合成孔径时间

%%参数——慢时间域
Ka=-2*（V*cos（theta））^2/lambda/Rc;            %doppler调频斜率
Ba=abs（Ka*Tsar）;                %doppler调频频宽
PRF=Ba;                         %脉冲重复频率
PRT=1/PRF;                      %脉冲重复时间
ds=PRT;                         %慢时域采样间距
                                %为了进行FFT，对PRT、PRF、ds进行运算并更新
Nslow=ceil（（2*X0+Lsar）/V/ds）;   %慢时域采样点数
Nslow=2^nextpow2（Nslow）;        %使采样点数为2的幂；N点

%Nslow=1024;

sn=linspace（-（X0+Lsar/2）/V（X0+Lsar/2）/VNslow）;%慢时域离散时间序列，1*N矩阵
PRT=（2*X0+Lsar）/V/Nslow;        %更新PRT、PRF、ds
PRF=1/PRT;
ds=PRT;
%%参数——快时域
Tr=5e-6;                        %Chirp脉冲持续时间：5 us
Br=50e6;                        %chirp信号调频带宽：50 MHz
Kr=Br/Tr;                       %chirp信号调频斜率
Fsr=2*Br;                       %快时域采样频率
dt=1/Fsr;                       %快时域采样间距
Rmin=sqrt（（Y-Y0）^2+H^2）;        %最小斜距
Rmax=sqrt（（Y+Y0）^2+H^2+（Lsar/2）^2）;     %最大斜距
Nfast=ceil（2*（Rmax-Rmin）/C/dt+Tr/dt）;   %快时域采样点数
Nfast=2^nextpow2（Nfast）;                %使采样点数为2的幂；M点
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）;%快时域离散时间序列，1*M矩阵
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;        %更新Fsr、dt
Fsr=1/dt;

%%参数——分辨率
DY=C/2/Br;                          %距离向分辨率:DY=C/2/Br;
DX=D/2;                             %方位向分辨率:DX=D/2;

%%参数——点目标
Ntarget=3;                          %点目标个数

%格式[方位向 距离向水平投影 目标散射系数]
%标准点位置
Ptarget=[0Y1
         62*DXY1
         75Y+25*DY1 ];
%分辨率点位置
%Ptarget=[0Y1
%         3*DXY1
%         0Y+2*DY1 ];

disp（‘参数:‘）
disp（‘快时域采样数：‘）;disp（Nfast）
disp（‘慢时域采样数：‘）;disp（Nslow）
disp（‘距离向分辨率：‘）;disp（DY）
disp（‘方位向分辨率：‘）;disp（DX）     
disp（‘合成孔径长度：‘）;disp（Lsar）     
%%========================================================
%%产生单、多点回拨回波信号，一个N*M矩阵
K=Ntarget;                       
N=Nslow;                         
M=Nfast;                        
T=Ptarget;                      
Srnm=zeros（NM）;
for k=1:1:K                      
    sigma=T（k3）;                 %获取点目标散射系数
    Dslow=sn*V-T（k1）;            %慢时域中雷达与目标方位向距离，1*N矩阵
    R=sqrt（Dslow.^2+T（k2）^2+H^2）;%斜距，1*N矩阵
    tau=2*R/C;                    %时间延迟，1*N矩阵
    Dfast=ones（N1）*tm-tau‘*ones（1M）;                  %[t-2*R/C]项，N*M

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

     文件       5366  2009-05-15 18:57  pangbanbi1.m

     文件       6315  2009-05-21 13:50  chengxiang.m

     文件       5324  2009-05-21 14:37  dingweixianshi.m

     文件         68  2009-04-10 13:32  ftx.m

     文件         72  2009-04-10 13:32  fty.m

     文件         71  2009-04-10 13:32  iftx.m

     文件         75  2009-04-10 13:32  ifty.m

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

                17291                    7