正向欧拉法的稳定性

计算科学 稳定 显式方法
2021-12-23 06:36:13

我试图了解正向欧拉方法的稳定性。我读到有一个模型问题来查看稳定性。

y(t)=λy(t)t(0,)
y(0)=1

然后这本书显示了这一点:

u0=1does this come from the initial condition?

然后

un+1=un(1+λh)=(1+λh)n+1,n0

我无法理解是如何出现的(1+λh)n+1

编辑

然后分析表明的唯一方法是|1+\lambda h| < 1(为什么它需要小于 1?)这给出(根据书) h < \frac{2}{|\lambda|}

limnun=0
|1+λh|<1
h<2|λ|

对于不等式|1+λh|<1我得到的是\frac{-2}{\lambda} < h但仍然不知道他们如何获得\lambda2λ<h的绝对值λ

这本书是如何得出这个结论的?我不希望得到一个非常复杂的答案(值得赞赏),而是一种简单的方法来理解未显示的中间步骤。

1个回答

好的。所以 , , , , ..., , 这将回答您问题的一部分。我不明白你问题的另一部分。为了评估稳定性,我们假设您可以稍后自己考虑其他可能性。微分方程的真解是当 t 趋于无穷大时,解趋于零。离散方程也必须如此。时,我们的离散方程的解将变为零所以u0=1u1=(1+λh)u2=(1+λh)2u3=(1+λh)3un=(1+λh)nun+1=(1+λh)n+1λ<0u0eλt(1+λh)>1h<2/λ=2/|λ|.

另一个提示:当你将小于一的东西与自身相乘时,你会得到小于一的东西。当它变大时,它会增长。

其它你可能感兴趣的问题