KNOWECON WIKI

ARTICLE

频率采样法

频率采样法 (Frequency Sampling Method) 频率采样法 (Frequency Sampling Method) 是数字信号处理中设计有限冲激响应 (FIR，Finite Impulse Response) 数字滤波器的一种经典方法。其核心思想是：在频域中对期望的频率响应进行采样，然后通过离散傅里叶逆变换 (IDFT) 求得滤波器的单位

浏览 0 更新 2025-11-09

频率采样法 (Frequency Sampling Method)

频率采样法 (Frequency Sampling Method) 是数字信号处理中设计有限冲激响应 (FIR，Finite Impulse Response) 数字滤波器的一种经典方法。其核心思想是：在频域中对期望的频率响应进行采样，然后通过离散傅里叶逆变换 (IDFT) 求得滤波器的单位冲激响应系数。该方法因直观、灵活而广泛应用于滤波器设计领域，尤其适用于设计具有任意幅频响应的 FIR 滤波器。

基本原理

频率采样法建立在离散傅里叶变换 (DFT) 的理论基础之上。对于一个长度为 N N 的 FIR 滤波器，其频率响应 H(ejω) H(e^{j\omega}) 与其单位冲激响应 h[n] h[n] 之间构成一对离散时间傅里叶变换 (DTFT) 关系：

H(ejω)=∑n=0N−1h[n]e−jωnH(e^{j\omega}) = \sum_{n=0}^{N-1} h[n] e^{-j\omega n}

频率采样法的基本策略是：先在 0 0 到 2π 2\pi 范围内等间隔地选取 N N 个频率点 ωk=2πkN(k=0,1,…,N−1) \omega_k = \frac{2\pi k}{N} (k = 0, 1, \dots, N-1) ，在这些频率点上指定期望的频率响应值 Hd[k]=Hd(ejωk) H_d[k] = H_d(e^{j\omega_k}) ，然后通过 N N 点 IDFT 求得滤波器系数：

h[n]=1N∑k=0N−1Hd[k]ej2πkn/N,n=0,1,…,N−1h[n] = \frac{1}{N} \sum_{k=0}^{N-1} H_d[k] e^{j 2\pi k n / N}, \quad n = 0, 1, \dots, N-1

由此得到的 h[n] h[n] 所对应的实际频率响应 H(ejω) H(e^{j\omega}) 在采样点上与期望值精确一致，即 H(ejωk)=Hd[k] H(e^{j\omega_k}) = H_d[k] ，而在采样点之间的频率响应则由内插公式决定。

数学表述与内插公式

给定 N N 个频率采样值 H[k](k=0,1,…,N−1) H[k] (k = 0, 1, \dots, N-1) ，滤波器的系统函数 H(z) H(z) 可以表示为：

H(z)=∑n=0N−1h[n]z−nH(z) = \sum_{n=0}^{N-1} h[n] z^{-n}

利用 DFT 的性质，可以推导出用频率采样值直接表示的系统函数：

H(z)=1−z−NN∑k=0N−1H[k]1−ej2πk/Nz−1H(z) = \frac{1 - z^{-N}}{N} \sum_{k=0}^{N-1} \frac{H[k]}{1 - e^{j 2\pi k / N} z^{-1}}

由此可得频率响应的内插公式：

H(ejω)=e−jω(N−1)/2N∑k=0N−1H[k]sin⁡[N(ω−2πk/N)/2]sin⁡[(ω−2πk/N)/2]H(e^{j\omega}) = \frac{e^{-j\omega (N-1)/2}}{N} \sum_{k=0}^{N-1} H[k] \frac{\sin[N(\omega - 2\pi k / N) / 2]}{\sin[(\omega - 2\pi k / N) / 2]}

内插公式揭示了一个重要性质：在采样点 ωk \omega_k 处，分子和分母同时趋于零，利用洛必达法则可得 H(ejωk)=H[k] H(e^{j\omega_k}) = H[k] ，即频率响应精确通过采样点；而在采样点之间，频率响应由所有采样值共同加权内插决定。这意味着采样点之间的响应并非完全可控，可能出现较大的波动——这正是频率采样法的主要局限性所在。

线性相位约束与采样值设计

在实际应用中，FIR 滤波器通常要求具有线性相位特性，以保证信号通过滤波器后的相位不失真。对于长度为 N N 、满足偶对称或奇对称条件的线性相位 FIR 滤波器，其频率采样值 H[k] H[k] 不能任意指定，而必须满足特定的幅度和相位约束。

类型 I：N N 为奇数，h[n] h[n] 偶对称

此时频率响应可写为：

H(ejω)=A(ω)e−jω(N−1)/2H(e^{j\omega}) = A(\omega) e^{-j\omega (N-1)/2}

其中 A(ω) A(\omega) 为实值幅度函数，满足 A(ω)=A(2π−ω) A(\omega) = A(2\pi - \omega) 。相应地，频率采样值 H[k] H[k] 的幅度 ∣H[k]∣ |H[k]| 和相位 θ[k] \theta[k] 必须满足：

∣H[k]∣=∣H[N−k]∣,θ[k]=−N−1Nπk|H[k]| = |H[N-k]|, \quad \theta[k] = -\frac{N-1}{N} \pi k

设计时只需要指定前 (N+1)/2 (N+1)/2 个采样点的幅度值，其余采样点由对称关系自动确定。

