ISAR 飞机模型成像

clear all;
clc;
C=3e8;               %光速
ima=sqrt（-1）;        %虚数单位
%%场景中心参数
R0=2e4;              %场景中心斜距
%%雷达参数
fc=10e9;             %工作频率10GHz
lamda=C/fc;          %波长
SNR=-10;
%%天线参数
D=4;                 %天线孔径长度
Lsar=lamda*R0/D;     %雷达合成孔径长度
%%距离向参数
Tr=10e-6;            %发射脉冲时宽
B=180e6;             %带宽
Kr=B/Tr;             %调频率
Fsr=200e6;           %距离采样率
nrn_Tr=Fsr*Tr;       %脉冲区间采样数
nrn=2048;            %距离向采样数
tnrn=2*R0/C+[-nrn/2:nrn/2-1]‘/Fsr;  %距离向快时间
%%方位向参数
Fa=100;              %方位采样率
PRF=500;             %脉冲重复频率
V=100;               %飞机速度
ar=5;                %飞机加速度
nan=512;             %方位向采样数
tnan=[-nan/2:nan/2-1]/PRF;  %方位向慢时间
%%分辨率参数
DY=C/2/B;           %距离向分辨率
DX=DY;              %为保持距离分辨率和方位分辨率相同，图像比例一致
wt=lamda/（2*DX*（nan/PRF））;   %方位向分辨率计算出应该转多大的角度wt，wt等于角速度*积累时间
%%目标参数，飞机模样
Ntarget=11;       	%点目标个数
deltaX=15;
deltaY=15;
points=zeros（Ntarget3）;              
points（:1）=[2*deltaX;deltaX;deltaX;0;0;0;0;0;-deltaX;-deltaX;-2*deltaX];
points（:2）=[0;0;-2*deltaY;2*deltaY;deltaY;0;-deltaY;-2*deltaY;0;-2*deltaY;0];
points（:3）=ones（Ntarget1）;
figure;
h=stem（points（:1）points（:2）‘fill‘）;
grid on;
set（h‘Linestyle‘‘none‘）;
axis（[-50 50 -50 50]）
title（‘点目标真实图‘）;xlabel（‘x‘）;ylabel（‘y‘）;

%雷达回波
s=zeros（nrnnan）;
for n=1:nan
Rs=R0+V*tnan（n）+1/2*ar*tnan（n）^2;  
for m=1:Ntarget
   rp=sqrt（points（m1）^2+points（m2）^2）;  
   phi=atan2（points（m2）points（m1））+wt*tnan（n）;
   R=sqrt（Rs^2+rp^2+2*Rs*rp*sin（phi））;
   td=2*R/C;
   win_r=（abs（tnrn-td）<=Tr/2）;     %距离包络（矩形窗函数）
   s_mn=points（m3）.*exp（ima*pi*Kr*（tnrn-td）.^2）.*win_r.*exp（-ima*4*pi*R/lamda）;
   s（:n）=s（:n）+s_mn;
end
s（: