切比雪夫II型带通滤波器

NANHUA University课程设计(论文)题目用切比雪夫Ⅱ型 IIR 设计带通（数字频带变换）滤波器学院名称 指导教师 班 学 级 号电 气 工 程 学 院 陈 忠 泽电子 091 [1**********]学生姓名潘星2012 年12月1

1．课程设计的内容和要求（包括原始数据、技术要求、工作要求等） ：1.设计内容： 根据自己在班里的学号20（最后两位）查表一得到一个四位数，0702，由该四位数 索引表二确定待设计数字滤波器的类型：等波纹FIR数字带通滤波器 2.滤波器的设计指标： ⑴阻带下截止频率  sl  e  0.2rad ⑵通带下截止频率  pl  e 50  0.3rad ⑶通带上截止频率  pu  eid 50id 50id 0.7rad⑷阻带上截止频率 su  e  0.8rad ⑸通带最大衰减错误！未找到引用源。 ⑹阻带最小衰减错误！未找到引用源。 其中，错误！未找到引用源。—你的学号的最后两位 3. 滤波器的初始设计通过手工计算完成； 4. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响 （至少选择两种以上合 适的滤波器结构进行分析）； 5. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 6. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 7. 课程设计结束时提交设计说明书。id 502

2．对课程设计成果的要求〔包括图表（或实物）等硬件要求〕 ：滤波器的初始设计通过手工计算完成； 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适 的滤波器结构进行分析）； 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 课程设计结束时提交设计说明书。3．主要参考文献：[1]高息全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8 [2]陈怀琛. 《数字信号处理教程——MATLAB 释义与实现》 [M].北京： 电子工业出版社， 2004.12 [3]张德丰.《详解 MATLAB 数字信号处理》[M].北京：电子工业出版社，2010.6 [4]飞思科技产品研发中心.《MATLAB7 辅助信号处理技术与应用》[M].北京：电子工业 出版社，2005.34．课程设计工作进度计划： 序号1 2 3 4 5 6起 迄 日 期2012.12.26-2013.12.31 2012.12.31-2013.1.3 2013.1.3-2013.1.5 2013.1.5-2013.1.7 2013.1.8-2013.1.9 2013.1.10-2013.1.13工 作 内 容接到题目，搜集资料 整理资料，构思设计方案 手工计算进行滤波器的初步设计 完善初步设计，学习 Matlab 软件操作 通过 Matlab 软件分析设计内容，逐步落实课题目标 上交课程设计，并做细节修改并完成设计主指导教师日期：3年月日

一 手工计算完成切比雪夫 2 型 IIR 带通滤波器初始设计1.设计要求滤波器的设计指标要求为 ⑴ 阻 带 下 截 止 频 率 s1  e  0.5  0.2radid 50 idid 50⑵ 通 带 下 截 止 频 率p1  e  0.5  0.3rad ⑶通带上截止频率, p 2  e 50  0.5  0.7rad ⑷阻带上截止频率 s 2  e 50  0.5  0.8rad ⑸通带最大衰减错误！未找到引用源。, ⑹阻带最小 衰减错误！未找到引用源。 其中 为我学号的后两位。我的学号：[1**********]，所以 =21。id由此计算性能指标（保留小数点后 2 位）得： 阻带下边沿频率为 s1 通带上截止频率 0.48rad ，通带下边沿频率为 p1  0.72rad ，p 2  1.68rad ,阻带上上截止频率 s 2  1.92rad ，通 带 最 大 衰 减 为 Ap  1dB ， 阻 带 最 小 衰 减 为 As  60 dB ， 采 样 频 率f s  2Hz2.手工计算（保留 2 位小数） 3.  s1  2 f s tans 22 0.96 rad / s p1  2 f s tanp12 1.52rad / s p 2  2 f s tanp 22 4.48rad / ss 2  2 f s tans 22 5.72 rad / sBW   p 2   p1  2.96ad / s因为带通滤波器的通带（阻带）边界频率关于中心频率 0 几何对称，即满 足 p1p2  s1s 2  0 ，而现在 p1p 2 s1s 2 ，顾增大  s1 ，使其满足24