类型 II：N N 为偶数，h[n] h[n] 偶对称

此类滤波器在 ω=π \omega = \pi 处幅度恒为零，因此在设计高通或带阻滤波器时需注意这一固有约束。相位约束与类型 I 相同。

通过合理设置频率采样点的幅度值，可以设计出低通、高通、带通、带阻等各类线性相位 FIR 滤波器。

过渡带采样点的优化

频率采样法的一个关键缺陷是：在通带截止频率和阻带起始频率之间存在较大的波动，即吉布斯现象 (Gibbs 现象)，导致阻带衰减不够理想。为改善这一状况，可以在过渡带中引入一个或多个非零的过渡采样点，以平滑滤波器的频率响应。

具体做法是：在理想频率响应的通带和阻带之间的过渡带放置一个或多个采样点，其幅度值在 0 和 1 之间自由调节，通过优化算法（如线性规划、梯度下降或网格搜索）寻找使阻带衰减最大化的最优过渡带采样值。常见的策略如下：

| 过渡带采样点数 | 近似阻带衰减 | |:---:|:---:| | 0 | -20 dB | | 1 | -40 dB 至 -50 dB | | 2 | -60 dB 至 -70 dB | | 3 | -80 dB 至 -90 dB |

增加过渡带采样点可以显著提升阻带衰减性能，但代价是增加了从通带到阻带的过渡带宽。设计者需要在过渡带宽和阻带衰减之间做出折中。

设计步骤

频率采样法设计 FIR 滤波器的标准流程如下：

  1. 确定滤波器规格：选定滤波器类型（低通、高通、带通、带阻），明确通带截止频率 ωp \omega_p 、阻带起始频率 ωs \omega_s 、通带最大衰减 Rp R_p 和阻带最小衰减 As A_s 。
  1. 选择滤波器长度 N N ：根据过渡带宽 Δω=ωs−ωp \Delta\omega = \omega_s - \omega_p 和阻带衰减要求，估算所需的滤波器阶数。一般而言，过渡带越窄，所需 N N 越大。
  1. 构造频率采样序列：在 N N 个频率采样点上，按照期望的幅频特性赋予幅度值（通带为 1，阻带为 0），并附加线性相位因子。如有必要，在过渡带放置优化采样点。
  1. 计算 IDFT：对构造的频率采样序列 H[k] H[k] 进行 N N 点 IDFT，得到滤波器系数 h[n] h[n] 。
  1. 验证滤波器性能：计算所得滤波器的实际频率响应，检查通带波动、阻带衰减和过渡带宽是否满足指标。若不满足，调整 N N 或过渡带采样值后重新设计。
  1. 实现滤波器：将求得的系数 h[n] h[n] 部署到目标硬件或软件平台。

与其他 FIR 设计方法的比较

频率采样法与窗函数法 (窗函数法) 和最优等波纹法 (Parks-McClellan 算法) 并列为 FIR 滤波器的三大经典设计方法，三者各有优劣：

| 特性 | 频率采样法 | 窗函数法 | Parks-McClellan | |:---|:---|:---|:---| | 设计灵活性 | 高（任意幅频响应） | 中 | 中（需分段常数） | | 最优性 | 非最优 | 非最优 | 最优（等波纹） | | 计算复杂度 | 低（一个 IDFT） | 低 | 高（迭代） | | 过渡带控制 | 间接（依赖 N N 和过渡采样） | 依赖窗函数类型 | 精确可控 | | 阻带衰减 | 可通过过渡采样大幅改善 | 受窗函数旁瓣限制 | 最佳 |

应用场景

频率采样法因其直观性和灵活性，在以下场景中尤为实用：

  • 任意幅频响应滤波器：当期望的频率响应不是标准形状（如多带滤波器、任意形状均衡器）时，频率采样法只需在采样点处指定增益即可，操作简便。
  • 快速原型设计：在设计初期需要快速验证滤波器概念时，频率采样法计算量小（仅需一次 IDFT），可迅速生成滤波器系数。
  • 频域直接设计：在某些应用中，设计者只有频域的采样数据（如从测量或仿真得到的期望频谱），频率采样法可以直接利用这些数据构造滤波器。
  • 联合时频域控制：频率采样法允许同时在频域（通过频率采样值）和时域（通过选择 N N 和冲激响应对称性）对滤波器施加约束。

注意事项与局限性

尽管频率采样法直观且易于编程实现，使用时需注意以下局限：

  1. 采样点间不可控：频率响应在非采样点处的行为完全由内插公式决定，采样点间可能出现超出预期的波动或陷波。
  2. 边缘频率对齐问题：当期望的截止频率与采样点位置不对齐时，实际截止频率会向最近的采样点偏移，需要增大 N N 以增加频率采样密度。
  3. 复系数问题：对于非对称的频率响应，IDFT 可能产生复数滤波器系数。实际应用中通常需强制系数的对称性以保证线性相位和实系数。
  4. 吉布斯效应：在频率响应跳变处，即使采样点精确匹配，采样点之间仍会出现由频域截断引起的过冲和波动，需要借助过渡带采样点加以抑制。

总体而言，频率采样法是 FIR 滤波器设计中一种简洁而有力的工具，在理解 DFT 与滤波器设计之间的内在联系、快速原型验证以及任意响应滤波器设计等场合具有不可替代的价值。