计算机仿真技术第三章86202.pptVIP

计算机仿真技术第三章 连续系统数值积分 主 要 内 容 3.1 连续系统数值积分法的基本原理 3.2 Runge-Kutta方法 3.3 线性多步法 3.4 数值积分法稳定性分析 3.5 数值积分法的选择与计算步距的确定 3.6 离散相似法 3.7 替换法 绪言 3.1 连续系统数值积分法基本原理 3.1.1 Euler法 3.1.2 改进的Euler法 * *研究各类系统（如天体、水流、飞行器等）的动态特性过程中，从理论上对其进行描述的基本手段就是数学模型，它们都可以转换成微分方程或状态空间方程，也可以由这两类模型导出对应的其它模型，如传递函数、权函数等等，将它们进行离散化处理又可得到离散化模型。其实建立系统数学模型的目的就是研究系统的运动规律，只有对建立的数学模型求解，并且将结果与实际工程的要求进行分析研究才能够达到掌握系统运行规律的目的。因此如何把已建立起来的数学模型转换成仿真运算模型，以便为分析解决实际问题服务，是系统仿真学科的一个重要研究内容。对于复杂系统的数学模型来说，求其解析解往往很烦琐、很困难，大多数情况求不出解析解，或者根本不存在解析解，因此必须借助于数值解法对连续系统进行仿真研究。一般情况下，实际工程中对连续系统的数学模型进行直接求解析解是一件很困难的事情。数值积分法（或称为数值解法），就是利用数值积分方法对较常见的常微分方程（组）建立离散化形式的数学模型 — 差分方程，并求其数值解。现代计算机提供了强大的数值计算功能，所以只要从算法上研究出了工程上可以接受的针对这些模型的计算理论，就可以实现对连续系统任何数学模型的数值求解。从形式上观察不难看出，连续系统常见的三类数学模型中最为简单的是形式为一阶微分方程的状态方程数学模型。我们先不考虑其系数为矩阵的情况，只把它作为一个简单的一阶微分方程看待。然后在有了求解算法的基础上再研究系数为矩阵时候的计算机算法设一阶向量微分方程及初值(3-1)其中 ，为 维状态向量和 维向量函数，积分得将积分区间进行分割，且设在每处的连续解为对于上面积分项不容易找到原函数，所以希望能找到近似公式上面为的近似值。所谓数值积分解法，就是寻求一级微分方程初值问题的真解在一系列离散点上的近似解（即数值解），相邻两个时间离散点的间隔称为计算步长，通常取（定值）。可见，数值积分法的主要问题归结为对函数的数值积分问题，即如何求出该函数定积分的近似解。为此，首先要把连续变量问题用数值积分方法转化为离散的差分方程的初值问题，然后根据已知的初始条件，逐步地递推计算后刻的数值解。所以，解初值问题的数值方法的共同特点是步进式的，采用不同的递推算法，就出现各种不同的数值积分方法。基本的方法有三类：单步法，多步法和预估-校正法。递推公式又有显式和隐式之分。采用不同的数值积分公式，对系统的求解精度、速度和稳定性均有不同的影响。Eu1er法是最简单的一种数值积分方法，也是所有数值积分方法的理论基础，虽然它的计算精度比较差，实际中也很少采用，但由于它导出简单，几何意义明显，便于理解，能说明构造数值解法一般计算方法的基本思想，所以通常用它引入所有这些理论研究中的基本概念。为了使读者对数值积分算法的基本思想有个透彻的了解，下面将从三个侧面推导出Euler法数值解近似公式。1. Taylor级数展开2. 矩形近似解法3. 切线近似法1. Taylor级数展开对于标量微分方程， 假定其解析解为 在 附近将其展成Taylor级数的形式 当 适当小时，可以近似地表示成 将上式写成差分方程2. JUxingfa对于长度 为的区间，在其上对前述微分方程求积分得 式子后面的积分是曲线下面的面积， 如果区间足够小，则可以用矩形面积 近似代替这个面积，即函数可 以近似地看成常数，即或写成差分方程 的形式3. 切线近似法在 的小邻域内，曲线可以用 处的切线表示， 处的斜率在 处的切线方程 取，使与 十分靠近， 利用切线方程可得处的近似 值 取，则可以写成差分方程的形式 例. 设系统方程，。用Euler法求数值解() 解：原方程化为，其中 一方面求得解析解 另一方面利用Euler公式得 数值解： 方程的解析解和数值解对照 0.4628 ... 0.647 0.6594 0.7519 0.819 0.9 1 数值解yn 0.5 ... 0.625 0.6667 0.7692 0.8333 0.9091 1 精确解y(t) 1.0 ... 0.5 0.4 0.3 0.2 .1 0 t用梯形面积代替每一小区间曲面面积，可提高精度。如图积分为曲边梯形的面积 直边梯形的面积为比较小时，直边梯形面积 近似 等于曲边梯形的面积 直接取值，得差分方程 或简记为 讨论前面式