[問題] 方程式準確性

看板Mathematica作者 (嘉峰)時間8年前 (2015/11/11 19:49), 8年前編輯推噓2(209)
留言11則, 4人參與, 最新討論串1/1
是這樣子的: 我有一個方程式,很醜,大致上像這樣 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,但其他是不同的參數空間), 我就疑惑,到底哪一個才可信.. 想請問有人遇到過類似的問題嗎?和處理他的方法 感謝,不知道這樣問得清不清楚。 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 103.5.101.7 ※ 文章網址: https://www.ptt.cc/bbs/Mathematica/M.1447242591.A.986.html ※ 編輯: a29788685 (103.5.101.7), 11/11/2015 19:55:13

11/12 00:36, , 1F
WorkingPrecision 調高一點?
11/12 00:36, 1F

11/13 11:44, , 2F
感謝回答。似乎是沒什麼用,當分母趨近0他數值會亂跑
11/13 11:44, 2F

11/13 11:44, , 3F
所以我跑去更早更早之前推導這方程式的地方直接下
11/13 11:44, 3F

11/13 11:45, , 4F
a>>b的條件,結果是好很多-.-
11/13 11:45, 4F

11/13 14:45, , 5F
我發現當分母很小,例如我設定b=0.1 , c=0.11的情況
11/13 14:45, 5F

11/13 14:45, , 6F
使用Expand or Simplify等等的,方程式出現的數值解
11/13 14:45, 6F

11/13 14:45, , 7F
會不同,這點可能要很小心
11/13 14:45, 7F

11/13 17:27, , 8F
你確定第一項有問題嗎? 看起來就只是a/(a+1)的形式啊
11/13 17:27, 8F

11/13 22:47, , 9F
他是有非常非常多項,形式都和第一項一樣
11/13 22:47, 9F

11/22 17:37, , 10F
轉換成另一個空間看看
11/22 17:37, 10F

11/22 22:05, , 11F
假設a>>b不夠,c也會移動pole位置,全疊起來就會變那樣
11/22 22:05, 11F
文章代碼(AID): #1MGojVc6 (Mathematica)
文章代碼(AID): #1MGojVc6 (Mathematica)