转发：常微分方程初值问题数值解法

• 初值问题：即满足初值条件的常微分方程的解

首先，常微分 方程得有解-- -有解条件----利普希茨条件---

 

两类问题：

①单步法 ---计算下一个点的值  只需要用到前面一个点的值 

②多步法---计算下一个点的值  需要用到前面  个点的值 

1、欧拉法---下一个点的计算值等于前一个点的计算值加上步长乘以前一个点的函数值

①初值条件：给出因变量在某个点上的值

• 具体过程

第一步：将 初值条件带入微分方程，得到在该点的导数值
第二步：在该点，用taylor进行展开， 舍去二次项，将一次函数近似函数y
第三：计算 第二个点在直线上的值，用这个值近似函数  的在第二个点的值,依此类推，直到迭代完成

一些批注：显式欧拉方程指下一步要计算的值，不在迭代方程中；隐式欧拉方程指下一步要计算的值，在迭代方程中。

• 显示欧拉法与隐式欧拉法

评价： ① 精度为1（整体） ②误差阶为2（局部）

• 梯形方法----将显式欧拉迭代方程与隐式欧拉迭代方程做一下加权平均，构造的计算公式.

• 改进的欧拉方法---

思想：因为梯形公式是隐式公式，将显式欧拉公式对下一步的计算值进行预估，用梯形公式对下一步的计算值进行校正.

评价： ① 精度为2(整体) ②误差阶为3（局部）

2、龙格-库塔方法

 

• 2阶-龙格-库塔方法----类似改进的欧拉法

 

• 4阶-龙格-库塔方法

k1为区间左端点的斜率；

k2为区间中点的斜率，通过欧拉法，用斜率k1确定函数f在中点上的斜率k2

k3也是区间中点的斜率，通过欧拉法，用斜率k2确定函数f在中点上的斜率k3

k4是区间右端点的斜率，通过欧拉法，用斜率k3确定函数f在端点上的斜率k4

4阶-龙格-库塔方法---- 局部误差  ;  整体误差为