clear;
clc;
c?=?3e8;
Tp?=?2e-6;????????????????%發(fā)射脈沖寬度
fc?=?10e9;????????????????%載頻10GHz
lambda?=?c/fc;??????????%波長(zhǎng)0.03m
B?=?200e6;???????????????%發(fā)射帶寬200MHz
gama?=?B/Tp;???????????%距離向調(diào)頻率γ=100MHz/s
Fs?=?280e6;?????????????%采樣頻率280MHz
N_r?=?2^（ceil（log2（Tp*Fs）））*2;%距離向采樣點(diǎn)數(shù)
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;%此仿真系統(tǒng)分辨率：0.75m
PRF?=?100;%脈沖重復(fù)頻率
PRT?=?1/PRF;%脈沖重復(fù)時(shí)間
w?=?0.02;%角速度
N_a?=?128;%方位向采樣點(diǎn)數(shù)
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（‘點(diǎn)目標(biāo)真實(shí)圖‘）;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））