拋物型方程的差分方法

關(guān)于拋物型方程的差分方法第一頁(yè)，共三十三頁(yè)，2022年，8月28日

其中,為平面上某一區(qū)域。眾所周知，一維線性拋物型方程的一般形式為第二頁(yè)，共三十三頁(yè)，2022年，8月28日通?？紤]的定解問題有：(1)初值問題

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

為給定的初始函數(shù)。第三頁(yè)，共三十三頁(yè)，2022年，8月28日(2)初邊值問題(或稱混合問題)

在區(qū)域上求函數(shù)，使?jié)M足第四頁(yè)，共三十三頁(yè)，2022年，8月28日

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

則在處對(duì)的一階偏導(dǎo)數(shù)有三個(gè)可能的近似:向后差商向前差商中心差商第六頁(yè)，共三十三頁(yè)，2022年，8月28日

顯然，用差商近似導(dǎo)數(shù)存在誤差，令則截?cái)嗾`差第七頁(yè)，共三十三頁(yè)，2022年，8月28日現(xiàn)記前差算子:,后差算子:,中心差算子:,為方向偏導(dǎo)數(shù)算子為方向位移算子，為方向平均算子，其中:,第八頁(yè)，共三十三頁(yè)，2022年，8月28日

建立差分算子和導(dǎo)數(shù)算子之間的關(guān)系由得或者同理有第九頁(yè)，共三十三頁(yè)，2022年，8月28日因?yàn)楣释硪驗(yàn)?/p>

則第十頁(yè)，共三十三頁(yè)，2022年，8月28日

利用這些關(guān)系式就可給出偏導(dǎo)數(shù)的差分表達(dá)式第十一頁(yè)，共三十三頁(yè)，2022年，8月28日又由可得二階偏導(dǎo)數(shù)的差分表達(dá)式第十二頁(yè)，共三十三頁(yè)，2022年，8月28日

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

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

在的前差表達(dá)式中取第一項(xiàng)，則有即截?cái)嗾`差階為。第十三頁(yè)，共三十三頁(yè)，2022年，8月28日

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

L

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

展開式推出第十四頁(yè)，共三十三頁(yè)，2022年，8月28日

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

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

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

假定在所考慮的區(qū)域保持有界，則古典顯式差分格式的截?cái)嗾`差階為。第二十二頁(yè)，共三十三頁(yè)，2022年，8月28日或者相應(yīng)的截?cái)嗾`差階為。通常，格式可用下圖表示。

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

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

我們也可通過直接用差分算子代替的方法，即代入微分方程，得到此格式。第二十七頁(yè)，共三十三頁(yè)，2022年，8月28日m,n+1m-1,n+1m+1,n+1m,n第二十八頁(yè)，共三十三頁(yè)，2022年，8月28日?qǐng)D方法第二十九頁(yè)，共三十三頁(yè)，2022年，8月28日

-4

-3

-2

-1

0

1

2

3

4

0

1

2

3