fir低通滤波器设计原理.docxVIP

fir低通滤波器设计原理

FIR（有限长单位冲激响应）低通滤波器的设计核心是通过构造有限长的单位冲激响应序列，让滤波器只允许低于截止频率的信号通过，同时抑制高于截止频率的信号，且避免非线性相位失真（可设计为线性相位）。

一、核心设计原理：从“理想低通”到“实际FIR”

1.理想低通滤波器的特性（理论基础）

理想低通滤波器的频率响应是“矩形窗”：

通带：在截止频率?fc??以下，增益为1（信号无衰减通过）；

阻带：在?fc??以上，增益为0（信号完全被抑制）；

相位：线性相位（相位随频率线性变化，避免信号失真）。

但其单位冲激响应?hideal?(n)?是无限长的sinc函数（hideal?(n)=fs?2fc???sinc(fs?2fc?n?)，其中?fs??为采样频率），无法物理实现——因此FIR低通设计的本质是“用有限长序列逼近无限长的理想冲激响应”。

2.关键步骤：截断与窗函数（解决“无限长”问题）

直接截断无限长的sinc函数会导致吉布斯现象（通带和阻带出现波动，过渡带变宽），因此需通过“窗函数”优化截断过程：

步骤1：确定理想冲激响应根据目标截止频率?fc?、采样频率?fs?，计算无限长的理想sinc响应?hideal?(n)。

步骤2：选择窗函数并截断用一个有限长的“窗函数”?w(n)（如矩形窗、汉宁窗、汉明窗）与?hideal?(n)?相乘，得到有限长的实际FIR冲激响应：h(n)=hideal?(n)?w(n)窗函数的作用是“平滑截断”，减少吉布斯现象，同时平衡过渡带宽度和阻带衰减（不同窗函数特性不同，需按需选择）。

步骤3：保证线性相位（可选但重要）若需避免信号相位失真，需让FIR的冲激响应满足对称条件：

偶对称：h(n)=h(N?1?n)（N为滤波器阶数，即冲激响应长度，N=阶数+1），适用于低通、带通、高通滤波器；满足对称条件后，频率响应的相位为??(ω)=?2(N?1)ω?（线性相位，仅延迟无失真）。

二、主流设计方法（3种核心实现方式）

不同方法的核心差异在于“如何逼近理想响应”，需根据对“过渡带、阻带衰减、复杂度”的需求选择：

设计方法

核心逻辑

优点

缺点

窗函数法

用窗函数截断理想sinc响应，简单直接

计算简单、易实现，可保证线性相位

过渡带固定（由窗函数类型决定），阻带衰减有限（如矩形窗阻带衰减仅20dB，汉明窗约40dB）

频率采样法

直接在频域设定“通带、过渡带、阻带”的目标采样点，反推时域冲激响应

可灵活控制频域各点的幅度（如通带允许微小波动，优化阻带衰减）

易出现时域响应“边缘振荡”，需额外加窗抑制，线性相位实现较复杂

等波纹法（Parks-McClellan）

用“切比雪夫逼近理论”，让通带和阻带的波动幅度均匀最小化（等波纹）

过渡带最窄（相同阶数下），阻带衰减高，资源利用率最优

计算复杂（需迭代优化），需借助工具（如MATLAB、Python的scipy库）实现

三、关键参数选择（影响滤波器性能）

设计时需明确4个核心参数，直接决定FIR低通的效果：

截止频率?fc?需根据信号需求设定（如要保留1kHz以下信号，fc?=1kHz），通常以“归一化频率”（ωc?=fs?2πfc??）计算，避免受采样频率影响。

滤波器阶数NN越大（冲激响应越长），过渡带越窄、阻带衰减越高，但计算复杂度和延迟也越大（需平衡性能与成本）。

窗函数类型

若需窄过渡带（如实时信号处理）：选矩形窗（过渡带最窄，但阻带衰减低）；

若需高阻带衰减（如噪声抑制）：选布莱克曼窗（阻带衰减约74dB，但过渡带宽）；

兼顾两者：选汉明窗（阻带衰减40dB，过渡带中等，最常用）。

采样频率?fs?需满足奈奎斯特准则（fs?2fmax?，fmax??为信号最高频率），且?fc??需远小于?fs?/2，避免过渡带与Nyquist频率（fs?/2）重叠。

四、设计流程总结（以窗函数法为例）

明确需求：确定通带截止频率?fc?、阻带衰减要求、采样频率?fs?；

计算理想冲激响应：根据?fc??和?fs?，生成无限长的sinc序列?hideal?(n)；

选择窗函数：根据阻带衰减和过渡带需求，选汉明窗/汉宁窗等，确定窗长（即滤波器阶数N）；

生成FIR冲激响应：将?hideal?(n)?与窗函数相乘，得到有限长的?h(n)；

验证性能：通过傅里叶变换计算?h(n)?的频率响应，检查通带波动、过渡带宽度、阻带衰减是否满足需求；

实现滤波：将?h(n)?用于信号的卷积运算（时域滤波），或通过FFT实现频域滤波。