数字butterworth滤波器的设计与实现.docx

课程设计任务书学生姓名： 专业班级：指导教师：工作单位：题 目: 数字Butterworth滤波器的设计初始条件：Matlab软件；课程设计辅导资料：“Matlab语言基础及使用入门”、“信号与系统”、“数字信号处理原理与实现”、“Matlab及在电子信息课程中的应用”等；先修课程：信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）课程设计时间：1周；课程设计内容：数字Butterworth滤波器的设计，具体包括：基本Butterworth滤波器的设计，数字高通、带通滤波器的设计，冲激响应不变法和双线性变换法的应用等；本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结，按要求进行实验演示和答辩等；课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括：目录；与设计题目相关的理论分析、归纳和总结；与设计内容相关的原理分析、建模、推导、可行性分析；程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结；课程设计的心得体会（至少500字）；参考文献（近5年的、不少于5篇）；其它必要内容等。时间安排：1周附——具体设计内容：低通巴特沃斯模拟滤波器设计。设计一个低通巴特沃斯模拟滤波器：指标如下：通带截止频率：=3400HZ, 通带最大衰减：=3dB 阻带截至频率：=4000HZ，阻带最小衰减：=40dB模拟低通转换为数字低通滤波器已知一模拟滤波器的系数函数为 分别用冲激响应不变法和双线性变换法将转换成数字滤波器系统函数，并图示和 的幅度相应曲线。分别取采样频率＝1000Hz和＝500Hz，分析冲激响应法中存在的频率混叠失真和双线性变换法存在的非线性频率失真等。3.设计一个10阶的Butterworth滤波器，通带为100～250Hz，采样频率为1000Hz，绘出滤波器的单位脉冲响应。指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目 录1 设计总体思路12 低通巴特沃斯模拟滤波器设计22.1 MATLAB提供的设计函数22.1.1巴特沃斯滤波器阶数选择函数22.1.2 完整的巴特沃斯滤波器滤波器设计函数22.2 设计步骤32.3 MATLAB程序及运行结果32.4 结果分析43 模拟低通转换为数字低通滤波器设计63.1冲激响应不变法63.2 双线性变换法63.3 MATLAB提供的设计函数73.4 MATLAB程序及运行结果73.5 结果分析94 设计10阶巴特沃斯滤波器104.1设计原理分析104.2 matlab 设计程序及运行结果104.3 结果分析115 心得体会126 参考文献131 设计总体思路IIR数字滤波器的设计，主要采用间接法，即：首先设计出低通模拟滤波器H(S)；进行频率变换，将其转换为高通、带通、带阻滤波器；再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤波器。对单极点的N阶H(S)用部分展开式：冲激响应不变法取H(S)的单阶极点Sk的指数函数作为H（Z）的极点双线性变换法是用代换H（S）中的S得到H(Z),双线性变换法可完全消除频率混叠失真但存在非线性频率失真，而冲激响应不变法存在混叠失真。在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的滤波器设计函数，这些函数代表了不同类型的逼近函数的滤波器，常用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。 模拟域 冲激响应不变法 频率变换 双线性变换法图1.1 设计整体思路2 低通巴特沃斯模拟滤波器设计2.1 MATLAB提供的设计函数2.1.1巴特沃斯滤波器阶数选择函数（1）求最小阶数N的函数buttord调用格式1：[N, Wn] = BUTTORD(Wp, Ws, Rp, Rs, s)功能：求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn说明：Wp, Ws,为通带截止频率和阻带截止频率如式所示；，Rp, Rs为通带最大衰减 Apass(dB)和阻带最小衰减Astop(dB)如式s对应模拟滤波器。（2）模拟低通滤波器设计函数buttap调用格式：[Z,P,K] = BUTTAP(N)功能：按给定的阶数N设计出巴特沃斯模拟低通滤波器说明：Z,P,K为返回的N阶模拟滤波器的零、极点和增益系数。如：其中k为增