对称5-对称矩阵的Cholesky分解 - 图1对称5-对称矩阵的Cholesky分解 - 图2阶实对称矩阵。设(正交)矩阵对称5-对称矩阵的Cholesky分解 - 图3使得对称5-对称矩阵的Cholesky分解 - 图4#card=math&code=P%5E%5Ctop%20A%20P%20%3D%20P%5E%7B-1%7DAP%20%3D%20%5CLambda_A%3D%5Cmathrm%7Bdiag%7D%28a_1%2C%20%5Cldots%2C%20a_n%29)是对角矩阵。

    定义矩阵对称5-对称矩阵的Cholesky分解 - 图5#card=math&code=%5Csqrt%7B%5CLambda_A%7D%3D%5CLambda_A%5E%7B1%2F2%7D%20%3D%20%5Cmathrm%7Bdiag%7D%28%5Csqrt%7Ba_1%7D%2C%20%5Cldots%2C%20%5Csqrt%7Ba_n%7D%29). 其中如果某个对称5-对称矩阵的Cholesky分解 - 图6, 则对称5-对称矩阵的Cholesky分解 - 图7.

    则有对称5-对称矩阵的Cholesky分解 - 图8%5E2%20%3D%20(P%5Csqrt%7B%5CLambda_A%7D%20%5C%2C%20P%5E%7B-1%7D)%5E2%20%3D%20P%5CLambda_A%20P%5E%7B-1%7D%3DA#card=math&code=%28P%5Csqrt%7B%5CLambda_A%7D%5C%2CP%5E%5Ctop%29%5E2%20%3D%20%28P%5Csqrt%7B%5CLambda_A%7D%20%5C%2C%20P%5E%7B-1%7D%29%5E2%20%3D%20P%5CLambda_A%20P%5E%7B-1%7D%3DA). 通常记矩阵对称5-对称矩阵的Cholesky分解 - 图9.

    性质1:如果对称5-对称矩阵的Cholesky分解 - 图10半正定,即所有特征值对称5-对称矩阵的Cholesky分解 - 图11, 则对称5-对称矩阵的Cholesky分解 - 图12也是半正定。

    虽然对称5-对称矩阵的Cholesky分解 - 图13%5E%5Ctop%20%3D%20P%5Csqrt%7B%5CLambda_A%7D%20P%5E%5Ctop#card=math&code=%28P%5Csqrt%7B%5CLambda_A%7D%20P%5E%5Ctop%29%5E%5Ctop%20%3D%20P%5Csqrt%7B%5CLambda_A%7D%20P%5E%5Ctop), 但这个矩阵”不够简单”。我们想找一个下三角矩阵对称5-对称矩阵的Cholesky分解 - 图14, 使得对称5-对称矩阵的Cholesky分解 - 图15.

    对称5-对称矩阵的Cholesky分解 - 图16对称5-对称矩阵的Cholesky分解 - 图17对称5-对称矩阵的Cholesky分解 - 图18分解。这里对称5-对称矩阵的Cholesky分解 - 图19是下三角矩阵,对角元都是对称5-对称矩阵的Cholesky分解 - 图20. 记对称5-对称矩阵的Cholesky分解 - 图21对称5-对称矩阵的Cholesky分解 - 图22的转置矩阵。则对称5-对称矩阵的Cholesky分解 - 图23 左边是对称的,右边是上三角的。这只能是等号两边的矩阵都是对角矩阵。令对称5-对称矩阵的Cholesky分解 - 图24, 则对称5-对称矩阵的Cholesky分解 - 图25. 这称为实对称矩阵的对称5-对称矩阵的Cholesky分解 - 图26分解。

    注:电子计算机在计算(正定矩阵)的对称5-对称矩阵的Cholesky分解 - 图27分解时,是将对称5-对称矩阵的Cholesky分解 - 图28的严格下三角元素存储在对称5-对称矩阵的Cholesky分解 - 图29的对应位置上,而将对称5-对称矩阵的Cholesky分解 - 图30的对角元存储在对称5-对称矩阵的Cholesky分解 - 图31的对应的对角位置上。

    如果令对称5-对称矩阵的Cholesky分解 - 图32, 这是下三角矩阵。则对称5-对称矩阵的Cholesky分解 - 图33, 由此就有对称5-对称矩阵的Cholesky分解 - 图34. 分解对称5-对称矩阵的Cholesky分解 - 图35称为实对称矩阵对称5-对称矩阵的Cholesky分解 - 图36的Cholesky分解。通常只考虑对称5-对称矩阵的Cholesky分解 - 图37是对称正定的,此时对称5-对称矩阵的Cholesky分解 - 图38都是实矩阵。

    注:使用Cholesky分解,当求解线性方程组对称5-对称矩阵的Cholesky分解 - 图39时,记对称5-对称矩阵的Cholesky分解 - 图40, 则只需求解对称5-对称矩阵的Cholesky分解 - 图41. 由于对称5-对称矩阵的Cholesky分解 - 图42是下三角矩阵,此时用forward substituion (即顺次求出对称5-对称矩阵的Cholesky分解 - 图43)就能求解。然后再用back substitution (即顺次求出对称5-对称矩阵的Cholesky分解 - 图44)就能解出原方程组。

    性质2:设对称5-对称矩阵的Cholesky分解 - 图45是实对称矩阵的Cholesky分解,又设对称5-对称矩阵的Cholesky分解 - 图46是奇异值(SVD)分解。则对称5-对称矩阵的Cholesky分解 - 图47.

    这是因为:对称5-对称矩阵的Cholesky分解 - 图48%5E2%3D%20(U%5E%5Ctop%20%5CSigma%20U)%5E2%20%3D%20U%5E%5Ctop%20%5CSigma%5E2%20U%20%3D%20U%5E%5Ctop%20%5CSigma%20V%20V%5E%5Ctop%20%5CSigma%20U%20%3D%20GG%5E%5Ctop%20%3D%20A#card=math&code=%28%5Csqrt%7BA%7D%29%5E2%3D%20%28U%5E%5Ctop%20%5CSigma%20U%29%5E2%20%3D%20U%5E%5Ctop%20%5CSigma%5E2%20U%20%3D%20U%5E%5Ctop%20%5CSigma%20V%20V%5E%5Ctop%20%5CSigma%20U%20%3D%20GG%5E%5Ctop%20%3D%20A).

    例子:矩阵对称5-对称矩阵的Cholesky分解 - 图49对称5-对称矩阵的Cholesky分解 - 图50分解是

    对称5-对称矩阵的Cholesky分解 - 图51%20%5Cbegin%7Bpmatrix%7D%201%20%26%200%20%26%201%20%5C%5C%200%20%26%201%20%26%20-2%20%5C%5C%200%20%26%200%20%26%201%5Cend%7Bpmatrix%7D%3B%0A#card=math&code=A%20%3D%20%5Cbegin%7Bpmatrix%7D%201%20%26%200%20%26%200%20%5C%5C%200%20%26%201%20%26%200%20%5C%5C%201%20%26%20-2%20%26%201%5Cend%7Bpmatrix%7D%5Cmathrm%7Bdiag%7D%283%2C%201%2C1%29%20%5Cbegin%7Bpmatrix%7D%201%20%26%200%20%26%201%20%5C%5C%200%20%26%201%20%26%20-2%20%5C%5C%200%20%26%200%20%26%201%5Cend%7Bpmatrix%7D%3B%0A)

    其Cholesky分解是

    对称5-对称矩阵的Cholesky分解 - 图52

    值得注意的是,矩阵对称5-对称矩阵的Cholesky分解 - 图53的特征多项式为对称5-对称矩阵的Cholesky分解 - 图54%20%3D%20%5Clambda%5E3-12%5Clambda%5E2%2B22%5Clambda%20-3#card=math&code=%5Cchi_A%28%5Clambda%29%20%3D%20%5Clambda%5E3-12%5Clambda%5E2%2B22%5Clambda%20-3), 这个多项式的根很难计算。不过从对称5-对称矩阵的Cholesky分解 - 图55对称5-对称矩阵的Cholesky分解 - 图56分解中的矩阵对称5-对称矩阵的Cholesky分解 - 图57, 马上可以看到对称5-对称矩阵的Cholesky分解 - 图58是正定矩阵。