揭开IIR滤波器中映射的面纱-以双线性变换法和数字频率转换为例

背景

         最近学习双线性变换法的时候，一直不太理解什么是预畸变？畸变发生在哪？s到z平面的映射中有畸变吗？既然双线性变化中 模拟频率和数字频率之间的关系为 Ω = 2/Ttan(w/2)，为了得到所需要的数字频率 ω，自然要按照 Ω = 2/Ttan(w/2）所对应的模拟频率来设计模拟滤波器，但为什么要把 Ω = 2/Ttan(w/2）称为预畸 变呢？在接下来的数字域频率转换中也有类似的疑惑，为啥到LPF的频率映射中可以任意指定一个起始频率呢？现有的教材中一般都给出了预畸变的概念，但概念解释不太清楚。我查阅了一些文献资料，最终将自己的理解整理如下：

双线性变换法的处理思路

        首先我们明确我们的目的：需要设计一个DF（以低通为例），我们首先：

        一，通过角频率的映射，将他映射成一个AF；

        二：设计一个AF；

        三：通过双线性变换法映射回DF，分别对应我下图中的步骤1，2，3

步骤

        一: 给定一个DF的设计指标ws,wp,As,Ap,我们首先设计一个AF，AF的性能指标有两种映射关系：

        ① Ω = w/T ，ws,wp映射到了B1，B2点；

        ② Ω = 2/T*tan(w/2)  ,ws,wp映射到了A1,A2点；

      

二: 分别设计两种情况下的AF。

三: 设计好了AF,我们拿到了H(s)，下面需要通过双线性变换法映射回z平面，拿到H(z),对于s与z平面而言，映射关系为：s = 2/T*(1-z^-1)/(1+z^-1),对应频率响应而言，即为Ω = 2/Ttan(w/2) 

至此，我们可以发现，畸变的原因是因为 s到z是非线性映射，无论角频率那里的映射是如何的，最后一步的映射都要经过这个非线性映射，因此，只有对频率做对应的非线性变换，才会使得H(s)到H(z)的映射是不失真的，否则必有失真。

如果我们拿着模拟指标直接设计一个AF，那么他经过H(s)到H(z)的映射后，一定会失真。

再来聊聊数字域频率转换中的映射

数字域频率转换的思路

见下图。

我们可以看到，如果首先忽略HLPF(z)的具体设计步骤，认为直接可以设计出来一个满足要求的HLPF(z)，则一个数字高通滤波器的设计步骤为：

        一，通过角频率的映射，将他映射成一个HLPF；

        二：设计一个HLPF；

        三：映射回H HPF（z）

我们可以发现，这和双线性变换法的步骤是一样的！！！

为什么可以任意指定一个或?

        到此，我们应该知道，步骤①中低通滤波器和高通滤波器之间的角频率映射关系，其实就是步骤⑤中，z和Z取单位圆的时候的映射关系而已，即:

对于高通滤波器，其中a满足：

我们可以发现，z和Z平面的映射关系中的a是和、有关的，而和之间的映射关系，又取决于a,因此，你可以任意指定一个，当你指定以后，a随之确定，z和Z平面的映射关系就随之确定，映射变为了唯一的映射；当你改变指定的，只会使得LPF左右移动一下而已，映射关系随之改变，根本不会影响到你设计的HPF。而对于双线性变换法，其中的映射关系是确知的，因此肯定是不能随便指定着玩的。