条件。则 s1  p1 p 2 s 2 1.19rad / s定义 BW   p 2   p1  1.19ad / s 为通带带宽， 并以此为参考频率对  轴作归 一化处理，即 s1  s1 =0.41  BW, s 2 s 2 =1.93  BW, p1  p1  BW=0.51,p 2 p 2  BW=1.512 0  p1p2 =6.81，归一化的02  p1 p2 =0.77再定义中心频率为02 2   与  的转换关系为：  p 2 -p1 2p 归一化低通滤波器的技术指标 由于p2 -p1  1 ，  p  1 所以，   2 -02 s22 -02 所以 s  =1.53 s 2通带最大衰减 Ap  1dB 阻带最小衰减 As  60dB 由切比雪夫二型滤波器的特征函数H a ( j)21-1  2  1   2C N ( ) p    10 p  1 k11   1965 2 . 0.1 As 10 10.1 A5

archk11 N 9 archs  100.1A 1  0.509s代入归一化系统函数的表达式得模拟低通滤波器的传递函数：Qa p） （ 8 6 0.0026p  0.0724p  0.5508 4  1.5514 2  1.4567 p p 9 8 7 0.0001  0.0028  0.0073  0.0777p 6 p p p 0.0811 5  0.5761p 4  0.2866p 3 p  1.5852 2  0.3138  1.4567 p p2 s 2  0 2.56  s 2 p  p  Bw s 0.4sH a s） Qa p） 2.56s 2 （  （ p0.4s由 H（z） H（s）   a s a2 1  z 1 得数字低通滤波器的传递函数： T 1  z 1利用数字频带转换公式： H BP ( z )  H ( z ) z 1  其中： d1  2k 1 k ， d2  1 k 1 kz 2  d1 z 1  d 2 d 2 z 2  d1 z 1  1s 2   s1    s1 2 其中： k  tan s  tan s 2 ，   cos 0  s 2  s1 2 2 cos 2cos得数字带通滤波器传递函数 H BP z  系数： den=[1 -5.[1**********]162 17.[1**********]18 -41.[1**********]416

75.[1**********]50 -113.[1**********]7 144.[1**********]7 -157.[1**********]1 149.[1**********]3 -123.[1**********]6 88.[1**********]22 -55.[1**********]43 30.[1**********]75 -13.[1**********]75 5.[1**********]342 -1.[1**********]976 0.[**************] -0.[**************]9 0.[**************]05] num=[0.[**************]11 -0.[**************]8 0.[**************]8 -0.[**************]6 0.[**************]2 -0.[**************]3 0.[**************]0 -0.[**************]0 0.[**************]5 8.[1**********]740e-14 -0.[**************]7 0.[**************]2 -0.[**************]4 0.[**************]6 -0.[**************]3 0.[**************]7 -0.[**************]8 0.[**************]2 -0.[**************]11]二 滤波器的结构不同对性能指标的影响1.利用直接型结构构建数字滤波器直接型的结构流图如图所示：0.00411 x(n) y(n)z -1 z -1-0.01628z -1-5.46581z -10.03287 … 0.01628z -117.78253 … -0.07785z -1-0.004110.00841选择 filter structure 选项框中的 Direct-Form I 选项，点击窗口下方的 Import Filter 按钮，构建直接 1 型结构的切比雪夫 2 型带通滤波器，结果如图 所示：7

读图可以得 Direct-Form I 结构的滤波器技术指标（ s1 ， p1 ，p 2 ， s 2 ，单位为 rad/sample；  s ,  p ,单位为 dB）如表所示：性能指标 初始设计指标 Direct-Form I △rad/samples10.15 0.23 0.53 0.61 60 1rad/sample0.138271 0.246126 0.508165 0.581445 58.75763 0.6252307 -0.011729 0.016126 -0.021835 -0.028855 -1.24237 -0.3747693p1p 2s 2sp由上图和上表可以看出： ⑴滤波器幅频曲线在通带和阻带内波动幅度不太均匀。 ⑵阻带最小衰减  s 比初始设计低 1.24237dB，通带最大衰减  p 比初始设计低 0.3747693dB。 ⑶ω pl 和ω su 分别比初始设计高 0.016126 和-0.028855 rad/sample，ω sl 和ω pu 分 别比初始设计低-0.011729 rad/sample和-0.021835 rad/sample，截止频率坡度 较初始设计更加平缓。8

