数字信号处理教学课件第五章.pptVIP

* 5.巴特沃思滤波器的MATLAB实现 在MATLAB信号处理工具箱中，与Butterworth滤波器 相关的有三个函数： （1）[n,wn]=buttord(wp,ws,Rp,Rs,‘s’)，求出给定通带 截止频率wp、阻带截止频率ws、通带最大衰减Rp、阻带最 小衰减Rs条件下，所需要的最小Butterworth模拟滤波器阶 数n和3dB截止频率wn。 （2）模拟低通原型滤波器设计函数buttap，调用格式 为[z,p,k]=buttap(n)，其中参数z、p、k分别为滤波器的零 点、极点和增益。 * （3）Butterworth滤波器完全设计函数，调用格式为 [b,a]=butter(n,wn,‘ftype’,‘s’)，其中’ftype’为滤波器类型，缺省时表示低通或带通滤波器。b、a分别为滤波器传输函数 分子、分母多项式系数向量。 例5-4 利用MATLAB实现例5-2所述巴特沃斯滤波器。 解 MATLAB实现程序如下： wp=2*pi*5000; ws=2*pi*10000; Rp=3; Rs=20; %计算满足性能指标的滤波器阶数n和3dB截止频率wn [n,wn]=buttord(wp,ws,Rp,Rs,s) * %设计模拟低通原型滤波器，传输函数为零点、极点和增益 形式 [z,p,k]=buttap(n) %采用另一种方法设计模拟低通滤波器，传输函数为分子、 分母多项式形式 [b,a]=butter(n,wn,s) %求模拟滤波器的频率响应 [db,mag,pha,w]=freqs_m(b,a,2*ws); %绘图 plot(w,db,k);axis([0,2*ws/pi,-50,1]);grid on; * 仿真曲线 * 5.3 脉冲响应不变法 1.变换原理 数字滤波器的单位脉冲响应h(n)模仿模拟滤波 器的冲激响应ha(t)，使h(n)等于ha(t)的采样值，即 h(n)＝ha(nT)。 设模拟滤波器的传输函数Ha(s)只有单阶极点， 且分母的阶数高于分子阶数，则 其拉普拉斯反变换 * 对ha(t)采样，得到数字滤波器的单位脉冲响应 h(n)为 对h(n)取z变换，即得到数字滤波器的系统函数 可见，s平面的极点si映射到z平面的极点 ， 而对应系数不变。 * 例5-5 利用脉冲响应不变法将模拟滤波器 变换为数字滤波器 ，采样周期 。 解：模拟滤波器的传输函数 极点： ， 。 因此，所求数字滤波器的系统函数为 * 1）h(n)的z变换与ha(t)的拉普拉斯变换之间的 关系 设模拟滤波器的单位冲激响应ha(t)经理想采样 后的采样信号为 则 2．s平面与z平面的映射关系 比较序列h(n)的Z变换 ,可得 * 由于 ，将 代 入，得 因此 s与z的 关系 或 * 2）关于z=esT的讨论 令 ， 代入上式，得 因此 讨论： （1）σ=0，r=1，表明s平面虚轴映射为z平面的单 位圆。 * 讨论： （2）σ0，r1；σ0，r1。表明s左半平面映射 为z平面的单位圆内部，而s右半平面则映射为z平 面单位圆外部。 * 讨论： （3）由于ω=ΩT Ω：0 ——±π/T ω：0 ——±π s平面上每一条宽为2π/T的横条，都将重叠的 映射到整个z平面上。脉冲响应不变法不是从s平面 到z平面的简单代数映射关系。 * 3．频率混叠效应 根据时域采样理论，采样序列h(n)的z变换与模 拟信号ha(t)的拉普拉斯变换之间满足如下关系 将 和 代入得数字滤波器的频率响应 数字滤波器的频响是模拟滤波器频响的周期延 拓。 * 图5-11 脉冲响应不变法中的频率混叠现象 * 只有当模拟滤波器的频响是带限的，且带限于 折叠频率以内，即 有 此时，数字滤波器的频率响应重现模拟滤波器 的频率响应而不产生混叠失真。 * 但是，任何一个实际的模拟滤波器，其频响都 不可能是真正带限的，这就不可避免的会产生混叠 失真。 答：当模拟滤波器的频响在折叠频率以上处衰 减越大时，混叠失真就越小。 问题：要求混叠 失真小，对模