由模拟滤波器设计IIR数字滤波器.docVIP

由模拟滤波器设计IIR数字滤波器 为了从模拟滤波器设计IIR数字滤波器，必须先设计一个满足技术指标的模拟滤波器，然后将其数字化，即从s平面映射到z平面，得到所需的数字滤波器。虽然IIR数字滤波器的设计本质上并不取决于连续时间滤波器的设计，但是因为在许多应用中，数字滤波器就是用来模仿模拟滤波器功能的，所以由模拟滤波器转化为数字滤波器是很自然的是。另外，模拟滤波器的设计技巧非常成熟，不仅有封闭形式的公式，而且设计系数已经表格化。因此，有模拟滤波器设计数字滤波器的方法准确、简便，是目前最普遍采用的方法。 在模拟滤波器的设计中，低通滤波器是最基本的。设计模拟滤波器的方法有多种，如巴特沃兹（Butterworth）型、切比雪夫型（Chebyshev）型、椭圆型（Elliptic）型滤波器。为了能从模拟滤波器的低通原型设计各种IIR DF,一般需如下四个步骤： 把要求的低通（LP）、高通（HP）、带通（BP）、或带阻（BS）的特征频率参数转化为模拟低通滤波器低通原型的设计参数。 用模拟逼近的方法获的巴特沃兹、切比雪夫或椭圆模拟低通原型的传递函数Hp(s)。 通过s平面到z平面的映射关系，由Hp(s)求出相应的数字低通的系统函数Hp(z)。 用数字域的频率变换，从Hp(z)求出所需的数字LP、HP、BP、或BS数字滤波器的系统函数H(z)。下面将对上述四个步骤分别加以介绍。 　 5.2.1 在模拟滤波器的设计中，巴特沃兹、切比雪夫以及椭圆滤波器的设计都是低通逼近。所以，如果设计的滤波器不是低通，就需要将HP、BP、或BS的频率参数变换为低通原型的相应参数。这个变换是在模拟域进行的，所以叫模拟频域变换。 低通原型的设计参数 设计一个低通滤波器需要给出4个参数：通带临界频率fp(Hz)，阻带临界频率fs(Hz)，通带最大衰耗αp(dB)，阻带最小衰耗αs(dB)。这4个参数构成的低通样板图如图5.2所示。图中横坐标代表频率f(Hz)，纵坐标是相对域零频的幅值衰耗 dB 从零频到fp称为通带，fs到称为∞阻带，fp到fs称为过渡带。参数αp是通带中允许的衰耗值波动边界。当|H(f)| 下降到|H(0)|的0.707（半功率点）时，A(f)=3dB，所以，αp最多取3dB。参数αs是阻带最小衰耗，一般在30dB～80dB之间。所设计的低通滤波器，只要器相对域零频的幅值衰耗曲线落在图中的空白出，该滤波器的性能就算合格了。 需要指出的是，任何实际的滤波器与图5.1的理想滤波器器都是由距离的。理想滤波器没有过渡带，相应的冲激响应h(n)为非因果序列，而非因果系统在工程上是不可实现的。为使设计的过渡带较窄，需要提高滤波器的阶数，这就相应提高了滤波器的成本和价格。因此，在滤波器的设计中，需要综合考虑滤波器的性能及价格等因素。 高通→低通的频率变换 设计一个高通滤波器，也应给出4个参数：通带临界频率fp，阻带临界频率fs，通带最大衰耗αp，阻带最小衰耗αs（见图5.3）。其中最主要的参数是fp，它代表滤波器的截止频率。 为了将给定的高通滤波器参数变换为低通原型的参数，需要将高通的阻带和通带的位置互换，并保持截止频率fp不变。这种变换是变量f的倒量变换，只有将f换成1/f，则零频变为∞，而∞变为零频，原来的高通也就成了低通。为了保持fp位置不变，可先将f除以fp，把频率点fp归一，然后球倒数，最后再乘以fp反归一，把频率点fp回归原位。所以整个过程是 (归一) 则变换后低通原型的通带频率为fp（不变），阻带临界频率为fsˊ=fp2/fs，αp和αs不变（见图5.4）。 带通→低通频率变换 设计一个带通滤波器，需要给出6个参数：通带上下临界频率fp1、fp2，上、下阻带临界频率fs1、fs2，通带最大衰耗αp和阻带最小衰耗αs(见图5.5)。 通带上、下临界几何中心称为中心频率f0。 f0= (5.13) 参照图5.6，可见带通到低通的频率变换目标是： 将通带中心频率f0变换到低通零频点：f0→0。 将带通的半频率轴f：（0,+∞）扩展为低通频率轴fˊ(-∞,∞),宽度扩大一倍。这样低通fˊ的宽度应是（fp2-fp1）的两倍，就是 即 (5.14) (3)f轴上的几何对称点变换为fˊ轴上对零点的镜像对称点。频率f的几何对称点是f02/f(f02/f与f的乘积为f02)，而fˊ镜像对称点是- fˊ。就是说，如果f→fˊ，则必有f02/f→- fˊ。 （4）相对位置不变，也就是指f轴上两几何对称点距离与总带宽之比应等于-fˊ轴上两像对称点的距离与总带宽之比，即 (5.15) 将式（5.14）代入式（5.15），可得 (5.16) 将f=fs2代入上式，得 设fp1fp2≈fs1fs2,则 因此可以得出结论，为了将给定得带通滤波器的参数变换为低通原型