基于GSL的IIR滤波器实例

以下是一个基于GSL的IIR滤波器实例，该实例使用了GSL库中的函数来实现一个二阶IIR低通滤波器：

#include <stdio.h>
#include <gsl/gsl_filter.h>

int main()
{
    gsl_filter_iir_coeffs(&b, &a, 0.5, 0.0, 0.0, 0.5, 0.2, 0.1, 0.0);
    gsl_filter_iir_workspace* w = gsl_filter_iir_alloc(2);
    gsl_filter_iir_reset(w);
    double input[10] = {1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0};
    double output[10];
    gsl_filter_iir(w, b, a, 10, input, output);
    for(int i=0; i<10; i++) {
        printf("%f\n", output[i]);
    }
    gsl_filter_iir_free(w);
    return 0;
}

在上面的代码中，我们首先使用 `gsl_filter_iir_coeffs` 函数来计算IIR滤波器的系数，然后使用 `gsl_filter_iir_alloc` 函数来分配IIR滤波器的工作空间。接着我们使用 `gsl_filter_iir_reset` 函数来重置IIR滤波器的状态，然后使用 `gsl_filter_iir` 函数来对输入数据进行滤波，并将结果存储在输出数组中。最后，我们使用 `gsl_filter_iir_free` 函数来释放IIR滤波器的工作空间。

值得注意的是，上面的代码中的 `b` 和 `a` 数组是IIR滤波器的系数数组，可以通过调用 `gsl_filter_iir_coeffs` 函数来计算。其中，第一个参数是指向 `b` 数组的指针，第二个参数是指向 `a` 数组的指针，后面的参数是IIR滤波器的设计参数。在这个例子中，我们实现了一个二阶IIR低通滤波器，其截止频率为0.1，通带增益为0.5，阻带增益为0.2。

总之，使用GSL库中的函数来实现IIR滤波器非常简单，只需要计算系数、分配工作空间、重置状态、滤波输入数据即可。