数值剖析ex12—13“数值剖析”习题课II.pptVIP

* 高斯消元法 矩阵的三角分解 雅可比迭代与赛德尔迭代 迭代法收敛定理 最速下降法 《数值分析》习题课 II Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 2/20 一、高斯消元法 = = 三角方程组解法、顺序消元法、列主元法、追赶法 二、矩阵的三角分解 矩阵的紧凑格式分解、改进平方根法 三、向量范数和矩阵范数 常用的三种向量范数、常用的三种矩阵范数、条件数 四、迭代法及收敛性分析 雅可比迭代法、高斯-赛德尔迭代、收敛定理、误差定理、初等变分原理 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 定理3.1 约化主元ak+1,k+1(k) ≠ 0 (k=0,1,···,n-1)的充分必要条件是 矩阵A的各阶顺序主子式不为零. Ex1.如果A是严格主对角占优矩阵, 则 det(A) ≠0. 证: 用反证法。设det(A) = 0, 则齐次方程组Ax=0有非零解 u =[u1, u2, ···, un ]T. 设 考虑Au =0的第k个等式 ? 3/20 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 4/20 两边约去 |uk|，得 这与主对角占优矛盾, 故det(A) ≠0。 Ex2.设A对称且a11≠ 0,经过高斯消元法一步后,A约化为 证明A2 也是对称矩阵。 证明:设 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 经高斯消元一步后,得 5/20 所以, A2 = A2T 思考: 1.若A是对称正定矩阵,经高斯消元一步后,右下角子矩阵A2也是对称正定矩阵; 2.若A为对角占优矩阵，经过高斯消元法一步后，右下角子矩阵A2也是对角占优矩阵。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Ex3.对任何一种矩阵的算子范数,证明矩阵A的谱半径与A的范数有关系:ρ(A) ≤ || A || 证:设 ? 是矩阵A任一特征值,x 是对应的特征向量,则 ? ? ? Ex4.若矩阵A是n阶对称矩阵, 则有 证：设? 是A的任一特征值,由于A对称,故?2 是矩阵ATA的特征值,即 6/20 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 7/20 由2-范数计算公式 Ex5.对任意x，y∈Rn，利用向量范数的三角形不等式证明： 证: || x || = || (x – y )+ y || ≤|| x – y || + || y || ? || x || – || y || ≤|| x – y || 同理, || y || – || x || ≤|| y – x || =|| x – y || ? || x || – || y ||≥ – || x – y || ? – || x – y || ≤ || x || – || y || ≤|| x – y || ? Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Jacobi 迭代法的迭代矩阵 8/20 Gauss-Seidel迭代法的矩阵: BG-S= (D – L)-1U Ax = b, 将矩阵分裂: A = D – U – L BJ = D-1(U