动力学与控制基础-共103页课件.ppt

1、动力学和控制基础动力学和控制基础陈 立 群上海大学力学系，上海市应用数学和力学所lqchenstaff.shu.edu提纲提纲1 引言引言2 工程背景工程背景3 系统建模系统建模4 系统分析系统分析5 系统控制系统控制6 总结总结引言动力学与控制动力学与控制研究随时间变化过程及其调节的力学分支学科。主要内容主要内容系统的建模、分析和控制2.1 车辆中的力学问题车辆中的力学问题设：车轮、车身作平面运动设：车轮、车身作平面运动车身和车轮的运动车身和车轮的运动不仅仅是平面运动不仅仅是平面运动 车身作什么运动？车身作什么运动？ 车轮作什么运动？车轮作什么运动？工程背景汽车的减振测试与疲劳测试汽车的减振

2、测试与疲劳测试越野赛车越野赛车问题的研究是：问题的研究是： 从简单到复杂；从简单到复杂； 从特殊到一般；从特殊到一般； 从单一到综合。从单一到综合。工程背景1893年生产的轿车年生产的轿车车轴与车体之间无减振器车轴与车体之间无减振器1904年生产的轿车年生产的轿车车轴与车体之间有减振器车轴与车体之间有减振器汽车发展的要求：汽车发展的要求：便捷、舒适、安全、环保便捷、舒适、安全、环保研究对象：研究对象：多个物体组成，结构更加复杂多个物体组成，结构更加复杂 减振结构：独立悬架减振结构：独立悬架工程背景现代研制的轿车、吉普车减振结构：现代研制的轿车、吉普车减振结构：独立悬架独立悬架独立悬架独立悬架共

3、轴式悬架共轴式悬架工程背景研制方法：研制方法：计算机的引入计算机的引入工程背景模拟实验和物理实验对比模拟实验和物理实验对比车辆碰撞的计算机模拟实验车辆碰撞的计算机模拟实验工程背景工程背景机器人骑自行车机器人骑自行车工程背景摩托车行驶失稳的现象摩托车行驶失稳的现象工程背景航天器对接航天器对接2.2 飞行器的动力学问题飞行器的动力学问题工程背景工程背景飞机起落架的飞机起落架的动力学仿真动力学仿真工程背景卫星太阳翻版展开的动力学仿真卫星太阳翻版展开的动力学仿真问题问题1：用什么方法建立系用什么方法建立系统的动力学方程便统的动力学方程便于编程计算？于编程计算？工程背景问题问题2：用什么方法定性和定量地

4、验证计算结果的正确性？用什么方法定性和定量地验证计算结果的正确性？系统建模3.1 概述概述建模是指导出时间变化系统控制方程的过程。建模先要对系统的构成进行分析，尽可能考虑系统的主要因素，再根据具体问题确定要采用的方法。完全根据力学理论的建模称为原理建模。主要包括能量法(Lagrange方程、变分原理等)和矢量法(Newton-Euler法)两大类。以实验作为主要方法的建模称为数据建模。主要包括模型辨识和参数估计等。系统建模3.2 Lagrange方程方程1736年年1月月25日生于意大利的都灵日生于意大利的都灵,1813年年4月月10日在法国巴黎去世，日在法国巴黎去世，19岁当数学教授，岁当数

5、学教授，为变分法奠定了理论基础，为变分法奠定了理论基础， 是当时欧洲公认的第一流的数学家。是当时欧洲公认的第一流的数学家。 他写的他写的分析力学分析力学（1788出版）一书运用出版）一书运用变分原理变分原理和和分析的方法，建立起完整和谐的力学体系，使力学分析化。分析的方法，建立起完整和谐的力学体系，使力学分析化。系统建模1769年蒸汽机车的模型年蒸汽机车的模型问题的引出问题的引出系统建模问题：问题：在不计摩擦的情况下，用什么方法建在不计摩擦的情况下，用什么方法建立系统的运动与主动力的关系立系统的运动与主动力的关系? ?质点系或刚体系的动力学问题质点系或刚体系的动力学问题系统建模设：质点系中第

