dsp 期末作业.docVIP

Yibin University DSP期末学科设计报告 题 目 第五章5-18、5-19 院 别 物理与电子工程学院 专 业 电子信息科学与技术 学生姓名 梁 杰 　 学 号 110303038 班级 11 级 3 班 2014年1月 6 日  实验目的 了解巴特沃斯方法 掌握脉冲响应不变法设计数字低通滤波器 熟悉各种滤波器特点 实验原理 脉冲不变法原理 脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。它特别适合于对滤波器的时域特性有一定要求的场合。 脉冲响应不变法的核心原理是通过对连续函数ha(t)等间隔采样得到离散序列ha(nT)。令h(n)= ha(nT) ，T为采样间隔。它是一种时域上的转换方法。 一个模拟滤波器的传递函数可以用有理分式表达式表示为： 通过反拉普拉斯变换我们就可以得到它的冲激相应： 脉冲响应不变法就是要保证脉冲响应不变，即： 对上式的冲激相应序列作变换，就可以得到数字滤波器的传递函数： 一般来说，的分母多项式阶次总是大于分子多项式的阶次。假定的没有多重极点，则式就可分解为： 式中，均为复数，是的极点。其拉氏变换为 一般来说，的分母多项式阶次总是大于分子多项式的阶次。数字滤波器的传递函数H(z)经过合并简化，成为一般形式的有理分式传递函数 在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道 按照z=的关系，每一个s平面上宽度为/T的水平条带将重迭映射到z平面上。因此脉冲响应不变法将s平面映射到z平面，不是一个简单的一一对应的关系。对于高采样频率(T小)的情况，数字滤波器在频域可能有极高的增益。为此我们采用 在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的()。同时，它可以保持脉冲响应不变， =。因此，这一方法往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。 2、双线性变换法设计数字滤波器原理 将S平面j轴压缩变换到s1平面j轴上的-/T到/T一段，可以采用以下变换关系： 这样变换到，变换到，可将上式写成 令，，解析延拓到整个s平面和s1平面，可得 再将平面通过以下变换关系映射到z平面，即 从而得到s平面和z平面的单值映射关系为 一般来说，为了使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系，可引入待定常数c， 将代入到上式，可得 在MATLAB中，双线性Z变换可以通过bilinear函数实现，其调用格式为： [Bz，Az]＝bilinear(B，A，Fs)； 根据解题思路，编写源程序如下： %截止频率wc=1Kz，采样频率fs=6.283 KHz fs=6.283*1000; wc=2*pi*1000; [B,A]=butter(3,wc,s); [num,den]=impinvar(B,A,fs); [db,mag,pha,grd,w]=freqz_m(num,den); f=0.5*fs*w/pi; plot(f,db,-);grid; xlabel(以HZ为单位的频率);ylabel(对数幅度/db); title(巴特沃斯低通滤波器); 图一、实验1仿真结果图 2、5-20题用双曲线变换法设计一个满足下面指标的数字带阻巴特沃思滤波器；通带上下