.
(2)设对称正定阵系数阵线方程组
二、数学原理
1、平方根法
解 n 阶线性方程组 Ax=b 的 choleskly 方法也叫做平方根法,这里对系数
矩阵 A 是有要求的,需要 A 是对称正定矩阵,根据数值分析的相关理论,如果
A 对称正定,那么系数矩阵就可以被分解为的 形式,其中 L 是下三角
矩阵,将其代入 Ax=b 中,可得:
进行如下分解:
那么就可先计算 y,再计算 x,由于 L 是下三角矩阵,是 上三角矩阵,这
样的计算比直接使用 A 计算简便,同时你应该也发现了工作量就转移到了矩阵
的分解上面,
那么对于对称正定矩阵 A 进行 Cholesky 分解,我再描述一下过程吧:
如果你对原理很清楚那么这一段可以直接跳过的。
设 ,即
其中
第 1 步,由矩阵乘法, 故求得
.页脚.