6、i 个质点的质量为 ;作用在其上的有主动力 ; 约束力 . 质点的惯性力为 imiFiNFi IF应用dAlembert原理：), 1( ,INniiii0FFF0)(1INniiiiiWrFFF应用虚位移原理：若质点系所受的约束为理想约束01NniiirF0)(1N1IniiiniiiirFrFF其中：iiim aFI), 1( , 0)(INniiiiirFFF0)(1IniiiirFF动力学普遍方程系统建模0)(1IniiiirFFLagrange形式的形式的dAlembert原理原理0)()()(111niiiiizniiiiiyniiiiixzzmFyymFxxmF iixiyizI

7、iiiiiiiiFFFmxmymzxyz FijkFijkrijk受有理想约束的质点系，在运动过程中，其上所受的主动力和惯性力在质点系的任何虚位移上所作的虚功之和为零。动力学普遍方程的直角坐标形式系统建模0)(1niiiiimraF的独立性，有：利用), 1(kjqj设：具有完整理想约束的非自由质点系有 k 个自由度 系统的广义坐标为：kqqq,21T 为系统的动能，可表示成：),(11tqqqqTTkkjQ为对应于广义坐标jq的广义力niikjjjWqQ11)(Fd,(1, )djjjTTQjktqqkjjjjjqQqTqTt10dd如何推导？系统建模再对广义速度jq 求偏导数，得1( ,)

8、,iikt qqrr 因对时间t求导数，1kiiilllqqtrrr得(9-4)iijjqqrr (9-5)n iijjqqrr证明:n个质点，s个完整约束，k3ns，12 (1,2)iik=q ,q ,.q ,ti,.,nrr系统建模ddiijjtqqrr 2)将式(9-4)两边对广义坐标jq证明:求偏导数，有1kiiilljljjqqq qqt rrr而1d()()dkiiilljljjqtqqqqtrrr比较以上两式，可得d()diijjtqqrr(9-6)系统建模jn*iiQii=1j-F=mqrvddiiiiiijjmmtqqvvvvddddiiiiiijjmmtqtqrrvv22d

