解線方程組的直接方法.ppt

2019/7/18,1,計算方法,解線性方程組的直接方法,2019/7/18,2,第五章 解線性方程組的直接方法 5.1 引言,解線性方程組的兩類方法： 直接法: 經(jīng)過有限次運算后可求得方程組精確解的 方法(不計舍入誤差) 迭代法：從解的某個近似值出發(fā)，通過構(gòu)造一個無窮序列去逼近精確解的方法。（一般有限步內(nèi)得不到精確解）,2019/7/18,3,n階線性方程組,2019/7/18,4,2019/7/18,5,5.2 高斯消去法,轉(zhuǎn)化為同解的方程組,2019/7/18,6,5.2.1 高斯消去法計算過程,2019/7/18,7,2019/7/18,8,2019/7/18,9,2019/7/18,10,2019/7/18,11,系數(shù)矩陣與常數(shù)項：,2019/7/18,12,回代過程:,2019/7/18,13,2019/7/18,14,消去第一列的 n-1 個系數(shù)要計算n*(n-1） 個乘法。,5.2.2 高斯消去法計算量,2019/7/18,15,每一步消去過程相當于左乘初等變換矩陣Lk,5.2.3 矩陣的三角分解,2019/7/18,16,2019/7/18,17,依次遞推,2019/7/18,18,定理7（矩陣的LU分解） 設(shè)A為n階矩陣，如果A的順序主子式 Di0（i=1，2，n-1），則A可分解為一個單位下 三角矩陣L和一個上三角矩陣U的乘積，且這種分解是 唯一的。,2019/7/18,19,2019/7/18,20,5.3 高斯主元素消去法,為避免此種情況的發(fā)生，可通過交換方程的次序，選取絕對值大的元素作主元。,5.3.1 列主元素消去法,2019/7/18,21,選取,或,2019/7/18,22,2019/7/18,23,m31= -0.0005,2019/7/18,24,定理8（列主元素的三角分解定理） 如果A為非奇異 矩陣，則存在排列矩陣P使 PA=LU 其中L為單位下三角陣，U為上三角陣。,2019/7/18,25,5.3.2 高斯若當消去法,2019/7/18,26,2019/7/18,27,5.4 矩陣三角分解法 5.4.1 直接三角分解法,將高斯消去法改寫為緊湊形式，可以直接從矩陣A的元素 得到計算L,U元素的遞推公式，而不需要任何中間步驟， 這就是直接三角分解法。,由于A=LU，求解Ax=b的問題就等價于求解兩個三角形方程組, Ly=b，求y；, Ux=y，求x.,2019/7/18,28,1、不選主元的三角分解法,A=LU,其中L為單位下三角陣，U為上三角陣,（4.1）,2019/7/18,29,一、直接計算 A 的 LU 分解(例),2019/7/18,30,2019/7/18,31,二、一般計算公式,2019/7/18,32,三、LU 分解求解線性方程組,2019/7/18,33,矩陣A的直接分解法稱為杜利特爾（Doolittle）分解,2019/7/18,34,例1：將方程組,的系數(shù)矩陣A作LU分解，并求方程組的解,2019/7/18,35,解,LU分解的緊湊格式為,2019/7/18,36,推出：,由Ly=b得,2019/7/18,37,由Ux=y，即,用回代法解得,即為線性方程組的解,2019/7/18,38,2、選主元的三角分解法,采用與列主元消去法類似的