資源描述：
《橢圓型方程的有限差分方法（I）ppt課件.ppt》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、1第五章橢圓型方程的差分方法(一)Poisson方程（二）差分格式的性質(zhì)（三）邊界條件的處理（四）變系數(shù)方程（五）雙調(diào)和方程（六）特征值問題2(一)Poisson方程區(qū)域以及邊界離散4內(nèi)點(diǎn)邊界點(diǎn)返回5返回1.五點(diǎn)差分格式利用Taylor級(jí)數(shù)展開有則8由于上差分方程中只出現(xiàn)u在(i,j)及其四個(gè)臨點(diǎn)上的值,故稱為五點(diǎn)差分格式Possion方程第一邊值問題的差分逼近為改寫為代數(shù)方程,僅考慮解：11先僅下標(biāo)j小的未知量放前面,相同的j的再按i從小到大順利排列11按自然順序排列網(wǎng)點(diǎn)(i,j)按自然順序排列網(wǎng)點(diǎn)(i,j)
2、定義向量于是差分方程為:12分析系數(shù)矩陣H對(duì)于第一個(gè)結(jié)點(diǎn)(1,1),未知節(jié)點(diǎn)分析系數(shù)矩陣H對(duì)于第二個(gè)結(jié)點(diǎn)(2,1),分析系數(shù)矩陣H對(duì)于前I個(gè)結(jié)點(diǎn),系數(shù)矩陣H為五點(diǎn)差分格式等價(jià)于B,I為I維方陣,H為I*J維方陣對(duì)于第一個(gè)結(jié)點(diǎn),即例：17同理,對(duì)于另外三個(gè)結(jié)點(diǎn),分別有聯(lián)立,有線性代數(shù)方程組18計(jì)算機(jī)實(shí)現(xiàn)五點(diǎn)差分格式的系數(shù)矩陣算法程序返回19五點(diǎn)差分格式的系數(shù)矩陣對(duì)possion方程第一邊值問題建立五點(diǎn)差分格式將雙下標(biāo)數(shù)列ui,j按自然順序排成單下標(biāo)uk,那么20這樣有代數(shù)方程組Hu=g,其中21這樣有代數(shù)方程組H
3、u=g,其中H是稀疏矩陣,每一行最多5個(gè)數(shù)不為0H是對(duì)角占優(yōu)矩陣H是對(duì)稱正定課堂練習(xí)23微分方程差分格式解u=u(x,y)解u=un實(shí)用性分析:相容性收斂性穩(wěn)定性24相容性所以五點(diǎn)差分格式是相容的.2.九點(diǎn)差分格式27四階精度3極坐標(biāo)形式的差分格式若求解域是圓環(huán)、環(huán)形域或扇形域，則采用極坐標(biāo)是方便的，此時(shí)Poisson方程形如任意點(diǎn)用中心差商公式分別取等步長亦可用有限體積法獲得P1411.用五點(diǎn)差分格式求解Poisson方程的邊值問題其中作業(yè)
