[数字信号处理]IIR滤波器的直接设计(C代码)

1.IIR滤波器的直接设计原理

        不利用模拟滤波器，直接进行数字滤波器的设计的方法，称为直接设计。回忆之前所说的IIR滤波器的直接设计，我们首先设计了巴特沃斯模拟滤波器，然后进行双线性变换，得到数字滤波器。我所使用的是巴特沃斯低通滤波器作为原型滤波器，其振幅特性如下所示。

        

        首先，我们先把数字滤波器的指标，根据下式转为模拟滤波器的指标。

        然后根据模拟滤波器的设计指标，计算次数N，然后计算极点，最后选择出稳定的极点，计算模拟滤波器的传递函数。然后，我们通过拉普拉斯反变换，再使用差分替代微分，得到了S平面-Z平面的对应关系，将其转为数字滤波器。

        详细的，戳链接  ---------[数字信号处理]IIR滤波器的间接设计（C代码）

        现在，我们将这个过程，稍微调整一下。

      1.1 巴特沃斯低通数字滤波器的设计

         首先，转换指标的公式，直接带入巴特沃斯低通模拟滤波器的振幅特性表达式，然后化简，我们可以得到一个新的振幅特性

同样的，这个振幅特性具有巴特沃斯低通模拟滤波器的一些特性。函数为单调递减的，通带与阻带都没有纹波。这个新的滤波器，称为巴特沃斯低通数字滤波器。

        与之前一样的，我们的首先先确定所需要的滤波器的次数N，根据次数计算出极点，选择稳定的极点来计算出传递函数，然后就可以得到滤波器的系数了。

        首先是次数的计算，次数的计算，还是要根据阻带衰减指标，根据下式计算。

将巴特沃斯低通数字滤波器的振幅特性带入，我们就可以计算出所需要的滤波器的次数N。

         下一步，我们是需要计算极点。这里与之前的间接设计不同，我们需要做一些变形。首先，将其振幅特性做平方，变为

然后，把 部分稍微做下变形，

然后，将 替换为Z，带入，则得到了如下式子

通过这个式子，就可以很方便的计算极点与零点。很容易的能看出，这个滤波器的零点是-1，并且为N重极点（这里是振幅特性的平方所以不是2N）。此时，分母多项式为

由于一步解开很麻烦，我们先将这个式子，关于 解开，我们可以得到如下解，

再解z，然后再我们就可以得到极点了，如下所示。

这里所求得的极点为2N个，为了所设计的滤波器是稳定的，我们需要选择出稳定的极点。在Z平面内，其摸小于1的极点，就是稳定的极点。或者来说，只要滤波器的所有极点均在Z平面的单位圆里&