r - Haar ウェーブレット - 部分信号再構成

Question

dwt を使用して時系列/ベクトルを概算しようとしています。完全な dwt および idwt 再構成を適用できますが、(すべてではない) わずかな係数に基づいて元の信号を再構成/近似する方法がわかりません。

私はこれを試しました：

library(wavelets)
x<-c(7, 5, 5, 3, 3, 3, 4, 6)
w <- dwt(x, filter="haar",n.levels = 3)
rec<-idwt(w)

係数の切り捨てについて読みましたが、その方法がわかりません。

たとえば、2 つの係数に基づく近似値が必要な場合はどうすればよいでしょうか。

どんな助けにも感謝します。

Answer 1

3行目と4行目の間でできること

str    (w)     # see what w looks like
w@W$W1 []  = 0 # set the 1st level detail coefficients to 0 

それからあなたのrec意志は

[1] 6 6 4 4 3 3 5 5

これは理にかなっています。最も細かいレベルの詳細係数がない場合、Haar 変換は、両方のポイントについて、ポイントの各ペアの平均を返します。