Question

我正在尝试近似时间序列/向量usind 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变换会返回两个点的每对点的平均值。