控制系统的数学模型 (2).ppt

《控制系统的数学模型 (2).ppt》由会员分享，可在线阅读，更多相关《控制系统的数学模型 (2).ppt（47页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、上海大学机电工程与自动化学院工 程 控 制 原 理2.数学模型与传递函数主讲：彭艳 办 公 室：机械楼205室 电子邮件：办公电话：56334137上海大学机电工程与自动化学院 建立控制系统的数学模型，并在此基础上对控制系统进行分析、综合，是机电控制工程的基本方法。如果将物理系统在信号传递过程中的动态特性用数学表达式描述出来，就得到了组成物理系统的数学模型。经典控制理论 经典控制理论采用的数学模型主要以传递函数 传递函数为基础。而现代控制理论 现代控制理论采用的数学模型主要以状态空间方程 状态空间方程为基础。而以物理定律及实验规律为依据的微分方程又是最基本的数学模型，是列写传递函数和状态空间方

2、程的基础。2.数学模型与传递函数上海大学机电工程与自动化学院2.1 控制系统的数学模型2.1.1 数学模型 对于一个复杂的物理系统，为了对系统的动态特性进行分析和综合，必须用数学表达式来描述该系统，这个表达式称为该系统的“数学模型”。由于动态过程中有关物理量都是时间的函数，所以，通常用微分方程 微分方程来描述系统。2.数学模型与传递函数数学模型 数学模型描 描述 述系 系统 统或 或元 元件 件的 的动 动态 态特 特性 性的 的数 数学 学表达式 表达式 深 深入 入了 了解 解元 元件 件及 及系 系统 统的 的动 动态 态特 特性 性，准确建立它们的数学模型 准确建立它们的数学模型 建

3、建 模 模上海大学机电工程与自动化学院2.1.1 数学模型 数学模型是描述物理系统的运动规律、特性和输入输出关系的一个或一组方程式。系统的数学模型可分为静态和动态数学模型。静 静态 态数 数学 学模 模型 型：反映系统处于平衡点（稳态）时，系统状态有关属性变量之间关系的数学模型。即只考虑同一时刻实际系统各物理量之间的数学关系，不管各变量随时间的演化，输出信号与过去的工作状态（历史）无关。因此，静 静态 态模 模型 型都 都是 是代数式，数学表达式中不含有时间变量。代数式，数学表达式中不含有时间变量。2.1 控制系统的数学模型上海大学机电工程与自动化学院2.1.1 数学模型 动态数学模型 动态数

4、学模型：描述动态系统瞬态与过渡态特性的模型。也可定义为描述实际系统各物理量随时间演化的数学表达式。动态系统的输出信号不仅取决于同时刻的激励信号，而且与它过去的工作状态有关。微分方程或差分方程常用作动态数学模型。对于给定的动态系统，数学模型不是唯一的。对于给定的动态系统，数学模型不是唯一的。工程上常用的数学模型包括：微分方程，传递函数和状态方程。对于线性系统，它们之间是等价的。针对具体问题，选择不同的数学模型。2.1 控制系统的数学模型上海大学机电工程与自动化学院2.1.2 建立控制系统数学模型的方法 物理系统往往比较复杂，因而必须作一些理想化的假设，获得简化的数学模型。理论分析法（解析法）：对

5、系统各部分的运动机理进行分析，依据系统本身所遵循的有关定律列写数学表达式，并在列写过程中进行必要的简化。实验研究法 根据系统对某些典型输入信号的响应或其它实验数据建立数学模型。即人为施加某种测试信号，记录基本输出响应。对于比较复杂的系统 复杂的系统，需要通过理论分析与实验研究结合起来，才能获得适用的数学模型。2.1 控制系统的数学模型上海大学机电工程与自动化学院2.1.3 线性系统与非线性系统(1)线性系统 若描述系统的微分方程是变量及其导数的一次有理整式，则此系统称为线性系统。线性系统可以用线性微分方程描述。如果方程的系数为常数，则为线性定常系统；如果方程的系数是时间 t 的函数，则为线性时

