脉冲压缩代码

%%雷达成像处理脉冲压缩的仿真
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