孟波，孟純青

（1．聯(lián)科應用研究所，山東 濟 南250031；2．山東省濟南中學，山東 濟 南250001）

0 引言與基本知識介紹

提出一種構建常微分方程數(shù)值方法的新思路，按照這種思路，構建一些具體的計算公式．先介紹一階和二階偏差分．

對于函數(shù)f（x，y）在點（xi，yi）處的的一階形式如下．

對于二階偏差分有以下形式．

最后給出兩個變量的兩種Taylor公式．

1 新方法的基本思想與相關公式的推導

設存在常微分方程初值問題f（x，y）關于y滿足Lipschitz條件，即對于y1、y2有｜f（x，y1）－f（x，y2）｜≤L｜y1－y2｜．設步長為h，這時有xi＋1＝xi＋h，xi－1＝xi－h(huán)．現(xiàn)在介紹新方法的基本思想．

如果假設yi＋1－yi＝ki，忽略余項，（1）式和（2）式可以寫成以下形式：

把（6）～（7）式跟一些傳統(tǒng)方法結合起來就得到新方法中的公式，這些公式的形式可參見文獻［1－3］．

隱式Euler方法的形式為：

如果把右邊的f（xi＋1，yi＋1）按照（6）～（7）式展開后得到

公式（A1）、（A2）就是兩種n階新方法中的公式，其中（A1）是按照（6）式推導出來的，稱為新方法中的單差分公式，（A2）式是按照（7）式推導出來的，稱為雙差分公式．選擇具體的偏差分形式可以構建不同的計算公式．下面就是一些具體公式．

如果只考慮到一階偏差分，對于（A1）式，設0≤r≤1，記：

有計算公式：

對于（A2）式，設0≤r1，r2≤1，記：

有計算公式：

設0＜s＜1，把以上思想用于加權梯形公式得到：

如果只考慮到一階偏差分，把（8）～（10）式代入得到以下公式：

如果只計算到二階偏差分，公式（A4）的一種具體形式為：

把以上思想和其他公式相結合，可以得到許多新方法中的公式．

與Adams隱式方法相結合可以得到許多公式，以下僅列出從三階、四階Adams隱式方法推導出來的新公式．

如果計算到二階偏差分，公式（B4）的一種具體形式為

從Gear方法中可以推出許多新方法中的公式．以下給出三步新方法中的公式

如果只計算到二階偏差分，公式（C2）的一種具體形式為：

用以上思想可以構造出許多新方法中的公式，這類公式數(shù)量極多，本文中不再一一列舉．

2 新方法中的加權平均公式和預估校正公式

按照筆者的思路，可以構造很多具體的加權平均公式和預估校正公式，限于篇幅分別給出一例．

設存在r1，r2，滿足0＜r1，r2＜1，則（A1）式和（A2）式的加權平均公式為：

使用具體的r1，r2，構建具體的加權平均公式，按照以上思路，可以構造出許多加權平均公式．

新方法中的顯式公式可以跟隱式公式相結合形成預估校正公式（A4．1）中一個的具體公式和梯形公式組成的預估校正公式為：

3 pi、qi值的確定

以上算法中出現(xiàn)了pi、qi兩個參數(shù)，現(xiàn)在來討論pi、qi的計算方法，計算方法可以有多種形式．

如果假設?xi＝dxi＝h，即dyi＝hf（xi，yi），則

如果存在常數(shù)r1，r2，滿足?yi＝r1（yi＋1－yi）＋r2（yi－yi－1），則：

對于qi則可以有以下選擇．

如果假設?yi＝y(tǒng)i＋1－yi，那么

如果假設?yi＝y(tǒng)i－yi－1，那么

如果存在常數(shù)r1，r2，滿足?yi＝r1（yi＋1－yi）＋r2（yi－yi－1），那么：

由于（E2）式、（E4）式、（E5）式出現(xiàn)了yi＋1，這就變成了隱式公式．為便于計算，可用一些近似方法．

一是直接計算方法：把某種顯式公式代入直接給出yi＋1的表達式；二是采用預估校正的方法得到相應的計算公式．以下給出幾個具體公式．

直接計算公式．如果假設yi＋1＝y(tǒng)i＋hf（xi，yi），則代入（E2）式、（E4）式得到

代入（E5）式得到

預估校正公式有很多，下面介紹兩種．

一階單差分預估校正公式：

二階雙差分預估校正公式：

直接計算公式和預估校正公式還可以有其它形式，限于篇幅本文只給出以上形式的公式．

4 遞進計算公式

實際計算時，為了提高精度，可以采用遞進計算公式，計算步驟如下：

