[問題]矩陣做LU分解
如標題
在寫C時遇到一個解不出來的矩陣
於是便想說把它丟到matlab試試看解不解得出來
(如下連結 )
http://ppt.cc/_b_P
而該矩陣為200*200方陣，且有大約98%為0
目前已知det(J)為0
且rank(J)為200
inv(J)也可求出且
由rank=200看似沒問題
但是做LU分解時
[L,U]=lu(J)
發現在下三角矩陣的地方(也就是L)
並不為下三角矩陣(對角線的上方仍有非零值 如下圖)
http://ppt.cc/RAJj
在L(9,9)的上面應該為0，可是出現了-0.2602的值
而在後面的對角線上方，也是有同樣出現不為0的值(似乎不為下三角矩陣!?)
想請問版友們是不是哪一個地方出現了問題(例如special case之類的)
PS：如果是因為我不懂一些限制的關係
例如矩陣本身的限制不能做LU分解或是其他計算什麼的
那我可能就PO錯版了，麻煩請各位鞭小力一點 @@
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.154.186
→
07/19 18:23, , 1F
07/19 18:23, 1F