第四章-控制系统数字仿真-课件.ppt

1、第四章 控制系统数字仿真 4.1第4章控制系统数字仿真 第四章 控制系统数字仿真 4.2第四章 控制系统数字仿真 4.3连续系统数值积分方法 连续系统的数学模型，一般都能以微分方程的形式给出，所以，连续系统数学仿真算法问题通常可归结为如何用计算机来求解微分方程的问题。数值积分法是解决该问题的重要方法之一，例如已知00d(,)d()()yf y tty tty t(4.1)求()y t第四章 控制系统数字仿真 4.4连续系统数值积分方法解：对式(4.1)两边积分，则00()()(,)dtty ty tf y tt当 1mtt1mtt0mtt11()()(,)dmmtmmty ty tf y tt

2、时 令 1(,)dmmtmtQf y tt则 1()()mmmy ty tQ(4.2)数值积分法是解决在已知初值的情况下，对(,)f y t进行近似积分，对()y t进行数值求解的方法，在数学上称为微分方程初值问题的数值方法。第四章 控制系统数字仿真 4.5连续系统数值积分方法一.欧拉法mQ表示函数(,)f y t在 mt和 1mt相邻两次采样时刻之间的积分。若将此定积分中的(,)f y t近似看成常数，即(,)(,)mmf y ty yt则(,)mmmQf yth由式(4.2)得1()()(,)mmmmy ty tf yth即1(,)mmmmyyf yth(4.3)第四章 控制系统数字仿真

3、4.6 连续系统数值积分方法式中，1mmhtt为计算步距。0m，1，2，(1)n 式为(4.3)著名的欧拉积分公式。欧拉公式计算简单，但精度较低，其原因是将 (,)f y t看成常数(,)mmf yt，从而用矩形面积代替准确的曲面面积 mQ，形成了较大的误差面积。为了提高计算精度产生了梯形法。第四章 控制系统数字仿真 4.7连续系统数值积分方法二.梯形法mQmQmQmQ梯形法是用梯形的面积近似代替定积分 ，即 式中，；，所以 即mQ1()2mmmhQff(,)mmmff yt111(,)mmmff yt11()()()2mmmmhy ty tff11()2mmmmhyyff(4.4)第四章 控

4、制系统数字仿真 4.8连续系统数值积分方法(0)1(0)11(,)()2mmmmmmmmyyhf ythyyff从式(4.4)可以看出，用梯形法进行数值积分，会出现一个问题：在计算 1my时，先要用1my去计算式(4.4)右端的 1mf，而此时 还未求出，1my显然，这是不可能实现的。所以，一般采用欧拉公式先预报一个(0)1my然后将预报 的代入式(4.4)进行校正，求出 ，即(0)1my1my(4.5)式中，(,)mmmff yt(0)(0)111(,)mmmff yt一般称式(4.5)为预报校正公式。显然，梯形法比欧拉法精度要高，同时，每前进一个步距，计算工作量也比欧拉法约多一倍。第四章

5、控制系统数字仿真 4.9连续系统数值积分方法 三.对于微分方程式(4.1)，若在其初值附近展开成泰勒级数，并只取前三项，则有2210200d1dd2dyyyyhhtttttt200001(,)()2ffyf y t hfhttty(4.6)设式(4.1)的解可以写成如下形式101 122()yya ka kh100(,)kf y t202101(,)kf yb k h tb h(4.7)第四章 控制系统数字仿真 4.10连续系统数值积分方法2k00200121(,)t tyyffkf y tbhb k hty将 、代入式(4.7)2k1k10100200121(,)(,)ffyya hf y

6、taf y tbhb k hhty(4.8)比较式(4.6)与式(4.8)，得122 12211212aaa ba b第四章 控制系统数字仿真 4.11连续系统数值积分方法 三个方程中有四个未知数，因而解不是唯一的。所以，限定三个方程中有四个未知数，因而解不是唯一的。所以，限定 ，则求，则求出出 12aa1212aa121bb。代入式(4.7)得1012()2hyykk100(,)kf y t2010(,)kf yk h th写成一般形式112()2mmhyykk1(,)mmkf yt21(,)mmkf yk h th(4.9)第四章 控制系统数字仿真 4.12连续系统数值积分方法 ()mmy

7、y t1mtt(1)my t(1)my t 1my为用式(4.9)求得的近似解，则它们之差 ，称为此时计算的截断误差。111()mmmy ty式(4.6)只取到泰勒级数展开式中 的二阶导数项，略去了三阶以上高阶导数项。为纪念提出该方法的德国数学家C.Runge和M.W.Kutta，称这种计算方法为二阶龙格-库塔法。其截断误差正比于步长 的三次方，同理若在式(4.6)的计算中，取到的 三阶或四阶导数项，则有相应的三阶或四阶龙格-库塔法，相应的截断误差也应正比于 或 。yhy4h5h第四章 控制系统数字仿真 4.13连续系统数值积分方法 一般在计算精度要求较高的情况下，多使用四阶龙格-库塔法。其计

