代码复现如下：

clc
close all
clear all%参数设置
B = 5.80e+6;        %信号带宽
T = 7.26e-6;        %脉冲持续时间
K = B/T;            %线性调频频率
alpha = 5;          %过采样率
F = alpha*B;        %采样频率
N = F*T;            %采样点数
dt = T/N;           %采样时间间隔%变量设置
t = -T/2:dt:T/2-dt;             %时间变量
%另一种变量设置方式
%t = linspace(-T/2,T/2-dt,N);%信号表达
st = exp(1j*pi*K*t.^2);          %chirp信号复数形式%其他参数
f = K*t;            %瞬时频率
phi = pi*K*t.^2;    %瞬时相位%绘图
figure
subplot(221),plot(t*1e+6,real(st),'r')
axis([-4 4,-1.2 1.2])
subplot(222),plot(t*1e+6,phi)
axis([-4 4,0 40])
subplot(223),plot(t*1e+6,imag(st),'r')
axis([-4 4,-1.2 1.2])
subplot(224),plot(t*1e+6,f*1e-6)
axis([-4 4,-3 3])