x=beatftr(1,:);
figure;
plot(x)
wpt=wpdec(x,4,'db3'); %小波包分解
plot(wpt);
for j=1:2^4 %4层分解,16组系数
E(j)=sum(abs(wprcoef(wpt,[4,j-1])).^2);
end
E1=sum(E);
dim=length(E);
for j=1:dim
p(j)=E(j)/E1;
end
Shang=-sum(p.*log(p));
figure;
plot(x)
wpt=wpdec(x,4,'db3'); %小波包分解
plot(wpt);
for j=1:2^4 %4层分解,16组系数
E(j)=sum(abs(wprcoef(wpt,[4,j-1])).^2);
end
E1=sum(E);
dim=length(E);
for j=1:dim
p(j)=E(j)/E1;
end
Shang=-sum(p.*log(p));