完整的巴特沃斯高通滤波器实现代码

以下是一个实现巴特沃斯高通滤波器的代码示例，使用Python编写，使用了scipy库中的signal模块来实现滤波器。

import numpy as np
from scipy import signal

def butter_highpass(cutoff, fs, order=5):
    # 计算正常化截止频率
    nyq = 0.5 * fs
    normal_cutoff = cutoff / nyq

    # 设计Butterworth高通滤波器
    b, a = signal.butter(order, normal_cutoff, btype='highpass', analog=False)

    return b, a

def butter_highpass_filter(data, cutoff, fs, order=5):
    # 获取Butterworth高通滤波器的系数
    b, a = butter_highpass(cutoff, fs, order=order)

    # 应用滤波器
    y = signal.filtfilt(b, a, data)

    return y

这个函数可以通过传递要滤波的数据、截止频率、采样率和滤波器阶数来使用。例如，要将一段音频信号从100Hz开始的高频噪声中滤波掉，可以这样调用该函数：

# 读取音频文件
data, fs = librosa.load('audio_file.wav')

# 设定截止频率和滤波器阶数
cutoff = 100
order = 5

# 应用高通滤波器
filtered_data = butter_highpass_filter(data, cutoff, fs, order=order)

值得注意的是，不同的滤波器类型和参数设置会产生不同的滤波效果，因此需要根据实际应用需求进行选择和调整。