8、算公式为11234(22)6mmhyykkkk1(,)mmkf yt21(,)22mmhhkf yk t32(,)22mmhhkf yk t43(,)mmkf yhk th(4.10)第四章 控制系统数字仿真 4.14连续系统数值积分方法 将欧拉公式与 的泰勒级数展开式比较可知，欧拉公式即泰勒级数展开式只取前两项。而将预报校正公式与二阶龙格库塔公式比较，发现它们完全是一样的。这样，通过对龙格库塔法的介绍，可将前面讲的几种数值积分方法统一起来，都可看成是 在初值附近展开成泰勒级数所产生的，只是其泰勒级数所取项数的多少不同。欧拉法只取前两项，梯形法与二阶龙格库塔法取了前三项，四阶龙格库塔法取了前五

9、项。从理论上说，取的项数越多，计算精度越高，但相应的计算公式更复杂，计算工作量也更大。y()y t第四章 控制系统数字仿真 4.15连续系统数值积分方法 以上介绍的几种数值积分公式，有一个共同的特点，在计算 时只用到 ，而不直接用 ，各项的值，即在本次计算中，仅仅用到前一次的计算结果，而不需要利用更前面各步的结果。这类计算方法称为单步法。单步法运算有下列优点：1mymy1my2my(1)需要存储的数据量少，占用的存储空间少。(2)只需要知道初值，就可启动递推公式进行运算，具有这种能力的计算方法称为自启动的计算方法。(3)容易实现变步长运算。与单步法相对应的还有一类数值积分方法，在它的数值积分公

10、式中，本次计算不仅利用前一次的计算结果，还必须利用更前面各次结果，此类方法称为多步法。如四阶阿达姆斯(Adams)积分公式第四章 控制系统数字仿真 4.16连续系统数值积分方法 1123(5559379)24mmmmmmhyyffff(4.11)1mymf1mf2mf3mf1k2k3k4k(,)f y t第四章 控制系统数字仿真 4.17连续系统数值积分方法 四.前面介绍的数值积分公式，都是以求解单个典型微分方程 进行介绍的。而工程实际系统中大量的仿真对象是以一阶微分方程组或矩阵微分方程的形式给出，如 d/d(,)ytf y t1112221212(,)(,)(,)nnnnnxf x xx t

11、xfx xx txfx xx t&L&LM&L或 (,)Ft&XX第四章 控制系统数字仿真 4.18连续系统数值积分方法 12,nx xxLTX11221212(,)(,)(,)(,)nnnnf x xx tfx xx tFtfx xx tLLMLX在这种情况下，各数值积分公式显然应采用相应的矩阵形式 第四章 控制系统数字仿真 4.19连续系统数值积分方法(1)欧拉公式1(,)mmmmhFtXXX(4.12)(2)梯形公式11()2mmmmhFFXX(4.13)(3)二阶龙格-库塔公式1111()2(,)(,)mmmmmmmmmmmhFFFF XtFFhF tXXX(4.14)第四章 控制系统

12、数字仿真 4.20连续系统数值积分方法11234(22)6mmhXXKKKK(4.15)112111(,)mnkkFtkMKX第四章 控制系统数字仿真 4.21连续系统数值积分方法 1222212(,)22mmnkkhhFtkMKXK1323323(,)22mmnkkhhF XtkMKK第四章 控制系统数字仿真 4.22连续系统数值积分方法1424434(,)mmnkkFhthkMKXK对一个 维向量 ，每前进一个步距 ，至少要求 个 之值。nXh4nijk对常见的线性定常系统 ，四阶龙格库塔法的四个可表示为&XAXBu第四章 控制系统数字仿真 4.23连续系统数值积分方法 ()mmu t1K

13、AXB21()22mmhhu tKA XKB32()22mmhhu tKA XKB43()mmAhu thKXKB第四章 控制系统数字仿真 4.24连续系统数值积分方法 对 维向量 ，取 是有四个分量的一维矢量，其中 再取 为一零向量，其中 ，。nXjh10h 32hh 4hh0K00ik1,2,inL则四个求jK的公式可合并为一个公式1()jmjjmjhu thKA XKB1,2,3,4j(4.16)在以后介绍的程序中，四阶龙格库塔法的计算将建立在上式的基础上。第四章 控制系统数字仿真 4.25连续系统数值积分方法 五.数值积分方法的计算稳定性 这里所说的数值积分方法的计算稳定性，在数学上是

