[討論] matlab判斷矩陣rank

[j0958322080 (Tidus)]

我自己有寫個程式去解一個 15 by 15 非稀疏矩陣的rank，

解出來是沒有問題，解到最後對角線元素都是大於一，rank為15。

但我用matlab解的時候他的答案只有12，

直接去解反矩陣會警告我說這個矩陣是deficient。

我的解反矩陣是用 LU 分解，聽說MATLAB也是，

雖然反矩陣的元素有到10^-42次方，

但是不知道為什麼RANK的部分會有差異??

--
！！！！！！！！！！！！！簽名檔破530000點擊率啦！！！！！！！！！！！！！！！
Fw: [問卦] 電影：決勝21點的機率問題  https://goo.gl/2BpbB7  #1MfN3FgZ (joke)
→ yeebon: chx64的1/2悖論真的很經典呢
！！！！！！！！！！！！！！簽名檔破530000點擊率啦！！！！！！！！！！！！！！

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 39.12.192.18
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1514705793.A.99E.html
推 sunev: numerical instability
→ j0958322080: 那matlab是如何算矩陣rank的??
推 sunev: matlab不是open source，所以看他說明文件怎麼寫囉
→ rex0707: 我也曾經遇過自己寫的程式沒問題 但是matlab會跳出警告