6、变系统。线性系统的线性性质是指系统满足叠加性 叠加性和齐次性 齐次性。v v叠 叠加 加性 性：指当几个激励信号同时作用于系统时，总的输出响应等于每个激励单独作用所产生的响应之和。v v齐 齐次 次性 性：指当输入信号乘以某常数时，响应也倍乘相同的常数。2.1 控制系统的数学模型上海大学机电工程与自动化学院(2)非线性系统 不是线性系统的系统称为非线性系统。非线性系统用非线性微分方程描述。非线性系统不满足叠加性 叠加性和齐次性 齐次性。系统中只要含有一个非线性性质的元件，就成为一个非线性系统。许多机械系统各物理量之间的关系都是非线性的，即使对所谓的线性系统来说，也只是在一定的工作范围内保持线性

7、关系。因此，研究机械系统的某些动态特性时，必须考虑系统中的非线性特征。2.1.3 线性系统与非线性系统上海大学机电工程与自动化学院 例其中：a，b，c，d 均为常数 线性定常系统线性时变系统非线性系统2.1.3 线性系统与非线性系统上海大学机电工程与自动化学院(2)非线性系统 许多机械系统、电气系统、液压及气动系统等，其物理量之间都包含有非线性关系。例如：在大输入信号作用下，元件的输出量可能饱和(即饱和非线性)；在小信号输入下，元件没有输出量(即死区非线性)；某些元件中可能存在着平方非线性关系。如下图所示。2.1.3 线性系统与非线性系统饱和非线性 死区非线性 平方律非线性各种非线性因素的特性

8、曲线 各种非线性因素的特性曲线上海大学机电工程与自动化学院(2)非线性系统 是否线性元件的判别 是否线性元件的判别：元件的输出与输入间关系 输出与输入间关系为 一次幂函数 线性元件二次或高次幂函数、周期函数或超越函数 非线性元件 判别系统数学模型微分方程是否非线性 判别系统数学模型微分方程是否非线性：其中的函数出现高于一次的项，或者函数导数项的系数是输出量的函数。2.1.3 线性系统与非线性系统上海大学机电工程与自动化学院 机械系统中常见的一些非线性问题 机械系统中常见的一些非线性问题 2.1.3 线性系统与非线性系统传动间隙 传动间隙：由齿轮及丝杠螺母副组成的机床进给传动系统中，经常存在有传

9、动间隙，使输入转角 x xi i 和输出位移 x xo o 之间有滞环关系。传动间隙oxixo死区xidxodto死区 死区：在死区范围内，系统有输入而没有输出。例如：负开口的液压伺服阀 负开口的液压伺服阀上海大学机电工程与自动化学院2.1.3 线性系统与非线性系统摩擦力 摩擦力：机械滑动运动副，如：机床滑动导轨运动副、主轴套筒运动副、活塞液压缸运动副等，在运动中都存在摩擦力。右图为干 干摩 摩擦 擦力 力（也称库伦摩擦力，其大小为 f f）与运动速度 d dx/x/d dt t 的关系。摩擦力 f f 总是与速度 d dx/x/d dt t 的方向相反。干摩擦力foFdxdt上海大学机电工程

10、与自动化学院2.1.3 线性系统与非线性系统摩擦力 摩擦力：实际运动副中的摩擦力与运动速度的大小有关，如右图。F F与d dx/x/d dt t 的曲线大致分为三个阶段：（F F与d dx/x/d dt t 的方向相反）粘性摩擦力dxdtoFq起动时的静动摩擦力；（摩擦力数值较大）q低速时的混合摩擦力；（摩擦力呈下降特性）q粘性摩擦力。（摩擦力随速度的增加而增加）上海大学机电工程与自动化学院2.1.4 系统非线性微分方程的线性化 工程上 工程上，绝对的线性元件和线性系统是不存在的；数学上 数学上，非线性微分方程的求解困难，目前只能采用数值解法，但也存在较大的数值误差。为了解决工程实际问题，必须

