代码复现如下：

clc
clear
close all% 参数设置
TBP = 100;                              % 时间带宽积
T = 7.2e-6;                             % 脉冲持续时间
t_0 = 1e-6;                             % 脉冲回波时延% 参数计算
B = TBP/T;                              % 信号带宽
K = B/T;                                % 线性调频频率
alpha_os = 1.25;                        % 过采样率，使用较高的过采样率是为了提高采样频率
F = alpha_os*B;                         % 采样频率
N = 2*ceil(F*T/2);                      % 采样点数
dt = T/N;                               % 采样时间间隔
df = F/N;                               % 采样频率间隔% 变量设置
t = -T/2+t_0:dt:T/2+t_0-dt;             % 时间变量
f = -F/2:df:F/2-df;                     % 频率变量%信号表达
st = exp(1j*pi*K*t.^2);                 %chirp信号复数表达式
srt = exp(1j*pi*K*(t-t_0).^2);          %chirp信号回波表达式
Srf = fft(srt);                         %chirp信号回波频域表达式
Hf = exp(1j*pi*f.^2/K);                 %匹配滤波器频域表达式，见公式3.35
S_out_f = Srf.*Hf;                      %匹配滤波后的信号频域表达式%绘图
set(figure,'position',[500,500,600,150])
subplot(121),plot(f*1e-6,real(S_out_f)),axis([-10 10,-15 15])
subplot(122),plot(f*1e-6,imag(S_out_f)),axis([-10 10,-15 15])