python import numpy as np from scipy import signal import matplotlib.pyplot as plt

t = np.linspace # 生成时候序列，长远度为1000个点 x = np.sin + 0.5 * np.sin # 生成包含10Hz和50Hz成分的信号

b, a = signal.butter, 'lowpass', fs=1000) # fs是采样频率

y = signal.filtfilt

plt.figure plt.plot plt.plot plt.xlabel plt.grid plt.legend plt.title plt.show

Fx = np.fft.fft Fy = np.fft.fft freqs = np.fft.fftfreq, d=1/1000) # d是时候间隔，即采样周期

plt.figure plt.plot, label='Original signal') plt.plot, label='Filtered signal') plt.xlabel plt.xlim # 管束频率轴的范围 plt.grid plt.legend plt.title plt.show