matlab中tf2par函数,《数字信号处理教程——MATLAB释义与实现》第六章.ppt

《数字信号处理教程——MATLAB释义与实现》第六章

第六章 数字滤波器结构 数字滤波器结构 6.1 概述 6.2 滤波器的信号流图结构和解法 6.3 FIR基本网络结构 6.4 IIR基本网络结构 6.5 一些其他类型的滤波器结构 6.6 状态变量分析法 6.7 各种结构形式的比较 6.1 概述 任何线性时不变集总参数离散系统都可用下列三种形式之一来表述其输入输出关系：差分方程，卷积公式和系统函数。这三种方法是等价的，从任何一个都能推导出其它两个。 即使用同一类表示方法，也存在着许多等价的算法结构，滤波器的工程实现要用计算机的硬件或软件。不同的算法要求不同的实现结构，也会影响系统的某些实际性能。 概述 滤波器实现中需要考虑许多问题，如 (1)。计算的效率：即完成整个滤波所需要的乘法和加法次数； (2)。需要的存储量； (3)。滤波器系数的量化影响； (4)。运算中的舍入和截断误差、饱和和溢出； 不同的算法在满足上述的要求方面是有很大差别的。在研究FFT算法时已经看到过，同样的DFT运算，改变计算结构对提高计算效率和节省存储量能起多大的作用。 6.2 滤波器的信号流图结构 信号流图由节点和有向支路组成。每个节点表示一个信号，有向支路表示信号流动方向和算法。在数字信号处理中只有右图所示的三种基本算法，即倍率、加法和单位时延。没有输入箭头的节点称为输入节点；没有输出箭头节点称为输出节点。 滤波器的信号流图结构 例6.2.1 右图表示一个信号流图的实例，可以依次对每个节点列出方程。五个节点可以列出五个联立方程： 滤波器的信号流图结构 用信号流图可以简明地表示系统的运算情况，同时也很容易列写它的方程。以下均用信号流图表示网络结构。 信号流图中的环路，是指信号顺着箭头流的方向能回到起点处的整条路径。有环路意味着存在反馈。上图中x1-x2-x3-x4-x1, x1-x3-x4-x1和x3-x4-x3都是环路。 不同的信号流图代表不同的运算方法，对于同一个系统函数，可以有很多种信号流图与之相对应。 基本信号流图结构 满足以下条件的信号流图，是在计算机上可实现的。称为基本信号流图。 (1)信号流图中所有支路进行的都是基本运算，即支路增益是常数或者是z -1； (2)信号流图中如果有环路，则环路中必须存在延迟支路，否则将出现信号xi(n)的计算要要依赖于xi(n)值的悖论状况，是数字系统无法实现的。没有延迟支路的环路称为代数环，代数环是构成系统时必须避免的。 (3)节点和支路的数目是有限的。 信号流图方程的解法 当系统内没有环路时，只要靠代入法消去中间变量就能解出。如果系统内存在环路，就要进行等式两端的移项合并，有些麻烦。当回路较多时，那是相当繁琐和容易出错的。 求解时，因变量可以用信号序列x(n)，也可以用其z变换X(z)。自变量宜用可实现的算子z-1。 一步求解的方法有梅森公式，它利用信号流图的拓朴关系。但计算很繁，不实用。 由于应用了MATLAB，本书作者创新地推出了它的计算机解法。 信号流图的代数解法例6.2.1 求图6.2.2的系统函数H(z)=Y(z)/U(z)。 解 对(6.2.1)式进行z变换，得到： 经过消元移项等代数处理，得到： 信号流图的计算机解法 设信号流图中有Ki个输入节点，K个中间节点和输出节点，分别代表输入信号ui(i=1,2,…Ki)和系统状态xj(j=1,2,…K)。令q= z-1，则信号流图中的任何节点变量xj可以表为输入变量ui和和其它节点变量xj的线性组合： 用矩阵表示可写成： 信号流图的计算机解法 其中：X = [x1, x2 , …, xK]T 为K阶状态变量列向量, U=[u1, u2, … , uKi]T 为Ki阶输入列向量，Q为K×K阶的连接矩阵，P为K×Ki阶的输入矩阵。 令W=X/U为以U向量为输入，X向量为输出的系统函数，把QX移到等式左方，不难得到： 信号流图的计算机解法 其中W为K×Ki的矩阵，也称传递矩阵。若只有一个输入，则W为K阶单列向量，表示输入为U，输出为X(其中含K个变量)的K个传递函数组成的向量。 这个简明的公式就等价于信号流图中的梅森公式。只要写出P和Q，任何复杂系统的传递矩阵都可用这个简单的式子求得。在没有计算工具的时候，计算(I-Q) –1涉及一个阶次很高的矩阵的求逆运算，所以虽然简明却不好使用。现在有了科学计算软件，这就不是问题了。 信号流图的计算机解法 用来解此题的计算软件必须有公式推导功能：由于待求的传递函数中存在着一个无法赋值的符号变量q= z-1，它无法为仅会做数值运算的软件所接受。幸而MATLAB已经引进了符号推理的功能。只要在程序的开始，规定某些变量是符号变量，系统就不会去追究它的值，而在运算中始终保留这个符号。 定义符号变量使用的MAT