9、11d22iiiijjmmtqq vvddjjTTtqq系统建模Lagrange方程几种形式), 2 , 1(kj1、当所有主动力均为有势力时jjjQqTqTtddjkjqqqVQ),(1设：LT-V （ Lagrange函数）jjjqVqTqTtdd0)(ddjjqVTqTt0ddjjqLqLtjjjjLTVTqqqq=0系统建模2、当部分主动力为有势力时1),(jjkjQqqqVQddjjjQqLqLt), 2 , 1(kjjjjQqTqTtdd设：LT-V （Lagrange函数）应用Lagrange方程建立系统动力学方程的基本步骤：1、确定系统的自由度和广义坐标；2、用广义速度和广义坐

10、标给出系统的动能和势能；3、给出系统的拉格朗日函数；4、确定系统的非有势力的广义力。系统建模 xg1mg2mAB0l)(tFk解：解：1、求系统的动能和势能 ( Lagrange函数 ) 例1：图示机构在铅垂面内运动，均质杆AB用光滑铰链与滑块连接，求系统运动微分方程。AB2LAv cCAvCAACAxvvvv , cosLxvCxsinLvCy22221212121CCAJvmvmT222222132cos)(21LmLxmxmm2221)cos1 (kxgLmV( , , )( , )( , , , )LT xV xL xx 系统建模 xg1mg2mAB0l)(tFk0sincos)2(3

11、1)(sincos)(222222221gLmxLmLmtFkxLmLmxmm xkiixWQxQ1)(F)(tFQx0QddjjjQqLqLt222222132cos)(21LmLxmxmmVTL2221)cos1 (kxgLm2 2、求、求非有势力非有势力的广义力的广义力3 3、建立系统运动微分方程、建立系统运动微分方程方程的物理意义？方程的物理意义？12qxq系统建模 对于具有完整理想约束的质点系，若系统的自由度为k，则系统的动力学方程为：其中：VTLT：为系统的动能，V：为系统的势能ddjjjQqLqLt), 1(kjjQ：为对应于广义坐标 的非有势力的广义力jq当系统为保守系统时，有

12、：1：若系统存在循环坐标 ，则：qconst.pqTqL2：若系统的拉格朗日函数不显含时间t，则：const.02VTT系统建模3.3 Hamilton方程简介方程简介), 2 , 1( ,), 2 , 1( ,kjqHpkjpHqjjjj对于具有理想约束的质点系，其主动力为有势力，若设：111,(1,2, ),(, )jjkjjkkjTpjkqHp qLH pp qq tHamilton方程是关于广义坐标和广义动量的一阶微分方程组。方程是关于广义坐标和广义动量的一阶微分方程组。则系统的Hamilton方程为：系统建模对于定常约束的质点系，若主动力为有势力，则哈密顿函数H就是系统的动能与势能之

13、和，即：2HTVHamilton方程为数值计算提供了很好的微分-代数结构，在此基础上建立的辛算法可保持长期数值计算的稳定性。20TTV2HTV1,(1,2, ),kjjjjjTpjkHp qLq1kjjjTHqLq212102()TTTTTV2101()kjjjTTTqLq对于定常约束的质点系，有：00T 系统建模 在研究星球的运动轨道时，太阳系可视为Hamilton系统，其动力学方程可表示成：), 2 , 1( ,), 2 , 1( ,kjqHpkjpHqjjjj问题：如何精确地计算行星的运动轨迹，准确预测行星位置?解决问题的方法：提高计算方法精度和速度、通过数值仿真预测行星的运动轨迹和位置

14、，从而估计小行星撞击地球的可能性。k=3n, n为行星的个数(=9大行星+近百个小行星)系统建模冯 康(1920.91993.8) 数学与物理学家、计算数学家。1944年毕业于重庆中央大学物理系。19511953年赴前苏联进修。 曾任中国数学会理事，计算数学分会副理事长，中国计算机学会副主任等职。 1980年被选为中国科学院学部委员(数学物理学部院士)。 在拓扑代数、广义函数和计算数学等领域取得多方面首创性成就，并对我国计算机事业的创建和发展做出了重要贡献。 20世纪80年代，他提出了Hamilton系统的辛算法。该算法可保持长期数值计算的稳定性。系统建模 20世纪90年代末，中国科学院计算数

15、学研究所，秦孟兆院士领导的课题组，用Hamiltom系统的辛算法，预测小行星撞击地球的可能性是： 50年内不会发生 50年后，即使有小行星撞击地球的可能性，那时侯人类的科技手段一定能够阻止灾难的发生。系统分析3.1 概述概述分析是确定系统状态变化定性或定量特征的过程。除采用实验方法进行分析以外，常用的理论分析方法包括：几何方法、解析方法和数值方法。 几何方法是对系统作定性分析的方法。不仅能得到直观的定性结果，而且可为其它研究方法提供理论依据。 解析方法是对系统作定量分析的方法。不仅能确定系统运动随时间变化的规律，而且还能得到运动特性对系统参数的依赖关系。 数值方法是对系统作定量计算的方法。通过

16、数值求解非线性微分方程得到系统在特定的参数条件和初始条件的运动规律。 系统分析3.2 单自由度线性系统分析单自由度线性系统分析问题：对已知F(t)，求解常微分方程初值问题 00; (0), (0)mxcxkxF txxxx思路：线性系统叠加原理000; (0), (0)mxcxkxxxxxn0n00ddd( )e(cossin)txxx txtt2ndn,12kcmmk 其中 ; (0)0, (0)0mxcxkxF txx将激振力 F(t) 看作一系列脉冲力的叠加系统分析对于脉冲激励情形，系统只有暂态响应而不存在稳态响应 单位脉冲力可利用狄拉克（Dirac）分布函数(t) 表示 函数也称为单位

17、脉冲函数，定义为： )()(0)( ttt且()d1tt)(t的图象用位于时刻、长度为 1 的有向线段表示 )(tt10系统分析函数： )()(0)( ttt1)(dtt)(t 是一个广义函数可以看作矩形脉冲、脉冲面积为 1 而脉冲宽度趋于零时的极限 即 )(lim)(0 tt=为其它ttt 0)(1其中)(t也可以定义为其它形状的面积为 1 的脉冲 量纲：1/秒 )(tt10)(tt01系统分析 函数的性质: ( ) ()d( )f tttf特别地，当时刻 = 0 时，有( ) ( )d(0)f tttf实际应用时，通常 f (t) 在0t时才有意义0( ) ()d( )tf tttf冲量为

18、0I的脉冲力可借助函数表示为 )()(0tItF 当 I0 =1 时，为单位脉冲力 因而 )(1)(ttF系统分析现求处于零初始条件下的系统对单位脉冲力的响应单位脉冲响应记： 0+、0- 为单位脉冲力的前后时刻 运动微分方程与初始条件可合写为 0)0(, 0)0()(1xxtkxxcxm 或脉冲响应 ddd( )dmx tcx tkx ttt乘dt dx dd0 x tx在脉冲力作用的瞬间，位移来不及变化，但速度可产生突变 0 xdx令：令：d( )dmxtt)(tt00-0+系统分析d( )dmxtt两边在区间-00t内对时间积分 0000( )ddttmx t)0()0(1xmxmmx1)

