第1章控制系统的状态空间表达式.ppt

现代控制理论,张涛,自动化专业学位课程,华北科技学院自动化系,1.1 状态变量及状态空间表达式,1.3 状态变量及状态空间表达式的建立(一),1.2 状态变量及状态空间表达式的模拟结构图,1.5 状态矢量的线性变换(坐标变换),1.4 状态变量及状态空间表达式的建立(二),1.8 时变系统和非线性系统的状态空间表达式,1.6 从状态空间表达式求传递函数阵,1.7 离散时间系统的状态空间表达式,1.1 状态变量及状态空间表达式,1.1.1 状态变量,状态变量确定系统状态的最小一组变量，如果知道这些变量在任意初始时刻t=t0 的值以及 tt0 的系统输入，便能够完整地确定系统在任意时刻 t 的状态。,状态动态系统的状态是一个可以确定该系统行为的信息集合。这些信息对于确定系统未来的行为是充分且必要的。,状态变量的特点：状态变量的选择不唯一，但相互独立，各个变量线性无关，其个数为微分方程的阶数，也为系统中独立储能元件的个数多组状态变量间存在某种非奇异变换关系状态变量可选任何变量，但通常选易测量或易观察的量，以便满足实现状态反馈、改善系统性能的需要,1.1 状态变量及状态空间表达式,如果 个状态变量用 表示，并把这些状态变量看作是矢量 的分量，则 就称为状态矢量，记作：,1.1.2 状态矢量,1.1.3 状态空间,以状态变量 为坐标轴所构成的 维空间，称为状态空间。,1.1.4 状态方程,反映系统内部状态变量和输入变量间因果关系的一阶微分（差分）方程组称为系统的状态方程。,用下图所示的 网络，说明如何用状态变量描述这一系统。,图1.1,此系统有两个独立储能元件，有两个状态变量，分别为uc和i，容易写出两个含有状态变量的一阶微分方程组：,式(1)就是图1.1系统的状态方程，式中若将状态变量用一般符号 表示，即令 并写成矢量矩阵式，则状态方程变为：,或,1.1.5 输出方程,在指定系统输出的情况下，该输出与状态变量间以及输入量之间的函数关系式，称为系统的输出方程。如在图1.1系统中，指定 作为输出，输出一般用y表示，则有：,式（3）就是图1.1系统的输出方程，它的矩阵表示式为：,式中,在经典控制理论中，用指定某个输出量的高阶微分方程来描述系统的动态过程。如上图一所示的系统，在以 作输出时，从式(1)消去中间变量i，得到二阶微分方程为：,其相应的传递函数为：,（6）,（5）,回到式（5）或式（6）的二阶系统，若改选 和 作为两个状态变量，即令 则得一阶微分方程组为：,1.1.6 状态空间表达式,1.1.6 状态空间表达式,即：,（8）,设单输入一单输出定常系统，其状态变量为 则状态方程的一般形式为：,输出方程式则有如下形式：,用矢量矩阵表示时的状态空间表达式则为：,多输入一多输出系统状态空间表达式的矢量矩阵形式为：,式中，x和A分别为n维状态矢量和nn系统矩阵；,（9）,（10）,1.1.7 状态空间表达式的系统框图,和经典控制理论相类似，可以用框图表示系统信号传递的关系。对于式(9)和式(10)所描述的系统，它们的框图分别如图a和b所示。,为了简便，下面除特别申明，在输出方程中，均不考虑输入矢量的直接传递，即令D=0。,1.2 状态变量及状态空间表达式的模拟结构图,状态空间表达式的框图可按如下步骤绘制：积分器的数目应等于状态变量数，将它们画在适当的位置，每个积分器的输出表示相应的某个状态变量，然后根据所给的状态方程和输出方程，画出相应的加法器和比例器，最后用箭头将这些元件连接起来。,对于一阶标量微分方程：,它的模拟结构图示于下图,再以三阶微分方程为例：,将最高阶导数留在等式左边，上式可改写成,它的模拟结构图示于下图,同样，已知状态空间表达式，也可画出相应的模拟结构图，下图是下列三阶系统的模拟结构图。,下例是二输入二输出的二阶系统的模拟结构图。,可见较复杂，所以多输入多输出系统的结构图多以矢量结构图的形式表示。,1.3 状态变量及状态空间表达式的建立(一),状态空间表达式一般可以从三个途径求得：一是由系统框图来建立，即根据系统各个环节的实际连接，写出相应的状态空间表达式；二是从系统的物理或化学的机理出发进行推导；三是由描述系统运动过程的高阶微分方程或传递函数予以演化而得。,1.3.1 从系统框图出发建立状态空间表达式,该法是首先将系统的各个环节，变成相应的模拟结构图，并把每个积分器的输出选作一个状态变量 其输入便是相应的 然后，由模拟图直接写出系统的状态方程和输出方程。,一、由系统框图建立状态空间描述,例1-1：系统框图如下：,关键：将积分部分单独表述出来，对结构图进行等效变换,等效变换如下：,图中有三个积分环节，三阶系统，取三个状态变量如上图（选择积分环节后的变量为状态变量）：则有：,写成矩阵形式：,系统,1.3.2 从系统的机理出发建立状态空间表达式,一般常见的控制系统，按其能量属性，可分为电气、机械、机电、气动液压、热力等系统。根据其物理规律，如基尔霍夫定律、牛顿定律、能量守恒定律等，即可建立系统的状态方程。当指定系统的输出时，很容易写出系统的输出方程。,步骤：1)根据系统的机理建立相应的微分方程或差分方程；2)选择有关的物理量作为状态变量；3)导出状态空间表达式。,电路如图所示。建立该电路以电压u1,u2为输入量，uA为输出量的状态空间表达式。,例1-2,解：1)选择状态变量 两个储能元件L1和L2，可以选择i1和i2为状态变量，且两者是独立的。,2）根据基尔霍夫电压定律，列写2个回路的微分方程：,整理得：,3）状态空间表达式为：,例1-3试列出在外力f作用下，以质量 的位移 为输出的状态空间描述。,解：该系统有四个独立的储能元件。取状态变量如下：,则有：,及：,将所选的状态变量,例1-4 建立右图所示机械系统的状态空间表达式（注：质量块 m 的重量已经和弹簧 k 的初始拉伸相抵消）,根据牛顿第二定律,即：,选择状态变量,则：,机械系统的系统方程为,该系统的状态图如下,例1-5 建立电枢控制直流他励电动机的状态空间表达式,电枢回路的电压方程为,系统运动方程式为,（式中，为电动势常数；为转矩常数；为折合到电动机轴上的转动惯量；为折合到电动机轴上的粘性摩擦系数。）,可选择电枢电流 和角速度 为状态变量，电动机的电枢电压 为输入量，角速度 为输出量。,状态空间表达式,状态图如下：,1.3.5 传递函数（矩阵）描述和状态空间描述的比较,1）传递函数是系统在初始松弛的假定下输入-输出间的关系描述，非初始松弛系统，不能应用这种描述；状态空间表达式即可以描述初始松弛系统，也可以描述非初始松弛系统。,2）传递函数仅适用于线性定常系统；而状态空间表达式可以在定常系统中应用，也可以在时变系统中应用。,3）对于数学模型不明的线性定常系统，难以建立状态空间表达式；用实验法获得频率特性，进而可以获得传递函数。,4）传递函数仅适用于单入单出系统；状态空间表达式可用于多入多出系统的描述。,5）传递函数只能给出系统的输出信息；而状态空间表达式不仅给出输出信息，还能够提供系统内部状态信息。,综上所示，传递函数（矩阵）和状态空间表达式这两种描述各有所长，在系统分析和设计中都得到广泛应用。,状态空间法具备如下优点：（1）在数字计算机上求解一阶微分方程组或者差分方程组，比求解与它相当的高阶微分方程或差分方程要容易。,（2）状态空间法引入了向量矩阵，大大简化了一阶微分方程组的数学表示法。,（3）在控制系统的分析中，系统的初始条件对经典法感到困难的问题，采用状态空间法就迎刃而解了。,（4）状态空间法能同时给出系统的全部独立变量的响应，不但反映了系统的输入输出外部特性，而且揭示了系统内部的结构特性，既适用单输入单输出系统又适用多输入多输出系统。,（5）状态空间法可利用计算机进行分析设计以及实时控制，所以可应用求解大量的非线性系统、时变系统、随机过程和采样系统。,（6）利用现代空间法进行系统综合时，是非常有利的。,1.4 状态变量及状态空间表达式的建立(二),对于给定的系统微分方程或传递函数，寻求对应的状态空间描述而不改变系统的输入-输出特性，称此状态空间描述是系统的一个状态空间实现。由于所选状态变量不同，其状态空间描述也不同，故其实现方法有多种。,考虑一个单变量线性定常系统，它的运动方程是一个 n 阶线性常系数微分方程：,相应的传递函数为:,1.4 状态变量及状态空表达式的建立(二),实现的存在条件是，当 时，d=0;当m=n时，d=bm，不为零，此时，上式可写成如下形式：,这里只研究不存在零、极点对消的情况，所求得的状态空间描述中，状态变量数量最少，各矩阵的维数最小，构造硬件系统时所需的积分器个数最少，称为最小实现。,1.4.1 传递函数中没有零点时的实现,在这种情况下，系统的微分方程为：,1.4 状态变量及状态空表达式的建立(二),相应的系统传递函数为：,上式的实现，可以有多种结构，常用的简便形式可由相应的模拟结构图(下图)导出。这种由中间变量到输入端的负反馈，是一种常见的结构形式，也是一种最易求得的结构形式。,将图中每个积分器的输出取作状态变量，有时称为相变量，它是输出 的各阶导数。至于每个积分器的输入，显然就是各状态变量的导数。,依据上图，容易列出系统的状态方程：,输出方程为：,表示成矩阵形式，则为：,顺便指出，当 A 矩阵具有式上矩阵的形式时，称为友矩阵，友矩阵的特点是主对角线上方的元素均为1；最后一行的元素可取任意值；而其余元素均为零。,下面看个例题。,例1 设 求（A，B，C，D）解：选,则：,状态空间表达式为,1.4.2 传递函数中有零点时的实现,此时，系统的微分方程为：,相应地，系统传递函数为：,（26）,设待实现的系统传递函数为：,因为 n=m,上式可变换为,关键是选取合适的结构，使状态方程中不包含输入函数的导数项。先从三阶微分方程出发，找出其实现规律，然后推广到n阶系统。,令,则,对上式求拉氏反变换，可得：,据此可得系统模拟结构图，如下图所示。(先画出本页最上方传递函数的模拟结构图，再添加。),每个积分器的输出为一个状态变量，可得系统的状态空间表达式：,或表示为：,推广到 n 阶系统，式(26)的实现可以为：,（28）,可见状态方程与无零点时是相同的，只与传递函数的分母有关。,由于实现非唯一，下图示出上例的另一种实现方法，从输入、输出的关系看，二者是等效的。,可见输出方程与无零点时不同。根据这个结论，由传递函数中分子分母多项式的系数，可写出系统的状态空间表达式。,从图中可以看出，将输入函数的各阶导数作适当的等效移动，可用下图表示。,将综合点等效地移到前面，得到等效模拟结构图如下图所示。,由上图可求得其对应的传递函数为：,（29）,为求得 令式（29）与式（26）相等，通过对 多项式系数的比较得：,也可将式(30)写成式(31)的形式，以便记忆。,（31）,将上图的每个积分器输出选作状态变量，如图所示，得这种结构下的状态空间表达式：,扩展到n阶系统，其状态空间表达式为：,（33）,式中：,或记为：,值得注意的是，这两种方法所选择的状态变量是不同的。这点从它们的模拟结构图可以清楚地看到。,1.4.3 多输入一多输出系统微分方程的实现,以双输入一双输出的三阶系统为例，设系统的微分方程为：,（35）,同单输入单输出系统一样，式(35)系统的实现也是非唯一的。现采用模拟结构图的方法，按高阶导数项求解：,1.4.3 多输入一多输出系统微分方程的实现,对每一个方程积分：,故得模拟结构图，如下图所示：,取每个积分器的输出为一个状态变量，如上图所示。则式（35）的一种实现为：,或表示为：,（36）,1.5 状态矢量的线性变换(坐标变换),1.5.1 系统状态空间表达式的非唯一性,对于一个给定的定常系统，可以选取许多种状态变量，相应地有许多种状态空间表达式描述同一系统，也就是说系统可以有多种结构形式。所选取的状态矢量之间，实际上是一种矢量的线性变换(或称坐标变换)。,1.5 状态矢量的线性变换(坐标变换),1.5.1 系统状态空间表达式的非唯一性,设给定系统为：,（37）,我们总可以找到任意一个非奇异矩阵 T,将原状态矢量X 作线性变换，得到另一状态矢量Z,设变换关系为：,即,代入式（37），得到新的状态空间表达式：,（38）,由于T为任意非奇异阵，故状态空间表达式非唯一。通常称T为变换矩阵。,1.5.2 系统特征值的不变性及系统的不变量,1.系统特征值,系统,系统特征值就是系统矩阵 A 的特征值，也即特征方程：,（43）,的根。方阵A共有n个特征值；实际物理系统中，为实数方阵，故特征值或为实数，或为成对共轭复数；如 为实对称方阵，则其特征值都是实数。,2系统的不变量与特征值的不变性,同一系统，经非奇异变换后，得：,式(43)与式(44)形式虽然不同，但实际是相等的，即系统的非奇异变换，其特征值是不变的。可以证明如下：,将特征方程写成多项式形式 由于特征值全由特征多项式的系数 唯一确定，而特征值经非奇异变换是不变的，那么这些系数 也是不变的量。所以称特征多项式的系数为系统的不变量。,一个 维矢量：经过以 作为变换阵的变换，得到一个新的矢量 即,如果此 即矢量，经 线性变换后，方向不变，仅长度变化 倍，则称 为 的对应于 的特征矢量，此时有,3特征矢量,1.5.3 状态空间表达式变换为约旦标准型,根据系统矩阵A 求其特征值，可以直接写出系统的约旦标准型矩阵J无重根时,有重根时,而欲得到变换的控制矩阵 和输出矩阵CT，则必须求出变换矩阵T。下面根据A阵形式及有无重根的情况，分别介绍几种求T 的方法。,1.A阵为任意形式,(1)A阵的特征值无重根时,设 是A的 个互异特征根，求出A的特征矢量 则变换矩阵由A的特征矢量 构成，即,2A阵为标准型，即,(1)A的特征值无重根时，其变换是一个范德蒙德(Vandermonde)矩阵，为：,(2)A特征值有重根时，以有 的三重根为例：,(3)有共轭复根时，以四阶系统其中有一对共轭复根为例，即,此时,已知系统传递函数：,（55）,现将式(55)展开成部分分式。由于系统的特征根有两种情况：一是所有根均是互异的，一是有重根。,3系统的并联型实现,讨论此系统：,也有一个q重极点：,分析：,既有互异极点：,实现方法：整理得,系数 为待定系数，其中，采用留数定理计算：,（1）对于互异极点部分：,令,拉氏反变换可得：,（2）对于重极点部分：,令,则：,联立上两式得：,拉氏反变换可得：,联立(1)、(2)、(4)可得：,由(3)、(6)、(7)可得状态空间描述为：,xn,xq+1,x11,x12,x1q,y(t),u(t),+,+,+,+,+,-1,-q+1,-n,-1,-1,c11,c12,c1q,cq+1,cn,约当标准型状态结构图,1.6 从状态空间表达式求传递函数阵,1.6.1 传递函数(阵),1单输入一单输出系统,已知系统的状态空间表达式：,式中，为 维状态矢量；和 为输出和输入，它们都是标量；A为 方阵;为 列阵；c为 行阵；d为标量，一般为零。,（62）,1.6 从状态空间表达式求传递函数阵,对式(62)进行拉氏变换，并假定初始条件为零，则有：,（63）,故UX间的传递函数为：,（64）,它是一个 的列阵函数。,间的传递函数为：,它是一个标量。,2多输入一多输出系统,已知系统的状态空间表达式：,（66）,同前，对式(66)作拉氏变换并认为初始条件为零，得：,间的传递函数为：,（69）,它是一个mr矩阵函数，即,其中各元素 都是标量函数，它表征第 个输入对第 个输出的传递关系。当 时，意味着不同标号的输入与输出有相互关联，称为有耦合关系，这正是多变量系统的特点。,式(69)还可以表示为：,可以看出，的分母，就是系统矩阵A的特征多项式，的分子是一个多项式矩阵。,应当指出，同一系统，尽管其状态空间表达式可以作各种非奇异变换而不是唯一的，但它的传递函数阵是不变的。对于已知系统如式(66)，其传递函数阵为式(69)。当做坐标变换，即令 时，则该系统的状态空间表达式为：,那么对应上式的传递函数阵 应为：,即同一系统，其传递函数阵是唯一的。,（71）,1.6.2 子系统在各种连接时的传递函数阵,实际的控制系统，往往由多个子系统组合而成，或并联，或串联，或形成反馈连接。现仅以两个子系统作各种连接为例，推导其等效的传递函数阵。,设系统1为：,（72）,简记为：,设系统2为：,简记为：,1并联连接,所谓并联连接，是指各子系统在相同输入下，组合系统的输出是各子系统输出的代数和，结构简图如下图所示。,由式(72)和式(73)，并考虑 得系统的状态空间表达式：,从而系统的传递函数阵为：,故子系统并联时，系统传递函数阵等于子系统传递函数阵的代数和。,2串联连接,串联连接如下图所示：,其串联连接传递函数阵为：,即子系统串联时，系统传递函数阵等于子系统传递函数阵之积。但应注意，传递函数阵相乘，先后次序不能颠倒。,3具有输出反馈的系统,如下图所示，由图可得：,即,从而系统的传递函数阵为：,可推导出：,同理也可求得：,本章到此结束！谢谢您的观看！请提宝贵意见！,