重庆邮电大学本科毕业设计（论文）
图1- 1滤波器的幅频特性
第三节 滤波器的技术指标
一、 滤波器的技术指标及表示方法
设计中确定了实际需求就需要确定技术指标。具体的，滤波器的技术指标可用以下几种方式来表示：
①幅度范围定义了通带的最小幅度错误！未找到引用源。和阻带的最大幅度错误！未找到引用源。。
②幅度容限值指定通带幅度减小最大值?p?1?MP和阻带幅度的最大值?s?1?Ms ③幅度纹波容限描述了通带内幅度变化最大量错误！未找到引用源。和阻带内幅度变化最大量错误！未找到引用源。。
④衰减范围用dB表示，指定通带内的最大衰减错误！未找到引用源。和阻带内的最小衰减错误！未找到引用源。。
⑤增益范围用dB表示，通带内最小增益Gp??Ap恶化阻带内最大增益Gs??As。
二、 四种选频滤波器的技术指标
①低通滤波器的技术指标： H(ej?)=1,0
②带通滤波器的技术指标： H(ej?)=0，0<错误！未找到引用源。<0.4错误！未找到引用源。；H(ej?)=1；0.4?
H(ej?)=0，0.6<错误！未找到引用源。
③高通滤波器的技术指标： H(ej?)=0，0
④带阻滤波器的技术指标： H(ej?) =1，0<错误！未找到引用源。<0.3?；H(ej?)=0；0.3错误！未找到引用源。<错误！未找到引用源。<0.7错误！未找到引用源。；
。 H(ej?)=1，0.7?
第四节 滤波器的系统特性
- 6 -
重庆邮电大学本科毕业设计（论文）
已知在滤波器的设计的过程中，我们需要考虑信号与系统中所要求的很多因素，例如稳定性等。现在来具体介绍数字滤波器的系统特性。
一、 稳定性
稳定系统是系统只要没有外部源的激励，就能保持静止状态没有输出。或者输入有界输入即输入值在任意时刻都小于某一有限值，输出也能呈现有界性。这样的系统被称作稳定系统。考虑到系统的稳定性提供了性能实用性等方面的影响，故它也是系统设计中的一个主要的问题，因此有必要通过确定特定系统参数值的范围来确定系统稳定性。
二、 时不变性
如果输入信号的时移会引起相同的输出信号的时移就叫做时不变系统。用数学表达的方式，即如果输入x(t)产生输出y(t)，则输入x(t?t0)产生y(t?t0)。也就是说，在这种系统中，任意时刻输入的信号，输出的都是同一个波形。也可以说，输入延迟
t0秒，输出也保持同样的延迟步调，这个系统就是时不变的。同理，对于离散时间系统而言，要想时不变就需要满足输入信号x(k?k0)产生对应输出信号y(k?k0)。时不变系统的特性和参数不随时间变化，这一点在差分方程中表现的尤为明显。相反，非时不变的系统称为时变系统。
三、 线性
首先介绍叠加原理：几个输入x1(t),x2(t),...,xN(t)同时作用于系统，其响应y(t)等于每个输入单独作用于系统的响应之和。具体的，如果系统满足：
①叠加性：输入x1(t)产生输出y1(t)，输入x2(t)产生输出y2(t)，那么对所有的输入x1(t)和x2(t)有，输入(x1(t) + x2(t))产生输出(y1(t)+y2(t))。
②齐次性：输入x(t)产生输出y(t)，同时对任何输入及任何常数a，都有输入ax(t)得到的输出为ay(t)。
同时满足叠加性和齐次性的系统被称作线性系统。线性时不变系统可以用它的单位脉冲响应来表征。其输入输出的描述法通常采用线性常系数差分方程。
实际研究的大部分系统都是线性时不变系统，可以缩写为LTI系统。这类系统对输入信号的响应可以通过求取输入系统的冲击响应的卷积得到。输入输出的关系变得相对简单。它的缺点在于：用卷积运算来求解LTI系统的响应时并不容易，同时系统的真实情况也不能很好地反映。但无论怎样，LTI系统在信号与系统的研究中仍然占有举足轻重的位置。
- 7 -
重庆邮电大学本科毕业设计（论文）
四、 因果性和最小相位系统
如果系统以时间为自变量，且输出仅与现在和过去的输入值有关，就叫做因果系统。即y(t)仅取决于x(t)在??t的值。相反，如果和未来有关或是在实际中是不可能实现的，这类系统被称作是非因果系统。
为了实现滤波器设计实际总就必须保证系统的可实现性，所以必须为因果系统。如果一个因果稳定的线性时不变系统，它的系统函数H(z)所有的零点都在单位圆内则称其为最小相位系统。相反，如果零点都在单位圆之外是最大相位系统。最小相位系统的相频特性，其相角变化范围是最小的。
第二章 模拟滤波器设计分析
第一节 滤波器逼近的一般方法
模拟滤波器基于连续时间系统，在频谱变换中采用s变换。已知滤波器设计的首要任务就是要确定满足技术指标要求的系统函数，系统函数必须是物理可实现的，也就是要满足因果性和稳定性。为此系统就必须满足一下几个条件：
①H?s?是s的实系数有理函数，也就是滤波器的单位脉冲响应h(t)是实函数 ②H?s?的极点必须全部位于s平面的左半平面。（稳定性要求） ③H?s?的分子多项式的阶数必须小于分母多项式的系数。 由上推论可得，滤波器的频率响应满足
|H?j??|2?H?j??H??j???H?s?H??s?|s?j?
（2- 1）
由于系统函数H(s)是s的有理函数，则H(s)H(-s)是s2的有理函数，所以|H(j?)|2必为?2的有理函数。并且在滤波器逼近函数问题中，一般由|H(j?)|2确定H(s)。
在模拟滤波器中通常用一些逼近函数进行设计，常见的逼近函数有：巴特沃思逼近，切比雪夫（I和II型）逼近，椭圆逼近，贝塞尔逼近。图2-1给出了几种逼近函数设计的低通滤波器幅频特性，他们具有相同的阶数和波纹要求。
- 8 -
重庆邮电大学本科毕业设计（论文）
|H(j?)|11|H(j?)|1|H(j?)|1|H(j?)|01234?01234?01234?01234? (a) 巴特沃思滤波器 (b)切比雪夫I型滤波器 (c) 切比雪夫II滤波器 (d)椭圆滤波器
图2- 1几种逼近函数的幅频特性
从图中观察可以得知，在相同的阶数和波纹要求之下，椭圆滤波器的选择性相对最优，它的过渡带也相对更窄，带外抑制更加陡峭。其次是切比雪夫滤波器，最后才是巴特沃思滤波器。切比雪夫的两种类型所变现出的波纹和单调型出现频带相反也从图中明显观察得知。
第二节 巴特沃思滤波器实现与分析
一、 巴特沃思逼近
巴特沃思（Butterworth）低通滤波器的平方幅度响应函数为：
|H(j?)|2?H????21 2N1?(?/?c)（2- 2）
式中： N为滤波器的阶数，通带频率取为1rad/s。
由2-2式可计算出：??0rad/s时，|H(j0)|?1；??1rad/s时，H(j1)?1/2；也就是说在通带频率处衰减约3dB。故有时也称巴特沃思滤波器3dB通带频率为1rad/s；???时，H(j?)?0。
阶数N值和其他参数的关系可被归纳为：
所示。
log(100.1As?1) N???2log?s（2- 3）
巴特沃思系统函数的特点可以被总结为以下几点：它的分母多项式A(s)如表2-1①阶数N越高，特性越接近于矩形，过渡带越窄。
②通带没有波纹，其频率响应随着频率增大而平滑单调下降。 ③H(j?)在通带内具有最大平坦特性。传递函数无零点。 ④缺点是从通带到阻带的过渡较大。
- 9 -
重庆邮电大学本科毕业设计（论文）
表2- 1巴特沃思低通滤波器分母多项式的系数A(s)?sN?a1sN?1???aN?1s?1N 2 3 4 5 6 a1 1.414214 2.000000 2.613126 3.236068 3.863703 a2 2.000000 3.414214 5.236068 7.464102 a3 2.613126 5.236068 9.141620 a4 3.236068 7.464102 a5 3.863703
二、 巴特沃思滤波器的实现与分析
在MATLAB信号处理工具箱中巴特沃思滤波器函数的调用格式如下：[z,p,k,]=buttap（n）式中：n为巴特沃思滤波器的阶数，z,p,k分别为滤波器的零点、极点和增益。设计中还涉及到[num,den]=zp2tf(z,p,k)函数，它是传递函数（transfer function），num为分子多项式的系数，而den是分母多项式的系数。
设计并产生一个25阶低通滤波器原型，其实现的MATLAB程序代码如下：
[Z,P,K]=buttap(25);[num,den]=zp2tf(Z,P,K);freqs(num,den);运行程序效果如图2-2示：
图2- 2巴特沃思模拟滤波器特性图（n=25） 图2- 3巴特沃思低通滤波器幅频响应曲线
从图2-2中可以明显观察得知，巴特沃思模拟滤波器不具备线性相位。再者设计并绘制阶数分别为3、5、10、15的巴特沃思滤波器的幅频响应曲线。其实现的MATLAB程序代码见附录代码2-1，运行程序幅频响应绘制曲线如图2-3所示。从图中可论证，滤波特性一般，且阶数N越大，滤波效果越接近理想值矩形过渡带越窄。这种滤波器无论是在通带还是阻带内都是频率单调的函数。可以发现，当通带边界处满足指标要求，通带内有余量，可以将精度均匀分布在整个通带内。为此，选择传递函数具有等波纹特性的滤波器，就能在阶数较低的情况下满足系统要求。
- 10 -