在一维心电信号中进行加噪的python代码
时间: 2024-05-12 16:18:43 浏览: 45
以下是一个简单的Python代码示例,用于在一维心电信号中添加高斯白噪声:
```python
import numpy as np
import matplotlib.pyplot as plt
# 生成原始信号
fs = 1000 # 采样频率
t = np.linspace(0, 1, fs) # 时间轴
f = 60 # 心率
ecg_signal = np.sin(2 * np.pi * f * t)
# 添加高斯白噪声
mean = 0 # 均值
std = 0.1 # 标准差
noise = np.random.normal(mean, std, size=len(ecg_signal))
noisy_ecg = ecg_signal + noise
# 绘制原始信号和加噪信号
plt.figure(figsize=(10, 6))
plt.plot(t, ecg_signal, label='Original ECG signal')
plt.plot(t, noisy_ecg, label='Noisy ECG signal')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')
plt.title('ECG signal with Gaussian white noise')
plt.legend()
plt.show()
```
在此示例中,我们首先生成一个心电信号作为原始信号。然后,我们使用`numpy.random.normal`函数生成一些高斯白噪声,并将其添加到原始信号中。最后,我们使用`matplotlib.pyplot`库绘制原始信号和加噪信号。在这个简单的例子中,我们使用了一个固定的均值和标准差来生成高斯白噪声,但实际上,您可能需要根据您的数据和应用程序来调整这些参数。
阅读全文