连续小波变换的频率和时间局部重构

重建神户地震资料的频域近似。从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范围内的频率数据。

次要情节(2,1,1)情节(科比)网格在标题(“原始数据”) subplot(2,1,2) plot(xrec) grid在标题(带通滤波重建[0.030 0.070]Hz)

在CWT中，还可以使用时间周期，而不是频率。加载厄尔尼诺数据，得到其CWT，指定时间周期，以年为单位。

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

得到第2年到第8年的CWT的倒数。

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

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

类图(xrec年(1/12))

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

图次plot(2,1,1) plot(nino)网格在标题(“原始数据”) subplot(2,1,2) plot(xrec) grid在标题(“厄尔尼诺数据-第2-8年”)

小波工具箱文档