基于月份的颤振滤波器列表

import numpy as np

def chatter_filter(signal, cutoff_freq, sample_rate):
    """
    应用颤振滤波器到输入信号
    :param signal: 输入信号数组
    :param cutoff_freq: 滤波器的截止频率
    :param sample_rate: 采样率
    :return: 滤波后的信号数组
    """
    nyquist_freq = 0.5 * sample_rate
    normalized_cutoff = cutoff_freq / nyquist_freq
    b, a = signal.butter(4, normalized_cutoff, btype='low')
    filtered_signal = signal.lfilter(b, a, signal)
    return filtered_signal

# 示例使用
import matplotlib.pyplot as plt

# 生成一个包含颤振的信号
t = np.linspace(0, 1, 1000, False)
signal = np.sin(2 * np.pi * 50 * t) + 0.5 * np.sin(2 * np.pi * 200 * t)

# 应用颤振滤波器
filtered_signal = chatter_filter(signal, cutoff_freq=150, sample_rate=1000)

# 绘制结果
plt.figure(figsize=(10, 6))
plt.plot(t, signal, label='原始信号')
plt.plot(t, filtered_signal, label='滤波后信号')
plt.legend()
plt.show()