计算机控制系统电子教案单元设计 (9).doc





《计算机控制系统电子教案单元设计 (9).doc》由会员分享,可在线阅读,更多相关《计算机控制系统电子教案单元设计 (9).doc(70页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第9章 线性离散系统的离散状态空间分析、设计法9.1 线性离散系统的离散状态空间分析法在连续控制系统中,状态空间分析法是分析、研究系统的有力工具,它解决了频率特性解决不了的问题,如多变量问题、时变问题等。计算机的广泛普及和应用为状态空间分析法提供了有力的手段。对于离散系统同样可以用离散状态空间分析法来研究和分析。离散状态空间分析法有以下的优点: 离散状态空间表达式适宜于计算机求解。 离散状态空间分析法对单变量和多变量系统允许用统一的表示法。 离散状态空间分析法能够应用于非线性系统和时变系统。9.1.1 线性离散系统的离散状态空间表达式在线性连续系统中,是用控制变量,状态变量和输出变量来表征系统
2、的动态特性的,如图9.1所示。状态变量是表征系统本身特性的变量。系统的状态变量是可以有多种选择方案,但是当系统确定时,状态变量的个数就确定 图9.1 线性连续系统的变量关系了,而且是最少的,它就是系统的阶数。 状态变量可以表示成列向量 (9.1)控制变量可以表示成列向量 (9.2)输出变量可以表示成列向量 (9.3)线性连续系统的状态空间表达式为 (9.4) (9.5)A,B,C,D是定常的系数矩阵。式(9.4)称为状态方程,式(9.5)称为输出方程。与线性连续系统类似,线性离散系统的离散状态空间表达式可以表示为 (9.6) (9.7)式(9.6)称为状态方程,式(9.7)称为输出方程。 F是
3、维矩阵,称为状态矩阵或系统矩阵。 G是维矩阵,称为输入矩阵或驱动矩阵。 C是维矩阵,称为输出矩阵。 D是维矩阵,称为直传矩阵或传输矩阵。图 9.2 线性离散系统的状态变量图 线性离散系统的状态变量图如图9.2所示。9.1.2 由差分方程导出离散状态空间表达式 对于单输入、单输出的线性离散系统,可以用阶差分方程来描述 (9.8)或表示成 (9.9) 为了得到状态空间表达式,应适当选择状态变量,将高阶的差分方程化为一阶差分方程组,然后表示成向量的形式,便可以得到离散状态空间表达式了。1. ,即控制变量不包含高于一阶的差分(9.10)选择状态变量 (9.11)由(9.11)式,可得(9.12)写成矩
4、阵形式 (9.13) (9.14)将式(9.13),式(9.14)表示成向量形式 (9.15)式中状态矩阵为 (9.16)输入矩阵为 (9.17)输出矩阵为 (9.18)直传矩阵为 (9.19)例 9.1 线性离散系统的差分方程为试导出离散状态空间表达式。 解 由差分方程知:(输出向量维数)。可知 2. ,即控制变量包含高于一阶的差分令,则差分方程为 (9.20)选择状态变量(9.21)作为一组状态变量,其中是待定常数。下面,我们来确定它们的值,使状态方程中不出现的高阶差分。 对式(9.21)中最后一个方程的两边前移一拍,得 (9.22)分别用乘式(9.21)中的第1,第2,第个方程,得(9.
5、23)将方程(9.22)和方程(9.23)中的个方程一齐加起来,得 (9.24)将方程(9.20) 代入方程(9.24),并加以整理,则得 (9.25)由上式可以看出,若令上式右边除最后一项外,其余各项系数为零,则状态方程中不会出现的高阶差分。据此可以解出 (9.26)这样一来,式(9.21)中的个状态变量的个待定常数就唯一地确定出来了。再将方程(9.25)中的系数简记为,连同(9.21)中后面个方程一起,便得到系统的离散状态空间表达式 (9.27)由式(9.21)的第一个方程,得系统输出方程为 (9.28)将式(9.27)和(9.28)写成矩阵形式 (9.29)其中 在这个表达式中,不论输入
6、函数是否含有高阶方程,系数矩阵都是一样的,而的高阶差分只会改变控制矩阵的各元。系统的初始条件可由下式决定(9.30) 例 9.2 设线性离散系统的差分方程为 试写出离散状态空间表达式。 解 设状态变量则系统的离散状态空间表达式为 初始条件可由下式决定此题也可以直接利用公式求出F,G,C,从而得到离散状态空间表达式。 线性离散系统的阶数为, 9.1.3 线性离散系统状态方程的求解 线性离散系统的离散状态方程是由高阶差分方程化为一阶方程组得到的,所以求解差分方程的方法可以适用于求解离散状态方程。通常离散状态方程的求解方法有迭代法和Z变换法。1. 迭代法设线性离散系统的离散状态空间表达式为 (9.3
7、1)初始值。以代入式(9.31)可得 (9.32)有了离散状态方程的解(9.32)式,便可得到输出方程 (9.33)令式(9.32)中的,称为线性离散系统的状态转移矩阵。状态转移矩阵具有如下的性质 (9.34)式中是单位矩阵。 可以看作是式(9.34)矩阵差分方程的唯一解。它描述了线性离散系统由的初始状态向任意时刻的状态转移的特性。因此称为线性离散系统的状态转移矩阵。 线性离散系统的解式(9.32)还可以用状态转移矩阵来表示,即 (9.35) 例 9.3 用迭代法求解线性离散系统的状态方程 解 令及初始条件代入离散状态空间表达式,可以得到 将以上计算数据列表9.1。表9.1迭代计算数据表012
8、311-0.3-0.69-10.30.691.327-10.30.691.327 用迭代法求解离散状态方程只能得到有限项时间序列,得不到状态变量和输出变量的数学解析式。2. Z变换法设线性离散系统的状态空间表达式为 (9.36)对(9.36)式作Z变换,可得 (9.37)对式(9.37)作Z反变换,可得 (9.38)对比式(9.32),下式成立 (9.39)例 9.4 用Z变换法求解线性离散系统的状态方程 解 取Z的反变换即可得到方程的解 用Z变换法求解离散状态方程,可以得到状态变量和输出变量的数学解析式。9.1.4 线性离散系统的Z传递矩阵设线性离散系统的状态空间表达式为 (9.40)式中
9、是维状态向量;是维状态向量;是维状态向量。对式(9.40)作Z变换,可得初始条件为零,即时,则有 因而 (9.41)称为线性离散系统的Z传递矩阵(维矩阵)。它反映了在初始静止的条件下,输出量的Z变换与输入量的Z变换之间的关系。 对于单输入、单输出系统,是维矩阵,即为脉冲传递函数(Z传递函数)。 例 9.5 设线性离散系统的状态空间表达式为初始条件为零。试求线性离散系统Z传递矩阵,并求出单位阶跃输入时的输出响应。 解 单位阶跃输入时对上式作Z反变换可得 9.1.5 线性离散系统的Z特征方程在线性连续系统中,用特征方程来表征系统的动态特性,同样在线性离散系统中引进Z特征方程的概念来描述一个线性离散
10、系统的动态特性。 设线性离散系统的状态方程为 (9.42)对式(9.42)作Z变换,可得仿照线性连续系统,令矩阵的行列式 (9.43)称式(9.43)称为线性离散系统的Z特征方程。 例 9.6设线性离散系统的状态矩阵试求线性离散系统的Z特征方程。 解 (9.44)式(9.44)即为线性离散系统的Z特征方程或称为矩阵的Z特征方程。Z特征方程的根也称为矩阵的特征值,就是线性离散系统的极点。 对于一个阶的系统,仅有个特征值。因为特征方程表征系统的动态特性,因此尽管一个系统的状态变量的选择不是唯一的,但是系统的Z特征方程是不变的。 例 9.7 设线性离散系统的Z传递函数为试导出离散状态方程,并求出Z特
11、征方程。解 用直接程序法求系统的离散状态方程(参阅第8章内容,直接程序法) (9.45)上式中,令 (9.46)由(9.46)式得出 (9.47) (9.48)(9.49) 图 9.3 例9.7状态变量框图根据式(9.47)和式(9.49)画出状态变量框图如图9.3所示。 由图9.3 可知 (9.50) (9.51)由上式可知则,Z特征方程为特征值 用并联程序法求系统离散状态方程令 , (9.52) (9.53)由(9.52)式和(9.53)式,画出系统状态变量框图如图9.4所示。 图 9.4 并联程序框图 由图9.4可以写出系统离散状态方程 (9.54) (9.55)由式(9.54)可知则,
12、Z特征方程为特征值 由本例可以清楚看到一个线性离散系统,用不同的方法可以得到不同的离散状态方程,如式(9.50)、(9.51)以及(9.54)、(9.55)中的状态矩阵是不同的。但是,它们的特征方程是相同的,因而它们的特征值也相同。而且特征值的个数就等于离散系统的阶数。9.1.6 计算机控制系统的离散状态空间表达式 一个计算机控制系统,通常由数字计算机和连续环节组成。计算机控制系统的离散状态空间表达式,可以通过连续环节的状态空间表达式离散化得到,也可以用Z变换得到。下面介绍用Z变换求得计算机控制系统的离散状态空间表达式。 图 9.5 计算机控制系统 例 9.8 已知计算机控制系统如图9.5所示
13、,初始静止。试求系统的离散状态空间表达式。 解 广义被控对象的传递函数广义对象的Z传递函数图 9.6 图9.5的等效方框图 (9.56)由式(9.56)可以建立系统的方框图,如图9.6所示。 选择状态变量 (9.57)由式(9.57)可得(9.58)由式(9.58)可得离散状态方程,由图9.6可得离散状态空间表达式。(9.59)当时,系统特征方程为 图9.7 计算机控制系统方框图当用连续部分的状态空间表达式离散化时,也即要求导出连续部分的离散状态方程。假设计算机控制系统的方框图如图9.7所示。由图9.7可见,计算机控制系统由连续部分和离散部分组成。是数字控制器,为离散部分;是保持器传递函数,它
14、与被控对象合在一起为,称为广义被控对象。绝大多数系统中,为零阶保持器,它的输出为阶 图 9.8 零阶保持器的输出梯信号,如图9.8所示。连续对象的动态特性可以用状态空间表达式表示 (9.60)式(9.60)的解为 (9.61) 当为阶梯信号时 (9.62) 初始条件为 积分上限为 积分下限为 ,且 由式(9.61)可得 (9.63)在积分区间内,输入是常数,而且积分对所有的都成立,作变量置换,则有 (9.64)式(9.64)是与有关的常数矩阵。将其代入(9.63)式,可得 (9.65)式(9.65)即为所求的离散状态方程,写成标准形式为 (9.66)式中 图9.9 计算机控制系统结构图 例 9
15、.9已知计算机控制系统如图9.9所示,试求计算机控制系统的离散状态空间表达式。解 由图9.9 可得对象的状态空间表达式可以写成下列形式由于系统采用零阶保持器,是阶梯信号,由此可求得对象的离散状态空间表达式。 由对象的状态空间表达式,可知则对象的离散状态空间表达式的状态矩阵对象的输入矩阵所以,离散状态方程为以代入上式 经整理得 系统特征方程 对照例9.8,对于同一个计算机控制系统,用不同的方法都能得出相同的Z特征方程。9.1.7 用离散状态空间法分析系统的稳定性 设线性离散系统的状态方程为 (9.67)用迭代法求(9.67)的解,得 (9.68)分析系统的稳定性只考虑状态方程的齐次解,即的情况,
16、则 (9.69) 设是矩阵,具有两两相异的特征值。根据西尔维斯特展开定理,可以展开成级数 (9.70)式中称为的相容矩阵或要素矩阵,且 (9.71)矩阵跟及其个特征值有关,但是是确定的,与无关。于是状态方程的解为 (9.72) 显然只有,当时,才收敛,系统才是稳定的。 由此可以得到线性离散系统的稳定判据:线性离散系统稳定的充分必要条件是系统的Z特征方程的所有特征根。特征值在Z平面上的分布与稳定性关系如图9.10所示。图 9.10 线性离散系统的稳定区 例 9.10 设线性离散系统如图9.5所示,。试判断采样周期时系统的稳定性。 解 由例9.8可知,系统中,当时,线性离散系统的状态空间表达式为系
17、统的状态矩阵为 当T=1sec时Z特征方程为,可知所以系统稳定。 当T=3sec时Z特征方程为,可知所以系统稳定。 当T=5sec时Z特征方程为,可知所以系统不稳定。从例9.9可以看到,可用Z特征根在Z平面上的分布,来判断线性系统的稳定性。同时也可以看到,在线性离散系统中,采样周期T是系统的一个重要参数:当T比较小时,系统稳定;当T增大时,特征根的模加大,在Z平面内向单位圆靠近;当T大到一定值时,特征根的模大于1,即在单位圆外,系统变得不稳定了。9.2 线性系统离散状态空间设计法离散状态空间设计法是利用离散状态空间表达式,根据性能指标要求,设计出满足要求的计算机控制系统。离散状态空间设计法的主
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机控制系统电子教案单元设计 9 计算机控制 系统 电子 教案 单元 设计

限制150内