C语言实现三阶Savitzky - Golay滤波器：代码详解与开发注意要点

```c #include <stdio.h> #include <stdlib.h> // 三阶Savitzky - Golay滤波器函数 void savitzky_golay_filter(double *input, double *output, int n, int window_size) { // 确保窗口大小是奇数 if (window_size % 2 == 0) { printf("窗口大小必须是奇数。\n"); return; } // 计算窗口的半宽 int half_window = window_size / 2; // 处理边界情况，这里简单采用复制边界值的方法 for (int i = 0; i < half_window; i++) { output[i] = input[i]; } // 计算滤波器系数 double coefficients[window_size]; for (int i = 0; i < window_size; i++) { int k = i - half_window; coefficients[i] = (3 * (3 * window_size * window_size - 7) - 30 * k * k) / (4 * window_size * (window_size * window_size - 4)); } // 应用滤波器 for (int i = half_window; i < n - half_window; i++) { output[i] = 0.0; for (int j = 0; j < window_size; j++) { output[i] += coefficients[j] * input[i - half_window + j]; } } // 处理另一侧边界情况 for (int i = n - half_window; i < n; i++) { output[i] = input[i]; } } int main() { // 示例输入数据 double input[] = {1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0}; int n = sizeof(input) / sizeof(input[0]); double *output = (double *)malloc(n * sizeof(double)); // 窗口大小 int window_size = 5; // 应用Savitzky - Golay滤波器 savitzky_golay_filter(input, output, n, window_size); // 输出滤波后的数据 for (int i = 0; i < n; i++) { printf("output[%d] = %f\n", i, output[i]); } // 释放内存 free(output); return 0; } ``` 这段C语言代码实现了一个三阶Savitzky - Golay滤波器。`savitzky_golay_filter`函数接收输入数组、输出数组、数组长度和窗口大小作为参数，在函数内部首先检查窗口大小是否为奇数，然后计算窗口半宽，接着计算滤波器系数，最后将滤波器应用到输入数据上。在`main`函数中，我们创建了示例输入数据，调用滤波函数并输出滤波后的结果。 我在以往的开发过程中，碰到过以下几点需要格外小心处理的地方： 1. 窗口大小必须为奇数，否则滤波器系数的计算会出错。 2. 边界处理很重要，如果不恰当处理边界，可能会导致滤波结果在边界处出现异常。 3. 要注意内存的分配和释放，避免出现内存泄漏的问题。 ######[AI写代码神器 | 764点数解答 | 2025-08-13 10:02:46]