c++实现python的filtfilt函数

filtfilt函数是一个数字信号处理函数，可以实现无相位延迟滤波器。以下是一个用C语言实现filtfilt函数的示例代码：

#include <stdio.h>
#include <stdlib.h>

void filtfilt(double *b, double *a, int nb, int na, double *x, int nx, double *y)
{
    int i, j;
    double *y1, *y2;
    y1 = (double *) malloc(sizeof(double) * nx);
    y2 = (double *) malloc(sizeof(double) * nx);
    // forward filtering
    y1[0] = b[0] * x[0];
    for (i = 1; i < nx; i++) {
        y1[i] = 0.0;
        for (j = 0; j < nb; j++) {
            if (i - j >= 0) {
                y1[i] += b[j] * x[i - j];
            }
        }
        for (j = 1; j < na; j++) {
            if (i - j >= 0) {
                y1[i] -= a[j] * y1[i - j];
            }
        }
    }
    // backward filtering
    y2[nx - 1] = y1[nx - 1];
    for (i = nx - 2; i >= 0; i--) {
        y2[i] = 0.0;
        for (j = 0; j < nb; j++) {
            if (i + j < nx) {
                y2[i] += b[j] * x[i + j];
            }
        }
        for (j = 1; j < na; j++) {
            if (i + j < nx) {
                y2[i] -= a[j] * y2[i + j];
            }
        }
    }
    // scaling
    for (i = 0; i < nx; i++) {
        y[i] = 2.0 * y2[i] - y1[i];
    }
    free(y1);
    free(y2);
}

int main()
{
    double b[] = {0.2, 0.2};
    double a[] = {1.0, -0.6};
    double x[] = {1.0, 2.0, 3.0, 4.0, 5.0};
    double y[5];
    filtfilt(b, a, 2, 2, x, 5, y);
    for (int i = 0; i < 5; i++) {
        printf("%f ", y[i]);
    }
    printf("\n");
    return 0;
}

该示例代码实现了一个双指数衰减低通滤波器，并对一个输入信号进行filtfilt滤波处理。可以根据需要修改滤波器系数和输入信号，以实现不同的滤波效果。