每日精学 120|FIR系数量化频响误差界

2026-09-25

雷达 · 通信 · 电子战 / 基础精学 120

FIR系数量化误差在频域表现为误差系数的傅里叶和,其全频带幅度可由误差绝对值之和界定。

01
模型与符号

只量化系数,不考虑输入、乘法和累加器误差;舍入到 F 位二进制小数且无系数溢出。

δ h[k]:系数量化误差(无量纲)

F:小数位数(bit)

L:系数数量(tap)

Δ H:频率响应误差(无量纲)

02
从模型到公式

误差分离

令

h_q[k]=h[k]+\delta h[k]

代入响应定义得到

H_q(\omega)-H(\omega)=\sum_{k=0}^{L-1}\delta h[k]e^{-\mathrm j\omega k}=\Delta H(\omega).

输入信号不参与此系数扰动恒等式。


一致误差界

因为单位圆上的复指数模为1,

|\Delta H(\omega)|\leq\sum_k|\delta h[k]|.

最近舍入满足 |δ h[k]|≤2−F−1,因此

\sup_\omega|\Delta H(\omega)|\leq L2^{-F-1}.

此结论是绝对幅度界,不是相对分贝界。


输出误差

同一输入下

\delta y[n]=\sum_k\delta h[k]x[n-k]

由三角不等式

|\delta y[n]|\leq X\sum_k|\delta h[k]|.

当理想阻带接近零,极小绝对扰动仍可能导致很大相对误差,故应同时报告绝对误差和阻带底限。


确定性误差界与平均误差的区别

设系数误差为 ek,则频响误差是误差序列的傅里叶变换。三角不等式给出逐频统一上界,而正交积分给出全频平均能量:

|\Delta H(e^{\mathrm j\omega})|\le\sum_k|e_k|, \qquad \frac1{2\pi}\int_{-\pi}^{\pi}|\Delta H|^2\,d\omega=\sum_k|e_k|^2.

二者分别回答最大可能偏差和平均平方偏差,不应混作同一指标。

若有 M 个抽头且最近舍入步长为 Δ,逐个误差至多半步长,因此

\|\Delta H\|_\infty\le M\Delta/2, \qquad \frac1{2\pi}\int|\Delta H|^2d\omega\le M\Delta^2/4.

平均平方界不能保证深阻带内的相对误差很小:理想响应接近零时,有限绝对误差对应任意大的相对比值。

因此通带可报告相对偏差,阻带仍应报告绝对幅度或衰减下限。

03
物理含义

量化误差系数构成一只独立的误差滤波器。绝对误差界对所有频率有效,但通常比实际峰值保守。

04
固定参数算例

对17抽头固定窗化低通分别使用6、10、14位小数;密集频率网格计算最大响应差,并比较实际绝对系数和界与通用长度界。

固定参数数值验证结果
量数值
taps17
F min6
F max14
measured error Q140.000128762022104
l1 bound Q140.000202485805146
generic bound Q140.000518798828125
Q6 to Q14 ratio310.737224419
图 1:FIR系数量化频响误差界的理论关系与固定参数验证图 1:FIR系数量化频响误差界的理论关系与固定参数验证

图 1 FIR系数量化频响误差界的理论关系与固定参数验证

表中数值和图中曲线由 MATLAB 实际运行得到,随机种子为 20261039;运行版本与完整结果另存为独立文件。

05
数字实现的边界

RTL实现误差滤波器 [1,−1,1] 的整数核,外部乘以系数误差量化单位即可形成指定误差通路;不等价于任意17抽头滤波器。 RTL 采用同步有效信号和异步低有效复位;

空拍保持状态且输出无效。除另有说明外,输入为16位有符号整数,结果为64位有符号整数,测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位;

未进行器件布局布线、时序收敛或功耗验证。

独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。

06
适用范围

系数量化界可用于把滤波器阻带预算分配给系数小数位。例如阻带绝对允许误差为某阈值时,可先由长度界保守选择位数,再通过实际误差系数的频率扫描收紧估计。最终验收还应把数据通路舍入误差与系数误差分开测量。

07
工程上需要注意

网格最大值只是连续频率上界的下界,严格保证来自解析不等式。

保持对称系数配对可保留线性相位,但不保证原阻带指标。

运行时舍入噪声不包含在系数误差界内。

FIR系数量化误差在频域表现为误差系数的傅里叶和,其全频带幅度可由误差绝对值之和界定。

内容依据:每日精学第 120 课《FIR系数量化频响误差界》。图表沿用原课固定参数数值结果,属于数值验证,不代表设备实测。

阅读8
分享
写评论...