常微分方程的数值方法

线性多步法

张瑞
中国科学技术大学数学科学学院

线性多步法(Linear multistep method)

线性多步法

通常的$k$步线性多步格式可以表示为

\[\begin{aligned} y_{i+1}=&\alpha_0y_i+\alpha_1y_{i-1}+\cdots+\alpha_{k-1}y_{i-k+1} \\ &+h(\beta_{-1}f_{i+1}+\beta_0f_{i}+\cdots+\beta_{k-1}f_{i-k+1} \end{aligned} \]

其中$f_j=f(x_j,y_j)$, $\alpha$, $\beta$为系数。

当$\beta_{-1}\neq 0$时，格式是显式的，否则是隐式的。

它的局部截断误差为（假定$y_{i},\cdots,y_{i-k+1}$是精确的），

\[LTE=y(x_{i+1})-(\sum_{j=0}^{k-1}\alpha_jy(x_{i-j})+h\sum_{j=-1}^{k-1}\beta_jf(x_{i-j},y(x_{i-j}))) \]

利用$f(x_{i-j},y(x_{i-j}))=y'(x_{i-j})$，把上式在$x_i$处Taylor展开，可以得到误差表达式。

\[\begin{aligned} LTE=&y(x_i)+hy'(x_i)+\frac{h^2}{2!}y''(x_i)+\cdots \\ & -\sum_{j=0}^{k-1}\alpha_j[y(x_i)+(-jh)y'(x_i)+\frac{(-jh)^2}{2!}y''(x_i)+\cdots ]\\ & -h\sum_{j=-1}^{k-1}\beta_j[y'(x_i)+(-jh)y''(x_i)+\frac{(-jh)^2}{2!}y'''(x_i)+\cdots ] \end{aligned} \]

记

\[\begin{cases} c_0=&\displaystyle 1-\sum_{j=0}^{k-1}\alpha_j \\ c_m=&\displaystyle\frac1{m!}-\frac1{m!}\sum_{j=0}^{k-1}(-j)^m\alpha_j-\frac1{(m-1)!}\sum_{j=-1}^{k-1}(-j)^{m-1}\beta_j , m=1,\cdots \end{cases} \]

待定系数法构造

命题 .当$c_i=0, i=0,1,\cdots,p$, $c_{p+1}\neq 0$时，格式的局部截断误差是$O(h^{p+1})$，是$p$阶格式。

要构造$p$阶的格式，取合适的$\alpha_i$, $\beta_i$，使得$c_i=0, i=0,1,\cdots,p$

例 1. 当k=1时，可以得到条件

\[\alpha_0=1, \beta_{-1}+\beta_0=1 \]
  • $\beta_{-1}=0$时，可以得到向前Euler方法 $y_{j+1}=y_j+hf_i$
  • $\beta_{-1}=1$时，可以得到向后Euler方法 $y_{j+1}=y_j+hf_{i+1}$
  • $\beta_{-1}=\beta_0=\frac12$时，可以得到梯形公式 $y_{j+1}=y_j+\frac{h}2(f_{i+1}+f_j)$

例 2. 推导形如

\[\begin{aligned} y_{i+1}=&\alpha_0y_i+\alpha_1y_{i-1}+\alpha_2y_{i-2} \\ &+h(\beta_{-1}f_{i+1}+\beta_0f_i+\beta_1f_{i-1}+\beta_2f_{i-2}) \end{aligned} \]

的隐式线性三步格式，确定系数，使格式为4阶方法。

解. 列出方程为

\[\begin{array}{l} \alpha_0+\alpha_1+\alpha_2=1 \\ \alpha_1+2\alpha_2-\beta_{-1}-\beta_0-\beta_1-\beta_2=-1 \\ \frac12\alpha_1+2\alpha_2+\beta_{-1}-\beta_1-2\beta_2=\frac12 \\ \frac16\alpha_1+\frac43\alpha_2-\frac12\beta_{-1}-\frac12-\beta_1-2\beta_2=-\frac16 \\ \frac1{24}\alpha_1+\frac23\alpha_2+\frac16\beta_{-1}-\frac16\beta_1-\frac43\beta_2=\frac1{24} \\ \frac1{120}\alpha_1+\frac4{15}\alpha_2-\frac1{24}\beta_{-1}-\frac1{24}\beta_1-\frac23\beta_2\neq-\frac1{120} \end{array} \]