14、指微分方程初值问题算法的数值稳定性或计算稳定性。由于在系统仿真运算阶段大量求解微分方程问题时，实际上是用不定数值积分方法求解微分方程的问题。如欧拉法就是在已知 时，用 完成对微分方程 的求解。习惯上将这些方法称为数值方法。0nyy1nnnyyhfd(,)dyf x yx从稳定性理论中，可以知道如何从系统的微分方程或传递函数去判断该系统的稳定性。那么，对于一个稳定的连续系统，当用某数值积分方法进行仿真计算时，是否仍然稳定呢？先看下面的例子。第四章 控制系统数字仿真 4.26连续系统数值积分方法 【例4.1】已知微分方程d30d1(0)3yyty 01.5t 其精确解为 301()e3ty t取

15、，用欧拉法和四阶龙格-库塔法计算 时的 ：0.1h 1.5t()y t第四章 控制系统数字仿真 4.27连续系统数值积分方法 欧拉法 4(1.5)1.09225 10y 四阶龙格-库塔法 1(1.5)3.95730 10y精确解 21(1.5)9.54173 10y第四章 控制系统数字仿真 4.28连续系统数值积分方法 显然，此时数值积分法计算的结果是错误的。为什么会出现这种情况呢？这是因为数值积分方法只是一种近似方法，它在反复的递推运算中将引入误差。若误差的积累越来越大，将使计算出现不稳定，从而得出错误的结果。所以，原系统的稳定性与数值积分法计算的稳定性是不同的两个概念。前者用原系统的微分方

16、程、传递函数来讨论，后者用逼近微分方程的差分方程来讨论。由于选用的数值积分法不同，即使对同一系统，差分方程也各不相同，稳定性也各不一样。如何分析数值积分法的数值计算稳定性呢？一般来讲，对高阶微分方程的数值计算稳定性作全面分析是比较困难的，通常用简单的一阶微分方程来考查其相应差分方程的计算稳定性。第四章 控制系统数字仿真 4.29连续系统数值积分方法 微分方程0dd(0)yytyy(4.17)称测试方程。据稳定性理论，当其特征方程的根在平面的左半平面，即根的实部 时，则原系统稳定。此时相应的数值积分法的计算稳定性如何呢？Re0第四章 控制系统数字仿真 4.30连续系统数值积分方法 1.欧拉法的计

17、算稳定性欧拉法的计算稳定性对式(4.17)按欧拉公式计算有1(1)mmmmyyh yhy0,1,(1)mnL(4.18)用式(4.18)计算的 ，并不是 时的真值，而只是真值包含了各种误差的近似值。随着递推次数的增加，此误差是否会不断扩大，使 完全不能表示此时的真值呢？这决定了此差分方程计算的稳定性。为简化问题的讨论，假设只在 时刻有误差引入，这样 的误差 仅由 的误差引起，所以1my(1)tmh1mytmh1my1mmy11()(1)()mmmmmmmmyyhyhy(4.19)第四章 控制系统数字仿真 4.31连续系统数值积分方法 式(4.19)减式(4.18)得误差方程1(1)mmh 同理

18、 221(1)(1)mmmhh (1)nm nmh 当 ，则 ，表明此算法若在计算中的某步引入了误差，随着计算步数的增加，这个误差将逐渐扩大，以致差分方程的解完全失真。11h1mm第四章 控制系统数字仿真 4.32连续系统数值积分方法 反之，当 ，则 随着计算步数的增加，误差 逐渐减小并趋于零或有界，在这种情况下，称此差分方程的算法是计算稳定的。显然，合理的选择步长使其满足 ，是保证欧拉法计算稳定的重要条件。11h1mm1m11h第四章 控制系统数字仿真 4.33连续系统数值积分方法 2.梯形法的计算稳定性梯形法的计算稳定性对测试方程按梯形法计算公式，则111()2mmmmyyyy1111()

19、()2mmmmmmmmhyyyy(4.20)111(1)/(1)22mmhh 所以(4.21)第四章 控制系统数字仿真 4.34连续系统数值积分方法 若原系统稳定，则根据稳定性理论，应是小于零的实数。对式(4.21)，取任何正数都可以保证下式成立：h1121112hh所以 1mm这就表明，某步计算引入的误差，将随着计算步数的增加而减小。也就是说，梯形法的计算在任何步长下都是稳定的，是一种绝对稳定的计算方法。第四章 控制系统数字仿真 4.35连续系统数值积分方法 六.数值方法的计算精度、速度、稳定性与步距的关系从上面的分析可以看出，一般数值积分方法的稳定性与所选用的计算步距有关。各种数值积分方法

