1. 引言
1.1微分方程是描述一个系统的状态随时间和空间变化的最基本的数学表达式,在物理,经济,工程等方面都有及其重要的作用。然而,微分方程大多数并不能直接求出解析式,因此在求解微分方程时更多采用数值计算。
1.2 数值求解微分方程的基本思路是将时间和空间离散化,然后将微分转化为差分,推导出差分递推式,从而利用边界条件和初始条件进行求解。
2 常微分方程的差分求解
一般的,一阶微分方程为
首先将连续变量y和t离散化,连续的函数y(t)和f(y,t)化为离散的
y n y_{_{n}} yn和 f ( y n , t n ) f(y_{n},t_{n}) f(yn,tn),上述微分方程可以化为差分方程
y n + 1 − y n Δ t = f ( y n , t n ) + f ( y n + 1 , t n + 1 ) 2 \frac{y_{n+1}-y_{n}}{\Delta t}=\frac{f(y_{n},t_{n})+f(y_{n+1},t_{n+1})}{2} Δtyn+1−yn=2f(yn,tn<