数字滤波器的网络结构课件.ppt
《数字滤波器的网络结构课件.ppt》由会员分享,可在线阅读,更多相关《数字滤波器的网络结构课件.ppt(63页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、数字滤波器的网络结数字滤波器的网络结构构第1页,此课件共63页哦5.1 5.1 概述概述一一.数字滤波器的基本原理数字滤波器的基本原理 是数字信号处理的一个重要技术分支;是数字信号处理的一个重要技术分支;是用一有限精度算法实现的是用一有限精度算法实现的离散时间系统离散时间系统;具有对信号进行具有对信号进行滤波处理滤波处理的功能的功能,即在形形色色的信号中即在形形色色的信号中提取提取所需要的信号所需要的信号,抑制抑制不需要的信号(噪声、干扰);不需要的信号(噪声、干扰);滤波器实质是一种滤波器实质是一种运算过程运算过程(差分方程的计算或卷积的计(差分方程的计算或卷积的计算)算)第2页,此课件共6
2、3页哦二二.数字滤波器数字滤波器的数学模型的数学模型时域离散系统时域离散系统时域中:时域中:)()()(10knyainxbnyNkkMii差分方程差分方程 mmnhmxnhnxny)()()(*)()(卷积公式卷积公式 Z域中:域中:NNMMzazaazbzbbzXzYzH110110)()()(系统函数系统函数第3页,此课件共63页哦 同一种数学模型,具有不同的算法。数字滤波器的工程实同一种数学模型,具有不同的算法。数字滤波器的工程实践要用计算机的硬件和软件来完成,不同的算法会影响系统的践要用计算机的硬件和软件来完成,不同的算法会影响系统的一些实际性能。一些实际性能。(1)(1)计算的效率
3、,即完成整个滤波所需要的乘法和加法次计算的效率,即完成整个滤波所需要的乘法和加法次数;数;(2)(2)需要的存储单元;需要的存储单元;(3)(3)计算机的位数对滤波系数的量化误差;计算机的位数对滤波系数的量化误差;本章着重研究系统结构的等效变换方法。本章着重研究系统结构的等效变换方法。第4页,此课件共63页哦 数字滤波器的实质是一个运算过程,对输入序列数字滤波器的实质是一个运算过程,对输入序列x(n)x(n)进行一定的运算操作,从而得到输出序列进行一定的运算操作,从而得到输出序列y(n)y(n)。三三.基本运算单元的信号流图表示基本运算单元的信号流图表示)()()(10knyainxbnyNk
4、kMii数字滤波器的差分方程模型:数字滤波器的差分方程模型:包含三种运算:加法、数乘、单位延时包含三种运算:加法、数乘、单位延时 第5页,此课件共63页哦 (1)(1)加法加法x x1 1(n)(n)x x2 2(n)(n)x x1 1(n)+x(n)+x2 2(n)(n)(2)(2)数乘数乘a)(nx)(nax (3)(3)单位延时单位延时)(nx)1(nxz-1第6页,此课件共63页哦四四.数字滤波器网络结构的分类数字滤波器网络结构的分类1.1.IIRIIR数字滤波器数字滤波器特点特点 网络结构存在输出到输入的反馈。网络结构存在输出到输入的反馈。单位脉冲响应单位脉冲响应h(n)h(n)是无
5、限长的。是无限长的。Infinite Impulse ResponseInfinite Impulse Responseaz z-1-1x(n)x(n)y(n)y(n)第7页,此课件共63页哦2.2.FIRFIR数字滤波器数字滤波器特点特点 网络结构没有输出到输入的反馈。网络结构没有输出到输入的反馈。单位脉冲响应单位脉冲响应h(n)h(n)是有限长的。是有限长的。Finite Impulse ResponseFinite Impulse Responseaz z-1-1x(n)x(n)y(n)y(n)第8页,此课件共63页哦5.2 IIR5.2 IIR滤波器的基本网络结构滤波器的基本网络结构I
6、IR滤波器的差分方程:滤波器的差分方程:)()()(10knyainxbnyNkkMiiIIR滤波器的系统函数:滤波器的系统函数:NNMMzazaazbzbbzXzYzH110110)()()(第9页,此课件共63页哦一一.直接型直接型 采用加法、数乘、单位延时三种基本运算单元,由差分方程采用加法、数乘、单位延时三种基本运算单元,由差分方程直接实现。直接实现。设设M2,N2,则则)2()1()2()1()()(21210nyanyanxbnxbnxbnyb1z z-1-1x(n)x(n)y(n)y(n)z z-1-1b2b0-a1-a2z z-1-1z z-1-1直接直接I型型第10页,此课件
7、共63页哦b1z z-1-1x(n)x(n)y(n)y(n)z z-1-1b2b0-a1-a2z z-1-1z z-1-1直接直接I型型NNMMzazaazbzbbzXzYzH110110)()()(系统函数系统函数 先实现系统函数先实现系统函数H(Z)的分子部分,后实现其分母部分,的分子部分,后实现其分母部分,然后再将两部分级联起来。然后再将两部分级联起来。第11页,此课件共63页哦b1z z-1-1x(n)x(n)y(n)y(n)z z-1-1b2b0-a1-a2z z-1-1z z-1-1直接直接I型型H1(z)H2(z)()()(21zHzHzH)()(12zHzH交换级联两部分的连接
8、次序,得交换级联两部分的连接次序,得 第12页,此课件共63页哦b1z z-1-1x(n)x(n)y(n)y(n)z z-1-1b2b0-a1-a2z z-1-1z z-1-1直接直接I型型H1(z)H2(z)x(n)x(n)-a1-a2z z-1-1z z-1-1b1z z-1-1z z-1-1b2b0y(n)y(n)第13页,此课件共63页哦x(n)x(n)-a1-a2z z-1-1z z-1-1b1z z-1-1z z-1-1b2b0y(n)y(n)w w1 1(n)(n)w w2 2(n)(n)可共用延时单元可共用延时单元x(n)x(n)-a1-a2z z-1-1z z-1-1b1b2
9、b0y(n)y(n)直接直接II型型第14页,此课件共63页哦2211022110)()()(zazaazbzbbzXzYzHx(n)x(n)-a1-a2z z-1-1z z-1-1b1b2b0y(n)y(n)直接直接II型型系统函数系统函数第15页,此课件共63页哦 直接直接II型与型与H(Z)的关系:的关系:H(Z)H(Z)的分子系数决定非反馈支路的系数,符号相同。的分子系数决定非反馈支路的系数,符号相同。H(Z)H(Z)的分母系数决定反馈支路的系数,但要注意:的分母系数决定反馈支路的系数,但要注意:符号相反符号相反;a0 0要归一化为要归一化为1 1。第16页,此课件共63页哦例例5.3
10、.1 IIR5.3.1 IIR数字滤波器的系统函数数字滤波器的系统函数H(z)H(z)为为12312384112()5311448zzzH zzzz画出该滤波器的直接型结构。画出该滤波器的直接型结构。x(n)y(n)z1z1z1 4811 2454381第17页,此课件共63页哦说 明 yfilter(b,a,x)在在matlab中,输入信号中,输入信号x(n)通过通过IIR直接型滤波器直接型滤波器的输出的输出y(n)可由可由filter函数求得,调用方式为:函数求得,调用方式为:bb0,b1,bM%系统函数的分子系数系统函数的分子系数aa0,a1,aN 系统函数的分母系数系统函数的分母系数
11、第18页,此课件共63页哦二二.级联型级联型1.1.基本原理基本原理NNMMzazaazbzbbzXzYzH110110)()()(将分子、分母多项式进行因式分解:将分子、分母多项式进行因式分解:12211111111111)1)(1()1()1)(1()1()(NkNkkkkMkkkMkkzpzpzdzqzqzcGzH第19页,此课件共63页哦12211111111111)1)(1()1()1)(1()1()(NkNkkkkMkkkMkkzpzpzdzqzqzcGzH式中:式中:c ck k、d dk k 为实数零、极点;为实数零、极点;为共轭复数零、极点;为共轭复数零、极点;kqkqkpk
12、p、G 为增益。为增益。把任意两个实数零、极点组合,每一对共轭零、极点把任意两个实数零、极点组合,每一对共轭零、极点也进行组合,这样,就得到一些实系数的二阶子系统。也进行组合,这样,就得到一些实系数的二阶子系统。第20页,此课件共63页哦H H(z z)分解为二阶子系统的乘积:)分解为二阶子系统的乘积:2211221111)(zzzzzHjjjjj每一个子系统具有如下形式:每一个子系统具有如下形式:kj,2,1)()()()(21zHzHzHGzHk第21页,此课件共63页哦x(n)y(n)z1x(n)y(n)z1z1(a)(b)j 0j 1j 2j 0j 1j 2j 1j 1j 0kj,2,
13、1)()()()(21zHzHzHGzHk每一个二阶子系统都用直接型网络结构实现:每一个二阶子系统都用直接型网络结构实现:2211221111)(zzzzzHjjjjj说 明、0 jj0都为都为1。第22页,此课件共63页哦2.2.用用Matlab求级联形式的各系数求级联形式的各系数NNMMzazaazbzbbzXzYzH110110)()()(2211221111)(zzzzzHiiiii)()()()(21zHzHzHGzHk第23页,此课件共63页哦sos,G=tf2sos(b,a)sos,G=tf2sos(b,a)b b H(z)H(z)分子多项式按分子多项式按z z1 1升幂排列的各
14、项系数;升幂排列的各项系数;a a H(z)H(z)分母多项式按分母多项式按z z1 1升幂排列的各项系数;升幂排列的各项系数;kkkkkkaaabbbaaabbbaaabbbsos210210221202221202211101211101 H H1 1(z)(z)H H2 2(z)(z)H Hk k(z)(z)G G 整个系统归一化以后的增益整个系统归一化以后的增益第24页,此课件共63页哦例例5.3.2 5.3.2 设系统函数设系统函数H(z)H(z)如下式:如下式:12312384112()1 1.250.750.125zzzH zzzz试画出其级联型网络结构。试画出其级联型网络结构。
15、解解:b=8,-4,11,-2;:b=8,-4,11,-2;a=1,-1.25,0.75,-0.125 a=1,-1.25,0.75,-0.125 sos,G=tf2sos(b,a)sos,G=tf2sos(b,a)第25页,此课件共63页哦)5.01)(25.01()3161.131.01)(19.01(8)(211211zzzzzzzHx(n)z12y(n)z14z1 0.3790.25 1.245.264 0.580.190.311.3161第26页,此课件共63页哦注意注意x(n)z12y(n)z14z1 0.3790.25 1.245.264 0.5)5.01)(25.01()264
16、.524.14)(38.02()(211211zzzzzzzH)5.01)(25.01()3161.131.01)(19.01(8)(211211zzzzzzzHx(n)z12y(n)z14z1 0.3790.25 1.245.264 0.580.190.311.3161教材上:教材上:第27页,此课件共63页哦b,a=sos2tfb,a=sos2tf(sos(sos,G)G)说 明由级联型到直接型的转换:由级联型到直接型的转换:3.3.求求信号通过级联形式滤波器的输出信号通过级联形式滤波器的输出y=G y=G*sosfilt(sos,x)sosfilt(sos,x)式中:式中:x x 输入序
17、列输入序列 sos sos 级联各系数级联各系数 G G 级联增益级联增益 y y 系统输出系统输出第28页,此课件共63页哦只需设计一个网孔(即二阶程序只需设计一个网孔(即二阶程序)调整零、极点方便,通过调整零极点配对调整零、极点方便,通过调整零极点配对,可减少量化误差。,可减少量化误差。4.4.级联型网络结构的优点级联型网络结构的优点级联结构中后面的网络输出不会再流到前级联结构中后面的网络输出不会再流到前面,运算误差的累积相对直接型要小。面,运算误差的累积相对直接型要小。第29页,此课件共63页哦三三.并联型并联型1.1.基本原理基本原理将级联型将级联型H(z)H(z)展成部分分式形式,就
18、得到展成部分分式形式,就得到IIRIIR滤波器的并联滤波器的并联型结构。型结构。)()()()(21zHzHzHzHk其中,每个子系统其中,每个子系统Hi(z)Hi(z)具有下列形式:具有下列形式:22111101)(zzzzHiiiii第30页,此课件共63页哦)()()()(21zHzHzHzHk其中,每个子系统其中,每个子系统Hi(z)Hi(z)具有下列形式:具有下列形式:每个子系统也用直接型网络结构实现:每个子系统也用直接型网络结构实现:x(n)x(n)-a1i-a2iz z-1-1z z-1-10iy(n)y(n)1i22111101)(zzzzHiiiii第31页,此课件共63页哦
19、2.2.用用matlabmatlab求并联型各系数求并联型各系数MatlabMatlab中没有直接把中没有直接把 b b、a a 转换成并联型结构的函数。转换成并联型结构的函数。1.r,p,c=residuez(b,a)1.r,p,c=residuez(b,a)求出留数、极点、直接项求出留数、极点、直接项2.B2.Bk k,A,Ak k=residuez(r=residuez(rk k,p,pk k,c,ck k)将共轭的极点留数或二个实极点留数进行合并,形成一个将共轭的极点留数或二个实极点留数进行合并,形成一个二阶子系统。二阶子系统。r rk k,p,pk k,c,ck k 共轭留数、极点、
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字滤波器 网络 结构 课件
限制150内