数值积分法稳定性分析

3.4数值积分法稳定性分析3.4.1数值解法稳定性含义考虑如下一阶系统采用Euler法求其数值解。设计算步长为,则Euler递推公式为(1)当时，，递推结果发散；(2)当时，数值解显等幅振荡趋势；(3)当时，递推结果收敛。10,010yyyh101110110110110nnnnnnhyyhyhyy　　　2.0h2.0h2.00h1101h所谓数值解的稳定性:指在扰动(初始误差、舍入误差、截断误差等)影响下，其计算过程中的累积误差不会随计算步数的增加而无限增增长。判断：不同的数值解法对应着不同的差分递推公式。一个数值法是否稳定取决于该差分方程的特征根是否满足稳定性要求。3.4.2稳定性分析以Euler法为例说明各种数值积分方法稳定性分析方法。Euler公式有以下三种形式：(1)前差公式(2)后差公式(3)梯形公式nnnhfyy111nnnhfyy112nnnnffhyy以检验方程为例进行稳定性讨论，。(1)前差公式为要使上述差分方程稳定，必须使当系统有实根时，为了保证计算稳定性，要求结论：步长必须小于系统时间常数的两倍。yy0,jnnnnyhhyyy1111h0，2hh(2)后差公式为差分方程的特征根结论：只要原方程稳定，那么利用后差公式获得的差分方程的特征根一定落在单位圆内，与步长无关。后差公式是恒稳定的。01111nnnnnyyhyhyy11111222hhhz(3)对于梯形公式，其差分方程特征根为也是恒稳定的。0122122122222，　　hhhhz此思想，也适用于其他数值积分方法。类似地可得RK法的绝对稳定域据此可得出各类RK公式的稳定条件。1!1!211:2rrhrhhhG表3.4RK方法的稳定区域表3.5Adams法的稳定域rrG1G2G3G4G0,20,20,51.20,78.20,20,10,1160,1030,0,0,60,3条件稳定算法，步长必须满足下列不等式其中为由积分方法确定的常数。相当于连续系统微分方程或状态方程的特征根或闭环系统的极点。MhhM