求解线性方组的直接解法.docVIP

求解线性方程组的直接解法 5.2?????LU分解 ①?? Gauss消去法实现了LU分解 顺序消元结束时的上三角矩阵U和所用的乘数,严格下三角矩阵。将下三角矩阵的对角元改成1,记为L,则有A=LU, 这事实是一般的,我们不难从消去的第k个元素时的矩阵k行及k列元素的历史得到这一点.因为从消元的历史有 ukj=akj-mk1u1j- mk2u2j -…- mk,k-1uk-1,j, j=k,k+1,…,n mik=(aik-mi1u1k- mi2u2k -…-mi,k-1uk-1,k)/ukk i=k+1,k+2,…,n 于是 akj=mk1u1j+mk2u2j+…+mk,k-1uk-1,j+ukj, j=k,k+1,…,n aik=mi1u1k+mi2u2k+…+mi,k-1uk-1,k+mikukk i=k+1,k+2,…,n 从前面两个式子我们可以直接计算L和U(见下段).将矩阵分解为单位下三角矩阵和上三角矩阵之积称为矩阵的LU分解.顺序消元实现了LU分解,同时还求出了g, Lg=b的解. ② ? 直接LU分解 上段我们得到(lij=mij) ukj=akj-lk1u1j-lk2u2j -…- lk,k-1uk-1,j, j=k,k+1,…,n lik=(aik-li1u1k-li2u2k -…-li,k-1uk-1,k)/ukk i=k+