连续小波变换的频域和时域重构

打开实时脚本

重建神户地震数据的频率局部近似。从CWT中提取频率范围为[0.030,0.070]Hz的信息。

负载科比

获取数据的CWT。

[wt,f] = cwt(神户，1);

重建地震数据，将信号均值重新加入转换后的数据中。

Xrec = icwt(wt，[]，f，[0.030 0.070]，“SignalMean”意思是(科比));

绘制并比较原始数据和频率范围为[0.030,0.070]Hz的数据。

Subplot (2,1,1) plot(神户)网格在标题(“原始数据”) subplot(2,1,2) plot(xrec)网格在标题(带通滤波重建[0.030 - 0.070]Hz）

图中包含2个轴对象。标题为原始数据的坐标轴对象1包含一个类型为line的对象。axis对象2的标题为Bandpass Filtered Reconstruction [0.030 0.070] Hz包含一个类型为line的对象。

对于CWT，您还可以使用时间段，而不是频率。加载El Nino数据并获得其CWT，指定以年为单位的时间段。

负载ninoairdata[cfs,period] = cwt(nino,years(1/12));

获得第2年到第8年的逆CWT。

Xrec = icwt(cfs，[]，period，[years(2) years(8)]);

绘制重建数据的CWT。注意在2到8年的周期带之外没有能量。

类图(xrec年(1/12))

ans = proplistener with properties:对象:{[1x1 Axes]}来源:{[1x1 matlab.graphics.internal.GraphicsMetaProperty]} EventName: 'PostSet'回调:@changeYTickLabels启用:1递归:0

图中包含一个轴对象。标题为Magnitude scalalogram的轴对象包含图像、线、区域类型的3个对象。

将2 - 8年的原始数据与重建数据进行比较。

图中子图(2,1,1)图(尼诺)格在标题(“原始数据”) subplot(2,1,2) plot(xrec)网格在标题(“厄尔尼诺数据- 2-8年”）

图中包含2个轴对象。标题为原始数据的坐标轴对象1包含一个类型为line的对象。标题为El Nino Data - Years 2-8的axis对象2包含一个类型为line的对象。