主要内容

icwt

一维连续小波逆变换

描述

例子

xrec= icwt (wt)倒置了连续小波变换(CWT)系数矩阵wt使用默认值。icwt假设您获得类使用使用默认分析莫尔斯(60)小波。这小波的对称性3和带宽的60。icwt还假定CWT使用默认的尺度。如果wt是一个二维矩阵,icwt假设CWT获得从一个实值信号。如果wt是一个三维矩阵,icwt假设CWT得到复值信号。对于一个三维矩阵,的第一个页面wt是积极的CWT(逆时针)组件和第二页的wt是消极的(顺时针)组件。页面代表的分析和anti-analytic部分类,分别。

例子

xrec= icwt (wt,wname)使用小波分析wname反变换。指定的小波必须使用的小波一样。有效的选择wname“莫尔斯”,“爱”,“撞”指定莫尔斯,Morlet,凹凸小波,分别。

例子

xrec= icwt (wt,f,freqrange)反转CWT /在指定的频率范围freqrangefscale-to-frequency转换获得吗

例子

xrec= icwt (wt,,periodrange)反转CWT在指定的时间范围periodrange是一个数组的时间获得的输入时间。的输出获得使用持续时间输入。必须增加,包含在周期范围内

例子

xrec= icwt (___,名称,值)返回逆CWT由一个或多个指定附加选项名称,值对参数。

例子

全部折叠

获得语音样本的类和转化的CWT使用默认分析莫尔斯波。

负载mtlb;wt = cwt (mtlb);xrec = icwt (wt);

获得连续小波变换的语音样本和重构示例使用凹凸小波而不是默认的莫尔斯波。

负载mtlbdt = 1 / f;t = 0: dt元素个数(mtlb) * dt-dt;

获取类。

bumpmtlb = cwt (mtlb Fs,“撞”);

获得反变换。信号的意思是添加到输出。

xrec = icwt (bumpmtlb,“撞”,“SignalMean”意思是(mtlb));

绘制原始和重建信号。

情节(t, mtlb)包含(“秒”)ylabel (“振幅”)举行情节(t, xrec“r”)传说(“原始”,“重建”)

玩和比较原始和重建信号。

p = audioplayer (mtlb Fs);播放(p)暂停(2)px = audioplayer (xrec Fs);玩(像素)

重建一个frequency-localized近似神户地震数据中提取信息的类。提取的信息对应的频率范围(0.030 - 0.070)赫兹。

负载科比;

获取类。然后,获得反变换和添加信号均值回重建数据。CWT不保留信号的意思。

(wt, f) = cwt(科比,1);xrec = icwt (wt f (0.030 - 0.070),“SignalMean”意思是(科比));

绘制原始和重建数据。

次要情节(211)情节(科比);网格标题(“原始数据”);ylabel (“振幅”)次要情节(212)情节(xrec);网格标题(带通滤波重建[0.030 - 0.070]赫兹的);包含(“频率”);ylabel (“振幅”);

图包含2轴。轴与标题1原始数据包含一个类型的对象。轴2标题带通滤波重建[0.030 - 0.070]赫兹包含一个类型的对象。

利用连续小波逆变换重建一个近似厄尔尼诺现象数据基于2 - 8年时间。

加载厄尔尼诺现象数据并获得其类。每月数据采样。获取时间,指定采样间隔为一年的1/12。

负载ninoairdata;[cfs,时期]= cwt(尼诺,年(1/12));

获得的逆CWT 2到8年的周期。

xrec = icwt (cfs,期间,(2)年(8)]);

情节的CWT重建数据和原始数据的类进行比较。

类(尼诺,年(1/12));标题(“原始数据”);

图包含一个轴。坐标轴标题原始数据包含3类型的对象形象,线,区域。

图;类(xrec年(1/12));标题(基于2 - 8年时间近似的);

图包含一个轴。坐标轴标题近似基于2 - 8年时间包含3图像类型的对象,线,区域。

比较原始数据和重建数据。

图;次要情节(211)情节(datayear,尼诺);网格;甘氨胆酸ax =;斧子。XTickLabel =;轴;标题(“原始数据”);次要情节(212)情节(datayear xrec);网格;轴;包含(“年”);标题(“厄尔尼诺数据- 2 - 8年时间”);

图包含2轴。轴与标题1原始数据包含一个类型的对象。与标题厄尔尼诺现象数据轴2——2 - 8年时间包含一个类型的对象。

添加一个趋势的连续小波变换复值数据集和重建。

