雷达 · 通信 · 电子战 / 基础精学 120
FIR系数量化误差在频域表现为误差系数的傅里叶和,其全频带幅度可由误差绝对值之和界定。
只量化系数,不考虑输入、乘法和累加器误差;舍入到 F 位二进制小数且无系数溢出。
δ h[k]:系数量化误差(无量纲)
F:小数位数(bit)
L:系数数量(tap)
Δ H:频率响应误差(无量纲)
令
代入响应定义得到
输入信号不参与此系数扰动恒等式。
因为单位圆上的复指数模为1,
![|\Delta H(\omega)|\leq\sum_k|\delta h[k]|.](http://mmo.508mallsys.com/image/loading/transparent.png)
最近舍入满足 |δ h[k]|≤2−F−1,因此
此结论是绝对幅度界,不是相对分贝界。
同一输入下
由三角不等式
当理想阻带接近零,极小绝对扰动仍可能导致很大相对误差,故应同时报告绝对误差和阻带底限。
设系数误差为 ek,则频响误差是误差序列的傅里叶变换。三角不等式给出逐频统一上界,而正交积分给出全频平均能量:
二者分别回答最大可能偏差和平均平方偏差,不应混作同一指标。
若有 M 个抽头且最近舍入步长为 Δ,逐个误差至多半步长,因此
平均平方界不能保证深阻带内的相对误差很小:理想响应接近零时,有限绝对误差对应任意大的相对比值。
因此通带可报告相对偏差,阻带仍应报告绝对幅度或衰减下限。
量化误差系数构成一只独立的误差滤波器。绝对误差界对所有频率有效,但通常比实际峰值保守。
对17抽头固定窗化低通分别使用6、10、14位小数;密集频率网格计算最大响应差,并比较实际绝对系数和界与通用长度界。
| 量 | 数值 |
|---|---|
| taps | 17 |
| F min | 6 |
| F max | 14 |
| measured error Q14 | 0.000128762022104 |
| l1 bound Q14 | 0.000202485805146 |
| generic bound Q14 | 0.000518798828125 |
| Q6 to Q14 ratio | 310.737224419 |
图 1 FIR系数量化频响误差界的理论关系与固定参数验证
表中数值和图中曲线由 MATLAB 实际运行得到,随机种子为 20261039;运行版本与完整结果另存为独立文件。
RTL实现误差滤波器 [1,−1,1] 的整数核,外部乘以系数误差量化单位即可形成指定误差通路;不等价于任意17抽头滤波器。 RTL 采用同步有效信号和异步低有效复位;
空拍保持状态且输出无效。除另有说明外,输入为16位有符号整数,结果为64位有符号整数,测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位;
未进行器件布局布线、时序收敛或功耗验证。
独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。
系数量化界可用于把滤波器阻带预算分配给系数小数位。例如阻带绝对允许误差为某阈值时,可先由长度界保守选择位数,再通过实际误差系数的频率扫描收紧估计。最终验收还应把数据通路舍入误差与系数误差分开测量。
网格最大值只是连续频率上界的下界,严格保证来自解析不等式。
保持对称系数配对可保留线性相位,但不保证原阻带指标。
运行时舍入噪声不包含在系数误差界内。
FIR系数量化误差在频域表现为误差系数的傅里叶和,其全频带幅度可由误差绝对值之和界定。
内容依据:每日精学第 120 课《FIR系数量化频响误差界》。图表沿用原课固定参数数值结果,属于数值验证,不代表设备实测。