11、对非线性微分方程进行线性化处理。本质非线性性质：在工作点附近存在着不连续直线、跳跃、折线，以及非单值关系等严重非线性性质的元件或系统。2.1 控制系统的数学模型上海大学机电工程与自动化学院 对于非 非本 本质 质非线性元件或系统，可以在工作点附近用切线来替代函数关系，这就是非线性数学模型的线性化方法之一（微小偏差法）。系统正常工作时，通常都有一个预定工作点，即系统处于某一平衡位置，对于自动调节系统或随动系统，只要系统的工作状态稍一偏离此平衡位置，整个系统就会立即作出反应，并力图恢复原来的平衡位置。2.1.4 系统非线性微分方程的线性化上海大学机电工程与自动化学院 具有连续变化的非线性函数的线性

12、化，可用切线法或微小偏差法。在预定工作点处一个小范围内，将非线性特性用一断直线来代替。一个变量的非线性函数f f(x x)，在 x x0 0处连续可微，则可将它在该点附件用Taylor级数展开2.1.4 系统非线性微分方程的线性化of(x)xf(x0)x0 xf(x)增量较小时，可略去其二次以上高次幂项，则有上海大学机电工程与自动化学院 两个变量的非线性函数 f f(x x,y y)，在(x x0 0,y y0 0)处连续可微，则也可将它在该点附件用Taylor级数展开2.1.4 系统非线性微分方程的线性化 增量较小时，可略去其二次以上高次幂项，则有上海大学机电工程与自动化学院单变量非线性函数

13、 f f(x x)的线性化数学模型为2.1.4 系统非线性微分方程的线性化双变量非线性函数 f f(x x,y y)的线性化数学模型为式中：；式中：；上海大学机电工程与自动化学院 非线性系统线性化时的注意事项 注意事项：(1)必须明确系统处于平衡状态的工作点，因为不同工作点所得到线性化方程的系数不同，即：非线性曲线上各点的斜率（导数）是不同的；(2)如果变量在较大范围内变化，则用这种线性化方法建立的数学模型，除工作点外的其他工况必定存在较大的误差。因此非线性系统线性化是有条件的：变量偏离预定工作点很小。(3)对于某些典型的本质非线性，如果非线性函数是不连续的，则在不连续点附近不能得到收敛的Ta

14、ylor级数，此时就不能线性化。(4)线性化后的微分方程是以增量为基础的增量方程。2.1.4 系统非线性微分方程的线性化上海大学机电工程与自动化学院 试把非线性方程 z z xy xy 在区域 5 5x x7 7、10 10y y12 12 上线性化。求用线性化方程来计算当 x x=5=5，y y=10=10 时 z z 值所产生的误差。解 解：由于研究的区域为5x7、10y12，故选择工作点x0=6，y0=11。于是 z0=x0y0=611=66。求在点 x0=6，y0=11，z0=66 附近非线性方程的线性化表达式。将非线性方程在点（x0,y0,z0）处展开成泰勒级数，并忽略其高阶项，则有

15、例 题上海大学机电工程与自动化学院因此，线性化方程式为：例 题 当 x x=5=5，y y=10=10 时，z z 的精确值为：z z xy xy=510=50由线性化方程求得的z z值为 z z=11x+6y-66=55+60-66=49 因此，误差为 50-49=1，表示成百分数为上海大学机电工程与自动化学院2.1.5 系统微分方程的建立 建立微分方程的步骤 建立微分方程的步骤是：1)分析系统的工作原理和系统中各变量间的关系，确定出待研究元件或系统的输入量和输出量；2)从输入端入手（闭环系统一般从比较环节入手），依据各元件所遵循的物理、化学、生物等规律，列写各自方程式，但要注意负载效应。所

