第三章解线性代数方程组得直接法讲述.pptVIP

越大，解的相对误差可能越大，A对求 解线性方程组来说就越可能呈现病态． 反之， 越小，解的相对误差越小， 越呈现良态． 下面给出两个与条件数有关的定理 定理 设 Ax = b，A为非奇异矩阵，b为非零 向量且 A 和b均有扰动． ，则 使得 非常小, 若A的扰动 定理 设Ax=b，A为非奇异矩阵，b为非零向量，则方程组近似解 的事后估计式为 其中 称为近似解 的余量，简称余量。 因为 * * 因为 * 3.2 直接三角分解法 矩阵的三角分解（LU分解） 高斯消元法的矩阵形式 为上三角矩阵 其中 为单位下三角 矩阵分解为单位下三角 和上三角矩阵的乘积 Doolittle分解法 。 杜利特尔分解 (Doolittle) 常用的两种三角分解 克洛特分解 (Crout) ?利用矩阵 分解求解方程组的思想 记 方程组可化为下面两个易求解的三角方程组 下三角阵L和上三角阵U的计算。由 根据矩阵乘法，比较等式两边可知： 这样，U的第一行和L的第一列被求出。 当k=2,3,…,n 时，有 由此得到LU分解的计算公式： 对k=1,2,…,n, 计算 计算中，上述两个公式交替使用，即从U的第一行和L的第一列开始，计算共分步，每步先计算U的一行，再计算L的相应一列，最后可求出L和U的全部元素，从而实现A的LU分解. 最后，求解两个三角形方程组。 求解Ux=y： 求解Ly=b： 3.3 解三对角方程组的三对角算法（追赶法） 特别的，如果 A 是三对角线性方程组 ，即 ?三对角矩阵 分解的计算公式： ?方程组求解的计算公式： ?解方程组 ?解方程组 “追”的过程 “赶”的过程 ?追赶法实现的一个充分条件 定理3.2 设 为前述三对角矩阵，且满足下列条件： ? ? 则 非奇异，且 例5：用追赶法求解三对角方程组 , 其中: 解： 注意到本例并不满足定理3.2的条件,但仍然可以 利用追赶法来求解.因此,定理3.2的条件仅是充分 条件. 求解方程组 求解方程组 3.4 平方根法和改进的平方根法 当A为对称正定矩阵时，A=LU分解将有更简单的对称形式，并且分解过程无需选主元，有良好的数值稳定性。 一、平方根法 定理3.3 对任意n阶对称正定矩阵 A，均存在下三角矩阵L， 使 A=LLT （1） 成立。进而, 如果规定 L 的对角元为正数，则 L 是唯一确定的，称（1）式为对称正定矩阵A的Cholesky分解. 设 = 由矩阵乘法， A的下三角元素 可由下式表出： 由上式，自左向右逐列计算L中待定元素lij ,公式如下 对于 j= 1, 2,…, n 计算 ，i = j+1, j+2,…,n 计算次序为 这样， 求解Ly=b 求解LTx = y 由 得 由此推出 因此在分解过程中L的元素的数量级不会增长，故平方根法通常是数值稳定的，不必选主元。 , k=1,2,…,j． 3.5 条件数与方程组的性态 一、 向量范数 ?非负性： ?齐次性： ?三角不等性： 且 则称 为 中向量 的范数。 非负实值函数 存在唯一实数 与之对应，且满足 定义：设 是 的一个映射，若对 ?常用的几种向量范数： 设 ? 1-范数： ? 2-范数： ? -范数： 上述3种向量范数统称为P-范数 例：求向量 的１，2和∞-范数。 解： （向量范数的等价性定理）设 和 为 上的任意两种向量范数，则存在两个与向量无关的正常数 c1和c2，使得下面的不等式成立 并称 和 为 上的等价范数。 定理 二、 矩阵范数 ?非负性： ?齐次性： ?三角不等性： 且 ? 定义：设 是 的一个映射，若对 ，存在唯一实数 与之对应，且满足 则称 为 中矩阵 的范数。 ?列范数： 记 ?行范数： ?谱范数： 其中 是 的最大特征值 谱半径 ?常用的几种矩阵范数： 设 求 例2 解： 引入条件数的原因 定义 如果线性方程组Ax=b中，A或b的元素的微小变化，就会引起方程组解的巨大变化，则称方程组为“病态”方程组，矩阵A称为“病态”矩阵．否则称方程组为“良态”方程组，矩阵A称为“良态”矩阵。 我们需要一种能刻画矩阵和方程组“病态”标准的量。 三、条件数