20、的差分方程是对原微分方程的近似逼近，存在明显的截断误差。且计算机的字长有限，计算只能限制在有限位，将引入舍入误差。这些误差都和计算步距密切相关，所以计算步距是影响计算精度、速度和稳定性的重要因素。若步距取得较大，截断误差就会相应地增大；反之，若步距取得较小，截断误差就会减小，但在给定时间范围内，计算次数必然增加，使舍入误差积累，相应地增加。第四章 控制系统数字仿真 4.36连续系统数值积分方法 这两种误差与步距的关系可用图4.1表示。从图中可知，两种误差对步距的要求是矛盾的，但两者之和有一个最小值，步距最好能选在最小值。然而，实际要做到这一点是很困难的。一般只能根据经验确定一个 附近的合理步长

21、区，如可将 限制在系统的最小时间常数数量级上。0hh图4.1 误差与步距关系第四章 控制系统数字仿真 4.37控制系统的结构及其描述 一.控制系统的典型结构形式 1.单输入-单输出开环控制结构控制系统结构如图4.2所示。其中控制器、控制对象等环节可用第3章所述任何一种数学模型描述,为参考输入量,为控制量,为输出量。()r t()u t()y t图4.2 单输入-单输出开环控制结构图第四章 控制系统数字仿真 4.38控制系统的结构及其描述 2.单输入-单输出前馈控制结构控制系统结构如图4.3所示。这种结构形式本身仍为开环控制形式，但为补偿输入引起的误差，在已知误差变化规律的情况下，加入补偿环节，

22、对误差作提前修正。补偿环节同样可用任一种数学模型描述。图4.3 单输入-单输出前馈控制结构图第四章 控制系统数字仿真 4.39控制系统的结构及其描述 3.单输入单输出闭环控制结构这是控制系统中应用最广泛的控制结构形式，大多数控制系统为保证其控制精度、性能指标要求，都采用这种闭环负反馈形式，如图4.4所示。该系统控制原理主要是根据偏差(或误差)确定控制量，使输出量按希望的精度变化，减小或消除偏差，即所谓的偏差控制系统。图中 为偏差量(或误差量)，为反馈量。反馈环节可用任一种数学模型描述。当反馈为单位反馈时，有()e t()b t()()b ty t则偏差()()()()()e tr tb tr

23、ty t第四章 控制系统数字仿真 4.40控制系统的结构及其描述 图4.4 单输入-单输出闭环控制结构图称作误差，而此时的闭环系统称作单位反馈系统。第四章 控制系统数字仿真 4.41控制系统的结构及其描述 凡是单输入-单输出控制系统结构均能方便地表示为图论中的拓扑结构形式，即使系统环节再增多，反馈和前馈联系再复杂，表述也很清晰，如图4.5所示。图中各环节用序号表示，每个环节都有自己的输入和输出变量 和 (1,2,.，6)，对应其相互关系的数学描述，可以是第3章所述形式的任一种。iuiy图4.5 控制系统拓扑结构图第四章 控制系统数字仿真 4.42控制系统的结构及其描述 4.多输入-多输出控制系

24、统结构若干个单输入-单输出控制系统结构形式通过一定方式组合在一起，构成了多输入-多输出控制系统结构形式，用来描述较复杂的多变量控制系统。根据不同的环节组合形式，也有不同的控制结构。如变量和环节有单向耦合作用关系的结构形式见图4.6(a)、交叉反馈耦合作用关系的结构形式见图4.6(b)等。第四章 控制系统数字仿真 4.43控制系统的结构及其描述 (a)单向耦合结构 (b)交叉反馈耦合结构图4.6 多输入-多输出控制结构图第四章 控制系统数字仿真 4.44控制系统的结构及其描述 二.控制系统的典型环节描述任何一个复杂的控制系统，都是由一些简单的不同类型的典型环节采用串联、并联、反馈等形式组合而成的

25、。若对常见的一些典型环节能准确地加以定量描述，则复杂系统的描述也只是复杂在各部分的相互连接关系上。采用前述系统关系描述，把这种连接关系用相应的数学关系表达出来，就可以得到在计算机上能方便运行的“仿真模型”。经典控制理论中常见的典型环节如下所述：第四章 控制系统数字仿真 4.45控制系统的结构及其描述 ()iiiiyG sKu()()()iY sG sKU s()1iiiiiyKG suTs()()()1Y sKG sU sTs(1)比例环节 (2)惯性环节 第四章 控制系统数字仿真 4.46控制系统的结构及其描述 ()y t()(1)()()1Y sKsG sU sTs()iiiiiyKG s