16、谓负载效应，就是考虑后一级对前一级的影响。3)将所有方程联解，消去中间变量，得出系统输入输出的标准方程。2.1 控制系统的数学模型上海大学机电工程与自动化学院 标准方程包含三方面的内容 标准方程包含三方面的内容：将与输入量有关的各项放在方程的右边；与输出量有关的各项放在方程的左边；各导数项按降幂排列。2.1.5 系统微分方程的建立上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：机械系统中部件的运动有直线和转动两种。机械系统中以各种形式出现的物理现象，都可简化为质量、弹簧和阻尼三个要素。列写其微分方程通常用达朗贝尔原理。即：作用于每一个质点上的合力，同质点惯性力

17、形成平衡力系。质量：质量：2.1.5 系统微分方程的建立mu(t)x(t)fm(t)上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：弹簧：弹簧：弹簧各点受力相同，但各点的变形量不同。2.1.5 系统微分方程的建立ku1(t)x1(t)fk(t)u2(t)x2(t)fk(t)上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：阻尼：阻尼：2.1.5 系统微分方程的建立Du1(t)x1(t)fD(t)u2(t)x2(t)fD(t)上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：2.1.5 系统微分方程

18、的建立 例1直线运动（机械平移系统）直线运动（机械平移系统）m0 xo(t)fi(t)k Dm0 xo(t)fi(t)fm(t)fk(t)fD(t)机械平衡系统及其力学模型静止（平衡）工作点作为零点，以消除重力的影响上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：2.1.5 系统微分方程的建立 例1直线运动（机械平移系统）直线运动（机械平移系统）式中：m、D、k 通常均为常数，故机械平移运动系统可以由上述二阶常系数微分方程描述。上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：2.1.5 系统微分方程的建立 例2弹簧弹簧-阻尼系

19、统阻尼系统 系统运动方程为一阶常系数微分方程。0 xo(t)fi(t)k D弹簧-阻尼系统上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：2.1.5 系统微分方程的建立 例3机械旋转系统机械旋转系统 图中：图中：J 旋转体转动惯量；k 扭转刚度系数；D 粘性阻尼系数。0i(t)kD齿轮o(t)0JTD(t)Tk(t)柔性轴粘性液体上海大学机电工程与自动化学院(1)机械系统微分方程的列写 机械系统微分方程的列写：2.1.5 系统微分方程的建立 例3机械旋转系统机械旋转系统 数学模型：扭矩平衡方程上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系

20、统微分方程的列写：电网络系统分析主要根据基尔霍夫电流定律和电压定律写出微分方程式，进而建立系统的数学模型。a)a)基 基尔 尔霍 霍夫 夫电 电流 流定 定律 律：汇聚到某节点的所有电流之代数和应等于0（即流出节点的电流之和等于所有流进节点的电流之和）。b)b)基 基尔 尔霍 霍夫 夫电 电压 压定 定律 律：电网络的闭合回路中电势的代数和等于沿回路的电压降的代数和。电网络系统三个基本元件：电阻、电容和电感 电阻、电容和电感。2.1.5 系统微分方程的建立上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系统微分方程的列写：电阻：电阻：2.1.5 系统微分方程的建立Ru(t)i

21、(t)电压：电容：电容：电压：Cu(t)i(t)电感：电感：电压：Lu(t)i(t)上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例1RR-LL-CC 无源电路网络无源电路网络 LRCi(t)ui(t)uo(t)基本物理规律：上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例1RR-LL-CC 无源电路网络无源电路网络 数学模型：一般情况下，R、L、C均为常数，上式为二阶常系数微分方程。若L0，则系统微分方程简化为上海大学机电工程与自动化学院(2)电网络

