Savitzky-Golay滤波器（2）
前几天写过一篇介绍 Savitzky-Golay滤波器的文章, 没想到最近做项目还真的用上了。
因此就顺便写了个 C 语言的自动计算生成 SG 滤波器系数的程序。利用这里的代码可以生成任意阶数的 SG 滤波器。对于一些需要动态调整 SG 滤波器系数的场合，这里的代码还是很实用的。 上次的文章中给出了计算 SG 滤波器系数的 matlab 代码。

/** * 计算 Savitzky-Golay 滤波器系数 * SG 滤波器的阶数为 2M+1，多项式的最高次数为 N */gsl_vector* SG_FilterCreate(int M, int N /* Poly Order */){    int i;    gsl_vector *x = gsl_vector_alloc(2 * M + 1);    gsl_vector *y = gsl_vector_alloc(2 * M + 1);    gsl_vector_set_zero(y);    gsl_vector_set(y, M, 1);    for(i= -M; i <= M; i++)    {        gsl_vector_set(x, i + M, i);    }    gsl_vector *c = PolyFit(x, y, N);    gsl_vector_free(x);    gsl_vector_free(y);    gsl_vector *fir = gsl_vector_alloc(2 * M + 1);    for(i = -M; i <= M; i++)    {        gsl_vector_set(fir, i + M, gsl_poly_eval (c->data, N + 1, i));    }    gsl_vector_free(c);    return fir;}