IIR数字滤波器设计与滤波

典型的模拟低通滤波器的指标有两种方法：

    第一种：Ωp,Ωs分别为通带频率和阻带频率，δp,δs分别为通带容限和阻带容限（峰波纹值）。

其中，在通带内要求1-δp<=|Ha(jΩ)|<=1。

有时候指标也可以这样表示：用通带最大衰减 αp 和阻带最小衰减 αs表示，它们的单位是dB：它的定义为αp=-20lg(1-δp);αs=-20lg(δs)。

    第二种：用参数ε，A来表示,它们和δp,δs的关系是：ε=[(1-δp)^(-2)-1]^(0.5),  A=1/δs。

以上是先对滤波器的一个了解，我们需要理解各个指标的意思，这样对下面的描述才能更理解，以下是我个人的理解：

我们要设计数字滤波器，总的归纳起来有以下步骤组成：

1，先根据指标，用函数得到滤波器的阶数和截止频率

2，根据所得到的阶数，用函数得到滤波器的零点，极点，增益因子

3，根据零点，极点，增益因子，用函数求得滤波器的传输函数的分子，分母系数。这时候我们得到的是模拟的低通滤波器。

4，根据需要，用以下函数把模拟低通滤波器转换成低通，高通，带通或带阻滤波器。

  低通：[B,A]=lp2lp(b,a,Wo)

  高通：[B,A]=lp2hp(b,a,Wo)

  带通：[B,A]=lp2bp(b,a,Wo,Bw)

  带阻：[B,A]=lp2bs(b,a,Wo,Bw)

到这已经差不多大功告成了，但是我们这时候求得的是s域上的，有实，虚轴。可是我们要的是数字滤波器，也是对离散信号进行滤波，于是我们就要把s域转换成z域。这时候我们就要用双线性变换了。

5，用双线性变换得到数字滤波器。在这里我们要注意，既然是双线性变换，我们在设计之前，应该先对频率进行