主要内容

墓场

小波同步压缩变换

描述

例子

风场=墓场(x返回小波同步压缩变换,风场,用于检查时频平面上的数据。与连续小波变换(CWT)相比,同步压缩变换减少了能量涂抹。输入,x,必须是一个至少有四个样本的一维实值信号。墓场利用解析Morlet小波计算同步压缩变换。

墓场函数对分析小波进行归一化以保留L1范数。有关更多信息,请参见算法

例子

风场f= wsst(x返回一个频率向量,f,以每个样本的周期为单位。频率对应于风场

例子

___= wsst(xfs使用指定的采样频率计算同步压缩变换,fs,单位为Hz,计算同步压缩变换。如果您指定f输出,墓场返回以Hz为单位的频率。您可以使用任何先前的输出值组合。

例子

___= wsst(xts使用一个持续时间ts用一个正的标量输入,作为采样间隔。持续时间的单位为年、天、小时、分钟和秒。如果你指定tsf输出,墓场返回频率f单位时间的周期,其中时间单位来源于指定的持续时间。

例子

___= wsst(___wav所指定的解析小波wav计算同步压缩变换。有效值为“爱”而且“撞”,分别指定解析Morlet小波和bump小波。

墓场(___在没有输出参数的情况下,将同步压缩变换绘制为时间和频率的函数。如果不指定采样频率,fs,或区间,ts,同步压缩变换以每个样本的周期为单位绘制。如果指定采样频率,则同步压缩变换以Hz为单位绘制。如果使用持续时间指定采样间隔,则图形以每单位时间的周期为单位。时间单位来源于持续时间。

___= wsst(___名称,值返回带有一个或多个指定的附加选项的同步压缩转换名称,值对参数。

例子

全部折叠

使用默认值获取语音样本的小波同步压缩变换。

负载mtlb;SST = wsst(mtlb);

得到语音信号的小波同步压缩变换,并对原始信号和重构信号进行比较。

加载语音信号,得到语音信号的同步压缩变换。

负载mtlbsoundsc(mtlb,Fs) dt = 1/Fs;T = 0:dt:数字(mtlb)*dt-dt;[sst,f] = wst (mtlb,Fs);

画出同步压缩变换。

pcolor (t、f、abs (sst)阴影插值函数包含(“秒”) ylabel (的频率(赫兹))标题(“Synchrosqueezed变换”

得到同步压缩逆变换,并播放重构后的语音信号。

Xrec = iwsst(sst);soundsc (xrec Fs)

得到并绘制二次啁啾的小波同步压缩变换。采样频率为1000赫兹。

负载quadchirp;[sst,f] = wsst(quadchirp,1000);HP = pcolor(tquad,f,abs(sst));惠普。EdgeColor =“没有”;标题(“小波同步变换”);包含(“时间”);ylabel (“赫兹”);

图中包含一个轴对象。标题为“小波同步压缩变换”的axis对象包含一个类型为surface的对象。

利用默认Morlet小波获得太阳黑子数据的小波同步压缩变换。指定采样间隔为一年。

负载太阳黑子墓场(太阳黑子(:,2)年(1))

图中包含一个轴对象。标题为“小波同步压缩变换”的axis对象包含一个类型为surface的对象。

利用凹凸小波对太阳黑子数据进行小波同步压缩变换。指定抽样间隔为1,即每年一次抽样。

负载太阳黑子墓场(太阳黑子(:,2)年(1)“撞”

图中包含一个轴对象。标题为“小波同步压缩变换”的axis对象包含一个类型为surface的对象。

输入参数

全部折叠

输入信号,指定为行或列向量。x必须是一个至少有四个样本的一维实值信号。

采样频率,指定为正标量。

采样间隔,也称采样周期,指定为a持续时间输入是正标量。有效期限为小时,分钟.不能使用日历持续时间(caldayscalweekscalmonthscalquarters,或calyears).你不能同时指定两者ts而且fs

例子:SST = wsst(x,hours(12))

分析小波用于计算同步压缩变换,指定为以下之一:

  • “爱”解析Morlet小波

  • “撞”-凹凸小波

名称-值参数

指定可选参数对为Name1 = Value1,…,以=家,在那里的名字参数名称和价值对应的值。名称-值参数必须出现在其他参数之后,但对的顺序无关紧要。

在R2021a之前,使用逗号分隔每个名称和值,并将其括起来的名字在报价。

例子:“VoicesPerOctave”,26岁

同步压缩转换中每八度使用的声音数,指定为由逗号分隔的对组成“VoicesPerOctave”一个从10到48的偶数。每个八度音阶的声音数与八度音阶数的乘积就是音阶数。八度音阶的数目取决于输入的大小x,是地板(log2(元素个数(x))) 1

选项,以对称地扩展输入信号,指定为由逗号分隔的对组成“ExtendSignal”,要么真正的.对称扩展信号可以缓解边界效应。如果你指定,则信号不扩展。如果你指定真正的,则信号被扩展。

输出参数

全部折叠

同步压缩变换,作为矩阵返回。默认情况下,同步压缩转换使用地板(log2(元素个数(x))) 1八度,每八度32个声音,以及解析Morlet小波。风场是一个Na——- - - - - -N矩阵Na是天平的多少,和N样品数量在吗x.默认的刻度数为32 *(地板(log2(元素个数(x))) 1)

同步压缩变换的频率,作为矢量返回。频率对应于风场.如果您没有指定fsts时,频率为每个样本的周期。如果你指定fs,频率单位为Hz。如果你指定ts,频率单位为周期/单位时间。频率向量的长度与的个数相同风场行。如果你指定ts作为抽样区间,ts用于计算刻度到频率的转换f

算法

墓场函数对分析小波进行归一化以保留L1范数。同样的说法是墓场不会将小波带通滤波器的傅里叶变换乘以尺度的平方根。乘以刻度的平方根将不平等地加权不同的带通贡献。

通过L1归一化,如果你在不同尺度的数据中有等幅振荡分量,它们在CWT中就会有相等的幅值。的函数也使用L1归一化。有关更多信息,请参见CWT的L1范数

参考文献

[1] Daubechies,英格丽,卢剑锋和吴昊婷。同步压缩小波变换:经验模态分解类工具应用与计算谐波分析30日,没有。2(2011年3月):243-61。https://doi.org/10.1016/j.acha.2010.08.002。

[2]塔库尔,高拉夫,Eugene Brevdo, Neven S. fuikkar和Hau-Tieng Wu。时变光谱分析的同步压缩算法:鲁棒性和新的古气候应用。信号处理93年,没有。5(2013年5月):1079-94。https://doi.org/10.1016/j.sigpro.2012.11.029。

版本历史

在R2016a中引入