[問題] 用solve計算線性方程錯誤

作者: loser113 (洨大魯蛇ㄍ)   2015-11-15 22:57:59
我用 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在做計算的時候 有可能有這種誤差嗎謝謝
作者: celestialgod (天)   2015-11-15 23:03:00
這種誤差是有可能的
作者: Wush978 (拒看低質媒體)   2015-11-15 23:21:00
10^-17... 可能是浮點數運算的誤差
作者: loser113 (洨大魯蛇ㄍ)   2015-11-15 23:30:00
謝謝有方法解決嗎
作者: celestialgod (天)   2015-11-16 09:18:00
沒有...

Links booklink

Contact Us: admin [ a t ] ucptt.com