ISAR成像仿真matlab

clear all;
clc;
close all;
C=3e8;              %光速
ima=sqrt（-1）;         %虚数单位
%%场景中心参数
R0=2e4;             %场景中心斜距
Wr=0.12e3;          %场景宽度
%%雷达参数
fc=10e9;             %工作频率10GHz
lamda=C/fc;          %波长
H=5000;             %高度 5000m
SNR=-10;
%%天线参数
D=4;                %天线孔径长度
Lsar=lamda*R0/D;     %雷达合成孔径长度
%%距离向参数
Tr=10e-6;            %发射脉冲时宽
B=180e6;            %带宽
Kr=B/Tr;            %调频率
Fsr=200e6;          %距离采样率
dt=1/Fsr;            %距离采样时间
nrn=round（（Tr+2*Wr/C）*Fsr/2）*2;        %距离向的取样数
nrn_Tr=Fsr*Tr;                        %脉冲区间的取样数
Nfast_signal=exp（ima*pi*Kr*（[-nrn_Tr/2:nrn_Tr/2-1].‘/Fsr）.^2）;   %参考信号
Nfast=2.^ceil（log2（nrn+nrn_Tr））;  %为了达到线性卷积的要求而取两者之和
r_signal=[zeros（round（Nfast/2-nrn_Tr/2）1）;Nfast_signal;zeros（round（Nfast/2-nrn_Tr/2）1）];
%%方位向参数
Fa=100;              %方位采样率
PRF=500;            %脉冲重复频率
PRT=1/PRF;          %脉冲重复时间
V=100;
ar=5;
wt=0.05;
nan=round（PRF*lamda*R0/D/V/2）*2;     %方位向取样数
%%分辨率参数
DY=C/2/B;           %距离向分辨率
DX=D/2;             %方位向分辨率
%%目标参数
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‘）;
%雷达回波
nrn=2048;
nan=512;
tnan=[-nan/2:nan/2-1]*PRT;
tnrn=2*R0/C+[-nrn/2:nrn/2-1]‘/Fsr;
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