【数值分析】choleskey分解，matlab实现

最新推荐文章于 2024-11-03 14:00:00 发布

你哥同学

最新推荐文章于 2024-11-03 14:00:00 发布

阅读量1.1k

点赞数 20

分类专栏：数值分析文章标签： matlab 开发语言 choleskey分解乔列斯基分解

本文链接：https://blog.csdn.net/lafea/article/details/135321790

版权

数值分析专栏收录该内容

21 篇文章 16 订阅

订阅专栏

平方根分解（Choleskey分解）

$A=GG^ \mathrm T \,\,,\,\, A对称正定$
$\begin{align*} A =LDM= LDL^ \mathrm T=(LD^{1/2})(LD^{1/2})^ \mathrm T=GG^ \mathrm T \end{align*}$
$\begin{cases} Gy=b \\ \\ G^ \mathrm Tx=y \end{cases}$
手算的话根据转置的性质直接把 ${G}$ 矩阵写出来就行，非常好写，编程直接参考同济《现代数值计算》算法2.2.3。
matlab编程实现：

%% Choleskey分解解线性方程组
function [x,G] = cholesSolve(A,b)
    n = size(A);
    for i = 1:n
        t = 0;
        for s = 1:i-1
            t = t+ G(i,s)^2;
        end
        G(i,i) = sqrt(A(i,i)-t);
        for k = i+1:n
            t = 0;
            for s = 1:i-1
                t = t+G(i,s)*G(k,s);
            end
            G(k,i) = (A(k,i)-t)/G(i,i);
        end
    end
    % 回代
    for i = 1:n
        t = 0;
        for j = 1:i-1
            t = t+G(i,j)*y(j);
        end
        y(i)=(b(i)-t)/G(i,i);
    end
    for i = n:-1:1
        t = 0;
        for j = i+1:n
            t = t+G(j,i)*x(j);
        end
        x(i) = (y(i)-t)/G(i,i);
    end
    x = x';
end