小波包分解求各频带能量问题？

xiangyu537

clc
fs=12000;
fid=fopen('259.mat','r');
N=1200;
%lev=4;
%nfft=1024;
xdata=fread(fid,N,'int32');
fclose(fid);
%xdata=X247_DE_time;
xdata=xdata(1:N);
xdata=(xdata-mean(xdata))/std(xdata,1);
figure(1);
plot(linspace(0,0.1,N),xdata);
xlabel('时间 t/s');
ylabel('电压 V/v');
[T]=wpdec(xdata,3,'db4');
for i=7:14
    x=wprcoef(T,i);
     E(i)=wenergy(x,i);%求小波包能量
end
  figure(2)
t=1:8;
bar(t,E)
ylabel('能量')
title('外圈故障轴承'）
必须把t=1:8;改为t=1:16才正确  不知为什么？？？？？

ximifly

E,t数组元素不匹配

kaiber

E(i),i=7:14，虽然i只有8个数，但E(i)这样表示就代表E有14个元素，前面6个元素为0。应该把 t 改为1:14才对吧