Documentation

ICWT

Inverse continuous 1-D wavelet transform

Syntax

xrec = icwt(wt)
XREC= icwt(wt,wname)
XREC= icwt(wt,f,freqrange)
XREC= icwt(wt,period,periodrange)
XREC= icwt(___,名称,值)

描述

example

XREC= icwt(wt)反转连续小波变换(CWT)系数矩阵wtusing default values.ICWTassumes that you obtained the CWT usingcwt带有默认的摩尔斯小波。该小波的对称性为3,时间带宽为60。ICWTalso assumes that the CWT uses default scales.

example

XREC= icwt(wt,wname)使用分析小波wname倒转CWT。指定的小波必须与使用的小波相同cwt。有效的选项wnameare'morse','amor', 和'撞', which specify the Morse, Morlet, and bump wavelet, respectively.

example

XREC= icwt(wt,f,freqrange)inverts the CWT over the two-element frequency range specified infreqrangef是从比例到频率转换cwt。这frequency range must be increasing and contained inf

example

XREC= icwt(wt,时期,周期)将CWT倒置在两元素范围的周期范围内周期。这时期是个cwtoutput obtained using aduration输入。周期范围必须增加并包含时期

example

XREC= icwt(___,名称,价值)返回逆CWT,其中一个或多个指定的其他选项名称,价值配对参数。

例子

全部收缩

Obtain the CWT of a speech sample and invert the CWT using the default analytic Morse wavelet.

loadMTLB;wt = cwt(mtlb);XREC = ICWT(WT);

Obtain the continuous wavelet transform of a speech sample and reconstruct the sample using the bump wavelet instead of the default Morse wavelet.

loadMTLB;dt = 1/fs;t = 0:dt:numel(mtlb)*dt-dt;

获取并绘制CWT。

[bumpmtlb,f] = cwt(mtlb,Fs,'撞');p1 = pcolor(t,f,abs(bumpmtlb)); p1.EdgeColor ='none';xlabel('Seconds');ylabel('赫兹');

获得逆CWT。

xrec = icwt(bumpmtlb,'撞',“SignalMean”,平均值(mTLB));

Plot the original and reconstructed signals.

plot(t,mtlb); xlabel('Seconds');ylabel('振幅');抓住;情节(t,xrec,'r');传奇('Original','重建');

Play and compare the original and reconstructed signals.

p = audioplayer(mtlb,Fs); play(p); pause(2); px = audioplayer(xrec,Fs); play(px);

通过从CWT中提取信息,重建对科比地震数据的频率定位近似。提取的信息对应于[0.030 0.070] Hz范围内的频率。

loadkobe;

获得CWT。然后,获取逆CWT并将信号均值添加回重建数据中。CWT不能保留信号均值。

[wt,f] = cwt(kobe,1); xrec = icwt(wt,f,[0.030 0.070],“SignalMean”,平均(神户));

绘制原始数据和重建数据。

subplot(211) plot(kobe); grid标题(“原始数据”);ylabel('振幅')子图(212)图(XREC);网格标题('Bandpass Filtered Reconstruction [0.030 0.070] Hz');xlabel('频率');ylabel('振幅');

使用逆连续小波变换根据2到8年的时间重建对El Nino数据的近似值。

Load the El Nino data and obtain its CWT. The data is sampled monthly. To obtain the periods in years, specify the sampling interval as 1/12 of a year.

loadNinoairdata;[cfs,period] = cwt(nino,years(1/12));

Obtain the inverse CWT for periods of 2 to 8 years.

XREC= icwt(cfs,period,[years(2) years(8)]);

Plot the CWT of the reconstructed data and compare it to the CWT of the original data.

cwt(nino,years(1/12)); title(“原始数据”);数字;CWT(XREC,年(1/12));标题(“基于2  -  8年的近似值”);

Compare the original data with the reconstructed data in time.

