matlab 前向欧拉法,前向后项差分和显式隐式欧拉法

摘要：

本文主要介绍前向后向差分，显式隐式欧拉法及其稳定性分析。

前向差分对应显式方法，后向差分对应隐式方法。显式欧拉法是比较流行的显式方法，隐式欧拉法是比较流行的隐式方法。

显式欧拉法条件稳定，对积分步长有要求，隐式欧拉法无条件稳定，对积分步长无要求(理论上如此，但实际使用中对积分步长仍然有要求，只是比显式欧拉宽松)。

1. 前向差分与后向差分

设一元函数函数离散化为一系列的点，，其中，, , 称为步长。

1.1 前向差分(Forward Difference)前向差分的定义：

之所以称之为前向差分，是因为当前时刻的差分，是下一时刻位置(时间向前前进一步)与当前时刻位置之差。

高阶前向差分

根据上式递归可得k阶前向差分公式为：

例如，二阶前向差分公式为：

1.2 后向差分(Backward Difference)后向差分定义：

由定义可以看出，当前时刻的差分，是当前时刻的位置与前一时刻的位置之差。

高阶后向差分

同上，k阶后向差分公式为：

因此，二阶后向差分公式为：

1.3 差商和导数差商就是差分除以步长，即，

在数值计算中，需要以差商代替导数，即，如果使用前向差分，则为显式方法，例如，

如果使用后向差分，则为隐式方法，例如，

2. 显式欧拉和隐式欧拉