数字信号处理第六章2012.pptVIP

若令3dB截止频率?c =1，得归一化传输函数。 巴特沃斯低通滤波器的设计 通带中允许的最大衰减?p和通带截止频率?p； 阻带中允许的最小衰减?s和阻带截止频率?s； （一）根据性能指标 设计： 确定阶数N: 滤波器的阶数 令 同理 可推出 巴特沃斯低通滤波器的设计 2. 3dB通带截止频率?c: 或 3. 确定归一化传递函数Ha(p)。 4. 去归一化。 阻带指标有富裕量 通带指标有富裕量 巴特沃斯低通滤波器的设计 巴特沃斯低通滤波器的设计 （二）根据通带边界频率 ，通带最小幅度 阻带边界频率 ，阻带最大波纹 设计 由幅度平方函数 可得： 由 可得： 由 可得： 巴特沃斯低通滤波器的设计 例6.2.1: 已知通带截止频率fp=5kHz,通带最大衰减?p=2dB,阻带截止频率fs=12kHz,阻带最小衰减?s=30dB,设计一巴特沃兹模拟低通滤波器。 解（1）确定阶数N 取N=5 （2）查表确定归一化原型函数 （3）确定3dB截止频率?c （4）去归一化 用matlab工具箱函数设计巴特沃斯滤波器 ?1. [Z,P,K]=buttap（N） 用于计算N阶巴特沃斯模拟低通原型滤波器系统函数的零、极点和增益因子。Z和P是长度为N的列向量。 ?2. [N,wc]=buttord(wp,ws,Rp,As) 用于计算巴特沃斯模拟低通原型滤波器的阶数N和3dB截止频率wc。 其中Wp和Ws分别是通带和阻带的截止频率，其值 ，当其值为1时，表示0.5 fs Rp,RS分别是通带和阻带区的衰减。 4.[B,A]=butter(N,wc,’ftype’) 计算N阶巴特沃斯滤波器系统函数的分子、分母多项式的系数向量B和A。 3.[N,wc]=buttord(wp,ws,Rp,As，’s’) 其中Wp和Ws模拟角频率。 5.[B,A]=butter(N,wc,’ftype’，’s’) 计算N阶模拟巴特沃斯滤波器系统函数的分子、分母多项式的系数向量B和A。 切比雪夫（ Chebyshev ）低通滤波器 切比雪夫-I型低通滤波器的频率响应其通带内是等波纹的 切比雪夫-II型低通滤波器的频率响应其阻带内是等波纹的 (巴特沃兹低通滤波器的频率响应在通带和阻带内都是单调的) 3.对相同的频带要求，切比雪夫低通滤波器比巴特沃兹低通滤波器需要更低的阶数（牺牲单调性） 相同的阶数情况下，同样的边界频率,切比雪夫滤波器具有更好的幅度指标 N 为滤波器的阶数,ε为通带内波纹因子 Ωp截止频率 ?切比雪夫Ⅰ型低通滤波器幅度平方函数 N 阶切比雪夫多项式 切比雪夫多项式 当 归纳出高阶切比雪夫多项式的递推公式： 不同N 值的切比雪夫多项式的曲线 切比雪夫多项式的零点在 区间内； 当 时，曲线具有等纹波特性； 当 时，曲线按双曲余弦函数单调上升。 当|x|?1， 即在[0，?p]范围内， 当??p， 切比雪夫滤波器的幅频特性 当 Ω=0 N 奇数： ；N 偶数 其幅度平方函数与三个参数有关：?，?p，N 切比雪夫滤波器的幅度平方函数 1. 确定? 定义 因为 ?是通带内允许的波动大小。 2. 确定阶数N 令 由 或 3. 确定3dB截止频率?c 由 通常 上式仅取正号，得 以上参数确定后，可以求出N个极点的位置，从而综合出传递函数Ha(s)。 求Ha(s) Ha(-s)的极点可以通过求 的根得到。 选择上述多项式在左半平面的根作为Ha(s)的极点。 其中 归一化切比雪夫滤波器Ha(p) 去归一化切比雪夫滤波器Ha(s) 式中 切比雪夫滤波器的设计步骤： 1. 确定技术要求 2. 确定? 3. 确定N 4.归一化切比雪夫滤波器Ha(p) 式中 5.去归一化 切比雪夫滤波器的设计步骤： 用matlab函数设计切比雪夫滤波器 ?1. [Z,P,K]=cheb1ap（N,Rp） 用于计算N阶切比雪夫1型模拟低通原型滤波器系统函数的零、极点和增益因子。Z和P是长度为N的列向量。 ?2. [N,wpo]=cheb1ord(wp,ws,Rp,As) 用于计算切比雪夫1型模拟低通原型滤波器的阶数N和通带截止频率wpo。 其中Wp和Ws分别是数字域通带和阻带的截止频率，其值 ，当其值为1时，表示0.5 fs Rp,RS分别