EMD分解去噪代码实例 python

以下是使用Python实现EMD分解去噪的代码示例：

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

# 生成测试数据
t = np.linspace(0, 1, 1000)
s = np.sin(20 * np.pi * t) + np.cos(10 * np.pi * t)

# 使用EMD分解去噪
emd = EMD()
IMFs = emd(s)
s_clean = IMFs[-1]

# 绘制原始信号和去噪后的信号
plt.subplot(2, 1, 1)
plt.plot(t, s)
plt.title('Original Signal')
plt.xlabel('Time')
plt.ylabel('Amplitude')
plt.subplot(2, 1, 2)
plt.plot(t, s_clean)
plt.title('Clean Signal')
plt.xlabel('Time')
plt.ylabel('Amplitude')
plt.show()

在上面的代码中，我们首先生成了一个包含正弦和余弦波组合的测试信号。然后，我们使用PyEMD库中的EMD类进行了EMD分解去噪。最后，我们将分解后的信号的最后一个分量作为去噪后的信号，并绘制了原始信号和去噪后的信号的图形。

请注意，这个示例只是演示了如何使用EMD分解去噪，实际应用中可能需要进行更多的调整和优化。