从高通到带阻：一文搞懂巴特沃斯数字滤波器的“变形记”（含Matlab/Octave实例）

从高通到带阻：一文搞懂巴特沃斯数字滤波器的“变形记”（含Matlab/Octave实例）

在音频均衡器设计或生物电信号处理中，工程师常面临一个核心挑战：如何将基础的低通滤波器转化为适应不同频段需求的高通、带通或带阻滤波器？巴特沃斯滤波器因其最大平坦特性成为首选，但类型转换的数学原理常让实践者望而生畏。本文将用工程视角拆解频率变换的底层逻辑，并通过可复现的代码示例展示从理论到实践的完整路径。

1. 巴特沃斯滤波器的设计哲学与频率变换原理

巴特沃斯滤波器的核心优势在于其极点均匀分布的特性——在复平面上，N阶滤波器的极点等距分布在单位圆的左半平面。这种几何对称性带来了通带内的最大平坦响应，但也意味着：

• 低通原型（Prototype）的传递函数形式最为简洁
• 其他类型滤波器可通过频率变量替换（Frequency Transformation）派生

关键变换公式（模拟域）：

% 低通到高通的频率变换
omega_c = 2*pi*fc;  % 截止频率(rad/s)
H_highpass = subs(H_lowpass, s, omega_c/s);  % 变量替换

这种变换的本质是复平面上的保角映射：

• 低通的通带（s=0附近）映射为高通的阻带（s→∞）
• 低通的阻带（s→∞）映射为高通的通带（s=0附近）

注意：数字滤波器需先通过双线性变换将模拟传递函数转换为数字域，再进行频率变换，否则会导致非线性频率畸变。

2. 四类滤波器的Matlab实现对比

2.1 低通滤波器设计基准

fs = 1000;  % 采样率(Hz)
fc = 100;   % 截止频率(Hz)
[