Вейвлеты Хаара - Частичная реконструкция сигнала

Я пытаюсь аппроксимировать временной ряд/вектор с помощью 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)

Я читал об усечении коэффициентов, но я не знаю, как это сделать.

Что, если мне нужно приближение, основанное, скажем, на двух коэффициентах?

Я буду признателен за любую помощь.


person Nelson    schedule 22.05.2016    source источник


Ответы (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

и это имеет смысл: без коэффициентов детализации на самом высоком уровне преобразование Хаара возвращает среднее значение каждой пары точек для обеих точек.

person alle_meije    schedule 31.05.2016