[问题] 用solve计算线性方程错误 loser113 PTT批踢踢实业坊

[问题] 用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

谢谢有方法解决吗

没有...

继续阅读

Re: [问题] 资料转换celestialgod Re: [问题] 资料转换obarisk [问题] 改善这只龟速的 image filtering 程式Neisseria Re: [问题] 资料转换OGCKM [问题] 资料转换Rose9305 [问题] function的输出laputaca Re: [问题] 除了RStudio 还有其他的IDE?Wush978 Re: [问题] 除了RStudio 还有其他的IDE?Neisseria [问题] 除了RStudio 还有其他的IDE?rn940111 [心得] 中文编码处理Wush978