C语言实现Savitzky - Golay滤波：系数用Matlab计算及开发注意要点

下面是一个使用C语言实现Savitzky - Golay滤波的程序示例，其中滤波系数假设已经通过Matlab计算得到。 ```c #include <stdio.h> // 定义Savitzky - Golay滤波函数 // input: 输入信号数组 // output: 输出滤波后的信号数组 // n: 输入信号的长度 // window_size: 滤波窗口的大小 // coefficients: Savitzky - Golay滤波系数数组 void savitzky_golay_filter(double *input, double *output, int n, int window_size, double *coefficients) { int half_window = window_size / 2; // 对每个数据点进行滤波 for (int i = 0; i < n; i++) { double sum = 0; // 计算当前点的滤波值 for (int j = -half_window; j <= half_window; j++) { int index = i + j; // 处理边界情况 if (index < 0) { index = 0; } else if (index >= n) { index = n - 1; } // 累加滤波系数与对应输入值的乘积 sum += input[index] * coefficients[j + half_window]; } // 将滤波结果存入输出数组 output[i] = sum; } } 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]); // 滤波窗口大小，假设为5 int window_size = 5; // 假设这些系数是通过Matlab计算得到的 double coefficients[] = {-3.0 / 35, 12.0 / 35, 17.0 / 35, 12.0 / 35, -3.0 / 35}; // 用于存储滤波后的信号 double output[10]; // 调用Savitzky - Golay滤波函数 savitzky_golay_filter(input, output, n, window_size, coefficients); // 输出滤波后的信号 for (int i = 0; i < n; i++) { printf("output[%d] = %f\n", i, output[i]); } return 0; } ``` 代码 1. `savitzky_golay_filter`函数实现了Savitzky - Golay滤波的核心逻辑。它接受输入信号数组、输出数组、信号长度、滤波窗口大小和滤波系数数组作为参数。 2. 在`main`函数中，我们定义了一个示例输入信号，指定了滤波窗口大小，并假设了一组通过Matlab计算得到的滤波系数。然后调用` savitzky_golay_filter`函数进行滤波，并输出滤波后的信号。 我在以往的开发过程中，碰到过以下几点需要格外小心处理的地方： 1. 滤波系数的计算：确保通过Matlab计算得到的滤波系数是正确的，不同的窗口大小和多项式阶数会得到不同的系数。 2. 边界处理：在处理边界数据时，需要根据具体需求选择合适的边界处理方法，这里采用的是简单的边界复制方法，可能不适用于所有情况。 3. 数组越界：在计算滤波值时，要确保不会访问到数组之外的元素，特别是在处理边界数据时要格外小心。 ######[AI写代码神器 | 831点数解答 | 2025-09-27 12:40:13]