資源描述：
《有限差分法求解偏微分方程matlab.doc》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在行業(yè)資料-天天文庫。
1、.南京理工大學課程考核論文課程名稱：高等數(shù)值分析論文題目：有限差分法求解偏微分方程姓名：羅晨學號：115104000545成績：任課教師評語：簽名：年月日有限差分法求解偏微分方程word范文.一、主要內(nèi)容1.有限差分法求解偏微分方程，偏微分方程如一般形式的一維拋物線型方程：具體求解的偏微分方程如下：2.推導五種差分格式、截斷誤差并分析其穩(wěn)定性；3.編寫MATLAB程序實現(xiàn)五種差分格式對偏微分方程的求解及誤差分析；4.結論及完成本次實驗報告的感想。二、推導幾種差分格式的過程：有限差分法（finite-differencemethods）
2、是一種數(shù)值方法通過有限個微分方程近似求導從而尋求微分方程的近似解。有限差分法的基本思想是把連續(xù)的定解區(qū)域用有限個離散點構成的網(wǎng)格來代替；把連續(xù)定解區(qū)域上的連續(xù)變量的函數(shù)用在網(wǎng)格上定義的離散變量函數(shù)來近似；把原方程和定解條件中的微商用差商來近似，積分用積分和來近似，于是原微分方程和定解條件就近似地代之以代數(shù)方程組，即有限差分方程組，解此方程組就可以得到原問題在離散點上的近似解。推導差分方程的過程中需要用到的泰勒展開公式如下：(2-1)求解區(qū)域的網(wǎng)格劃分步長參數(shù)如下：(2-2)word范文.2.1古典顯格式2.1.1古典顯格式的推導由泰勒
3、展開公式將對時間展開得(2-3)當時有(2-4)得到對時間的一階偏導數(shù)(2-5)由泰勒展開公式將對位置展開得(2-6)當時，代入式(2-6)得(2-7)因為，代入上式得(2-8)得到對位置的二階偏導數(shù)word范文.(2-9)將式(2-5)、(2-9)代入一般形式的拋物線型偏微分方程得(2-10)為了方便我們可以將式(2-10)寫成(2-11)(2-12)最后得到古典顯格式的差分格式為(2-13)，古典顯格式的差分格式的截斷誤差是。2.1.2古典顯格式穩(wěn)定性分析古典顯格式(2-13)寫成矩陣形式為(2-14)上面的C矩陣的特征值是：wo
4、rd范文.(2-15)使，即結論：當時，所以古典顯格式是穩(wěn)定的。2.2古典隱格式2.2.1古典隱格式的推導將代入式(2-3)得(2-16)(2-17)得到對時間的一階偏導數(shù)(2-18)將式(2-9)、(2-18)原方程得到(2-19)為了方便把(2-19)寫成(2-20)word范文.(2-21)最后得到古典隱格式的差分格式為(2-22)，古典隱格式的差分格式的截斷誤差是。2.2.2古典隱格式穩(wěn)定性分析將古典隱格式(2-22)寫成矩陣形式如下(2-23)誤差傳播方程(2-24)所以誤差方程的系數(shù)矩陣為使，顯然恒成立。結論：對于，即任意
5、網(wǎng)格比下，古典隱格式是絕對穩(wěn)定的。word范文.2.3Richardson格式2.3.1Richardson格式的推導將，代入式(2-3)得(2-25)即(2-26)由此得到可得(2-27)將式(2-9)、(2-27)代入原方程得到下式(2-28)為了方便可以把式(2-28)寫成(2-29)即(2-30)最后得到Richardson顯格式的差分格式為(2-31)word范文.，古典顯格式的差分格式的截斷誤差是。2.3.2Richardson穩(wěn)定性分析將Richardson顯格式(2-31)寫成如下矩陣形式(2-32)誤差傳播方程矩陣形
6、式(2-33)再將上面的方程組寫成矩陣形式(2-34)系數(shù)矩陣的特征值是(2-35)解得特征值為(2-36)(恒成立)(2-37)結論：上式對任意的網(wǎng)比都恒成立，即Richardson格式是絕對不穩(wěn)定的。4.Crank-Nicholson格式3.4.1Crank-Nicholson格式的推導將代入式(2-9)得word范文.(2-40)即(2-41)得到如下方程(2-42)所以處的一階偏導數(shù)可以用下式表示：(2-43)將,代入式(2-6)可以得到式(2-9)；同理,代入式(2-6)可以得到(2-44)所以，處的二階偏導數(shù)用式(2-6)
7、、(2-44)表示：word范文.(2-45)所以，處的函數(shù)值可用下式表示：(2-46)原方程變?yōu)椋?2-47)將差分格式代入上式得：(2-48)為了方便寫成：(2-49)最后得到Crank-Nicholson格式的差分格式為(2-50)，Crank-Nicholson格式的差分格式的截斷誤差是。3.4.1Crank-Nicholson穩(wěn)定性分析Crank-Nicholson格式寫成矩陣形式如下：word范文.(2-51)誤差傳播方程是：(2-52)可以得到：(2-53)(2-54)使即(2-55)(2-56)(2-57)(2-58)
8、上式恒成立。結論：Crank-Nicholson格式對任意網(wǎng)格比也是絕對穩(wěn)定的。word范文.5.DuFortFrankle格式（Richardson格式的改進）將代入式(2-31)并化簡得到DuFortFrankle：
