xiangyu537 发表于 2009-9-28 17:12

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

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');
=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 发表于 2010-5-10 18:42

E,t数组元素不匹配

kaiber 发表于 2010-5-26 16:32

E(i),i=7:14,虽然i只有8个数,但E(i)这样表示就代表E有14个元素,前面6个元素为0。应该把 t 改为1:14才对吧
页: [1]
查看完整版本: 小波包分解求各频带能量问题?