clear;
clc;
n=5;
wpname='sym8';
% [b,a]=butter(8,[5/100 99/100]);
% load a1_1-1;
x=xlsread('dbxyp.xls');
Data=x(:,1);
% Data=filter(b,a,Data);
Data=Data-mean(Data);
Data=detrend(Data);
wpt1=wpdec(Data,5,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
评论0