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

显然这种矩阵的范数||A||v依赖于向量范数||x||v的具体含义. 也就是说, 当给出一种具体的向量范数||x||v时, 相应地就得到了一种矩阵范数||A||v. 定理18 设x?Rn, A aij ?Rn×n, 则有常用的算子范数 （称为A的行范数）, （称为A的列范数）, （称为A的2-范数）. 其中?max ATA 表示ATA的最大特征值, 由于矩阵2-范数与ATA 的特征值有关，所以又称为A的谱范数. 证明见书p202. 由定理18看出, 计算一个矩阵的||A||∞, ||x||1还是比较容易的, 而矩阵的2-范数在||x||2计算上不方便, 但是矩阵的2-范数具有许多好的性质, 它在理论上是非常有用的. 我们指出, 对于复矩阵 即A?Cn×n 定理18中1,2显然也成立, 对于3应改为 例7 设 则 求相应的各种范数. 解 因为 令 得ATA的特征值 故 可见 参见书p203-例7. 定义6 设A?Rn×n的特征值为?i i 1,2,?,n , 称 为A的谱半径. 例子 设 , 求A的谱半径. 解 得A的特征值 故得A的谱半径为 定理19 特征值上界 设A?Rn×n, 则ρ A ≤||A||,即A的谱半径不超过A的任何一种算子范数 对||A||F亦成立 . 证明 设?是A的任一特征值, x为相应的特征向量,则Ax ?x, 由 5.7 得 注意到||x||≠0, 即得 |?|≤||A||. 定理20 如果A?Rn×n为对称矩阵, 则||A||2 ρ A . 证明留作练习题或查有关参考书. 定理21 如果||B|| 1, 则I±B为非奇异矩阵, 且 其中||·||是指矩阵的算子范数. 证明 用反证法. 若det I-B 0, 则 I-B x 0有非零解, 即存在x0≠0使Bx0 x0, ||Bx0||/||x0|| 1, 故||B||≥1, 这与假设矛盾. 又由 I-B I-B -1 I, 有 I-B -1 I+B I-B -1, 从而 类似对加法有 I+B I+B -1 I, I+B -1 I-B I+B -1, 得 求‖A‖1, ‖A‖2 , ‖A‖∞ , ρ A . 例子 设 , 解: 显然‖A‖1 4,‖A‖∞ 4 其中lii 0 i 1,2,?,n . 由矩阵乘法及ljk 0 当j k时 ,得 于是得到解对称正定方程组Ax b的平方根法计算公式 对于 j 1,2,?,n 求解Ax b，即求解两个三角形方程组 1 Ly b，求y； 2 LTx y，求x. 由计算公式知道 所以有 于是 上面分析说明，分析过程中元素ljk的数量级不会增长且对角元素ljj恒为正数. 于是有 结论 不选主元素的平方根法是一个数值稳定的方法. 当求出L的第j列元素时，LT的第j行元素亦算出.所以平方根约需n3/6次乘除法，大约为一般直接LU分解法计算量的一半. 例题 用平方根法求解对称正定方程组 解 首先对A进行Cholesky分解 求解Ly b，得 y1 2, y2 3.5, y3 1. 求解LTx y，得 x1 1, x2 1, x3 1. 由公式 4.7 看出，用平方根法解对称正定方程组时，计算L的元素ljj需要用到开方运算. 为了避免开方，我们下面用定理10的分解式 A LDLT. 即 由矩阵乘法，并注意ljj 1,ljk 0 j k ，得 于是得到L的元素及D的对角元素公式： 为了避免重复计算，我们引进 tij lij dj . 由 4.9 得到按行计算L, T T LD 元素的公式： d1 a11 求解Ly b及DLTx y的计算公式为 计算公式 4.10 , 4.11 称为改进的平方根法. 5.4.3 追赶法 在一些实际问题中，例如解常微分方程边值问题，解热传导方程以及船体数学放样中建立三次样条插值函数，都会要求解系数矩阵为对角占优的三对角线方程组 Ax f，即: 其中, 当|i-j| 1时, aij 0, 且满足如下的对角占优条件: a |b1| |c1| 0; b |bi|≥|ai|+|ci|, ai,ci≠0, i 2,3,…,n-1 . c |bn| |an| 0. 我们利用矩阵的直接三角分解法来推导解三对角线方程组 4.12 的计算公式. 由系数阵A的特点，可以将A分解为两个三角矩阵的乘积，即 A LU, 其中取L下三角阵, 取U为单位上三角阵, 这样求解方程组Ax f的方法称为追赶法. 设 其中 为待定系数，比较 4.13 两边即得 求解Ax f等价于求解两个三角形方程组. 1 Ly f, 求y; 2 Ux y, 求x. 从而得到解三对角线方程组的