声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5604|回复: 17

[综合讨论] 信号小波包分解各频带能量比问题

  [复制链接]
发表于 2012-5-30 16:17 | 显示全部楼层 |阅读模式

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

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

x
用matlab对信号进行傅里叶变换并进行5层小波包分解,用wenergy()得到的各频段能量比为啥与信号的频谱图上该频率的幅值不是对应的呢,比如频谱在150kHZ上幅值很大,有尖峰,可是为什么在频段125~140.625  140.625~156.25  156.25 ~171.875 的能量比都很小。求各位高手指点啊
回复
分享到:

使用道具 举报

头像被屏蔽
发表于 2012-5-30 16:21 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2012-5-31 22:10 | 显示全部楼层
楼主,我也碰到这样的问题,发现好像图形对的起来,但是图谱就对不起来了。望高手路过
 楼主| 发表于 2012-6-4 15:15 | 显示全部楼层
回复 3 # wingxy_123 的帖子

图形对的起来是指什么图形啊
 楼主| 发表于 2012-6-5 19:52 | 显示全部楼层
有没有高手帮忙解答一下啊
发表于 2012-6-10 21:38 | 显示全部楼层
额,我的问题解决了,是采样频率的问题,你对比一下fft幅值谱与小波包能量谱。我的发现少了一块所以对不起来了

评分

1

查看全部评分

 楼主| 发表于 2012-6-13 15:17 | 显示全部楼层
我检查了好多遍了。。。可能当局者迷吧。。。这是我的程序,希望大家帮我纠纠错。
clear;
fs=1000000;
fid=fopen('2012-06-12-10-55-59-1.txt','r');
N=200000;
x=fscanf(fid,'%f',N);
fclose(fid);
figure(1);
subplot(3,1,1);
plot(x);
ylabel('幅值 A');
xlabel('样本序号 n');
nfft=length(x);
y1=abs(fft(x,nfft));
y=y1/N*2;
f=(0:nfft-1)*fs/nfft/1000;
subplot(3,1,2);
plot(f,y);
axis([0,300,min(y),max(y)]);
ylabel('幅值 A');
xlabel('频率 f/kHz');
%------------------------------------------------------
wpt=wpdec(x,5,'db10','shannon');
E = wenergy(wpt);
e=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];
subplot(3,1,3);
bar(E(e),'group');
E(e)
发表于 2012-6-15 08:52 | 显示全部楼层
楼主,你在写毕业论文吧,能不能加个qq,或者msn啥的聊,这样实在太慢了
 楼主| 发表于 2012-6-19 09:33 | 显示全部楼层
回复 8 # wingxy_123 的帖子

加了你Q:@)
发表于 2012-6-19 11:31 | 显示全部楼层
恩恩,加了

点评

希望能将你的讨论结果与大家分享!  发表于 2012-7-5 09:46
发表于 2012-9-30 10:22 | 显示全部楼层
那个e向量是什么?是表示数据的列吗
发表于 2012-10-14 09:34 | 显示全部楼层
楼主的结果呢?为什么用wenergy求出的能量比例,与自己用范数、求和、相除计算得到的不一样?
发表于 2012-10-18 11:46 | 显示全部楼层

wenergy是采用小波系数平方,然后求和,并进行归一化得到的
发表于 2012-10-25 10:33 | 显示全部楼层
happy 发表于 2012-10-18 11:46
wenergy是采用小波系数平方,然后求和,并进行归一化得到的

那么,应该采用哪一个?wenergy?
发表于 2012-11-20 22:34 | 显示全部楼层
学习过一段时间,程序没错,但是好像频段的选择是个问题
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 00:25 , Processed in 0.089805 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表