方程组同直接三角分解法.pptVIP

第四章方程组的直接解法 4.2 直接三角分解法 4.2.3 平方根法 4.2.1 一般矩阵的直接三角分解法 4.2.2 三对角方程组的追赶法 4.2 直接三角分解法 4.2.1 一般矩阵的直接三角分解法 　　本节讨论矩阵A的三角分解法的直接计算以及直接利用A的三角分解式来求解方程组。 1.不选主元的三角分解法 　　设A=LU，记 　　其中L为单位下三角阵，U为上三角阵。我们可直接给出L和U的元素的计算公式。 　　由A的第1行和第1列可计算出U的第1行和L的第1列，即 （4.2.1） （4.2.2） 如果U的第1至k-1列和L的第1至k-1列已经算出，则由 可得U的第k行元素 同理，由 ukj =akj - ,j =k,k+1, ···,n。 （4.2.3） akj = ,i=k+1,k+2,···,n, 可得L的第k列元素 　　交替使用（4.2.3）和 （4.2.4），就能逐次计算出U（按行）和L（按列）的全部元素，而且可以把它们存放在矩阵A对应的位置上（L的对角线元素不必存放