如何解决C语言指数滤波器输出因浮点精度无法收敛到目标值的问题

问题根因 #

你遇到的现象是两个原因共同作用的结果：

• 指数滤波本身属于渐进收敛算法，理论上输出只会无限趋近输入值，永远不会严格等于输入
• IEEE 754 32位浮点数的精度有限，当输出与输入的差值小到低于当前数值对应的最小可表示步进时，(差值 * 收敛系数)的结果会小到与输入相加后，舍入结果和原输出完全一致，迭代就此卡住，数值再也不会变化。

另外你示例里的循环用Out != 0.0f作为终止条件属于典型的浮点数误用，工程代码中除非是判断手动赋值的精确0值，否则永远不要用浮点数的精确相等/不等判断。

最优解决方案（无额外逻辑开销） #

直接修改终止判断逻辑，用阈值替代精确相等判断，阈值可以根据你的应用精度要求选择：

// 例如要求输出小于1e-6即认为收敛到0
for (float Out = 1.0f; fabsf(Out) > 1e-6f;)
{
  Exp_Filt(0.0f, &Out, 0.1f, 1.0f);
}

这个方案不需要修改滤波函数本身，完全符合工程使用逻辑——你实际使用中根本不需要输出严格等于0，只要小到你的应用可接受的误差范围内即可。

轻量修改滤波函数方案（仅1条判断语句） #

如果你需要滤波函数本身能保证最终输出精确等于输入值，只需要加一行极低开销的判断，完全不属于你说的「大量条件语句」的情况：

#include <math.h>
void Exp_Filt(float In, float *Out, float time_base, float time_constant) {
  const float alpha = expf(-time_base / time_constant);
  *Out = In + ((*Out - In) * alpha);
  // 差值小于阈值时直接赋值为输入，阈值可按需调整
  if (fabsf(*Out - In) < 1e-7f) *Out = In;
}

这个判断只有在输出快收敛的时候才会触发，平时几乎没有额外开销。

无判断精度优化方案 #

如果连一条判断都不想加，可以用FMA（融合乘加）指令提升计算精度，大幅延后迭代卡住的阈值，大多数场景下足够让输出降到你需要的精度范围内：

#include <math.h>
void Exp_Filt(float In, float *Out, float time_base, float time_constant) {
  const float alpha = expf(-time_base / time_constant);
  // fmaf会把乘加合并为一条指令，中间不做舍入，精度比分开的乘加高得多
  *Out = fmaf(*Out - In, alpha, In);
}

目前主流的x86、ARM处理器都硬件支持FMA指令，运算速度和普通乘加完全一致，没有额外开销。

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