切比雪夫II型带通滤波器.doc

NANHUA University课程设计(论文)题 目 用切比雪夫型IIR设计带通（数字频带变换）滤波器 学院名称 电 气 工 程 学 院 指导教师 陈 忠 泽 班 级 电子091 学 号 20094470121 学生姓名 潘星 2012年 12 月1课程设计的内容和要求（包括原始数据、技术要求、工作要求等）：1.设计内容：根据自己在班里的学号20（最后两位）查表一得到一个四位数，0702，由该四位数索引表二确定待设计数字滤波器的类型：等波纹FIR数字带通滤波器2.滤波器的设计指标：阻带下截止频率通带下截止频率通带上截止频率阻带上截止频率通带最大衰减阻带最小衰减其中，你的学号的最后两位3. 滤波器的初始设计通过手工计算完成；4. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器结构进行分析）；5. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响；6. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表；7. 课程设计结束时提交设计说明书。 2对课程设计成果的要求包括图表（或实物）等硬件要求：滤波器的初始设计通过手工计算完成；在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器结构进行分析）；在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响；以上各项要有理论分析和推导、原程序以及表示计算结果的图表；课程设计结束时提交设计说明书。3主要参考文献：1高息全 丁美玉.数字信号处理M.西安：西安电子科技大学出版社，2008.82陈怀琛.数字信号处理教程MATLAB释义与实现M.北京：电子工业出版社，2004.123张德丰.详解MATLAB数字信号处理M.北京：电子工业出版社，2010.64飞思科技产品研发中心.MATLAB7辅助信号处理技术与应用M.北京：电子工业出版社，2005.34课程设计工作进度计划：序号起 迄 日 期工 作 内 容12012.12.26-2013.12.31接到题目，搜集资料22012.12.31-2013.1.3整理资料，构思设计方案32013.1.3-2013.1.5手工计算进行滤波器的初步设计42013.1.5-2013.1.7完善初步设计，学习Matlab软件操作52013.1.8-2013.1.9通过Matlab软件分析设计内容，逐步落实课题目标62013.1.10-2013.1.13上交课程设计，并做细节修改并完成设计主指导教师日期： 年 月 日 一 手工计算完成切比雪夫2型IIR带通滤波器初始设计1.设计要求滤波器的设计指标要求为阻带下截止频率 通带下截止频率 通带上截止频率, 阻带上截止频率 通带最大衰减, 阻带最小衰减其中为我学号的后两位。我的学号：20094470121，所以=21。由此计算性能指标（保留小数点后2位）得：阻带下边沿频率为，通带下边沿频率为，通带上截止频率,阻带上上截止频率，通带最大衰减为，阻带最小衰减为，采样频率2. 手工计算（保留2位小数）3. 因为带通滤波器的通带（阻带）边界频率关于中心频率几何对称，即满足，而现在，顾增大，使其满足条件。则 定义 为通带带宽，并以此为参考频率对轴作归一化处理，即=0.41 , =1.93 , =0.51 , =1.51再定义中心频率为 =6.81，归一化的=0.77与的转换关系为：归一化低通滤波器的技术指标由于， 所以， 所以=1.53 通带最大衰减阻带最小衰减 由切比雪夫二型滤波器的特征函数 代入归一化系统函数的表达式得模拟低通滤波器的传递函数：由得数字低通滤波器的传递函数：利用数字频带转换公式：其中：，其中：，得数字带通滤波器传递函数系数：den=1-5.4658111553316217.7825351855018-41.227345342314175.5263259691950-113.614161413897144.624984603567-157.884313921591149.527541834463-123.25796900412688.6819841889922-55.482520642344330.0604635686775-13.92705924515755.44206528681342-1.737801575149760.437024902060527-0.07785646114950190.00841264650675005num=0.00411590260300511-0.01628627686930280.0328729018367188-0.05023194250783660.0687699378963542-0.07823879138932330.0700901835449720-0.05379986881091100.03261278709786558.15796163854740e-14-0.03261278709801870.0537998688110372-0.07009018354506240.0782387913893786-0.06876993789638230.0502319425078477-0.03287290183672180.0162862768693032-0.00411590260300511 二 滤波器的结构不同对性能指标的影响1.利用直接型结构构建数字滤波器直接型的结构流图如图所示： 0.00411 x(n) y(n) -0.01628 -5.46581 0.03287 17.78253 0.01628 -0.07785 -0.00411 0.00841 选择filter structure选项框中的 Direct-Form I选项，点击窗口下方的Import Filter按钮，构建直接1型结构的切比雪夫2型带通滤波器，结果如图所示： 读图可以得Direct-Form I结构的滤波器技术指标（，单位为；,单位为dB）如表所示：性能指标初始设计指标Direct-Form I0.150.138271-0.0117290.230.2461260.0161260.530.508165-0.0218350.610.581445-0.0288556058.75763-1.2423710.6252307-0.3747693由上图和上表可以看出：滤波器幅频曲线在通带和阻带内波动幅度不太均匀。阻带最小衰减比初始设计低1.24237dB，通带最大衰减比初始设计低0.3747693dB。pl和su分别比初始设计高0.016126和-0.028855，sl和pu分别比初始设计低-0.011729和-0.021835，截止频率坡度较初始设计更加平缓。2利用级联型结构构建数字滤波器直接型的结构流图如图所示：x(n) y(n) 选择Edit下拉菜单中点击 Convert to Second-order Sections选项，将构建好的Direct-Form I结构的切比雪夫2型带通滤波器转换为级联滤波器，结果如图所示:读图可以得级联型结构的滤波器技术指标（，单位为；,单位为dB）如表2所示：性能指标初始设计指标级联型0.150.1478271-0.00137290.230.222168-0.0016320.530.5220947-0.00010530.610.59814450.00134456058.75001-1.2499910.6252766-0.3747234表2 级联结构滤波器对性能指标的影响由图和表中可以看出：滤波器幅频曲线在通带和阻带内波动比直接型结构滤波器幅频曲线更均匀。阻带最小衰减比初始设计值低1.24999dB，通带最大衰减比初始设计低0.3747234dB。比初始设计高0.0013445；，分别比初始设计低了0.0013729、0.001632和0.0001053，滤波器的截止频率坡度较初始设计更加平缓。3、两种滤波器结构对性能指标影响的比较与总结比较两种结构发现在性能指标方面均有误差。由于直接型滤波器的系数不是直接决定单个零极点，不能很好的进行滤波器性能的控制，且直接型滤波器的极点对参数的变化过于敏感，从而使得系统的频率响应对参数的变化也特别敏感，容易出现不稳定或产生较大误差；而级联型滤波器每个二阶系数单独控制一对零、极点，有利于控制频率响应。因此直接型误差比级联型更大，受有限参数字长影响更大，主要表现在直接型的（，和与设计要求的相应性能间的差的绝对值普遍大于级联型。此外，级联型的的幅频响应曲线的通带的波动稳定性要稍好于直接型。所以，参数字长有限的情况下，级联结构型滤波器对参数变化的反应要比直接结构型的更小，性能指标误差更小，滤波效果更好，更能符合设计指标的要求。三 参数字长对性能指标的影响在实际的数字滤波器的设计中，由于计算机或DSP芯片等的字长和存储空间有限，所以也只能对设计参数取有限的字长进行设计。然而，如果字长太短，则设计的滤波器误差就会太大，造成滤波效果不佳。下面研究不同参数字长对性能指标的影响。将计算获得的切比雪夫带通IIR数字滤波器的系数输入FDAtool中的filter coefficients工具中，并点击Import Filter按钮，生成数字滤波器。运用Set quantization parameters按钮，在filter arithmetic下拉菜单下选择Fixed Point选项，进入如图7所示的界面。1、参数字长取8位对性能指标的影响将coefficient word length的值改为8，点击下方的Apply按钮，此时设计的滤波器幅频响应曲线和性能指标如图7所示：参数字长取2位时的滤波器幅频响应曲线图图中的虚线为供参考的理想字长下生成的滤波器的幅频响应曲线，图中实线为参数字长取为8位时的滤波器幅频响应曲线。从图中可以看出：字长为8位时，滤波器的各项性能指标离设计指标偏差很大，滤波器失真明显，滤波效果很差，并且系统部稳定远远不能满足设计指标的要求。2参数字长取16位对性能指标的影响将coefficient word length的值改为16，点击下方的Apply按钮，此时设计的滤波器幅频响应曲线和性能指标如图所示： 图9 参数字长取4位时的滤波器幅频响应曲线图由图可以看出，当参数字长取为16位时，幅频曲线已经初具带通滤波器的形态，失真度虽然比8位时有明显改进，但仍很明显，滤波效果很差，通带最大衰减频率和阻带最小衰减频率与设计指间标误差还是很大，并且系统仍然不稳定离设计指标差距仍然很大。 3参数字长取24位对性能指标的影响将coefficient word length的值改为24，点击下方的Apply按钮，此时设计的滤波器幅频响应曲线和性能指标如图所示：性能指标初始设计指标Direct-Form I0.150.134561-0.034790.230.2461260.0161260.530.508165-0.0218350.610.581445-0.0288556058.75763-1.2423710.6252307-0.3747693由图和表可以看出，当参数字长取为24位及以上时，幅频曲线失真几乎为零，设计的曲线与要求的曲线几乎完全重合，截止频率与实际要求几乎完全相同，通带最大衰减频率和阻带最小衰减频率与设计指间相差也很小，设计的滤波器的各项性能指标达到了设计要求。6结论总结发现：参数字