关闭

关闭

封号提示

内容

首页 巴特沃斯高通数字滤波器.doc

巴特沃斯高通数字滤波器.doc

巴特沃斯高通数字滤波器.doc

上传者: 李勤力 2017-10-19 评分 5 0 202 28 920 暂无简介 简介 举报

简介:本文档为《巴特沃斯高通数字滤波器doc》,可适用于IT/计算机领域,主题内容包含巴特沃斯高通数字滤波器巴特沃斯高通数字滤波器数字信号处理课程设计老师陈忠泽老师学院电气工程学院班级电子信息工程班学号姓名何依阳二一一年五月目录:一、符等。

巴特沃斯高通数字滤波器巴特沃斯高通数字滤波器数字信号处理课程设计老师陈忠泽老师学院电气工程学院班级电子信息工程班学号姓名何依阳二一一年五月目录:一、IIR数字高通滤波器的设计、数字滤波器的概述、数字滤波器的设计步骤、设计方法、IIR巴特沃斯数字高通滤波器的实例计算二、软件仿真工具及实现环境简介、计算机辅助设计方法、MATLAB直接设计IIR巴特沃斯数字高通滤波器三、滤波器结构对数字滤波器性能指标的影响分析、IIR系统的基本网络结构()直接型()级联型()并联型四、有限字长运算在网络结构中对数字滤波器的影响、运算量化效应对数字滤波器的影响、参数的字长对数字滤波器性能指标的影响、系数量化对数字滤波器的影响五、运用MATLAB的辅助工具FDATOOL画出系统函数图像六、设计心得IIR数字高通滤波器的设计一、IIR数字高通滤波器的设计、数字滤波器的概述所谓数字滤波器是指输入、输出均为数字信号通过数值运算处理改变输入信号所含频率成分的相对比例或者滤除某些频率成分的数字器件或程序。、数字滤波器的设计步骤设计一个IIR数字滤波器主要包括下面个步骤:()确定滤波器要求的规范指标。()选择合适的滤波器系数的计算(如图一流程图所示)。()用一个适当的结构来表示滤波器(实现结构)。()有限字长效应对滤波器性能的影响分析。()用软件或硬件来实现滤波器。本次设计的IIR数字滤波器系数的计算是根据已知的模拟滤波器的特性转换到等价的数字滤波器。本次设计用双线性变换法得到数字滤波器。而且双线性变换法得到的数字滤波器保留了模拟滤波器的幅度响应特性。、设计方法频率变换法设计思想:、从归一化模拟低通原型出发先在模拟域内经频率变换成为所需类型的模拟滤波器然后进行双线性变换由S域变换到Z域而得到所需类型的数字滤波器。、先进行双线性变换将模拟低通原型滤波器变换成数字低通滤波器然后在Z域内经数字频率变换为所需类型的数字滤波器。、IIR巴特沃斯数字高通滤波器的实例计算()设计要求设计一个巴特沃斯数字高通滤波器,要求通带截止频率减不大于dB,阻带截止频率prad,通带衰prad阻带衰减不小于dB()确定数字高通的技术指标:通带截止频率paprad通带最小衰减阻带截止频率dBradp阻带最大衰减asdBphtanshprads,pdBtansrads,sdB。()模拟低通滤波器的技术指标计算如下:对通带边界频率归一化由于本设计的为dB截止频率所以pc,pdB把pssh带入和pphsh求得归一化巴特沃斯低通滤波器的阻带截止频率为:spphshsdB()设计归一化巴特沃斯模拟滤波器G(p)。ksp:spNsp::lgksplgsp所以取N=根据巴特沃斯归一化低通滤波器参数表(见附录)可得归一化模拟低通原型系统函数G(p)为:G(p)pppp(pp)(pp)()利用频率变换公式ppphspphs将G(p)转换成模拟高通HHP(s):Ha(s)G(p)|p把SphphSphSphSSph代入此式可得:sssss用双线性变换法将模拟高通Ha(s)转换成数字高通H(z):令STz,zzzzH(z)Ha(s)|zszzzzz二、软件仿真工具及实现环境简介、计算机辅助设计方法在优秀科技应用软件MATLAB的信号处理工具箱中提供了一整套模拟数字滤波器的设计命令和运算函数方便准确简单容行使得设计人员除了可按上述传统设计步骤快速的进行较复杂高阶选频滤波器的计算、分析外还可通过原型变换直接进行各种典型数字滤波器设计即应用MATLAB设计工具从模拟原型直接变换成满足原定频域指标要求有数字滤波器。、MATLAB直接设计IIR巴特沃斯数字高通滤波器MATLAB编程如下:fs=wp=*fsws=*fsrp=rs=Nn=N,wn=buttord(wp,ws,rp,rs)b,a=butter(N,wn,#high#)freqz(b,a,Nn,fs)MATLAB运算结果如下:N=wn=b=a=三、滤波器结构对数字滤波器性能指标的影响分析、IIR系统的基本网络结构IIR系统的基本网络结构有三种即直接型、级联型和并联型。()直接型N阶差分方程如下:MNy(n)bix(ni)aiy(ni)ii对应的系统函数为H(z)ibziMaiziiiN设M=N=按照差分方程可以直接画出网络结构如下图(a)所示。图中第一部分系统函数用H(z)表示第二部分用H(z)表示那么H(z)H(z)H(z)当然也可以写成H(z)H(z)H(z),按照该式相当于将下图(a)中两部分流图交换位置如下图(b)所示。该图中节点变量w=w因此前后两部分的延时支路可以合并形成如下图(c)所示的网络结构流图我们将下图(c)所示的这类流图称为IIR直接型网络结构。M=N=时的系统函数为bbzbzH(z)azaz接按照H(z)画出直接型结构流图。对照下图(c)的各支路的增益系数与H(z)分母分子多项式的系数可见可以直IIR网络直接型结构由bz和az写出数字滤波器系统函数:zzzzH(z)zzzz有H(z)写出差分方程如下:y(n)y(n)y(n)y(n)y(n)x(n)x(n)x(n)x(n)x(n)直接型网络结构如下图:直接型网络结构图在后面的分析中我们将发现直接型系统对滤波器的性能控制作用不明显极点对系数的变化不灵敏易出现不稳定或较大误差而且运算的累积误差较大。因此在设计时一般不选用直接型。()级联型级联型结构是将系统传递函数H(z)写成具有实系数的二阶节的乘积。将分子和分母多项式分解为各自的根然后将一对复数共轭根(或者任意两个实数根)组合成二阶多项式。在直接型表示的系统函数H(z)中分子、分母均为多项式且多项式的系数一般为实数。现将分子、分母多项式分别进行因式分解得到:H(z)A(C(drrNMrz)z)r上式中,A是常数Cr和dr分别表示H(z)的零点和极点。由于多项式的系数是实数Cr和dr是实数或者是共轭成对的复数将共轭成对的零点(极点)放在一起形成一个二阶多项式其系数仍为实数再将分子、分母均为实系数的二阶多项式放在一起形成一个二阶网络Hj(z)。Hj(z)如下式:Hj(z)上式中Hj(z)jjzjzajzajzHj(z)表示一个一阶或二阶的数字网络的子系统函数每个的网络结构均采用前面介绍的直接型网络结构如下图所示H(z)则由k个子系统级联构成。一阶和二阶直接型网络结构级联型的MATLAB的表示与实现:将数字滤波器系统函数H(z):zzzzH(z)zzzz用MATLAB转换成级联型的程序如下:bz=az=S,G=tfsos(bz,az)freqz(bz,az)说明:S,G=tfsos(bz,az):实现直接型到级联型的变换。B和A分别为直接型系统函数的分子和分母多项式系数向量。返回L级二阶级联型结构的系数矩阵S和增益常数G。MATLAB运算结果如下:S=G=由S和G写出数字滤波器级联型系统函数:zzzzH(z)()()zzzz级联型网络结构如下图:在级联型结构中每一级分子的系数确定一对零点分母的系数确定一对极点因为子网络的零极点也即整体网络的零极点所以整个系统的零极点都可以准确的由每一级的系数来调整和控制这样便于调整滤波器的频率响应性能其灵敏度特性优于直接型和正准型结构。其次级联结构具有最少的存储器。并联级联型网络结构图支路的极点也是整个网络的极点而并联支路的零点却不是整个网络的零点因此并联网络能独立的调整系统的极点位置但不能控制零点。并联结构的灵敏度由于直接型和正准型运算累积误差比级联型小。()并联型如果将级联形式的H(z)展成部分分式形式则得到:H(z)H(z)H(z)Hk(z)对应的网络结构为这k个子系统并联。上式中Hi(z)通常为一阶网络或二阶网络网络系统均为实数。二阶网络的系统函数一般为iizHi(z)iziz式中i、i、i和ii都是实数。如果i=i=则构成一阶网络。由上式其输出Y(z)表示为Y(z)H(z)X(z)H(z)X(z)Hk(z)X(z)上式表明将x(n)送入每个二阶(包括一阶)网络后将所有输出加起来得到输出y(n)。在并联型结构中每一个一阶网络决定一个实数极点每一个二阶网络决定一对共轭极点因此调整极点位置方便但调整零点位置不如级联型方便。另外各个基本网络是并联的产生的运算误差互不影响不像直接型和级联型那样有误差积累因此并联形式运算误差最小。由于基本网络并联可同时对输入信号进行运算因此并联型结构与直接型和级联型比较其运算速度最高。四、有限字长运算在网络结构中对数字滤波器的影响、运算量化效应对数字滤波器的影响在实现数字滤波器时将遇到相乘与求和运算。在定点制运算中每一次乘法运算之后都要作一次舍入(或截留)处理研究定点实现相乘运算的流图如下图所示。图(a)表示无限精度乘积y(n)图(b)表示有限精度乘积y(n),表示舍入运算。采用统计分析方法时可以将舍入误差作为独立噪声e(n)叠加在信号上如图(c)所示。x(n)y(n)x(n)y(n)(n)yx(n)(n)ya理想相乘b实际相乘的非线性流图c统计模型的线性流图定点相乘运算的流图表示显然采用统计分析方法后实际的输出可以表示为:对于舍入处理e(n)的均值为零方差为:e(n),y(n)e(n)yq现在以一个一阶IIR巴特沃斯数字高通滤波器为例来讨论分析方法。表示其输入与输出关系的差分方程为y(n)ay(n)x(n)n式中|a|。它含有乘积项ay(n)这将引入一个舍入噪声其统计分析流图示于下图。一阶IIR滤波器的幅频特性整个系统可以当作线性系统来处理。输出噪声ef(n)是由噪声源e(n)造成的输出误差可由量化误差通过线性系统的方法求得输出噪声。由于e(n)叠加在输入端因此式中h(n)=anu(n)是一阶系统的单位脉冲响应,由式ef(n)e(n)h(n)e(n)anu(n)fh(m)h(l)(ml)femlemh(m)(f为输出噪声方差)和dzH(z)H(z)cjze可求得输出噪声的方差mafem式中H(z)为一阶IIR系统的传递函数即zH(z)za有此可以求得qbbqa(a)(a)(,)fe由此可见字长b越大数字滤波器输出端的噪声越小。、参数的字长对数字滤波器性能指标的影响、系数量化对数字滤波器的影响系统对输入信号进行处理时需要若干参数或者称为系数这些系数都要存储在有限位数的寄存器中因此存在系数的量化效应。系数的量化效应误差直接影响系统函数的零、极点位置如果发生了偏移会使系统的频率响应偏移理论设计的频率响应不满足实际需要。下面分析系数量化误差对极零点位置的影响。如果极零点位置改变了严重时不仅IIR系统的频率响应会发生变化还会影响系统的稳定性。因此研究极点位置的改变更加重要。为了表示系数量化对极点位置的影响引入极点位置灵敏度的概念所谓极点灵敏度,是指每个极点对系数偏差的敏感程度。相应的还有零点位置灵敏度分析方法相同。下面讨论系数量化对极点位置的影响。现分析一个N阶直接型结构的IIR滤波器的传递函数H(z)kbzkM上式表示了一个N阶直接型结构的IIR数字滤波器的系统函数该滤波器的极点kakzkkNB(z)A(z)都在单位圆内聚集在z=附近。系数ak和bk必须用有限位二进制数进行量化存储在有限长的寄存器中经过量化后的系数误差为ar和br量化后的系数用ar和br即rararabbbrrr则实际的系统函数可表示为:(z)HzkbkkzkakkNM显然系数量化后的频率响应已不同于原来设计的频率响应。用直接型结构来实现该滤波器时系数ak和bk都将直接出现在信号流程图中其中ak影响着极点的位置。当由于系数量化误差使一个极点从单位圆内移动到单位圆上或单位圆外时滤波器的稳定性即受到破坏。所以只要有一个系数由于量化产生很微小的误差就有可能使系统失去稳定。反馈支路的阶次N越高使滤波器失去稳定的系数量化误差的绝对值就越小则越容易使滤波器变得不稳定。rarar和设滤波器的传输函数H(z)系数ak和bk经舍入量化后为abbbrrr这里ar和br是量化误差。分母多项式有N个极点用的传输函数为:pi(i=,,…N)表示。这样实际的滤波器H(z)zbkkNMkzkbkk(pipi)ziNMkzkak上式中pi是第i个极点的偏差称为极点误差它应该和各个系数偏差都有关它和各系数偏差的关系用下式表示:pipiak,kakNi,,Npipiaap上式中ak的大小决定着系数k的偏差k对极点偏差i的影响程度。akpipap越大k对i影响也越大。称ak为极点i对系数变化的灵敏度。五、运用MATLAB的辅助工具FDATOOL画出系统函数图像系统函数系数向量经过位二进制舍入量化前后幅频特性和相位响应系统函数系数向量经过位二进制舍入量化前后幅频特性和相位响应系统函数系数向量经过位二进制舍入量化前后零极点特性图系统函数系数向量经过位二进制舍入量化前后零极点特性图由MATLAB运算位和位二进制舍入量化后结果可以看出因为系数的量化使极点位置发生变化系数量化的相对误差(ppq)不到极点位置的相对误差ap达到了。问题不但是数量的变化算一下极点的模可以发现所有根的模都变大了说明量化后的极点离单位圆稍远一些如系统函数系数向量经过和位二进制舍入量化前后的零极点特性图如图(c)、(d)所示。这致使数字高通滤波器的幅频特性降低运用MATLAB的辅助工具FDATOOL画出位和位量化墙后的幅频特性曲线分别如图(a)、(b)中的实线和虚线所示这说明由于系数量化效应使极点位置发生了变化从而改变了原来设计的频率响应特性。另外我们还注意到位二进制舍入量化后极点都在单位圆内部但在系统函数系数向量经过位二进制舍入量化后原来较小的系数相对误差变化较大使滤波器性能偏离原设计指标要求使本来稳定的系统变成了不稳定滤波器。从以上分析可以看出系数的量化效应误差直接影响系统函数的零、极点位置如果发生了偏移会使系统的频率响应偏移理论设计的频率响应不满足实际需要。因此在设计滤波器的时候应尽量选择合理的系数量化以减小量化误差以便设计出较为理想的实用滤波器。六、设计心得:在课设之前我对MATLAB软件特别是滤波器设计中的函数基本上处于一种模糊状态。但是通过与同学不断的交流最后完成了这次课设对滤波器的设计有了比较清楚的了解。这次课设的完成真的不太容易。我的理论基础不是很扎实所以完成起来要比其他同学费劲。在课程设计的过程中我学到了很多东西比如设计滤波器的一些基本函数的用法各种模拟滤波器的特性设计滤波器的一些基本方法。但更为重要的是我对于解决一个问题的思路更加清晰找到了属于自己的方法。当然在设计的过程中不可能避免的遇到了很多问题如刚开始思路比较混乱没有明确的方向。也由于一些基本的概念了解的不够清楚比如模拟频带变换我理解为是做模拟滤波器而老师的课程设计题上有一句话似乎要求是要做数字滤波器。所以在这个问题上走了很多弯路。但通过与同学交流终于弄清楚数字滤波器而不是模拟滤波器。模拟频带变换只不过是一种变化方式或者说一种变化过程。是由模拟低通变为模拟高通再由模拟高通变为数字高通的过程。总的来说这次课程设计让我对MATLAB有了更深刻的了解对数字滤波器的设计流程有了大致的了解掌握了一些设计滤波器的基本方法提高了理论用于实践的能力掌握了更多专业相关的使用知识与技能。同时也暴露了我很多的不足在以后的学习中将进一步发现并克服缺点。

类似资料

该用户的其他资料

心好累人生感悟的句子图片带字[权威资料].doc

[笔记]海底捞的员工福利.doc

2011JMS中国营销科学学术年会征稿通知 2011年《营销科学学报》学术.doc

Java占位符应用.doc

彩钢板吊顶、隔墙、门窗施工.doc

职业精品

精彩专题

用户评论

0/200
    暂无评论
上传我的资料

精选资料

热门资料排行换一换

  • 纯粹理性批判.pdf

  • 《战国中山国灵寿城》1975-1…

  • 航空电子系统.pdf

  • 经方躬行录.pdf

  • 春秋左传研究.pdf

  • 歌德文集.第02卷.威廉.麦斯特…

  • 与神对话.pdf

  • 猛醒的时刻:十一届三中全会十周年…

  • 钛与钛合金.pdf

  • 资料评价:

    / 17
    所需积分:0 立即下载

    意见
    反馈

    返回
    顶部