代码参数说明：Sf1为书中公式3.19，Sf2为时域信号快速傅里叶表达式，两种频谱表达式所做出的图可看出其区别
代码如下：
clc
clear all
close all
%参数设置
TBP = 720; %时间带宽积
T = 10e-6; %脉冲持续时间
%参数计算
B = TBP/T; %信号带宽
K = B/T; %线性调频频率
alpha = 1.25; %过采样率
F = alpha*B; %采样率
N = 2*ceil(F*T/2); %采样点数
dt = T/N; %采样时间间隔
df = F/N; %采样频率间隔
%变量设置
t = -T/2:dt:T/2-dt; %时间变量
f = -F/2:df:F/2-df; %频率变量
%信号表达
st = exp(1j*pi*K*t.^2); %chirp信号复数表达式
Sf1 = exp(-1j*pi*f.^2/K); %chirp信号频谱表达式
Sf2 = fftshift(fft(fftshift(st))); %chirp信号快速傅里叶变换
%绘图
%Sf1
figure
subplot(221),plot(f*1e-6,real(Sf1),'r')
axis([-10 10,-1.2 1.2])
subplot(222),plot(f*1e-6,abs(Sf1))
axis([-10 10,-1.2 1.2])
subplot(223),plot(f*1e-6,imag(Sf1),'r')
axis([-10 10,-1.2 1.2])
subplot(224),plot(f*1e-6,unwrap(angle(Sf1)))
axis([-28 28,0 400])
%Sf2
figure
subplot(221),plot(f*1e-6,real(Sf2),'r')
axis([-10 10,-40 40])
subplot(222),plot(f*1e-6,abs(Sf2))
axis([-50 50,-0 40])
subplot(223),plot(f*1e-6,imag(Sf2),'r')
axis([-10 10,-40 40])
subplot(224),plot(f*1e-6,unwrap(angle(Sf2)))
axis([-50 50,0 900])
1. N = 2*ceil(F*T/2);
ceil: 向正无穷舍入，保证N是偶数
2. plot(f*1e-6,unwrap(angle(Sf1)))
Q = unwrap(P) 展开向量 P 中的弧度相位角。每当连续相位角之间的跳跃大于或等于 π 弧度时，unwrap 就会通过增加 ±2π 的整数倍来平移相位角，直到跳跃小于 π。如果 P 是矩阵，unwrap 将按列运算。如果 P 是多维数组，unwrap 将对大小大于 1 的第一个维度进行运算。
举个例子：一般在我们计算一个系统相频特性时，就要用到反正切函数提取相位，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四象限的角度为0～-pi。但实际得到的结果会发生相位跳变，跳变幅度为2pi，这就叫相位的卷绕。unwrap函数的作用就是解卷绕，使相位在pi处不发生跳变，从而反应出真实的相位变化。
