最新解线方程组的直接方法课件

解线方程组的直接方法解线方程组的直接方法第五章解线性方程组的直接方法

§5.1引言解线性方程组的两类方法：直接法:经过有限次运算后可求得方程组精确解的方法(不计舍入误差)迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。（一般有限步内得不到精确解）12/12/20222第五章解线性方程组的直接方法

§5.1引言解线性方程最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件12/12/2022912/10/2022912/12/20221012/10/202210系数矩阵与常数项：12/12/202211系数矩阵与常数项：12/10/202211回代过程:12/12/202212回代过程:12/10/20221212/12/20221312/10/202213消去第一列的n-1个系数要计算n*(n-1）

个乘法。§5.2.2高斯消去法计算量12/12/202214消去第一列的n-1个系数要计算n*(n-1）个乘法。§每一步消去过程相当于左乘初等变换矩阵Lk§5.2.3矩阵的三角分解12/12/202215每一步消去过程相当于左乘初等变换矩阵Lk§5.2.3矩阵的12/12/20221612/10/202216i+1行

i+1行依次递推12/12/202217i+1行i+1行依次递推12/10/202217定理7（矩阵的LU分解）设A为n阶矩阵，如果A的顺序主子式Di≠0（i=1，2，…，n-1），则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，且这种分解是

唯一的。12/12/202218定理7（矩阵的LU分解）设A为n阶矩阵，如果A的顺序主子式12/12/20221912/10/202219§5.3高斯主元素消去法为避免此种情况的发生，可通过交换方程的次序，选取绝对值大的元素作主元。§5.3.1列主元素消去法12/12/202220§5.3高斯主元素消去法为避免此种情况的发生，可通过交换方选取或称此方法为全主元素高斯消去法称此方法为列主元素高斯消去法12/12/202221选取或称此方法为全主元素高斯消去法称此方法为列主元素高斯消去3m31=-200012/12/2022223m31=-200012/10/202222m31=-0.000512/12/202223m31=-0.000512/10/202223定理8（列主元素的三角分解定理）如果A为非奇异矩阵，则存在排列矩阵P使

PA=LU其中L为单位下三角阵，U为上三角阵。12/12/202224定理8（列主元素的三角分解定理）如果A为非奇异12/10/§5.3.2高斯—若当消去法12/12/202225§5.3.2高斯—若当消去法12/10/20222512/12/20222612/10/202226§5.4矩阵三角分解法

§5.4.1直接三角分解法将高斯消去法改写为紧凑形式，可以直接从矩阵A的元素得到计算L,U元素的递推公式，而不需要任何中间步骤，这就是直接三角分解法。由于A=LU，求解Ax=b的问题就等价于求解两个三角形方程组①Ly=b，求y；②Ux=y，求x.12/12/202227§5.4矩阵三角分解法

§5.4.1直接三角分解法1、不选主元的三角分解法A=LU其中L为单位下三角阵，U为上三角阵（4.1）12/12/2022281、不选主元的三角分解法A=LU其中L为单位下三角阵，U为上一、直接计算A的LU分解(例)12/12/202229一、直接计算A的LU分解(例)12/10/2022212/12/20223012/10/202230二、一般计算公式12/12/202231二、一般计算公式12/10/202231三、LU分解求解线性方程组12/12/202232三、LU分解求解线性方程组12/10/202232矩阵A的直接分解法称为杜利特尔（Doolittle）分解12/12/202233矩阵A的直接分解法称为杜利特尔（Doolittle）分解12例1：将方程组的系数矩阵A作LU分解，并求方程组的解12/12/202234例1：将方程组的系数矩阵A作LU分解，并求方程组的解12/1解LU分解的紧凑格式为12/12/202235解LU分解的紧凑格式为12/10/202235推出：由Ly=b得12/12/202236推出：由Ly=b得12/10/202236由Ux=y，即用回代法解得即为线性方程组的解12/12/202237由Ux=y，即用回代法解得即为线性方程组的解12/10/202、选主元的三角分解法采用与列主元消去法类似的方法，通过交换A的行实现矩阵PA的LU分解。12/12/2022382、选主元的三角分解法采用与列主元消去法类似的方法，12/1本章作业1、设Ax=b，用LU分解法求此方程组2、P230712/12/202239本章作业1、设Ax=b，用LU分解法求此方程组2、

结束语谢谢大家聆听！！！40

结束语谢谢大家聆听！！！40解线方程组的直接方法解线方程组的直接方法第五章解线性方程组的直接方法

§5.1引言解线性方程组的两类方法：直接法:经过有限次运算后可求得方程组精确解的方法(不计舍入误差)迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。（一般有限步内得不到精确解）12/12/202242第五章解线性方程组的直接方法

§5.1引言解线性方程最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件最新解线方程组的直接方法课件12/12/20224912/10/2022912/12/20225012/10/202210系数矩阵与常数项：12/12/202251系数矩阵与常数项：12/10/202211回代过程:12/12/202252回代过程:12/10/20221212/12/20225312/10/202213消去第一列的n-1个系数要计算n*(n-1）

个乘法。§5.2.2高斯消去法计算量12/12/202254消去第一列的n-1个系数要计算n*(n-1）个乘法。§每一步消去过程相当于左乘初等变换矩阵Lk§5.2.3矩阵的三角分解12/12/202255每一步消去过程相当于左乘初等变换矩阵Lk§5.2.3矩阵的12/12/20225612/10/202216i+1行

i+1行依次递推12/12/202257i+1行i+1行依次递推12/10/202217定理7（矩阵的LU分解）设A为n阶矩阵，如果A的顺序主子式Di≠0（i=1，2，…，n-1），则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，且这种分解是

唯一的。12/12/202258定理7（矩阵的LU分解）设A为n阶矩阵，如果A的顺序主子式12/12/20225912/10/202219§5.3高斯主元素消去法为避免此种情况的发生，可通过交换方程的次序，选取绝对值大的元素作主元。§5.3.1列主元素消去法12/12/202260§5.3高斯主元素消去法为避免此种情况的发生，可通过交换方选取或称此方法为全主元素高斯消去法称此方法为列主元素高斯消去法12/12/202261选取或称此方法为全主元素高斯消去法称此方法为列主元素高斯消去3m31=-200012/12/2022623m31=-200012/10/202222m31=-0.000512/12/202263m31=-0.000512/10/202223定理8（列主元素的三角分解定理）如果A为非奇异矩阵，则存在排列矩阵P使

PA=LU其中L为单位下三角阵，U为上三角阵。12/12/202264定理8（列主元素的三角分解定理）如果A为非奇异12/10/§5.3.2高斯—若当消去法12/12/202265§5.3.2高斯—若当消去法12/10/20222512/12/20226612/10/202226§5.4矩阵三角分解法

§5.4.1直接三角分解法将高斯消去法改写为紧凑形式，可以直接从矩阵A的元素得到计算L,U元素的递推公式，而不需要任何中间步骤，这就是直接三角分解法。由于A=LU，求解Ax=b的问题就等价于求解两个三角形方程组①Ly=b，求y；②Ux=y，求x.12/12/202267§5.4矩阵三角分解法

§5.4.1直接三角分解法1、不选主元的三角分解法A=LU其中L为单位下三角阵，U为上三角阵（4.1）12/12/2022681、不选主元的三角分解法A=LU其中L为单位下三角阵，U为上一、直接计算A的LU分解(例)12/12/202269一、直接计算A的LU分解(例)12/10/2022212/12/20227012/10/202230二、一般计算公式12/12/202271二、一般计算公式12/10/202231三、LU分解求解线性方程组12/12/202272三、LU分解求解线性方程组12/10/202232矩阵A的直接分解法称为杜利特尔（Doolittle）分解12/12/202273矩阵A的直接分解法称为杜利特尔（Doolittle）分解12例1：将方程组的系数矩阵A作LU分解，并求方程组的解12/12/202274例1：将方程组的系数矩阵A作LU分解，并求方程组的解12/1解LU分解的紧凑格式为12/12/202275解LU分解的紧凑格式为12/10/202235推出：由Ly=b得12/12/202276推出：由Ly=b得12/10/202236由Ux=y，即用回代法解得即为线性方程组的解12/12/202277由Ux