石头王石头 发表于 2012-9-17 10:42

关于峭度的程序

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

ddy41 发表于 2012-9-17 11:05

端点发散？两端数据截去了没？

石头王石头 发表于 2012-9-17 17:18

ddy41 发表于 2012-9-17 11:05 static/image/common/back.gif
端点发散？两端数据截去了没？

没有，应用的hht工具包是从本论坛下载的，据说里面是应用镜像法处理的端点效应，不知道里面是否已经应用

ddy41 发表于 2012-9-17 17:33

石头王石头 发表于 2012-9-17 17:18 static/image/common/back.gif
没有，应用的hht工具包是从本论坛下载的，据说里面是应用镜像法处理的端点效应，不知道里面是否已经应用

一定要截的！镜像法是最基本的抑制方法，效果一般

石头王石头 发表于 2012-9-17 19:23

可我不会截呀，明天我问问我师兄吧

石头王石头 发表于 2012-9-20 16:53

石头王石头 发表于 2012-9-17 19:23 static/image/common/back.gif
可我不会截呀，明天我问问我师兄吧

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);
=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 static/image/common/back.gif
一定要截的！镜像法是最基本的抑制方法，效果一般

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);
=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

顶一下，顶一下
页: [1]
查看完整版本: 关于峭度的程序