在系列(一)中我们已经说过,几乎所有滤波器的设计都是以低通模拟原型滤波器为基础的,那么在本文中,主要针对低通模拟原型滤波器来进行讨论分析。
低通模拟原型滤波器主要有以下五种:butterworth滤波器、chebyshevI型滤波器、chebyshevII型滤波器、椭圆滤波器、bessel滤波器。下面分别对每种滤波器的设计进行讨论。进行讨论之前,我们假设各滤波器的阶数n已知。
1、butterworth滤波器
[z,p,k]=buttap(n);
[b,a]=zp2tf(z,p,k);
[H,W]=freqs(b,a,w);
magH2=(abs(H)).^2;
其中n为巴特沃斯低通模拟原型滤波器阶数,buttap为滤波器设计函数,生成滤波器对应的是滤波器的零极点模型,换句话说,就是用零极点模型来表示巴特沃斯模拟低通滤波器的。那么第一行表达式左边的三个参数z,p,k分别对应着滤波器的零点、极点和增益。第二个表达式是将上式所对应的零极点模型转换为传递函数形式,其中b为分子多项式系数,a为分母多项式系数。如果学过随机信号处理的话,可以将此处的传递函数与一般的AR模型、MA模型以及ARMA模型联系起来,这些内容我们以后再详细讨论。第三个表达式中,w为频率点,就相当于横轴坐标,只是该坐标为频率,因为要做出频率响应图形的。而函数freqs的作用就是求出该滤波器的复数频率响应。第四个表达式很显然就是求出该原型滤波器的幅频图了。
下面简单讨论其他四种滤波器即可,因为参数意义是基本一样的。
2、chebyshevI型滤波器
[z,p,k]=cheb1ap(n,Rp);
Rp为通带波纹,单位为dB,一般为1~5dB。其他参数类似于巴特沃斯原型滤波器。后面三个表达式与巴特沃斯一致,此处不再赘述。
3、chebyshevII型滤波器
[z,p,k]=cheb2ap(n,Rs);
Rs为阻带波纹,单位为dB,一般大于16dB。其他不再赘述。
4、椭圆滤波器
[z,p,k]=ellipap(n,Rp,Rs);
Rp一般1~5dB,Rs一般大于15dB。其他不再赘述。
5、bessel滤波器
[z,p,k]=besselap(n);
此处n一般小于25。其他不再赘述。
由以上讨论我们可以得出这五种低通模拟原型滤波器,需要说明的是,模拟低通滤波器的截止频率为1。下一讲我们讨论如何由模拟低通原型滤波器得出低通、高通、带通、带阻这四种类型的模拟滤波器。