SS数值积分法仿真.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 一般算法的稳定性分析 根据上例可得数值积分方法稳定域的一般方法。 设系统测试方程为： 而数值积分公式为： 只有当 时，算法稳定。