19、0(在单位脉冲力的作用下，系统的速度发生了突变，但在这一瞬间，位移则来不及有改变，即有：x(0+) = x(0-) 又当 t 0+ 时，脉冲力作用已经结束，所以 t 0+ 时，有 mxxkxxcxm1)0(, 0)0(0 质量越大， 越小)0( x 质量越小， 越大)0( x 系统分析 mxxkxxcxm1)0(, 0)0(0 系统的单位脉冲响应即初始位移为零，而初始速度为 1/m 的自由振动 记为 h(t) nd1( )( )esin tdx th ttm无阻尼系统 nn( )( )sin x th ttm若单位脉冲力不是作用在时刻 t = 0，而是作用在 t =时刻 n()ddesin()

20、 th ttm响应响应如果系统在 t=时刻受到冲量为 I0 的任意脉冲力作用，则系统暂态响应可用脉冲响应函数表示为 tthItx),()(0n0n00ddd( )e(cossin)txxx txtt系统分析当处于零初始条件的系统受到任意激振力时，可以将激振力 F(t) 看作一系列脉冲力的叠加 对于时刻 t =的脉冲力系统受脉冲作用后产生速度增量mdF)( 并引起 t 各个时刻的响应 系统的脉冲响应 d( ) ()dxFh tdt)(tF其冲量为( )dF由线性系统的叠加原理，系统对任意激振力的响应应等于系统在时间区间t0内各个脉冲响应的总和 0( )( ) ()dtx tFh tDuhamel

21、积分 n()d0d1( )esin()dttFtmn()dd1esin() th ttm系统分析0( )( ) ()dtx tFh t利用卷积性质 0( )() ( )dtx tF th初始条件nn0n 00ddd()d0d( )e(cossin)1( )esin()d tttxxx txttFtm若阻尼为零00nnn0nn1( )(cossin)( )sin()txx tx t tF t dmn()d0d1( )esin()dttFtm系统分析例2：有阻尼系统在过渡阶段对简谐激励的响应 000)0(,)0(sinxxxxtFkxxcxm nn0n00dddnddd( )e(cossin)es

22、incos( sincos )sinsin() ttxxx txttBtstBtnkmkmc22dn1nskFB0222)2()1 (1ss122tan1ss初始条件响应自由伴随振动强迫响应计算Duhamel积分：系统分析经过充分长时间后，作为瞬态响应的前两种振动都将消失，只剩稳态强迫振动 0)(txt0 x强迫响应强迫响应全响应全响应)sin(sin)cossin(cossin)(00tBtsteBtxdddt0)0(x0)0(x 对于零初始条件系统分析3.3 多自由度无阻尼线性系统分析多自由度无阻尼线性系统分析作用力方程：作用力方程：)(tPKXXM nRXnnRKM、nRt )(P固有振

23、动方程：固有振动方程：（自由振动方程）（自由振动方程）0KXXM 在考虑系统的固有振动时，最感兴趣的是系统的在考虑系统的固有振动时，最感兴趣的是系统的同步振动同步振动，即系统在各个坐标上除了运动幅值不相同外，随时间变化的即系统在各个坐标上除了运动幅值不相同外，随时间变化的规律都相同的运动规律都相同的运动 )(tfX1)(Rtf 运动规律的时间函数运动规律的时间函数 常数列向量常数列向量 Tn21 Tnxxx21 X多自由度系统的固有频率系统分析0KXXM )(tf X代入，并左乘代入，并左乘 ：T0KM )()(tftfTT MKTTtftf)()( ：常数：常数M 正定，正定，K 正定或半正

24、定正定或半正定 对于非零列向量对于非零列向量 ： 0MT0KT20令：令：对于半正定系统，有对于半正定系统，有 0对于正定系统必有对于正定系统必有 02nRXnR系统分析2)()( MKTTtftf 0)()(2tftf 0 ,)(0),sin()(battftatfa、b、 为常数为常数0KXXM )(tf X（1）正定系统）正定系统 只可能出现形如只可能出现形如 的同步运动的同步运动)sin( taX系统在各个坐标上都是按相同频率及初相位作简谐振动系统在各个坐标上都是按相同频率及初相位作简谐振动 （2）半正定系统）半正定系统 可能出现形如可能出现形如 的同步运动的同步运动)sin( taX

25、也可能出现形如也可能出现形如 的同步运动的同步运动)(bat X（不发生弹性变形（不发生弹性变形 ）主振动主振动00系统分析首先讨论正定系统的主振动首先讨论正定系统的主振动 M 正定，正定，K 正定正定0主振动：主振动：)sin(taX正定系统：正定系统：0KXXM nRX将常数将常数 a 并入并入 中中)sin(tXTn21 代入振动方程：代入振动方程： 0)(2MK有非零解的充分必要条件：有非零解的充分必要条件：0MK2特征方程特征方程 0222212122222222212211211221211211nnnnnnnnnnnnmkmkmkmkmkmkmkmkmk系统分析022221212

26、2222222212211211221211211nnnnnnnnnnnnmkmkmkmkmkmkmkmkmk021)1(212nnnnaaa解出解出 n 个值，按升序排列为：个值，按升序排列为： 222210ni：第第 i 阶固有频率阶固有频率频率方程频率方程或特征多项式或特征多项式仅取决于系统本身的刚度、质量等物理参数仅取决于系统本身的刚度、质量等物理参数 1：基频：基频系统分析多自由度系统的模态（主振型）正定系统0主振动)sin( taX0KXXM nRXnnRKM、nR0MK)(2特征值问题特征值特征向量 n 自由度系统（固有频率）（模态）i)(i一一对应ni11)()(1)(nini

27、iR0MK)(2)(ii)(ii、代入系统分析0MK)(2)(iiTinii)()(1)( 当 不是特征多项式的重根时，上式 n 个方程中有且只有一个是不独立的 i设最后一个方程不独立，把它划去，并且把含有 的某个元素（例如 ）的项全部移到等号右端 )(i)(in)(, 12, 1)(11, 121, 1)(11 , 121 , 1)(121)(11, 121, 1)(111211)()()( )()()(innninninnninninininnininniniimkmkmkmkmkmk若这个方程组左端的系数行列式不为零，则可解出用 表示的 )(in)(1)(2)(1inii，)(i否则应把

28、含 的另一个元素的项移到等号右端，再解方程组 系统分析)(iia将 ， 代入主振动方程ii并将改为第 i 阶主振动 )sin()()(iiiiitaXTiniixx)()(1)(XTinii)()(1)( 系统在各个坐标上都将以第 i 阶固有频率 做简谐振动，并且同时通过静平衡位置 i)()()(2)(2)(1)(1ininiiiixxx 比值第 i 阶特征向量 中的一列元素，就是系统做第 i 阶主振动时各个坐标上位移（或振幅）的相对比值 )(i系统分析虽然各坐标上振幅的精确值并没有确定，但是所表现的系统振动形态已确定 描述了系统做第 i 阶主振动时具有的振动形态，称为第 i 阶主振型，或第

29、i 阶模态)(i主振动仅取决于系统的 M 阵、K 阵等物理参数。这一重要概念是单自由度系统所没有的 系统的固有振动niiiiinnnntatataat1)()(222)2(111)1()sin( )sin()sin()sin()(Xn个主振动的叠加 模态叠加法模态叠加法由于各个主振动的固有频率不相同，多自由度系统的固有振动一般不是简谐振动，甚至不是周期振动 ：)1(,niaii初始条件决定初始条件决定系统分析模态的正交性，主质量和主刚度)(ii)( jj两式相减)(2)()(2)(jjjiiiMKMK)()(2)()(jTiijTiMK)( j转置右乘转置右乘Ti)(左乘左乘)()(2)()(

30、jTijjTiMK0)()()(22 jTijiMji ji若若 时，时， 0)()( jTiM0)()( jTiK模态关于质量的正交性模态关于刚度的正交性均满足当当 ij 时时piiTim )()(MpiiTik )()(K第 i 阶模态主质量第 i 阶模态主刚度)(i第第 i 阶主模态阶主模态恒成立系统分析 00)()()()(jTijTiKM模态关于质量的正交性模态关于刚度的正交性当 ij 时piiTim )()(M主质量piiTik )()(K主刚度ji 当 时利用 Kronecker 符号 piijjTipiijjTikm)()()()(KMjijiij01第 i 阶固有频率)1(n

31、imkpipii 系统分析另一种模态：正则模态定义：全部主质量皆为1的主模态 )(iN1)()( iNTiNpimMni1)()(iiiNc令12)()(2)()( piiiTiiiNTiNmccMM)()(1ipiiNm 正则模态和主模态之间的关系：)(iN相对于 的主刚度2)()()()(1ipipiiTipiiNTiNmkm KKpiimc1 正则模态的正交性条件：正则模态的正交性条件：2)()()()(iijjNTiNijjNTiNKMjijiij01 piijjTipiijjTikm)()()()(KM主模态的正交性条件：主模态的正交性条件：系统分析)()(iTipimM )()(i

32、TipikK 第 i 阶模态主质量第 i 阶模态主刚度)(i第 i 阶主模态主模态：主模态：ni11)()(iNTiNM2)()(iiNTiNK主质量为1固有频率的平方)(iN第 i 阶正则模态正则模态：正则模态：ni1)1()(nii将 组成矩阵模态矩阵)()1(n nnR ppnpTppnpTkkdiagmmdiagKKMM),(),(11主质量矩阵主刚度矩阵系统分析)()1 ()()1 (nTnTMM ppnpTmmdiagMM ),(1推导：)()1()()1(nTnTM )()1()()1(nTnT MM )()()1()()()1()1()1(nTnTnnTTMMMM pnpmm0

33、01系统分析正则模态)1()(niiN将 组成矩阵正则模态矩阵)()1(nNNN nnRKIMNTNNTN单位矩阵谱矩阵ni1正交性条件2)()()()(iijjNTiNijjNTiNKM221n系统分析)(2)(iiiMK特征值问题：依次取 ，得到的 n 个方程，可合写为：ni1 piijjTipiijjTikm)()()()(KMni1主模态正交性条件：MK 左乘 TppMK ppKM1系统分析求解无阻尼系统对初始条件的响应求解无阻尼系统对初始条件的响应采用正则模态坐标自由振动方程 00)0()0(XXXX0KXXM ，nR XnnRKM、Tnxxx)0(,),0(),0(210 XTnx

34、xx)0(,),0(),0(210 XIM NTNK NTNNNXX 坐标变换NX：正则模态坐标N ：正则模态矩阵0XKXMNTNNTN 代入，并左乘 TN0XXINN 模态坐标初始条件01)0(XXNN01)0(XXNNTNnNNNxxx)0(,),0(),0(210 XTNnNNNxxx)0(,),0(),0(210 X系统分析NNXX 坐标变换0101)0()0(XXXX0XXI NNNNNN，)1(, 02nixxNiiNi )1(,sin)0(cos)0(nitxtxxiiNiiNiNi 在求得 后，可利用 式求得原系统的解 )1(nixNi NNXX TNnNNNxxx)0(,),

