ToolPythonIntermediate
Filtering with scipy.signal: mains hum and noise out of an ECG
- y = signal.sosfiltfilt(sos, signal.filtfilt(b, a, x))
- part_out = signal.sosfiltfilt(sos, signal.filtfilt(b, a, part)) # filters are linear: shares add up
- noise_left = signal.sosfiltfilt(sos, signal.filtfilt(b, a, noise))
- ecg_left = signal.sosfiltfilt(sos, signal.filtfilt(b, a, ecg))
- hum_left = signal.filtfilt(b, a, hum)
- print(f"same, method=\"gust\" {np.abs(signal.filtfilt(b, a, hum, method='gust')[:n01]).max():.3f} mV")