电子科大MATLAB第5节 线性方程组求解的直接法.pptVIP

正定对称矩阵 若A为n阶满秩方阵，则有 A T A为对称正定矩阵 1、 A T A为对称矩阵 2、 根据正定矩阵定义： 由于A满秩，即方程 A x = 0 只有0解，所以任意任意非零向量有： 对称正定矩阵Cholesky分解： Cholesky 分解：A 为对称正定矩阵，则有： 令： 则： A 对称，则 上三角 下三角 所以有： 需要引理 对称正定矩阵Cholesky分解 说明：上一页的证明过程是一个构造式的证明方法，不但得到了结论，也给出了LU分解和Cholesky矩阵分解的关系： 已知LU分解求cholesky分解： C = D-1U， 已知cholesky分解，求LU分解： U = DC， L = CTD-1 。 如何计算D？ U对角线元素开方，或C的对角线元素。 优势：cholesky方法误差可控，不需选主元。 对称正定矩阵Cholesky分解 例： 对称正定矩阵Cholesky分解 clear all; close all; clc % CHOL分解 display(CHOL分解) A = [4, 1, 0; 1, 4, 1; 0, 1, 4] display(C1 = chol(A)) C1 = chol(A) display([L, U] = lu(A)) [L, U] = lu(A) for iii = 1 : 3 for jjj = 1 :