当前位置：首页 > news >正文

别再只用巴特沃斯了！用MATLAB的cheby1函数快速搞定带通滤波器设计（附完整代码）

news 2026/5/27 4:48:28

别再只用巴特沃斯了！用MATLAB的cheby1函数快速搞定带通滤波器设计（附完整代码）

在信号处理领域，滤波器设计是每个工程师都绕不开的基础技能。提到IIR滤波器，很多人的第一反应就是巴特沃斯（Butterworth）——它简单、稳定，通带平坦，似乎是个"不会出错"的选择。但当你面对一个需要从嘈杂的传感器信号中提取特定频段的实际工程问题时，是否考虑过：同样的性能指标下，是否存在计算效率更高的方案？

契比雪夫I型滤波器（Chebyshev Type I）就是这样一个常被低估的利器。与巴特沃斯相比，它在通带内允许等波纹波动，却能以更低的滤波器阶数实现相同的阻带衰减——这意味着更少的计算资源消耗和更快的实时处理速度。本文将带你用MATLAB的cheby1函数家族，从参数计算到效果验证，完整实现一个高性能带通滤波器设计。

1. 为什么选择契比雪夫I型滤波器？

1.1 巴特沃斯的局限

巴特沃斯滤波器的幅频响应在通带和阻带都是单调变化的，这种特性带来了两个固有缺陷：

过渡带平缓：要达到理想的阻带衰减，往往需要较高的滤波器阶数
计算成本高：高阶滤波器意味着更多的乘加运算，对实时系统构成压力

下表对比了两种滤波器在相同指标下的阶数需求：

指标要求	巴特沃斯阶数	契比雪夫I型阶数
通带波纹1dB	8	5
阻带衰减40dB
过渡带宽度100Hz

1.2 契比雪夫的优势

契比雪夫I型滤波器通过允许通带内等波纹波动，换取了更陡峭的过渡带特性。这种设计哲学特别适合以下场景：

对通带平坦度要求不严苛的场合（如语音信号处理）
需要降低处理器负载的嵌入式系统
宽带信号中提取窄带成分的应用

注意：当信号对通带相位响应敏感时，可能需要考虑椭圆滤波器或FIR方案。

2. 设计流程与MATLAB工具链

2.1 核心函数三剑客

MATLAB提供了完整的契比雪夫I型设计工具链：

cheb1ord- 计算最小阶数和截止频率
```
[n, Wn] = cheb1ord(Wp, Ws, Rp, Rs)
```
- Wp：通带边缘频率（归一化）
- Ws：阻带起始频率（归一化）
- Rp：通带最大波纹（dB）
- Rs：阻带最小衰减（dB）
cheby1- 生成滤波器系数
```
[b, a] = cheby1(n, Rp, Wn, 'bandpass')
```
freqz- 验证频率响应
```
freqz(b, a, 1024, Fs)
```

2.2 参数选择实战技巧

频率归一化：所有频率参数应除以奈奎斯特频率（采样率/2）
通带波纹：通常设置在0.1-1dB之间，过大会导致信号失真
阻带衰减：根据干扰强度决定，一般不少于40dB

3. 完整带通滤波器设计案例

假设我们需要从采样率为10kHz的ECG信号中提取80-120Hz的心率成分，设计指标如下：

通带：80-120Hz (±3dB)
阻带：<60Hz和>140Hz处衰减≥50dB
采样率：10kHz

3.1 步骤实现

Fs = 10000; % 采样率 Wp = [80 120]/(Fs/2); % 通带边缘频率(归一化) Ws = [60 140]/(Fs/2); % 阻带边缘频率(归一化) Rp = 3; % 通带波纹(dB) Rs = 50; % 阻带衰减(dB) % 计算最小阶数 [n, Wn] = cheb1ord(Wp, Ws, Rp, Rs); % 生成带通滤波器 [b, a] = cheby1(n, Rp, Wn, 'bandpass'); % 绘制频率响应 freqz(b, a, 1024, Fs); title(['Chebyshev Type I Bandpass Filter: Order = ' num2str(n)]);

3.2 效果优化技巧

如果阶数过高，可以适当放宽通带波纹要求
使用fvtool交互式工具调整参数：
```
fvtool(b, a)
```

对于实时处理，可将系数导出为C头文件：

fid = fopen('filter_coeffs.h','w'); fprintf(fid, 'const float b[] = {%f, %f, %f};\n', b); fprintf(fid, 'const float a[] = {%f, %f, %f};\n', a); fclose(fid);

4. 进阶应用与性能调优

4.1 多级滤波器设计

对于特别严格的指标要求，可以采用级联设计：

% 第一级：低通 [n1, Wn1] = cheb1ord(120/(Fs/2), 140/(Fs/2), Rp, Rs/2); [b1, a1] = cheby1(n1, Rp, Wn1); % 第二级：高通 [n2, Wn2] = cheb1ord(80/(Fs/2), 60/(Fs/2), Rp, Rs/2); [b2, a2] = cheby1(n2, Rp, Wn2, 'high'); % 级联响应 H1 = tf(b1, a1, 1/Fs); H2 = tf(b2, a2, 1/Fs); Hcas = H1 * H2; bode(Hcas)