k***s 发帖数: 277 | 1 【 以下文字转载自 Mathematics 讨论区 】
发信人: kenjs (乔峰.Leon), 信区: Mathematics
标 题: 一个矩阵分解的问题(对称实矩阵)
发信站: BBS 未名空间站 (Thu Nov 9 18:11:19 2006), 转信
cholesky decomposition (from matlab)
可以把实对称矩阵X分解为
R' * R = X
可是现在我需要的是把X分解为R,R也是upper trianglular matrix
R * R' = X
请问有没有什么现成的方法?(最好不要自己从头到尾coding) |
l***t 发帖数: 81 | 2 what do you mean ?
they are the same.
【在 k***s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 【 以下文字转载自 Mathematics 讨论区 】 : 发信人: kenjs (乔峰.Leon), 信区: Mathematics : 标 题: 一个矩阵分解的问题(对称实矩阵) : 发信站: BBS 未名空间站 (Thu Nov 9 18:11:19 2006), 转信 : cholesky decomposition (from matlab) : 可以把实对称矩阵X分解为 : R' * R = X : 可是现在我需要的是把X分解为R,R也是upper trianglular matrix : R * R' = X : 请问有没有什么现成的方法?(最好不要自己从头到尾coding)
|
k***s 发帖数: 277 | 3 no,
R is upper triangular, then R' is lower triangular
【在 l***t 的大作中提到】![](/moin_static193/solenoid/img/up.png) : what do you mean ? : they are the same.
|
p***o 发帖数: 1252 | 4 transpose X first
【在 k***s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : no, : R is upper triangular, then R' is lower triangular
|
t****t 发帖数: 6806 | 5 X is symmetric. go math board and read answer. this is not programming trick
at all.
【在 p***o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : transpose X first
|