%%雷達成像處理脈沖壓縮的仿真
close?all;
clear;
clc;
T=10e-6;
B=60e6;
k=B/T;
fs=100e6;
N=fs*T;
deltat=1e-6;
deltaN=deltat*fs;
t_ser=（0:T/N:T*（1-1/N））-T/2;????%碼元持續時間序列
tstart?=?5000;
tlen?=?16384;???????%回波信號長度
techo_ser=（1:tlen）/fs;

w_echo=（0:fs/tlen:fs-1/tlen）-fs/2;
w_ser=（0:fs/N:fs-1/N）-fs/2;
xlfm=exp（pi*1i*k*t_ser.^2）;
%?figure;
%?plot（t_serabs（xlfm））;
%?title（‘LFM信號時域包絡‘）;
%?axis（[-T/2?T/2?0?2]）;
%?figure;
%?plot（t_serphase（xlfm））;
%?title（‘LFM信號相位‘）;
%?
%?figure;
%?plot（t_serreal（xlfm））;
%?title（‘LFM信號實部‘）;
%?
%?figure;
%?plot（t_serimag（xlfm））;
%?title（‘LFM信號虛部‘）;

fft_lfm=?fft（xlfm）;
fft_lfm=fftshift（fft_lfm）;

%?figure;
%?plot（w_ser20*log（abs（fft_lfm）））;
%?title（‘LFM信號頻譜幅頻特性‘）;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘幅度/dB‘）;
%?
%?figure;
%?plot（w_serphase（fft_lfm））;
%?title（‘LFM信號頻譜相頻特性‘）;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘相位‘）;
%?
%?
%?figure;
%?plot（w_ser20*log（real（fft_lfm）））;
%?title（‘LFM信號頻譜實部‘）;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘幅度/dB‘）;
%?
%?
%?figure;
%?plot（w_ser20*log（imag（fft_lfm）））;
%?title（‘LFM信號頻譜虛部‘）;
%?xlabel（‘頻率/Hz‘）;
%?ylabel（‘幅度/dB‘）;


echo?=zeros（1tlen）;
echo（1tstart:tstart+N-1）=echo（1tstart:tstart+N-1）+xlfm;
echo（1tstart+deltaN:tstart+deltaN+N-1）=echo（1tstart+deltaN:tstart+deltaN+N-1）+0.5*xlfm;
%?figure
%?plot（techo_serabs（ech