主要内容

墓场

小波同步压缩变换

描述

例子

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

例子

风场f= wsst()x返回频率向量,f,单位为周期/样本。频率对应于的行风场

例子

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

例子

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

例子

___= wsst()___wav使用指定的解析小波wav来计算同步压缩变换。有效值是“爱”“撞”,分别指定了解析Morlet小波和碰撞小波。

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

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

例子

全部折叠

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

负载mtlb;SST = wsst(mtlb);

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

加载语音信号并获得其同步压缩变换。

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

绘制同步压缩变换。

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

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

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

得到并绘制了二次啁啾的小波同步压缩变换。啁啾以1000hz采样。

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

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

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

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

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

利用碰撞小波获得并绘制了太阳黑子数据的小波同步压缩变换。指定抽样间隔为1,每年一个样本。

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

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

输入参数

全部折叠

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

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

采样间隔,又称采样周期,指定为a持续时间有正的标量输入。有效期限为小时,分钟.不能使用日历持续时间(caldayscalweekscalmonthscalquarters,或calyears)。不能同时指定tsfs

例子: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

参考文献

[10]多贝西,卢杰,吴海涛。同步压缩小波变换:一种类似经验模态分解的工具。应用与计算谐波分析.第30卷第2期,2011年,第243-261页。

[10]吴海涛,吴海涛,吴海涛。时变光谱分析的同步压缩算法:鲁棒性和新的古气候应用。信号处理.第93卷第5期,2013年,第1079-1094页。

版本历史

在R2016a中引入