一个设计低通巴特沃斯数字滤波器的实例

本人本科渣渣一个，前两天导师让我设计一个数字滤波器。由于本人基本没有数字信号处理基础，于是只能依靠百度和matlab，折腾了半天总算是摸索明白了。百度上有一些文章不靠谱，很容易误导别人，故在此发一篇博客。

滤波器设计目标：设计一个1Hz截止频率的2阶低通巴特沃斯数字滤波器，并转化成C语言函数。（国标里提的要求）

滤波器指标：指标：截止频率Fc = 1Hz，阶数N=2，低通巴特沃斯滤波器，采样频率Fs = 15Hz。

一、Matlab计算滤波器系数

Matlab计算巴特沃斯低通滤波器系数过程如下：

①根据给定的通带截止频率、通带截止增益、阻带截止频率、阻带截止增益，利用buttord函数计算巴特沃斯滤波器所需的最小阶数和截止频率。（由于要求中已给出截止频率，故这步省去）

（图是网上扒拉的，指标与本设计不符）

②根据上述计算得到的阶数，利用buttap函数计算出巴特沃斯滤波器原型。

③利用lp2lp函数，将原型滤波器转换成目标截止频率的滤波器。

④利用脉冲响应不变法（impinvar函数）或是双线性变换法（bilinear函数）将模拟滤波器转换为数字滤波器。数字滤波器形式为z域有理函数，分子分母系数即为滤波器系数。

我这里选用的是脉冲响应不变法，因为计算得到的滤波器比较简单，运算速度比较快。