用C语言写一个抑制603Hz和1745Hz声音信号的滤波器，编程采用Linux风格

好的，以下是一个简单的抑制603Hz和1745Hz声音信号的IIR滤波器的C语言实现，采用Linux风格编程：

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

#define PI 3.14159265358979323846

// 603Hz和1745Hz的角频率
double w1 = 2 * PI * 603 / 44100;
double w2 = 2 * PI * 1745 / 44100;

// 滤波器系数
double b0 = 1.0;
double b1 = -2 * cos(w1);
double b2 = 1.0;
double a0 = 1.0;
double a1 = -2 * cos(w2);
double a2 = 1.0;

// IIR滤波器函数
double IIR_filter(double x)
{
    static double x1 = 0.0, x2 = 0.0, y1 = 0.0, y2 = 0.0;
    double y = (b0 / a0) * x + (b1 / a0) * x1 + (b2 / a0) * x2 - (a1 / a0) * y1 - (a2 / a0) * y2;
    x2 = x1;
    x1 = x;
    y2 = y1;
    y1 = y;
    return y;
}

int main()
{
    // 读入音频数据
    FILE *fp_in = fopen("input.pcm", "rb");
    if (fp_in == NULL) {
        printf("Failed to open input file!\n");
        return -1;
    }
    // 创建输出文件
    FILE *fp_out = fopen("output.pcm", "wb");
    if (fp_out == NULL) {
        printf("Failed to create output file!\n");
        return -1;
    }

    // 逐个样本进行滤波
    short sample;
    while (fread(&sample, sizeof(sample), 1, fp_in)) {
        double sample_double = (double)sample / 32768.0; // 转换为[-1,1]范围内的浮点数
        double filtered_sample_double = IIR_filter(sample_double);
        short filtered_sample = (short)(filtered_sample_double * 32768.0); // 转换回[-32768,32767]范围内的整数
        fwrite(&filtered_sample, sizeof(filtered_sample), 1, fp_out);
    }

    fclose(fp_in);
    fclose(fp_out);
    return 0;
}

注意，这只是一个简单的实现，实际应用中可能需要更复杂的滤波器设计和优化，以满足特定的要求。