35、0(),0(210 XTNnNNNxxx)0(,),0(),0(210 X系统分析任意外部激励的响应任意外部激励的响应 n 自由度系统 )(0tFKXXM )0()0(XX，nRXNNXX )(tNFXXINN TNnNNTNNFFF,210FFnitFxxNiNiiNi1),(2 )0(01XXNN）（正则坐标初始条件)0(01XXNN）（ tiNiiiiNiiNiNidtFtxtxx0)(sin)(1sin)0(cos)0(：正则模态矩阵 N在得到 后，利用 得出原系统的解NXNNXX )1(ni 模态广义力系统分析小结小结物理空间 tMXKXF耦合主模态空间pXX pXX iTpipip

36、ipim xk xtF解耦物理空间 tMXKXF耦合正则模态空间NNXX NNXX 2 iTNiiNiNxxtF解耦系统分析3.4 稳定性浅说稳定性浅说问题：如何分析平衡位置的稳定性？系统分析稳定性的萌芽思想 2000年前 ，汉朝的淮南王刘安 淮南子说山训 ：“下轻上重，其覆必易”； 宋朝沈括在 梦溪笔谈中把这种观察到的事实付诸于应用 ，他在忘怀录 中指出：“安车车轮不欲高，高则摇” ；类似稳定，至少可以追溯1500年前到晋书上所述“行人安稳，布帆无恙” ；西方“stable”源出于拉丁文“stabilis” ，表示坚持、保持的意思；以上说法与观念表现了对稳定这一概念的最初理解。系统分析稳定性

37、科学概念的发展 18世纪下半叶到19世纪末 ，发生了一些具有深远影响的事件，从中人们可以看到稳定性理论产生的必然性。 J. Watt 1765改进了T. Newcomen 发明的蒸气机 ，引发了工业革命； J. L. Lagrange 1780年出版 分析力学，科学地讨论了平衡位置的稳定性； C. Hermite 1856年建立了关于多项式对根交错的理论；J. C. Maxwell 1868年发表的“论调节器” ，讨论了蒸气机自动调速器与时钟机构的运动稳定性；系统分析A.L. Cauchy 在19世纪给出了关于极限描述的，N语言；H. Poincare在微分方程定义的积分曲线和天体力学方面作出

38、了贡献；G. Peano，I. Bendixson和G. Darboux微分方程解对初值及参数连续依赖性的研究。 上述这些重要事件及相关科学的进展促成了19世纪末稳定性理论的两个主要学派的形成。稳定性科学概念的发展 系统分析Routh-Hurwitz （1875，1895）通过判断系统的特征根是否在左半平面判定系统是否稳定；A.M. Lyapunov 1892发表著名的博士论文运动稳定性一般问题，通过考察系统能量是否衰减来判定稳定性。两个主要学派 系统分析 平衡的稳定性(stability of equilibrium)：质点系处于某一平衡位置，若受到微小干扰偏离平衡位置后总不超出平衡位置邻近

39、的某个微小区域，则称质点系在该位置的平衡是稳定的(stable)，否则是不稳定的(unstable)。) *, *(:),(:11nnqqqq其平衡位置为设系统的广义坐标为成立使得若201021)(,*)(:0,0, 0tqqtqiii，否则是不稳定的则称平衡位置是稳定的成立时，有当,*)(0iiqtqtt系统分析质点系在势力场中平衡的充分必要条件：), 1( , 0kjqVj ABC定理：质点系在势力场中的平衡位置是稳定的充分必要条件是系统在平衡位置的势能为极小值。系统分析平衡是如何实现的？系统分析3.5 混沌浅说混沌浅说混沌是非线性系统特有的一种运动形式，产生确定性系统的敏感依赖于初始条件

40、的往复性非周期运动，类似于随机振动而具有长期不可预测性。Edward N. Lorenz系统分析例例2 2：系统如图示，设均质杆系统如图示，设均质杆OA长为长为 L 质量为质量为m，其阻力系数为，其阻力系数为 c，扭，扭簧刚簧刚度系数为度系数为 k。当。当 OA BC 时，时，扭扭簧无簧无变形，已知变形，已知BC杆的转动：杆的转动： , , 试建立试建立OA杆的运动微分方程。杆的运动微分方程。tbcos OABCgm系统分析解：研究OA杆，进行受力分析kcmgOMMMJ cMc)( kMk niiOOMtJ122)(ddFOAgmxFyFkMcMsin2LmgMmg tkbmgLkcmLcos

41、sin21312 常系数非齐次非线性微分方程系统分析OABCgmtkbmgLkcmLcossin21312 rad/s,0 . 0rad,0 . 0d0.65Nms/ra0.07rad,Nm/rad,0 .14m,1kg,300cbkLm在一定参数条件下，响应周期与激励周期相同在一定参数条件下，响应周期与激励周期相同系统分析系统分析rad/s,0 . 0rad,0 . 0d0.65Nms/ra0.22rad,Nm/rad,0 .14m,1kg,300cbkLm在一定参数条件下，响应周期与激励周期不同在一定参数条件下，响应周期与激励周期不同系统分析系统分析s/trad/rad/s,0 . 0ra

42、d,001. 0:(2)rad/s,0 . 0rad,0 . 0: ) 1 (d0.64Nms/ra0.02rad,Nm/rad,0 .14m,1kg,30000cbkLm系统分析系统分析 xg1mg2mAB0lk0102030405060708090 100-0.500.511.522.5s/ t050100150200250300-10123456m/)(txrad/)(t摆杆的运动摆杆的运动滑块的运动滑块的运动例例3：例：例1的数值仿真结果的数值仿真结果212222222()cossin01(2 )cossin03mmxm Lm LkxmLm Lxm gL系统分析00.10.13.2 0

43、.4m/1xrad/2x初值在初值在红色红色区域内系区域内系统将产生统将产生 混沌运动混沌运动2x1xg1mg2mAB0l0k1k系统控制4.1 概述概述控制是调节系统状态变化定性或定量特征的过程。控制理论是门内涵丰富的技术科学。它研究如何按被控对象和环境特性、通过能动地采集和运用信息、施加控制作用而使系统正常运行且具有预定的功能。 控制理论的主要方法包括状态空间方法和频域方法。 状态空间方法是通过一组适当选取的状态变量，将系统中的动态过程及静态过程用常微分方程组及代数方程组的形式表示。这里系统的方程可反映系统的每个状态变量的变化过程，故是种内部的描述。在现代控制理论中，常采用状态空间方法。本

44、节仅介绍状态空间方法。 系统控制频域方法是通过传递函数描述系统或系统的各部分的输入和输出的关系。由于它没有直接涉及系统内部的众多变量，而仅涉及最主要的变量：输入和输出，故它是种外部描述。传统上频域方法主要适用于线性系统，近来也被发展到非线性系统，主要是通过广义频率响应函数或非线性传递函数的分析确定非线性系统某些频率特性。频域方法物理意义明确，便于实验，所设计的控制器具有鲁棒性。因此，经典控制理论中，频域方法是最常用的，并取得良好效果。 系统控制4.2 控制系统控制系统一般的控制系统由受控系统、测量装置和执行机构组成。 控制系统数学描述包括状态方程和输出方程，状态方程表示系统由于输入所引起内部状

45、态的变化，输出方程表示系统输出与系统状态和系统输入之间的关系。状态方程与输出方程的组合称为状态空间表达式，又称为动态方程。 , , ,.ttxfx uyh x u其中为xRn状态变量，yRl为输出变量，uRm为控制信号，又称为系统输入，实数t为时间。 控制系统分为定常系统和非定常系统，前者不显含时间t，后者显含时间t。控制系统还可以分为线性系统和非线性系统，若状态方程和输出方程的右端均为状态、输入和时间的线性函数，则为线性系统，否则为非线性系统。 系统控制4.3 局部线性化与线性系统局部线性化与线性系统局部线性化是种近似的线性化方法。该方法基于将非线性函数在工作状态附近展开成Taylor级数并

46、仅保留线性项。因为忽略了Taylor级数展开式中的高阶项，故这些被忽略的高阶项必须很小，即状态变量仅能对工作状态有微小偏离。局部线性化是非线性系统的一种最简单的简化方式，适用于许多非线性系统。该方法存在误差，而且误差随偏离理想状态的增大而急剧增大，故该线性化方法仅在工作状态的小邻域内研究控制精度要求不是很高的问题才有效。当工作状态的可能变化范围很大或需要控制精度很高时，局部线性化方法不适用。为保证局部线性化的适用性，有必要对基于该方法得到的结论如控制律代入原来的非线性模型进行理论分析或数值仿真。 系统控制0000ttttttxuxuxfxxf ufyhxxh uhDDDDDD其中各Jacobi

47、矩阵 11,iiitjknn mn nppptjkl mll nfffxuthhhxutxuxufffhhhDDDDDD系统控制例例4：倒立摆的控制：倒立摆的控制单摆倒立在可水平移动的小车上。小车和单摆均在同一平面内运动。车的质量为M，摆球的质量为m，摆杆的质量不计，长度为l。控制力u作用在小车上。摆在有干扰时保持摆垂直，且在每一控制过程结束时小车都返回初始位置。导出受控系统数学模型，进行局部线性化。 小车和摆球均可模型化为质点。用小车离开初始位置的距离x和摆杆偏离垂直位置的角度描述系统的运动。系统的Lagrange函数 22211cossincos22LMm xxllmgl系统控制应用Lag

48、range方程可导出系统运动方程 2cossincossin0Mm xmlmlumlmlxmgl引入状态变量1234,xxxxxx受控系统的数学模型 122112113421141sinsinsin1cossinsin1cosxxmgxxxugxxlMmxlxxmgxxxuxMmx系统控制分别在x=0计算式右端关于x和u的Jacobi矩阵，导出受控系统数学模型的局部线性化 11223344010001000000101000 xxMmxxMlMluxxxxmgMM系统控制4.4 控制设计控制设计设计问题是确定系统的反馈控制律使系统具有稳定性等性能指标，又称为综合问题。设计问题的一个特殊情形是参

49、数选择问题。假设控制律具有某种给定的形式，确定其中的待定参数。 控制系统的任务可以分为两类：镇定(又称调节)和跟踪(又称伺服)。镇定问题是要计控制系统使得受控系统的状态稳定在一个平衡点的周围。跟踪问题是要设计控制系统使得受控系统的输出跟踪一个给定的时间信号。严格而论，镇定问题是指渐近镇定问题。对于满足关系 , ,ttfh0 0000求控制律u使从任意允许区域开始的状态有输出变量满足 limtty= 0系统控制跟踪问题也是指渐近跟踪问题。对于系统和任意给定的期望输出yd(t)，求控制律u，使得从允许区域中任意状态开始，跟踪误差e(t)=y(t)-yd(t)满足 dlimtttyy= 0求解跟踪问

50、题一般比求解镇定问题困难，因为在跟踪问题中控制器不仅要保持整个状态的稳定而且还要驱使系统的输出逼近期望的输出。在理论上，两者是互相联系的。跟踪问题等价于以跟踪误差为状态变量的镇定问题，而镇定问题是以系统平衡点为期望输出的跟踪问题。 tuu具体实现控制有两种基本方案：开环和闭环。在开环控制律中，控制信号仅是时间的函数 系统控制在闭环控制中，控制信号还与状态变量和输出变量有关。开环控制又称为前馈控制，闭环控制又称反馈控制。根据控制变量的不同，反馈控制又可分为输出反馈和状态反馈。在输出反馈中，控制信号是输出变量和时间的函数 ,tuu y在状态反馈中，控制信号是状态变量和时间的函数 ,tuu x所引入