matlab filtfilt c,matlab filter与filtfilt函数实现，C语言实现

嗯，算法非常简单，就是网上搜不到C代码实现。filter是个很万能的数字滤波器函数，只要有滤波器的差分方程系数，IIR呀FIR呀都能通过它实现。在MATLAB里面，filter最常用的格式是这两个：

[y,zf] = filter(b,a,X)

[y,zf] = filter(b,a,X,zi)

其中b和a就是差分方程的系数，X是输入向量，zi是“初始状态”。可能这么说明还是不很清晰，那么请看图(注意，a(1)为1，这个可以通过差分方程所有系数除以a(1)所得)：

嗯，这样子很轻松地就能把各个y值给算出来了，哦注意上面式子里面的n是“向量a或者b中最长的长度”，在实际编程的时候，如果a和b长度不一样，短者显然是要用0补齐的。对于那个初始状态zi，忽略的时候，比如(顺便提醒一句MATLAB的下标从1开始)

y(1)=b(1)x(1)

y(2)=b(1)x(2)+Z1(1)= b(1)x(2) + b(2)x(1) - a(2)y(1)

不忽略的时候

y(1)=b(1)x(1) + Z1(0)

y(2)=b(1)x(2)+Z1(1)= b(1)x(2) + b(2)x(1) - a(2)y(1) + Z2(0)

因为实际程序中自己定义的东西比较多(=,=|||这也是没办法的事情不是)，而filtfilt这个超级无敌的“零相移滤波函数”更是复杂到稍微调用了一下自己写的矩阵运算函数，所以代码全部贴上来实在是太乱。等这次工程做完会大概地把代码打包发一下。现在就先贴点代码片段好了……但愿我的代码风格没那么难懂……