获得NPG2006数据集的类。

负载npg2006.matwt = cwt (npg2006.cx);

创建一个来自数据的时变趋势。

趋势= smoothdata (npg2006.cx,“movmean”,100);

获得反变换和增加的趋势。画出原始数据和重建数据。

xrec = icwt (wt,“SignalMean”、趋势);情节([真正的(xrec)”(npg2006.cx)])网格传奇(“趋势”,“原始”)

图包含一个轴。轴包含2线类型的对象。这些对象代表趋势,原创。

图绘制([图像放大(xrec)的图像放大(npg2006.cx)])网格传奇(“趋势”,“原始”)

图包含一个轴。轴包含2线类型的对象。这些对象代表趋势,原创。

输入参数

全部折叠

连续小波变换系数,指定为一个矩阵复杂的值。wt的输出是什么函数。

数据类型:
复数的支持:万博1manbetx是的

分析小波用于反变换,指定为“莫尔斯”,“爱”,或“撞”。这些特征向量指定分析莫尔斯,Morlet,和凹凸小波分别。指定的小波必须得到使用的相同类型的小波

默认的莫尔斯对称小波使用默认参数, γ 3和有一个默认的时间带宽60。

波变换频率,指定为一个向量。频率向量中的元素的数量必须等于输入变换系数矩阵的行数,wt。如果您指定f,您还必须指定freqrange

频率范围为返回逆连续小波变换值,指定为一个双元素向量或2×2矩阵。如果wt是一个二维矩阵,freqrange必须是一个双元素向量。如果wt是一个三维矩阵,freqrange可以是一个双元素向量或一个2×2的矩阵。如果wt是一个三维矩阵和freqrange执行是一个矢量,反演在同一频率范围的积极(分析)和消极(anti-analytic)的组件wt。如果freqrange是一个2×2矩阵第一行包含积极的频率范围的一部分吗wt(第一页)和第二行包含负面的一部分的频率范围wt(第二页)。一个向量的元素freqrange必须严格增加和包含在范围的频率向量f。为一个矩阵,每一行freqrange必须严格增加,包含的范围ff在CWT scale-to-frequency转换获得。对于复值信号的反演,可以指定的一行freqrange作为一个零的向量。如果第一行freqrange0是一个矢量,只有负(anti-analytic部分)是用于反演。如果您指定freqrange,您还必须指定f

例如[0 0;1/10的1/4)颠倒的负频率范围(顺时针)组件(1/10 1/4)。积极的(逆时针)组件是执行反转之前首先设置为0。同样的,[1/10的1/4;0 0]反转CWT通过选择频率范围(1/10 1/4)从正(逆时针)组件和设置负所有零组件。

时间对应的行波变换系数矩阵wt,指定为一个向量。的输出是当CWT获得使用持续时间输入。

期范围的返回逆连续小波变换值,指定为一个双元素向量或2×2矩阵。如果wt是一个二维矩阵,periodrange必须是一个双元素向量的持续时间。如果wt是一个三维矩阵,periodrange可以双元素向量的时间或时间的2×2矩阵。如果periodrange是一个向量的持续时间和wt执行3-Dmatrix,反演在同一频率范围的积极的(分析)和消极(anti-analytic)的组件wt。如果periodrange是持续时间的2×2矩阵第一行包含的时间范围的积极的部分wt(第一页)和第二行包含负面的周期范围的一部分wt(第二页)。一个向量的元素periodrange必须严格增加和包含在范围的向量。的元素periodrange必须有相同的单位。为一个矩阵,每一行periodrange必须严格增加和包含在范围的矢量p .复值信号的反演,可以指定一行的periodrange作为一个向量的零时间。如果第一行periodrange是一个向量的零时间,只有负面(anti-analytic部分)是用于反演。如果您指定periodrange,您还必须指定

