C语言 SG滤波平滑函数

以下是一个简单的 C 语言函数，可以实现 SG 滤波平滑：

void smooth(float *data, int n, int window_size, int degree) {
    int i, j, k;
    int half_window = window_size / 2;
    float sum, coeff;

    float *smoothed_data = malloc(n * sizeof(float));

    for (i = 0; i < n; i++) {
        sum = 0.0;
        for (j = -half_window; j <= half_window; j++) {
            if ((i + j >= 0) && (i + j < n)) {
                coeff = 0.0;
                for (k = 0; k <= degree; k++) {
                    coeff += pow(j, k);
                }
                sum += coeff * data[i + j];
            }
        }
        smoothed_data[i] = sum;
    }

    for (i = 0; i < n; i++) {
        data[i] = smoothed_data[i];
    }

    free(smoothed_data);
}

这个函数接受一个 float 类型的数组 data，数组大小为 n。window_size 参数定义了滤波窗口大小，degree 参数定义了滤波器的阶数。函数通过计算每个数据点的加权平均值来实现平滑。加权系数由 SG 滤波器的多项式计算而来，具体实现可以参考上面的代码。