第10講收斂性與穩(wěn)定性方程組與高階方程

1、第五章 常微分方程的差分方法5.3線性多步法一、教學(xué)目標(biāo)及基本要求通過對本節(jié)課的學(xué)習(xí)，使學(xué)生掌握常微分方程、常微分方程方程組的線性多 步法。二、教學(xué)內(nèi)容及學(xué)時分配本節(jié)課主要介紹常微分方程的數(shù)值解法。具體內(nèi)容如下：講授內(nèi)容：歐拉公式、改進(jìn)的歐拉公式。三、教學(xué)重點難點1 教學(xué)重點：開型求解公式，閉型求解公式。2.教學(xué)難點：收斂性與穩(wěn)定性。四、教學(xué)中應(yīng)注意的問題多媒體課堂教學(xué)為主。適當(dāng)提問，加深學(xué)生對概念的理解五、正文線性多步法及其收斂性與穩(wěn)定性、方程組與高階方程1引言收斂性問題微分方程數(shù)值解法的基本思想是：通過某種離散化手段，將微分方程轉(zhuǎn)化為 差分方程（代數(shù)方程）來求解。這種轉(zhuǎn)化是否合理，還要看差

2、分問題的解yn，當(dāng)h 0時是否會收斂到微分方程的準(zhǔn)確解y（Xn），需要注意的是，如果只考慮 ho，那么節(jié)點X Xo nh對固定的n將趨向于Xo，這時討論收斂性是沒有意 義的，因此，當(dāng)h0時，同時n時才合理。定義：若一種數(shù)值方法對于任意固定的XXo nh，當(dāng)h 0 （同時n ：） 時，有yn * y（Xn）,則稱該方法是收斂的考察歐拉公式 yn 1 = yn hf(Xn,yn)()設(shè)yn 1為在yn二y(Xn)條件下按歐拉公式計算的結(jié)果,yn 1 =y(Xn)hf(Xn,y(Xn)( 2)y(Xn.l) -yn1即為局部截斷誤差。h2 .Tn “y(Xn 1) yn 2 y()，存在常數(shù) C使

3、y(Xn 1) -石:Ch2(3)考慮整體截斷誤差編舟=$匕卄)_%y(Xn 弭一yn 0時宀 0，歐拉公式是收斂的。進(jìn)一步考察一般的單步法：所謂單步法，就是在計算yn1時只用到它前一步的信息yn。顯式單步法的共同特征是，它們都是將yn加上某種形式的增量得出% 1 ,其計算公式的形式為：yn 1二yn h (冷，yn，h),-(冷，稱為增量函數(shù), 不同的單步法，對應(yīng)不同的增量函數(shù)。定理：單步法滿足條件(x,y，h)(x,y，h)|n |yy| (李普希茲條件)， 且設(shè)初值yo是準(zhǔn)確的，即yo =y(Xo)，則該單步法是收斂的。2穩(wěn)定性問題對于一個數(shù)值方法，即使是收斂的，由于初始值一般都帶有誤差

4、，同時，在 計算過程中還常常產(chǎn)生舍入誤差，這些誤差又必然會傳播下去，對后續(xù)的計算結(jié) 果都將產(chǎn)生影響，數(shù)值穩(wěn)定性問題是討論這種誤差的積累和傳播能否得到控制的 問題。定義若用某一數(shù)值方法計算yn時，所得到的實際計算結(jié)果為yn，且由擾 動；n =| yn -輪|引起以后各節(jié)點ym(m n)的擾動為，如果總有丨m |J 5 |,則稱 該方法是穩(wěn)定的。一種數(shù)值方法是否穩(wěn)定，不僅與該數(shù)值方法本身有關(guān)，而且還與微分方程的 右端函數(shù)f(x，y)，以及步長h有關(guān)，因此穩(wěn)定性問題比較復(fù)雜。為了簡化討論 只考慮模型方程z o y() = y0歐拉公式穩(wěn)定性：yn (1 h )ynyn處有擾動n，它的傳播使節(jié)點 齊1

5、產(chǎn)生擾動，假設(shè)歐拉公式計算中不再引入新誤差，則n(1 h，y n如果原差分方程yn 1 = (1 h )yn的解不增長，即有|yn 1口興|，就能保證歐 拉方法的穩(wěn)定性。yn (1 h )yn的解不增長，h需要充分小，使|11o故歐拉方法是條件穩(wěn)定的隱式歐拉公式穩(wěn)定性:yn 1yn 1，從而|yn1 Iyn I ,隱式歐拉公式是恒穩(wěn)定的。3方程組與高階方程(1)一階方程組直接推廣各種算法到方程組y = f (x, y,z), y(x) =y如 z =g(x, y,z),z(x) =z令x x nh，yn,Zn表示節(jié)點xn上的近似解改進(jìn)的歐拉公式為：ynyn - hf (Xn, 丫.冃)預(yù)報 Z

6、ni % hg(xn,yn,zn)校正rhyn 1 = yn - - f (Xn, yn,Zn) f (Xn 1, yn 1, Zn 1)hZn 1 二 Zng(xn,yn,zn) g(Xn 1,n 1,厶 JL2四階龍格一庫塔方法為：hYnYn -K1 2K2 2K3 K46Zn1 二厶L1 2L2 2L3 L40 二 f (Xn,yn,Zn),Li 二 9(心幾，召)= g(xn羊，* +沙,召+加)n 2 2 2心=胡,幾+尹1,召+2）丄2n 2 2 2hhhh&厲嚴(yán)n+尹2尹）丄3弋甘尹2產(chǎn)） 二 f (Xn 1, yn hK3,厶 1), L4 二 g(Xn 1, y. h,厶 hL?)（2）化高階方程為一階方程組y =f(x,y,y)對.y(X0)=y0,y(X0)= y0，引入新變量Z = y即可化為一階方程組:二 Z, y(Xo) = yoI二 f (x, y, z),z(Xo) = y