控制工程基础-第二章

1、设计控制系统应完成哪些工作？控制对象运动规律的描述控制对象运动规律定性分析控制对象运动规律定量分析控制系统的设计与综合控制对象和控制系统的数学模型本章任务8/6/20221北京科技大学信息工程学院自动化系2、控制系统的数学模型2.1 控制系统的运动方程2.2 线性系统的频域模型2.3 方框图与信号流图2.4 状态空间与状态空间表达式2.5 控制系统不同模型间的关系小 结8/6/20222北京科技大学信息工程学院自动化系本章学习要点简单物理系统的微分方程的列写；非线性模型的线性化方法；传递函数和传递函数矩阵的概念；结构图和信号流图的变换与化简； 状态空间与状态空间表达式；控制系统不同模型形式及其

2、之间的转换。8/6/20223北京科技大学信息工程学院自动化系2.1 控制系统的运动方程例2.1.1研究RLC电路，试找出输出电压uc(t)随输入电压ur(t)变化的规律。解R、C、L以及初始uc(0)确定时,已知ur(t)就可以确定uc(t)8/6/20224北京科技大学信息工程学院自动化系2.1 控制系统的运动方程例2.1.2如图：由质量为m的木块、弹性系数为K的弹簧和阻尼系数为B的系统，试找出木块的位移x(t)与外力F(t)之间的关系。解m、K、B以及初始x(0)确定时,已知f(t)就可以确定x(t)8/6/20225北京科技大学信息工程学院自动化系2.1 控制系统的运动方程 直流他励电

