ISAR转台成像matlab程序

%%%%%%参数设置%%%%%%
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=[-3 2.5 0 0];%散射点横坐标
% Ry=[80 80 78.5 16];%散射点纵坐标
 Ry=[17 13.5 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=linspace（0Tpndistance）;%脉内时间
sr=zeros（npulsenfast）;
ss=zeros（npulsenfast）;
s0=zeros（npulsenfast）;
for n=1:npulse