用iircase5设计一个巴特沃斯滤波器写出代码

以下是使用iircase5设计巴特沃斯滤波器的代码示例，其中包含了系数计算和滤波操作两部分：

import scipy.signal as signal

# 设计巴特沃斯滤波器
fs = 1000  # 采样频率
fc = 50  # 截止频率
b, a = signal.butter(8, fc/(fs/2), 'lowpass')

# 模拟输入信号
import numpy as np
t = np.arange(0, 1, 1/fs)
x = np.sin(2*np.pi*100*t) + np.sin(2*np.pi*200*t) + np.sin(2*np.pi*300*t)

# 滤波操作
y = signal.filtfilt(b, a, x)

# 绘制输入信号和滤波输出信号的时域波形
import matplotlib.pyplot as plt
plt.figure()
plt.plot(t, x, label='Input')
plt.plot(t, y, label='Output')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')
plt.legend()
plt.show()

需要注意的是，以上代码中的系数计算使用了`scipy.signal.butter`函数，该函数的第一个参数表示滤波器的阶数，第二个参数表示截止频率，第三个参数表示滤波器类型（例如`lowpass`表示低通滤波器）。滤波操作使用了`scipy.signal.filtfilt`函数，该函数可以对输入信号进行前向和后向滤波，避免了滤波器的相位延迟。最后，绘制了输入信号和滤波输出信号的时域波形，以便于观察滤波效果。