2.5平方根法追赶法.pdfVIP

§5 平方根法 一、教学设计 1．教学内容：对称正定矩阵的 Cholesky 分解法、三对角线 矩阵分解的追赶法。 2 ．重点难点： Cholesky 分解法、追赶法。 3 ．教学目标：掌握对称正定矩阵的 Cholesky 分解的计算过 程，掌握三对角线矩阵分解的追赶法。 4 ．教学方法：讲授与讨论。 二、教学过程 §5 平方根法 在工程计算中，常遇到求解解对称再正定线性方程组问题， 如应用有限元法解结构力学问题，应用差分方法解椭圆型偏 微分方程等，最后都归结为求解系数矩阵为对称正定阵的线 性方程组。根据系数矩阵的特殊性，是否有更好的解决方案 （在存贮空间上的好处是显而易见的） ，算法上是否有所简 化？ 5 －0 对称正定矩阵及性质复习 定义 ：设 A Rn n ，如果 A 满足条件 （1） AT A ；（2 ）对任意非零向量 x Rn ，有 xT Ax 0 ，则称 A 为对称正定矩阵。 定理 1 （对称正定矩阵的性质）如果 A Rn n 为对称正定矩 阵，则 （1）A 为非奇异阵，且 A 1 亦是对称正定阵； （2 ）记 Ak 为 A 的顺序主子阵，则 Ak 亦是对称正定阵 (k 1,2, , n) ； （3 ）A 的特征值 i ( A ) 0 (i 1,2, , n) ； （4 ）A 的顺序主子式都大于零，即 det( Ak ) 0 (k 1,2, , n) 。 定理 2 设 A Rn n 为对称矩阵（判据） （1）若 A 的特征值 i ( A ) 0 (i 1,2, , n) ，则 A 为对称正定 矩阵； （2 ）若 A 的顺序主子式都大于零， 即 det( Ak ) 0 (k 1,2, ,n) ， 则 A 为对称正定阵。 5 －1 对称正定矩阵的三角分解 由 前 述 定 理 3.1 知 ， 若 n 阶 方 阵 A 的 顺 序 主 子 式 d e tA(k ) ( k 1,2, , n 1) 均不为零，则 A 有唯一的三角分解 A LU ，其中 L 为单位下三角阵， U 为上三角阵。 n 阶对称正 定阵 A 的顺序主子式都大于零，当然有 LU 分解，进一步地， 此 时 L,U 之 间 有 什 么 关 系 ？ 这 对 解 方 程 组 有 用 处 。 由 T T T A U L A LU 及分解的唯一性，想到若 U 的主对角元素 皆为 1，就有可能获得一些结果。为此，再将 U 分解 u12 u1n 1 u u u u 11 12 1n 11 u u 11 11 u22 u2 n u22 u2n U