计算机仿真技术第九次课.pptVIP

第九次课 3.4数值积分法稳定性分析 3.4.2 稳定性分析 3.5 数值积分法的选择与计算步距的确定 3.5.2 积分步长的确定 3.5.3 误差估计与步长控制 * 3.4.1 数值解法稳定性含义 考虑如下一阶系统: 其解析解为： 采用Euler法求其数值解。设计算步长为 则Euler递推公式为： 其时间常数是0.1 当 时， 递推结果显然是发散的； 数值解等幅振荡； (2) 当 递推结果是收敛的。 (3) 当 上述的临界步长 ,是系统时间常数的2倍。 时， 时， 以检验方程 为例进行讨论。 前向欧拉公式 要是上述差分方程稳定，必须使下式成立 或 前向欧拉法稳定区域 当 时，为了使计算稳定要求 ，即积分步长 必须小于系统时间常数的2倍。 （2） 后向欧拉公式 要是上述差分方程稳定，必须使下式成立 只要原系统稳定，此不等式必然成立。 因此，后向欧拉法是恒稳定的。 （3）梯形公式， 对于梯形公式，其稳定条件为 只要原系统稳定，此不等式必然成立。 因而，梯形数值积分公式也是恒稳定的。 RK(Runge-Kutta)法的稳定条件为 1-4 Adams 法的稳定区域 数值积分算法稳定性的结论是以检验方程为例得出的。在实际应用中，绝对稳定的积分算法，并不意味着计算结果不发散。它可以说明算法稳定性的大小或好坏。 3.5.1 积分方法的选择 1.精度要求 2.计算速