2010第3讲解线性方程组的直接法.pptVIP

数值计算方法B 教师：唐剑梅 Email：tjm112@163.com 第3章 解线性方程组的直接法 3．1 高斯消元法 3．1．1 三角形方程组 3．1．2 顺序高斯消元法 3．1．3 选主元高斯消元法 3．2 系数矩阵的三角分解法 3．2．1 高斯消元法的矩阵形式 3．2．2 LU分解法 3．2．3 解三对角线方程组的追赶法 3．2．4 平方根法与LDLT分解法 以三阶为例： 令Ux=y，由Ly=b，求出yi，再由Ux=y即可求出xi。问题的关键在于首先求出L、U中的各个元素lij、uij。其步骤如下： Doolittle(杜利特尔)分解法? 方法说明：将待求线性方程组Ax=b的系数矩阵A分解为两个矩阵的乘积A=LU，其中L为单位下三角矩阵，U为上三角矩阵。 1、计算U中的第一行元素（考虑A中的第一行元素） 由计算矩阵乘积的方法，A中的a11可由如下的计算式而得： 同理 2、计算L中的第一列元素（考虑A中的第一列元素） 由计算矩阵乘积的方法，A中的a21可由如下的计算式而得： 同理可得： 3、计算U中第二行元素（考虑A中第二行元素） 由计算矩阵乘积的方法，A中a22、a23可由如下的计算式而得： 一般地： 4、计算L中的第2列元素（考虑A中的第2列元素） 由计算矩阵乘积的方法，A中a32、a42可由如下的计算式而得： 同理 计算U中第k行元素（考