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
全站熱搜
留言列表