基于GSL的FIR滤波器实例

以下是一个基于GNU Scientific Library (GSL)的FIR滤波器实例：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_fft_complex.h>
#include <gsl/gsl_wavelet.h>

#define N 512
#define L 5

int main(int argc, char *argv[]) {
    double data[N];
    gsl_wavelet_workspace *work = gsl_wavelet_workspace_alloc(N);

    // 生成输入数据
    for (size_t i = 0; i < N; i++) {
        data[i] = sin(2 * M_PI * i / (double) N) + 0.1 * rand() / (double) RAND_MAX;
    }

    // 初始化低通滤波器
    gsl_wavelet *w = gsl_wavelet_alloc(gsl_wavelet_daubechies, L);
    gsl_wavelet_transform_forward(w, data, 1, N, work);

    // 提取低频分量
    for (size_t i = 0; i < N; i++) {
        if (i >= N / 2) {
            data[i] = 0.0;
        }
    }

    // 反向变换
    gsl_wavelet_transform_inverse(w, data, 1, N, work);

    gsl_wavelet_free(w);
    gsl_wavelet_workspace_free(work);

    // 输出结果
    for (size_t i = 0; i < N; i++) {
        printf("%f\n", data[i]);
    }

    return 0;
}

在这个例子中，我们使用Daubechies小波进行FIR滤波。我们首先生成一个长度为512的正弦信号，并添加一些随机噪声。然后我们初始化一个低通滤波器，并进行一次正向小波变换。接着，我们提取出低频分量，将高频分量置为0，并进行一次反向变换。最后，我们输出结果。

注意，在这个例子中，我们使用了GNU Scientific Library (GSL)提供的小波变换函数。GSL是一个广泛使用的数值库，提供了各种数值分析和科学计算函数。如果您想了解更多关于GSL的信息，请参考官方文档：https://www.gnu.org/software/gsl/