3、动机电枢电路，取电枢电压ua为输入量,电动机角速度m为输出量,讨论它们之间的关系。电枢回路电压平衡方程：电磁转矩方程：ammtiCM)(=电动机轴上的转矩平衡方程：例2.1.3解电枢反电势是电枢电流产生的电动转矩是电动机转矩系数mMmC是折合到电动机轴上的总负载转矩(t)McJm:电动机和负载折合到电动机轴上的转动惯量；fm:电动机和负载折合到电动机轴上的黏性摩擦系数；8/6/20226北京科技大学信息工程学院自动化系2.1 控制系统的运动方程注意观察三个示例的微分方程可以通过求解得到ur(t)uc(t)，f(t)x(t)之间内在运动的关联关系、分析系统的运动特性。进而改造系统-选择适当的R、

4、L、C和m、B、K得到希望的运动规律。许多表面上看来似乎毫无共同之处的控制系统，其物理背景可能完全一样， 可以用一个运动方程来表示，我们可以不单独地去研究具体系统而只分析其数学表达式，即它们具有相同的数学模型。这类系统被称为相似系统。8/6/20227北京科技大学信息工程学院自动化系2.1 控制系统的运动方程控制系统的运动对系统施加控制（即输入控制信号），从而得到系统输出量（即受控量）随时间的变化规律（即输出响应信号）。控制系统的运动方程根据描述系统特性的物理学定律，如机械，电气，热力，液压等方面的基本定律写出。展示系统在运动过程中各变量之间的相互关系，既定性又定量地描述整个系统的运动过程。

5、数学模型描述系统内部物理量（或变量）之间的数学表达式，是分析和设计自动控制系统的基础。静态模型：在静态条件下（即变量不随时间变化），描述变量之间关系的代数方程(组)。动态模型：描述变量各阶导数之间关系的微分方程(组)。8/6/20228北京科技大学信息工程学院自动化系2.1 控制系统的运动方程建立数学模型的方法解析法 依据描述系统运动规律的运动定律来得到微分方程的方法。实验法 基于系统输入输出的实验数据来建立数学模型的方法。数学模型的形式时域模型微分方程、差分方程和状态方程；复频域模型传递函数、结构图、频率特性。8/6/20229北京科技大学信息工程学院自动化系2.1 控制系统的运动方程问题：

6、从严格意义上讲，绝大多数系统的数学模型都不是线性模型（即系统并非是线性系统）。事实上，任何一个元件总是存在一定程度的非线性。即使假设具有线性的特性，也是局限在一定的范围内。几种常见的非线性8/6/202210北京科技大学信息工程学院自动化系2.1 控制系统的运动方程两类非线性系统具有连续变化的非线性系统动态：y(n)=f(t;y,y(1),y(n-1),x,x(1),x(m)静态：y=f(x)本质非线性系统8/6/202211北京科技大学信息工程学院自动化系2.1 控制系统的运动方程 非线性微分方程的求解很困难。在一定条件下，近似地转化为线性微分方程，可以使系统的动态特性的分析大为简化。实践证

7、明，这样做能够圆满地解决许多工程问题，有很大的实际意义。线性化的方法忽略弱非线性环节：如果元件的非线性因素较弱或者不在系统线性工作范围以内，则它们对系统的影响很小，就可以忽略。台劳级数展开法(小偏差法，切线法，增量线性化法)：适用前提假设在控制系统的整个调节过程中，各个元件的输入和输出量只是在平衡点附近作微小变化。8/6/202212北京科技大学信息工程学院自动化系2.1 控制系统的运动方程忽略二次以上的各项，上式可以写成：A(x0,y0)平衡点，函数在平衡点处连续可微，则可将函数在平衡点附近展开成台劳级数：其中：非线性元件的线性化数学模型8/6/202213北京科技大学信息工程学院自动化系2

8、.1 控制系统的运动方程平均斜率法：如果一非线性元件输入输出关系如下图所示，此时不能台劳级数展开法，可用平均斜率法得线性化方程为 ：其中：8/6/202214北京科技大学信息工程学院自动化系2.1 控制系统的运动方程注意：这几种方法只适用于一些非线性程度较低的系统，对于某些严重的非线性（本质非线性)不能作线性化处理，一般用相平面法及描述函数法进行分析。(此部分超出本课程的内容，可参考非线性控制的章节或教材。)8/6/202215北京科技大学信息工程学院自动化系2.1 控制系统的运动方程例2.1.4 水位自动控制系统输入量为Q1 ,输出量为水位H，求水箱的微分方程。水箱的横截面积为C，R表示流阻

9、。8/6/202216北京科技大学信息工程学院自动化系根据托里拆利定理，出水量与水位高度平方根成正比，则有：其中 为比例系数。水箱的线性化微分方程:整理得水箱的标准线性化微分方程为:其中：显然这个式子为非线性关系, 在工作点 附近进行台劳级数展开。取一次项得：2.1 控制系统的运动方程解在 时间中，水箱内流体变化量 .则：注意： 是常数； 。8/6/202217北京科技大学信息工程学院自动化系2.1 控制系统的运动方程说明本质非线性系统一般不可线性化。多变量情况处理类似。工作点不同，所得线性化方程的线性化系数不同，即线性化方程不同。非线性系统的线性化方程只在工作点附近才成立 。8/6/2022

10、18北京科技大学信息工程学院自动化系2.2 线性系统的复数域模型2.2.1 拉普拉斯变换2.2.2 传递函数2.2.3 传递函数矩阵2.2.4 典型元部件及典型环节的传递函数问题： 1）微分方程求解比较困难，不利于工程实现； 2）有时分析控制系统的性质时不必求解方程； 是否有更方便的形式描述系统？8/6/202219北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型 拉氏变换法是一种数学积分变换，其核心是把时间函数f(t)与复变函数F(s)联系起来，把时域问题通过数学变换为复频域问题，把时间域的高阶微分方程变换为复频域的代数方程以便求解。对应 时域函数f(t)(原函数)复频域函数F(

11、s)(象函数)s为复频率2.2.1 拉普拉斯变换8/6/202220北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型正变换反变换正变换反变换象函数F(s) 用大写字母表示,如I(s)，U(s)。原函数f(t) 用小写字母表示，如 i(t), u(t)。12象函数F(s) 存在的条件：拉氏变换的定义t 0 , f(t)=08/6/202221北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型如果存在有限常数M和c使函数f(t)满足：则 总可以找到一个合适的s值使上式积分为有限值，即f(t)的拉氏变换式F(s)总存在。8/6/202222北京科技大学信息工程学院自动化系2.2

12、 线性系统的复频域模型典型函数的拉氏变换 (1)单位阶跃函数的象函数(2)单位冲激函数的象函数(3)指数函数的象函数8/6/202223北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型(4)正弦函数的象函数(5)余弦函数的象函数8/6/202224北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型拉普拉斯变换的基本性质线性性质时间比例性质（相似定理）其中为实常数8/6/202225北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型微分性质时域导数性质频域导数性质8/6/202226北京科技大学信息工程学院自动化系积分性质延迟性质频域延迟时域延迟在时间域的平移变换

13、在复数域有对应的衰减变换。时间信号f(t)在时间域的指数衰减，其拉氏变换在复数域有对应的坐标平移。2.2 线性系统的复频域模型8/6/202227北京科技大学信息工程学院自动化系初值定理f(t)和 的拉氏变换存在， 也存在，则终值定理f(t)和 的拉氏变换存在， ，并且除在原点处唯一的极点外，SF(S)在包含j轴的右半平面是解析的（即t时，f(t)为常数），则时域函数的初值，可以由变换域求得。时域函数的终值，也可以由变换域求得。2.2 线性系统的复频域模型8/6/202228北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.1已知微分方程如下，试求初值皆为零时输出量的拉氏变

14、换与输入的拉氏变换之比。解初值皆为零有由微分性质对微分方程作拉氏变换得：8/6/202229北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型拉普拉斯反变换的求法(1)按定义(2)对简单形式的F(s)可以查拉氏变换表得原函数(P28)f(t)F(s)f(t)F(s)(t)1Sint1(t)1/sCostt1/(s+a)8/6/202230北京科技大学信息工程学院自动化系(4)把F(S)分解为简单项的组合部分分式展开法(3)利用拉氏变换的性质解由延迟性质知：思考的原函数2.2 线性系统的复频域模型例2.2.28/6/202231北京科技大学信息工程学院自动化系利用部分分式可将F(s)分

15、解为：象函数的一般形式：待定常数12.2 线性系统的复频域模型8/6/202232北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型待定常数的确定：方法1方法2求极限的方法8/6/202233北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.3求如下象函数的原函数。解解法1解法2原函数的一般形式：8/6/202234北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型一对共轭复根为一分解单元,设：28/6/202235北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.4解8/6/202236北京科技大学信息工程学院自动化系2.2 线性系

16、统的复频域模型其中：38/6/202237北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.5解8/6/202238北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型 n =m 时将F(s)化成真分式和多项式之和小结:由F(s)求f(t)的步骤求真分式分母的根，确定分解单元将真分式展开成部分分式，求各部分分式的系数对每个部分分式和多项式逐项求拉氏反变换8/6/202239北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.6解8/6/202240北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型2.2.2 系统的传递函数（1）定义： 单

17、输入单输出线性定常动态对象的传递函数G(S)是零初值下该对象的输出量的拉普拉斯变换Y(S)数与输入量的拉普拉斯变换R(S)之比。 回答本节开始的问题8/6/202241北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型RLC电路取ur为输入，uc为输出，得:拉氏变换得：则传递函数为：例2.2.7解8/6/202242北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.8解根据牛顿第二定律，得取外力f(t)为输入；位移x(t)为输出得微分方程：拉氏变换后得：传递函数为：8/6/202243北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型 一般有nm 。同一

18、个系统，当输入量和输出量的选择不相同时，可能会有不同的传递函数。不同的物理系统可以有相同的传递函数。传递函数表示系统传递输入信号的能力，反映系统本身的动态性能。它只与系统的结构和参数有关，与外部作用等条件无关。（2）传递函数的性质 G(s)与系统的微分方程有直接联系。G(s)是系统单位脉冲响应的拉氏变换8/6/202244北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型（3）传递函数的常用表示形式 时间常数形式根的形式8/6/202245北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型8/6/202246北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型（4）

19、传递函数局限 G(s)原则上不反映y(0)0时的系统的全部运动规律. G(s)只适用于单输入，单输出系统。 G(s)只适用于线性定常系统由于拉氏变换是一种线性变换.8/6/202247北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型特征多项式：G(s)的分母多项式D(s)特征方程：D(s)=0极点/特征根：D(s)=0的根零点：N(s)=0的根零极点对消系统的阶数：max(n,m)，（一般nm）系统的类型放大系数与上述传递函数有关的几个重要概念：系统的放大系数 K根轨迹放大系数 Kg零极点图8/6/202248北京科技大学信息工程学院自动化系G(s) 的零点、极点表示在S平面上零极

20、点图G(s)G(s)零极点分布图 系统性能 G(s)2.2 线性系统的复频域模型8/6/202249北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型2.2.3 传递函数矩阵将传递函数的概念推广到多输入多输出系统，传递函数G（s）推广为传递函数矩阵G（s）。设系统有p个输入量、q个输出量如下图。u1u2up。y1y2yq。8/6/202250北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型例2.2.9如图，直流他励电动机；ua是外加的输入变量电枢电压（伏），m表示电动机的角转速（弧度/秒），为输出量。讨论它们之间的关系。8/6/202251北京科技大学信息工程学院自动化系

21、2.2 线性系统的复频域模型由例2.1.3系统运动方程为：解拉氏变换得：整理得：8/6/202252北京科技大学信息工程学院自动化系 比例环节 控制系统通常由若干个基本部件组合而成，这些基本部件称为典型环节。 包括：比例环节、微分环节、积分环节、比例微分环节、一阶惯性环节、二阶振荡环节和延迟环节。（1） 典型环节2.2 线性系统的复频域模型2.2.4 典型元部件及典型环节的传递函数8/6/202253北京科技大学信息工程学院自动化系比例环节的单位阶跃响应当 时2.2 线性系统的复频域模型8/6/202254北京科技大学信息工程学院自动化系 一阶惯性环节当 时，微分方程是一阶的,且输出响应需一定

22、的时间才能达到稳态值。其中T为惯性环节的时间常数。2.2 线性系统的复频域模型8/6/202255北京科技大学信息工程学院自动化系惯性环节的单位阶跃响应求拉氏反变换得 2.2 线性系统的复频域模型8/6/202256北京科技大学信息工程学院自动化系 积分环节其中K=1/T ， T为积分环节的时间常数，表示积分的快慢程度。积分环节的单位阶跃响应2.2 线性系统的复频域模型8/6/202257北京科技大学信息工程学院自动化系 微分环节其中K为微分环节的时间常数，表示微分速率的大小。2.2 线性系统的复频域模型 理想微分环节 一阶微分环节（又称比例微分环节、实用微分环节） 8/6/202258北京科

23、技大学信息工程学院自动化系 二阶振荡环节这种环节包括有两个储能元件，当输入量发生变化时，两种储能元件的能量相互交换。在阶跃函数作用下，其暂态响应可能作周期性的变化。 式中： 自然振荡角频率 阻尼比由二阶微分方程描述的系统。2.2 线性系统的复频域模型8/6/202259北京科技大学信息工程学院自动化系当输入量为阶跃函数时，输出量的拉氏变换为：当 时，上式特征方程的根为共轭复数。因式分解得：振荡环节的单位阶跃响应：输出量为 ：2.2 线性系统的复频域模型8/6/202260北京科技大学信息工程学院自动化系 延迟/时滞环节2.2 线性系统的复频域模型带钢厚度检测环节写成一般形式 :零初始条件下，拉

24、氏变换为 传递函数为 例8/6/202261北京科技大学信息工程学院自动化系时滞环节的输出量时滞环节的传递函数 对于时滞时间很小的时滞环节，常把它展开成泰勒级数，并略去高次项，得：时滞环节在一定条件下可近似为惯性环节！2.2 线性系统的复频域模型8/6/202262北京科技大学信息工程学院自动化系2.2 线性系统的复频域模型2.2.4 典型元部件及典型环节的传递函数元部件名称传递函数电位器测速电机电加热炉单容水槽双容水槽（2） 典型元部件(有纯延迟)(也可有延迟，略)8/6/202263北京科技大学信息工程学院自动化系2.3 方框图与信号流图控制系统的结构图是由许多对信号进行单向运算的方框和一

25、些信号流向线组成，它包含4种基本单元。2.3.1 系统动态结构图 1）信号线 2）引出点（或测量点）3）比较点（或综合点）4）方框（或环节）8/6/202264北京科技大学信息工程学院自动化系2.3 方框图与信号流图8/6/202265北京科技大学信息工程学院自动化系2.3 方框图与信号流图思考:将两部分电路分开分别讨论然后在结合到一起结果和前面得到的是否相同?8/6/202266北京科技大学信息工程学院自动化系2.3 方框图与信号流图速度控制系统例2.3.18/6/202267北京科技大学信息工程学院自动化系2.3 方框图与信号流图解（1）比较环节和速度调节器环节式中：式中：8/6/2022

26、68北京科技大学信息工程学院自动化系2.3 方框图与信号流图式中整理得 8/6/202269北京科技大学信息工程学院自动化系2.3 方框图与信号流图（2）速度反馈的传递函数 式中： 为速度反馈系数 8/6/202270北京科技大学信息工程学院自动化系2.3 方框图与信号流图（3）电动机及功率放大装置8/6/202271北京科技大学信息工程学院自动化系2.3 方框图与信号流图（4）系统的动态结构图 8/6/202272北京科技大学信息工程学院自动化系2.3 方框图与信号流图2.3.2 系统的等效变换（1）典型连接的等效传递函数 串联8/6/202273北京科技大学信息工程学院自动化系2.3 方框

27、图与信号流图并联8/6/202274北京科技大学信息工程学院自动化系2.3 方框图与信号流图反馈连接8/6/202275北京科技大学信息工程学院自动化系2.3 方框图与信号流图（2）相加点及分支点的换位运算原则:换位前后的输入/输出信号间关系不变。 相加点后移8/6/202276北京科技大学信息工程学院自动化系2.3 方框图与信号流图相加点前移分支点后移8/6/202277北京科技大学信息工程学院自动化系2.3 方框图与信号流图分支点前移 分支点换位 8/6/202278北京科技大学信息工程学院自动化系2.3 方框图与信号流图相加点变位8/6/202279北京科技大学信息工程学院自动化系2.3

28、 方框图与信号流图相加点和分支点一般不能变位 8/6/202280北京科技大学信息工程学院自动化系2.3 方框图与信号流图(3)系统开环传递函数 定义：闭环系统反馈信号的拉氏变换与偏差信号的拉氏变换之比（反馈通道断开），定义为系统的开环传递函数，用 表示。 8/6/202281北京科技大学信息工程学院自动化系2.3 方框图与信号流图系统的开环传递函数是正向通道传递函数与反向通道传递函数的乘积。 正向通道传递函数 反向通道传递函数8/6/202282北京科技大学信息工程学院自动化系2.3 方框图与信号流图利用方块图变换法则 (a) 相加点A前移，分支点D后移8/6/202283北京科技大学信息工

29、程学院自动化系2.3 方框图与信号流图(b) 消除局部反馈回路 8/6/202284北京科技大学信息工程学院自动化系2.3 方框图与信号流图（C） 消除主反馈回路 方块图的化简方法不是唯一的，应充分地利用各种变换技巧，选择最简捷的路径，以达到省力省时的目。 8/6/202285北京科技大学信息工程学院自动化系2.3 方框图与信号流图无交叉局部反馈系统 例2.3.2解8/6/202286北京科技大学信息工程学院自动化系2.3 方框图与信号流图有交叉局部反馈系统 例2.3.3解8/6/202287北京科技大学信息工程学院自动化系2.3 方框图与信号流图(4)系统闭环传递函数 定义： 在初始条件为零

30、时，系统的输出量与输入量的拉氏变换之比称为系统的闭环传递函数，用 表示。 对于单位反馈系统，有 8/6/202288北京科技大学信息工程学院自动化系2.3 方框图与信号流图(5)系统对给定作用和扰动作用的传递函数 原则：对于线性系统来说，可以运用叠加原理，即对每一个输入量分别求出输出量，然后再进行叠加，就得到系统的输出量。8/6/202289北京科技大学信息工程学院自动化系2.3 方框图与信号流图只有给定作用 只有扰动作用 8/6/202290北京科技大学信息工程学院自动化系2.3 方框图与信号流图两个输入量同时作用于系统 8/6/202291北京科技大学信息工程学院自动化系2.3 方框图与信

31、号流图信号流图是一种用图线表示线性系统方程组的方法。设一组线性方程式如下：信流图的表示形式2.3.3 信号流图8/6/202292北京科技大学信息工程学院自动化系2.3 方框图与信号流图信号流图中的术语: (1)源点 ：只有输出支路的节点称为源点或称为输入节点。它一般表示系统的输入变量。(2)汇点 ：只有输入支路的节点称为汇点或称为输出节点。它一般表示系统的输出变量。 (3)混合节点 ：既有输入支点又有输出支点的节点称为混合节点。(4)通路：从某一节点开始，沿支路箭头方向经过各相连支路到另一节点（或同一节点）构成的路径，称为通路。通路中各支路传输的乘积称为通路传输（通路增益）。8/6/2022

32、93北京科技大学信息工程学院自动化系2.3 方框图与信号流图(5)开通路：与任一节点相交不多于一次的通路称为开通路。 (6)闭通路：如果通路的终点就是通路的起点，并且与任何其他节点相交不多于一次的通路称为闭通路或称为回环。(7)回环增益：回环中各支路传输的乘积称为回环增益（或传输）。(8)前向通路：是指从源点开始并终止于汇点且与其他节点相交不多于一次的通路，该通路的各传输乘积称为前向通路增益。 (9)不接触回环：如果一信号流图有多个回环，各回环之间没有任何公共节点，就称为不接触回环，反之称为接触回环。 8/6/202294北京科技大学信息工程学院自动化系2.3 方框图与信号流图梅逊增益公式：

33、式中： P 系统的总传输增益；pk 第k 条前向通道的传输增益；n 从输入节点到输出节点的前向通路数；信号流图的特征式。k与第k条前向通道不接触的那部分信号流图的；8/6/202295北京科技大学信息工程学院自动化系2.3 方框图与信号流图特征式的意义为：信号流图中所有不同回环的传输之和；信号流图中每两个互不接触回环的传输乘积之和；m个互不接触回环的传输乘积之和；称为第k条通路特征式的余因子，是在中除去第k条前向通路相接触的各回环传输（即将其置零）。 8/6/202296北京科技大学信息工程学院自动化系2.3 方框图与信号流图两个前向通道：P1，P2三个回路：La，Lb，Lc一对不相交回路：L

34、a，Lc例2.3.4解8/6/202297北京科技大学信息工程学院自动化系2.3 方框图与信号流图8/6/202298北京科技大学信息工程学院自动化系2.3 方框图与信号流图求：例2.3.5解（1）（2）8/6/202299北京科技大学信息工程学院自动化系2.3 方框图与信号流图（3）（4）8/6/2022100北京科技大学信息工程学院自动化系R(s)C(s)L1= G1 H1L2= G3 H3L3= G1G2G3H3H1L4= G4G3L5 = G1G2G3L1L2= (G1H1) (G3H3) = G1G3H1H3L1L4=(G1H1)(G4G3)=G1G3G4H1 G4(s) H1(s)

35、H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s)H3(s) G2(s) G3(s) G4(s) H1

