[问题] 方程式准确性

楼主: a29788685 (嘉峰)   2015-11-11 19:49:48
是这样子的:
我有一个方程式,很丑,大致上像这样
Func[a,b,c] = (abc)/(a-b+c)+ Func'[a,b,c]
Func'[a,b,c]包含了类似的项,大概只是a or c等等的次方不同,
我固定其他变量,对其中一个,例如a画图,他出现了很多锯齿状
(参数空间代a>>b, for example: a=10^5,b=10^-5),
就算改变plotpoint也没救-.-
http://imgur.com/z3K0O8Y
另外,我有一个条件是:
a >> b 所以我就带入
Func[a,b,c]/.{a-b -> a}
锯齿状的问题就解决了,但方程式的值(a>>b)和原本具有锯齿状的方程不同,
甚至连方程式画出来没有锯齿状的点也不同(一样保持a>>b,但其他是不同的参数空间),
我就疑惑,到底哪一个才可信..
想请问有人遇到过类似的问题吗?和处理他的方法
感谢,不知道这样问得清不清楚。
作者: sunev (Veritas)   2015-11-12 00:36:00
WorkingPrecision 调高一点?
楼主: a29788685 (嘉峰)   2015-11-13 11:44:00
感谢回答。似乎是没什么用,当分母趋近0他数值会乱跑所以我跑去更早更早之前推导这方程式的地方直接下a>>b的条件,结果是好很多-.-我发现当分母很小,例如我设定b=0.1 , c=0.11的情况使用Expand or Simplify等等的,方程式出现的数值解会不同,这点可能要很小心
作者: AmibaGelos (Amiba Gelos)   2015-11-13 17:27:00
你确定第一项有问题吗? 看起来就只是a/(a+1)的形式啊
楼主: a29788685 (嘉峰)   2015-11-13 22:47:00
他是有非常非常多项,形式都和第一项一样
作者: pig030 (FEBUR.PHEIX)   2015-11-22 17:37:00
转换成另一个空间看看
作者: AmibaGelos (Amiba Gelos)   2015-11-22 22:05:00
假设a>>b不够,c也会移动pole位置,全叠起来就会变那样

Links booklink

Contact Us: admin [ a t ] ucptt.com