[問題] 方程式準確性

作者: 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