代码复现如下:
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])