對于顯式公式，當計算第i＋1步的數(shù)值時，可以使用第i，i－1，i－2，…，i－j步時的數(shù)據(jù)，這樣可以按照實際的計算步數(shù)來確定計算公式的階數(shù)．當計算第一步數(shù)據(jù)時，j＝1，當計算第二步數(shù)據(jù)時j＝2，當計算第k步數(shù)據(jù)時，j最多可以等于k，這樣層層遞進從而得到一個遞進計算公式．遞進計算公式實際上是一個變階計算公式．

5 相容性、收斂性和穩(wěn)定性分析

以下討論相容性，收斂性和穩(wěn)定性．這涉及到qi的表達式的選擇，由于yi＋1可以預估而yi，yi－1已知，這樣每一個qi也可以得到一個預估值．所以下面的討論將每一個qi當作常數(shù)．

先討論相容性．限于篇幅，僅討論（A2．1）式，為了簡便起見，設，對于（A2．1）式，取，顯然當h→0時，均有偏差分?xf（xi，yi）→0，?yf（xi，yi）→0，于是ψ（x，y，0）＝f（x，y）．因此方法具有相容性．

下面討論收斂性．對于（A2．1）式，為了簡便起見，設r1＝0，r2＝0．對于y1，y2，有

設｜y1－δ－y2－δ｜＝w｜y1－y2｜，則上式滿足

因此ψ（x，y，h）關于y滿足Lipschitz條件，再考慮到相容性，故算法是收斂的．可用類似的方法討論其他新方法的相容性和收斂性．

最后討論穩(wěn)定性．限于篇幅以下僅討論（A3）式和（A4）式的穩(wěn)定性．設存在以下實驗方程：

代入（A3）式或（A4）式，這時兩個公式一樣．現(xiàn)給出代入（A4）式的情形．設s＝0．5．

限于篇幅，以下僅討論到n＝1時的情況．如果取?yyi＝?yi＝y(tǒng)i－yi－1，qi＝1代入（11）式有：

移項后得到：

相應的特征方程為：

如果取?yyi＝?yi＝y(tǒng)i＋1－yi，代入（11）式有：

移項后得到：

6 實例計算

例：求解下列常微分方程初值問題，y′＝x－y＋1，y0＝1，h＝0．1，0≤x≤1．

本例有精確解y＝x＋e－x，其傳統(tǒng)方法的計算結果可見參考文獻［4］．為了比較各方法的精確度，以下采用多種方法進行求解，以便比較各種方法的精確度，同一個公式使用pi＝1，而使用不同的qi分別計算，以便觀察計算精度．

現(xiàn)把一些相關的計算公式列舉如下：公式（B4．1），公式（C2．1），（A4．2）．以上公式分別根據(jù)不同的qi的值或者表達式進行計算．

表2給出了公式（B4．1），公式（C2．1）的計算結果，其中前3個公式為（B4．1）的計算結果，后3個公式為（C2．1）的計算結果．6個公式的參數(shù)分別為：

在計算中由于一些公式是高階公式，所以前兩個或3個值使用精確解的數(shù)值．前3個精確值為：y0＝1，y1＝1．004 837 418，y2＝1．018 730 753．計算結果如下：

表1 傳統(tǒng)公式計算結果及精確解

表2 公式（B4．1）、（C2．1）的計算結果

表3 公式（A4．2）的計算結果

從以上計算結果看出：以上公式的計算結果明顯好于隱式Euler公式，（C2．1）式的計算結果與兩步Gear公式相當，精度高低跟選擇的具體參數(shù)有關，（B4．1）式的3個公式計算結果比隱式Adams公式略差，（A4．2）式的計算結果跟梯形公式大致相當，精度的高低跟具體的參數(shù)選擇有關，第5個公式的計算結果跟精確解非常接近．

新差分方法具有以下特點：可以把隱式公式變成顯式公式，這大大減少了計算難度；二是可根據(jù)不同的差分形式得到種類繁多的具體計算公式．可以根據(jù)具體情況酌情使用，這對實際計算非常有利；三是可以構造出種類繁多的加權平均公式．筆者認為，新方法是一種值得認真研究的常微分方程數(shù)值方法，具有極大的理論價值和應用價值．

［1］李榮華，劉播．微分方程數(shù)值解法［M］．4版．北京：高等教育出版社，2009：1－40．

［2］戴嘉尊．微分方程數(shù)值解法［M］．南京：東南大學出版社，2012：1－45．

［3］倪興．常微分方程數(shù)值解法及其應用［D］．合肥：中國科學技術大學，2010．

［4］王立秋，魏煥彩，周學圣．工程數(shù)值分析［M］．濟南：山東大學出版社，2002：350－370．