matlab中的低通freqs,MATLAB中滤波器设计系列（二）低通模拟原型滤波器

在系列(一)中我们已经说过，几乎所有滤波器的设计都是以低通模拟原型滤波器为基础的，那么在本文中，主要针对低通模拟原型滤波器来进行讨论分析。

低通模拟原型滤波器主要有以下五种：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。下一讲我们讨论如何由模拟低通原型滤波器得出低通、高通、带通、带阻这四种类型的模拟滤波器。