與PDE數(shù)值解相關的線性代數(shù)方程組求解.pdf

1、在偏微分方程數(shù)值解中,常采用有限差分法、有限元法以及有限體積法等方法,其實質(zhì)都是將連續(xù)的問題離散化,最后問題常歸結為方程組的求解,比較常見的是大型線性方程組,非線性方程組也往往要轉換成線性代數(shù)方程組的求解.因此對于數(shù)值求解偏微分方程而言‘,線代數(shù)方程組求解是非常重要的一環(huán).方程組求解方法主要是數(shù)值代數(shù)研究的領域,對于偏微分方程數(shù)值解研究方向的科技人員而言,往往并不親自去研究方程組的求解方法以及理論分析,而注重的是算法的使用、算法的結論以

2、及算法的計算效果(計算精度、快慢等).本文寫作目的即是對大型線代數(shù)方程組求解的各種方法進行討論,理論與實踐相結合,對一些尚不完善的理論問題進行深入分析,并選擇Poisson方程差分格式對各種方法進行綜合實驗比對,提供求解方程組時算法選擇的依據(jù).在直接法部分,首先應該指出在條件數(shù)大小與病態(tài)/良態(tài)問題的關系、判斷一個具體問題是否病態(tài)的理論分析方法等還有許多值得細致研究的問題,本文只就通常的數(shù)值分析教材中不明確或易混淆的一些地方進行了說明,并

3、對一些理論問題進行了深入探討,提出一些疑問.其次針對高斯列主元消去法進行了具體的理論與實踐相結合的分析:對稠密矩陣進行了數(shù)值實驗,對該算法適應方程組類型、可求解方程組的規(guī)模、計算速度、計算結果精度等問題給出了答案;分析了列主元消去法的舍入誤差理論,給出了因舍入誤差影響導致計算失效的一個例子,指出了舍入誤差理論不完善的地方;針對典型的偏微分方程數(shù)值解中的Poisson方程問題進行了進一步實驗,了解將連續(xù)問題經(jīng)由差分格式離散化后產(chǎn)生的大型稀

4、疏型方程組用直接法來編程求解的過程以及最后的計算效果.最后結合一維兩點邊值問題數(shù)值實驗討論了三對角方程組求解的追趕法.在迭代法一節(jié),分析了各種基本迭代算法的收斂條件、收斂速度,為比較算法之間計算效果,仍采用Poisson方程問題對各種算法進行了數(shù)值實驗綜合比對.進一步地,針對Jacobi方法收斂速度慢的問題,本文給出了Jacobi方法結合Chebyshev半迭代算法,數(shù)值實驗驗證了它的計算速度和精度的提高;針對逐次超松弛(SOR)迭代法