25g传递函数阵ppt课件.ppt

目录(1/1)目 录•概述•2.1 状态和状态空间模型•2.2 根据系统机理建立状态空间模型•2.3 根据系统的输入输出关系建立状态空间模型 •2.4 状态空间模型的线性变换和约旦规范型•2.5 传递函数阵•2.6 线性离散系统的状态空间描述 •2.7 Matlab问题 •本章小结传递函数阵传递函数阵(1/1(1/1) )2.5 传递函数阵•对于SISO线性定常系统,标量传递函数表达了系统输入与输出间的信息动态传递关系–对于MIMO线性定常系统,将每个输入通道至每个输出通道之间的标量传递函数按序排列成的矩阵函数,•即传递函数阵,可用来表达系统输入与输出间的信息动态传递关系–下面将从状态空间模型出发,分别讨论MIMO系统的•传递函数阵的定义•由状态空间表达式建立系统的传递函数阵,以及•组合系统的状态空间模型和传递函数阵 传递函数阵的定义传递函数阵的定义(1/2)(1/2)2.5.1 传递函数阵的定义•在引入传递函数阵概念之前,需将标量函数拉氏变换的定义扩展到向量函数和矩阵函数–为此,定义对向量函数和矩阵函数的拉氏变换为分别对该向量函数和矩阵函数的各个元素求相应的拉氏变换,–那么我们可对矩阵函数和向量函数进行拉氏变换及其拉氏反变换。

传递函数阵传递函数阵(2/2)(2/2)•对r维输入、m维输出的MIMO系统,若其输入输出的拉氏变换分别为U(s)和Y(s),则系统的输入输出间的动态关系可表示为Y(s)=G(s)U(s)其中G(s)称为传递函数阵,其每个元素为标量传递函数–G(s)的形式为其中Gij(s)描述了第i个输出与第j个输入之间的动态传递关系由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵 (1(1/1)/1)2.5.2 由状态空间表达式求传递函数阵•前面已经介绍了SISO系统从传递函数求系统的状态空间表达式,下面将介绍其逆问题,即怎样从状态空间表达式求系统的传递函数阵–主要内容有:•传递函数矩阵的推导•函数矩阵(sI-A)的逆矩阵的快速计算传递函数矩阵的推导传递函数矩阵的推导(1(1/7)/7)1. 传递函数矩阵的推导•前面已经介绍了SISO系统从传递函数求系统的状态空间表达式,下面将介绍其逆问题,–即怎样从状态空间表达式求系统的传递函数阵•已知MIMO线性定常系统的状态空间表达式为其中x为n维状态向量;u为r维输入向量;y为m维输出向量传递函数矩阵的推导传递函数矩阵的推导(2(2/7)/7)•对上式取拉氏变换,有其中X(s)、U(s)和Y(s)分别为x(t)、u(t)和y(t)的拉氏变换;x(0)为x(t)的在初始时刻t=0的值。

q由于传递函数阵描述的是系统输入输出间动态传递关系,不考虑系统初始条件的影响Ø因此令x(0)=0,于是由状态方程的拉氏变换式有X(s)=(sI-A)-1BU(s)传递函数矩阵的推导传递函数矩阵的推导(3(3/7)/7)——例例2-122-12–将上述X(s)代入式(2-60)中的输出方程,有Y(s)=[C(sI-A)-1B+D]U(s)–因此,可得线性定常连续系统的传递函数阵为G(s)=C(sI-A)-1B+D–若对于输入与输出间无直接关联项(即D=0)的系统,则有G(s)=C(sI-A)-1B•例2-12 求如下系统的传递函数传递函数矩阵的推导传递函数矩阵的推导(4(4/7)/7)——例例2-122-12•解 (1) 先计算逆矩阵C(sI-A)-1B代数余子式代数余子式传递函数矩阵的推导传递函数矩阵的推导( (5/7)5/7)——例例2-122-12•(2) 由传递函数计算公式可得q由于状态变换仅对状态变量进行,保持系统的输入和输出变量及它们间的动静态关系不变因此,有如下结论:Ø描述系统输入与输出间动态传递关系的传递函数阵对状态变换具有不变性函数矩阵函数矩阵(sI-A)的逆矩阵的快速计的逆矩阵的快速计算算(1(1/4)/4)2. 函数矩阵(sI-A)的逆矩阵的快速计算•在传递函数矩阵的许多分析与计算问题中,涉及函数矩阵函数矩阵(sI-A)的逆矩阵(sI-A)-1的计算问题。

–当系统的阶数较高(方阵A的维数较大)时,应用此方法将会遇到多项式矩阵函数的逆阵的计算量大,计算困难问题–下面介绍一种计算 的实用递推算法,其证明可从相关的《矩阵分析》的书籍中找到 函数矩阵函数矩阵(sI-A)的逆矩阵的快速的逆矩阵的快速计算计算(2/4)•设A为nn维的矩阵,则sI-A的逆阵为(sI-A)-1=adj(sI-A)/|sI-A|其中adj(·)为伴随矩阵,|sI-A|为如下特征多项式:|sI-A|=sn+a1sn-1+…+an-1s+an–由伴随矩阵的定义可知,adj(sI-A)可表示为如下多项式矩阵函数adj(sI-A)=sn-1I+sn-2B2+…+sBn-1+Bn其中矩阵Bi(i=2,3,…,n)为nn维的矩阵函数矩阵函数矩阵(sI-A)的逆矩阵的快速的逆矩阵的快速计算计算(3/4)–可以证明,特征多项式的系数ai和伴随矩阵的系数矩阵Bj满足如下递推计算关系式:其中tr(M)表示矩阵M的迹,即M中的主对角线上各元素之代数和函数矩阵函数矩阵(sI-A)的逆矩阵的快速的逆矩阵的快速计算计算(4/4)–递推计算sI-A的伴随矩阵和A的特征多项式后,由(2-63)可计算得(sI-A)-1。