26、uTs()()()Y sKG sU sTs(3)惯性比例环节(4)积分环节 第四章 控制系统数字仿真 4.47控制系统的结构及其描述 (1)()iiiiiiyKsG suTs()(1)()()Y sKsG sU sTs22222()221iiniiinniiyKKG sussT sTs222()()()2nnnY sG sU sss(5)比例积分环节(6)二阶振荡环节 式中，为阻尼比，01时，系统为欠阻尼二阶振荡情况；为无阻尼自然振荡频率。n第四章 控制系统数字仿真 4.48控制系统的结构及其描述 ()eisiiiyG su()()e()sY sG sU s0101()mimmininnyb

27、sbsbG sua sasaL LL L110110()()()mmnnnb sbsb sbY sG sU sa sasa saLL(7)纯滞后环节 式中，(8)高阶线性环节为纯滞后时间。第四章 控制系统数字仿真 4.49控制系统的结构及其描述 式中，,为分母多项式系数，亦称特征多项式系数；,为分子多项式系数；它们均为反映高阶线性环节动态响应性能的重要参数。0a1ana0b1bmb以上典型环节，无论其阶次高低，均可应用MATLAB语言方便地求出输入量作用下输出量的响应性能。和 之间的数学表达关系可以由传递函数 方便地转化为其他各种形式，如状态方程、零极点增益、部分分式等。YU()G s第四章

28、控制系统数字仿真 4.50控制系统的结构及其描述 三.控制系统的连接矩阵 对图4.7所表示的线性系统，各环节均为线性的，在各自的输入量 作用下，给出各自的输出量 ，这种作用关系是通过各环节的数学描述体现出来的，其数学关系可以是(4.2.2)节所述的任何一种。但各环节之间存在相互作用，不是孤立的，只要与其他环节有连接关系，就要受到相应 变化的影响。因此，要完整地将系统描述出来，还应该分析各环节输出 对其他环节有无输入作用，才能完整地进行仿真分析。iuiyiuiyiy第四章 控制系统数字仿真 4.51 控制系统的结构及其描述图4.7 控制系统拓扑结构第四章 控制系统数字仿真 4.52控制系统的结构

29、及其描述 根据图4.7中 、的连接关系，可逐个写出每个环节输入 受哪些环节输出 的影响，现列写如下：iuiyiuiy1144522122324353uryyuryyuyuyuy(4.22)第四章 控制系统数字仿真 4.53控制系统的结构及其描述 引入向量 T12,r rR0UW RWYT12345,y yyyyY则式(4.22)可表示为11422213324455100001011000000100000001000000100uyuyruyruyuy 第四章 控制系统数字仿真 4.54控制系统的结构及其描述 即 0UW RWY式中，01001000000W420001100001000001

