MATLAB生成巴特沃斯滤波器+C语言实现之二（C语言篇）

上一个文章具体介绍了怎么用MATLAB生成滤波器参数，这些参数在C语言中使用可实现滤波器设计效果。这里就主要介绍这些参数怎么在C语言中使用。
使用之前得到的参数[b,a]：
b =
1.0e-03 *
0.0231 0.1388 0.3470 0.4627 0.3470 0.1388 0.0231
a =
1.0000 -4.5450 8.7464 -9.0984 5.3862 -1.7181 0.2305
在C语言中，首先声明数组，存放这些参数：
float b[7]={0.0000231,0.0001388,0.0003470,0.0004627,0.0003470,0.0001388,0.0000231};
float a[7]={1,-4.5450,8.7464,-9.0984,5.3862,-1.7181,0.2305};
再声明另两个数组：
float xBuf[7]={0,0,0,0,0,0,0};
float yBuf[7]={0,0,0,0,0,0,0};
然后使用函数：
float IIRLowPass(float x)
{
int i;
//运算之前Buf向前移动一个位置，以保存之前Buf的数据；
for(i=6; i>0; i–)
{
yBuf[i] = yBuf[i-1];
xBuf[i] = xBuf[i-1];
}
xBuf[0] = x;
yBuf[0] = 0;
for(i=1;i<7;i++)
{
yBuf[0] = yBuf[0] + b2[i]*xBuf[i];
yBuf[0] = yBuf[0] - a2[i]*yBuf[i];
}
yBuf[0] = yBuf[0] + b2[0]*xBuf[0];
return yBuf[0];
}

最终实现C语言滤波器设计。