现代控制理论-第2章-状态空间分析法ppt课件.ppt

第一部分 线性系统理论经典控制理论描述系统数学模型的方法： 外部描述：时域内为高阶微分方程、复频域内为输入输 出关系的传递函数； 电机现代控制理论描述系统数学模型的方法： 内部描述：一阶微分方程（时域） 从传递函数的零点、极点分布得出系统定性特性，并已建立起一整套图解分析设计法，至今仍得到广泛成功地应用。电机内部工作原理点击观看 利用状态分析法，对系统进行一系列特性分析，来设计状态反馈和输出反馈。线性系统理论的主要内容：状态空间分析法 线性系统内部特性线性系统状态空间 的综合设计经典控制理论的传递函数描述方法的不足之处： 系统模型为单输入单输出系统； 忽略初始条件的影响； 不包含系统的所有信息； 无法利用系统的内部信息来改变系统的性能。第二章 状态空间分析法 复杂的时变、非线性、多输入多输出系统的问题，需要用对系统内部进行描述的新方法状态空间分析法。本章主要内容本章主要内容 2.1 状态空间描述的基本概念状态空间描述的基本概念 2.2 线性定常连续系统动态方程的建立线性定常连续系统动态方程的建立 2.3 线性定常连续系统状态方程的解线性定常连续系统状态方程的解 2.4 动态方程与传递函数矩阵动态方程与传递函数矩阵 2.5 线性离散系统的动态方程及其解线性离散系统的动态方程及其解2.1 状态空间描述的基本概念状态空间描述的基本概念一一 状态变量状态变量 状态变量指描述系统运动的一组独立（数目最少的）变量。当系统能用最少的n个变量 完全确定系统状态时，则称这个变量为系统的状态变量状态变量。 12,nx txtxt状态变量选取的特点： 状态变量的选取具有非唯一性：即可用某一组， 也可用另一组数目最少的变量。状态变量个数的选取具有唯一性：二二 状态向量状态向量 把描述系统状态的n个状态变量 看作向量X(t)的分量，则X(t)称为状态向量，记以 , 上标T为矩阵转置记号。 若状态向量由n个分量组成，则称n维状态向量。一旦给定 时的初始状态向量 及 的输入向量 ，则 的状态由状态向量 唯一确定。 1nx txt， ， 1Tntxtxtx， ，0tt 0tx0tt tu0tt tx三 状态空间 以n个状态变量作为坐标轴所组成的n维空间称状态空间。系统在任一时刻的状态由状态空间中一点表示，例如二阶系统的状态可由 轴、 轴组成的状态平面（即相平面）中一点表示；三阶系统的状态可由 轴、 轴、 轴组成的三维状态空间中一点来表示；n阶系统的状态则由轴 ， 轴组成的n维状态空间中一点来表示。 初始时刻 的状态 在状态空间中为一初始点；随着时间推移，系统状态在变化，便在状态空间中描绘出一条轨迹，称状态轨迹。1x2x3x1x2x1xnx0t 0tx四四 状态方程状态方程 状态变量的一阶导数与状态变量、输入变量关系的数学表达式称为状态方程。状态方程一阶微分方程或差分方程。状态方程是状态空间分析法的基本数学方程。故系统的状态方程具有非唯一性。 一般形式的状态方程：式中常系数 与系统特性有关。 111,nnnaabb， ，； ， 111 112211221 1222221 1122nnnnnnnnnnnx ta x ta xta xtbuxta x ta xta xtb uxta x taxta xtb u (2-1) 称系统矩阵（系数矩阵，状态阵），称输入矩阵（在此为列矩阵）。12nxxxx12nxxxx111212122212nnnnnnaaaaaaaaaA12nbbbb式中 方程（21）可写成向量矩阵形式： ttu txAxb(2-2)多输入（含p个输入变量）线性定常连续系统的状态方程一般表达式为： 111 1111 11221 1221 121 11 1nnppnnppnnnnnnnppx ta x ta xtb ub uxta x ta xtb ub uxta x ta xtb ub u(2-3)方程(2-3)的向量矩阵形式为 ttxAxbu(2-4)式中u为p维列向量，B为 输入矩阵，或称控制系数矩阵，有np111212122212nnnnnaaaaaaaaaA111212122212ppnnnpbbbbbbbbbB 12nx txttxtx12puuuu五五 输出方程输出方程 系统输出量与状态变量、输入变量关系的数学表达式称输出方程，它是一个代数方程。单输出定常连续系统的输出方程一般形式为：式中常系数 与系统特性有关。可写成向量矩阵形式： 1nccd， ，； ( )y ttdu tcx(2-6) 1 122nny tc x tc xtc xtdu t(2-5)式中 为输出矩阵（在此为行矩阵），d为直接联系输入量、输出量的前向传递（前馈）系数，又称前馈系数。 12ncccc， ，多输入多输出（含q个输出变量）线性定常连续系统的输出方程一般表达形式为：111 1111 111 11 1nnppqqqnnqqppyc xc xd ud uyc xc xd ud u(2-7)其向量矩阵形式为yCxDu(2-8)式中12qyyyy1112112122212nqqqncccccccccC1112121222111ppqqqpdddddddddD12uuupuC为 输出矩阵，D为 前馈矩阵。()qp()qn六六 状态空间表达式状态空间表达式 状态方程、输出方程的组合称为状态空间表达式，简称动态方程。状态空间法用状态方程、输出方程来表达输入输出关系，提示了系统内部状态对系统性能的影响。单输入单输出系统动态方程一般形式为式中 为 维状态向量，u与y为标量，A为n阶方阵，b为 向量，c为 向量，d为标量。xnnp(1)nu yduxAx bcx，（2-9）多输入多输出系统动态方程一般形式为xAxBuyCxDu，（2-10）式中x为 向量，u为 向量，y为 向量，A为n阶方阵，B为 矩阵，C为 矩阵，D为 矩阵。由于 完整地表征了系统动态特性，故有时把一个指定的系统简称为系统 。1n(1)p1qnpqpqnABCD、 、 、ABCD、 、 、动态方程的结构图表示见图21，各方块的输入输出关系规定为: 输出向量（方块所示矩阵）（输入向量）注意到在向量、矩阵的乘法运算中，相乘顺序不允许任意颠倒。图21 动态方程的结构图表示七七 状态空间分析法状态空间分析法 以状态向量描述、分析系统性能的方法称为状态空间分析法。它具有下列优越之处：便于在数字计算机上求解；容易考虑初始条件；能了解并利用处于系统内部的状态信息；数学描述简化；适于描述多输入多输出、时变、非线性、随机、离散等各类系统，是最优控制、最优估计、辨识、自适应控制等现代控制系统的基本描述方法。倒立摆控制系统航天器控制系统机器人控制系统导弹控制系统2.2 线性定常连续系统动态方程的建立线性定常连续系统的动态方程的形式： 一般形式 典型形式xAxBuyCxDu，一 物理系统动态方程的建立实际物理系统动态方程的建立的原则：根据所含元件遵循的物理、化学定律，列写其微分方程；选择可以量测的物理量作为状态变量。例2-1 设机械位移系统如图2-2所示。力F及阻尼器汽缸速度v为两种外作用，给定输出量为质量块的位移x及其速度 、加速度 。图中m、k、f分别为质量、弹簧刚度、阻尼系数。试求该双输入-三输出系统的动态方程。x x 图22 双输入-三输出机 械位移系统点击观看解解 据牛顿力学，故有显见为二阶系统，若已知质量块的初始位移及初始速度，该微分方程在输入作用下的解便唯一确定，故选 和 作为状态变量。设 ，三个输出量为 ，可由微分方程导出下列动态方程： mxfxvkxFxx 12xxxx ，123yx yx yx，12221112232111xxxxfxvkxFmyxyxyfxvkxFm其向量-矩阵形式为xAxBuyCxDu，式中1200001xFkffxvmmmmxuAB123100001001yyykffmmmmyCD例2-2 设空间飞行器如图2-3所示。利用本体坐标系和飞行器本地垂线参考坐标系，试求空间飞行器的动态方程。 图23 空间飞行器点击观看解：空间飞行器相对于参考坐标系进行姿态定向，用一组旋转Euler角即俯仰角、偏航角和滚动角可以唯一的确定飞行器的定向。 利用动力矩定理和动量定理，同时考虑姿态偏移小、速度低、动量小及忽略惯量直积的情况下，可得俯仰轴方向的线性化方程为 ：2222222222010030010001nIuhh 而滚动轴和偏航轴方向的线性化方程为 ：其中11332111113211331133010000000100003000000000001100000100000001nnunnunIIhnhhnh 2311III 3122III 1233III 状态变量图状态变量图 将状态方程中的每个一阶微分方程用图解来表示，即每个一阶微分方程的右端诸项之和，构成了状态变量的导数，经积分可得该状态变量，最终按照系统中各状态变量的关系连接成封闭的图形，便是状态变量图。 它便于在模拟计算机上进行仿真，是向量-矩阵形式状态方程的展开图形，揭示了系统的详细的内部结构。 状态变量图中仅含积分器、加法器、比例器三种元件及一些连接线。积分器的输出均为状态变量。输出量可根据输出方程在状态变量图中形成和引出。例1-1的状态变量图见图1-3，图中 为拉普拉斯算子。 s图2-4 例2-1状态变量图二二 由微分非常或传递函数建立动态方程由微分非常或传递函数建立动态方程 1 实现： 对于给定的系统微分方程或系统传递函数，寻求对应的动态方程而不改变系统的输入-输出特性，称此动态方程是系统的一个状态空间实现。 由于状态变量的选择不唯一，所以状态空间实现也不唯一，最小实现也不唯一。设单输入-输出线性定常连续系统的微分方程具有下列一般形式： 1(2)1210(1)(2)1210nnnnnnnnnyayaya y a yuuuu(2-11)式中y为系统输出量，u为系统输入量，其系统传递函数为 121210121210( )( )nnnnnnnnnN sy ssssG sD su ssasasa sa(2-12) 2 典型实现：1. 能观测标准形实现 设11,1niiiixyxxa yuin，(2-13)其展开式为111112122112223344232212222212231211112nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnxxa yuy a yuxxa yuy a yu a yuxxa yuya yua yua yuxxayuya yua y 3211nnuayu考虑式(2-11)可得10000nxa yua xu 故有状态方程：10021111222111nnnnnnnnnnnnxa xuxxa xuxxaxuxxaxu (2-14)输出方程为nyx(2-15)其向量-矩阵形式为uyxAxbcx，(2-16)式中 0121000100010001naaaaA0121n b123nxxxx x001c式(2-16)所示动态方程，称能观测标准形实现。 2能控标准形式实现将式(2-12)所示传递函数 分解为两部分相串联，并引入中间变量 ，见下图所示图所示 ： G s z s由第一个方块可导出以u作为输入、z作为输出的不含输入导数项的微分方程，由第二个方块可导出系统输出量y可表为z及其导数的线性组合，即 11101110nnnnnza zaz a zuyzzz(2-17)定义如下一组状态变量1120nxz xzxz，(2-18)可得状态方程1223101101121nnnnnxxxxxa za zazua xa xaxu (2-19)输出方程为0 11 21nnyxxx (2-20)其向量-矩阵形式为uyxAxbcx，(2-21)式中0121010000100001naaaa A0001 b121nnxxxx x011nc式(2-21)所示动态方程，称能控标准形实现。o注意到能控、能观测两种准形实现动态方程中诸矩阵存在下列关系： (2-22)式中下标表示能控标准形， 表示能观测标准形，T为转置记号。式(2-22)所示关系称为对偶关系。000TTTcccAAbccb，o对偶关系：3G(s)的对角形实现设D(s)的因式分解为 (2-23)式中为系统的相异实极点，则G(s)可展开成部分分式之和，即12( )()()()nD ssss 1niiiy sN scG su sD ss(2-24) 其拉氏反变换结果有 1niiiiiix tx tu ty tc x t，(2-28) 若令状态变量为 11iix su sins， ，(2-27) 称为极点 的留数。且有ici 1niiicy su ss(2-26) 式中 iiisN scsD s(2-25) 其向量-矩阵形式为uyxAxbcx，(2-31) 1 122nnyc xc xc x(2-30) 展开可得111222nnnxxuxxuxxu(2-29) 式中 1200nA111 b12ncccc4 的约当形实现 G s式中 为 重实极点， 为相异实极点，则 可展成下列部分分式之和，即 1k1kn， ， Gs当 不仅含有相异实极点，还含有相同实极点时，除了可化为能控、能观测标准形实现以外，还可化为约当形实现，其A阵是一个含约当块的矩阵。设 的因式分解为 G s D s 11kknD ssss(2-32) 1112111111nkikki kiy sN sccccG su sD sssss (2-33) 式中 1111111121 !iiisikiisN scsiknD sN sdcsikidsD s， ， ， ，(2-34) 且 1112111111nkikki kiccccy su su su su sssss (2-35) 取状态变量 为1 iixx及 11111212111nkkiii ky sc xsc xsc xsc xs A 则(2-37) 11112111111111kkkiixsu ssxsu ssxsu ssxsu sikns， ，(2-36) 由式(2-36)有 111211213111111111111111kkkkknnxsxssxsxssxsxssxsussxsussxsuss(2-38) 故有状态方程 11 1112121111kkkknnyc xc xx xscxc x111 1112121 1213111 1,1111 1111kkkkkkkknnnxxxxxxxxxxxuxxuxxu输出方程为(2-40) 其向量-矩阵形式为uxAxbycx，(2-41) (2-39) 1111kknccccc111110110knA 00111 b111211kknxxxxxx式中当系统传递函数为应用综合除法有 11101110nnnnnnny sb sbsb sbG su ssasa sa(2-42) 121210111( )( )nnnnnnnnncssbs bN sG sbbsa sas aDs(2-43) 式中 是直接联系输入、输出量的前馈系数， 是严格有理真分式，其系数同综合除法得nb N sD s000111111nnnnnnba bba bbab(2-44) 其动态方程为nuyb uxAxbcx，(2-45) 22yyyTuu 2212y sTsG su sss试确定可控标准形、可观测标准形动态方程；分别确定状态变量与系统输入量 、系统输出量 的关系；画出状态变量图。解 该系统传递函数 为uy G s可控标准形动态方程各矩阵为12cccxxx例2-3 设二阶系统微分方程为2012cA01c b1cTc 与 的关系可如下导出：将 串联分解并引入中间变量 ，有cxuy、 G sz22zzzuyTzz，则 212yTzzT zTzTu令 ，可得所选状态变量为12ccxzxz ，2122222211212112ccxTyT yT uTTxyTyTuTT可观测标准形动态方程各矩阵为01002xxx20012A01T b001c所选状态变量由式(2-13)可得01022xyyTuxy图2-5(a)、(b)分别示出可控及可观测标准形实现的状态变量图。(a) 可控标准形实现的状态变量图图2-5(b) 可观测标准形实现的状态变量图三三 动态方程的线性变换动态方程的线性变换 设系统的动态方程为式中P为非奇异线性变换矩阵。令 xP x(2-47) uyxAxbcx，(2-46) 11AP APbP bccP，式中 u yyxAxbcx，(2-48) 可得变换后系统动态方程为四四 线性化动态方程的建立线性化动态方程的建立实际的物理系统通常含有非线性因素，其一阶微分方程组的一般形式为输出方程的一般形式为1112121212,;,;,npnnnpxfx xx u uuxfx xx u uu， ， ，(2-49) 1112121212,;,;,npqqnpygx xx u uuygx xx u uu， ， ，(2-50) 其向量-矩阵形式为xf x,uyg x,u(2-51) 式中、g的诸元是 和 的某类非线性函数。当 和 限制在工作点或平衡点附近的小偏差范围内工作时，系统非线性方程能足够精确的用线性化方程来描述。1nxx， ，1puu， ，ixiu ,ttxf xu0(2-52) 故有按理想弹道飞行的导弹00(,)x u00(,)x u 即为平衡点。将式(2-51)在 附近展开成台劳级数并略去二次及其以上各项有00TT x ,ux ,u00000ffxxf x ,uxuxu0TT x ,ux ,u000000ggyyg x ,uxuxu可解得0000,0f x u0yg x ,u，(2-53) 式(2-54)和式(2-55)即为线性化动态方程，式中00111122221212,nnTnnnnx ufffxxxfffxxxfffxxx00 x ufAx00111122221212,ppTnnnpx ufffuuufffuuufffuuu00 x ufBuTT 0000 x ,ux ,uggyxuCxDuxu(2-55) 故 TT 0000 x,ux,uffxxuAxBuxu(2-54) 00111122221212,nnTqqqnx ugggxxxgggxxxgggxxx00 x ugCx00111122221212,ppTqqqpx uggguuuggguuuggguuu00 x ugDu2.3 线性定常连续系统状态方程的解 一 齐次状态方程的解 称为齐次状态方程，其解描述的是无控情况下在初始状态作用下系统的自由运动。xAx 无控情况下倒立摆的齐次状态方程的解是不稳定的。1. 幂级数法xAxt 2012kkttttxbbbb0,kx bb，n设的解是的向量幂级数式中都是维向量，则 12120122kkkkttk ttttxbbbA bbbb由对应项系数相等条件，有10200121!kk2kbAbbA bbA b 齐次状态方程的解法是将标量齐次微分方程的解法推广到向量微分方程中去。常见如下两种解法。且 时， 故0t 00 xb 2 21102!k kttttkxIAAAx2 201112!tk kk kkeIttttkk AAAAA定义 (2-56) 0tteAxx(2-57)则众所周知，标量微分方程 的解为 称为指数函数；而式(2-57)所示向量微分方程的解，在形式上是相似的，故把 称为知指数函数；简称矩阵指数。x ax 0 ,atatx te xeteA由于 系 同转移而来， 又有状态转移矩阵之称，并记以 ，即 tx 0 xteA ttetA(2-58)2. 拉普拉斯变换法将 取拉氏变换有xAx 0sssxAxx 10ss xIAx故 110tsxI Ax取拉氏反变换有 (2-59) 11()0tesAIAx(2-60) 不论不论A是否非奇异是否非奇异， 总是非奇异的。由于()sI A 22310kksssssIAAAxx22312222kkkkkksIssssssssss IA AAAA AAA AAII可验证故 一定存在，即 总是非奇异的。系统自由运动的性质完全由状态转移矩 阵确定，故有必要研究 的运算性质 1sIAsIA t t由逆矩阵定义，有21231kksssssIAAAI A(2-61) 二二 状态转移矩阵的运算性质状态转移矩阵的运算性质 意为时刻零的状态即为初始状态。意为时刻零的状态即为初始状态。重写 的幂级数表达式 22112!tkktetttkAIAAA t(2-62) 具有下列性质 t0I1(2-63) 且有 (2-65)0A3 (2-66) 121221tttttt2 tttA A (2-64) 4 11tttt， (2-67)5 00ttttxx(2-68) 000ttxx 10000000tttttttttttxxxxx式(2-68)意为 转移至 的状态转移矩阵为 。由于则 0tx tx0tt6202110tttttt (2-69)2200ttttxx 1100ttttxx 221121100tttttttttxxx由于 又 故式(2-69)成立，意为 至 的状态转移过程可分解为 至 及 至 的分段转移过程。0t1t2t0t1t2t 故式(2-70)成立。 以上七条运算性质与标量指数运算性质相同。7ktk t ktkktk tteeekt AAAA (2-70)由于8 ttttttttttee ee eee ee eA BABBAA BABBAABBAABBA (2-71)9. 设 的状态转移矩阵为 ，则引入非奇异线性变 换 后的状态转移矩阵 为xAx txPx t 11ttteAP PPP(2-72) 10两种常见的状态转移矩阵设 ，即A为对角阵且具有互异元素时，有1ndiagA 1200nttteete 设A为 约当阵，即 ()n n1010A 21221 !2 !0nttttnttttttteteeenteteetntee则有 (2-73)三三 非齐次状态方程的解非齐次状态方程的解 称为非齐次状态方程，其解描述控制作用下系统的强迫运动。uxAxB 有控情况下倒立摆的齐次状态方程的解就会稳定的。非齐次状态方程常见如下两种解法。1 积分法,且有uxAxB设tteeuAAxAxBttttdeeeedt AAAAxAxxxAx由于式中第一项是对初始状态的响应分量，第二项是对控制输入的响应分量，适当地选择控制作用，可使系统响应过程按预定性能指标变化。 ttdeeudtAAxB故积分可得 00tteteudAAxxB 00tttteeudAAxxB故 00ttttudxxB或 (2-75) (2-74) 当取 作为初始时刻，则积分可得0t 000tttteteteudAAAxxB 000tt tttteteudAAxxB 000tttttttudxxB故或(2-77) (2-76) 2拉普拉斯变换法 将ux Ax B取拉氏变换有 0sssu sxAxxB 0ssu sIA xxB 110sssu sxIAxIAB则 故 11110tssu sxIAxIAB取拉氏反变换有取拉氏变换卷积定理 112120tsstdFFff可见与积分法得出相同结果。在此 视为 视为 ，故1sIA 1suFB， 2sF 1200tttetdAxxff(2-78) 00ttttudxxB故(2-79) 例24 试求下列状态方程在 及 作用下的解。 12000Txx x 1u tt1122010231xxuxx 0000ttttBdtdxxxB解解 已知 及 ，为简化计算，可引入变量置换，即令 ，故 00ttttudxxB 1ut 计算 t 11tsIA001102323sssssIA12111311121222212121212sadj ssssssssssssssIAIAIA 2222222tttttttteeeeteeee 2222002201112222ttttttteeeeeetddeeeeeeB 222122221102220222ttttttttttttxeeeeeetxeeeeeex故 一组状态曲线（点击观看）2.4 动态方程与传递函数矩阵一 由动态方程求传递函数矩阵 sssu sxAxB 1ssu sxIAB设系统动态方程为 ,uxAxByCxDu令初始条件为零，取拉氏变换有 (2-80) 称为系统传递函数矩阵，它表示初始条件为零时，输出向量与输入向量拉氏变换式之间的传递关系， 为 矩阵。式(2-81)的展开式为 sG()q p sG 1sssssyCIABD uGu故 (2-81) 1ssGCIABD式中 (2-82) 1111211221222212( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )ppqqqqppy sgsgsgsu sysgsgsgsusysgsgsgsus (2-83) 对于单输入单输出系统，传递函数矩阵蜕变为传递函数。传递函数矩阵是系统的外部描述（即输入输出描述），与状态变量无关。当 时， 是方阵。若 是对角方阵， pq G s G s表示该系统是解耦系统，整个系统同 个独立的子系统组成。 p q( )( ) ( )iiiiy sg s u s则 (2-85) 1,1,iijjy sgsiq jpu s， ；，式中 (2-84) 二 闭环系统中的传递矩阵 zHyHGe定义为开环传递矩阵，它确定偏差向量反馈向量间的传递关系。由于设多输入多输出系统结构图如图如图26所示。图中 分别为系统的输入、偏差、输出、反馈向量， 分别为前向通路、反馈通路传递矩阵。由图可见 u eyz、 、 、GH、yGeG uzGuGHy图26 系统结构图 定义 为偏差传递矩阵，它确定输入向量至偏差向量间的传递关系。e1uy1GHGu 11GHG故式中(2-87) (2-86) 定义 为闭环传递矩阵。由于euzuHGe1ee1HGu u故式中1e1HG(2-88) (2-89) 三 由传递函数矩阵求动态方程给定一传递函数矩阵 ，求系统的 诸矩阵使下式成立，即 (2-90)并称系统( )是 的一个实现。 第二节中研究了单输入单输出系统的传递函数的实现，鉴于多输入多输出系统传递函数矩阵的实现问题比较复杂，这里只研究单输入多输出、多输入单输出系统、前者的传递函数矩阵是列向量，后者的传递函数矩阵是行向量。 sGA B C D、 、 、 1ssC I AB D GA B C D、 、 、 sG 1单输入多输出系统传递矩阵的实现单输入多输出系统传递矩阵的实现 系统由q个独立的子系统组成，传递矩阵 为 11111( )qqqqqgsdgsdgsssgsdgsdgs GdG(2-91) 设最小公分母为 1110nnnD ssasa sa(2-92) 则 的一般形式为 sG 11,1111012,121201,1101nnnnnq nqqsssssD sssG(2-93) 将 进行串联分解，并引入中间变量 ，令 sGz (1)11nnz su sxzxzD s，(2-94) 式中 是 个子传递函数的公共部分。由于单输入，其输入矩阵为一列；由于多输出，其输出矩阵 行，故可导出可控标准形实现的状态方程： 1D sqq01210100000100000101nuuaaaa xxAxb(2-95) 其向量矩阵形式为110111,110,1nqqqq njydxuuydyCxd(2-97) 各子系统的输出方程均为 及其各阶导数的线性组合， 有 z11,111 210 11,11 20 1nnqq nnqqqyxxxd uyxxxd u(2-96) 3.多输入单输出系统传递矩阵的实现多输入单输出系统传递矩阵的实现 系统由p个独立子系统组成，传递矩阵为 11111pppppsg sgsdg sdgsddg sgssGdG(2-98) 同理设 的最小公分母为式(2-92)，则 1,pgsgs， 1111,11101,1101nnpnnppsddssssD sG(2-99) 由于多输入，其输入矩阵为 列；由于单输出矩阵为一行，故可导出可观测标准形实现的动态方程，即p1020010112112112222321,12,1,11000100010001pppnnp npnuauauauaxxAxBu(2-100) 120001pydddxu cxdu(2-101) 例2-5 试求下列单输入-双输出系统传递函数矩阵的可控标准 形实现及对角形实现。 31241ssssssG解解 由于单输入，其输入矩阵为一列；由于双输出，其输出矩阵 为二行。将 化为严格有理真分布： 3301212133111ssssssssss GdG sG求 各元的最小公分母 ： sG D s 21232D sssss故 301321sssD s G则可控标准形动态方程为010231uu xAxbx310163uu ycxdx由 确定系统的极点为-1、-2，它们是对角形状态阵的元素。由于输入矩阵只有一列， 式中所示 阵是元素全为1的列向量，故对角形实现的状态方程为 0D s b101021uu xAxbx其输出矩阵有二行，式(1-31)示出 阵由对应极点的留数确定。 在极点-1的留数 为c sG1c 11132113232sssssss cG 在极点-2的留数 为 sG2c 22231123201ssssssscG故对角形实现的输出方为12210301uccuu ycxdxdx 2.5 线性离散系统的动态方程及其解 有些系统是完全离散的，基输入量、中间传递的信号、输出量等都是离散信息；有些系统是局部离散的，其输入量、受控对象所传送的信号、输出量等都是连续信息。 计算机传送处理离散信号，这时，连续部分在采样点上的数据才是有用信息，故需将连续部分离散化，为研究方便，不论完全的或局部的离散系统，均假定采样是等间隔的；在采样间隔内，其变量均保持常值。导弹传输系统 经典控制理论中，线性离散系统的动力学方程是用标量差分方程或脉冲传递函数来描述的。而对于离散系统或离散化系统的状态模型，它们的解法则是一样的。一 由差分方程或脉冲传递函数建立动态方程单输入-单输出线性定常离散系统差分方程的一般形式为 1101101111nnny knay kna y ka y kb u knb u knbu kb u k (2-102) 式中 表示 时刻， 为采样周期； 为 时刻的输出量， 为时刻 的输入量； 是与系统特性有关的常系数。kkTT y k u k,iiabkTkT初始条件为零时，离散函数的 变换关系为z ,izy kyzzy kiz yz (2-103) 对式(1-102)进行 变换，整理为z 1110111011101110nnnnnnnnnnnnnny zb zbzb zbG zu zzaza zaN zzzbbzaza zaD z(2-104) 在 的串联分解实现中，引入中间变量 ，则有 /NzD z Q z 11101110nnnnnz Q za z Q zazQ zaQ zu zy zz Q zzQ zQ z(2-105) 而连续系统动态方程的建立方法，对离散系统是同样适用的。定义如下一组状态变量： 12111nnnxzQ zxzzQ zzxzxzzQ zzxz(2-106) 于是 0 11 21nnnnzQ zzx zax zax za x zu z 0 11 21nny zx zx zx z(2-107) (2-108) 由式(2-106)式(2-108)可得动态方程 122310 11 211111nnnnnx kx kx kx kxkx kx kax kax ka x ku k (2-110) 利用 反变换关系 11,1iiiiZx zx kZzx zx kz(2-109) 011nnykxkb uk (2-113) 其向量-矩阵形式为 1122110121( )010010( )00101000011( )01( )1nnnnnx kx kx kxku kxkxkaaaaxkx k (2-112) 0 11 21nny kx kx kx k(2-111) 与连续系统的情况相类似，单输入-单输出线性定常离散系统的动态方程的形式可推广到多输入-多输出系统，有简记为 1kku kkkdu kxGxhycx(2-114) 1kkkkkkxG xH uyC xD u (2-115) 二 定常连续动态方程的离散化已知定常连续系统状态方程xAxBu在 及 作用下的解为 0tx u t 000tttttttdxxBu (2-116) 令 ，有 ；令 ，有 在 有 常数，于是其解为0tkT 0tkTkxxx1tkT11kTkxx,1tk k 1kkuu 111kTkTkTkkTdkxxBu (2-117) 11kTkTTkTdGB记 (2-118) 为便于计算 ，引入下列变量置换，即令 TG1kT 00TTTBdBdG(2-119) 则 故离散化系统状态方程为 1kTkTkxxGu (2-120) 式中 由连续系统的状态转移矩阵 导出，有 T t t TTt (2-121) 离散化系统输出方程为 kkkyCxDu (2-122) 令式(2-120)中 ，可得到 时刻的状态，即0,1,1kk,2 ,TTkT三 定常离散动态方程的解 232121100:1001:2110012:32200121:11001210kkkkkkiikTTkTTTTTTkTTTTTTTTkkkTkTkTTTTTTkTkTTTi xxGuxxGuxGuGuxxGuxGuGuGuxxGuxGuGuGuGuxGu(2-123)式(2-123)为离散化状态方程的解，又称为离散状态转移方程。当 时，有 00,11iiku， ， 000kkTkTkxxxx (2-124) 式中 称为离散化系统状态转移矩阵。 k离散化系统输出方程为 1100kkkiikku kCTTTik yCxDxCGuDu(2-125) 对于离散动态方程式(2-125)，其解为 11011000kkkiikkkiikxikxik xGGHuyCGCGHuDu(2-126) 式中 表示 个 自乘。kGkG例2-7 求下列连续状态方程的离散化状态方程。设采样周期 秒。1T 010231 xxu解解 先求该连续系统的状态转移矩阵 :由例2-4已知 为 t t 22222222tttttttteeeeteeee 10.60040.23250.46510.0972i TTt 故 计算 ： G T 22200211222TTTTTTeeeeTt BddeeeeG 10.19980.2325TTG故 得离散化状态方程为 0.60040.23250.199810.46510.09720.2325kTkTkkkxxGuxu