30、0000100W第四章 控制系统数字仿真 4.55控制系统的结构及其描述 式中，W为连接矩阵，nn型，矩阵中元素清楚地表示出各环节之间的连接关系，W0为输入连接矩阵，nm型(当参考输入为m维时)，阵中元素表示环节与参考输入之间的连接关系。仔细研究连接矩阵W，可从其元素值直接看出各环节之间的连接情况，即wij0，环节j不与环节i相连；wij 0，环节j与环节i有连接关系；wij 0，环节j与环节i直接相连(wij=1)或通过比例系数相连(wij为任意正实数)；wij 0，环节j与环节i直接负反馈相连(wij=-1或通过比例系数负反馈相连(wij为任意负实数)；wii 0，环节i单位自反馈(wii

31、=l或wii=-1)或通过比例系数自反馈(wii为任意实数)。第四章 控制系统数字仿真 4.56面向结构图的数字仿真 一.典型闭环系统的数字仿真 1.典型闭环系统结构形式控制系统最常见的典型闭环系统结构如图4.8所示。其中0101()()()mmmnnb sbsbY sG sU sa sasaLL11101110()()()mmmmnnnnb sbsb sbY sG sU sa sasa saLL第四章 控制系统数字仿真 4.57面向结构图的数字仿真 图4.8 典型闭环系统结构图第四章 控制系统数字仿真 4.58面向结构图的数字仿真 2.系统仿真模型与求解思路所谓仿真模型是指经一定方式把数学模

32、型转化为便于在计算机上运行的表达形式。这种表达形式往往是一些适合于具体编程实现的数学关系式。对图4.8所示系统的开环传递函数，可按照能控标准型写出其开环状态方程UYAXBCXX X(4.23)式中,01010000101naaLLLLOLLA0001 MB0100mbbbLLC第四章 控制系统数字仿真 4.59面向结构图的数字仿真 其中，矩阵中，(1,2,，0,1,，,)为式(3.4)首一化后分母、分子各系数，即 ，,且 。ACjaibj ni mjjnaaaiimbbb1mn由图又知，控制量 ，代入式(4.23)得UrvY()rvY&XAXB再由，Y CX，则()vrr&XAB C XBAX

33、B(4.24)即得系统闭环状态方程，其中vA A B C为系统闭环系数矩阵，而输入矩阵B和输出矩阵C不变。这就是图4.8系统的仿真模 第四章 控制系统数字仿真 4.60面向结构图的数字仿真 仿真模型一旦确立，就可以考虑求解与编程实现。观察式(4.24)可知，该式其实为一个一阶微分方程组的矩阵表达形式，而数值积分法最适宜解一阶微分方程，当采用四阶龙格库塔法求解此闭环状态方程时，其步骤如下：r&XAXB(,)f tXrXAB对 知此为对应 个状态变量的一阶导数 的n维向量表达式。其中，r为随时间变化的已知输入函数。于是，当求 解过程进行到 时刻，欲求 时刻各量，需先求n&XT12,nx xxLXk

34、tt1kt第四章 控制系统数字仿真 4.61面向结构图的数字仿真 1(,)()kkkkKf tr tXAXB211(,)()()2222kkkkhhhhKf tKKr tXA XB312(,)()()2222kkkkhhhhKf tKKr tXA XB433(,)()()kkkkKf thhKhKr thXA XB以上各系数 、可认为是对应维状态变量在时刻的四组斜率，每组为维斜率向量，即每组个龙格库塔系数。1K2K3K4K第四章 控制系统数字仿真 4.62面向结构图的数字仿真 最后，再由 11234(22)6kkhKKKKXX求得 时刻状态 ，立即可得输出相应时刻值1kt1kX11kky CX

35、按以上算式，取k0，1，2，N不断递推，即求得所需时刻各点的状态变量 和输出量 。()ktX()kY t第四章 控制系统数字仿真 4.63面向结构图的数字仿真 3仿真程序框图与实现基于以上仿真模型和求解思路，我们可以着手考虑面向这类典型闭环控制系统结构图的仿真程序的编制与实现。作为系统仿真程序，使用时应尽可能方便，使用者只要将开环传递函数 的分母、分子各系数 ，和反馈系数输入计算机，计算机就掌握了有关该系统的基本信息模型参数。然后求取首一化表达形式，以及形成开、闭环状态方程各矩阵等步骤均由仿真程序自动完成，无需人工干预。因此，程序应有输入数据模块和初始化程序模块。()G s(1,2,)jajn

36、L(1,2,1)ib inL第四章 控制系统数字仿真 4.64面向结构图的数字仿真 在参考输入函数 的作用下，系统输出 开始随时间变化，仿真程序应能按照给定的计算步长，采用已确定的数值算法，对系统中各状态变量和输出逐点变化情况进行求解运算。这部分模块是整个仿真程序的核心部分，计算速度、精度误差均取决于它。通常称为运行程序模块。在使用者规定的时间范围内，将计算数据按照一定要求存储，并在仿真结束时，按使用者指定格式输出仿真结果，以便对系统进行分析研究，这就是所谓的输出程序模块。综上所述，构成一个完整的仿真程序必须至少建立：输入数据块、初始化块、运行计算块、输出结果块四个子程序模块，才能正常地完成仿

37、真任务。()r t()y t第四章 控制系统数字仿真 4.65面向结构图的数字仿真 1)程序框图程序框图如图4.9所示。图4.9 典型闭环系统的仿真程序框图第四章 控制系统数字仿真 4.66面向结构图的数字仿真 2)程序语句用MATLAB语言编制的程序语句如下：EXP4-1.M%输入数据 aao,a1，an；%n1维分母系数向量 bb0,b1,.,bm；%ml维分子系数向量 XOx10,x20,.,xn0；%状态向量初值 VVo；%反馈系数 nno；%系统阶次 T0=t0；%起始时间 Tf=tf；%终止时间 h=h0；%计算步长 R=r%阶跃输入函数的幅值%形成开、闭环系数矩阵第四章 控制系统

