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

看板R_Language作者 (洨大魯蛇ㄍ)時間10年前 (2015/11/15 22:58), 編輯推噓1(104)
留言5則, 3人參與, 最新討論串1/1
我用 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://www.ptt.cc/bbs/R_Language/M.1447599481.A.855.html

11/15 23:03, , 1F
這種誤差是有可能的
11/15 23:03, 1F

11/15 23:21, , 2F
10^-17... 可能是浮點數運算的誤差
11/15 23:21, 2F

11/15 23:30, , 3F
謝謝
11/15 23:30, 3F

11/15 23:30, , 4F
有方法解決嗎
11/15 23:30, 4F

11/16 09:18, , 5F
沒有...
11/16 09:18, 5F
文章代碼(AID): #1MI9rvXL (R_Language)
文章代碼(AID): #1MI9rvXL (R_Language)