自动控制数学模型精选PPT.ppt
《自动控制数学模型精选PPT.ppt》由会员分享,可在线阅读,更多相关《自动控制数学模型精选PPT.ppt(94页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、04.04.2023第二章数学模型 关于自动控制数学模型第1页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.1数学模型数学模型概述概述为了从理论上对自动控制系统进行定性分析和定量计算,首先需要建为了从理论上对自动控制系统进行定性分析和定量计算,首先需要建立系统的数学模型。立系统的数学模型。系统的数学模型:描述系统输入、输出变量以及内部各变量之系统的数学模型:描述系统输入、输出变量以及内部各变量之间关系的数学表达式。间关系的数学表达式。常用的动态数学模型有常用的动态数学模型有微分方程、传递函数及微分方程、传递函数及动态结构图。动态结构图。系统数学模型的建立,一般采用解析法或
2、实验法。系统数学模型的建立,一般采用解析法或实验法。解析法:依据系统及元件各变量之间所遵循的物理、化学定律,列写解析法:依据系统及元件各变量之间所遵循的物理、化学定律,列写出变量间的数学表达式,从而建立数学模型。出变量间的数学表达式,从而建立数学模型。本章仅讨论解析法,关于实验法将在后面的章节进行介绍。本章仅讨论解析法,关于实验法将在后面的章节进行介绍。第2页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.1.1线性系统的微分方程模型线性系统的微分方程模型很很多多常常见见的的元元件件或或系系统统的的输输出出量量和和输输入入量量之之间间的的关关系系都都可可以以用用一一个个微微
3、分分方方程程表表示示。微微分分方方程程的的阶阶数数一一般般是是指指方方程程中中最最高高导导数数项项的的阶数阶数,又称为系统的阶数。又称为系统的阶数。如图机械系统,由牛顿定理得到以下关系:如图机械系统,由牛顿定理得到以下关系:第3页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 如图如图RLC网络,由电路定律可得网络,由电路定律可得:不同的物理系统可能得到相似的数学表达式。如果它们对应的不同的物理系统可能得到相似的数学表达式。如果它们对应的系数和初始条件相同,则它们的解将完全相同。这样系数和初始条件相同,则它们的解将完全相同。这样就可以撇开系就可以撇开系统的具体物理属性,研究这些
4、系统的运动过程的共同规律。统的具体物理属性,研究这些系统的运动过程的共同规律。有了数学表达式,就可从理论上进行普遍意义上的分析。有了数学表达式,就可从理论上进行普遍意义上的分析。第4页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 机械系统中,设外力机械系统中,设外力F1,质量,质量m2,弹性系数,弹性系数k1,若阻尼系,若阻尼系数较小数较小1,则发生震荡,若阻,则发生震荡,若阻尼系数较大尼系数较大10,不会产生震荡。,不会产生震荡。但无但无论阻尼大小如何,最终物体将下降一个单位长度,新增的弹力正论阻尼大小如何,最终物体将下降一个单位长度,新增的弹力正好和外力相抵,系统进入一个
5、新的平衡点。好和外力相抵,系统进入一个新的平衡点。总之,建立合理的数学模型,是至关重要的问题。许多系统,事件及总之,建立合理的数学模型,是至关重要的问题。许多系统,事件及项目就是因为无法建立合理的数学模型而不能加以预测和控制。项目就是因为无法建立合理的数学模型而不能加以预测和控制。第5页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.1.2 2.1.2 列写微分方程的一般方法列写微分方程的一般方法用解析法列写系统或元件微分方程的一般步骤是:用解析法列写系统或元件微分方程的一般步骤是:1 1根根据据实实际际工工作作情情况况,将将系系统统划划分分为为多多个个独独立立的的环环节节
6、,标标出出各各环环节节的输入、输出变量。的输入、输出变量。各环节之间无负载效应。各环节之间无负载效应。2从输入端开始,按照信号的传递顺序,依据各环节所遵循的从输入端开始,按照信号的传递顺序,依据各环节所遵循的物理定律,列写的动态方程,一般为微分方程组。物理定律,列写的动态方程,一般为微分方程组。3消去中间变量,写出系统输入、输出变量的微分方程。消去中间变量,写出系统输入、输出变量的微分方程。4标准化。即将与输入有关的各项放在等号右侧,与输出有关标准化。即将与输入有关的各项放在等号右侧,与输出有关的各项放在等号左侧,并按降幂排列。最后将系数归化为具有一定物的各项放在等号左侧,并按降幂排列。最后将
7、系数归化为具有一定物理意义的形式。理意义的形式。第6页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 例例2.1 2.1 列写如图所示列写如图所示RCRC滤波电路的微分方程。(假设电路的输入电滤波电路的微分方程。(假设电路的输入电源的内阻为零,输出接的负载具有无限大阻抗)源的内阻为零,输出接的负载具有无限大阻抗)解解 根据基尔霍夫定律得:根据基尔霍夫定律得:消除中间变量,得到滤波网络的微分方程式为消除中间变量,得到滤波网络的微分方程式为 :第7页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 若若撇撇开开具具体体系系统统的的物物理理属属性性,令令r(t)r(t)
8、为为输输入入,c(t)c(t)为为输输出出。线线性性n阶阶系统的输入输出微分方程式的一般表达式可写为系统的输入输出微分方程式的一般表达式可写为 式中式中均为由系统结构参数决定的常系数,均为由系统结构参数决定的常系数,且有且有nm。令令 则上式可改写为:则上式可改写为:第8页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.1.3 非线性数学模型的线性化非线性数学模型的线性化 在在建建立立控控制制系系统统的的数数学学模模型型时时,常常常常遇遇到到非非线线性性的的问问题题。严严格格地地讲讲,实实际际的的物物理理系系统统都都包包含含着着不不同同程程度度的的非非线线性性因因素素。但但
9、是是,许许多多非非线线性性系系统统在在一定的条件下可以近似地视作线性系统。一定的条件下可以近似地视作线性系统。若若控控制制系系统统在在工工作作点点的的附附近近微微小小运运动动,则则可可将将非非线线性性函函数数展展开开为为泰泰勒勒级级数数,并并忽忽略略级级数数展展开开式式中中的的高高次次项项,从从而而得得到到只只含含一一次次项项的的线线性性化化方方程。即用工作点的切线代替非线性曲线。程。即用工作点的切线代替非线性曲线。对对于于一一般般的的非非线线性性系系统统,假假设设其其输输入入量量为为r r,输输出出量量为为c c,并并设设在在给给定定工工作作点点处处c0 0=f(r 0 0),各各阶阶导导数
10、数均均存存在在,则则可可在在 的的邻邻域域展展开泰勒级数,即开泰勒级数,即第9页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 当当(rr 0 0),很小时,可以忽略上式中二阶以上各项,得很小时,可以忽略上式中二阶以上各项,得 或或 在处理非线性问题时,应注意以下几点:在处理非线性问题时,应注意以下几点:1 1线线性性化化是是在在输输入入、输输出出量量围围绕绕平平衡衡点点作作小小范范围围变变化化的的假假设设下下进行的。一般取零误差状态作为平衡工作状态。进行的。一般取零误差状态作为平衡工作状态。2 2线线性性化化以以切切线线代代替替曲曲线线,是是一一种种近近似似处处理理。系系统统
11、的的实实际际变变化化量量如如果果很很大大,则则采采用用小小偏偏差差线线性性模模型型将将会会带带来来较较大大的的计计算算误误差。差。3.3.对对于于某某些些严严重重的的典典型型非非线线性性,不不能能进进行行求求导导运运算算,因因此此原原则则上上不不能能用小偏差法进行线性化用小偏差法进行线性化第10页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 例例2.22.2 图图示示为为一一个个单单摆摆系系统统,输输入入量量M为为零零(不不加加外外力力矩矩),),输输出出量量为为摆摆幅幅(t)。摆摆锤锤的的质质量量为为m,摆摆杆杆长长度度为为l,l,空空气气阻阻尼尼系系数数为为,重重力加速度
12、为力加速度为g。试建立系统的近似线性运动方程。试建立系统的近似线性运动方程。解解 对对于于图图示示的的单单摆摆系系统统,根根据据牛牛顿顿运运动动定定律律可可以以直直接接推推出出如如下下系系统统运动方程运动方程:显显然然方方程程是是一一个个二二阶阶的的非非线线性性微微分分方方程程(因因为为含含有有sinsin),),但但是是在在摆摆幅幅较较小小的的情情况况下下,将其线性化处理:将其线性化处理:第11页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 令非线性函数令非线性函数sin()=f,则,则工作点在工作点在0=0,f0=0。线线性化:性化:即单摆系统的近似线性化动态方程为:即单
13、摆系统的近似线性化动态方程为:第12页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.2传传递递函函数数2.2.1拉氏变换拉氏变换 1.1.拉氏变换的定义拉氏变换的定义 将时间函数将时间函数f(t)乘上指数函数乘上指数函数e-st(其中其中s=+j是一个复数是一个复数),),并且并且在在0,+0,+上对上对t积分积分,称为称为f(t)的拉氏变换的拉氏变换,并用并用L f(t)表示。表示。拉拉氏氏变变换换将将原原来来的的时时间间函函数数f(t)转转化化为为复复变变量量函函数数F(s)。通通常常将将F(s)称作称作f(t)的的象函数象函数,将将f(t)称作称作F(s)的的原函数
14、原函数。传递函数是对微分方程取传递函数是对微分方程取拉氏变换后推导出来的概念。拉氏变换后推导出来的概念。第13页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.2.拉氏变换的计算拉氏变换的计算1)1)根据定义积分计算,各典型函数的拉氏变换根据定义积分计算,各典型函数的拉氏变换见下表见下表。第14页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2)MATLAB计算计算symsst;Ft=1-sin(t)Fs=laplace(Ft,t,s)执行结果:执行结果:Fs=1/s-1/(s2+1)3.拉氏反变换拉氏反变换已知时间函数的象函数通过拉氏反变换求出其时间函数
15、:已知时间函数的象函数通过拉氏反变换求出其时间函数:第15页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 1)部分分式法部分分式法将将F(s)展开成多个典型函数的象函数之代数和,查表。展开成多个典型函数的象函数之代数和,查表。例例2.3F(s)含单极点和重极点时的拉氏反变换。含单极点和重极点时的拉氏反变换。解解:第16页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2)MATLAB拉氏反变换指令:拉氏反变换指令:ilaplace(Fs,s,t)例例2.3的的MATLAB求解程序:求解程序:symss,t;ilaplace(1/s*(s+3)*(s+1)2)计
16、算结果与手算结果完全一样。计算结果与手算结果完全一样。第17页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 例例2.4 F(s)含有共轭复极点时的反变换。解:解:用用MATLAB求解:求解:symsst;ft=ilaplace(s+1)/s*(s2+s+1);pretty(ft)%将符号表达式写成易读形式将符号表达式写成易读形式与手算结果一样与手算结果一样第18页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 4.4.拉氏变换的基本定理拉氏变换的基本定理1)1)线性定理线性定理两个函数和的拉氏变换两个函数和的拉氏变换,等于每个函数拉氏变换的和等于每个函数拉氏变
17、换的和,即即函数放大函数放大k倍的拉氏变换等于该函数拉氏变换的倍的拉氏变换等于该函数拉氏变换的k倍倍,即即第19页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2)2)微分定理微分定理成立成立,则有则有如果初始条件如果初始条件第20页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 3)3)终值定理终值定理 函函数数 f(t)在在 t+时时的的函函数数值值(即即稳稳定定值值)可可以以通通过过 f(t)的的拉拉氏变换氏变换F(s)乘以乘以 s 取取 s0时的极限而得到时的极限而得到,即即总结:微分方程通过拉氏变换变成代数方程,解代数方程可求总结:微分方程通过拉氏变
18、换变成代数方程,解代数方程可求出输出的象函数,对象函数取拉反变换,可求出微分方程的解。出输出的象函数,对象函数取拉反变换,可求出微分方程的解。第21页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.2.2 2.2.2 传递函数的定义和特点传递函数的定义和特点 1.1.传递函数的定义传递函数的定义 线线性性定定常常系系统统的的传传递递函函数数,定定义义为为零零初初始始条条件件下下,系系统统输输出出量量的的拉拉氏变换与输入量的拉氏变换之比。氏变换与输入量的拉氏变换之比。设输入量为设输入量为r(t)(t);输出量为;输出量为 c(t),定义传递函数为:,定义传递函数为:一般线性定
19、常系统由下面的一般线性定常系统由下面的n阶线性常微分方程描述阶线性常微分方程描述:第22页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 如如果果r(t)和和c(t)及及其其各各阶阶导导数数在在t=0时时的的值值均均为为零零,则则根根据据拉拉氏氏变变换换的的定定义和性质义和性质,对微分方程进行拉氏变换对微分方程进行拉氏变换,可得可得由传递函数的定义可得系统的多项式形式的传递函数为由传递函数的定义可得系统的多项式形式的传递函数为第23页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 用用MATLAB指令:指令:Gs=tf(b0,b1,bm,a0,a1,an)或者或
20、者s=tf(s);Gs=关于关于s的多项式的多项式构造多项式形式的传递函数后,可以用构造多项式形式的传递函数后,可以用MATLAB的各种控制系统的各种控制系统指令分析系统。指令分析系统。传递函数的零极点形式传递函数的零极点形式zi(i=1,2,m)和和pj(j=1,2,n)分别称为传递函数的分别称为传递函数的零点和极点零点和极点,K1称为称为传递函数的增益或传递函数的增益或根轨迹增益根轨迹增益。第24页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 i(i=1,2,=1,2,m)和和Tj(j=1,2,=1,2,n)为为系系统统中中各各环环节节的的时时间间常常数数,K为为系系统统
21、的的放大倍数放大倍数。用用MATLAB指令:指令:Gs=zpk(z0,z1,zm,p0,p1,pn,K)或者或者s=tf(s);Gs=关于关于s的因式的因式可构造零极点形式的传递函数。可构造零极点形式的传递函数。传递函数的参数形式传递函数的参数形式第25页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 使用使用Gtf=tf(Gzpk)或者或者Gzpkzpk(Gtf)可实现传递函数在零极点形式可实现传递函数在零极点形式和多项式形式之间的互换。即可将传递函数进行展开和因式分解。和多项式形式之间的互换。即可将传递函数进行展开和因式分解。例例2.4求传递函数求传递函数的零极点形式。的零
22、极点形式。解解G=tf(26,4,1,14,63,90);F=zpk(G)执行结果:执行结果:Zero/pole/gain:2(s+2)(s+1)-(s+6)(s+5)(s+3)第26页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 2.2.传递函数的特点传递函数的特点 (1)(1)传传递递函函数数的的概概念念适适用用于于线线性性定定常常系系统统,传传递递函函数数的的结结构构和和各各项项系系数数(包包括括常常数数项项)完完全全取取决决于于系系统统本本身身结结构构,因因此此,它它是是系系统统的的动动态态数数学学模模型型,而而与与输输入入信信号号的的具具体体形形式式和和大大小小无无
23、关关,也也不不反反映映系系统统的的任任何何内内部部信信息息。同同一一个个系系统统若若选选择择不不同同的的量量作作为为输输入入量量和和输输出出量量,所所得得到到的的传传递函数可能不同。所以谈到传递函数递函数可能不同。所以谈到传递函数,必须指明输入量和输出量必须指明输入量和输出量。已知传递函数,可求任意输入已知传递函数,可求任意输入R(s)R(s)下的输出下的输出C(s)C(s):第27页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 (2)(2)传传递递函函数数是是在在零零初初始始条条件件下下定定义义的的。但但是是,对对输输入入量量加加于于系统之前系统之前,系统处于系统处于稳定工
24、作状态稳定工作状态的情况同样适用。的情况同样适用。(3)(3)对对于于实实际际的的物物理理系系统统和和元元件件而而言言,传传递递函函数数的的分分子子多多项项式式的的阶阶次次总总是是小小于于分分母母多多项项式式的的阶阶次次,即即mn。它它反反映映了了一一个个基基本本事事实实:一一个个物物理理系系统统的的输输出出不不可可能能立立即即复复现现输输入入信信号号,只只有有经经过过一一段段时时间间后后,输出量才能达到输入量所要求的数值。输出量才能达到输入量所要求的数值。第28页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 (4)(4)传传递递函函数数与与线线性性常常微微分分方方程程一一一
25、一对对应应。将将传传递递函函数数展展开开并并取取拉拉氏反变换可得到微分方程。例如氏反变换可得到微分方程。例如,由传递函数由传递函数可得可得s的代数方程的代数方程(a0s2+a1s+a2)C(s)=(b1s+b2)R(s)对方程两端取拉氏反变换,对方程两端取拉氏反变换,便得到相应的微分方程便得到相应的微分方程第29页,讲稿共94张,创作于星期三04.04.2023第二章数学模型 (5)(5)传传递递函函数数不不能能反反映映系系统统或或元元件件的的学学科科属属性性和和物物理理性性质质。物物理理性性质质和和学学科科类类别别截截然然不不同同的的系系统统可可能能具具有有完完全全相相同同的的传传递递函函数
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 自动控制 数学模型 精选 PPT
限制150内