我有下面的数据和代码来创建 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]