滤波器是无线信号处理里最常用的一把刀:通带保留要的信号、阻带挡掉不要的频段。本节先立 FIR/IIR 两大体系的取舍,再实操窗函数法设计一只低通,最后把"阶数、滚降、延迟"这把三角天平讲透,并落到它在调制成形、均衡、抗干扰里的实际位置。
如果给一串波形的第一站挑一件工具,十有八九是滤波器。它不是那种"偶尔亮一下相"的配角,而是从头到尾都捏在自己手里的常备:信号过采样前端要挡高频、中频要选台、基带要整形,进到接收端还要均衡、平滑导频。正因为如此,本册开头那些"把波形搬上发射天线"的动作,其实每一步都有滤波器的影子。我们先不急着叩响复杂公式,先给它立两套派系、再亲手做一只,让你往后走到哪一图,都能认出这把刀、知道它此时在替你挡谁、放谁。
阅读完本节,你应当能够:
打开任何一张无线框图,滤波器一定无处不在:前端抗混叠、中频选频道、基带整形、均衡抽头、导频平滑。它的原理其实就一句话——一个 LTI 系统,放某些频率、抑某些频率。可这句"取与舍"落到实际,比听起来麻烦得多,因为它要在"滚得多快、延时多久、算得多累"三个方向里反复打太极。
滤波器分两大体系,各有山头:
FIR(有限冲激响应)。系数就是冲激响应,输出是输入和系数的有限长卷积。优点:天然稳定,可以做成严格线性相位(波形不畸变——通信里最看重的就是"过一趟不画歪"),窗函数法设计直观。缺点:达到同样滚降需要较长滤波器,系数多、计算量大、群延迟大。
IIR(无限冲激响应)。输出还要用一部分过去的输出,因此冲激响应无限拖尾巴。优点:阶数低、同样指标计算量小(比切比雪夫、巴特沃斯等经典原型),延迟小。缺点:不一定稳定(极点必须在单位圆内)、相位非线性明显、无法做成精确线性相位。
| 维度 | FIR | IIR |
|---|---|---|
| 稳定性 | 天然稳定 | 需极点落单位圆内 |
| 线性相位 | 可精确 | 一般非线性 |
| 同指标阶数 | 高、计算量大 | 低、省算了 |
| 群延迟 | 大 | 小 |
| 典型用途 | 脉冲成形、线性相位滤波 | 语音、粗略频选 |
通信处理看重线性相位,所以基带侧钟爱 FIR;射频前端讲究省算省延时,IIR 常被请上台。没有谁总赢,只有"这笔该不该为波的形状买单"。
FIR 最友好的入场式是窗函数法:先构造理想低通的无限冲激响应,再乘上一个有限长窗把它截成 FIR 系数。理想低通的冲激响应是正弦与分母的形式,截它时窗的形状决定滚降与旁瓣。下面直接落地:
import numpy as np def fir_lowpass(cutoff, fs, ntaps): # cutoff: 截止频率;fs:采样率;ntaps:阶数(奇数) t = np.arange(ntaps) - (ntaps-1)//2 h = np.sin(2*np.pi*cutoff/fs*t) / (np.pi*t) h[t == 0] = 2*cutoff/fs # 补上直流点 return h * np.hanning(ntaps) # 加汉宁窗收旁瓣 fs = 8000.0 h = fir_lowpass(cutoff=1000.0, fs=fs, ntaps=51) print("滤波器系数若干:", np.round(h[:5], 4)) # 频响在通带近似1、阻带压低;窗越长滚降越陡 for nt in (21, 51, 101): ht = fir_lowpass(1000.0, fs, nt) print("阶数", nt, "窗长占比", len(ht))
跑几次对比阶数 21/51/101,你会看到滚降越来越陡、但滤波器长度(底延迟)变长。把"用更长系数换更陡滚降"记在心里,就是滤波器设计全程的主演。
拦下来的代价不止延迟,还有过渡带那一截"半吞不吐"。理想的滤波器希望通带到阻带一瞬截止,物理上做不到,落下来必然有一段过渡带宽。窗函数估算过渡带的经验很好用:过渡带宽度大约与"窗主瓣宽度"成正比,而主瓣宽度又约等于 4·π/N。意思是一句话:要过渡带窄,就得把阶数 N 做大,这和滚降陡是完全同一件事的两样说法。设计者通常是先量好"我能容忍多宽的过渡带",反手推回 N,再去验旁瓣够不够低——两条约束一夹,滤波器设计就从"拍脑袋"变成了"查表反算"。
滤波器设计其实在三个标尺上做三角交易:
这句话基本能定案:通信基带里,宁可多花一点滤波器长度换陡滚降与线性相位,也少在相位和波形上妥协;射频里正好反过来,能省一点算力最是要紧。这跟第3.2节"滚降系数买鲁棒性"是同一条哲学在换主次位置。
窗的选择也藏着一条明文:矩形窗前截会导致吉布斯振荡与高的旁瓣,汉宁窗主瓣略宽但旁瓣低得多,凯泽窗用一个滑动的 β 参数在"主瓣宽度"与"旁瓣高度"之间自由取舍。工程上常见做法是先定允许的最大旁瓣(对应带外泄漏约束),再反推需要的窗型与阶数——这又是一次"先定指标、再定资源"的教科书式交易。
| 应用 | 想要的性质 | 为啥要它 |
|---|---|---|
| 脉冲成形 | 线性相位、精确频响 | 收带宽、削 ISI |
| 均衡 | 逆信道响应 | 消除多径畸变 |
| 抗混叠 | 陡阻带 | 采前挡高频 |
再看清波形的"镜"——傅里叶变换与频谱分析怎么把时域看见频域。