平方根法和改进平方根法求解线性方程组例题与程序.pdfVIP

平方根法和改进平方根法求解线性方程组例题 与程序 2、数学原理 1、平方根法解 n阶线性方程组 Ax=b的 choleskly方法也叫 做 平方根法，这里对系数矩阵 A是有要求的，需要 A是对称正定 矩 阵，根据数值分析的相关理论，如果 A对称正定，那么系数矩 阵就 可以被分解为的形式，其中 L是下三角矩阵，将其代入 Ax二 b 中， 可得：进行如下分解：那么就可先计算 y,再计算 x,由于L是下三角 矩阵，是上三角矩阵，这样的计算比直接使用 A计算简 便，同时你 应该也发现了工作量就转移到了矩阵的分解上面，那 么对于对称正 定矩阵 A进行 Cholesky分解，我再描述一下过程 吧：如果你对原理 很清楚那么这一段可以直接跳过的。设，即其 中第 1步，由矩阵乘 法，故求得一般的，设矩阵 L的前 k-l列元 素已经求出第 k步，由 矩阵乘法得于是 2、 改进平方根法在平方根的基础上，为了避免开方运算，所 以用计算；其中，；得按行计算的元素及对元素公式对于、、 计算 出的第行元素后，存放在的第行相置，然后再计算的第行元 素，存 放在的第行、的对角元素存放在的相应位置、对称正定矩 阵按分解 和按分解计算量差不多，但分解不需要开放计算。求解, 的计算