2．利用级联型结构构建数字滤波器直接型的结构流图如图所示：x(n) …y(n)1 j1 j1 j1 j2 j选择 Edit 下拉菜单中点击 Convert to Second-order Sections 选项，将 构建好的 Direct-Form I 结构的切比雪夫 2 型带通滤波器转换为级联滤波器，结 果如图所示:读图可以得级联型结构的滤波器技术指标（ s1 ， p1 ， p 2 ， s 2 ，单位为r a d/sa m p ；  s ,  p ,单位为 dB）如表 2 所示： le性能指标 初始设计指标 0.15 0.23 级联型 0.1478271 0.222168 △ -0.0013729 -0.001632s1p19

p 2s 20.53 0.61 60 1 表20.5220947 0.5981445 58.75001 0.6252766-0.0001053 0.0013445 -1.24999 -0.3747234sp级联结构滤波器对性能指标的影响由图和表中可以看出： ⑴滤波器幅频曲线在通带和阻带内波动比直接型结构滤波器幅频曲线更均匀。 ⑵阻带最小衰减  s 比初始设计值低 1.24999dB，通带最大衰减  p 比初始设计低 0.3747234dB。 ⑶ s 2 比初始设计高 0.0013445 rad/sample； s1 ， p1 ， p 2 ，分别比初始设计低 了 0.0013729 rad/sample、 0.001632 rad/sample和 0.0001053 rad/sample，滤波 器的截止频率坡度较初始设计更加平缓。3、两种滤波器结构对性能指标影响的比较与总结 比较两种结构发现在性能指标方面均有误差。 由于直接型滤波器的系数不是直接决定单个零极点， 不能很好的进行滤波器性能 的控制， 且直接型滤波器的极点对参数的变化过于敏感，从而使得系统的频率响 应对参数的变化也特别敏感， 容易出现不稳定或产生较大误差；而级联型滤波器 每个二阶系数单独控制一对零、极点，有利于控制频率响应。因此直接型误差比 级联型更大，受有限参数字长影响更大，主要表现在直接型的（ s1 ， p1 ，p 2 ，s 2 ， s 和  p 与设计要求的相应性能间的差的绝对值普遍大于级联型。此外，级联型的的幅频响应曲线的通带的波动稳定性要稍好于直接型。 所以， 参数字长有限的情况下，级联结构型滤波器对参数变化的反应要比直 接结构型的更小， 性能指标误差更小， 滤波效果更好， 更能符合设计指标的要求。三 参数字长对性能指标的影响在实际的数字滤波器的设计中， 由于计算机或 DSP 芯片等的字长和存储空间10

有限，所以也只能对设计参数取有限的字长进行设计。然而，如果字长太短，则设计的滤波器误差就会太大，造成滤波效果不佳。下面研究不同参数字长对性能指标的影响。

将计算获得的切比雪夫带通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按钮，此时设计的滤波器幅频响应曲线和性能指标如图所示：

由图和表可以看出，当参数字长取为24位及以上时，幅频曲线失真几乎为零，设计的曲线与要求的曲线几乎完全重合，截止频率与实际要求几乎完全相同，通带最大衰减频率和阻带最小衰减频率与设计指间相差也很小，设计的滤波器的各项性能指标达到了设计要求。 6．结论

总结发现：参数字长越长，设计出的滤波器就越符合设计指标要求，误差越小，稳定性越好，滤波性能越好。当参数字长达到14及位及以上时，设计的滤波器便可达到设计性能指标。

参考文献

[1]高息全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8 [2]陈怀琛.《数字信号处理教程——MATLAB释义与实现》[M].北京：电子工业出版社，2004.12

[3]张德丰.《详解MATLAB数字信号处理》[M].北京：电子工业出版社，2010.6 [4]飞思科技产品研发中心.《MATLAB7辅助信号处理技术与应用》[M].北京：电子工业出版社，2005.3