36、(s)H3(s) G1(s) G2(s) G3(s) G1(s) G2(s) G3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s)G4(s)G3(s)梅逊公式例 H1(s)H3(s) G1(s) G2(s) G4(s) H1(s)H3(s) G1(s) G2(s) G3(s)P2= G4G3P1=G1G2G31=12=1+G1H18/6/2022101北京科技大学信息工程学院自动化系2.3 方框图与信号流图用梅森公式法求C(s)/R(s)、 E(s)/R(s)

37、、C(S)/N(S)。例2.3.6解(1)求C(s)/R(s) ：8/6/2022102北京科技大学信息工程学院自动化系2.3 方框图与信号流图两两互不接触回路有L1L2(2)若以E(s) 为输出，R(s) 为输入，求E(s)/R(s)：两两互不接触回路仍为L1L2无论输入输出是什么，回路是不变的，所以不变.8/6/2022103北京科技大学信息工程学院自动化系2.3 方框图与信号流图(3)求C(s)/N(s)：(4).若求R(s),N(s) 同时作用下的总输出：8/6/2022104北京科技大学信息工程学院自动化系2.3 方框图与信号流图例：利用梅逊公式，求：C(s)/R(s)解：画出该系统

38、的信号流程图 8/6/2022105北京科技大学信息工程学院自动化系2.3 方框图与信号流图该系统中有四个独立的回路： L1 = -G4H1 L2 = -G2G7H2 L3 = -G6G4G5H2 L4 = -G2G3G4G5H2互不接触的回路有一个L1 L2。所以，特征式 =1-（L1 + L2 + L3 + L4）+ L1L2该系统的前向通道有三个： P1= G1G2G3G4G5 1=1 P2= G1G6G4G5 2=1 P3= G1G2G7 3=1-L1 因此，系统的闭环系统传递函数C(s) / R(s)为8/6/2022106北京科技大学信息工程学院自动化系2.3 方框图与信号流图例：

