11讲：ch4-3平方根法.pptVIP

第4节 上一讲 要 点 回 顾 LU分解法—直接三角分解法 上一讲 要 点 回 顾 消元法的过程解析 上一讲 要 点 回 顾 令 故消元过程是把系数矩阵A分解成下三角矩阵与上三角矩阵的乘积的过程—LU分解 上一讲 要 点 回 顾 假定Ax=b AX=（LU)x=L( U x ) = b 令 U x=y,则原线性方程组 Ax=b LU分解的应用 上一讲 要 点 回 顾 三对角线性方程组： 四、追赶法与三对角线性方程组 定理 设矩阵A满足下列条件： 则它可分解为 追赶法与三对角线性方程组 追赶法与三对角线性方程组 平方根法与改进的平方根法 第4 节 假设线性方程组: 其系数矩阵A对称正定，则A的各阶顺序主子式和 全部特征值均大于零 （对称正定矩阵的LU分解形式更加简单，平方根法就是针对正定矩阵的LU分解法） 对角线元素未必是1 定理: 设A是对称正定矩阵,则存在唯一的非奇异 下三角阵L,使得 且L的对角元素皆为正. 平方根法 证明： 因A对称正定，其各阶顺序主子式均大于零， 设 A=LU， 令 由LU分解的唯一性，得 对任意非零向量x， 存在性（构造性证明） 定理证明 于是由A正定得： 即D正定，所以D的对角元素均为正数。 定理证明 平方根(Cholesky分解法）法 由 其中 由矩阵乘法及 (当jk时), 得 也可通过紧凑格式来观察公式 平方根(Cholesky分解法）法 这里规定 计算顺序是按列进行, 即 当矩阵A完成平方根分解后,求解Ax=b, 即求解两个 三角形方程组. 平方根(Cholesky分解法）法 比较 推 论 改进平方根法 其中 由比较得 对于 改进平方根法 求解线性方程组的这一方法称为改进的平方根法。 求解 计算公式 对于正定矩阵A按 分解和按 分解计算量 差不多, 但 分解不需要开方计算. 改进平方根法 改进平方根法 故最后解方程组 可求出方程组的解。 事实上， 即为L、U分解中的U，故改进的平方根法也有紧凑格式 改进平方根法 由公式 可得紧凑格式方法 顺序：先求行后求列 改进平方根法 例1、用改进的平方根法解方程组 解： 改进平方根法举例 改进平方根法举例 L={{1,0,0},{2,1,0},{1,-2,1}} DD=DiagonalMatrix[{1,1,9}] L.DD.Transpose[L]; MatrixForm[%] 验证 改进平方根法举例 由公式得： 改进平方根法举例 例2用改进的平方根求解方程组 (与LU分解法基本相同) 解：系数矩阵为对称正定矩阵(用紧凑格式计算) 改进平方根法举例 改进平方根法举例 同解方程组为： 改进平方根法举例 改进平方根法举例