38、数字仿真 4.67面向结构图的数字仿真B=b/a(1)；a=a/a(1)；A=a(2:n+1)%首一化处理A=rot90(rot90(eye(n-1,n)；-fliplr(A)；%形成能控标准形A阵B=zeros(1,n-1)，1；%形成输入矩阵B(n维列向量)m1=length(b)；%分子系数向量维数m+1C=fliplr(b),zeros(1,n-m1)；%形成输出矩阵C(n维行向量)=A-B*C*V；%形成闭换系数矩阵X=X0；y=0；t=T0；%设初值，准备开始递推运算%运算求解N=round(Tf-T0)/h)；%确定输出点数For i=1:N%四阶龙格-库塔法 K1=*X+B*R

39、 K2=*(X+h*K1/2)+B*R；K3=*(X+h*K2/2)+B*R；K4=*(X+h*K3)+B*R%求各系数K X=X+h*(K1+2*K2+2*K3+K4)/6%求状态 y=y,C*X；%求输出并以向量形式保存 t=t,t(i)+h；%输出对应时刻以列向量形式保存end%输出结果 t,y%输出数据形式结果 plot(t,y)%输出曲线形式结果 第四章 控制系统数字仿真 4.68面向结构图的数字仿真以上即为采用MATLAB语言实现面向典型闭环系统数字仿真的程序。上面程序中用到了几个MATLAB语言的特殊功能，使得程序较为简练。简要说明如下：eye(m,n)产生mn型 阵，即主对角元

40、素均为1，其余元素为零的长方阵length()求取括号内向量的维数zeros(m,n)产生mn型零阵rot90()将括号内矩阵左转90fliplr()将括号内矩阵或向量左右翻转round()对括号内运算取整plot(t,y)将打印为横轴，为纵轴的曲线 第四章 控制系统数字仿真 4.69面向结构图的数字仿真【例4.2】求图4.10所示系统的阶跃响应 的数值解。()y t图4.10 控制系统结构图第四章 控制系统数字仿真 4.70面向结构图的数字仿真 解：该系统结构形式为典型闭环控制系统，用EXP4-1.M求解过程如下：(1)取开环放大系数k=1，反馈系数v=1，阶跃输入幅值r=1。(2)利用MA

41、TLAB语言中conv()卷积函数功能，先将系统开环传递函数 化为要求的传递函数形式的分母、分子多项式，系数向量为 和 。()G s01,na aaL01,mb bbL(3)设系统状态向量初值 均为零。(4)系统运行参数 ，。10200,nxxxL04n 00t 10ft 00.25h 第四章 控制系统数字仿真 4.71面向结构图的数字仿真(5)按以上步骤和参数，在MATLAB语言环境下，输入以下命令语句：k=1；a=conv(1 0 0,conv(0.25 1,0.25 1)；a=0.0625 0.5000 1.0000 0.0000 0.0000 b=2*k k b=2 1X0=0 0 0

42、 0v=1；n=4；T0=0；Tf=10；h=0.25；R=1；EXP4-1.M%运行EXP4-1.M仿真程序Plot(t,y)%输出曲线形式结果(见图4.11)第四章 控制系统数字仿真 4.72面向结构图的数字仿真 图4.11 输出响应 曲线()y t第四章 控制系统数字仿真 4.73面向结构图的数字仿真 二.复杂连接的闭环系统数字仿真 实际工程中常常遇到的是复杂结构形式的控制系统，这些控制系统是由若干典型环节按照一定规律连接而成的，可用前述的形式将其描述为复杂连接的闭环系统结构图。要采用上节所述方法仿真，必须先将复杂结构图简化成为图4.8的典型结构式，求出开环传递函数，然后运用上节程序进行

43、仿真分析。但存在如下问题：系统结构复杂，内部存在交叉耦合、局部回环等情况时，化简并非易事，尤其通过手工化简使工作量陡增；有时在分析中，还需要得知结构图中某些环节的输出变量情况，若经过化简消去这些环节，则不便进行观察分析；在分析中常常需要改变某个参数，观察其对输出的影响，但改动一个参数值，就有可能需将所有开环传递函数分子、分母系数统统重新计算再次输入计算机，很不方便；对实际系统中存在的非线性情况，更是无法加以考虑。鉴于以上原因，编制面向复杂连接的闭环系统结构图的系统仿真程序，必须使其能克服上节仿真程序的不足之处且应该具有以下特点：(1)可按照系统结构图输入各环节参数，对应关系明确，改变参数方便。

44、(2)可方便地观察各环节输出动态响应。(3)各环节存在非线性特性时易于处理。第四章 控制系统数字仿真 4.74面向结构图的数字仿真 1典型环节的二次模型化复杂连接闭环系统数字仿真程序的基本思路是：与实际系统的结构图相对应，在计算机程序中也应能方便地给出表示各实际环节的典型环节，并将环节之间的连接关系输入计算机，由计算机程序自动形成闭环状态方程，运用数值积分方法求解响应。因此，选定典型环节是很重要的，要使其既具有代表性，又不致于造成输入数据复杂烦琐。考虑控制系统典型环节有以下几种情况。比例环节：()G sK积分环节：()KG ss第四章 控制系统数字仿真 4.75面向结构图的数字仿真 比例积分：

45、12()K sKG ss惯性环节：()1KG sTs一阶超前滞后环节：121()1T sG sKT s二阶振荡环节：22()21KG sT sTs第四章 控制系统数字仿真 4.76面向结构图的数字仿真 可见，除二阶振荡环节外，都是一阶环节，完全可用一个通用一阶环节(见图4.12)表示，即图4.12 典型一阶环节()()()iiiiiiY sCD sG sU sAB s(1,2,)inL(4.25)式中，为第i个环节的输出；为第i个环节的输入；n为系统中的环节数(也就是系统的阶次)。而二阶振荡环节可以化为如图4.13方式连接而成的等效结构图。()iy s()iu s第四章 控制系统数字仿真 4.

46、77面向结构图的数字仿真 由图可见，二阶振荡环节完全可以用一阶环节等效连接得到。因此，选定式(4.25)作为组成复杂系统仿真模型的典型环节是适合的。图4.13 二阶振荡环节的等效结构图第四章 控制系统数字仿真 4.78面向结构图的数字仿真 设：输入向量 ，其中各分量表示各环节输入量；T12,nu uuLU输出向量 ，其中各分量表示各环节输出量；T12,ny yyLY模型参数矩阵为1200nAAAOA1200nBBBOB1200nCCCOC1200nDDDOD于是系统中所有环节输出、输入关系用矩阵表示如下:()()ssAB YCD U第四章 控制系统数字仿真 4.79面向结构图的数字仿真 2系统

47、的连接矩阵与仿真求解1)连接矩阵仅把系统中各环节描述出来还不够，要进行数值积分求解，还必须把各环节之间的相互作用关系清楚地表达出来，这种表达是通过建立本章第4.2节所述的连接矩阵实现的。为说明问题，现举例说明如下：图4.14 某闭环系统结构图第四章 控制系统数字仿真 4.80面向结构图的数字仿真 设系统结构如图4.14所示，图中编号1、2、3表示三个典型环节，、为环节间直接作用的比例系数。各环节输入 与输出 的关系为23iuiy10uv33y 2122uyy32uy整理为矩阵形式得0UWYW V(4.27)第四章 控制系统数字仿真 4.81面向结构图的数字仿真 320010010W0100 W

48、式中,(4.28)第四章 控制系统数字仿真 4.82面向结构图的数字仿真 2)求系统的解将式(4.27)代入式(4.26)，则0()()()ssAB YCDWYW V整理，得00()()&BDW YCWA YCW VDW V为避免输入函数为阶跃函数时,初始点导数 ,为阶跃函数)为无穷大，要求使 ，则d(0divtt iv00DW0()()&BDW YCWA YCW V即0&QYRYCW V 第四章 控制系统数字仿真 4.83面向结构图的数字仿真 式中，QBDWRCWA若Q矩阵的逆矩阵 存在，则式(4.29)两边同时左乘 ，得1Q1Q110&YQ RYQ CW V即%&%YAYBV1%AQ R1

49、0%BQ CW式中，(4.30)第四章 控制系统数字仿真 4.84面向结构图的数字仿真 式(4.30)是典型的一阶微分方程组矩阵形式，利用4.3.1节介绍的求解方法，可方便地求出各环节的输出响应。然而，建立系统仿真模型还应注意以下两点：(1)存在的条件。为了分析 存在的条件，分析图4.15所示例子。1Q1Q图4.15 系统结构图第四章 控制系统数字仿真 4.85面向结构图的数字仿真 211AA200BB120CCC300DD0100 W001100010W230000000BDQBDW第四章 控制系统数字仿真 4.86面向结构图的数字仿真 由于Q阵中有两列全为零，所以 不存在。其原因在于它的第

50、一个环节为纯比例环节，第三个环节为纯微分环节，按典型环节给出的方程分别为1Q11 1233yC uyD u&显然不是微分方程而是代数方程，因而使 不存在。为避免系统中出现代数方程，必须保证系统中所有环节的 ，从而使环节能用微分方程给出。从数学上说，因为 ，保证了B非奇异，即使此时 ，也可保证Q非奇异。所以当系统中出现纯比例、纯微分这类使 的环节时，则应设法与其他环节合并处理，或设法化为系统可接受的环节。1Q0iB QBDW0iB 0DW0iB 第四章 控制系统数字仿真 4.87面向结构图的数字仿真 (2)00DW的条件 当被仿真系统的输入函数有阶跃函数时，阶跃函数在t=0点的 ,此时若 ,将使