二阶系统的时间响应课件.ppt
《二阶系统的时间响应课件.ppt》由会员分享,可在线阅读,更多相关《二阶系统的时间响应课件.ppt(188页珍藏版)》请在三一办公上搜索。
1、第2章 数 学 模 型,为了从理论上对控制系统进 行性能分析,首先要建立系统的数学模型。系统的数学模型,是描述系统输入、输出量以及内部各变量之间关系的数学表达式,它揭示了系统结构及其参数与其性能之间的内在关系。系统数学模型有多种形式,这取决于变量和坐标系统的选择。在时间域,通常采用微分方程或一阶微分方程组的形式;在复数域则采用传递函数形式;而在频率域采用频率特性形式。必须指出,建立合理的数学模型,对于系统的分析和研究极为重要。由于不可能将系统实际的错综复杂的物理现象完全表达出来,因而要对模型的简洁性与精确性进行折衷的考虑。一般是根据系统的实际结构参数和系统分析所要求的精度,忽略一些次要因素,建
2、立既能反映系统内在本质特性,又能简化分析计算工作的模型。,建立系统数学模型,一般采用解析法或实验法。所谓解析法建模,即依据系统及元件各变量之间所遵循的物理学定律,理论推导出变量间的数学关系式,从而建立数学模型。本章仅讨论解析建模方法,关于实验法建模将在后面的章节进行介绍。2.1 控制系统的运动微分方程 2.1.1 建立数学模型的一般步骤 用解析法列写系统或元件微分方程的一般步骤是:(1)分析系统的工作原理和信号传递变换的过程,确定系统和各元件的输入、输出量。(2)从系统的输入端开始,按照信号传递变换过程,依据各变量 所遵循的物理学定律,依次列写出各元件、部件动态微分方程。(3)消去中间变量,得
3、到一个描述元件或系统输入、输出变量之间关系的微分方程。(4)写成标准化形式。将与输入有关的项放在等式右侧,与输出有关的项放在等式的左侧,且各阶导数项按降幂排列。,2.1.2 控制系统微分方程的列写,1机械系统 任何机械系统的数学模型都可以应用牛顿定律来建立。机械系统中以各种形式出现的物理现象,都可以使用质量、弹性和阻尼三个要素来描述。(1)机械平移系统 图2.1所示为常见的质量-弹簧-阻尼系统,图中的、分别表示质量、弹簧刚度和粘性阻尼系数。以系统在静止平衡时的那一点为零点,即平衡工作点,这样的零位选择消除了重力的影响。设系统的输入量为外作用力,输出量为质量块的位移。现研究外力 与位移 之间的关
4、系。在输入力 的作用下,质量块 将有加速度,从而产生速度和位移。质量块的速度和位移使阻尼器和弹簧产生粘性阻尼力 和弹性力。这两个力反馈作用于质量块上,影响输入 的作用效果,从而使质量块的速度和位移随时间发,图2.1 机械平移系统力学模型,生变化,产生动态过程。根据牛顿第二定律,有 点击观看公式推导 由阻尼器、弹簧的特性,可写出 由以上三个式子,消去 和,并写成标准形式,得 一般、均为常数,故式(2.1)为二阶常系数线性微分方程。它描述了输入 和输出 之间的动态关系。方程的系数取决于系统的结构参数;而方程的阶次等于系统中独,(2.1),立的储能元件(惯性质量、弹簧)的数量。当质量 很小可忽略不计
5、时,系统由并联的弹簧和阻尼器组成,如图2.2所示。此时,系统的运动方程为一阶常系数微分方程 这说明,同一系统由于简化程度的不同,可以有不同的数学模型。,图2.2 弹簧-阻尼系统力学模型,(2)机械旋转系统 包含定轴旋转的机械系统用途极其广泛。其建模方法与平移系统非常相似。只是这里将质量、弹簧、阻尼分别变成转动惯量、扭转弹簧、旋转阻尼。图2.3所示为一机械旋转系统,旋转体通过柔性轴(用扭转弹簧 表示)与齿轮连接。旋转体在粘性介质中旋转,因而承受与旋转速度成正比的阻尼力矩。设齿轮转角 为系统输入量,旋转体转角 为系统输出量,据此建立系统的运动微分方程(忽略轴承上的摩擦)。扭转弹簧左、右端的转角分别
6、为、,设它加给旋转体的扭矩为(当 时,弹簧的扭矩为零),则 旋转体上除了受弹簧的扭矩外,也受阻尼扭矩 作用,因而有扭矩平衡方程,和旋转阻尼特性方程 由以上三式整理可得机械旋转系统运动微分方程,图2.3 机械旋转系统力学模型,(2.2),2电气系统 电阻、电感 和电容器 是电路中的三个基本元件。通常利用基尔霍夫定律来建立电气系统的数学模型。电气系统数学模型 无源电路网络如图2.4所示,设输入端电压 为系统输入量。电容器 两端电压 为系统输出量。现研究输入电压 和输出电压 之间的关系。电路中的电流为中间变量。,图2.4 无源电路网络,根据基尔霍夫定律,有 点击观看公式推导 消去中间变量,稍加整理,
7、即得 一般假定、都是常数,则上式为二阶常系数线性微分方程。若,系统也可简化为一阶常微分方程 有源电路网络如图2.5所示,设电压 为系统输入量,电压 为系统输出量。现建立 与 之间的关系式。,(2.3),(2.4),图2.5 有源电路网络,图中 点为运算放大器的反相输入端,为运算放大器的开环放大倍数。因为 且一般 值很大,所以 点电位 运算放大器的输入阻抗一般都很高,故而可认为 因此,可以得到 即,(2.5),3流体系统 流体系统比较复杂,但经过适当简化也可以用微分方程加以描述。图2.6所示为一简单的液位控制系统。在此系统中,箱体通过输出端的 节流阀对外供液。设流入箱体的流量 为系统输入量,液面
8、高度 为输出量,下面列写液位波动的运动微分方程。,图2.6 液位控制系统,根据流体连续方程,可得 式中 箱体的截面积。设液体是不可压缩的,通过节流阀的液流是紊流,则其流量公式为 式中 由节流阀通流面积和通流口结构形式决定的系数,通流面积不变时 为常数。消去中间变量 得液位波动方程为 显然,式(2.8)是一个非线性微分方程。4模型分析 将上述系统模型进行比较,可清楚地看到,物理本质不同的,(2.6),(2.7),(2.8),系统,可以有相同的数学模型。反之,同一数学模型可以描述物理性质完全不同的系统。因此,从控制理论来说,可抛开系统的物理属性,用同一方法进行普遍意义的分析研究,这就是信息方法,从
9、信息在系统中传递、转换的方面来研究系统的功能。而从动态性能来看,在相同形式的输入作用下,数学模型相同而物理本质不同的系统其输出响应相似,若方程系数等值则响应完全一样,这样就有可能利用电系统来模拟其它系统,进行实验研究。这就是控制理论中的功能模拟方法的基础。分析上述系统模型还可以看出,描述系统运动的微分方程的系数都是系统的结构参数及其组合,这就说明系统的动态特性是系统的固有特性,取决于系统结构及其参数。用线性微分方程描述的系统,称为线性系统。如果方程的系数为常数,则称为线性定常系统;如果方程的系数不是常数,而是时间 的函数,则称为线性时变系统。线性系统的特点是具有线性性质,即服从叠加原理。这个原
10、理是说,多个输入同时作用,于线性系统的总响应,等于各输入单独作用时产生的响应之和。用非线性微分方程描述的系统称为非线性系统,如前述的液位控制系统。在工程实践中,可实现的线性定常系统,均能用 阶常系数线性微分方程来描述其运动特性。设系统的输入量为,系统的输出量为,则单输入、单输出 阶系统常系数线性微分方程有如下的一般形式(2.9)式中,和,由系统结 构参数决定的实常数。由于实际系统中总含有惯性元件以及受到能源能量的限制,所以总是,2.2 拉氏变换与反变换,机电控制工程所涉及的数学问题较多,经常要解算一些线性微分方程。按照一般方法解算比较麻烦,如果用拉普拉斯变换求解线性微分方程,可将经典数学中的微
11、积分运算转化为代数运算,又能够单独地表明初始条件的影响,并有变换表可查找,因而是一种较为简便的工程数学方法。2.2.1 拉普拉斯变换的定义 如果有一个以时间 为自变量的实变函数,它的定义 域是,那么 的拉普拉斯变换定义为 式中,是复变数,(、均为实数),称为拉普拉斯积分;是函数 的拉普拉斯变换,它是一个复变函数,通常也称 为 的象函数,而称 为 的原函数;L是表示进行拉普拉斯变换的符号。,(2.10),式(2.10)表明:拉氏变换是这样一种变换,即在一定条件下,它能把一实数域中的实变函数变换为一个在复数域内与之等价的复变函数。2.2.2 几种典型函数的拉氏变换 1.单位阶跃函数 的拉氏变换 单
12、位阶跃函数是机电控制中最常用的典型输入信号之一,常以它作为评价系统性能的标准输入,这一函数定义为 单位阶跃函数如图2.7所示,它表示在 时刻突然作用于系统一个幅值为1的不变量。单位阶跃函数的拉氏变换式为 当,则。,所以 2.指数函数 的拉氏变换 指数函数也是控制理论中经常用到的函数,其中 是常数。令 则与求单位阶跃函数同理,就可求得,(2.11),(2.12),图2.7 单位阶跃函数,3.正弦函数与余弦函数的拉氏变换 设,则 由欧拉公式,有 所以,(2.13),同理 4.单位脉冲函数(t)的拉氏变换 单位脉冲函数是在持续时间 期间幅值为 的矩形波。其幅值和作用时间的乘积等于1,即。如图2.8所
13、示。单位脉冲函数的数学表达式为,图2.8单位脉冲函数,其拉氏变换式为此处因为 时,故积分限变为,(2.15),5.单位速度函数的拉氏变换 单位速度函数,又称单位斜坡函数,其数学表达式为见图2.9所示。单位速度函数的拉氏变换式为,图2.9 单位速度函数,利用分部积分法 令 则 所以当 时,,则,(2.16),6.单位加速度函数的拉氏变换 单位加速度函数的数学表达式为 如图2.10所示。其拉氏变换式为 通常并不根据定义来求解象函数和原函数,而可从拉氏变换表(见附录A)中直接查出。,图2.10 单位加速度函数,(2.17),2.2.3 拉氏变换的主要定理 根据拉氏变换定义或查表能对一些标准的函数进行
14、拉氏变换和反变换,但利用以下的定理,则对一般的函数可使运算简化。1.叠加定理 拉氏变换也服从线性函数的齐次性和叠加性。1)齐次性 设,则 式中 常数。(2)叠加性 设,则 两者结合起来,就有 这说明拉氏变换是线性变换。,(2.18),(2.19),2.微分定理 设 则 式中 函数 在 时刻的值,即初始值。同样,可得 的各阶导数的拉氏变换是,(2.20),式中,原函数各阶导数在 时刻的值。如果函数 及其各阶导数的初始值均为零(称为零初始条件),则 各阶导数的拉氏变换为 3.复微分定理 若 可以进行拉氏变换,则除了在 的极点以外,式中,。同样有,(2.21),(2.22),一般地,有 4.积分定理
15、设,则 式中 积分 在 时刻的值。当初始条件为零时,对多重积分是,(2.23),(2.24),(2.25),(2.26),当初始条件为零时,则 5.延迟定理 设,且 时,则 函数为原函数 沿时间轴延迟了,如图2.11所示。,(2.27),(2.28),图2.11 函数,6.位移定理 在控制理论中,经常遇到 一类的函数,它的象函数只需把 用 代替即可,这相当于在复数 坐标中,有一位移。设,则 例如 的象函数,则 的象函数为 7.初值定理 它表明原函数在 时的数值。即原函数的初值等于 乘以象函数的终值。,(2.29),(2.30),8.终值定理 设,并且 存在,则 即原函数的终值等于 乘以象函数的
16、初值。这一定理对于求瞬态响应的稳态值是很有用的。9.卷积定理 设,则有即两个原函数的卷积分的拉氏变换等于它们象函数的乘积。式(2.32)中,为卷积分的数学表示,定义为 10.时间比例尺的改变,(2.31),(2.32),式中 比例系数 例如,的象函数,则 的象函数为 11.拉氏变换的积分下限 在某些情况下,在 处有一个脉冲函数。这时必须明确拉普拉斯积分的下限是 还是,因为对于这两种下限,的拉氏变换是不同的。为此,可采用如下符号予以区分:,(2.33),若在 处 包含一个脉冲函数,则 因为在这种情况下显然,如果 在 处没有脉冲函数,则有2.2.4 拉普拉斯反变换 拉普拉斯反变换的公式为 式中 表
17、示拉普拉斯反变换的符号 通常用部分分式展开法将复杂函数展开成有理分式函数之和,然后由拉氏变换表一一查出对应的反变换函数,即得所求的原函数。,(2.36),1.部分分式展开法 在控制理论中,常遇到的象函数是 的有理分式 为了将 写成部分分式,首先将 的分母因式分解,则有 式中,是 的根的负值,称为 的极点,按照这些根的性质,可分为以下几种情况来研究。2.的极点为各不相同的实数时的拉氏反变换,,,,,(2.37),式中,是待定系数,它是 处的留数,其求法如下再根据拉氏变换的迭加定理,求原函数 例 2.1 求 的原函数。解:首先将 的分母因式分解,则有,(2.38),即得 3.含有共轭复数极点时的拉
18、氏反变换 如果 有一对共轭复数极点,其余极点均为各不相同的实数极点。将 展成,式中,和 可按下式求解即因为(或)是复数,故式(2.39)两边都应是复数,令等号两边的实部、虚部分别相等,得两个方程式,联立求解,即得、两个常数。例 2.2 已知,试求其部分分式。解:因为,(2.39),、,(2.40),含有一对共轭复数极点,和一个极点,故可将 式(2.40)因式分解成以下求系数、和。由式(2.40)和式(2.41)相等,有用 乘以上式两边,并令,得到,(2.41),(2.42),上式可进一步写成由上式两边实部和虚部分别相等,可得联立以上两式,可求得为了求出系数,用 乘方程(2.42)两边,并令,将
19、 代入,得将所求得的、值代入(2.41),并整理后得 的部分分式,查拉氏变换表便得,结果见式(3.16)。例 2.3 已知 求。解:将 的分母因式分解,得,利用方程两边实部、虚部分别相等得 解得,所以,,,这种形式再作适当变换:查拉氏变换表得,4.中含有重极点的拉氏反变换 设 有 个重根,则 将上式展开成部分分式 式中,的求法与单实数极点情况下相同。,的求法如下:,(2.43),例 2.4 设,试求 的部分分式。解:已知 含有2个重极点,可将式(2.45)的分母因式分解得 以下求系数、和。,(2.44),(2.45),(2.46),、,将所求得的、值代入式(2.46),即得 的部分分式 查拉氏
20、变换表可得。例 2.5 求 的拉氏反变换。解:将 展开为部分分式,上式中各项系数为 于是查拉氏变换表,得,应当指出,对于在 分母中包含有较高阶次多项式的复杂函数,用人工算法进 行部分分式展开则相当费时费力。这种情况下,采用MATLAB工具就方便多了。5.用MATLAB展开部分分式(1)概述 MATLAB是美国Math Works公司的软件产品,是一个高级的数值分析、处理与计算的软件,其强大的矩阵运算能力和完美的图形可视化功能,使得它成为国际控制界应用最广的首选计算机工具。SIMULINK是基于模型化图形的动态系统仿真软件,是MATLAB的一个工具箱,它使系统分析进入一个崭新的阶段,它不需要过多
21、地了解数值问题,而是侧重于系统的建模、分析与设计。其良好的人机界面及周到的帮助功能使得它广为科技界和工程界所采用。(2)用MATLAB进行部分分式展开,MATLAB有一个命令用于求B(s)/A(s)的部分分式展开式。设s 的有理分式为式中(i=)和(j=)的某些值可能为零。在MATLAB的行向量中,num和den分别表示F(s)分子和分母的系数,即 num=den=1 命令 MATLAB将按下式给出F(s)部分分式展开式中的留数、极点和余项:,r,p,k=residue(num,den),上式与式(2.37)比较,显然有p(1)=-,p(2)=-,p(n)=-;r(1)=,r(2)=,r(n)
22、=;k(s)是余项。例2.6 试求下列函数的部分分式展开式 解:对此函数有 num=1 11 39 52 26den=1 10 35 50 24 命令于是得到下列结果 r,p,k=residue(num,den)r=1.0000 2.5000-3.0000,r,p,k=residue(num,den),0.5000 p=-4.0000-3.0000-2.0000-1.0000 k=1 则得如果F(s)中含重极点,则部分分式展开式将包括下列诸项式中,p(j)为一个q重极点。例2.7 试将下列函数展开成部分分式,解:对于该函数有 num=0 1 4 6 den=1 3 3 1 命令 r,p,k=r
23、esidue(num,den)将得到如下结果:r,p,k=residue(num,den)r=1.0000 2.0000 3.0000 p=-1.0000-1.0000-1.0000k=,所以可得 注意,本例的余项k为零。2.2.5 应用拉氏变换解线性微分方程 应用拉氏变换解线性微分方程时,采用下列步骤:(1)对线性微分方程中每一项进行拉氏变换,使微分方程变为 的代数方程;(2)解代数方程,得到有关变量的拉氏变换表达式;(3)用拉氏反变换得到微分方程的时域解。整个求解过程如图2.12所示。,图2.12 应用拉氏变换法求解线性微分方程的过程,设系统微分方程为若,初始条件分别为、,试求。解:对微分
24、方程左边进行拉氏变换 利用迭加定理将上式逐项相加,即得方程左边的拉氏变换 对方程右边进行拉氏变换,例2.8,得 写成一般形式 应该强调指出 是微分方程的特征方程,也是该系统的特征方程。利用部分分式将 展开为,求待定系数、:代入原式得,查拉氏变换表得当初始条件为零时,得 2.3 传 递 函 数 在控制工程中,直接求解系统微分方程是研究分析系统的基本方法。系统方程的解就是系统的输出响应,通过方程的表达式,可以分析系统的动态特性,可以绘出输出响应曲线,直观地反映系统的动态过程。但是,由于求解过程较为繁琐,计算复杂费时,而且难以直接从微分方程本身研究和判断系统的动态性能,因此,这种方法有很大的局限性。
25、显然,仅用微分方程这一数学模型来进行系统分析设计,显得十分不便。,对于线性定常系统,传递函数是常用的一种数学模型,它是 在拉氏变换的基础上建立的。用传递函数描述系统可以免去求解微分方程的麻烦,间接地分析 系统结构及参数与系统性能的关系,并且可以根据传递函数在复平面上的形状直接判断系统的动态性能,找出改善系统品质的方法。因此,传递函数是经典控制理论的基础,是一个极其重要的基本概念。2.3.1 传递函数的概念和定义 对于线性定常系统,在零初始条件下,系统输出量的拉氏变换与引起该 输出的输入量的拉氏变换之比,称为系统的传递函数。图2.1所示质量-弹簧-阻尼系统,由二阶微分方程式(2.1)来描述它的动
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 系统 时间 响应 课件
链接地址:https://www.31ppt.com/p-3827687.html