- 22
- 0
- 约1.97千字
- 约 18页
- 2023-12-01 发布于上海
- 举报
解微分方程的欧拉法,龙格-库塔法简单实例比较
欧拉方法(Euler method)用以对给定初值的常微分方程(即初值问题)求解分为前EULER 法、后退 EULER 法、改进的 EULER 法。
缺点:
欧拉法简单地取切线的端点作为下一步的起点进行计算,当步数增多时,误差会因积累而越来越大。因此欧拉格式一般不用于实际计算。
改进欧拉格式(向前欧拉公式):
为提高精度,需要在欧拉格式的基础上进行改进。采用区间两端的斜率的平均值作为直线方程的斜率。改进欧拉法的精度为二阶。
算法:
微分方程的本质特征是方程中含有导数项,数值解法的第一步就是设法消除其导数值。对于常微分方程:
dy ? f (x, y) x ?[a, b]
dx
y(a) ? y
0
可以将区间[a, b] 分成n 段,那么方程在第x
i
点有 y (x
i
) ? f (x
i
, y(x
i
)) ,再用
向前差商近似代替导数则为:
( y(x
i
?1) ? y(x ))
i
h
? f (x
i
, y(x ))
i
在这里,h 是步长,即相邻两个结点间的距离。因此可以根据 x
i
点和 y 的数
i
值计算出 y 来:
i?1
y ? y
i ?1 i
h ? f (x , y
i i
) i ? 0,1,2,
L这就是向前欧拉公式。
L
改进的欧拉公式 :
将向前欧拉公式中的
原创力文档

文档评论(0)