3.1连续系统数值积分法基本原理数值积分法-Read.pptVIP

第三章 连续系统数值积分仿真方法学 　 　本章主要内容： 连续系统数值积分方法的基本原理 有关积分方法的一些基本概念 几种最常用的数值积分方法 3.1 连续系统数值积分法基本原理 数值积分法： 　 连续系统的数值积分法（或称为数值解法），就是利用数值积分方法对常微分方程（组）建立离散化形式的数学模型－差分方程，并求其数值解。 　　特点： 在数字计算机上构造数字积分器的方法就是数值积分法，而且数字机的硬件特点决定了这种积分运算： 离散的、串行的、近似的、递推的。 数值积分的基本原理： 若一阶向量微分方程及初值为 　为　维状态向量，　　为　维向量函数。 设方程在 　　 　　　　处的形式上的连续解为： 希望能找到一个近似公式 为准确解 的近似值 是准确积分值的近似值 。 所谓数值解法，就是寻求初值问题的真在 一系列离散点 上的近似解 。 相邻两个时间离散点的间隔 称为计算步长， 积分方法分类： 基本的方法有三类： 单步法，多步法和预估－校正法。 递推公式有显式和隐式之分。 采用不同的数值积分公式，对系统的求解精度、速度和稳定性均有不同的影响。 Euler法 特点：导出简单，几何意义明显，便于理解，能说明构造数值解法一般计算公式的基本思想。通常用它来说明有关的基本概念。 从三个侧面出Euler方法数值解近似公式。 Taylor 级数展开 考虑如下标量微分方程 假定 为其解析解。将 展成Taylor级数 从而 将上式写成差分方程 矩形近似解法 对(3.5)在区间 上求积分，得 用矩形面积近似地代替该区间上的曲线面积 于是 即 此方法的实质：将曲线 看成阶梯函数。 切线近似 在 的一个小邻域 内，曲线 可以用 处的切线来表示， 在 处的斜率为： 则 在 处的切线方程为 Euler法的优点是方法简单，计算量小，缺点是精度低。 例3.1 设系统方程为 用Euler法求其数值解（取步长 , ) 解：先将原方程变换为状态方程形式为 则 递推公式为 已知方程的解析解为 　　　 精确解和解析解作比较： 误差在 数量级，精度较差。 改进Euler法 　　梯形面积代替每一个小区间的曲线面积，则可提高精度。 根据定积分的几何意义，积分为曲边梯形的面积 　 直边梯形的面积为 　特点： 上式称为梯形公式或二阶隐式Adams公式。右端是隐函数形式，故梯形法不能自行起步运算，而需有其它算法的帮助。 每次计算可用Euler法求出初值. 步骤： 1用Euler法算出 的近似值 　 2 将 代入原微分方程右端，计算导函数 　 的近似值 3最后利用梯形公式求出修正后的 。 改进的Euler法为预估－校正法: 第一式称为预估公式，第二式称为校正公式。 改进Euler公式的思想 改进Euler公式实质上是采用了 和 两点斜率平均值的结果。 ? 在实际计算中通常采取加权平均的方法，在每个积分步中取几个点，分别求出其斜率 , 取不同的权值，公式为 数值积分法的几个基本概念 算法自启动 有一些数值积分方法只要知道了方程的初值 ,便可以利用算法公式求得后续一系列点 上的近似解，而无需其它算法的辅助。 具有这种性质的数值积分方法可以自启动，如Eul