[信息与通信]数字信号处理.pptVIP

作业 1. 概述模拟滤波器的产生与发展； 2. 概述数字滤波器的产生与发展 作业 1.编程设计一个N=4阶，Ωc=0.4π的巴特沃思滤波器，画出幅频特性曲线 2.编程设计一个N=6阶， Ωc=0.4π的且比雪夫滤波器，画出幅频特性曲线 4．系统函数 5 4．系统函数 5 例5.2 求三阶巴特沃斯低通滤波器的系统函数，Ωc=2rad/s， 画出频率响应图。 解 方法1：按幅度平方函数求解。因为 则 所有的极点为 选出左半平面的三个极点： = 例 5 解 方法2：查表得知归一化?三阶巴特沃斯低通滤波器的 = = 用s/2代替式中的s得到 查表法 5 例5.3 设计一个模拟低通巴特沃斯滤波器，且满足以下要求：通带截止频率 ，通带波纹 阻带截止频率 阻带波纹 例 5 解：明确要求谁？ 谁未知？N和Ωc，设计的实质是这两个参数的设计 续例 5 解方程组： N取大于等于所求结果的最小整数，所以取N=7 续例 5 现在，在上面两个数之间可任选Ωc值。现选Ωc=0.7，这样就必须设计一个N=7和Ωc=0.7 的巴特沃斯滤波器。过程同上例题。 例5.3巴特沃思低通滤波器频率响应 续例 5 MATLAB程序为 rp=1; rs=17; wp=0.2*pi; ws=0.3*pi; %求取滤波器阶数和3dB截止频率点 [N1,wn1]=buttord(wp,ws,rp ,rs,s); %巴特沃斯模拟低通原型滤波器设计 [b1,a1] = butter(N1,wn1,s) ; 续例 5 %巴特沃斯模拟低通滤波器频率响应 [h1,w1]=freqs(b1,a1); plot(w1/(2*pi),20*log(abs(h1))) grid on 运行结果为 滤波器阶数N1=7， 3dB截止频率wn1=0.6283 分子和分母多项式系数 b1 =0 0 0 0 0 0 0 0.0942 a1 =1.0000 3.2070 5.1424 5.3029 3.7842 1.8688 0.5935 0.0942 续例 5 切比雪夫滤波器的幅度特性是在通带或阻带具有等波纹特性。幅度特性在通带中是等波纹的，在阻带中是单调的，称为切比雪夫Ⅰ型。幅度特性在通带内是单调下降的，在阻带内是等波纹的，称为切比雪夫Ⅱ型。 5.3 切比雪夫滤波器 5 切比雪夫Ⅰ型和巴特沃斯低通滤波器 5.3 切比雪夫滤波器 5 切比雪夫Ⅱ型 式中，Ωc为通带截止频率，不一定是3dB衰减点；ε表示通带波纹的大小，0ε1，其值越大，通带波动越大；N为滤波器的阶数；VN(x)为N阶切比雪夫多项式。 幅度平方函数 5 N阶切比雪夫多项式 5 设x=Ω/Ωc为归一化频率，则切比雪夫滤波器VN(x)为N阶切多项式为 V V V V V V V V V V V 1 x 2x2-1 4x3-3x 8x4-8x2+1 16x5-20x3+5x 32x6-48x4+18x2-1 0 1 2 3 4 5 6 VN(x) N 由切比雪夫幅度平方函数 分析其函数特点为 2. 幅度函数特点 5 N为2、5、7切比雪夫幅度平方函数特性 当Ω=0时， 当Ω=Ωc时， 当ΩΩc时，即在通带内，|Ha(jΩ)|在(1, )之间等波纹起伏； 当ΩΩc时，即在通带外，随着Ω的增大，|Ha(jΩ)|迅速单调下降趋向0。 2. 幅度函数特点 5 2. 幅度函数特点 5 3. N阶特性 5 阶数N等于通带内最大和最小值个数的总和。 可由幅频特性中看出N阶数。且当： N=奇数，则Ω=0处有一最大值， N=偶数，则Ω=0处有一最小值。 N值如果未知，则是根据阻带的边界条件来确定。 切比雪夫幅度平方公式 3. N阶特性 5 思考：为什么不用通带截止频率？ N不是整数怎么处理？ 3. N阶特性 5 切比雪夫滤波器特性参数有三个ε、Ωc和N，通常通带截止频率Ωc预先给定。通带波纹Ap(以dB表示)的定义为 可以推导出Ap与ε关系 4. 性能指标 5 由切比雪夫幅度平方函数(归一化)可得 用s代替jΩ，可以得到 5. 系统函数 5 可以证明， 的极点分布在一个椭圆上，椭圆的短轴半径为a，长轴半径为b，取左半平面的N个极点构成，设这N个极点为 5. 系统函数 5 根据极点和其幅度平方函数可以得到切比雪夫滤波器的系统函数为 5. 系统函数 5 总结以上内容切比雪夫