ISAR 对4个点目标进行仿真成像

%%%%%%参数设置%%%%%%
%--对4目标成像-----
%---采用解线频调方法成像-----

clc
clear all
Tp=4e-6;    %脉冲宽度
B=300e6;      %信号带宽
k=B/Tp;     %调频斜率
Fs=30e6;        %采样频率
prf=50;         %脉冲周期
c=3.0e8;
lamda=0.05;     %波长
w=0.05;             %转速
Rx=[-13 5 9 -10];        %散射点横坐标
 Ry=[5 10 7 -10];        %散射点纵坐标
R=sqrt（Rx.^2+Ry.^2）;%散射点到中心点距离
theta0=atan2（RyRx）;%散射点角度
Rmax=10100;     %距离门最大距离
Rmin=9900;      %距离门最小距离
Rref=10000;     %参考距离
pr=c/2/B;       %距离分辨率
pa=pr/2;        %要求的方位分辨率
delta_theta=lamda/2/pa %要求的转角
t_n=delta_theta/w%要求的观测时间
npulse=ceil（t_n*prf）;%要求的积累脉冲数
nfast=round（（2*Rmax/c-2*Rmin/c+Tp）*Fs）;%距离门采样数
ndistance=round（Tp*Fs）;%目标距离向采样数
t=linspace（2*Rmin/c-Tp/22*Rmax/c+Tp/2nfast）;%距离门内时间
tn=lin