声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1311|回复: 5

[小波] 求助_振动信号的功率谱计算。

[复制链接]
发表于 2009-4-21 13:40 | 显示全部楼层 |阅读模式

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

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

x
下面这段是我对振动信号进行的去噪处理:
load zhendong.m
s=zhendong(1:10000)
figure(1)
sbuplot(2,2,1);plot(s);
title('原始振动信号');

[thr,sorh,keepapp,cirt]=ddencmp('den','wp',s);
[c,treed,perf0,perfl2]=wpdencmp[s,sorh,3,'db2',crit,thr,keepapp];
subplot(2,2,2);plot(c);
title('默认阀值去噪信号');

thr=thr+15;
……
thr=thr+30;
……

请问,我怎样才能对c信号(默认阀值去噪信号)进行功率谱计算,然后对计算结果作图,获得二维[P,f]和三维[P,f,t]功率谱图
求matlab程序实现。。。
附件是我做的振动波形图。。。
zhendong1-4.jpg

zhendong1.m

59.73 KB, 下载次数: 28

原始数据m文件

回复
分享到:

使用道具 举报

 楼主| 发表于 2009-4-21 16:19 | 显示全部楼层
请问,如果采用PSD生成功率谱,怎么编程才能实现三维功率谱?
发表于 2009-4-21 22:49 | 显示全部楼层
可以用连续小波生成功率谱   1/N*sum(wt(c)^2)/fb   
其中fb为小波带宽
 楼主| 发表于 2009-4-22 13:53 | 显示全部楼层
不明白。。。哎,急死我了
发表于 2009-4-22 14:55 | 显示全部楼层

回复 板凳 aini2009 的帖子

这样做x轴是小波尺度?
还有小波带宽怎么取呢?

[ 本帖最后由 sinada 于 2009-4-22 14:57 编辑 ]
发表于 2009-4-30 10:32 | 显示全部楼层
x轴坐标可以根据尺度和频率的关系转化
带宽选取的标准我也没找到,但是可以自己设定,因为带宽也随尺度成比例变化
所以在计算功率谱的时候在纵坐标上一个等间距的问题。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-26 06:32 , Processed in 0.067570 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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