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個節(jié)點的系數

E(i)=norm(wpcoef(wpt1,[n,i-1]),2);%求第i個節(jié)點的范數平方，其實也就是平方和

end

% disp(‘每個節(jié)點的能量E(i)’);

% E

% disp(‘小波包分解總能量E_total’);

E_total=sum(E); %求總能量

for i=1:2^n

pfir(i)= E(i)/E_total;%求每個節(jié)點的概率

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個節(jié)點的系數

E(i)=norm(wpcoef(wpt2,[n,i-1]),2);%求第i個節(jié)點的范數平方，其實也就是平方和

end

E_total=sum(E); %求總能量

for i=1:2^n

psec(i)= E(i)/E_total;%求每個節(jié)點的概率

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個節(jié)點的系數

E(i)=norm(wpcoef(wpt3,[n,i-1]),2);%求第i個節(jié)點的范數平方，其實也就是平方和

end

E_total=sum(E); %求總能量

for i=1:2^n

pthi(i)= E(i)/E_total;%求每個節(jié)點的概率

end

pthi

總結

以上是生活随笔為你收集整理的python小波包分解_Matlab小波包分解后如何求各频带信号的能量值？的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。