wuyunfei001 发表于 2010-4-8 00:31

小波包分解频带能量

load wusun
E=zeros(16,19) %0矩阵
for j=1:19
a=wusun(:,j);
depth=4;%分解层数
=wpdec(a,depth,'db4');%小波包分解
for i=1:2^depth
x=wprcoef(T,i);%重构最后一层小波包系数
E(i,j)=wenergy(x,i);%求小波包能量
end
end

我想依次分析19列信号，按照wenergy是能量百分比，分解结果每列
能量相加应该等于100，但是实际上相差很远，
，另外提出任意一列信号，用以上命令来做，此时for j=1:1
分析出来的能量也和图形操作分析出来的能量结果相差很远，
尝试了用每一介能量E除以总能量所得的比例也不是图形操作分析结果的能量比例，
请求大家讨论解决，指出我哪里错了，这么来改？，谢谢！

sparrow 发表于 2010-4-8 02:29

load wusun
E=zeros(16,19)% 0矩阵
for j=1:19
a=wusun(:,j);
depth=4;%分解层数
=wpdec(a,depth,'db4');   %小波包分解
E( :,j)=(wenergy(T))';   %求底层小波包能量
end

或者
load wusun
E=zeros(16,19)% 0矩阵
for j=1:19
a=wusun(:,j);
depth=4;%分解层数
=wpdec(a,depth,'db4');   %小波包分解
for i=1:2^depth
x=wprcoef(T,); %重构最后一层小波包系数
E(i,j)=sum(x.^2); %求小波包能量
end
end
E=100*E./repmat(sum(E),16,1);

WENERGY 直接读取系数计算能量，因为卷积的原因在边界有溢出，算出来的能量应该比WPRCOEF算出来的要大一点。（我没在MATLAB上验证，LOL）

wuyunfei001 发表于 2010-4-8 10:18

对比第二列信号分析结果，wenergy函数结果，最大能量88.486。利用小波交换系数计算能量在计算，最大能量比96.957。
图形窗口操作结果，最大能量bi87.64,到底我以那儿为准确的能量，下一步分析用？欢迎交流，
file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/EL(O6$M1QB$$%7DQBL98ZX5YM.jpg
file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/D%60AYZ0~JYXMSSL7SJO4H4LR.jpgfile:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/0_)JZVM@1P_6.jpg

wuyunfei001 发表于 2010-4-8 10:20

能量比，接上一贴

wuyunfei001 发表于 2010-4-8 10:24

4

e

jingjing1811 发表于 2010-7-14 17:47

初学：
在MATLAB中运行出错，怎么回事？:@)
?? Error using ==> load
Unable to read file wusun: No such file or directory.

Error in ==> Untitled3 at 1
load wusun

ChaChing 发表于 2010-7-14 21:14

回复 6楼 jingjing1811 的帖子

No such file or directory
没有这个档案wusun
页: [1]
查看完整版本: 小波包分解频带能量