關(guān)于拋物型方程的差分方法第1頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

其中,為平面上某一區(qū)域。眾所周知，一維線性拋物型方程的一般形式為第2頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四通?？紤]的定解問(wèn)題有：(1)初值問(wèn)題

在區(qū)域上求函數(shù)，使?jié)M足

為給定的初始函數(shù)。第3頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四(2)初邊值問(wèn)題(或稱混合問(wèn)題)

在區(qū)域上求函數(shù)，使?jié)M足第4頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

為了構(gòu)造微分方程的有限差分逼近，首先將求解區(qū)域用二組平行于軸和軸的直線構(gòu)成的網(wǎng)格覆蓋，網(wǎng)格邊長(zhǎng)在方向?yàn)?，在方向?yàn)?。分別稱為空間方向和時(shí)間方向的步長(zhǎng)，網(wǎng)格線的交點(diǎn)稱為網(wǎng)格的結(jié)點(diǎn)。差分格式的建立第5頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四由Taylor展開(kāi)，有

則在處對(duì)的一階偏導(dǎo)數(shù)有三個(gè)可能的近似:向后差商向前差商中心差商第6頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

顯然，用差商近似導(dǎo)數(shù)存在誤差，令則截?cái)嗾`差第7頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四現(xiàn)記前差算子:,后差算子:,中心差算子:,為方向偏導(dǎo)數(shù)算子為方向位移算子，為方向平均算子，其中:,第8頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

建立差分算子和導(dǎo)數(shù)算子之間的關(guān)系由得或者同理有第9頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四因?yàn)楣释硪驗(yàn)?/p>

則第10頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

利用這些關(guān)系式就可給出偏導(dǎo)數(shù)的差分表達(dá)式第11頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四又由可得二階偏導(dǎo)數(shù)的差分表達(dá)式第12頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

從以上這些偏導(dǎo)數(shù)的差分表達(dá)式，我們可以得到偏導(dǎo)數(shù)的各種精度的近似表達(dá)式。且

又由二階導(dǎo)數(shù)的前差表達(dá)式，得因此

在的前差表達(dá)式中取第一項(xiàng)，則有即截?cái)嗾`差階為。第13頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

現(xiàn)在研究構(gòu)造微分方程的差分方程的方法，為此記微分方程為

L

是關(guān)于的線性算子，。包括二個(gè)相鄰時(shí)間層的網(wǎng)格結(jié)點(diǎn)的差分方程可以從Talor

展開(kāi)式推出第14頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

現(xiàn)在，對(duì)拋物型方程的幾種特殊情況，從方程出發(fā)，構(gòu)造微分方程的有限差分近似。

首先考慮一維熱傳導(dǎo)方程的差分近似。顯式格式第15頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四由，方程為代入則其中為步長(zhǎng)比。第16頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四在上式中，如果僅僅保留二階中心差分，且設(shè)為相應(yīng)差分方程解在結(jié)點(diǎn)(mh,nk)上的值，則代入的表達(dá)式，則得差分方程將格式應(yīng)用于解初值問(wèn)題第17頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四此差分格式也可簡(jiǎn)單地由導(dǎo)數(shù)的差商近似表達(dá)式得到代入微分方程，并令差分方程解為即可。雖然在邊界結(jié)點(diǎn)上，差分方程和微分方程具有相同的初值或者初邊值條件，但是，一般而言，結(jié)點(diǎn)上微分方程的精確解和古典顯式差分格式的精確解不相等。記第18頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

假定具有下面推導(dǎo)中所需要的有界偏導(dǎo)數(shù)，則由展開(kāi)，有截?cái)嗾`差42第19頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四則那么得從而有第20頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四或第21頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四從而，上式右邊量描寫(xiě)了古典顯式差分格式在點(diǎn)對(duì)微分方程的近似程度，將其定義為差分格式在點(diǎn)的截?cái)嗾`差，記為，即

假定在所考慮的區(qū)域保持有界，則古典顯式差分格式的截?cái)嗾`差階為。第22頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四或者相應(yīng)的截?cái)嗾`差階為。通常，格式可用下圖表示。

為了提高截?cái)嗾`差的階，我們也可用在式中保留四階中心差分項(xiàng)的辦法達(dá)到，這時(shí)有差分格式第23頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四m,n+1m-2,nm-1,nm,nm+1,nm+2,nm,n+1m-1,nm,nm+1,n第24頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四隱式格式隱式差分格式特點(diǎn)：1.具有二個(gè)或二個(gè)以上結(jié)點(diǎn)處的值未知；2.計(jì)算工作量較大；3.穩(wěn)定性較好。第25頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四得由推導(dǎo)其最簡(jiǎn)單的隱式差分逼近─古典隱式格式。

現(xiàn)在對(duì)熱傳導(dǎo)方程第26頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四格式用下圖表示，其截?cái)嗾`差階為，與古典顯式差分格式相同?；蛘弑Ａ舳A導(dǎo)數(shù)項(xiàng)，且以替代，則得差分格式

我們也可通過(guò)直接用差分算子代替的方法，即代入微分方程，得到此格式。第27頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四m,n+1m-1,n+1m+1,n+1m,n第28頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四圖方法第29頁(yè)，共33頁(yè)，2022年，5月20日，22點(diǎn)19分，星期四

-4

-3

-2

-1

0

1

2

3

4

0

1

2

3

4