数值分析第二章解线性方程组的直接方法.ppt

如何求PA=LU中的P, L, U 上机作业 第61页至62页任选一题. 范数等价定理: Rn中任意两个范数等价. ? Rn中范数的重要性质: 范数等价定理 例 1-范数, 2-范数和 ?-范数是两两等价的. ? 当不需要指明使用哪一种向量范数时，就用记号||.|| 泛指任何一种向量范数. ? 有了向量的范数就可以用它来衡量向量的大小和表示向量的误差. ? 设 x 为Ax=b 的精确解，x*为其近似解 ? 绝对误差 ? 相对误差 ? 矩阵范数是用于定义矩阵“大小”的量，类似于向量范数，可以定义 n 阶方阵A的范数. 定义(矩阵范数) 设A为n 阶方阵，按照一定规则有一实数与之对应，记为||A||, 若||A||满足： 1) ||A||?0, ||A||=0当且仅当A=0时； 2) 对任意实数?， ||?A||=|?| ||A||； 3) 对任意两个n阶方阵A, B, 都有 ||A+B|| ? ||A||+||B||; 4) ||AB|| ? ||A|| ||B|| (相容性条件) 则称||A||为矩阵A的范数. ? 矩阵范数 ? 常用的三种矩阵范数 ? 1-范数或列范数 ? 2-范数或谱范数 ? ?-范数或行范数 ? 用矩阵范数的定义来验证, 即验证它们满足矩阵范数定义中的四个性质. ? 三对角矩阵三角分解中LU的求解次序 对 k=3, …, n ? 三对角矩阵三角分解中LU的求解 ? 右端用矩阵乘法展开, 比较两边元素得 乘除运算量 2n－2 ? 解三对角方程组的追赶法 先求Ly=d, 得 y; 再求 Ux=y, 得x. 追：消元过程 赶：回代过程 ? Ly=d ? 乘除运算量为 n－1 ? Ux=y ? 乘除运算量为 2n－1 追赶法总的乘除运算量 5n－4 ? 追赶法的实质就是Gauss消元法, 只是由于系数中出现了大量的零, 在计算过程中将它们撇开, 从而使计算公式大大简化, 也大大减少了计算量. ? 为节省计算机存储单元, 计算得到的 lk, uk 分别存放在 ak, bk 的存储单元内, 而 yk, xk 存放在 dk 的存储单元内. ? 当系数矩阵为满足定理条件的严格对角占优阵时,追赶法具有良好的数值稳定性. §4 平方根法与改进的平方根法 求解对称正定方程组, 即 其中A为对称正定阵. ? 利用对称性, 平方根法的计算量是Gauss消去法的一半. ? 对称正定阵的Cholesky分解 定理 若 A 是对称正定阵, 则存在唯一的非奇异下三角阵L, 使得 且 L的对角元素皆为正数, 即 lii 0 ( i=1, 2, …, n). 证明 ? A可进行Doolittle分解. A为对称正定阵, 它的n个顺序主子式均大于0, 设 其中 为单位下三角阵, U为上三角矩阵. 用反证法 ? 对U进一步分解. ? A对称 ? A正定?对角矩阵D的对角元均为正. 令 则 其中 为非奇异下三角阵, 且对角元均为正. ? A=LLT. 对称 ? A为3阶对称正定阵, A=L LT, 怎样求L? 对称 例 对下列矩阵进行Cholesky分解 例 对下列矩阵进行Cholesky分解 Matlab解法： A=[4 -2 6;-2 17 5;6 5 22] R=chol(A) L=R’ ? A为n阶对称正定阵, A=L LT, 怎样求L? ? L中元素的求解次序 ? 依次求L的第一列, 第二列, …, 第n列. ………... 对称 ? A为n阶对称正定阵, A=L LT, 怎样求L? li 为 L的第 i 个行向量 对称 ? A为n阶对称正定阵, A=L LT, 怎样求L? ? 求 L 的第一列 ? 求 L 的第二列 ? A为n阶对称正定阵, A=L LT, 怎样求L? ? A为n阶对称正定阵, A=L LT, 怎样求L? ? 设已经求得 L 的前k－1列, 现求L 的第 k 列 ( k=3, 4, …, n) Matlab函数：chol 例 求正定阵 的Cholesky分解. 解 例 求正定阵 的Cholesky分解. 解 先求 Ly=b, 得 y, 再求 LTx=y, 得 x. ? 解正定线性方程组的平方根法或Cholesky分解法. 平方根法或Cholesky分解法 ? 设A为对称正定阵 定理 若 A 是对称正定阵, 则存在唯一的单位下三角阵L和对角阵D, 使得 且 D的对角元素皆