NANHUA University课程设计(论文)题目用切比雪夫Ⅱ型 IIR 设计带通（数字频带变换）滤波器学院名称 指导教师 班 学 级 号电 气 工 程 学 院 陈 忠 泽电子 091 [1**********]学生姓名潘星2012 年12月1

1．课程设计的内容和要求（包括原始数据、技术要求、工作要求等） ：1.设计内容： 根据自己在班里的学号20（最后两位）查表一得到一个四位数，0702，由该四位数 索引表二确定待设计数字滤波器的类型：等波纹FIR数字带通滤波器 2.滤波器的设计指标： ⑴阻带下截止频率  sl  e  0.2rad ⑵通带下截止频率  pl  e 50  0.3rad ⑶通带上截止频率  pu  eid 50id 50id 0.7rad⑷阻带上截止频率 su  e  0.8rad ⑸通带最大衰减错误！未找到引用源。 ⑹阻带最小衰减错误！未找到引用源。 其中，错误！未找到引用源。—你的学号的最后两位 3. 滤波器的初始设计通过手工计算完成； 4. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响 （至少选择两种以上合 适的滤波器结构进行分析）； 5. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 6. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 7. 课程设计结束时提交设计说明书。id 502

2．对课程设计成果的要求〔包括图表（或实物）等硬件要求〕 ：滤波器的初始设计通过手工计算完成； 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适 的滤波器结构进行分析）； 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 课程设计结束时提交设计说明书。3．主要参考文献：[1]高息全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8 [2]陈怀琛. 《数字信号处理教程——MATLAB 释义与实现》 [M].北京： 电子工业出版社， 2004.12 [3]张德丰.《详解 MATLAB 数字信号处理》[M].北京：电子工业出版社，2010.6 [4]飞思科技产品研发中心.《MATLAB7 辅助信号处理技术与应用》[M].北京：电子工业 出版社，2005.34．课程设计工作进度计划： 序号1 2 3 4 5 6起 迄 日 期2012.12.26-2013.12.31 2012.12.31-2013.1.3 2013.1.3-2013.1.5 2013.1.5-2013.1.7 2013.1.8-2013.1.9 2013.1.10-2013.1.13工 作 内 容接到题目，搜集资料 整理资料，构思设计方案 手工计算进行滤波器的初步设计 完善初步设计，学习 Matlab 软件操作 通过 Matlab 软件分析设计内容，逐步落实课题目标 上交课程设计，并做细节修改并完成设计主指导教师日期：3年月日

一 手工计算完成切比雪夫 2 型 IIR 带通滤波器初始设计1.设计要求滤波器的设计指标要求为 ⑴ 阻 带 下 截 止 频 率 s1  e  0.5  0.2radid 50 idid 50⑵ 通 带 下 截 止 频 率p1  e  0.5  0.3rad ⑶通带上截止频率, p 2  e 50  0.5  0.7rad ⑷阻带上截止频率 s 2  e 50  0.5  0.8rad ⑸通带最大衰减错误！未找到引用源。, ⑹阻带最小 衰减错误！未找到引用源。 其中 为我学号的后两位。我的学号：[1**********]，所以 =21。id由此计算性能指标（保留小数点后 2 位）得： 阻带下边沿频率为 s1 通带上截止频率 0.48rad ，通带下边沿频率为 p1  0.72rad ，p 2  1.68rad ,阻带上上截止频率 s 2  1.92rad ，通 带 最 大 衰 减 为 Ap  1dB ， 阻 带 最 小 衰 减 为 As  60 dB ， 采 样 频 率f s  2Hz2.手工计算（保留 2 位小数） 3.  s1  2 f s tans 22 0.96 rad / s p1  2 f s tanp12 1.52rad / s p 2  2 f s tanp 22 4.48rad / ss 2  2 f s tans 22 5.72 rad / sBW   p 2   p1  2.96ad / s因为带通滤波器的通带（阻带）边界频率关于中心频率 0 几何对称，即满 足 p1p2  s1s 2  0 ，而现在 p1p 2 s1s 2 ，顾增大  s1 ，使其满足24

