import wfdb import pywt import numpy as np from matplotlib import pyplot as plt from scipy import signal # 加载数据 record = wfdb.rdrecord('D:/下载/ptb-xl-a-large-publicly-available-electrocardiography-dataset-1.0.3/records100/00000/00001_lr') signal = record.p_signal[:, 0] # 取第一列信号 fs = record.fs # 采样率 # 小波滤波 w = 'db4' # 选用Daubechies4小波 level = 6 # 小波分解层数 coeffs = pywt.wavedec(signal, w, level=level) for i in range(1, level+1): coeffs[i] = pywt.threshold(coeffs[i], np.std(coeffs[i])/2)# 阈值为标准差的一半 filtered_signal = pywt.waverec(coeffs, w) # 去除基线漂移 baseline = signal - filtered_signal filtered_signal += np.mean(baseline) # 去除工频干扰 notch_freq_Hz = np.array([50.0]) # 工频干扰频率为50Hz for freq_Hz in np.nditer(notch_freq_Hz): bp_stop_Hz = freq_Hz + 3.0 * np.array([-1, 1]) b, a = signal.butter(3, bp_stop_Hz / (fs / 2.0), 'bandstop') filtered_signal = signal.lfilter(b, a, filtered_signal) # 绘制图像 plt.figure(figsize=(10, 5)) plt.plot(signal, label='Original') plt.plot(filtered_signal, label='Filtered') plt.title('Filtered Signal') plt.xlabel('Sample')代码如上，出现错误

您未在使用`signal`模块的函数时使用正确的完整名称。在您的代码中，您将`signal`模块导入为`from scipy import signal`，但在代码中您使用了`signal.butter`和`signal.lfilter`，这将导致名称冲突并引发错误。

为了解决这个问题，您可以尝试以下两个解决方案之一：

1. 在使用`signal`模块的函数时使用完整的模块名称。例如，您可以将以下代码：

   b, a = signal.butter(3, bp_stop_Hz / (fs / 2.0), 'bandstop')
   filtered_signal = signal.lfilter(b, a, filtered_signal)

更改为：

   b, a = scipy.signal.butter(3, bp_stop_Hz / (fs / 2.0), 'bandstop')
   filtered_signal = scipy.signal.lfilter(b, a, filtered_signal)

2. 重命名您的导入语句。例如，您可以将以下代码：

   from scipy import signal

更改为：

   import scipy.signal as sig

然后在使用`signal`模块的函数时使用`sig`作为前缀。例如，您可以将以下代码：

   b, a = signal.butter(3, bp_stop_Hz / (fs / 2.0), 'bandstop')
   filtered_signal = signal.lfilter(b, a, filtered_signal)

更改为：

   b, a = sig.butter(3, bp_stop_Hz / (fs / 2.0), 'bandstop')
   filtered_signal = sig.lfilter(b, a, filtered_signal)

希望这些解决方案可以帮助您解决问题。