基于GSL的FIR滤波器

GSL (GNU Scientific Library)是一个开源的科学计算库，其中包括了FIR滤波器的实现。下面是一个基于GSL的FIR滤波器的示例代码：

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

#define N 1000
#define L 51

int main(void)
{
    double data[N], output[N];
    gsl_filter_fir_workspace *w = gsl_filter_fir_alloc(L);
    gsl_filter_fir_init(w);

    // Generate input signal
    for (int i = 0; i < N; i++) {
        data[i] = sin(i * 0.1);
    }

    // Apply FIR filter
    gsl_filter_fir(GSL_FILTER_END_PADVALUE, 0, data, 1, output, 1, N, w);

    // Print output signal
    for (int i = 0; i < N; i++) {
        printf("%g\n", output[i]);
    }

    gsl_filter_fir_free(w);
    return 0;
}

这段代码演示了如何使用GSL实现一个FIR滤波器。前面的部分定义了滤波器的长度L和输入信号的长度N，以及分配了一个FIR滤波器所需的工作区间。接下来生成了一个简单的输入信号并初始化了FIR滤波器。最后，调用gsl_filter_fir函数进行滤波，并输出结果。