39、画出信号流图，并利用梅逊公式求取它的传递函数C(s) / R(s)。信号流图：8/6/2022107北京科技大学信息工程学院自动化系2.3 方框图与信号流图注意：图中C位于比较点的前面，为了引出C处的信号要用一个传输为1的支路把C、D的信号分开。单独回路有L1、L2和L3，互不接触回路有 L1L2，即 ：前向通路只有一条，即 所以 8/6/2022108北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式问题：（1）系统运动方程是时变的传递函数不存在？（2）多输入多输出系统传递函数矩阵不好处理？（3）能否反映系统的中间变量的变化？（4）非线性系统如何描述？高阶线性微分方程求解方法：

40、（1）拉氏变换（2）一阶线性微分方程组8/6/2022109北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式例2.4.1研究RLC电路，试找出输出电压uc(t)随输入电压ur(t)变化的规律。解8/6/2022110北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式例2.4.2如图：由质量为m的木块、弹性系数为K的弹簧和阻尼系数为B的系统，试找出木块的位移x(t)与外力f(t)之间的关系。解8/6/2022111北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式在例2.4.2中取：x1(t)=x(t), x2(t)=v(t), u(t)=f(t),y

41、(t)=x(t)在例2.4.1中取：x1(t)=uc(t), x2(t)=i(t)，u(t)=ui(t),y=uc8/6/2022112北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式状态变量一组能够完全表征系统运动状态的相互独立的最小个数的变量。 x1(t), x2(t), xn(t) 状态向量以状态变量为分量构成的向量，维数与状态变量的个数相同，一般等于系统中储能元件的个数。 xT(t)=(x1(t), x2(t), xn(t)状态空间以状态变量x1(t), x2(t), xn(t)为坐标轴构成的欧氏空间。状态空间的基本概念8/6/2022113北京科技大学信息工程学院自动

42、化系2.4 状态空间与状态空间表达式状态变量的特点独立性：状态变量之间线性独立；多样性：状态变量的选取并不唯一，实际上存在无穷多种方案；等价性：两个状态向量之间只差一个非奇异线性变换；现实性：状态变量通常取为涵义明确的物理量；抽象性：状态变量可以没有直观的物理意义。8/6/2022114北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式状态方程和输出方程-状态空间表达式状态方程输出方程线性系统状态空间表达式n维系统的状态方程的标准形式写成向量形式8/6/2022115北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式输出方程的标准形式写成向量形式：y=g（x，u，t

43、）状态空间表达式的标准形式8/6/2022116北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式与状态空间表达式相关的术语8/6/2022117北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式加法器： 积分器：放大器：状态空间表达式的方框图表示与模拟结构图表示8/6/2022118北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式 B(t)C(t)A(t)D(t)表示系统信号的传递关系+8/6/2022119北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式模拟结构图用来反映系统状态变量间的信息传递关系uy-8/6/2022120北京

44、科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式状态空间表达式建立方法由物理机理分析直接建立由微分方程建立由传递函数建立（例2.4.1，2.4.2）（1）由物理机理分析直接建立的步骤由牛顿定理，基尔霍夫定律等写出运动方程选择适当状态变量写出状态空间表达式8/6/2022121北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式（2）由微分方程建立的方法微分方程不含有输入项的导数项微分方程含有输入项的导数项8/6/2022122北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式微分方程不含有输入项的导数项整理得能控标准型8/6/2022123北京科技大学信息

45、工程学院自动化系2.4 状态空间与状态空间表达式另一种取状态变量的方法所以有8/6/2022124北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式整理得能观标准型8/6/2022125北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式微分方程含有输入项的导数项8/6/2022126北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式uyx1x2x38/6/2022127北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式系统的状态空间表达式需要确定系数8/6/2022128北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式用代

46、数方法确定待定系数。其中x4辅助状态带入方程后 整理有解此方程即可得 8/6/2022129北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式状态空间表达式为：推广到n维：8/6/2022130北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式由微分方程所描述的系统，试写出其状态空间表达式。a2=18 , a1=192 , a0=640b3=0 , b2=0 , b1=160 , b0=640例2.4.2解3=b3=02=b2-a23=01=b1-a22-a13=1600=b0-a21-a12-a03 =640-18*160 =-2240系统的状态空间表达式为：8/

47、6/2022131北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式另一种实现形式：8/6/2022132北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式能控标准型：能观标准型：8/6/2022133北京科技大学信息工程学院自动化系2.4 状态空间与状态空间表达式（3）由传递函数建立传递函数微分方程并联分解(极点互异的情况)当 p1pn两两互异时，传函可分解为：其中 8/6/2022134北京科技大学信息工程学院自动化系uy此时，系统的状态模拟图可绘为：2.4 状态空间与状态空间表达式8/6/2022135北京科技大学信息工程学院自动化系如图所选状态，则有： p1

48、c1pncnxnx1uy2.4 状态空间与状态空间表达式约当标准型8/6/2022136北京科技大学信息工程学院自动化系已知系统传函为： 求状态空间表达式。 将G(s)分解为： 求以上4个系数： 取 解2.4 状态空间与状态空间表达式例2.4.48/6/2022137北京科技大学信息工程学院自动化系整理有： 2.4 状态空间与状态空间表达式约当标准型8/6/2022138北京科技大学信息工程学院自动化系现在考虑运动对象的状态方程与传递函数之间的关系。设某p输入q输出的对象的状态方程和输出方程如下： 其中A为nn的矩阵。在零初值条件下，对式（2.5.1）取拉普拉斯变换，得sX(s)=AX(s)+BU(s). 移项，得(sIn-A) X(s)=BU(s).其中多项式矩阵（sIn-A）必为非奇异。2.5 控制系统不同模型间的转换 8/6/2022139北京科技大学信息工程学院自动化系可求出 X(s)= (sIn-A)-1BU(s).对式（2.5.2）取拉普拉斯变换后以上式代入，即得 Y(s)= (C(sIn-A)-1B+D)U(s),与式Y(s)=G(s)U(s)比较，就知道该系统的传递函数是 G(s)=C(sIn-A)-1B+D. （2.5.3）许多情况下D=0，则有 G(s)=C(sIn-A)-1B. （2.5.4）式（2.5.3）和式（2.5.4）给出了状态空间描述下的传递