Matlab 课程设计任务书,学生姓名： 管行 专业班级： 通信 0606 指导教师：王晟 工作单位： 信息工程学院 题 目: 利用 MATLAB 仿真软件系统结合双线性变换法设计一个数字巴特沃斯 高通 IIR 滤波器 基础强化训练目的 理论目的 数字信号处理课程设计的目的之一是为了巩固课堂理论学习，并能用所学理论知 识正确分析数字信号处理的基本问题和解释数字信号处理的基本现象。 实践目的 数字信号处理课程设计的目的之二是通过设计具体的各种滤波器掌握滤波器设计 方法和步骤。 训练内容和要求 利用 MATLAB 仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通 IIR 滤波器 在数字信号处理平台上（PC 机MATLAB 仿真软件系统和 TC+编程环境）进行软件仿真设 计，并进行调试和数据分析。课程设计题目由指导教师提供，每人隶属一组完成任务，每 组完成的内容不能雷同（按学号分组） 初始条件 MATLAB 软件 数字信号处理与图像处理基础知识 时间安排： 第 21 周，安排任务（鉴 3-204，7 月 14 日） 第 21 周，仿真设计（鉴主 13 楼计算机实验室） 第 21 周，完成（答辩，提交报告，演示） 指 导 教 师 签 名 ： 2008 年 7 月 14 日 系主任（或责任教师）签名： 年 月 日 1,1,目 录 摘要 . 1 Abstract . 2 1 设计项目要求与说明 . 3 2 系统设计 . 3 2.1 设计思路 . 3 2.2 设计方法对比 . 4 2.3 典型模拟滤波器比较 . 5 2.4 设计步骤 . 5 3 仿真程序的设计与调试 . 6 3.1 数字域指标变换成模拟域指标 . 6 3.2 数字域频率进行预畸变 . 6 3.3 模拟滤波器的设计 . 7 3.4 模拟滤波器变成数字滤波器 . 8 3.5 理论计算数字滤波器的仿真 . 11 4.程序调试中出现的问题 . 12 5. 总结与体会 . 13 参考文献 . 14 附录一 总程序如下 . 15 附录二 设计数字滤波器函数总结 . 17,1,武汉理工大学数字信号处理报告,摘要,此报告重点介绍了用双线性不变法设计 IIR 数字滤波器的基本流程，比较了各种设计 方法的优缺点，总结了模拟滤波器的性能特征。最后以双线性不变法设计了一个高通巴特 沃斯 IIR 数字滤波器，介绍了设计步骤，然后在 Matlab 环境下进行了仿真与调试，实现 了设计目标。 关键词：Matlab 双线性不变法 IIR 数字滤波器 巴特沃斯 高通,2,武汉理工大学数字信号处理报告,Abstract,This report introduced with emphasis of the basic flow of designing the IIR digit filter by the bilinear political reform, compared with each kind of design methods good and bad points, summarized analog filters performance characteristic. Finally design one by the bilinear political reform to pass Butterworth high IIR digit filter, introduced the design procedure, then has carried on the simulation and the debugging under the Matlab environment, has achieved the project objective. Key word: Matlab bilinearity political reform the IIR digital filterpass highbutterworth,武汉理工大学数字信号处理报告,Matlab 课程设计,设计一个数字巴特沃斯高通 IIR 滤波器,1 设计项目要求与说明,课题要求设计一个 IIR 数字滤波器，高通，采用双线性变换法，用巴特沃斯实现，用 matlab 软件对其进行仿真与调试。本设计将先说明用双线性法设计IIR 数字滤波器的原理， 然后写出基于 matlab 的软件设计流程。在对设计进行调试，分析实验数据。,2 系统设计,此部分将详细介绍 IIR 数字滤波器的设计流程，比较各种设计方案的优劣。 2.1 设计思路 IIR滤波器设计的主要方法是先设计低通模拟滤波器，然后转换为高通、带通或带阻 数字滤波器。对于其他如高通，带通，则通过频率变换转换为设计相应的高通，带通等。 在设计的全过程的各个步骤，matlab都提供相应的工具箱函数，使得IIR数字滤波器设计 变得非常简单。总的来说，我的设计思路主要有以下两种： 思路一：从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟 滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。,模拟域,冲激响应不变法,频 率 变 换 双 线 性 变 换 法,图2-1 先频率变换再离散 思路二：先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在 Z域内经数字频率变换为所需类型的数字滤波器。,数字域,双线性变换法,频率变换,图2-2 先离散再频率变换,归一化模拟 低通原型,模拟高,带通 或带阻,数字高,带通 或带阻,归一化模拟 低通原型,数字原型低 通,数字高,带通 或带阻,3,4,武汉理工大学数字信号处理报告 以上两种思路都可以，我最后选择了第一种思路进行设计，即先在模拟域内经频率变 换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型 的数字滤波器。 2.2 设计方法对比 方案一：冲激响应不变法 冲激相应不变法是从时域出发，要求数字滤波器的激响应 h(n)对应于模拟滤波器ha(t) 的等间隔抽样，h(n)=ha(nT) ，其中 T 是抽样周期，因此时域逼近良好。 优点：1 h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好 2 线性相位模拟滤波器转变为线性相位数字滤波器 缺点：1 对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的 频率响应与原来模拟滤波器的频率响应相差很大， 2 不能用来设计高通和带阻滤波器。只适用于限带的低通、带通滤波器 方案二：双线性变换法 双线性变换法是从频域出发，使 DF 的频率响应与 AF 的频率响应相似的一种变换法。 直接使数字滤波器的频率响应，逼近模拟滤波器的频率响应,进而求得 H(z)。 优点：1 避免了频率响应的混迭现象 2 在特定 AF 和特定 DF 处，频率响应是严格相等的，它可以较准确地控制截止 频率的位置。 3 它是一种简单的代数关系，设计十分方便。 缺点：1 除了零频率附近， 与 之间严重非线性，即线性相位模拟滤波器变为非线 性相位数字滤波器 2 要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变 3 对于分段常数型 AF 滤波器，经双线性变换后，仍得到幅频特性为分段常数 的 DF.但在各个分段边缘的临界频率点产生畸变，这种频率的畸变，可通过频率预畸变加 以校正。 方案三：频率变换法 设计思想：1 从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的 模拟滤波器；然后进行双线性变换，由 S 域变换到 Z 域，而得到所需类型的数字滤波器。,武汉理工大学数字信号处理报告 2 先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然 后在 Z 域内经数字频率变换为所需类型的数字滤波器。 综上所述，频率变换法是基于双线性变换法或冲激响应不变法（用于频带变换），此 处选择双线性变换法主要是基于要获得严格的频率响应，以及较准确地控制截止频率的位 置。而由于此种方法是一种简单的代数关系，设计也十分方便。,典型模拟滤波器比较 1、Butterworth 巴特沃斯滤波器：它具有单调下降的幅频特性；即最平幅度。 2、Chebyshev 切比雪夫滤波器：在通带或阻带等波纹，可提高选择性。 Bessel 贝塞尔滤波器：在通带内有较好的线性相位特性。 Ellipse 椭圆滤波器：其选择性相对前三种是最好的。 此处选择巴特沃斯主要是想获得最平稳的幅频响应。而不计较相位特性，而用双线 性变换法也会将线性相位变为非线性相位。 设计步骤 如设计一个数字低通滤波器，其技术指标为： 通带临界频率fp ，通带内衰减小于rp； 阻带临界频率fs，阻带内衰减大于s；采样频率为FS 1 将指标变为角频率,wp=fp*2*pi;ws= fs*2*pi;,（式2-1）,2 将数字滤波器的频率指标Wk由wk=(2/T)tan(Wk/2)转换为模拟滤波器的频率指标wk, 由于是用双线性不变法设计，故先采取预畸变。,（式2-2）,p 2 tg wp ; s 2 tg ws ; T2T2 3 将高通指标转换为低通指标，进而设计高通的s域模型,p, p,1,5,p,4 归一化处理 p 1； p , 1；,（式2-3）,武汉理工大学数字信号处理报告,s,p,1,ss,s, ； ；,（式2-4）,1,10S /10 1,N lg/ lg( ) ；（式2-5） 210 p /10 1S,由式2-3，2-4，2-5计算出N，查表可得模拟低通滤波器的