小波系数去噪可以用wavethresh包，里面有去噪的函数，不过这个包中的小波分解不能设置分解层数

我自己根据无偏风险估计阈值编了个函数，可以进行小波系数去噪
for (j in 1:n) {
  for (i in names(DWT[[j]]@W)[1:length(names(DWT[[j]]@W))]) {
    DWTi <- DWT[[j]]@W[[i]]
    ni <- length(DWTi)
    factor <- median(abs(DWTi))/0.6745
    xi <- DWT[[j]]@W[[i]]/factor
    sxi <- sort(abs(xi))^2
    s <- cumsum(sxi) + ((ni - 1):0) * sxi
    risk <- (ni - (2 * (1:ni)) + s)/ni
    surethresh <- sqrt(sxi[order(risk)[1]])
    ###软阈值
    DWT[[j]]@W[[i]]<-sign(DWTi)*(abs(DWTi)-factor*surethresh)*(abs(xi) > surethresh)
   
  }
}
你可以参考一下