4、FIR滤波器设计（窗函数法）：矩形窗、汉宁窗、海明窗、布莱克曼窗设计实例
各位同学，咱们今天聊聊FIR滤波器设计里最经典、也最直观的方法——窗函数法。说实话，我刚入行那会儿，觉得这玩意儿不就是拿个窗函数截一下理想滤波器嘛，有啥好学的？后来在项目里吃过亏才明白，选错一个窗，整个系统可能就白干了。今天咱们就把矩形窗、汉宁窗、海明窗、布莱克曼窗这四位“老兄”请出来，好好盘一盘。
4.1 窗函数法的核心思想
先说说基本原理。你想想看，理想的低通滤波器在频域是个矩形，但时域是无限长的sinc函数。这玩意儿没法实现啊！怎么办？我们只能截取一段，就像用窗户看风景，只能看到窗户框住的那部分。这个“窗户”就是窗函数。
说白了，窗函数法就是三步走：
- 构造理想滤波器的单位脉冲响应 hd(n)，这玩意儿通常是无限长的。
- 选一个窗函数 w(n)，长度 N 自己定。
- 两者相乘：h(n) = hd(n) · w(n)，得到实际可实现的 FIR 系数。
嗯，这里要注意：窗函数的选择直接决定了滤波器的过渡带宽度和阻带衰减。我在项目中遇到过，有人为了追求窄过渡带，选了矩形窗，结果阻带衰减只有-21dB，带外干扰压不下去，整个系统指标全废了。所以，选窗这事儿，真得慎重。
4.2 四种经典窗函数对比
咱们直接上干货，看看这四种窗函数到底有啥区别。
| 窗函数类型 | 主瓣宽度 | 旁瓣峰值衰减 (dB) | 阻带最小衰减 (dB) | 过渡带宽度 (近似) |
|---|---|---|---|---|
| 矩形窗 | 4π/N | -13 | -21 | 0.9π/N |
| 汉宁窗 | 8π/N | -31 | -44 | 3.1π/N |
| 海明窗 | 8π/N | -41 | -53 | 3.3π/N |
| 布莱克曼窗 | 12π/N | -57 | -74 | 5.5π/N |
看到这个表，你可能会问：矩形窗过渡带最窄，为啥不用它？我告诉你，矩形窗的旁瓣只衰减-13dB，这意味着阻带里会有大量能量泄漏。我曾经在一个雷达脉冲压缩项目里试过矩形窗，结果距离旁瓣把弱目标全淹了，那叫一个惨。所以，除非你对阻带衰减完全没要求，否则别碰矩形窗。
4.3 设计实例：低通滤波器
咱们来点实际的。假设要设计一个低通FIR滤波器，指标如下：
- 采样频率 fs = 1000 Hz
- 通带截止频率 fp = 100 Hz
- 阻带起始频率 fst = 150 Hz
- 阻带衰减 ≥ 50 dB
根据阻带衰减要求，海明窗（-53dB）或布莱克曼窗（-74dB）都能满足。我个人习惯先试试海明窗，因为它的过渡带比布莱克曼窗窄，计算量也小一些。
滤波器阶数 N 怎么定？有个经验公式：N ≈ (A - 8) / (2.285 × Δf)，其中 A 是阻带衰减（dB），Δf 是归一化过渡带宽度。对于海明窗，Δf = (150-100)/1000 = 0.05，算下来 N ≈ (53-8)/(2.285×0.05) ≈ 394。嗯，取整到奇数，N=395。
下面是MATLAB代码示例：
% 参数设置
fs = 1000; % 采样频率
fp = 100; % 通带截止频率
fst = 150; % 阻带起始频率
fc = (fp + fst)/2; % 截止频率
N = 395; % 滤波器阶数
% 理想低通滤波器系数
n = 0:N-1;
hd = 2*fc/fs * sinc(2*fc/fs*(n - (N-1)/2));
% 海明窗
w = hamming(N)';
% 实际滤波器系数
h = hd .* w;
% 频率响应
[H, f] = freqz(h, 1, 1024, fs);
plot(f, 20*log10(abs(H)));
grid on;
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');
title('海明窗设计低通滤波器');4.4 四种窗函数的实际效果对比
咱们用上面的例子，分别用四种窗函数设计滤波器，看看效果。
| 窗函数 | 阶数 N | 过渡带宽度 (Hz) | 阻带衰减 (dB) | 通带波纹 (dB) |
|---|---|---|---|---|
| 矩形窗 | 395 | 约 4.5 | -21 | ±0.74 |
| 汉宁窗 | 395 | 约 15.5 | -44 | ±0.055 |
| 海明窗 | 395 | 约 16.5 | -53 | ±0.019 |
| 布莱克曼窗 | 395 | 约 27.5 | -74 | ±0.001 |
看到没？矩形窗过渡带最窄，但阻带衰减惨不忍睹。布莱克曼窗阻带衰减最好，但过渡带宽得吓人。说白了，这就是个“鱼和熊掌”的取舍问题。
4.5 如何选择窗函数？
我个人总结了一套选择策略，供你参考：
- 阻带衰减要求 < -30 dB： 别用矩形窗，至少用汉宁窗。
- 阻带衰减要求 -30 ~ -50 dB： 汉宁窗或海明窗都行。海明窗的旁瓣更低，但主瓣稍宽。
- 阻带衰减要求 -50 ~ -70 dB： 海明窗是首选，性价比最高。
- 阻带衰减要求 > -70 dB： 上布莱克曼窗，或者考虑凯泽窗（咱们后面章节会讲）。
你想想看，如果系统对带外抑制要求极高，比如雷达的脉冲压缩滤波器，那布莱克曼窗几乎是标配。但代价就是阶数要翻倍，FPGA资源消耗也大。所以，设计时一定要在性能和资源之间找到平衡点。
4.6 小结
今天咱们把窗函数法的四种经典窗过了一遍。记住几个关键点：
- 矩形窗过渡带最窄，但阻带衰减最差，慎用。
- 汉宁窗和海明窗是“万金油”，适合大多数场景。
- 布莱克曼窗阻带衰减最好，但阶数要求高。
- 选窗时，先看阻带衰减要求，再算阶数，最后评估资源。
嗯，下节课咱们聊聊凯泽窗，它有个可调的β参数，能让你在过渡带和阻带衰减之间自由调节，比这四种固定窗灵活多了。到时候咱们再细聊。