dt=0001 dt=9019 0.05 0.05 让我们看一下向前欧拉方法是如何解决刚性力的问题的。(见下一页的附图) 结果清楚地显示了同前面相同的稳定性趋势。由|λΔl≤2,我们得出稳定性 解(相应于△01002),但是我们同样可以由1<1a2(相应于 0001<Mt≤002)得出震荡解。 发生什么情况了?能够通过我们的直觉找到关于这个问题的一个可能 解吗?让我们回到一般的问题的标量方程: f(v,r),其中,v(0) dt 向前欧拉问题是: f(v",t") 这是一个外推法,它由t=时一外推得到v+ dt让我们看一下向前欧拉方法是如何解决刚性力的问题的。(见下一页的附图) 结果清楚地显示了同前面相同的稳定性趋势。由 λΔt ≤ 2,我们得出稳定性 解(相应于 Δt≤0.002 ),但是我们同样可以由 1< λ+t ≤ 2 (相应于 0.001 0.002 <Δ ≤t )得出震荡解。 发生什么情况了?能够通过我们的直觉找到关于这个问题的一个可能 解吗?让我们回到一般的问题的标量方程: (,) dv f v t dt = ,其中, 0 v v (0) = 向前欧拉问题是: 1 (,) n n v v n n f v t t + − = + 这是一个外推法,它由 时n t t = dv dt 外推得到 n 1 v + :