条件。则 s1  p1 p 2 s 2 1.19rad / s定义 BW   p 2   p1  1.19ad / s 为通带带宽， 并以此为参考频率对  轴作归 一化处理，即 s1  s1 =0.41  BW, s 2 s 2 =1.93  BW, p1  p1  BW=0.51,p 2 p 2  BW=1.512 0  p1p2 =6.81，归一化的02  p1 p2 =0.77再定义中心频率为02 2   与  的转换关系为：  p 2 -p1 2p 归一化低通滤波器的技术指标 由于p2 -p1  1 ，  p  1 所以，   2 -02 s22 -02 所以 s  =1.53 s 2通带最大衰减 Ap  1dB 阻带最小衰减 As  60dB 由切比雪夫二型滤波器的特征函数H a ( j)21-1  2  1   2C N ( ) p    10 p  1 k11   1965 2 . 0.1 As 10 10.1 A5

archk11 N 9 archs  100.1A 1  0.509s代入归一化系统函数的表达式得模拟低通滤波器的传递函数：Qa p） （ 8 6 0.0026p  0.0724p  0.5508 4  1.5514 2  1.4567 p p 9 8 7 0.0001  0.0028  0.0073  0.0777p 6 p p p 0.0811 5  0.5761p 4  0.2866p 3 p  1.5852 2  0.3138  1.4567 p p2 s 2  0 2.56  s 2 p  p  Bw s 0.4sH a s） Qa p） 2.56s 2 （  （ p0.4s由 H（z） H（s）   a s a2 1  z 1 得数字低通滤波器的传递函数： T 1  z 1利用数字频带转换公式： H BP ( z )  H ( z ) z 1  其中： d1  2k 1 k ， d2  1 k 1 kz 2  d1 z 1  d 2 d 2 z 2  d1 z 1  1s 2   s1    s1 2 其中： k  tan s  tan s 2 ，   cos 0  s 2  s1 2 2 cos 2cos得数字带通滤波器传递函数 H BP z  系数： den=[1 -5.[1**********]162 17.[1**********]18 -41.[1**********]416

75.[1**********]50 -113.[1**********]7 144.[1**********]7 -157.[1**********]1 149.[1**********]3 -123.[1**********]6 88.[1**********]22 -55.[1**********]43 30.[1**********]75 -13.[1**********]75 5.[1**********]342 -1.[1**********]976 0.[**************] -0.[**************]9 0.[**************]05] num=[0.[**************]11 -0.[**************]8 0.[**************]8 -0.[**************]6 0.[**************]2 -0.[**************]3 0.[**************]0 -0.[**************]0 0.[**************]5 8.[1**********]740e-14 -0.[**************]7 0.[**************]2 -0.[**************]4 0.[**************]6 -0.[**************]3 0.[**************]7 -0.[**************]8 0.[**************]2 -0.[**************]11]二 滤波器的结构不同对性能指标的影响1.利用直接型结构构建数字滤波器直接型的结构流图如图所示：0.00411 x(n) y(n)z -1 z -1-0.01628z -1-5.46581z -10.03287 … 0.01628z -117.78253 … -0.07785z -1-0.004110.00841选择 filter structure 选项框中的 Direct-Form I 选项，点击窗口下方的 Import Filter 按钮，构建直接 1 型结构的切比雪夫 2 型带通滤波器，结果如图 所示：7

读图可以得 Direct-Form I 结构的滤波器技术指标（ s1 ， p1 ，p 2 ， s 2 ，单位为 rad/sample；  s ,  p ,单位为 dB）如表所示：性能指标 初始设计指标 Direct-Form I △rad/samples10.15 0.23 0.53 0.61 60 1rad/sample0.138271 0.246126 0.508165 0.581445 58.75763 0.6252307 -0.011729 0.016126 -0.021835 -0.028855 -1.24237 -0.3747693p1p 2s 2sp由上图和上表可以看出： ⑴滤波器幅频曲线在通带和阻带内波动幅度不太均匀。 ⑵阻带最小衰减  s 比初始设计低 1.24237dB，通带最大衰减  p 比初始设计低 0.3747693dB。 ⑶ω pl 和ω su 分别比初始设计高 0.016126 和-0.028855 rad/sample，ω sl 和ω pu 分 别比初始设计低-0.011729 rad/sample和-0.021835 rad/sample，截止频率坡度 较初始设计更加平缓。8