7个未知数，5个方程。有无穷组解。

Simpson格式

\[y_{i+1}=y_{i-1}+\frac{h}3(f_{i+1}+4f_i+f_{i-1}) \]

Hamming格式

\[y_{i+1}=\frac18(8y_i-y_{i-2})+\frac{3h}8(f_{i+1}+2f_i-f_{i-1}) \]

4阶Adams隐格式

\[y_{i+1}=y_i+\frac{h}{24}(9f_{i+1}+19f_i-5f_{i-1}+f_{i-2}) \]
\[\begin{aligned} &\alpha_0+\alpha_1+\alpha_2=1 \\ &\alpha_1+2\alpha_2-\beta_{-1}-\beta_0-\beta_1-\beta_2=-1 \\ &\frac12\alpha_1+2\alpha_2+\beta_{-1}-\beta_1-2\beta_2=\frac12 \\ &\frac16\alpha_1+\frac43\alpha_2-\frac12\beta_{-1}-\frac12-\beta_1-2\beta_2=-\frac16 \\ &\frac1{24}\alpha_1+\frac23\alpha_2+\frac16\beta_{-1}-\frac16\beta_1-\frac43\beta_2=\frac1{24} \\ &\frac1{120}\alpha_1+\frac4{15}\alpha_2-\frac1{24}\beta_{-1}-\frac1{24}\beta_1-\frac23\beta_2\neq-\frac1{120} \end{aligned} \]
\[\begin{array}{l} \alpha_0+\alpha_1+\alpha_2=1 \\ \alpha_1+2\alpha_2-\beta_{-1}-\beta_0-\beta_1-\beta_2=-1 \\ \frac12\alpha_1+2\alpha_2+\beta_{-1}-\beta_1-2\beta_2=\frac12 \\ \frac16\alpha_1+\frac43\alpha_2-\frac12\beta_{-1}-\frac12-\beta_1-2\beta_2=-\frac16 \\ \frac1{24}\alpha_1+\frac23\alpha_2+\frac16\beta_{-1}-\frac16\beta_1-\frac43\beta_2=\frac1{24} \\ \frac1{120}\alpha_1+\frac4{15}\alpha_2-\frac1{24}\beta_{-1}-\frac1{24}\beta_1-\frac23\beta_2\neq-\frac1{120} \end{array} \]

7个未知数，5个方程。有无穷组解。

基于数值积分的构造方法

线性多步法的相容性

当线性多步格式

\[\begin{aligned} y_{i+1}=&\alpha_0y_i+\alpha_1y_{i-1}+\cdots+\alpha_{k-1}y_{i-k+1} \\ &+h(\beta_{-1}f_{i+1}+\beta_0f_{i}+\cdots+\beta_{k-1}f_{i-k+1} \end{aligned} \]

的局部截断误差是步长$h$的高阶无穷小量时，格式是相容（Consistency）的。利用Taylor展开，可以证明

定理 1.
线性多步格式是相容的充要条件是

\[\sum_{j=0}^{k-1}\alpha_{j}=1 , \sum_{j=-1}^{k-1}\beta_j=k-\sum_{j=0}^{k-1}(k-1-j)\alpha_j \]

线性多步法的稳定性

通常的$k$步线性多步格式

\[\begin{aligned} y_{i+1}=&\alpha_0y_i+\alpha_1y_{i-1}+\cdots+\alpha_{k-1}y_{i-k+1} \\ &+h(\beta_{-1}f_{i+1}+\beta_0f_{i}+\cdots+\beta_{k-1}f_{i-k+1} \end{aligned} \]

称$P(\lambda)=\lambda^k-\alpha_0\lambda^{k-1}-\alpha_1\lambda^{k-2}\cdots-\alpha_{k-1}$ 为该$k$步法的特征多项式。设$\lambda_1$, $\lambda_2$, $\cdots$, $\lambda_k$是特征方程$P(\lambda)=0$的根。 则有

定理 2.
$k$步法格式稳定的充要条件是，$|\lambda_i|\leq 1, \forall i=1,2,\cdots,k$，且模为$1$的都是单根。

目录

本节读完

例 3.

[#ex9-1-0].