Управление коэффициентами в дискретном вейвлет-преобразовании пакета в R с помощью wavelim
В настоящее время я работаю с вейвлет-преобразованием в R с пакетом 'waveslim'. Я использую дискретный пакетный вейвлет-преобразование dwpt для восстановления определенной частотной части временного ряда. Там следует выбрать соответствующие коэффициенты (обычно первые n коэффициентов для самых низких частот). Если я выберу первые n коэффициентов, восстановленный сигнал изменится для одного и того же временного ряда после каждого прогона, и график будет неубедительным - он, похоже, нечастая часть этого сигнала. Вот мой фрагмент кода:
WaveletFilter_Waveslim = function(timeSeries, level = 2){
WaveletDecompositionVector <- dwpt(timeSeries, "la8", level)
print(WaveletDecompositionVector)
timeSeries.basis <- basis(WaveletDecompositionVector, c("w2.2", "w2.3"))
for(i in 1:length(WaveletDecompositionVector))
WaveletDecompositionVector[[i]] <- timeSeries.basis[i] * WaveletDecompositionVector[[i]]
filtered <- idwpt(WaveletDecompositionVector, timeSeries.basis)
dev.new()
par(mfrow=c(2,1), mar=c(5-1,4,4-1,2))
plot.ts(x, xlab="", ylab="", main="Original Series")
plot.ts(filtered, xlab="", ylab="", main="Reconstructed Series")
}Преобразование дискретного вейвлет-пакета с использованием алгоритма Маллата.