如标题 在写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