[問題]矩陣做LU分解

[smartfayfay (狒狒)]

如標題

在寫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
→ smartfayfay:對了 一開始目地是要解聯立方程組