例如(秒秒(0)(0);秒(1/10)秒(1/4)颠倒的负面(顺时针)组件在此期间范围内(秒(1/10)秒(1/4))。积极的(逆时针)组件是执行反转之前首先设置为0。同样的,(秒(1/10)秒(1/4);秒秒(0)(0))反转CWT通过选择时间范围(1/10 1/4)从正(逆时针)组件和设置负所有零组件。

名称-值对的观点

指定可选的逗号分隔条名称,值参数。的名字参数名称和吗价值相应的价值。的名字必须出现在引号。您可以指定几个名称和值对参数在任何顺序Name1, Value1,…,的家

例子:“TimeBandwidth”, 45设置时间带宽45。

莫尔斯波的时间带宽,指定为逗号分隔两人组成的“TimeBandwidth”和一个标量大于3小于或等于120。必须在指定的时间带宽与带宽使用的时间相同。莫尔斯的对称小波被认为是3。

你不能指定“TimeBandwidth”“WaveletParameters”。如果您指定“TimeBandwidth”,你就不能指定“WaveletParameters”。指定两个对称和时间带宽,使用“WaveletParameters”代替。

莫尔斯波的对称性和时间带宽,指定为逗号分隔组成的“WaveletParameters”双元素向量和一个标量。向量的第一个元素是对称的, γ ,第二个元素是时间带宽。指定的小波参数必须与使用的参数相同

你不能指定“WaveletParameters”“TimeBandwidth”。如果您指定“WaveletParameters”,你就不能指定“TimeBandwidth”。指定的时间对称带宽和使用默认值3,使用“TimeBandwidth”代替。

信号的意思是添加的icwt输出指定为逗号分隔组成的“SignalMean”和一个标量或矢量。如果信号的意思是一个向量,它必须相同长度作为小波系数矩阵的列的大小。如果wt是一个二维矩阵,意味着必须是实值标量或矢量信号。如果wt是一个三维矩阵,必须复值信号的意思是标量或矢量。因为不保留信号的意思是,icwt不包含信号的意思是默认。添加一个非零频率信号的意思是——或period-limited重建将零频组件添加到重建。

缩放系数用于逆类,指定为一个真正的或复数向量,获得作为一个可选的输出

比例系数的输出只支持莫尔斯和小波分析Morlet小波。万博1manbetx的长度ScalingCoefficients等于列的大小wt

你不能指定“SignalMean”“ScalingCoefficients”名称-值对。

数量每八度的声音,指定为逗号分隔组成的“VoicesPerOctave”甚至一个整数从4至48。CWT尺度离散使用指定的数量每八度的声音。每个八度的声音数必须相同的数量每八度用来获得CWT的声音。

你不能指定的数量每八度的声音如果指定的频率,f或持续时间,

输出参数

全部折叠

逆一维连续小波变换,作为一个真正的返回或复数的行向量。

数据类型:

更多关于

全部折叠

连续小波逆变换

icwt计算逆CWT基于离散由于Morlet版本的单积分公式[5]。简要描述的理论基础为单一的积分公式,明白了连续小波逆变换。额外的理论信息,参见2.4节[6]。提出了离散版本的积分[7]

引用

[1]莉莉,j . M。,S. C. Olhede. "Generalized Morse Wavelets as a Superfamily of Analytic Wavelets."IEEE信号处理60,不。(2012年11月11日):6036 - 41。https://doi.org/10.1109/TSP.2012.2210890。

[2]莉莉,J.M.,S.C. Olhede. "Higher-Order Properties of Analytic Wavelets."IEEE信号处理57岁的没有。1(2009年1月):146 - 60。https://doi.org/10.1109/TSP.2008.2007607。

[3]莉莉,j . M。为Matlab jLab:数据分析包1.6.2版本。2016年。http://www.jmlilly.net/jmlsoft.html。

[4]莉莉,j . M。,J.-C。Gascard。“时变椭圆信号的小波脊诊断与应用程序到一个海洋艾迪。”在地球物理非线性过程13日,没有。5(2006年9月14日):467 - 83。https://doi.org/10.5194/npg - 13 - 467 - 2006。

[5]Duval-Destin, M。,M. A. Muschietti, and B. Torresani. “Continuous Wavelet Decompositions, Multiresolution, and Contrast Analysis.”暹罗在数学分析》杂志上24日,没有。3(1993年5月):739 - 55。https://doi.org/10.1137/0524045。

[6]Daubechies,英格丽德。十个专题小波。61年应用数学CBMS-NSF地区会议系列。费城,宾夕法尼亚州:工业与应用数学学会,1992年。

[7]一如,克里斯托弗·吉尔伯特p .混合涂料。“小波分析的实用指南。”美国气象学会的公告79年,没有。1(1998年1月1日):61 - 78。https://doi.org/10.1175/1520 - 0477 (1998) 079 < 0061: APGTWA > 2.0.CO; 2。

[8]Holschneider, M。,Ph. Tchamitchian. “Pointwise Analysis of Riemann’s 'Nondifferentiable' Function.”发明Mathematicae105年,没有。1(1991年12月):157 - 75。https://doi.org/10.1007/BF01232261。

介绍了R2016b