偏微分方程數(shù)值解法答案

1、1. 課本有證明2. 課本有說明3. 課本有說明4. Rit2法，設是u的n維子空間，是的一組基底，中的任一元素可表為，則是的二次函數(shù)，令，從而得到滿足，通過解線性方程組，求的，代入，從而得到近似解的過程稱為Rit2法 簡而言之，Rit2法：為得到偏微分方程的有窮維解，構(gòu)造了一個近似解，利用確定，求得近似解的過程Galerkin法：為求得形式的近似解，在系數(shù)使關于，滿足，對任意或（取）的情況下確定，從而得到近似解的過程稱Galerkin法為Rit2-Galerkin法方程：5. 有限元法：將偏微分方程轉(zhuǎn)化為變分形式，選定單元的形狀，對求解域作剖分，進而構(gòu)造基函數(shù)或單元形狀函數(shù)，形成有限元空間，

2、將偏微分方程轉(zhuǎn)化成了有限元方程，利用有效的有限元方程的解法，給出偏微分方程近似解的過程稱為有限元法。6. 解：對求解區(qū)間進行網(wǎng)格剖分，節(jié)點得到相鄰節(jié)點之間的小區(qū)間，由節(jié)點上的一組值，按線性插值公式 確定試探空間，令把變到軸上的參考但愿0,1令則：，將帶入該函數(shù)得到：帶入可得 令其中從而得到的線性方程組！7.矩形剖分假定區(qū)域C1可以分割成為有限個互不重疊的矩形的和，且每個小矩形的邊和坐標軸平行，任意兩個矩形或者不相交或者有公共的邊和公共的頂點，成如此的分割為矩形剖分基函數(shù)的取法其中是以為頂點的矩形單元 為的底和商的長度。8. 何為三角剖分，基函數(shù)怎樣??？三角剖分：設G是多邊形域（否則可用多邊形域

3、逼近它），將G分割成有限個三角形之和，使不同三角形無重疊的內(nèi)部，且任一三角形的頂點不屬于其他三角形的內(nèi)部，這樣就把G分割成三角形網(wǎng)，稱為G的三角剖分?；瘮?shù)的取法：通過構(gòu)造Lagrange型插值公式可以得到基函數(shù)的取法。不妨以是一次多項式為例，得到，其中L1是相應于節(jié)點1的基函數(shù)在上的限制（具體的過程，可參考課本：P57 P58）9.題，參考課后習題P92的第一題，具體過程可參考積分插值的推導過程10，11題不會。 在此將14題推導過程介紹如下：12. 對Possion方程，建立五點差分格式，并估計截斷誤差。取定沿x軸和y軸方向的步長h1和h2，沿x,y方向分別用二階中心差商代替，則 （五點差

4、分格式）式中表示節(jié)點(i,j)上的網(wǎng)函數(shù)。令 利用Taylor展式有截斷誤差為13. 對possion方程建立，極坐標形式的差分格式poission 方程的極坐標形式為 - 其中 利用中心差商公式 - - 將 兩式代入式得 即poission 方程極坐標形式的差分方程。14. 解：將按照Taylor在處展開整理得到其截斷誤差為在Richardson格式（）中以代入，便得Du Fort Frankel格式： - - - 得 （省去了的商階無窮小）從而得到了微分方程左邊的誤差同理可得微分方程右邊的誤差：從而得到 15.用Fourier方法討論向前差分格式的穩(wěn)定性。解：向前差分格式。以代入得消去則知

5、增長因子由于在0，中分布稠密，(隨)為使?jié)M足von Neu-Mann條件，必須且只須網(wǎng)比所以向前差分格式的穩(wěn)定性條件是16. 用Fourier方法討論向后差分格式的穩(wěn)定性。解：對向后差分方程利用Fourier方法分析器穩(wěn)定性，整理得：。令，將代入得到：消去。則增長因子為。所以向后差分方程是恒穩(wěn)定的。17. 用Fourier方法討論六點對稱格式的穩(wěn)定性。解：六點對稱方程的格式為。令代入得= 。消去得增長因子為。所以六點對稱格式是無條件恒定的。18.證明：利用Fourier方程將兩端同時做變換。得)消去exp(ixjh)得增長因子為.即差分格式的充要條件是19.討論三維熱傳導方程向前差分格式的穩(wěn)定

6、性解：三位熱傳導方程為（向前差分格式）. 取通項代到上式消去公因子得。從而增長因子為為使|=1+O()必須且只須。當時三維熱傳導方程的向前差分格式穩(wěn)定。20. 討論三維熱傳導方程向前后分格式的穩(wěn)定性解：三維熱傳導方程的向后差分格式為：取通項=exp(i(+),=,=,=,帶入上式，消去共因子得：。恒成立所以 三維傳導方程向后差分格式是無條件穩(wěn)定的。21三維傳導方程的PR格式為： = （1） = (2) = (3)(1)(2)(3)合稱PR格式。22.將= = =將=exp帶入上式得= 對任何r01. 所以=絕對收斂。23解： (1) (2) (3) (4)(1)+(2) 得=+（3）+（4）得=+所以 其截斷誤差為 。24. 證