2．利用级联型结构构建数字滤波器直接型的结构流图如图所示：x(n) …y(n)1 j1 j1 j1 j2 j选择 Edit 下拉菜单中点击 Convert to Second-order Sections 选项，将 构建好的 Direct-Form I 结构的切比雪夫 2 型带通滤波器转换为级联滤波器，结 果如图所示:读图可以得级联型结构的滤波器技术指标（ s1 ， p1 ， p 2 ， s 2 ，单位为r a d/sa m p ；  s ,  p ,单位为 dB）如表 2 所示： le性能指标 初始设计指标 0.15 0.23 级联型 0.1478271 0.222168 △ -0.0013729 -0.001632s1p19

p 2s 20.53 0.61 60 1 表20.5220947 0.5981445 58.75001 0.6252766-0.0001053 0.0013445 -1.24999 -0.3747234sp级联结构滤波器对性能指标的影响由图和表中可以看出： ⑴滤波器幅频曲线在通带和阻带内波动比直接型结构滤波器幅频曲线更均匀。 ⑵阻带最小衰减  s 比初始设计值低 1.24999dB，通带最大衰减  p 比初始设计低 0.3747234dB。 ⑶ s 2 比初始设计高 0.0013445 rad/sample； s1 ， p1 ， p 2 ，分别比初始设计低 了 0.0013729 rad/sample、 0.001632 rad/sample和 0.0001053 rad/sample，滤波 器的截止频率坡度较初始设计更加平缓。3、两种滤波器结构对性能指标影响的比较与总结 比较两种结构发现在性能指标方面均有误差。 由于直接型滤波器的系数不是直接决定单个零极点， 不能很好的进行滤波器性能 的控制， 且直接型滤波器的极点对参数的变化过于敏感，从而使得系统的频率响 应对参数的变化也特别敏感， 容易出现不稳定或产生较大误差；而级联型滤波器 每个二阶系数单独控制一对零、极点，有利于控制频率响应。因此直接型误差比 级联型更大，受有限参数字长影响更大，主要表现在直接型的（ s1 ， p1 ，p 2 ，s 2 ， s 和  p 与设计要求的相应性能间的差的绝对值普遍大于级联型。此外，级联型的的幅频响应曲线的通带的波动稳定性要稍好于直接型。 所以， 参数字长有限的情况下，级联结构型滤波器对参数变化的反应要比直 接结构型的更小， 性能指标误差更小， 滤波效果更好， 更能符合设计指标的要求。三 参数字长对性能指标的影响在实际的数字滤波器的设计中， 由于计算机或 DSP 芯片等的字长和存储空间10

有限，所以也只能对设计参数取有限的字长进行设计。然而，如果字长太短，则设计的滤波器误差就会太大，造成滤波效果不佳。下面研究不同参数字长对性能指标的影响。

将计算获得的切比雪夫带通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按钮，此时设计的滤波器幅频响应曲线和性能指标如图所示：

由图和表可以看出，当参数字长取为24位及以上时，幅频曲线失真几乎为零，设计的曲线与要求的曲线几乎完全重合，截止频率与实际要求几乎完全相同，通带最大衰减频率和阻带最小衰减频率与设计指间相差也很小，设计的滤波器的各项性能指标达到了设计要求。 6．结论

总结发现：参数字长越长，设计出的滤波器就越符合设计指标要求，误差越小，稳定性越好，滤波性能越好。当参数字长达到14及位及以上时，设计的滤波器便可达到设计性能指标。

参考文献

[1]高息全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8 [2]陈怀琛.《数字信号处理教程——MATLAB释义与实现》[M].北京：电子工业出版社，2004.12

[3]张德丰.《详解MATLAB数字信号处理》[M].北京：电子工业出版社，2010.6 [4]飞思科技产品研发中心.《MATLAB7辅助信号处理技术与应用》[M].北京：电子工业出版社，2005.3