计算方法课件 第三章 线性代数方程组数值解法.ppt

（2）用全主元高斯-若当消去法 故得解为 （3）用列主元高斯消去法 回代解得 3.3 解线性方程组的矩阵分解法 一、 非对称矩阵的三角分解法 对于给定的线性方程组 矩阵分解法的基本思想是： （1） 分解 可逆下三角矩阵 可逆上三角矩阵 显见S是一个可逆的下三角阵 ——解两个三角形方程组。 Crout分解（以四阶为例） 2.利用三角分解法解方程组 例1. 试用克洛特分解法解线性方程组 例2 试用克洛特分解法解线性方程组 解 3.3.2 解三对角型线性方程组的追赶法 1.用LU分解矩阵A 3.3.3 对称正定矩阵的三角分解 定义 3.1 若n 阶方矩阵 A 具有性质 且对任何n 维向量 成立 ，则称 A 为对称正定矩阵。 定理3.4 若A 为对称正定矩阵，则 (1) A的k阶顺序主子式 (2)有且仅有一个单位下三角矩阵L和对角矩阵D 使得 （3-16） 这称为矩阵的乔里斯基（Cholesky）分解。 (3)有且仅有一个下三角矩阵 ，使 （3-17） 这称为分解矩阵的平方根法。 * 第三章 线性代数方程组的数 值解法 3.1 引言 3.2 解线性方程组的消去法 3.3 解线性方程组的矩阵分解法 3.4 解线性方程组的迭代法 3.1 引言 给定一个线性方程组 求解向量 x。 第一类是直接法。即按求精确解的方法运算求解。 第二类是迭代法。其思想是首先把线性方程组(3-1)等价变换为如下形式的方程组： 数值解法主要有两大类： 然后构造迭代格式 这称为一阶定常迭代格式，M 称为迭代矩阵。 3.2 解线性方程组的消去法 3.2.1 高斯消去法与高斯若当消去法 例1 第一步：先将方程（1）中未知数 的系数2除（1）的两 边，得到下列方程组： 再将第二个方程减去第一个方程的4倍，第三个方程减去 第一个方程的2倍。 第二步：将方程 中第二个方程的两边除以 的系数4 将第三个方程减去第二个方程： 第三步：为了一致期见，将第三个方程中的 系数变为1， 除以 我们来分析一下上述过程：整个过程分两大步。一是用 逐次消去未知数的方法，把原来的方程组化为与其等价的三 角形方程组。用矩阵的观点来看，就是用初等变换的方法将 方程组的系数矩阵进行初等变换，即 这样就将系数阵化为单位三角阵，这个过程称为“消元 过程”。二是解三角形方程组，称为“回代过程”，整个过程 称为“有回代过程的顺序消元法”。 下面我们来讨论一般的解n阶方程组的高斯消去法，且 就矩阵的形式来介绍这种新的过程： 一般地，第k步：即将矩阵变为如下 第n步：得到： 经过上述n步过程后，原系数矩阵A变为一个单位上三角 矩阵，即原方程组化为一个和它完全等价的三角形方程组， 即 高斯消去法： （1）消元过程： 对k=1,2, …, n 依次计算 (2) 回代过程： 例3.1 试用高斯消去法求解线性方程组 消元过程为 解 即把原方程组等价约化为 据之回代解得 为了避免回代的计算，我们可在消元过程中直接把系数矩阵A约化为单位矩阵I，从而得到解，即 相应地，计算公式可表述为： 从而得到解 这一无回代的消去法称为高斯-若当（Jordan）消去法 二、高斯-若当（Jordan）消去法 解 例 2 试用高斯-若当消去法求解例3.1的线性方程组。 因为 解 一般公式： 高斯约当消去法是一个具有消去过程而无回代过程的算法。 以上两种消去法都是沿系数矩阵的主对角线元素 进行的，即第k次消元是用经过前k-1次消元之后的系数阵位于（k,k)位置的元素作除数，这时的(k,k)位置上的元素可能为0或非常小，这就可能引起过程中断或溢出停机。因此： 3.2.2 消去法的可行性和计算工作量 定理 3.1 如果的各阶顺序主子式均不为零，即有 即消去法可行。 推论 若系数矩阵严格对角占优，即有 定理 3.2 求解 n 阶线性方程组 (3-1) 的高斯消去法的乘除工作量约为 ，加减工作量约为 ；而高斯-若当消去法的乘除工作量约为