该结果的STFT显示效果，参数不同，显示的结果也不同。
此处只演示线性关系的效果。
clc;
clear;
fs = 16000;
t = 0:2*fs;
phase = 2*pi*(0.05*t+1000)/fs.*t;
x = sin(phase);
spectrogram(x,256,250,256,fs,'yaxis');
clear;
fs = 16000;
t = 0:2*fs;
phase = 2*pi*(0.05*t+1000)/fs.*t;
x = sin(phase);
spectrogram(x,256,250,256,fs,'yaxis');
当然也可以是其他曲线。
只要把握好瞬时频率与相位之间的关系即可。
参考文献：https://en.wikipedia.org/wiki/Chirp