计算机仿真技术PPT件第三章 连续系统数值积分仿真方法学.ppt

只用三个方程 二阶二步显式算法的一般形式： 二阶二步Adams显式算法 为Adams法 为向后微分公式 显然是隐式公式 若要求向后微分公式为K阶算法，则 当K=1 为向后Eular法 K=1，2，….6的向后微分公式系数值 1 1 -1 1 2 2/3 1/3 -4/3 1 3 6/11 -2/11 9/11 -18/11 1 4 12/25 3/25 -16/25 36/25 -48/25 1 5 1 6 1 k 计算 时，仅已知 和相应的导数值 称为线性多步法的起步值，需用单步法来求出。 若线性多步法是P阶算法，计算起步值的算法不应低于P阶，否则影响计算精度。 仿真步长 h 的选取，是否影响仿真结果？先看一个例子： 第四节 数值积分法稳定性分析 系统稳定，精确解： 用前向欧拉法求解： ch3_6.m 由此可见，仿真步长 h 的选取，会影响仿真结果。 用前向欧拉法求解，当 h ≧ 0.2时不稳定，是由于步长太大，从而截断误差太大造成的。 一、数值解法稳定性的含义 数值解的稳定性：在扰动（初始误差、舍入误差、截断误差）的影响下，计算过程中的累积误差不会随计算步数的增加而无限增长。 微分方程的数值积分方法，实质是微分方程的差分化，然后从初始条件递推迭代。不同数值解法对应着不同的差分方程，是否稳定取决于