c++实现python中的filtfilt函数

filtfilt函数用于对信号进行前向和后向滤波，以消除滤波器响应对信号的影响。以下是一个用C实现filtfilt函数的示例代码：

#include <stdio.h>

void filtfilt(double *b, double *a, int nb, int na, double *x, int nx, double *y) {
    int i, j;
    double z[nx + na - 1];

    // Forward filter
    for (i = 0; i < nx + na - 1; i++) {
        z[i] = 0;
        for (j = 0; j < na; j++) {
            if (i >= j && i - j < nx) {
                z[i] += a[j] * x[i - j];
            }
        }
        for (j = 0; j < nb; j++) {
            if (i >= j && i - j < nx) {
                z[i] += b[j] * y[i - j];
            }
        }
    }

    // Backward filter
    for (i = nx + na - 2; i >= 0; i--) {
        y[i] = 0;
        for (j = 0; j < na; j++) {
            if (i + j < nx + na - 1) {
                y[i] -= a[j] * z[i + j];
            }
        }
        for (j = 0; j < nb; j++) {
            if (i + j < nx + na - 1) {
                y[i] += b[j] * y[i + j];
            }
        }
    }
}

int main() {
    double b[] = {1, 1, 1};
    double a[] = {1, 0.5};
    double x[] = {1, 2, 3, 4, 5};
    double y[5];
    int nb = sizeof(b) / sizeof(double);
    int na = sizeof(a) / sizeof(double);
    int nx = sizeof(x) / sizeof(double);
    filtfilt(b, a, nb, na, x, nx, y);
    for (int i = 0; i < nx; i++) {
        printf("%f ", y[i]);
    }
    printf("\n");
    return 0;
}

该代码实现了一个简单的filtfilt函数，用于对输入信号x进行前向和后向滤波，并返回输出信号y。其中b和a是低通滤波器的系数，nb和na是系数的长度，x是输入信号，nx是信号的长度，y是输出信号。在代码中，先对输入信号进行前向滤波，然后再对前向滤波的结果进行后向滤波，最终得到输出信号。