3

Morlet のウェーブレット変換から元の時系列を再構築しようとしています。私はR、パッケージRwave、関数cwtで作業しています。この関数の結果は、複素数値を含む n*m (n=期間、m=時間) の行列です。

信号を再構築するために、 Torrence & Compo classic textの式 (11) を使用しましたが、結果は元の信号とは何の関係もありません。私はウェーブレット変換の実部とスケールの間の分割に特に関心があります。このステップは結果を完全に歪めます。一方、すべてのスケールで実数部分を合計すると、結果は元の時系列に非常に似ていますが、値がわずかに広くなります (元の系列の範囲は ~ [-0.2, 0.5]、再構築された系列の範囲は~ [-0.4,0.7])。

元の時系列を再構築するための実用的な手順、式、またはアルゴリズムについて誰かが教えてくれるかどうか疑問に思っています。Torrence と Compo (1998)、Farge (1992) などの論文を読んだことがありますが、それらはすべて異なる公式で書かれていますが、実際に助けてくれる人は誰もいません。

4

1 に答える 1