系统仿真技术Chapter2经典的连续系统仿真建模方法学.ppt

控制系统仿真技术,盛 立 中国石油大学自动化系,Chapter 2 第2章经典的连续系统仿真建模方法学,对下面的控制系统描述，需要放在计算机上求解 常微分方程 传递函数 状态空间描述,方法一：ODE23,ODE45可解一阶微分方程组， 状态空间描述是一阶微分方程组 常微分方程，传递函数 状态空间表达式,求解？ ODE23,ODE45可解一阶微分方程组， 原理是什么？,对下面的控制系统描述，需要放在计算机上求解 常微分方程 传递函数 状态空间描述,方法一：ODE23,ODE45可解一阶微分方程组， 状态空间描述是一阶微分方程组 常微分方程，传递函数 状态空间表达式,原理: 一阶微分方程（线性，非线性）数值求解方法,数值求解方法,欧拉法,梯形法,龙格库塔法,RK2,RK4,2.1 离散化原理及要求,问题：数字计算机在数值及时间上的离散性-被仿真系统数值及时间上的连续性? 连续系统的仿真，从本质上：对原连续系统从时间、数值两个方面对原系统进行离散化并选择合适的数值计算方法来近似积分运算 离散模型原连续模型？,相似原理,设系统模型为： ，其中u(t)为输入变量，y(t)为系统变量；令仿真时间间隔为h，离散化后的输入变量为 ，系统变量为 ，其中 表示t=nh。 如果 , 且 即 , （对所有n=0,1,2,） ,则可认为两模型等价。,原连续模型,仿真模型,相似原理,对仿真建模方法三个基本要求,（1）稳定性：若原连续系统是稳定的，则离散化后得到的仿真模型也应是稳定的。 （2）准确性：有不同的准确性评价准则，最基本的准则是： 绝对误差准则： 相对误差准则： 其中 规定精度的误差量。,（3）快速性：若第k步计算对应的系统时间间隔为 计算机由计算需要的时间为 ，若 Tn=hn 称为实时仿真，Tnhn称为超实时仿真, Tnhn 称为亚实时仿真。,系统仿真中最常用、最基本的求解常微分方程数值解的方法主要是数值积分法。 设系统常微分方程为： （2-1） 为包含有时间t和函数y的表达式，y0为函数y在初始时刻t0时的对应初值。我们将求解方程（2-1）中函数 的问题称为常微分方程数值求解问题。,2.2 数值积分法,1欧拉公式的推导 将（2-1）式在小区间上进行积分可得：,其几何意义是把,在,区间内的曲边面积,用矩形面积近似代替，如图2-1所示。,2.2.1 欧拉法Euler,欧拉法Euler,当h很小时，可以认为造成的误差是允许的。所以有：,称之为欧拉公式。,截断误差正比于,欧拉Euler公式：,截断误差正比于,2 . 欧拉法具备以下特点： （1）欧拉法实际上是采用折线代替了实际曲线，也称之为折线法。 （2）欧拉法计算简单，容易实现。由前一点值仅一步递推就可以求出后一点值，所以称为单步法。 （3）欧拉法计算只要给定初始值，即可开始进行递推运算，不需要其它信息，因此它属于自启动模式。 （4）欧拉法是一种近似的处理，存在计算误差，所以系统的计算精度较低。,欧拉法Euler,数值求解方法,欧拉法,梯形法,龙格库塔法,RK2,RK4,截断误差正比于,1梯形公式 为了弥补欧拉法计算精度较低的不足，可以采用梯形面积公式来代替曲线下的定积分计算，如图2-2所示。 依然对式（2-1）进行求解，采用梯形法作相应近似处理之后，其输出为：,称为梯形积分公式 。,2.2.2 梯形法,梯形法,从中可以看到，在计算 时，其右端函数中也含有 ，这种公式称为隐式公式，不能靠自身解决，需要采用迭代方法来启动，称之为多步法。可以先采用欧拉公式进行预报，再利用梯形公式进行校正。即梯形法的预报校正公式 ：,梯形法,2. 梯形法具备以下特点： （1）采用梯形代替欧拉法的矩形来计算积分面积，其计算精度要高于欧拉法。 （2）采用预报校正公式，每求一个 ，计算量要比欧拉法多一倍。因此计算速度较慢。 （3）梯形公式中的右端函数含有未知数，不能直接计算左端的变量值，这是一种隐式处理，要利用迭代法求解。即梯形法不能自启动，要靠多步法来实现计算。,梯形法,数值求解方法,欧拉法,梯形法,龙格库塔法,RK2,RK4,截断误差正比于 ，记为,截断误差正比于 ，记为,欧拉Euler公式：,截断误差正比于,梯形法公式,截断误差正比于,2.2.3 龙格库塔法,1 龙格-库塔法基本原理 对 的数值求解：称作“右端函数”计算问题。 将 在 附近展开Taylor级数，只保留 项，则有：,若令：,则有,假设这个解可以写成如下形式： 其中 对 式右端的函数展成Taylor级数，保留h项，可得： 代入，则有：,龙格-库塔法基本原理（续）,进行比较，可得： 四个未知数 但只有三个方程，因此有无穷多个解。 若限定 ，则 计算公式： 其中,龙格-库塔法基本原理（续）,若写成一般递推形式，即为： 其中 截断误差正比于h3，称为二阶龙格-库塔法（简称RK-2）。,二阶龙格-库塔公式,四阶龙格库塔公式 ：,四阶龙格库塔（RungeKutta）法,截断误差正比于,（1）为单步法，并且可自启动。 （2）改变仿真步长比较方便，可根据精度要求而定。 （3）仿真计算量与仿真步长h的大小密切相关，h值越小计算精度越高，但所需仿真时间也就越长。 （4）用泰勒级数展开龙格库塔法计算公式时，只取h的一次项，即为欧拉法计算公式；若取到h2项，则为二阶龙格库塔法计算公式；若取到h4项，则为四阶龙格库塔法计算公式。,龙格库塔法特点,【例2.1】 已知一阶系统的微分方程为： ，初始条件 ，取仿真步长h=0.1，分别用欧拉法、梯形法和龙格库塔法计算该系统仿真第一步的值。,解：原方程可变为:,即,2.2.4 数值积分公式应用,（1）用欧拉法计算 根据欧拉公式，将函数表达式及其初始值代入后，可得该系统仿真第一步的值：,数值积分公式应用,（2）用梯形法计算： 根据预报校正公式，将函数表达式及其初始值代入后，可得仿真第一步的值。 用预报公式求起始值：,数值积分公式应用,再用校正公式得到系统仿真第一步的值：,数值积分公式应用,二阶龙格-库塔公式,（3）用二阶龙格库塔法计算 根据公式先计算出两个系数，再计算仿真第一步的值：,数值积分公式应用,则系统仿真第一步的值为：,数值积分公式应用,（4）用四阶龙格库塔公式计算 根据公式先计算出4个系数，再计算仿真第一步的值：,数值积分公式应用,四阶龙格库塔（RungeKutta）法,数值积分公式应用,则系统仿真第一步的值为：,数值积分公式应用,从上述结果可以看出: 对于同一个系统进行仿真计算时，其值的精度是随着数值积分公式的变化而改变的，其中欧拉法计算精度最低，其次为梯形法和二阶龙格库塔法，四阶龙格库塔法计算精度最高。,数值积分公式应用,数值积分公式在状态方程中应用,2.3.1 仿真精度与系统稳定性 1. 仿真过程的误差 （1）初始误差:现场采集数据不一定很准，会造成仿真过程中产生误差，称为初始误差。应对现场数据进行准确的检测，也可多次采集，以其平均值作为参考初始数据。 （2）舍入误差:由于不同档次的计算机其计算结果的有效值不一致，导致仿真过程出现舍入误差。 应选择挡次高的计算机，其字长越长，仿真数值结果尾数的舍入误差就越小。 （3）截断误差:仿真步距确定后，数值积分公式的阶次将导致系统仿真时产生截断误差，阶次越高，截断误差越小。仿真时多采用四阶龙格库塔法，其截断误差较小。,2.3 数值积分法性能分析,2. 仿真过程的稳定性 计算结果对系统仿真的计算误差反应不敏感，称之为算法稳定，否则称算法不稳定。对于不稳定的算法，误差会不断积累，最终可能导致仿真计算达不到系统要求而失败。 （1）系统的稳定性与仿真步长的关系 一个数值解是否稳定，取决于该系统微分方程的特征根是否满足稳定性要求，而不同的数值积分公式具有不同的稳定区域，在仿真时要保证稳定就要合理选择仿真步长，使微分方程的解处于稳定区域之中。,数值积分法性能分析,（2）积分步长的选择 由于积分步长直接与系统的仿真精度和稳定性密切相关，所以应合理地选择积分步长h的值。 通常遵循两个原则： 使仿真系统的算法稳定。 使仿真系统具备一定的计算精度。 一般掌握的原则是：在保证计算稳定性及计算精度的要求下，尽可能选较大的仿真步长。,数值积分法性能分析,由于工程系统的仿真处理采用四阶龙格库塔法居多，所以选择仿真积分步长可参考以下公式： 时域内： ；其中ts 为系统过渡过程调节时间 频域内： ；其中 为系统的开环截止频率,数值积分法性能分析,3.速度与精度 四阶方法的h可以比二阶方法的h大10倍，每步计算量仅比二阶方法大一倍 高于四阶的方法由于每步计算量将增加较多，而精度提高不快。,数值积分法性能分析,仿真步长与稳定性关系,2.4 稳定性分析,仿真方法选择的基本要求：仿真计算不改变原系统的绝对稳定性。 原系统是稳定的。观察欧拉法仿真递推公式 故有 (i) yn (n=0,1,2,)为它的一个仿真解，,稳定性分析（续）,设 为其准确解，即 （ii） 用(ii)式减去(i)式，可得： 即 特征方程为,稳定性分析（续）,特征方程为 显然，为了使扰动序列n不随n增加而增长，必须要求： 我们称它所对应的域就是该算法的稳定域： h 21/， 即h小于等于系统时间常数的两倍。,确定数值积分法稳定域的一般方法,测试方程： 数值积分公式 其中 是一个关于 高阶多项式函数， 则只有当 时，算法才稳定。,二阶