close

參照wiki和課本

裡面的LDLt和Cholesky

LDLt Factorization
先算出矩陣A的上三角和下三角
L  就是 lower triangular
Lt就是 lower triangular的transpose
D就是 U的對角那一行

A=L*D*Lt

Cholesky Factorization
http://ccjou.wordpress.com/2010/09/16/cholesky-%E5%88%86%E8%A7%A3/

和上面那很像 不過他是 newL * newLt

這裡的newL是等於是上面 L* (D)^1/2   (L 乘 D開根號)

然後在用新算出來的newL去做 transpose  = newLt

A=newL*newLt

arrow
arrow
    全站熱搜
    創作者介紹
    創作者 小穎 的頭像
    小穎

    屬於自己世界 ﹋ 穎〃    Black×Sky,,

    小穎 發表在 痞客邦 留言(0) 人氣()