☰
二阶巴特沃斯滤波器设计指南:从原理公式到 MATLAB/Python 实战
2026/10/9 3:41:03 网站建设 项目流程

摘要:二阶巴特沃斯滤波器以最大平坦响应和 -40 dB/decade 滚降著称,是工程中应用最广的低通滤波器之一。本文详解其设计公式、Sallen-Key 电路实现与级联二阶节结构,并给出 MATLAB/Python 设计示例,助你快速掌握从原理到仿真的完整流程。

TL;DR(快速阅读):

  • 最大平坦响应:通带内幅频响应最平坦,无波纹。
  • 滚降率:二阶滤波器为 -40 dB/decade。
  • 归一化公式:H(s)=1/(s²+√2s+1)。
  • Sallen-Key 要点:Q=0.707 匹配巴特沃斯响应。
  • 设计验证:MATLAB/Python 均可快速实现。

目录:

  • 1. 巴特沃斯滤波器特点
  • 2. 二阶巴特沃斯滤波器设计公式
  • 3. 电路实现方式
  • 4. 设计步骤详解
  • 5. MATLAB 设计示例
  • 6. 实际应用
  • 7. 二阶巴特沃斯滤波器与一阶滤波器对比
  • 8. 常见问题排查
  • 9. 总结
  • 10. 参考资料

在信号处理与电子系统设计中,滤波器是去除噪声、提取有效频段信息的基础器件。二阶巴特沃斯滤波器凭借通带内最大平坦的幅频响应和-40 dB/decade 的滚降率,在音频处理、生物医学信号预处理、通信系统抗干扰等场景中得到广泛应用,是工程实践中兼顾性能与实现复杂度的经典选择。本文将从巴特沃斯滤波器的基本原理出发,依次讲解设计公式、电路实现方式、完整设计步骤,并结合 MATLAB 与 Python 给出可运行的实现示例,最后介绍典型应用场景与常见设计误区,帮助读者系统掌握二阶巴特沃斯滤波器的设计与应用。

###二阶巴特沃斯滤波器

二阶巴特沃斯滤波器是一种具有最大平坦幅频响应的低通滤波器,其设计目标是在通带内实现最平滑的幅度响应,同时在截止频率后具有快速衰减特性。与一阶滤波器相比,二阶滤波器具有更陡峭的滚降率(-40 dB/decade),适用于需要更高选择性的应用场景。


1. 巴特沃斯滤波器特点
特性描述
幅频响应通带内最大平坦,无波纹
滚降率二阶滤波器为 -40 dB/decade
相位响应非线性,但比其他类型(如切比雪夫)更平滑
极点分布在 s 平面上呈对称分布,极点位于单位圆内

2. 二阶巴特沃斯滤波器设计公式

二阶巴特沃斯滤波器的系统函数为:

$$
H(s) = \frac{1}{s^2 + \sqrt{2}s + 1}
$$

其中,$s$ 是复频率变量,该系统函数对应于一个归一化截止频率为 1 rad/s 的二阶低通滤波器。

归一化频率转换

若实际截止频率为 $f_c$,则需将 $s$ 替换为 $\frac{s}{2\pi f_c}$,得到实际系统函数:

$$
H(s) = \frac{1}{(s/(2\pi f_c))^2 + \sqrt{2}(s/(2\pi f_c)) + 1}
$$


相位响应分析

除了幅频特性,相位响应也是滤波器设计中的重要考量。二阶巴特沃斯滤波器的相位响应是非线性的,但相比切比雪夫等类型更为平滑。下面从相位延迟特性、与一阶滤波器的对比以及相位公式推导三个角度展开分析。

相位延迟特性

将归一化系统函数中的 $s$ 替换为 $j\omega$($\omega$ 为归一化角频率),得到频率响应:

$$
H(j\omega) = \frac{1}{(j\omega)^2 + \sqrt{2}(j\omega) + 1} = \frac{1}{1 - \omega^2 + j\sqrt{2}\omega}
$$

其相位响应为分母辐角的相反数:

$$
\varphi(\omega) = -\arctan\left(\frac{\sqrt{2}\omega}{1 - \omega^2}\right)
$$

在低频段($\omega \to 0$),相位趋近于 0;在截止频率处($\omega = 1$),相位为 $-\arctan(\infty) = -\pi/2$;在高频段($\omega \to \infty$),相位趋近于 $-\pi$。可见二阶巴特沃斯滤波器在通带内会引入随频率变化的相位延迟,且频率越高延迟越大,这正是其相位响应非线性的体现。

