声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1492|回复: 2

[HHT] 关于hhspectrum的结果成图。

[复制链接]
发表于 2009-11-3 16:47 | 显示全部楼层 |阅读模式

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

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

x
我用小波分析信号间断点的相关信息,让后想用hhspectrun分析间断点的时频信息;
新图像:(有四张图的)表示的小波分析的第1,2 层近似系数和1,2 层的细节系数,看出间断点的位置在》50左右的位置。
两张图的是对第2层细节系数(四张图的最后一个)的振幅和频率的分析。我想得到的图像是横坐标是频率,纵坐标是振幅,可图像却像第三张所示,不知道是什么东西;
s=leleccum(1575:1720);
[c,l]=wavedec(s,5,'db3');
a1=wrcoef('a',c,l,'db3',1);
a2=wrcoef('a',c,l,'db3',2);
d1=wrcoef('d',c,l,'db3',1);
d2=wrcoef('d',c,l,'db3',2);
subplot(411);plot(a1);
subplot(412);plot(a2);
subplot(413);plot(d1);
subplot(414);plot(d2);a1,a2,d1,d2为第1,2 层近似系数和1,2层   [A,f,tt]=hhspectrum(d2);细节系数
然后就是画出A,f的图像
plot(f,A);(最后一张的命令);为什么不正确,请大侠致电一下
新图像.JPG
untitled.jpg
1.jpg
回复
分享到:

使用道具 举报

发表于 2009-11-4 14:42 | 显示全部楼层
plot(f,abs(A));试一下
 楼主| 发表于 2009-11-5 15:18 | 显示全部楼层
看了hhspectrum的源代码,A本来就是abs函数求的,所以plot(f,A)和plot(f,abs(A))是一样的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 19:04 , Processed in 0.069249 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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