加载在取样的语音信号 $$F_{\mathrm{年代}}=7418Hz$$。该文件包含一段女性声音的录音，她说“MATLAB®”。计算信号的傅里叶同步压缩变换。

负载MTLB％听，声音类型（MTLB，FS）[SST中，f] = FSST（MTLB，FS）;

反转变换以重构该信号。绘制原和重构信号，以及它们之间的区别。

xrec = ifsst (sst);t =(0:长度(mtlb) 1) / Fs;xrec mtlb情节(t, t, t, mtlb-xrec)包含('时间（s）')传说(“原始”,“重建”,“差异”)

通过计算检查重建的准确性 $${\ell }_{\infty }$$原信号和逆变换之间的差的范数。

LINF =范数（ABS（MTLB-xrec），天道酬勤）

Linf = 1.9762 e-14

%听，输入声音(mtlb-xrec,Fs)

产生在1024赫兹采样2秒钟的信号。

nSamp = 2048;Fs = 1024;t = (0: nSamp-1) / Fs;

在第一第二，该信号包括一个400赫兹正弦波和凹形二次啁啾的。指定一个线性调频脉冲是对称的，在250赫兹的频率的间隔中点，开始和结束，并达到最小150Hz的。

T1 = T（1：nSamp / 2）;X11 = SIN（2 * PI * 400 * T1）;X12 =啁啾（T1-T1（nSamp / 4），150，nSamp / Fs的，1750年，“二次”）;X1 = X11 + X12;

该信号的其余部分由减小的频率的两个线性调频脉冲的。一个线性调频脉冲具有250赫兹其降低至100Hz的初始频率。另一个线性调频脉冲具有400赫兹其降低至250Hz的初始频率。

t2 = t (nSamp / 2 + 1: nSamp);x21 =唧唧声(t2,400 nSamp / Fs 100);将=唧唧声(t2,550 nSamp / Fs 250);x2 = x21 +将;

计算信号的傅里叶同步压缩变换。使用形状参数指定一个256个样本的Kaiser窗口β= 100。使用的绘制功能FSST来显示结果。

SIG = [X1; X2];风=凯泽（256120）;[sigtr，FTR，TTR] = FSST（SIG，FS，风）;FSST（SIG，FS，风，“Y轴”)

反变换重建的功能。绘制原和反相信号以及它们之间的差异。

x = ifsst (sigtr、风);团体情节(t, t, x, t, x-sig)传说(“原始”,“重建”,“差异”)

diffnorm =规范（X-SIG）

diffnorm = 3.9026e-13

产生由两个啁啾组成的信号。信号以3千赫采样一秒。第一个啁啾的初始频率为400hz，在采样结束时达到800hz。第二个唧唧声从500hz开始，到最后达到1000hz。第二个啁啾的振幅是第一个啁啾的两倍。

FS = 3000;t = 0时：1 / FS：1-1 / FS;X1 =啁啾（吨，400，T（端部），800）;X2 = 2 *啁啾（吨，500，T（端部），1000）;

计算和绘制信号的傅里叶同步压缩变换。显示器上显示时间x轴以及在所述频率y设在。

[SST中，f] = FSST（X1 + X2，FS）;FSST（X1 + X2，FS“Y轴”)

提取信号的高能量分量对应的脊线，即振幅较大的啁啾。利用脊线重建信号。

[~, iridge] = tfridge (sst, f);[], xrec = ifsst (sst iridge);

绘制高能量组分的谱图。将组件划分为256个样本部分，并指定255个样本的重叠。使用512 DFT点和一个矩形窗口。

光谱图(xrec rectwin(256), 255512年,fs,“Y轴”)

要提取第二个唧唧声，请指定它tfridge搜索两个脊。的输出的第二列是信号的较低能量部分。

[~, iridge] = tfridge (sst, f,“NumRidges”，2）;xrec = ifsst（SST，[]，iridge（：，2））;光谱图(xrec rectwin(256), 255512年,fs,“Y轴”)