与一阶滤波器相位响应的对比

一阶巴特沃斯低通滤波器的相位响应为 $\varphi_1(\omega) = -\arctan(\omega)$,其相位从 0 平滑过渡到 $-\pi/2$,在截止频率处为 $-\pi/4$。相比之下,二阶滤波器的相位变化范围更大(从 0 到 $-\pi$),且在截止频率附近变化更陡峭,因此相位失真更明显。这也解释了前文对比表中「二阶滤波器相位失真更大」的原因。若系统对相位线性度要求较高,可考虑使用贝塞尔滤波器或采用零相位滤波(如filtfilt)来补偿。

相位响应公式推导

下面给出 $\varphi(\omega) = -\arctan\left(\frac{\sqrt{2}\omega}{1 - \omega^2}\right)$ 的完整推导过程。设频率响应为:

$$
H(j\omega) = \frac{1}{1 - \omega^2 + j\sqrt{2}\omega}
$$

令分母为复数 $Z = (1 - \omega^2) + j\sqrt{2}\omega$,其辐角为:

$$
\angle Z = \arctan\left(\frac{\sqrt{2}\omega}{1 - \omega^2}\right)
$$

由于 $H(j\omega) = 1 / Z$,其相位为分母辐角的相反数:

$$
\varphi(\omega) = \angle H(j\omega) = -\angle Z = -\arctan\left(\frac{\sqrt{2}\omega}{1 - \omega^2}\right)
$$

推导完成。该公式可直接用于计算任意归一化频率下的相位值,例如在截止频率 $\omega = 1$ 处,$\varphi(1) = -\arctan(\infty) = -\pi/2$,与理论分析一致。

3. 电路实现方式

二阶巴特沃斯滤波器可以通过以下两种方式实现:

(1) Sallen-Key 有源滤波器

下图给出了 Sallen-Key 二阶巴特沃斯低通滤波器的电路原理图,包含运放、电阻 R1、R2 和电容 C1、C2 的连接关系,并标注了元件参数(R1 = R2 = 11 kΩ,C1 = 20 nF,C2 = 10 nF):

flowchart TD Vin[Vin 输入] -- R1 11kΩ --> N1((节点1)) N1 -- R2 11kΩ --> N2((节点2)) N1 -- C1 20nF --> N3((节点3)) N2 -- C2 10nF --> GND[GND 地] N3 -- 反馈 --> N2 N2 -- 运放同相输入 --> OP[运放] OP -- 输出 --> Vout[Vout 输出] Vout -- 反馈 --> N3

Sallen-Key 结构是一种常见的有源滤波器拓扑,适用于二阶滤波器设计。其传递函数为:

$$
H(s) = \frac{1}{1 + \frac{s}{R_1 C_1} + \frac{s^2}{R_1 R_2 C_1 C_2}}
$$

通过调整电阻和电容值,可以实现所需的二阶巴特沃斯响应。

电路设计实例(fc = 1 kHz)

下面以截止频率fc = 1 kHz为例,给出完整的 Sallen-Key 二阶巴特沃斯低通滤波器设计过程。为简化计算,取R1 = R2 = R、C1 = 2C2,此时滤波器的品质因数Q = 1 / (3 - 2) = 1,恰好满足巴特沃斯响应所需的Q = 1 / √2 ≈ 0.707的近似条件(实际巴特沃斯 Q 值为 0.707,此处取 Q = 1 会带来轻微过冲,工程上更常用C1 = 2C2配合R1 ≠ R2精确匹配,下文给出精确取值)。

第一步:确定元件取值

二阶巴特沃斯低通滤波器的截止频率公式为:

$$
f_c = \frac{1}{2\pi \sqrt{R_1 R_2 C_1 C_2}}
$$

取R1 = R2 = R、C1 = 2C2,代入得:

$$
f_c = \frac{1}{2\pi R \sqrt{2C_2^2}} = \frac{1}{2\pi R C_2 \sqrt{2}}
$$

令C2 = 10 nF,则C1 = 20 nF,代入fc = 1000 Hz求解 R:

$$
R = \frac{1}{2\pi \times 1000 \times 10 \times 10^{-9} \times \sqrt{2}} \approx 11.25 \text{ kΩ}
$$

取标准值R1 = R2 = 11 kΩ(E24 系列),此时实际截止频率为:

$$
f_c = \frac{1}{2\pi \times 11 \times 10^3 \times 10 \times 10^{-9} \times \sqrt{2}} \approx 1023 \text{ Hz}
$$