数字;subplot(211) plot(datayear,nino); grid;ax = gca;ax.xticklabel ='';轴紧的;标题(“原始数据”);子图(212)图(Datayear,XREC);网格;轴紧的;xlabel(“年”);标题('El Nino Data - 2-8 Year Periods');

输入参数

全部收缩

Continuous wavelet transform coefficients, specified as a matrix of complex values.wt是个output from thecwtfunction.

数据类型:双倍的
Complex Number Support:Yes

分析小波用于倒转CWT,指定为'morse','amor', or'撞'。这些角色向量分别指定分析性摩尔斯科,莫雷特或凹凸小波。指定的小波必须是用于获得的小波cwt

默认摩尔斯小波使用默认对称参数, γ , that is 3 and has a default time bandwidth of 60.

CWT frequencies, specified as a vector. The number of elements in the frequency vector must equal to the number of rows in the input CWT coefficient matrix,wt。如果指定f, you must also specifyfreqrange

返回逆连续小波变换值的频率范围,指定为两元素向量。这freqrange值必须在频率值的范围内f。第一个和第二个要素freqrange分别定义频率范围的开始和结尾。第一个元素必须小于第二个元素。如果指定freqrange, you must also specifyf

Time periods corresponding to the rows of CWT coefficient matrixwt,指定为矢量。时期是输出cwt, when the CWT is obtained using aduration输入。

返回逆连续小波变换值的周期范围,指定为两元素向量。这周期values must be in the range of the period values in时期。第一个和第二个要素周期分别定义频率范围的开始和结尾。第一个元素必须小于第二个元素。如果指定周期, you must also specify时期

名称值对参数

Specify optional comma-separated pairs of名称,价值参数。Name是个argument name and价值是相应的值。Name必须出现在单引号中(' ')。您可以按任何顺序指定几个名称和值对参数NAME1,Value1,...,Namen,Valuen

例子:'TimeBandwidth',45将时间带宽设置为45。

全部收缩

Time bandwidth of the Morse wavelet, specified as a comma-separated pair consisting of'TimeBandwidth'标量大于3且小于或等于120。指定的时间带宽必须与在该中使用的时间带宽相同cwt。这symmetry of the Morse wavelet is assumed to be 3.

您不能同时指定the'TimeBandwidth''WaveletParameters'。如果指定'TimeBandwidth', you cannot specify'WaveletParameters'。To specify both the symmetry and time bandwidth, use'WaveletParameters'instead.

摩尔斯小波的对称和时间带宽,指定为逗号分隔对'WaveletParameters'和a two-element vector of scalars. The first element of the vector is the symmetry, γ ,第二个元素是时间带宽。指定的小波参数必须与cwt

您不能同时指定'WaveletParameters''TimeBandwidth'。如果指定'WaveletParameters', you cannot specify'TimeBandwidth'。To specify the time bandwidth and use the default symmetry value of 3, use'TimeBandwidth'instead.

信号是指添加到iwct输出,指定为逗号分隔对,由“SignalMean”和标量。因为cwt不保留信号的平均值,ICWTdoes not contain the signal mean by default.

Number of voices per octave, specified as the comma-separated pair consisting of“ VoicesPeroCtave”并且一个从4到48的整数。使用每个八度的指定数量的声音将CWT量表离散。每个八度的声音数必须与用于获得CWT的每个八度的声音数相同。如果指定频率,则不能指定每个八度的声音数,f, or duration,时期

Output Arguments

全部收缩

逆1-D连续小波变换,返回为行向量。

数据类型:双倍的

参考

[1] Lilly,J。M.和S. C. Olhede。“广义的摩尔斯小波是分析小波的超家族。”信号处理的IEEE交易。卷。60,第11页,2012年,第6036–6041页。

[2] Lilly,J。M.和S. C. Olhede。“分析波的高阶特性”。信号处理的IEEE交易。卷。57,第1期,2009年,第146-160页。

[3] Lilly, J. M.JLAB:MATLAB的数据分析软件包,版本1.6.2。2016. http://www.jmlilly.net/jmlsoft.html。

Introduced in R2016b

这个话题有帮助吗?