第5章解线性方程组的直接方法.pptVIP

显然这种矩阵的范数||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||∞, ||A||1还是比较容易的, 而矩阵的2-范数在||A||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 其中lii0(i=1,2,?,n). 由矩阵乘法及ljk=0(当jk时),得 于是得到解对称正定方程组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(jk)，得 于是得到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|