第三章 线性代数方程组的直接解法分析.pptVIP

（保带状结构定理） 设 为上半带宽为 ，下半带宽为 的带状矩阵， 且其顺序主子式 ，则 有唯一的三角分解 ， 其中 是下半带宽 为 的单位下三角阵， 是上半带宽为 的上三角阵。 证明可根据前面讲过的三角分解公式 保带状结构定理说明：矩阵的三角分解中， 和 带外元素为零，因此不必计算，且不必参加求和运算 ? 三对角线性方程组的三对角算法（追赶法） 三对角线性方程组 其中 根据保带状结构定理，系数矩阵可作如下三角分解： ?三对角矩阵 分解的计算公式： ?方程组求解的计算公式： ?解方程组 ?解方程组 “追”的过程 “赶”的过程 ?追赶法实现的一个充分条件(补充) 设 为前述三对角矩阵，且满足下列条件： ? ? 则 非奇异，且 特殊情况：如果三对角矩阵 为严格对角占优矩阵，则可以采用追赶法求解。 例2：用追赶法求解三对角方程组 , 其中: 解： 注意到本例并不满足定理3.1.3’的条件,但仍然可以利用追赶法来求解.因此,定理3.1.3’的条件仅是充分条件. 求解方程组 求解方程组 §3.2 选主元三角分解 ?选主元三角分解的思想 三角分解过程中存在的问题 ?Gauss消元法完成的条件是矩阵的各阶顺序主子式 (n=1,2,…,n-1)均不为零. ?三角分解过程中的除法运算要求分母不 能太小,否则 将可能产生不稳定情况. 选主元的目的就是为了完成消元且避免不稳定情况的发生 例3：在8位制计算机上解方程组 要求用三角分解方法计算。 8个 解： 小主元 可能导致计算失败 交换方程组的两行 8个 ? Gauss全主元三角分解法 交换单位矩阵 的第 列(行)和第 列(行)得到的矩 (初等置换矩阵) 阵 ,称之为初等置换矩阵. 列 列 Step 1(k=1)：第1步选择主元 寻求 和 满足 然后交换矩阵 的第 行和 行，第 列和 列 设给定 阶矩阵 记 然后按照前面讨论的方法进行三角分解. 用矩阵表示: 其中, 为初等置换矩阵. * 在没有舍入误差的情况下，经过有限次 运算可以得到方程组的精确解的方法。 第三章 线性方程组的直接解法 /*Direct Method for Solving Linear Systems*/ 求解 Cramer法则: 所需乘除法的运算量大约为(n+1)!+n n=20时，每秒1亿次运算速度的计算机要算30多万年！ 直接法 §3.1 三角形方程组和三角分解 一、 三角形方程组的解法 ?考虑下三角形方程组 的计算公式为： 算法3.1.1 下三角形方程组的前代法： for end ?考虑上三角形方程组 的计算公式为： 算法3.1.2 上三角形方程组的回代法： for end 两种算法的工作量(加减乘除运算次数之和)均为 ?三角分解法的基本思想: 记 方程组可化为下面两个易求解的三角方程组 设已知方程组系数矩阵的三角分解为 其中, 为下三角矩阵, 为上三角矩阵. 二、 高斯(Gauss)变换 取下三角形矩阵 则 可表示为 其中为 单位矩阵, 称下三角形阵 为高斯(Gauss)变换, 为高斯向量. ?Gauss变换的定义 ?高斯(Gauss)变换的性质 性质1 设向量 且 则存在唯一的下三角阵 ,满足 证明： 寻找满足条件的初等下三角阵 记 写成分量形式： 唯一确定 性质2 性质3 性质4 若记 ,则有 即单位下三角阵 可以分解为一系列初等下三角阵的乘积 三、 三角分解的计算 ? Gauss消去法 设给定矩阵 取Gauss变换矩阵 则有 再取Gauss变换矩阵 其中 设给定 阶矩阵 记 ? Gauss消去法的矩阵表示 令 Step 1：如果 高斯变换 取 其中 记 类似地，对 中的 部分重复以上做法 Step k：第k步消元过程的计算公式 整个消元过程的矩阵表示 上三角矩阵 计算 经过n-1次消元，并将 存放在矩阵零元素位置 for for for Gauss消去法的消元过程算法 Gauss消去法工作量为 三角分解的计算过程: Step1 Step2 Step3 S