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

打开实时脚本

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

负载科比

获取数据的CWT。

[wt，f]=cwt（神户，1）；

重建地震数据，将信号平均值添加回转换后的数据中。

xrec=icwt（重量，f，[0.030 0.070]，“信号平均值”，平均（神户）；

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

子地块（2,1,1）地块（神户）网格在…上头衔(“原始数据”)子地块（2,1,2）绘图（xrec）网格在…上头衔(‘带通滤波重建[0.030 0.070]Hz’)

图中包含2个轴对象。具有标题原始数据的轴对象1包含line类型的对象。具有标题带通滤波重建[0.030 0.070]Hz的轴对象2包含line类型的对象。

你也可以在CWT中使用时间段，而不是频率。加载厄尔尼诺数据并获取其CWT，以年为单位指定时间段。

负载花鸟[cfs，期间]=cwt（尼诺，年（1/12））；

获得第2年至第8年的反向CWT。

xrec=icwt（cfs，期限，[年（2）年（8）]）；

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

图cwt（xrec，年份（1/12））

图中包含一个axes对象。标题为Magnitude Scalegram的axes对象包含3个类型为image、line、area的对象。

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

图子地块（2,1,1）图（nino）网格在…上头衔(“原始数据”)子地块（2,1,2）绘图（xrec）网格在…上头衔(“厄尔尼诺数据——2-8年”)

图中包含2个轴对象。标题为原始数据的轴对象1包含线条类型的对象。标题为厄尔尼诺数据的轴对象2-2-8年包含线条类型的对象。

小波工具箱文档