求解微分方程的線性多步法

時(shí)間：2022-10-05 17:57:02 數(shù)學(xué)畢業(yè)論文我要投稿

相關(guān)推薦

求解微分方程的線性多步法

　　求解微分方程的線性多步法

　　摘要: 微分方程初值問題的求解一直是人們所關(guān)注的熱點(diǎn)問題,本文針對(duì)微分方程初值問題的線性多步法,利用它們的局部截?cái)嗾`差的主項(xiàng)系數(shù)的某種組合對(duì)它們進(jìn)行修正,提出了預(yù)測(cè)―校正算法,并給出了幾種常用的預(yù)測(cè)―校正算法。

　　關(guān)鍵詞: 微分方程初值問題線性多步法

　　1.微分方程初值問題的基本解法

　　=f(t,u)u(t)=u(1)

　　其中f為t和u的已知函數(shù),u為給定的初值。我們假設(shè)函數(shù)f(t,u)在區(qū)域:t≤t≤T,|u|<∞內(nèi)連續(xù),并且u滿足Lipschitz條件,即存在常數(shù)L,對(duì)所有t∈[t,T]和u,u,有|f(t,u)-f(t,u)|≤L|u-u|。在上述基本假設(shè)下,初值問題(1)在區(qū)間[t,T]上有唯一解,并且u(t)為連續(xù)可微的。

　　將區(qū)間[0,T]作N等分,小區(qū)間的長(zhǎng)度h=T/N稱為步長(zhǎng),點(diǎn)列t=nh(n=0,1,…,N)稱為節(jié)點(diǎn),t=0。由已知初值u(t)=u,可算出u(t)在t=t的導(dǎo)數(shù)值u′(t)=f(t,u(t))=f(t,u)。利用Taylor展式得u=u+hf(t,u),u就是u(t)的近似值。利用u又可算出u,如此下去可算出u在所有結(jié)點(diǎn)上的值,一般遞推公式為u=u+hf(t,u),n=0,1,…,N-1。

　　由文獻(xiàn)[1]知,求解該類方程分為單步法和多步法。只要利用h,t和u即可算出u的算法稱為單步方法,單步法包括Euler法、梯形法、Runge-kutta法,其一般形式為u=u+hφ(t,u;h),函數(shù)φ稱為增量函數(shù),它依賴于所給的微分方程。用單步法求出{φ},m=1,2,…,N只需要一個(gè)初值u。

　　需要用到h,t,t,…,t和h,u,u,…,u才能求出u(k>1)的算法稱為多步方法,線性多步法包括數(shù)值積分法和待定系數(shù)法,其一般形式為αu=hβf,這里α,β(j=0,1,K,k)是常數(shù),它需要有k個(gè)初值u,u,…,u才能得到整個(gè)序列{u},m=1,2,…,N。

　　2.線性多步法的改進(jìn)

　　常微分方程初值問題只給我們提供了一個(gè)初值,但線性k步方法需要k個(gè)初始值u,u,…,u,所以這其余的k-1個(gè)初始值就要通過其他方法來先期得到。最容易想到的是采用簡(jiǎn)單的單步方法,如Euler方法,或改進(jìn)的Euler法,但它們的收斂階較低。為保證整個(gè)計(jì)算的精度,我們應(yīng)使初始值的計(jì)算精度至少與所用的線性多步方法同階。

　　因此,若使用的線性多步方法為q階[2],[3]的,可以用q階Taylor展開方法或Runge-kutta方法來求u,u,…,u,然后進(jìn)入正式的求解過程。用Taylor展開確定初始值的同時(shí)還能得到關(guān)于選擇h的信息。如要求計(jì)算誤差不超過ε,則應(yīng)使

　　h|u(t)|≤εh|u(t)|>ε

　　同時(shí)成立。當(dāng)微商次數(shù)增大時(shí),計(jì)算量也將急速增大,因此利用低階微商構(gòu)造高精度的方法對(duì)實(shí)際計(jì)算是很有意義的,對(duì)u(t)=u(t)+?蘩u′(t)dt利用u′(t)在點(diǎn)t,t處的函數(shù)值及其各階導(dǎo)數(shù)值構(gòu)造Hernite型插值多項(xiàng)式,再代積分,舍去誤差項(xiàng)即可得到帶導(dǎo)數(shù)值的計(jì)算公式。

　　對(duì)隱式的計(jì)算格式一般采用迭代解法,將隱式線性多步方法改寫成

　　u-hβf(t,u)=-(αu-hβf)(2)

　　由于F是已知的,于是(2)式可用迭代法求解。

　　迭代法都需要有u,即u的迭代初值。選取一個(gè)好的u無疑是十分重要的,很自然的想法是用顯格式算出u――這稱為預(yù)測(cè),再用隱格式進(jìn)行迭代求u――這稱為校正,因此整個(gè)過程稱為預(yù)測(cè)―校正算法[4],簡(jiǎn)稱PC算法。顯然,預(yù)測(cè)的P式與校正的C式的階一般應(yīng)取成相同的。

　　預(yù)測(cè)―校正算法求u的過程一般如下:

　　P∶u=αu-hβf,

　　E∶f=f(t,u)C∶u=-αu+h(∑βf+βfs=1,2,…,n,

　　u=u,f=f(t,u)。

　　由于使用了同階的預(yù)測(cè)公式求得初始近似u,因此一般說來校正次數(shù)不會(huì)太多,大約兩三次就夠了。若校正步數(shù)太多則往往提示步長(zhǎng)過大,應(yīng)適當(dāng)減小步長(zhǎng)后再進(jìn)行計(jì)算。實(shí)際計(jì)算中經(jīng)常利用算式與算式同階的特點(diǎn),對(duì)(3)式再作一些修正。用線性組合的方法消去每個(gè)格式中的截?cái)嗾`差主項(xiàng)。僅花費(fèi)很小的代價(jià)就使方法的精度再提高一階,從而不需要再通過迭代來改善計(jì)算的近似值。

