声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1838|回复: 7

[综合讨论] 关于峭度的程序

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

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

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

x
我看了一遍文章,用eemd进行峭度分解,分解出来的波形挺好的,我就是了一把,用的数据是西储大学的轴承数据,结果算出来的结果不对,差太远了,到底怎么回事?请教高手
回复
分享到:

使用道具 举报

发表于 2012-9-17 11:05 | 显示全部楼层
端点发散?两端数据截去了没?
 楼主| 发表于 2012-9-17 17:18 | 显示全部楼层

没有,应用的hht工具包是从本论坛下载的,据说里面是应用镜像法处理的端点效应,不知道里面是否已经应用
发表于 2012-9-17 17:33 | 显示全部楼层
石头王石头 发表于 2012-9-17 17:18
没有,应用的hht工具包是从本论坛下载的,据说里面是应用镜像法处理的端点效应,不知道里面是否已经应用

一定要截的!镜像法是最基本的抑制方法,效果一般
 楼主| 发表于 2012-9-17 19:23 | 显示全部楼层
可我不会截呀,明天我问问我师兄吧
 楼主| 发表于 2012-9-20 16:53 | 显示全部楼层
石头王石头 发表于 2012-9-17 19:23
可我不会截呀,明天我问问我师兄吧

clear all;
close all;
load('IR014_1.mat','y')
load('zhengchang.mat','z')
subplot(211);plot(z);xlabel('点数');ylabel('幅值');title('正常信号');
subplot(212);plot(y);xlabel('点数');ylabel('幅值');title('故障信号');
fs=12000;
N=2400;
n=0:N-1;
t=n/fs;
f=1./(2400/12000)*n;
c=emd(y);
[r,s]=size(c);
figure;
for q=1:r
    subplot(r,1,q);plot(t,c(q,:));xlabel('时间/t');ylabel('幅值');
end
for q=1:r-1
  a(q)=kurtosis(c(q,101:2300));
end
峭度是这样写吗?怎么总感觉运行不对呀,请高手指点
 楼主| 发表于 2012-9-25 10:05 | 显示全部楼层
ddy41 发表于 2012-9-17 17:33
一定要截的!镜像法是最基本的抑制方法,效果一般

clear all;
close all;
load('IR014_1.mat','y')
load('zhengchang.mat','z')
subplot(211);plot(z);xlabel('点数');ylabel('幅值');title('正常信号');
subplot(212);plot(y);xlabel('点数');ylabel('幅值');title('故障信号');
fs=12000;
N=2400;
n=0:N-1;
t=n/fs;
f=1./(2400/12000)*n;
c=emd(y);
[r,s]=size(c);
figure;
for q=1:r
    subplot(r,1,q);plot(t,c(q,:));xlabel('时间/t');ylabel('幅值');
end
for q=1:r-1
  a(q)=kurtosis(c(q,101:2300));
end
峭度是这样写吗?怎么总感觉运行不对呀,请指点
还有我不会截取端点,怎么截呀?

 楼主| 发表于 2013-7-17 15:24 | 显示全部楼层
顶一下,顶一下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-10 08:32 , Processed in 0.073949 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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