ISAR成像matlab仿真

clear;
clc;
c = 3e8;
Tp = 2e-6;                %发射脉冲宽度
fc = 10e9;                %载频10GHz
lambda = c/fc;          %波长0.03m
B = 200e6;               %发射带宽200MHz
gama = B/Tp;           %距离向调频率γ=100MHz/s
Fs = 280e6;             %采样频率280MHz
N_r = 2^（ceil（log2（Tp*Fs）））*2;%距离向采样点数
pulsewidth = floor（Tp*Fs/2）*2;
t_pulse_r = （-pulsewidth/2+1:pulsewidth/2）/Fs;
t_N = （-N_r/2+1:N_r/2）/Fs;%此仿真系统分辨率：0.75m
PRF = 100;%脉冲重复频率
PRT = 1/PRF;%脉冲重复时间
w = 0.02;%角速度
N_a = 128;%方位向采样点数
t_pulse_a = （1:N_a）*PRT;

%waveform generation
target_num = 5;
Rs = 1e4;
r_target_x = [-10-1001010];
r_target_y = [Rs-10 Rs+10 Rs Rs-10 Rs+10];
figure（1）;
h = stem（r_target_xr_target_y‘fill‘）;
grid on;
set（h‘Linestyle‘‘none‘）;
axis（[-20 20 9980 10020]）;
title（‘点目标真实图‘）;xlabel（‘x‘）;ylabel（‘y‘）;

s = zeros（N_rN_a）;
for m = 1:N_a
    for n = 1:target_num
        tn = 2*（r_target_x（n）*w*t_pulse_a（m）+r_target_y（n）-Rs）/c;
        nn = floor（tn*Fs/2）*2;
        y_1 = exp（1i*pi*gama*t_pulse_r.^2）.*exp（1i*2*pi*fc*（t_pulse_r-tn））;
        y1 = zeros（N_r1）;
        y1（end/2+nn/2-pulsewidth/2+1:end/2+nn/2+pulsewidth/2） = y_1;
        s（:m） = s（:m） + y1;
    end
end
s = s.*exp（-1i*2*pi*fc*t_N.‘）;
figure（2）;
contour（abs（s））;
title（‘回波波形‘）;xlabel（‘方位向‘）;ylabel（‘距离单元‘）;

%距离压缩
s_ref = exp（1i*pi*gama*t_pulse_r.^2）;
s_ref_f = zeros（N_r1）;
s_ref_f（end/2-pulsewidth/2+1:end/2+pulsewidth/2） = s_ref;
s1 = zeros（N_rN_a）;
for m = 1:N_a
    s1（:m） = ifftshift（ifft（ifftshift（fftshift（fftshift（fft（s（:m））