如何在C语言中实现IIR滤波器？关于y[n-1]实现的困惑求助

看起来你卡在了IIR滤波器最关键的部分——维护历史输出的状态缓存上，这很正常，毕竟FIR只需要管理输入的历史数据，而IIR还要依赖之前的输出结果。先帮你梳理下现有代码的问题，再给出可运行的修正方案。

现有代码的核心问题 #

你的滤波器公式是对的，但处理y[n-1]相关的逻辑完全错误：

• 你直接操作*output--，但output数组初始阶段都是未初始化的垃圾值，就算后续有计算结果，这个指针操作也没有正确指向之前生成的y[n-1], y[n-2]。
• 没有专门的缓存保存历史输出值，IIR滤波器需要持续记录最近的N个输出（对应系数A的长度），供下一次计算使用，就像你用insamp保存历史输入一样。

修正思路：添加输出延迟线（状态缓存） #

我们需要一个专门的数组（比如叫delay_line_y）来保存最近的filterLength-1个输出值（因为系数A是从a1开始的，对应y[n-1]到y[n-(filterLength-1)]）。每次计算完当前的y[n]后，要更新这个延迟线：把旧的历史值往前移，新的y[n]放到末尾。

另外要注意：你的公式里默认a0=1（也就是公式是y(n) = b0x(n)+...+bM-1x(n-M+1) - a1y(n-1)-...-aN y(n-N)），如果你的系数A里包含a0，那还要把最终结果除以a0，这点要根据实际系数调整。

修正后的代码 #

#include <string.h> // 用于memmove操作

void IIRFloat(double *coeffs_B, double *coeffs_A, double *input, double *output, int length, int filterLength) {
    double b_acc = 0.0;
    double a_acc = 0.0;
    // 输出延迟线：保存最近的filterLength-1个输出值（y[n-1], y[n-2], ..., y[n-(filterLength-1)]）
    double delay_line_y[filterLength-1];
    // 初始化延迟线为0（系统初始静止，无历史输出）
    memset(delay_line_y, 0, sizeof(delay_line_y));

    for (int n = 0; n < length; n++) {
        // ---------------- 计算B部分：b0*x[n] + b1*x[n-1] + ... + bM-1*x[n-M+1] ----------------
        b_acc = 0.0;
        double *inputp = &input[n]; // 当前输入是input[n]，往前取历史输入
        for (int k = 0; k < filterLength; k++) {
            // 前filterLength个样本中，n-k会小于0，对应历史输入取0
            if (n - k >= 0) {
                b_acc += coeffs_B[k] * inputp[-k];
            }
            // 如果你的insamp是维护好的滑动窗口缓存，可替换为：
            // b_acc += coeffs_B[k] * insamp[filterLength-1 + n - k];
        }

        // ---------------- 计算A部分：a1*y[n-1] + a2*y[n-2] + ... + aN*y[n-N] ----------------
        a_acc = 0.0;
        for (int k = 1; k < filterLength; k++) {
            // 延迟线第k-1位对应y[n-k]（delay_line_y[0]是y[n-1]，以此类推）
            a_acc += coeffs_A[k] * delay_line_y[k-1];
        }

        // 计算当前输出y[n]
        output[n] = b_acc - a_acc;
        // 如果系数A包含a0（公式为a0*y[n] = ...），需改为：
        // output[n] = (b_acc - a_acc) / coeffs_A[0];

        // ---------------- 更新输出延迟线：将新输出加入，旧历史值左移 ----------------
        memmove(delay_line_y, delay_line_y + 1, sizeof(delay_line_y) - sizeof(double));
        // 把当前输出放到延迟线末尾，作为下一次计算的y[n-1]
        delay_line_y[filterLength-2] = output[n];
    }
}

关键细节解释 #

1. 延迟线初始化：系统初始静止，没有历史输出，所以把延迟线全部设为0，符合实际滤波场景的初始条件。
2. 输入历史处理：如果input数组只包含当前输入序列，前filterLength个样本的历史输入（n-k<0时）取0；如果你的insamp是已经维护好的滑动窗口缓存，可以替换成你原来的访问逻辑。
3. 延迟线更新：每次计算完y[n]后，用memmove将延迟线元素左移，再把新输出放到末尾，确保下一次计算时能拿到正确的历史输出值。
4. 系数A的a0处理：如果滤波器公式是a0y(n) + a1y(n-1) + ... = b0x(n) + ...，需要将最终结果除以coeffs_A[0]，代码中已添加注释提示。

额外建议 #

• 输入的历史缓存（比如你提到的insamp）也要用类似延迟线的方式维护，确保每次计算能正确获取x[n], x[n-1], ..., x[n-M+1]。
• 如果滤波器长度固定，可以将延迟线改为静态数组（或动态分配内存），避免每次调用函数都重新初始化，但静态数组要注意线程安全问题。

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