数值积分法仿真幻灯片.pptVIP

第三章 数值积分法仿真 数值积分的基本性能 基本性能 数值积分算法的性能包含： 定性特征：稳定性 时间粒度：计算速度 空间粒度：计算精度 不同的数值积分方法的具有不同的稳定性。 同一个模型采用不同的积分算法和不同的积分步长h，稳定性不同。 计算速度和计算精度 各种数值积分方法的差分方程是对原微分方程的近似逼近，并且因为计算机的字长有限，存在明显的截断误差。 这些误差都和计算步距h密切相关，所以计算步距是影响计算精度、速度和稳定性的重要因素。 h取得较大，计算时间少，截断误差大； h取得较小，截断误差就会减小，但在给定时间范围内，计算次数必然增加，使误差积累增加。 截断误差、累计舍入误差与步长h 截断误差、累计舍入误差与步长h关系如图。 3.1 预报公式 3.2 校正公式 3.3 Adams公式 根据前面的分析，我们可以将预报和校正公式统一写成： 显式Adams系数 隐式Adams系数 3.4 多步法的特点 与单步法相比，相同精度下，使用过去多步信息，计算量小。 隐式法的精度高，稳定性好，但在计算y(n+k)时需要用到f[y(n+k),t(n+k)]，只能采用迭代法计算。 缺点之一是不能自启动，需用单步法计算初始值才能启动计算。 4.2 一般算法的稳定性分析 根据上例可得数值积分方法稳定域的一般方法。 设系统测试方程为： 而数值积分公式为： 只有当