Matlab中滤波器系数量化引发异常幅值衰减的原因咨询

异常大幅衰减的来源 #

你在filterDesigner中看到的图2幅值响应，是未量化双精度浮点系数对应的理论响应，和你实际调用filter(Hd,y)时运行的量化后滤波器的真实响应完全不一致，衰减来自两个非常典型的系数量化错误：

• 滤波器结构选型问题：如果导出的是直接型结构的IIR滤波器，本身系数量化敏感度极高，尤其是量化字长设置较窄时，极点位置会发生严重偏移，直接导致通带内频点增益大幅下跌；换成二阶节（SOS）级联结构可以大幅降低量化敏感度。
• 量化定标错误：要么是量化字长给得过窄（比如小于8bit），留给系数的小数位精度不足，实际量化后的系数和理论设计值偏差过大；要么是SOS结构的各级缩放因子（Scale Value）在量化时被错误截断、饱和，导致通带总增益大幅降低。
  直接在命令行执行fvtool(Hd)就能看到量化后滤波器的真实幅频响应，200Hz位置的增益必然远低于0dB，和你看到的图2响应完全不符。

图1 滤波效果
图2 幅值响应

输出幅值恰好为0.52186的成因 #

这个数值不是随机误差，是固定量化规则下的确定计算结果：
你当前的量化配置（大概率是3~4bit的极窄字长系数量化，或是滤波运算过程中乘加结果错误截断了高位），会让200Hz频点对应的滤波器复数增益的模值固定为0.52186 / 1000 ≈ 0.00052186。
典型复现场景：500Hz截止、2kHz采样率的2阶巴特沃斯直接型低通滤波器，若将系数量化为3bit有符号定点数（可表示值仅为-1、-0.5、0、0.5四个），代入频响计算公式计算200Hz（归一化频率0.2π）处的增益，结果和该值几乎完全一致。这个数值是所有量化后的系数在200Hz频点的幅频贡献叠加得到的固定值，只要滤波器结构、量化配置不变，这个值就不会发生变化。

附：问题提供的测试代码 #

N = 2000;
f0=200;
f1=800;
fs=2000;
n=1:N;
y =1000*sin(2*pi*f0.*n/fs)+500*sin(2*pi*f1.*n/fs);
Hd = filter_below500Hz;
d = filter(Hd,y);
figure(2)
subplot(2,1,1)
plot(y);
title({['y，before filtering']});
subplot(2,1,2)
plot(d);
title({['d，after filtering']});

快速排查步骤 #

• 执行fvtool(Hd)查看量化后滤波器200Hz处的实际增益，先区分故障来自系数量化错误还是运算过程定标错误
• 将滤波器结构修改为二阶节级联形式，重新开展系数量化，给系数分配至少12bit以上的小数位精度
• 核对SOS结构各级缩放因子的量化结果，确认通带总增益和设计值一致

内容的提问来源于stack exchange，提问作者xxxfpga_so_difficult