《第十章模拟与数字滤波器讲义资料》由会员分享，可在线阅读，更多相关《第十章模拟与数字滤波器讲义资料（72页珍藏版）》请在金锄头文库上搜索。
1、第十章 模拟与数字滤波器一、滤波器的功能 滤波器的功能是对输入信号进行滤波以增强所需信号部分，抑制不要的部分。二、模拟滤波器（AF）是指输入输出均为模拟信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。10.1 引言三、数字滤波器（DF）是指输入输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。( )x t( )y t数字滤波器分类低通高通带通带阻n五、用数字滤波方法处理模拟信号n六、主要内容nAF的设计nDF的设计限带滤波A/D数字滤波D/A平滑滤波Xa(t)PrFADCDSPDACPoFYa(t)10.4 模拟
2、滤波器的设计n模拟滤波器n巴特沃斯 Butterworth 滤波器n切比雪夫 Chebyshev 滤波器n频率响应特性n及性能指标n系统函数一、实际模拟低通滤波器的特性：通带内允许误差：阻带衰减：通带截止频率：阻带截止频率二、Ha(s) Ha(-s)的零极点分布 （1） 的极点在s平面的左半平面；确定（2） 的零点在s平面的左半平面；（3）虚轴上的零点取一半；（4）系数为原系数开平方根。例：解：极点：零点： （二阶）零点：的极点：三、Butterworth滤波器1、逼近函数：称Wc为Butterworth低通滤波器的3分贝带宽N为滤波器的阶数Wc为通带截止频率2.幅频特性n 通带内幅度特性平坦
3、，单调减小； 过渡带及阻带内快速单调减小； 当WWs（阻带截止频率）时，衰减的d2为阻带最小衰减N越大，越接近理想特性Butterworth滤波器是一个全极点滤波器，其极点： 3.幅度平方函数的极点分布： 极点在s平面呈象限对称，分布在Buttterworth圆上， 共2N点;极点间的角度间隔为 极点不落在虚轴上，前N个极点落在s平面的半平面 N为奇数，实轴上有极点，N为偶数，实轴上无极点Ha(s) Ha(-s)的零极点分布(a) N=3 (三阶） (b)N=4 （四阶）4、滤波器系统函数Ha(s) 的确定：Ha(s) 的确定的方法：方法一：直接求取方法二：根据N查P186表101，得 ，5、
4、Butterworth滤波器设计步骤：n(2)根据 求出滤波器阶数N：n(1)确定技术指标：n(3)根据阻带衰减 确定 ：(4)根据N查P186表101，得 ，N：滤波器的阶数Wc ：截止频率，不一定为3dB带宽0e1，表示通带波纹大小，e越大，波纹越大TN(x) ：N阶Chebyshev多项式Type I Chebyshev四、Chebyshev滤波器1、逼近函数：切比雪夫多项式特点： (1)当|x|1时，|TN(x)|1,在|x|1时，TN(x) 随x单调上升; (3) (4)2、幅频特性：，N为偶数，N为奇数n n 通带内：在1和 间等波纹起伏n 通带外：迅速单调下降趋向0（1）（2）（
5、3）（4） 3.幅度平方函数的极点分布：令：记：则： 极点不落在虚轴上，前N个极点落在s平面的左半平面 极点在s平面呈象限对称，分布在椭圆上，共2N点其中：所以：4、滤波器系统函数Ha(s) 的确定：Ha(s) 的确定的方法：方法一：直接求取方法二：根据N、 查P196表104，得 ，5、Chebyshev滤波器设计步骤：n(3)根据 求出滤波器阶数N：n(1)确定技术指标：n(2) 根据 确定 ：(4)根据N、 查P196表104，得 ，n设计思想： s 平面 z 平面模拟系统 数字系统 设计步骤： 1、数字滤波器的性能指标 模拟滤波器的性能指标 2、设计模拟原型滤波器 3、 4、通过算法实
6、现H(z)10.6 无限冲激响应（IIR）数字滤波器两种设计方法：- 冲激响应不变法- 双线性变换法1、原理数字滤波器的单位冲激响应h(n) 模仿模拟滤波器的单位冲激响应ha(t)T抽样周期一、冲激响应不变法n系数相同：n极点：s 平面 z 平面k2、稳定性与逼近程度1）、稳定性不变AF稳定:DF稳定：2）、逼近程度当T 很小时，数字滤波器增益很大，易溢出，需修正3、修正4、设计步骤DF性能指标AF性能指标（1）指标转换（3）（2）设计AF试用冲激响应不变法，设计IIR数字滤波器例：设模拟滤波器的系统函数为解：设T = 1s，则n优点：n缺点：n保持线性关系：w=WT线性相位模拟滤波器转变为线
7、性相位数字滤波器n频率响应混迭只适用于限带的低通、带通滤波器nh(n)完全模仿模拟滤波器的单位抽样响应ha(t) 时域逼近良好5、优缺点二、 双线性变换法q冲激响应不变法：时域模仿逼近缺点是产生频率响应的混叠失真q为了克服这一缺点，采用双线性变换法。q使数字滤波器的频率响应与模拟滤波器的频率响应相似n设计思想1、脉冲响应不变法的主要缺点：对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大，而且不能用来设计高通和带阻滤波器。原因：从S平面到Z平面的映射是多值映射关系脉冲响应不变法映射过程图示 2、双线性变换法的改进：为避免频率的“混叠效应
8、” ，分两步完成S平面到Z平面的映射。 将S平面压缩到某一中介的S1平面的一条横带域 通过标准的变换将此横带域映射到整个Z平面上去。双线性变换的映射过程图示 3、双线性变换法的基本思路：从频率响应出发，直接使数字滤波器的频率响应 ，逼近模拟滤波器的频率响应 ,进而求得H(z)。n“双线性变换法”设计方法 通过Z变换： 将1映射到Z平面的单位圆上。 通过正切变换：将S平面的j轴压缩到S1平面的j1轴上的 内。 将正切变换延拓到整个将正切变换延拓到整个S S平面，得到平面，得到S S平面到平面到S S1 1平面的平面的映射关系映射关系： 将S1平面按 映射到Z平面得到： 令 或以保低频处有较确切的
9、对应关系：4、稳定与逼近情况（1）逼近情况s平面虚轴z平面单位圆（2）稳定情况左半平面单位圆内s平面z平面右半平面单位圆外虚轴单位圆上代入5、优缺点n优点：避免了频率响应的混迭现象2）双线性变换后，频率响应产生了畸变，但通带阻带仍为常数特性，不受影响，只是临界频率点产生畸变。1）线性相位模拟滤波器 非线性相位数字滤波器q 缺点：除了零频率附近，W与w之间严重非线性预畸变给定数字滤波器的截止频率 wc，则按Wc设计模拟滤波器，正好映射到wc为截止频率的数字滤波器6、设计步骤DF性能指标AF性能指标（1）指标转换（3）（2）设计AFPn16110.7 有限冲激响应（FIR）数字滤波器n线性相位FI
10、R滤波器的特点n窗函数设计法IIR数字滤波器：可以利用模拟滤波器设计但相位非线性FIR数字滤波器： 可以严格线性相位因果稳定系统可用FFT计算但阶次比IIR滤波器要高得多 FIR滤波器的单位冲激响应：系统函数：在 z 平面有N 1 个零点在 z = 0 处是N 1 阶极点 一线性相位特性1.无失真传输条件单位冲激响应 为实数，且关于奇对称或偶对称,即2.获取线性相位的条件（1）h(n)偶对称，N为奇数其中:g(2) h(n)偶对称，N为偶数其中， gg(3)h(n)奇对称，N为奇数其中,gg其中：(4) h(n)奇对称，N为偶数gg一、设计思想 先给定理想filter的频响 ，所要求设计一个F
11、IR的filter的频响为 ，使 逼近二、设计原理 设计是在时域进行的，先用傅氏反变换求出理想filter的单位抽样响应 ，然后加时间窗对 截断，以求得FIR filter的单位冲激响应h(n)。7-3 窗函数设计法或三、窗函数 1、基本概念（1）窗谱：窗函数的频响的幅度函数亦称作窗谱。（2）对窗函数要求 a）希望窗谱主瓣尽量窄，以获得较陡的过渡带，这 是因为过渡带等于主瓣宽度。 b）尽量减少窗谱最大旁瓣的相对幅度，这样可使肩峰 和波纹减少。 2、矩形窗 时域表达式： 频域表达式（频谱）：幅度函数：主瓣宽度3、三角形（Bartlett）窗时域表达式：10 1 2 3 4 频谱： 第一对零点为 ，即 ，所以主瓣宽度 ，比矩形宽一倍。4、汉宁窗（升余弦窗）旁瓣较大程度地互相抵消，但主瓣宽度为5、海明窗，（改进升余弦窗） 其窗函数为频率响应的幅度函数为 其主瓣宽度仍为 ，（旁瓣峰值/主瓣峰值）1%有99.963%的能量集中在主瓣内。四、设计步骤（1）给定频响函数（2）求出单位冲激响应（3）根据过渡带宽度和阻带最小衰减，查P342表7-3， 确定窗函数及N的大小：由阻带衰减率 确定 窗函数，由过渡带宽确定N。 （4）求 （5）求 ，检验是否达到所要求指标 例1：设计具有线性相位的FIR 低通滤波器，具体要求：阻带衰减 -20dB,过渡带宽度为
