matlab中wenergy函数,Python中一维小波的能量

我想知道在Python中是否有一个一维小波能量的实现，与Matlab的'[Ea，Ed]=wenergy(C，L)'相同。

我试着自己写一篇，但我不确定：

公式为：

其中Dj是细节向量，j=1,2，…，ld和N1是分解层的数据长度。在import json

import pywt

f=open('DataFile.txt','r')

D=json.load(f)

f.close()

#create the wavelet function

db1 = pywt.Wavelet('db13')

#calculate the number of necessary decompositions

NbrDecomp= pywt.dwt_max_level(len(D), db1)+1

#Initialize an empty list to receive the Detail and Approximation

Vector = [None] * NbrDecomp

#we use the Wavelet decomposition in the pywt module

Vector = pywt.wavedec(D, db1)

#Now Vector = [Approxiamtion N, Details N, Details N-1,.....]

#Where N would be the number of decompositions

根据定义，k级能量为：

^{pr2}$

实施是否正确？在