控制数学模型上课讲义

控制数学模型精品文档第二章控制系统的数学模型2—1数字模型在控制系统的分析和设计中，首先要建立系统的数学模型。电气的、自动控制系统：相同的数学模型进行描述，研究自动控机械的、制系统其内在共性运动规律。系统的数学模型，是描述系统内部各物理量之间动态关系的数学达式。微(差)分方程常用的数学模型有：传递函数(脉冲传递函数研究线性离散系统的数学模型)经典控制理论数学模型解析法：依据系统及元件各变量之间所遵循的物理、化学定的建立方法律，一般应尽可能采用线性定常数学模型描述控制系统。如果描述系统的数学模型是线性微分方程，则称该系统为线性系统，若方程中的系数是常数，则称其为线性定常系统。线性系统的最重要特性是可以应用叠加原理，在动态研究中，如果系统在多个输入作用下的输出等于各输入单独作用下的输出和(可加性)，而且当输入增大倍数时，输出相应增大同样倍数(均匀性)，就满足叠加原理，因而系统可以看成线性系统。如果描述系统的数学模收集于网络，如有侵权请联系管理员删除精品文档型是非线性微分方程，则相应系统称为非线性系统，其特性是不能应用叠加原理。建立系统数学模型的主要目的，是为了分析系统的性能。由数学模型求取系统性能指标的主要途径如图2—1所示。由图可见，傅里叶变换和拉普拉斯变换是分析和设计线性定常连续控制系统的主要数学工具。求解时间响应观察线性常微分方程性能指标傅拉氏变换拉氏反变换估算估算氏传递函数变换s=jω计算频率特性频率响应图2-1求取性能指标的主要途径2－2运用微分方程建立数学模型控制系统中的输出量和输入量通常都是时间的函数。很多常见的元件或系统的输出量和输入量之间的关系都可以用一个微分方程表示，方程中含有输出量、输人量及它们各自对时间的导数或积分。这种微分方程又称为动态方程、运动方程或动力学方程。微分方程的阶数一般是指方程中最高导数项的阶数，又称为系统的阶数。建立系统微分方程的一般步骤或方法是：1)根据研究问题的需要，确定系统的输入和输出。2)对实际系统进行适当的简化，如将分布参数集中化、将非线性因素线性化等。收集于网络，如有侵权请联系管理员删除精品文档3)根据系统、输入和输出三者之间动态关系的原理或定律，列写系统的微分方程。若系统比较复杂，则需分段列写微分方程，在这种情况下，必须注意各分段之间的负载效应问题。4)消去中间变量，将方程整理成标准形式，即将与输出有关的项列在等号左边，而将与输入有关的项列在等号右边，且各阶导数按降幂排列。列写微分方程的关键是元件或系统所属学科领域的有关规律而不是数学本身。但求解微分方程需要数学工具。下面分别以电路系统和机械系统为例，说明如何列写系统或元件的微分方程式。2－2－1电路系统电路系统的基本要素是电阻、电容和电感，而建立数学模型的基本定律是基尔霍夫电流定律∑i=0,以及基尔霍夫电压定律∑u=0。元件与电压电流的关系电阻：uRi电感：uLdi1电容：uidtdtC以下举例说明电路系统方程的建立。例2—1如图2—2所示为一个RLC串联电路，试求其数学模型。解设输入信号x(t)ui(t)输出信号y(t)u0(t)。按照基尔霍夫电压定律得图2—2RLC电路uiuRuLu0,uRRi收集于网络，如有侵权请联系管理员删除精品文档uldiLdtu01idtC消去中间变量i得系统的微分方程为：LCd2u0RCdu0u0ui（２－１）dt2dt令T1＝LC，2＝RC，同时将x(t)i与y(t)0(t)代人可得Tu(t)uT1d2y(t)T2dy(t)y(t)x(t)（２－２）dt2dt这是一个典型的二阶线性常系数微分方程，对应的系统也称为二阶线性定常系统。例2—2如图2-3所示为由两个ＲＣ电路串联而成的滤波网络，试建立输入电压ui和输出电压u。之间动态关系的微分方程。解设回路电流i1，和i2为中间变量。根据基尔霍夫电压定律对前一回路、后一回路有：uiR1i11i2)dt(i1C111图2—3两个RC串联网络i2dt(i1i2)dtR2i2C1C21i2dtu0C2由上三式消去中间变量i1，和i2，整理即得ui和u0之间动态关系的微分方程：d2u0(R1C1R2C2du0u0ui（２－３）R1C1R2C22R1C2)dtdt收集于网络，如有侵权请联系管理员删除精品文档由上例明显看出，系统中后一部分对前一部分的负载效应,反映在流过前一回路电容C１的电流上，没有后一回路时为i1，而当串联上后一回路则为i1-i2。从能量的角度看，负载效应就是后一回路带走了前一回路的一部分能量。从信息传递的角度看，负载效应就是系统的两个部分之间所存在的信息的内部直接反馈作用。2－2－2机械系统机械系统指的是存在机械运动的装置，常用的基本要素是质量、弹簧和阻尼器。它们遵循物理学的力学定律。机械运动包括直线运动(相应的位移称为线位移)和转动(相应的位移称为角位移)两种。做直线运动的物体要遵循的基本力学定律是牛顿第二定律md2ydt2式中F为物体所受到的力，m为物体质量，y是线位移，t是时间。转动的物体要遵循如下的牛顿转动定律Jd22dt式中T为物体所受到的力矩，J为物体的转动惯量，θ为角位移。例2－３如图2—4所示为一个，求其数学模型。图2-4带阻尼的质量弹簧系统收集于网络，如有侵权请联系管理员删除精品文档解设输入量为x(t)F，位移输出量为y(t)s。由牛顿定律得：FFmFfFkFfdsfdtFkksFmMd2sdt2代人力平衡方程式后得FMd2sfdsks（２－dt2dt４）令T1M/f,T2f/k，K1/k并将x(t)，y(t)代入上式得该机械运动系统的数学模型：T1T2d2y(t)T2dy(t)（２－dt2y(t)Kx(t)dt５）该系统是二阶线性定常系统。例２－4图2－５所示为一机械旋转系统。转动惯量为Ｊ的圆柱体，在转矩Ｔ的作用下产生角位移θ，求该系统的输入—输出描述。图2-5机械旋转系统(a)原理图(b)分离体解假定圆柱体的质量分布均匀，质心位于旋转轴上，而且惯性主轴和旋转主轴线相重合，则其运动方程可写成：收集于网络，如有侵权请联系管理员删除精品文档2TffdJdf2TTfTkdtdtTkk式中f──粘性摩擦系数，常数ω──角速度ｋ──弹性扭转变形系数，常数就得到输入与输出关系的微分方程：Jd2fdkT（２－dt2dt６）由以上描述的数学模型可以看出，系统的数学模型由其结构、参量及基本定律决定。还有如机电系统、热工系统、化工系统，都可以通过其物理、化学机理找到其数学模型。2-2-3线性系统微分方程的通用形式在一般线性系统，描述系统动态方程的标准形式为andny(t)an1dn1y(t)a1dy(t)a0y(t)dtndtn1dtbmdmx(t)dm1x(t)dx(t)（２－７）bm1b1b0x(t)dtmdtm1dt式中：x(t)为系统输入信号；y(t)为系统输出信号；ai(i＝0，1，2，，n)、bj(j＝0，1，2，，rn)为系数，n为输出信号的最高求导次数；m为输入信号的最高求导次数。若ai和bj均为常数时，上式为常系数线性微分方程，所描述的系统为定常线性系统。收集于网络，如有侵权请联系管理员删除精品文档2－3线性系统的传递函数微分方程：时间域；微分积分求解；环节增减分析不便，阶数高求解繁难不同的初始条件，输出响应不同传递函数：复数域；代数运算求解；环节增减分析方便，阶数高求解因式分解初始条件必须为零，研究动态特性，经典控制理论最基本数学方法微分方程与传递函数：连续系统利用传递函数还可研究系统参数变化或结构变化对动态过程的影响，因而使分析系统的问题大为简化。另外，还可以把对系统性能的要求转化为对系统传递函数的要求，使综合设计的问题易于实现。2－3－1传递函数的概念传递函数是描述线性定常系统输入－输出关系的一种最常用的表达式。引入微分算子：sd,则1dt。dts系统的传递函数可以定义为：在所有初始条件均为零时，系统输出的拉氏变换与系统输入的拉氏变换之比：G(s)Y(s)。（２－８）X(s)设有一线性定常系统，其微分方程表达式为2-7式。假定初始条件均为零，．．．．．．．．．前式的拉氏变换可写为：(ansnan1sn1a1sa0)Y(s)(bmsmbm1sm1b1sb0)X(s)由此可得系统的传递函数为：G(s)Y(s)bmsmbm1sm1b1sb0X(s)ansnan1sn1a1sa0(2-9)举例说明：例２－5由例2－1的RLC电路，求其传递函数。收集于网络，如有侵权请联系管理员删除精品文档解１由式(2-2)RLC电路的微分方程：T1d2y(t)T2dy(t)y(t)x(t)dt2dt初始条件为零，对上式进行拉氏变换得：(T1s2T2s1)Y(s)X(s)∴传递函数为：G(s)Y(s)1X(s)T1s2T2s1解２在推导电网络的传递函数时，对于无源元件电感Ｌ、电容C和电阻Ｒ，分别用它们的复阻抗求解往往是比较简便的。令Z１＝Ｒ+Ｌｓ，为电阻和电感的复数阻抗之和；Z２=1／Cs为电容的复数阻抗。则：U0(s)Z21/Cs11G(s)Z1Z2RLs1/CsLCs2RCs1T1s2T2s1Ui(s)另外例课本2-10,2-11,2-12。２－３－２传递函数的性质1)传递函数的定义，只是对线性系统而言，严格地说，还只是对定常系统而言。2)传递函数通常是复变量ｓ的有理分式，其分子、分母多项式各项系数均为实数，这些系数均由系统的物理参数所确定，且nm。3)传递函数表征了系统本身的特性，它是系统动态性能的解析描述，它与输入激励无关，也与初始条件无关。4)传递函数并不是系统具体物理结构的描述，所以对于许多物理性质截然不同的系统，如机械系统、电子系统、热传导系统，都可以具有相同的传递函数。收集于网络，如有侵权请联系管理员删除精品文档5)传递函数的分母多项式：D(s)ansnan1sn1a1sa0（２－１０）就是系统的特征多项式，它的阶次，也就代表了系统的阶次。6)传递函数的拉氏反变换是系统的脉冲响应。7)传递函数可以更形象地在复数平面上描述系统的动态特性。对于实际的元件和系统，传递函数是复变量s的有理分式，其分子N(s)和分母D(s)都是s的有理多项式，即它们的各项系数均是实数。传递函数因式分解后，可以写成：(a)零极点形式：K1(sz1)(sz2)(szm)G(s)（２－１１）(sp1)(sp2)(spn)式中，z１，z２，，zｍ称为传递函数的零点；p1，p2，pｎ称为传递函数的极点。(b)时间常数形式：mK(is1)G(s)i1，i,Tj为时间常数nj1(Tjs1)m(i)时间常数和零极点的关系：i1/zi，Tj1/pjKK1in1(Tj)j1由于传递函数中，分子及分母各项系数均为实数，因此传递函数若具有复数零、极点，则其复数零、极点必然是共轭的。’注意，传递函数分子分母多项式若有公因子可以消去，传递函数变成最简单的分式，这称为零极点相消。因此，只有当式(2－９)中的分子及分母多项式收集于网络，如有侵权请联系管理员删除精品文档间没有公因子时，传递函数的零、极点才会和系统的零、极点完全相同；分母多项式的阶次才代表系统的阶次。将传递函数G(s)的零点和极点同时表示在复数平面上的图形，称为传递函数的零极点分布图。（如课本P21图2-10）通常用“Ο”表示传递函数的零点，用“×”表示传递函数的极点。(s2)2(s3)G(s)5)(s22s2)s(s其零极点分布如图：传递函数的零极点分布图2－3－3典型环节的传递函数一个物理系统是由许多元件组合而成的。虽然各种元件的具体结构和作用原理是多种多样的。但若抛开其具体结构和物理特点，研究其运动规律和数学模型的共性，就可以划分成为数不多的几种典型环节。这些典型环节是：比例环节惯性环节积分环节微分环节振荡环节滞后环节。收集于网络，如有侵权请联系管理员删除精品文档典型环节是按照数学模型的共性划分的，它和具体元件不一定是一一对应的。其传递函数相同，动态特性就必然相同。如传递函数中有零值极点、共轭复数零点和极点时，传递函数式(2-9)改写为：m1m2K(is1)(2ks22kks1)G(s)i1k1（２－１２）n1n2s(Tjs1)(T2ls22lTls1)j1l1由此系统的传递函数是由若干典型环节组合而成。1．比例环节(又称放大环节)比例环节的输出量以一定比例复现输入信号，微分方程为：y(ｔ)＝Kx(ｔ)其传递函数为：G(s)Y(s)(2—１３)KX(s)实例如：理想运算放大器，齿轮传动变速箱、电子放大器、电路分压器和机械杠杆等。2．惯性环节(周期环节)Tdy(t)其微分方程为：yt)Kxt(()dt式中，T为时间常数，K为惯性环节增益，其传递函数为：收集于网络，如有侵权请联系管理员删除精品文档Y(s)KG(s)(2－１４)X(s)Ts1例：是由运算放大器构成的惯性环节。其传递函数为：G(s)Y(s)KX(s)Ts1式中，K＝R2／R1，T＝R2C，负号表示输出与输入反向。另有单容液位系统，电热炉炉温随电压变化系统和单容充放气系统也可视为惯性环节。3．积分环节(无差环节)积分环节的输入与输出的微分方程为：Tdy(t)x(t)dtY(s)1(2－15)其传递函数为：G(s)TsX(s)如图是由一个运算放大器构成的积分环节。其传递函数为：Y(s)1G(s)TsX(s)T=RC积分时间收集于网络，如有侵权请联系管理员删除精品文档实际工程中的电子积分器、水槽液位、烤箱温度、电动机转速等系统都属于积分环节。4．微分环节（超前环节）微分环节：理论微分环节实际微分环节(1)理论微分环节：是指仅在理论上存在，而在实际工程中不能单独实现的环节，包括纯微分环节、一阶微分环节和二阶微分环节。纯微分环节传递函数：G(s)Y(s)（2-16）TsX(s)一阶微分环节的传递函数为：G(s)Ts1（2-17）二阶微分环节的传递函数为：G(s)T2s22Ts1（2-18）它们不满足线性系统传递函数的基本性质，即n≥m，所以在实际中不会单独使用。(2)实际微分环节(也称复合微分环节)：如图(a)RC电路，其传递函数为：Ts（2-19）G(s)Ts1式中，T＝RC称为时间常数。如图(b)所示RC电路的传递函数为：K(T1s1)（2-20）G(s)1T2s收集于网络，如有侵权请联系管理员删除精品文档式中：KR2为放大系数；T1＝R1C1；T2R1R2C1为时间常数。R1R2R1R2由于它们满足n≥m的基本条件，所以可以付诸实际使用。5．振荡环节振荡环节的微分方程为：T2d2y(t)2Tdy(t)y(t)x(t)dt2dt式中：T为时间常数；ζ为阻尼比。振荡环节的传递函数为：Y(s)12G(s)n（2-21）X(s)T2s22Ts1s22nsn2式中n1为无阻尼振荡频率。有两个储能元件的系统属二阶系统，如TRLC电路、弹簧质量阻尼系统等可以用振荡环节描述。例RLC电路的传递函数为：G(s)1(2-22)LCs2RCs16．滞后环节滞后环节是输入信号加入后，其输出端隔一时间才能复现输入信号的环节。也称延迟环节或延时环节，其微分方程为：y(t)x(t)式中，τ为滞后时间。滞后环节的传递函数为：G(s)Y(s)es(2-23)X(s)2．4控制系统的结构图及其简化结构图又称为方框图，它是在传递函数的基础上建立的，是描述元、部件动态特性的图示模型。其中的每一个方块代表系统的一个组成环节。必须注．．．．．．收集于网络，如有侵权请联系管理员删除精品文档意，出现在方块图中的环节是以无负载效应为前提的，信息传递关系具有单向性。方块图中的一个环节不一定与实际系统中的一个元件相对应。用方块图表示系统的优点是：可以清楚地表明系统内部信号流动的情况和各环节各变量之间的关系；可以揭示各环节对系统性能的影响；可以根据信号的流向，将各环节的方块图连接起来，得到整个系统的方块图，从而较易写出整个系统的传递函数。2－4－1方块图的组成方块图的主要组成部分包括：函数方框(把一个环节的传递函数写在一个方块里)相加点(加减运算，“+”号表示相加，“-”号表示相减)分支点(同一分支点引出的信号，其信号和数值完全相同)信号流线(方块外面画上带箭头的线段，信息流向单向性)负反馈系统2－４－2绘制方块图的步骤首先写出每个环节的运动方程。特别要注意的是：在列写运动方程时，一定要考虑到相互连接部分间的负载效应，当一个方块的输出与下一个方块的输入相连接，并基本上不受下一方收集于网络，如有侵权请联系管理员删除精品文档块影响时，就认为是没有负载效应的。只有在没有负载效应时，前一方块中的传递函数才可以按独立方块进行计算。根据运动方程，写出传递函数。根据传递函数，画出各环节的方块图。根据信号流向，将各方块联结起来．便得出系统的方块图。例2—13如图2—19所示的RC电路，试绘制其系统结构图。解(1)电路的微分方程组为：uu1R1ii1ii21i1dtu1C1两级RC电路u1yR2i21iidtyC2(2)将上面各式取拉氏变换。取零初始条件，并整理成因果关系式：．．I(s)[U(s)U1(s)]/RI1(s)I(s)I2(s)U1(s)1I1(s)C1sI2(s)[U1(s)Y(s)]/R2Y(s)I2(s)1C2s收集于网络，如有侵权请联系管理员删除精品文档作出相应的方块图，如图(a)所示。将各元件方块图按信号流向联结起来，便得到两级RC网络的方块图，如图(b)所示。注意：图(b)并不等于两个RC网络的方块图的串联，因为两级RC电路之间有负载效应。(再例课本P252-13题)2－４－3方块图的简化利用方块图分析和设计系统时，常常要对方块图的结构进行适当的改动。用方块图求系统的传递函数时，总是要对方块图进行简化。这些统称为方块图的变换或运算。串联环节的简化多个串联方框的总传递函数等于各方框传递函数之乘积，如图。并联环节的简化多个并联方框的等效方框，等于各个方框传递函数的代数和。收集于网络，如有侵权请联系管理员删除精品文档（b）反馈回路的简化如果传递函数分别为G(s)和H(s)的两个方框如图(a)形式连接，就称为反馈连接。“+”为正反馈，“－”则表示相减，为负反馈。由图(a)所知：Y(s)＝G(s)E(s)，B(s)=H(s)Y(s)，E(s)=X(s)±B(s)E(s)——偏差信号，或称比较信号；B(s)——反馈信号。消去E(s)和B(s)后得：Ｙ(s)=G(s)[X(s)±H(s)Y(s)]整理得：(s)Y(s)G(s)X(s)1G(s)H(s)称Φ(s)为闭环传递函数。为了更方便地分析控制系统的传递函数，特引入以下的几个基本概念：(1)G(s)：前向通道传递函数。（由输入信号开始，沿信号流向到输出信号所经过的路径为系统的前向通道，前向通道中所有传递函数之积）(2)H(s)：反向通道传递函数。(由输出信号反向将信号馈送到输入端的通道为系统反向通道，其反向通道的所有传递函数之积)。收集于网络，如有侵权请联系管理员删除精品文档(3)GK(s)：开环传递函数。将一闭环系统从其反馈信号处打开后，以比较环节的输出为输入信号；以反馈信号为输出信号的传递函数。B(s)GK(s)G(s)H(s)E(s)(4)特征方程：若令闭环传递函数的分母等于零，则称其为闭环系统的特征方程。(5)单位反馈：如果反馈通道中的传递函数为1，则称为单位反馈系统。4）比较点移动(相加点移动)如图a为点后移，串接一个与比较点所越过的方框具有相同传递函数的函数方框；如图ｂ为点前移，串接一个与比较点所越过的方框具有相同传递函数成倒数的函数方框。以上比较点移动后是等效的。收集于网络，如有侵权请联系管理员删除精品文档5）分支点移动如图ａ为分支点前移；如图b为分支点后移。常用系统结构图等效变换的基本规则列于表2-1中。应用以上各项基本规则将包含许多反馈回路的复杂方框图简化。例２－1４简化例2－1３中图(b)，求出系统的传递函数。解利用方块图简化规则，可将图(b)的方块图按下列步骤求取传递函数：第一步：将信号线Ｉ2的分支点移到Ｙ处，信号取出线上应增加特性为Ｃ２ｓ的环节；将信号2的相加点移到U处，信号线上应增加特性为Ｒ１的环节。变换结果如图(a)所示。第二步：将串联环节1和1合并，消去其局部反馈；将串联环节1和1R1C1sR2C2s合并，消去其局部反馈；将R1和C2s合并，结果如图(b)所示。第三步：对图(b)的局部反馈进行简化，如图(c)所示。收集于网络，如有侵权请联系管理员删除精品文档第四步：对图(c)的前向通道进行合并，如图(d)所示。第五步：在图(d)中，消去反馈，简化后如图(e)所示。并据此写出系统的传递函数。Y(s)1U(s)R1R2C1C2s2(R1C1R2C2R1C2)s1收集于网络，如有侵权请联系管理员删除精品文档RLC电路结构图及简化：收集于网络，如有侵权请联系管理员删除精品文档对P27图2-21RC电路系统的结构图简化收集于网络，如有侵权请联系管理员删除精品文档收集于网络，如有侵权请联系管理员删除