系统时间响应及其仿真.ppt
《系统时间响应及其仿真.ppt》由会员分享,可在线阅读,更多相关《系统时间响应及其仿真.ppt(31页珍藏版)》请在三一办公上搜索。
1、第七章 系统时间响应及其仿真,本章主要内容如下:7.1 仿真算法7.2 系统仿真的MATLAB函数,系统的时间响应是指系统在输入信号或初始状态作用下,系统输出随时间变化的情况。系统的时间响应反映了系统的特征和性能,如快速性、稳定性等。对系统时间响应的分析是我们设计、校正系统的基础。,第七章 系统时间响应及其仿真,7.1 仿真算法 对系统的时间响应进行动态仿真,采用什么样的仿真算法是一个至关重要的问题。对连续时间系统进行数字动态仿真,主要是两种方法:基于数值积分的仿真方法;基于离散相似法的仿真方法。,由于后者涉及到离散控制系统理论,因此本节重点介绍基于数字积分的连续系统仿真方法。,7.1 仿真算
2、法,7.1.1 数值积分的基本原理设微分方程为 则求解方程中函数y(t)问题,就是已知初值的常微分方程求解问题。所谓数值求解就是要在时间区间a,b中取若干离散点 求出微分方程在这些时刻的近似值常微分方程数值求解的基本方法是数值积分法。,7.1 仿真算法,7.1.2 数值积分的基本原理 数值积分法:积分区间的划分,将区间a,b分成 N个小区间,时间间隔h()也称为积分步长,在第k个间隔 t=tk,tk+1内积分:,则可用yk(k=0,1,N)作为解y(t)的近似值,如图所示。,数值积分法:,7.1 仿真算法,7.1.2 数值积分的基本原理泰勒展开法:为避免(2)式中的积分项,将y在tk,以h为增
3、量展开成Taylor级数:,式(3)是一个递推公式。积分值与实际微分方程解的误差取决于步长h和计算所用的阶数,它是数值积分的基础。,7.1 仿真算法,7.1.2 数值积分的基本原理有关概念单步法和多步法 单步法指计算yk+1值只需利用tk时刻的信息,也称为自启动算法;多步法在计算yk+1值时,则需利用tk,tk-1,时刻的信息。,显示法和隐式法 显示法在计算yk+1时所需数据均已算出;隐式法在计算yk+1时需用到tk+1时刻的数据,该算法必须借助予估公式。,定步长和变步长 定步长为积分步长在仿真运行过程中始终不变;变步长指在仿真运行过程中自动修改步长。,7.1 仿真算法,7.1.3 数值积分的
4、几个算法欧拉算法 在(3)式中取前两项:可得欧拉算法:,【说明】欧拉法是用一条过各点的切线取代曲线来逼近精确解。该算法简单,计算量小,但精度较低。,7.1.3 数值积分的几个算法梯度法 梯度法是欧拉法的改进。与欧拉法相比,梯度法是用两个点(tm,ym)、(tm+1,ym+1)的斜率的平均值来确定下一点的y值。由于上式计算时需要用到ym+1的值,而 ym+1不能预先知道,故梯度法需要和欧拉法结合使用,即用欧拉法对ym+1 进行予估,再由梯度法计算ym+1,7.1 仿真算法,7.1.3 数值积分的几个算法龙格-库塔法龙格-库塔法的基本思想欧拉算法的精度较低,主要是其微分方程解 y 的 Taylor
5、 展开式所取的项数太少。显然为了提高计算精度,应当取泰勒公式(3)更高阶项。虽然增加高阶项可提高计算精度,但也同时带来了需要计算高阶导数的困难。龙格-库塔法的关键是利用低阶导数构成的曲线去拟合含有高阶导数的曲线,从而避免了计算高阶导数的问题。,7.1 仿真算法,7.1.3 数值积分的几个算法龙格-库塔法二阶龙格-库塔(RK)法 取(3)式的前三项,则有,设原微分方程(1)式解具有以下形式:,式中,a1,a2,b1,b2为待定系数。,将(8)式中K2按二元函数展开成泰勒级数,并取前三项,将K1,K2代入(8)式:,比较(10)、(7)式:,显然由(11)式并不能唯一确定a1,a2,b1,b2,因
6、为只有三个方程。因此对于同一种算法可以有不同的表现形式。,【说明】由于该算法只取到泰勒展开式的二阶导数项,所以称为二阶龙格-库塔法。但由(8)(12)式可知,算法并没有用y的二阶导数。,若设a1=a2,则,即二阶RK法公式为,7.1 仿真算法,7.1.3 数值积分的几个算法龙格-库塔法龙格-库塔(RK)法的一般形式,式中,i为待定权系数,ai,bij为待定系数,r为使用Ki的个数(即级数),Ki为所取各点导数f的值。Ki的个数与yk+1泰勒展开式所取的项数有关(即RK算法的阶数),同时还与计算区间内所取导数值的点数有关。,7.1 仿真算法,7.1.3 数值积分的几个算法龙格-库塔法四阶RK公式
7、 四阶RK公式用到了y的泰勒展开式的四阶导数。在RK算法的一般公式(13)中,取r=4可得:,由于(14)式在同级的RK算法中,计算精度较高,计算量较少,而在系统仿真的数值积分中应用十分广泛。称之为四阶四级RK公式。,7.1 仿真算法,7.1.3 数值积分的几个算法Gear算法“病态”常微分方程(刚性方程)的系数矩阵A的特征值具有如下特征:,则称为“病态”方程。,7.1 仿真算法,7.1.3 数值积分的几个算法Gear算法控制系统仿真中的“病态”问题病态系统中绝对值最大的特征值对应于系统动态性能解中瞬态分量衰减最快的部分,它反映了系统的动态响应和系统的反应灵敏度。一般与系统中具有最小时间常数T
8、min的环节有关,要求计算步长h取得很小。病态系统中绝对值最小的特征值对应于系统动态性能解中瞬态分量衰减最慢的部分,它决定了整个系统的动态过渡过程时间的长短。一般与系统中具有最大时间常数Tmax的环节有关,要求计算步长h取得很大。对于病态问题的仿真需要寻求更加合理的算法,以解决病态系统带来的选取计算步长与计算精度、计算时间之间的矛盾。,7.1 仿真算法,7.1.3 数值积分的几个算法Gear算法Gear算法Gear算法适用于病态系统的仿真,该算法类似于四阶RK算法,7.1 仿真算法,7.1.4 数值积分方法的选择 在选择积分方法时应考虑以下几个问题。计算精度 数值积分方法所得到的离散数值解只是
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 系统 时间 响应 及其 仿真
链接地址:https://www.31ppt.com/p-5373370.html