小波包反向DWT无法重建系列？

Question

我正在使用wavelets包，并注意到我尝试

时

library("wavelets")
x <- rnorm(100)
y <- idwt(dwt(x))
plot(x, y)

重建y显然不等于原始x。

这是预期的吗？

对于某些情况，我正在尝试使用一系列系列的小波变换进行（正则化）逻辑回归。然后，我想将回归系数映射回原始时间序列空间，以查看歧视中使用的时间点。

但我似乎无法重建原始系列。我可能完全误解了事情，有人可以对此有所了解吗？

Answer 1

在帮助文件?dwt之后，您可以修改脚本，例如：

library(wavelets)
set.seed(42)
x <- rnorm(100)
y <- idwt(dwt(x, n.levels=3, boundary="reflection", fast=FALSE))
plot(x, y)
abline(0,1)