IIR滤波算法c语言

下面是一个简单的基于IIR滤波器的C语言实现：

#define N 3 //IIR滤波器的阶数

float a[N+1] = {1.0, -1.9, 0.95, -0.9025}; //IIR滤波器的反馈系数
float b[N+1] = {0.0483, 0.1449, 0.1449, 0.0483}; //IIR滤波器的前馈系数
float x[N+1] = {0.0}; //输入信号的缓冲区
float y[N+1] = {0.0}; //输出信号的缓冲区

void iir(float input, float *output)
{
    int i;

    //将当前输入值存储到缓冲区的第一个位置
    x[0] = input;
    y[0] = 0.0;

    //计算输出值
    for (i = 0; i <= N; i++) 
    {
        y[0] += b[i] * x[i];
    }

    for (i = 1; i <= N; i++) 
    {
        y[0] -= a[i] * y[i];
    }

    //将缓冲区中的值向后移动
    for (i = N; i > 0; i--)
    {
        x[i] = x[i - 1];
        y[i] = y[i - 1];
    }

    //将输出值保存到输出指针所指向的位置
    *output = y[0];
}

这个代码中，我们假设IIR滤波器的阶数为3，反馈系数存储在a数组中，前馈系数存储在b数组中。输入信号存储在x数组中，输出信号存储在y数组中。每次新的输入值存储在x[0]中。输出信号保存在y[0]变量中，并且向后移动缓冲区中的值，以便下一次的运算。

我们可以在程序中调用iir函数，传递输入值和输出指针作为参数。下面是一个示例：

int main()
{
    float input = 1.0;
    float output;

    iir(input, &output);

    return 0;
}

这个示例中，我们假设输入信号为1.0。我们调用iir函数，并将输出指针传递给它。当函数返回时，输出指针所指向的位置将保存输出信号的值。