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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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 · 发表于 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才对吧

账号		自动登录	找回密码
密码			我要加入

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。