EMD降噪python代码

以下是一个简单的EMD降噪Python代码示例：

import numpy as np
import matplotlib.pyplot as plt
from PyEMD import EMD

# 生成信号
t = np.linspace(0, 1, 1000)
x = np.sin(2*np.pi*5*t) + 0.2*np.sin(2*np.pi*20*t)

# EMD降噪
emd = EMD()
IMFs = emd(x)
residual = IMFs[-1]
clean_signal = x - residual

# 绘制结果
plt.subplot(2, 1, 1)
plt.plot(t, x)
plt.title('Original Signal')
plt.subplot(2, 1, 2)
plt.plot(t, clean_signal)
plt.title('Clean Signal')
plt.tight_layout()
plt.show()

这个代码示例使用PyEMD库中的EMD类来实现EMD降噪。首先，我们生成一个包含5Hz和20Hz正弦波的信号。然后，我们使用EMD类将信号分解为一组内在模态函数（IMFs）。最后，我们从IMFs中获取残差，并将其从原始信号中减去以获得降噪信号。

请注意，这只是一个简单的示例，可能需要进行更多的调整和优化，以适用于更复杂的信号和降噪任务。