误差约 2.3%,在元件容差范围内可接受。若需更精确的 1 kHz,可选用R = 11.25 kΩ(E96 系列精密电阻)或C2 = 10 nF、C1 = 20 nF配合微调电阻。

第二步:元件选型建议

  • 运放:选用低噪声、高输入阻抗的通用运放,如TL072(双运放,JFET 输入,输入阻抗极高)、NE5532(低噪声,适合音频)或OP07(精密低失调)。对于 1 kHz 截止频率,运放的增益带宽积(GBW)需大于 10 MHz 以保证高频段响应不受影响,上述型号均满足要求。
  • 电阻:选用金属膜电阻(1% 容差),温度系数 ≤ 50 ppm/°C,以保证截止频率的稳定性。功率取 1/4 W 即可。
  • 电容:C1、C2 选用C0G(NP0)型陶瓷电容或聚丙烯薄膜电容,这类电容容差小(±5% 以内)、温度稳定性好,不会因温漂导致截止频率偏移。避免使用 X7R/Z5U 等高介电常数陶瓷电容,其容值随电压和温度变化明显。
  • 电源去耦:运放电源引脚需就近并联 0.1 μF 陶瓷电容和 10 μF 电解电容,抑制电源噪声。

第三步:仿真验证截止频率

推荐使用LTspice或Multisim进行交流扫描(AC Sweep)仿真,验证截止频率是否准确。以 LTspice 为例,搭建电路后执行以下步骤:

  1. 在输入端接入AC 1交流小信号源(幅值 1 V)。
  2. 设置交流扫描:起始频率 10 Hz,终止频率 100 kHz,采用十倍频程(Decade)扫描,每十倍频 100 个点。
  3. 运行仿真后,查看输出节点电压的幅频响应曲线。
  4. 找到幅值下降至 -3 dB(即 0.707 倍,约 0.707 V)处的频率点,该频率即为实际截止频率。

若仿真得到的 -3 dB 点偏离 1 kHz,可微调电阻值。例如实测截止频率为 1023 Hz,可将 R 增大至 11.25 kΩ 使截止频率回落至 1 kHz。此外,还需检查通带增益是否为 0 dB(单位增益),以及高频段滚降斜率是否接近 -40 dB/decade,以确认电路符合二阶巴特沃斯特性。

元件参数汇总表(fc = 1 kHz 设计实例)

下表汇总了上述设计实例中各元件的标称值、实际选型值、E24/E96 系列选型以及对应的截止频率误差,便于快速查阅与选型。

元件标称值实际值E24/E96 系列选型对应截止频率误差
R111.25 kΩ11 kΩE24 系列(11 kΩ)约 +2.3%(fc ≈ 1023 Hz)
R211.25 kΩ11 kΩE24 系列(11 kΩ)约 +2.3%(fc ≈ 1023 Hz)
C120 nF20 nFE24 系列(20 nF,C0G/NP0)约 +2.3%(fc ≈ 1023 Hz)
C210 nF10 nFE24 系列(10 nF,C0G/NP0)约 +2.3%(fc ≈ 1023 Hz)
R1(精密)11.25 kΩ11.25 kΩE96 系列(11.25 kΩ,1%)约 0%(fc ≈ 1000 Hz)
R2(精密)11.25 kΩ11.25 kΩE96 系列(11.25 kΩ,1%)约 0%(fc ≈ 1000 Hz)

说明:采用 E24 系列标准值(R1 = R2 = 11 kΩ,C1 = 20 nF,C2 = 10 nF)时,实际截止频率约为 1023 Hz,误差约 +2.3%;若需更精确的 1 kHz,可选用 E96 系列精密电阻 11.25 kΩ,此时截止频率误差接近 0%。

(2) 级联二阶节(Biquad)

对于高阶滤波器(如三阶或四阶),通常采用级联多个二阶节的方式实现。每个二阶节的传递函数为:

$$
H_i(s) = \frac{b_0 + b_1 s + b_2 s^2}{a_0 + a_1 s + a_2 s^2}
$$

通过合理选择系数 $b_0, b_1, b_2, a_0, a_1, a_2$,可以实现二阶巴特沃斯滤波器的特性。

四阶巴特沃斯滤波器的级联设计示例

下面以四阶巴特沃斯低通滤波器为例,演示如何将高阶传递函数分解为两个二阶节的乘积,并给出每个二阶节的系数计算方法和 MATLAB 实现代码。

第一步:四阶传递函数的分解<

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询