R小波分解
install.packages("wavelets")
library(wavelets)
X1 <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16)
#wt combine them and compute DWT 小波分解2层
wt <- dwt(X1, n.levels=2, boundary="reflection", fast=FALSE)
#将第二层细节函数赋值为0
wt@W$W2 <- as.matrix(data.frame(rep(0,8)))
#还原序列
idwt(wt)