C语言生成50HZ陷波滤波器代码

以下是一个简单的C语言程序，用于生成50Hz陷波滤波器代码：

#include <stdio.h>
#include <math.h>

#define PI 3.14159265358979323846

int main()
{
    double f0 = 50.0; // 陷波中心频率
    double Q = 10.0; // 陷波品质因数
    double fs = 1000.0; // 采样率
    double w0 = 2.0 * PI * f0 / fs; // 角频率
    double alpha = sin(w0) / (2.0 * Q); // 滤波器增益

    double b0 = 1.0;
    double b1 = -2.0 * cos(w0);
    double b2 = 1.0;
    double a0 = 1.0 + alpha;
    double a1 = -2.0 * cos(w0);
    double a2 = 1.0 - alpha;

    printf("陷波滤波器系数：\n");
    printf("b0 = %f\n", b0 / a0);
    printf("b1 = %f\n", b1 / a0);
    printf("b2 = %f\n", b2 / a0);
    printf("a1 = %f\n", a1 / a0);
    printf("a2 = %f\n", a2 / a0);

    return 0;
}

这个程序使用了公式：

w0 = 2 * pi * f0 / fs
alpha = sin(w0) / (2 * Q)
b0 = 1
b1 = -2 * cos(w0)
b2 = 1
a0 = 1 + alpha
a1 = -2 * cos(w0)
a2 = 1 - alpha

其中，`f0` 是陷波中心频率，`Q` 是品质因数，`fs` 是采样率，`w0` 是角频率，`alpha` 是滤波器增益，`b0`、`b1`、`b2`、`a0`、`a1`、`a2` 是滤波器系数。

注意：这里的滤波器是IIR（无限脉冲响应）滤波器，需要注意滤波器稳定性和溢出问题。