第5章 面向微分方程的数值积分法仿真 火电厂仿真课件.ppt

由此可见，简单的欧拉法是1阶精度；改进的欧拉法由于采用了平均斜率，相当于取了泰勒级数的前3项，因此为2阶精度。 分析欧拉法截断误差的思想，同样也适用于其它数值积分方法。 5、数值解的稳定性问题 显然，要使该差分方程是稳定的，必须使下式成立。 * * 第 5 章 面向微分方程的数值积分法仿真 数值积分是数值分析的一个基本问题。 也是复杂计算问题中的一个基本组成部分。 数值积分往往用极简单的方法就能较好地得出对所求解的具体数值问题的解答。 但数值积分的难点在于计算时间有时会过长，有时会出现数值不稳定现象。 另外，数值积分的理论性较强。其理论和方法都已经比较成熟，计算精度也比较高。 5.1 仿真中研究数值积分法的意义 数值解的一种近似方法。对于连续系统的高阶微分方程，可化为若干个一阶微分方程组成的方程组。 数值积分法是求解微分方程： 例如：下式所示的状态方程 可以化为一个一阶微分方程组 （5-1） （5-2） 所以，连续系统的仿真就是从给定的初始条件出发，对描述系统动态特性的常微分方程或常微分方程组进行求解，从而得到系统在一定输入作用下的变化过程。 在求解这些微分方程时，最常用、也是最有效的一种方法就是数值积分法。 （5-3） 5.2 数值积分法仿真的基本原理 对微分方程（5-1）两端同时取积分，可得 当 时，上式变为 ： （5-4） （5-5） 将积分项拆成两项 （5-6） 故上式可写为： （5-7） ——此式是方程（5-1）在tn+1时刻的精确解。 在数值解法中，希望用近似解 ： 代替准确解，其中 ： （5-8） 为 为 为 的近似值 令： ——称为计算步长或步距 式（5-8）是从常微分方程（5-1）出发建立的离散数学模型——差分方程。 由此可见，数值积分法就是在已知微分方程初值的情况下，求解该方程在一系列离散点 处的近似值，其特点是步进式——根据初始值逐步递推地计算出以后各时刻的值。 从式（5-8）可知，数值积分法的主要问题归结为如何对f（t，y）进行数值积分——求出f（t，y）在区间[tn， tn+1]上定积分的近似值Qn。 采用不同的方法求Qn，就出现了各种各样的数值积分方法。不同的数值积分将对求解的精度、速度和数值稳定性会产生不同的影响，这将在下述内容中具体介绍。 数值积分法种类繁多，在此从实用角度介绍几种基本的方法 5.3 欧拉(Euler)法 5.3.1 简单欧拉法 欧拉法是一种最简单的数值积分法，对于方程： 在区间[tn， tn+1]上求积分，得到： 若区间[tn， tn+1]足够小，则[tn， tn+1]上的f（t，y）可近似地看成常数f（tn，yn） 。故可用矩形面积近似代替 即： tn tn+1 f（tn，yn） 于是有： （5-9） 将此式写成差分方程为： （5-10） ——著名的欧拉公式 5.3.2 改进的欧拉法 如果用梯形面积而不是矩形面积来代替每一个小区间上的曲线积分，就可以提高计算精度，梯形法的计算公式为： （5-11） 式中的右端含有待求量yn+1，因而它是隐函数形式。这种方法不能自行启动运算，需要依赖其它算法的帮助。 每次计算都用欧拉法算出y（t n+1 ）的近似值 ，以此计算近似值 ，然后利用梯形公式（5-11）求出修正后的 。即有： 帮助方法： 预估式 校正式 （5-12） ——改进的欧拉公式 5.3.3 几个基本概念 简单的欧拉法是用前一时刻tn的yn求出后一时刻的yn+1，这种方法称为单步法，它是一种自行启动的算法。如果求yn+1时需要tn , tn-1 , tn-2 ……时刻yn , yn-1 , yn-2 ……的值，这种方法为多步法（改进的欧拉法为两步法），它是一种不能自行启动的算法。 1、单步法与多步法 简单的欧拉法表达式的右端，计算 所用的数据均已求出，这种公式称为显式公式。 改进的欧拉法表达式的右端，有待求量 ，这种公式称为隐式公式。 隐式公式不能自行启动，需要用预估-校正法。 单步法和显式在计算上比多步法和隐式方便，但有时为了满足精度、稳定性等方面的要求，需要采用隐式算法。 2、显式与隐式 3、截断误差 这里用泰勒级数为工具来分析数值积分公式的精度。 假定yn是精确的，用泰勒级数表示 处的精确解，即： 显然，简单的欧拉法是从以上精确解中取前两项之和来近似计算，每一步由这种方法引入的误差称为局部截断误差，简称截断误差。简单的欧拉法的截断误差为： 不同的数值方法有不同的截断误差。一般若截断误差为