作者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/cn.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