数值分析课件第五章4-5.ppt

5.3.2平方根法 应用有限元法解结构力学问题时，最后归结为求解线性方程组，系数矩阵大多具有对称正定性质。所谓平方根法，就是利用对称正定矩阵的三角分解而得到的求解对称正定方程组的一种有效方法，目前在计算机上广泛应用平方根法解此类方程组。 设 为对称矩阵，且 的所有顺序主子式均不为零，由本章定理7知， 可唯一分解为如（3.1）式的形式。 为了利用 的对称性，将 再分解为 其中 为对角矩阵， 为单位上三角矩阵， 于是 又 由分解的唯一性即得 代入（3.6）得到对称矩阵 的分解式 总结上述讨论有下面定理。 定理9(对称阵的三角分解定理)设 为 阶对称矩阵，且 的所有顺序主子式均不为零，则 可唯一分解为 其中 为单位下三角矩阵， 为对角矩阵。 现设 为对称正定矩阵，首先说明 的分解式 中 的对角元素 全身正数。 事实上，由 的对称正定性，定理6的推论成立，即 于是 由定理6得到 其中 为下三角矩阵。 定理10（对称正定矩阵的三角分解或楚列斯基（cholesky）分解）如果 为 阶对称正定矩阵，则存在一个实的非奇异下三角矩阵 使 ，当限定 的对角元素为正时，这种分解是唯一的。 下面我们用直接分解方法来确定计算 元素的递推公式。因为 其中 由矩阵乘法及 （当 时）， 得 于是得到解对称正定方程组 的平方根法计算公式： 对于 求解 ，即求解两个三角方程组： 由计算公式（3.7）知 于是 上面分析说明，分解过程中元素 的数量级不会增长且对角元素 恒为正数。于是不选主元素的平方根法式一个数值稳定的方法。 当求出 的第 列元素时， 的第 列元素亦算出。 所以平方根法约需 次乘除法，大约为一般直接 分解法计算的一半。 由于 为对称矩阵，因此在计算机实现时只需存储 的下三角部分，共需要存储 个元素，可用一维数组存放，即 矩阵元素 一维数组的表示为 ， 的元素存放在 的相应位置。 由公式（3.7）看出，用平方根法解对称正定方程组时，计算 的元素 需要用到平方运算。为了避免开方，我们下面用定理9的分解式 即 由矩阵乘法，并注意 ，得 于是得到计算 的元素及 的对角公式； 对于 为了避免重复计算，我们引进 ， 由（3.9）式得到按行计算 元素的公式： 对于 计算出 的第 行元素 后，存放在 的第 行相应位置，然后再计算 的第 行元素，存放在 的第 行。 的对角元素存放在 的相应位置。例如 对称正定矩阵 按 分解和按 分解计算量差不多，但 分解不需要开方计算。 求解 计算公式： 计算公式（3.10），（3.11）称为改进的平方根法。 5.3.3 追赶法 在一些实际问题中，例如解常微分方程边值问题，解热传导方程以及船体数学放样中建立三次样条函数等，都会要求解系数矩阵为对角占优的三对角线方程组 简记为 。其中，当 时， ，且： 我们利用矩阵的直接三角分解法来推导解三对角线方程组（3.12）的计算公式。由系数矩阵 的特点，可以将 分解为两个三角矩阵的乘积，即 其中 为下三角矩阵， 为单位上三角矩阵。下面我们来说明这种分解式可能的。设 其中 为待定系数。比较分解式（3.13）两边即得 由 ，得 。下面我们用归纳法证明 即 ，从而由（3.14）式可求出 。 （