22、系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例2有源电路网络有源电路网络 aRCi1(t)ui(t)uo(t)i2(t)-+即：上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例3直流电动机的伺服控制系统直流电动机的伺服控制系统 RaLaT(t)o(t)JDia(t)ei(t)em(t)if=常数激磁电流 激磁电流电机电枢 电机电枢输入电压 输入电压电机输出 电机输出转角 转角电机感应 电机感应电动势 电动势电机绕组电感 电机绕组电感 绕组电阻 绕组电阻电机扭矩上海大学机电工程与自动化学院(

23、2)电网络系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例3数学模型：基尔霍夫定律电磁感应定律扭矩平衡方程直流电动机的伺服控制系统直流电动机的伺服控制系统 上海大学机电工程与自动化学院(2)电网络系统微分方程的列写 电网络系统微分方程的列写：2.1.5 系统微分方程的建立 例3电枢控制式直流电动机的控制系统动态数学模型 当电枢电感很小时，通常可以忽略不计，则系统微分方程可简化为直流电动机的伺服控制系统直流电动机的伺服控制系统 上海大学机电工程与自动化学院流入 流入流出 流出例题：右图为四边伺服阀及液压缸组成的液压伺服控制系统。2.1 控制系统的数学模型cmmxi

24、xoxvp1p2液压伺服系统 活塞杆固定，液压缸与伺服阀阀体连为一体，工作时缸体带动工作台(质量为m)移动，负载包含有惯性和粘性负载。试建立该系统的数学模型。上海大学机电工程与自动化学院解：给阀杆一位移量 xi，此时阀口打开，液压油推动缸体并带动负载产生位移 xo，阀体有此位移相当于给阀口开度一个反馈量 xo，即 xv=xi-xo。该系统是输入为 xi、输出为 xo 的闭环系统。假设油液不可压缩，并忽略系统中液体的泄漏。系统动力学方程 系统动力学方程：2.1 控制系统的数学模型 例题例题质量守恒定律(连续性方程)p=p1-p2：负载压降粘性阻尼系数负载流量牛顿第二定律(力平衡方程)活塞有效面积

25、负载质量上海大学机电工程与自动化学院流体流经微小开口的伺服阀时，其流量关系 流体流经微小开口的伺服阀时，其流量关系：2.1 控制系统的数学模型 例题例题是一非线性函数关系。由上述三个关系式得出的系统数学模型也是非线性方程。因此必须对式(a)进行线性化处理。(a)写成增量方程-Kc：流量-压力系数Kq：流量增益 将式(a)在预定工作点(x xv0 v0,p p0 0)附近进行微小偏差线性化上海大学机电工程与自动化学院增量方程：2.1 控制系统的数学模型 例题例题(b)K Kq q、K Kc c可以根据工作点的不同，从阀的特性曲线求得。（不同的工作点，有不同的K Kq q、K Kc c）由于随着负

26、载压降 p p 的增大，负载流量 q q 将减小，因而q q/p p总是负值。为使流量-压力系数 K Kc c 为正，定义 若预定工作点选在阀的零位，即：x xv0 v0=0=0，p p0 0=0=0，q q(x xv0 v0,p p0 0)=0=0；这时 q q=q q，x xv v=x xv v，p p=p p。上海大学机电工程与自动化学院此时增量方程式(b)变成：2.1 控制系统的数学模型 例题例题此乃预定工作点为伺服阀零位时的线性化方程。(c)pqaaaaxv=2xvxv=xvxv=0 xv=-xvxv=-2xvo线性关系 右图为q q与p p的线性化关系曲线，q q与 x xv v之间也成线性化关系。上海大学机电工程与自动化学院 将式(c)、质量守恒方程代入力平衡方程，消去中间变量q q与p p后，得到系统的阀口开度x xv v与位移输出x xo o关系的线性化微分方程：2.1 控制系统的数学模型 例题例题将 xv=xi-xo 代入上式，得到闭环系统的线性微分方程 闭环系统的线性微分方程 若考虑油液的压 压缩 缩性 性，上式将变成三阶方程；如再考虑油液泄漏 泄漏，方程将更为精确。上海大学机电工程与自动化学院作业：P.362-12.1 控制系统的数学模型