作者loser113 (洨大魯蛇ㄍ)
看板R_Language
標題[問題] 計算線性方程錯誤
時間Sun Nov 15 22:57:59 2015
我用 slove(A,b) 計算 A=nxn矩陣 =b 時 我的解不該出現負值
x1~xn >0
n最多到7千多
code檢查多遍無誤
發現在結尾Xn和前幾個Xk 解出數值小於10^-17 (|Xk|<10^-17)
就會有可能是負值 (像Xk=-2.357...*10^-25)
R在做計算的時候 有可能有這種誤差嗎謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 114.42.139.52
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1447599481.A.855.html
1F：→ celestialgod: 這種誤差是有可能的 11/15 23:03
2F：推 Wush978: 10^-17... 可能是浮點數運算的誤差 11/15 23:21
3F：→ loser113: 謝謝 11/15 23:30
4F：→ loser113: 有方法解決嗎 11/15 23:30
5F：→ celestialgod: 沒有... 11/16 09:18