數(shù)值分析線性多步法

第五節(jié)線性多步法

龍格-庫塔方法是一類重要算法，但這類算法在每一步都需要先預(yù)報幾個點上的斜率值，計算量比較大?？紤]到計算yi+1之前已得出一系列節(jié)點上的斜率值，能否利用這些已知值來減少計算量呢？這就是線性多步法的設(shè)計思想。

構(gòu)造線性多步法有多種途徑，這里介紹兩種：

基于數(shù)值積分的構(gòu)造方法；

基于Taylor展開的構(gòu)造方法。)...(...110111101rnrnnnrnrnnnffffhyyyy--+---+++++++++=bbbbaaa線性多步法的通式可寫為：當(dāng)10時，為隱式公式;1=0則為顯式公式。5.5.1基于數(shù)值積分的構(gòu)造法將在上積分，得到只要近似地算出右邊的積分，則可通過近似y(xk+1)。而選用不同近似式I，可得到不同的計算公式。例如利用左矩形積分公式得到歐拉公式；梯形積分公式得到梯形公式。一般地，利用插值原理所建立的一系列數(shù)值積分方法也可以導(dǎo)出解微分方程的一系列計算公式。運用插值方法的關(guān)鍵在于選取合適的插值節(jié)點。假設(shè)已構(gòu)造出f(x,y(x))的插值多項式Pr(x),則5.5.2亞當(dāng)姆斯顯式公式利用r+1個節(jié)點上的被積函數(shù)值構(gòu)造r階牛頓后插多項式,有Newton插值余項Adams顯式公式j(luò)j0111/225/1233/8實際計算時，將差分展開

局部截斷誤差為：例：常用的是r=3的4階亞當(dāng)姆斯顯式公式Adams顯式公式用作為插值節(jié)點，在求積區(qū)間[xn,xn+1]上用插值函數(shù)Nr(x)近似f(x,y(x))，而xn+1不在插值節(jié)點內(nèi)，因此是一個外推的過程。雖然公式是顯式的，便于計算，但效果并不理想，比如穩(wěn)定性較差等。因此改用通過r+1個節(jié)點的插值多項式Nr(x)近似f(x,y(x))，由于xn+1是其中一個插值點，因此是內(nèi)插多項式，但導(dǎo)出的公式是隱式的。5.5.3亞當(dāng)姆斯隱式公式利用r+1個節(jié)點上的被積函數(shù)值fn+1

,fn,…,fnr+1

構(gòu)造r階牛頓前插多項式。與顯式多項式完全類似地可得到一系列隱式公式。j011-1/22-1/123-1/24局部截斷誤差為：常用的是k=3的4階亞當(dāng)姆斯隱式公式較同階顯式穩(wěn)定例：5.5.4亞當(dāng)姆斯預(yù)測-校正系統(tǒng)Step1:用Runge-Kutta法計算前i

個初值；Step2:用Adams顯式計算預(yù)測值；Step3:用同階Adams隱式計算校正值。注意：三步所用公式的精度必須相同。通常用經(jīng)典Runge-Kutta法配合4階Adams公式。4階Adams隱式公式的截斷誤差為4階Adams顯式公式的截斷誤差為預(yù)測值與校正值的事后誤差估計Adams顯隱預(yù)測-校正-改進系統(tǒng)5.5.5基于泰勒展開的構(gòu)造法)...(...110111101rnrnnnrnrnnnffffhyyyy--+---+++++++++=bbbbaaa

將通式中的右端各項yn1,…,ynk;fn+1,fn1,…,fnk(y’n+1,…,y’nk)分別在

xn點作泰勒展開，與精確解y(xn+1)在xn點的泰勒展開作比較。通過令同類項系數(shù)相等，得到足以確定待定系數(shù)0,…,r;

1,0,…,r

的等式，則可構(gòu)造出線性多步法的公式。當(dāng)10時，為隱式公式;1=0則為顯式公式。4階Adams顯式公式4階Adams隱式公式5.5.6Milne（米爾尼）格式可導(dǎo)出另一類四階方法，包含著名的米爾尼格式：4階Adams隱式公式Simpson格式用Simpson格式與米爾尼格式匹配，構(gòu)成預(yù)估校正系統(tǒng)：Milne-Simpson系統(tǒng)的缺點是穩(wěn)定性差，為改善穩(wěn)定性，考慮另一種隱式校正公式：要求公式具有4階精度。通過泰勒展開，可得到個等式，從中解出個未知數(shù)，則有個自由度。561取1=1得Simpson公式哈明