该博客探讨了一阶常微分方程的数值解法，包括二阶显式和隐式Adams公式以及Milne方法。通过实例比较，展示了三种方法在τ=0.2时的数值解，并与精确解和Euler方法进行对比。结果显示，二阶隐式Adams法和Milne方法的拟合度优于二阶显式Adams法，而Milne方法具有最高的精度。
一阶常微分方程的数值解法
这里我们介绍二阶显式、隐式 Adams 公式及 Milne 方法求解方程。
题目： 对初值问题
u ′ = u − t 2 , 0 ≤ t ≤ 1 , u ( 0 ) = 0 u^\prime=u-t^2,\ \ \ 0\le t\le1 ,u\left(0\right)=0 u′=u−t2, 0≤t≤1,u(0)=0
试分别用二阶显式，隐式 Adams 公式及 Milne 方法求数值解，取 τ = 0.2，并与精确解和 Euler 比较。
二阶显式的Adams 公式为：
u n + 1 = u n + τ 2 ( 3 f ( t n , u n ) − f ( t n − 1 , u n − 1 ) ) u_{n+1}=u_n+\frac{\tau}{2}\left(3f\left(t_n,u_n\right)-f\left(t_{n-1},u_{n-1}\right)\right) un+1=u
最低0.47元/天 解锁文章
扫一扫
2638
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
