MATLAB代码：小波包分解与重构、能量谱

两个matlab函数：一个是小波包分解与重构；一个是能量谱。可以自己更改成一个程序，可以达到能量特征提取的目的。

function E = wavelet_energy_spectrum（ wptn ）
%% 计算每一层每一个节点的能量
%  wpt-wavelet packet tree
%  n-第n层能量
% 
% Author hubery_zhang
% Date  20170714

%%
% 求第n层第i个节点的系数
E（1:2^n ）=0;
for i=1:2^n 
E（i） = norm（wpcoef（wpt[ni-1]）2）; %求第i个节点的范数平方，其实也就是平方和
end
%求每个节点的概率
E_total=sum（E）; 
for i=1:2^n
p_node（i）= 100*E（i）/E_total;
end
% E = wenergy（wpt）; only get the last layer
figure;
x=1:2^n;
bar（xp_node）;
title（[‘第‘num2str（n）‘层‘]）;
axis（[0 2^n 0 100]）;
xlabel（‘结点‘）;
ylabel（‘能量百分比/%‘）;
for j=1:2^n
text（x（j）p_node（i）num2str（p_node（j）‘%0.2f‘）...
    ‘HorizontalAlignment‘‘center‘...
    ‘VerticalAlignment‘‘bottom‘）
end

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        725  2017-10-19 20:32  wavelet_energy_spectrum.m

     文件       1486  2017-10-19 20:33  wavelet_packetdecomposition_reconstruct.m

----------- ---------  ---------- -----  ----

                 2211                    2