ifsst

加载一个语音信号在采样 $$F_{\mathrm{年代}}=7418Hz$$。文件包含一个记录一个女声说“MATLAB®。”这个词Compute the Fourier synchrosqueezed transform of the signal.

负载mtlb%听,声音类型(mtlb Fs)海温,[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 Hz 2秒。

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

在第一第二,信号由一个400 Hz正弦信号和凹二次唧唧喳喳。指定一个对称区间中点的唧唧声,开始和结束在250赫兹的频率,和至少达到150 Hz。

t1 = t (1: nSamp / 2);x11 =罪(2 *π* 400 * t1);x12 =唧唧声(t1-t1 (nSamp / 4), 150年,nSamp / Fs, 1750年,“二次”);x1 = x11 + x12;

其余的由两个线性啁啾信号降低频率。一个唧唧喳喳的初始频率250赫兹,减少到100赫兹。其他唧唧喳喳的初始频率400赫兹,减少到250赫兹。

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

计算信号的傅里叶synchrosqueezed变换。指定一个256 -样本Kaiser窗形状参数β= 100。使用的绘图功能fsst显示的结果。

sig = (x1, x2);风=凯瑟(256120);[sigtr,功能处理量,竞技场队伍]= fsst(团体、Fs、风);fsst(团体、Fs、风、“桠溪”)

反变换重构功能。画出原始和反向信号以及它们之间的区别。

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

diffnorm =规范(x-sig)

diffnorm = 3.9026 e-13

产生一个信号,包括两个啾啾。信号采样3千赫一秒钟。第一个唧唧喳喳的初始频率400赫兹和达到800 Hz的抽样。第二啁啾起价500 Hz,最后达到1000 Hz。第二个唧唧喳喳的短促声波的两倍。

fs = 3000;t = 0:1 / fs: 1 - 1 / f;x1 =唧唧喳喳(800 t, 400 t(结束),);x2 = 2 *唧唧喳喳(1000 t, 500 t(结束),);

计算信号的傅里叶synchrosqueezed变换和阴谋。显示时间x设在和频率y设在。

海温,[f] = fsst (x1 + x2, fs);fsst (x1 + x2, fs,“桠溪”)

提取脊对应信号的组成部分,即啁啾与更大的振幅。使用脊来重构信号。

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

图的谱图的高能组件。将组件划分为256 -示例部分和指定一个重叠的255个样本。用DFT 512点和一个矩形窗口。

光谱图(xrec rectwin(256), 255512年,fs,“桠溪”)

提取第二个唧唧喳喳,指定tfridge搜索两个山脊。第二列的输出信号的低能的组件。

[~,iridge] = tfridge (sst, f,“NumRidges”2);[],xrec = ifsst (sst iridge (:, 2));光谱图(xrec rectwin(256), 255512年,fs,“桠溪”)