Matlab小波包分解後如何求各頻帶信號的能量值？

Matlab小波包分解後如何求各頻帶信號的能量值？ [轉]

clear;
clc;
n=3;
wpname='db3';
% [b,a]=butter(8,[5/100 99/100]);
% load a1_1-1;
load k301_1-4;
% Data=filter(b,a,Data);
Data=Data-mean(Data);
Data=detrend(Data);
wpt1=wpdec(Data,n,wpname); %對數據進行小波包分解
for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n層第i個節點的係數
E(i)=norm(wpcoef(wpt1,[n,i-1]),2);%求第i個節點的範數平方，其實也就是平方和
end
% disp('每個節點的能量E(i)');
% E
% disp('小波包分解總能量E_total');
E_total=sum(E); %求總能量
for i=1:2^n
pfir(i)= E(i)/E_total;%求每個節點的概率
end
pfir
% load a1_1-2;
load k302_1-4;
% Data=filter(b,a,Data);
Data=Data-mean(Data);
Data=detrend(Data);
wpt2=wpdec(Data,n,wpname); %對數據進行小波包分解
for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n層第i個節點的係數
E(i)=norm(wpcoef(wpt2,[n,i-1]),2);%求第i個節點的範數平方，其實也就是平方和
end
E_total=sum(E); %求總能量
for i=1:2^n
psec(i)= E(i)/E_total;%求每個節點的概率
end
psec
% load a1_1-4;
load k303_1-4;
% Data=filter(b,a,Data);
Data=Data-mean(Data);
Data=detrend(Data);
wpt3=wpdec(Data,n,wpname); %對數據進行小波包分解
for i=1:2^n %wpcoef(wpt1,[n,i-1])是求第n層第i個節點的係數
E(i)=norm(wpcoef(wpt3,[n,i-1]),2);%求第i個節點的範數平方，其實也就是平方和
end
E_total=sum(E); %求總能量
for i=1:2^n
pthi(i)= E(i)/E_total;%求每個節點的概率
end
pthi

格式：

y = detrend(x) % 消除時間序列中的線性趨勢項
y = detrend(x,'constant') % 消除時間序列中的均值
y = detrend(x,'linear',bp) % 分段消除時間序列中的線性趨勢項,bp爲分段點向量

如

detrend的一個實例

 sig = [0 1 -2 1 0 1 -2 1 0]; % 無線性趨勢的信號
trend = [0 1 2 3 4 3 2 1 0]; % 有兩段線性的趨勢
x = sig+trend; % 將上面趨勢疊加到信號上
y = detrend(x,'linear',5) % 根據指定的分段點去除兩段線性趨勢


y =


   -0.0000    1.0000   -2.0000    1.0000         0    1.0000   -2.0000    1.0000    0.0000