–求得(sI-A)-1,再经拉氏反变换,即可求得矩阵指数函数eAt•上述计算方法采用逐次迭代求得多项式系数ai和伴随矩阵系数矩阵Bj,避免了直接求矩阵函数逆矩阵(sI-A)-1的困难由状态空间表达式求传递函数阵由状态空间表达式求传递函数阵 (1(1/1)/1)2.5.3 组合系统的状态空间模型和传递函数阵•对于许多复杂的生产过程与设备,其系统结构可以等效为多个子系统的组合结构,这些组合结构可以由–并联、–串联和–反馈3种基本组合联结形式表示–下面讨论的由这3种基本组合联结形式构成的组合系统的状态空间模型和传递函数阵并联联结并联联结(1(1/4)/4)1. 并联联结图图2-15 并联联接组合系统方块结构图并联联接组合系统方块结构图并联联结并联联结( (2/4)2/4)•设对应于图2-15示的并联联结的组合系统的两个子系统的传递函数阵为其对应的状态空间表达式分别为其对应的状态空间表达式分别为并联联结并联联结( (3/4)3/4)•从图2-15可知u1=u2=u y1+y2=y故可导出并联联结组合系统的状态空间模型为并联联结并联联结( (4/4)4/4)–因此,由上述状态空间表达式可知,并联组合系统的状态变量的维数为子系统的状态变量的维数之和。

–由组合系统的状态空间表达式可求得组合系统的传递函数阵为•因此,并联组合系统的传递函数阵为各并联子系统的传递函数阵之和串联联结串联联结(1(1/5)/5)2. 串联联结图图2-16 串联联接组合系统方块结构图串联联接组合系统方块结构图q设图设图2-16所示的串联联结的组合系统的两个子系统的传递函所示的串联联结的组合系统的两个子系统的传递函数阵分别和并联连结的结构相同数阵分别和并联连结的结构相同,其对应的状态空间表达式也其对应的状态空间表达式也分别相同分别相同串联联结串联联结( (2/5)2/5)•从图2-16可知 u1=u u2=y1 y2=y因此可导出串联组合系统的状态空间方程为串联联结串联联结( (3/5)3/5)–相应的输出方程为即有串联联结串联联结( (4/5)4/5)–因此,由上述状态空间模型可知,串联连接组合系统的状态变量的维数为子系统的状态变量的维数之和•由串联组合系统的状态空间模型可求得组合系统的传递函数阵为串联联结串联联结( (5/5)5/5)–因此,串联联结组合系统的传递函数阵为串联系统各子系统的传递函数阵的顺序乘积•应当注意,由于矩阵不满足乘法交换律,故在上式中G1(s)和G2(s)的位置不能颠倒,它们的顺序与它们在系统中的串联联结顺序一致。

反馈联结反馈联结(1(1/5)/5)3. 反馈联结图图2-17 反馈联接组合系统方块结构图反馈联接组合系统方块结构图反馈联结反馈联结( (2/5)2/5)•设对应于图2-17所示的反馈联结组合系统的两个子系统的传递函数阵为其对应的状态空间模型分别为其对应的状态空间模型分别为反馈联结反馈联结( (3/5)3/5)•从图2-17可知u1=u-y2 u2=y1=y因此可导出反馈组合系统的状态空间模型为反馈联结反馈联结( (4/5)4/5)即有故故反反馈馈联联结结组组合合系系统统的的状状态态变变量量的的维维数数为为子子系系统统的的状状态态变变量量的的维数之和维数之和反馈联结反馈联结( (5/5)5/5)Y(s)=G0(s)U1(s)=G0(s)[U(s)-Y2(s)]=G0(s)[U(s)-F(s)Y(s)]故[I+G0(s)F(s)]Y(s)=G0(s)U(s)或Y(s)=[I+G0(s)F(s)]-1G0(s)U(s)因此,反馈联结组合系统的传递函数为G(s)=[I+G0(s)F(s)]-1G0(s)q由反馈联结组合系统的联结由反馈联结组合系统的联结图图2-17可知可知反馈联结反馈联结( (6/5)6/5)U(s)=Y2(s)+U1(s)=F(s)G0(s)U1(s)+U1(s)=[I+F(s)G0(s)]U1(s)=[I+F(s)G0(s)]Y(s)故Y(s)=G0(s)[I+F(s)G0(s)]-1U(s)因此,反馈联结组合系统的传递函数又可写为G(s)=G0(s)[I+F(s)G0(s)]-1q按图按图2-17,还可作如下推导还可作如下推导。