s*******2
发帖数: 499 | 1 我用R算一个20X20的对称矩阵的逆矩阵。可是我发现算出来的矩阵有问题。因为这个逆
矩阵和逆矩阵的转置差别较大。
可是，如果我用原矩阵乘以逆矩阵，结果基本上是个单位矩阵。
请问这是怎么回事？谢谢。 |
s**c
发帖数: 1247 | 2 应该是数值近似计算的问题
【在 s*******2 的大作中提到】
: 我用R算一个20X20的对称矩阵的逆矩阵。可是我发现算出来的矩阵有问题。因为这个逆
: 矩阵和逆矩阵的转置差别较大。
: 可是，如果我用原矩阵乘以逆矩阵，结果基本上是个单位矩阵。
: 请问这是怎么回事？谢谢。
|
s*******2
发帖数: 499 | 3 那有算得准确的软件吗？matlab怎么样？
【在 s**c 的大作中提到】
: 应该是数值近似计算的问题
|
l*********s
发帖数: 5409 | 4 Make sure your matrix is non-singular. 20 is too small; precision of numeric
computation is unlikely to be the culprit.
【在 s*******2 的大作中提到】
: 我用R算一个20X20的对称矩阵的逆矩阵。可是我发现算出来的矩阵有问题。因为这个逆
: 矩阵和逆矩阵的转置差别较大。
: 可是，如果我用原矩阵乘以逆矩阵，结果基本上是个单位矩阵。
: 请问这是怎么回事？谢谢。
|
s*******2
发帖数: 499 | 5 你说的对，谢谢。
此前矩阵的det特别小，接近于0.
我把这个矩阵乘以一个常数，再求逆矩阵，结果有改善，可是仍然没法彻底解决问题。
numeric
【在 l*********s 的大作中提到】
: Make sure your matrix is non-singular. 20 is too small; precision of numeric
: computation is unlikely to be the culprit.
|
s*r
发帖数: 2757 | 6 不是full rank以后，inverse matrix就不一定要是对称阵了?
numeric
【在 l*********s 的大作中提到】
: Make sure your matrix is non-singular. 20 is too small; precision of numeric
: computation is unlikely to be the culprit.
|
s*****n
发帖数: 2174 | 7 如果没记错的话, matlab算逆矩阵远远好于R.
我本科刚接触R的时候, 最不爽的就是这点.
【在 s*******2 的大作中提到】
: 那有算得准确的软件吗？matlab怎么样？
|
c******r
发帖数: 300 | 8 Try the eigen-decomposition to invert the symmetric matrix
【在 s*******2 的大作中提到】
: 我用R算一个20X20的对称矩阵的逆矩阵。可是我发现算出来的矩阵有问题。因为这个逆
: 矩阵和逆矩阵的转置差别较大。
: 可是，如果我用原矩阵乘以逆矩阵，结果基本上是个单位矩阵。
: 请问这是怎么回事？谢谢。
|
