《2.5平方根法追赶法.docVIP

§5　平方根法　 一、教学设计 1．教学内容：对称正定矩阵的Cholesky分解法、三对角线矩阵分解的追赶法。 2．重点难点：Cholesky分解法、追赶法。 3．教学目标：掌握对称正定矩阵的Cholesky分解的计算过程，掌握三对角线矩阵分解的追赶法。 4．教学方法：讲授与讨论。 二、教学过程 §5　平方根法 在工程计算中，常遇到求解解对称再正定线性方程组问题，如应用有限元法解结构力学问题，应用差分方法解椭圆型偏微分方程等，最后都归结为求解系数矩阵为对称正定阵的线性方程组。根据系数矩阵的特殊性，是否有更好的解决方案（在存贮空间上的好处是显而易见的），算法上是否有所简化？ 5－0对称正定矩阵及性质复习 定义：设，如果满足条件 （1）；（2）对任意非零向量，有，则称A为对称正定矩阵。 定理1　（对称正定矩阵的性质）如果为对称正定矩阵，则 （1）A为非奇异阵，且亦是对称正定阵； （2）记为A的顺序主子阵，则亦是对称正定阵； （3）A的特征值； （4）A的顺序主子式都大于零，即。 定理2　设为对称矩阵（判据） （1）若A的特征值，则A为对称正定矩阵； （2）若A的顺序主子式都大于零，即，则A为对称正定阵。 5－1　对称正定矩阵的三角分解 由前述定理3.1知，若阶方阵A的顺序主子式均不为零，则A有唯一的三角分解，其中为单位下三角阵，为上三角阵。阶对称正定阵A的顺序主子式都大于零，当然有分解，进一步地，此时之间有什么关系？这对解方程组有用处。由及分解的唯一性，想到若的主对角元素皆为1，就有可能获得一些结果。为此，再将分解 易知　（用分别记矩阵的阶顺序主子阵，容易验证于是 ） 于是，所以 ， 即 由分解的唯一性知：，，于是 自然地，若记 则，其中是对角元为正数的下三角阵。 定理5.1　（Cholesky分解） 设A是阶对称正定矩阵，则存在唯一的对角线元素全是正数的下三角形矩阵，使得。称这种分解为Cholesky分解。 有了这种分解后，解线性方程组等价于解以下两个三角方程组，，这将带来一些简便。下面讨论如何计算的元素。 5－2　平方根法 设阶对称正定矩阵A有如下Cholesky分解 比较等式两边的第1列对应元素，得到 于是得到的第1列元素 比较等式两边的第2列（对角元及其以下）对应元素，得到 于是 得到的第2列元素 一般地，比较等式两边的第列（对角元以下）对应元素，得到 于是得到的第列元素 5.8 5.9 电算时，可只存A的下三角部分于二维数组A中，且注意到当计算出L的某一列元素后，A的相应列元素在后续计算过程中不再起作用，故可将求得的L的元素保存在A的相应位置上。假设已得到L的前列元素，从下图中注意观察L的第列元素的计算过程。 完成矩阵A的Cholesky分解后，解，化为解以下两个三角形方程组 （1），　 （5.10） （2），　 （5.11） 称由公式（5.8）~（5.11）解对称正定线性方程组的方法为平方根法。 例：用平方根法解对称正定方程组 A= 4 -1 1 -1 17/4 11/4 1 11/4 7/2 计算L的第1列 2 0 0 -1/2 17/4 0 1/2 11/4 7/2 计算L的第2列 2 0 0 -1/2 2 0 1/2 3/2 7/2 计算L的第3列 2 0 0 -1/2 2 0 1/2 3/2 1 Ly=b 2 0 0 0 -1/2 2 0 1 1/2 3/2 1 0