一维时间序列的经验模态分解(以轴承振动信号为例,Python)

import numpy as np
from pyhht import EMD
import matplotlib.pyplot as plt
from pyhht.visualization import plot_imfs
from time import time
import lee_dataset_CWRU
from lee_dataset_CWRU import *
DE_H1, FE_H1, t_DE_H1, t_FE_H1, RPM_H1, samples_s_DE_H1, samples_s_FE_H1 = lee_dataset('97.mat')
x = DE_H1
t = t_DE_H1
plt.figure(figsize=(7, 4))
plt.plot(t, x)

time0 = time()
decomposer = EMD(x)
imfs = decomposer.decompose()
elapsed_time = time() - time0
Mode 2, Iteration 600
2 minima > 0; 19594 maxima < 0.
Mode 2, Iteration 800
1 minima > 0; 19640 maxima < 0.
Mode 2, Iteration 1000
1 minima > 0; 19677 maxima < 0.
Mode 4, Iteration 1200
3 minima > 0; 7410 maxima < 0.
Mode 5, Iteration 1400
2 minima > 0; 4185 maxima < 0.
Mode 6, Iteration
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

哥廷根数学学派

码字不易,且行且珍惜

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值