　　設(shè)預(yù)測(cè)算式和校正算式的局部截?cái)嗾`差分別為

　　Δ[u(t),h]=chu(t)+o(h)(3)

　　Δ[u(t),h]=chu(t)+o(h)(4)

　　整理后分別得到:

　　u(t)-[u+(u-u)]=o(h)(5)

　　u(t)-[u+(u-u)]=o(h)(6)

　　這兩個(gè)式子告訴我們,在求(u和u)后,再利用它們的局部截?cái)嗾`差[5],[6]的主項(xiàng)系數(shù)的某種組合對(duì)它們進(jìn)行修正,可使得修正后的局部截?cái)嗾`差比原來高一階,即達(dá)到o(h)。這里和被稱為修正系數(shù)。

　　直接按照(5)式對(duì)u進(jìn)行修正是不行的,因?yàn)榇藭r(shí)的u還沒有算出。但是,在計(jì)算u的近似值時(shí),u已經(jīng)算出來了,因此u和u也都已經(jīng)有了,這時(shí)可以用u-u去代替(5)中的u-u。

　　綜上所述,我們可以將帶修正的預(yù)測(cè)―校正算法的一般形式表示如下:

　　P∶u=-α+hβfM∶+(u-u)E∶=f(t,)C∶u=-α+hβ+hβfM∶u=u+(u-u)E∶f=f(t,u)(7)

　　在多步方法中,改善步長(zhǎng)將帶來一個(gè)新問題,就是結(jié)點(diǎn)變了。在預(yù)期結(jié)點(diǎn)處已經(jīng)算出的近似值,對(duì)于計(jì)算當(dāng)前結(jié)點(diǎn)來說,很多都用不上了,而在以新步長(zhǎng)為單位的結(jié)點(diǎn)處,往往并沒有計(jì)算過相應(yīng)的近似函數(shù)值。在使用的方法中,一般采用插值的方法去補(bǔ)上所需的這些點(diǎn)上的近似值再進(jìn)行計(jì)算。

　　下面給出幾種常用的帶修正的預(yù)測(cè)―校正算法。

　　(1)Milne預(yù)校算法。Hamming針對(duì)Milne算法絕對(duì)不穩(wěn)定,提出了如下修正:

　　P∶ρ(λ)=λ-1,σ(λ)=(2λ-λ-λ);

　　C∶ρ(λ)=λ-1,σ(λ)=(λ+4λ+1)。

　　P和C均為四階方法,M和M中的修正系數(shù)分別為和-。

　　(2)Hamming預(yù)校算法。該預(yù)校算法的P式和C式分別為:P∶與Milne預(yù)校算法相同。C∶ρ(λ)=λ-λ+,σ(λ)=(λ+2λ+λ)。

　　P和C也均為四階方法,M和M中的修正系數(shù)分別為和-。

　　(3)Adams四階預(yù)校算法。此算法中的P算式和C算式分別為:P∶Adams四步四階外插公式。P∶Adams三步四階內(nèi)插公式。

　　M和M中的修正系數(shù)分別為和-。理論研究表明,這個(gè)算法的穩(wěn)定性比前兩個(gè)算法好,是常用的預(yù)校算法。

　　3.結(jié)語

　　針對(duì)微分方程初值問題的線性多步法,我們利用局部截?cái)嗾`差的主項(xiàng)系數(shù)的某種組合對(duì)其進(jìn)行修正,使修正后的局部截?cái)嗾`差比原來高一階,即達(dá)到o(h),并給出了帶修正的預(yù)測(cè)―校正算法的一般形式,見文中(7)式。

　　參考文獻(xiàn):

　　[1]李立康,於崇華,朱政華編著.微分方程數(shù)值解法[M].復(fù)旦大學(xué)出版社,2005.

　　[2]胡健偉.微分方程數(shù)值方法[M].科學(xué)出版社,1999.

　　[3]李榮華.微分方程數(shù)值解法(第三版)[M].高等教育出版社,1996.

　　[4]Dennis G.Zill編著.微分方程與邊界值問題[M].機(jī)械工業(yè)出版社,2003.

　　[5]黃振侃.編著.數(shù)值計(jì)算――微分方程數(shù)值解[M].北京工業(yè)大學(xué)出版社,2006.

　　[6][美]茲爾著.陳啟宏等譯.微分方程與邊界值問題[M].華章數(shù)學(xué)譯叢.機(jī)械工業(yè)出版社,2005.

【求解微分方程的線性多步法】相關(guān)文章：

多邊形內(nèi)角的求解技巧10-26

簡(jiǎn)單的線性規(guī)劃教學(xué)教案10-08

乒乓球常采用的步法05-05

打羽毛球的基本步法10-05