matlab中filter滤波器fdatool转换为C语言

在命令区中输入fdatool,按enter键，出现下图Filter Design&Analysis Tool界面：

在其中选择滤波器类型（低通，带通，高通），选择IIR，设置好采样频率及通频之后，点击Design Filter.

点击Edit中Convert to Single Section，即得到filter函数中的b,a参数：

点击Targets中的Generate C header...即可得到b,a参数的一维数组，下面即是C语言代码

#include <stdio.h>
#include <math.h>
#include <memory.h>
#include <stdlib.h>
#define EPS 0.000001

 

double filter(const double* x, double* y, int xlen, double* a, double* b, int nfilt)//nfilt为系数数组长度
    {
        double tmp;
        int i, j;

        //normalizati