1

我有下面的数据和代码来创建 DWT 小波,但我不知道如何反转它以使其与原始数据对齐。下面的代码可以创建系数,但后来我被困在反转它。

我不想使用诸如 pywavelets 之类的包,而是从头开始。

data = [56, 40, 8, 24, 48, 48, 40, 16]

def discreteHaarWaveletTransform(x):
    N = len(x)
    output = [0.0]*N

    length = N >> 1
    while True:
        for i in range(0,length):
            summ = x[i * 2] + x[i * 2 + 1]
            difference = x[i * 2] - x[i * 2 + 1]
            output[i] = summ/2
            output[length + i] = difference/2

        if length == 1:
            return output

        #Swap arrays to do next iteration
        x = output[:length << 1]
        length >>= 1

输出:[35.0, -3.0, 16.0, 10.0, 8.0, -8.0, 0.0, 12.0]

4

0 回答 0