PTT數位生活區 / Mathematica

[問題] NDSolve 數值上的問題

看板Mathematica作者 (Kuma)時間4年前 (2019/08/04 21:24), 編輯推噓3(3推 0噓 23→)
留言26則, 4人參與, 最新討論串1/1
版上的大神們好， 小弟最近利用MMA跑聯立微分方程式， 但會遇到數值精確度跑掉的問題（如圖所示） https://i.imgur.com/lsCGEyi.jpg
麻煩大神們撥冗幫助小弟，謝謝。 ----- Sent from JPTT on my iPhone -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 49.216.19.47 (臺灣) ※ 文章網址: https://www.ptt.cc/bbs/Mathematica/M.1564925058.A.2AA.html
推
08/05 09:30, 4年前 , 1F
數值方法會跑到這麼大的數出現可能有些發散性問題在
08/05 09:30, 1F
→
08/05 09:30, 4年前 , 2F
這比較不是精確度問題: 注意到錯誤訊息中的那些數
08/05 09:30, 2F
→
08/05 09:31, 4年前 , 3F
是 10 的 307 次方這種超大數
08/05 09:31, 3F
→
08/05 09:32, 4年前 , 4F
之所以問題不在倒數結果的原因是, 看起來倒數是方程一部份
08/05 09:32, 4F
→
08/05 09:33, 4年前 , 5F
(MMA 裡除法跟倒數是同樣的東西) 那這就代表這個超大數
08/05 09:33, 5F
→
08/05 09:33, 4年前 , 6F
確實是數值計算過程中的中間值, 那就要考慮為何會有這種值
08/05 09:33, 6F
→
08/05 10:33, 4年前 , 7F
您好，這個值確實是我要的，但希望MMA可以繼續算下去，請
08/05 10:33, 7F
→
08/05 10:33, 4年前 , 8F
問有什麼辦法嗎？
08/05 10:33, 8F
→
08/05 11:05, 4年前 , 9F
調working precision開big# 記得常數也要是big# 如有
08/05 11:05, 9F
→
08/05 11:05, 4年前 , 10F
理數 或改用exp寫該# 除法改減法 稍不精確但至少能跑
08/05 11:05, 10F
推
08/13 06:08, 4年前 , 11F
我有點好奇是什麼樣的數值微分方程會需要(可能暫時性地)
08/13 06:08, 11F
→
08/13 06:08, 4年前 , 12F
出現這種數量級的數...
08/13 06:08, 12F
→
08/13 06:10, 4年前 , 13F
就我對一些數學/物理/etc 方面的理解來看, 10^307 這數量級
08/13 06:10, 13F
→
08/13 06:10, 4年前 , 14F
大到不太可能會有一些應用上的意義
08/13 06:10, 14F
→
08/13 06:11, 4年前 , 15F
而數值模擬微分方程是一個很應用方向的東西
08/13 06:11, 15F
→
08/13 06:12, 4年前 , 16F
因此就算是中間值都會有某種程度上的應用意義在
08/13 06:12, 16F
→
08/13 06:13, 4年前 , 17F
這兩件事情是有點矛盾的, 所以才會認為問題是方程有發散性
08/13 06:13, 17F
→
08/13 06:14, 4年前 , 18F
而若發散性是可能狀況之一的話, 發散到多大就不重要了
08/13 06:14, 18F
→
08/16 00:47, 4年前 , 19F
抱歉，忘了回。我處理的系統是需要算高階的correlation f
08/16 00:47, 19F
→
08/16 00:47, 4年前 , 20F
unction, 所以當correlation function 在很高階的時候，
08/16 00:47, 20F
→
08/16 00:48, 4年前 , 21F
就會有這種數量級，或許我可能改用其他的方法來算我感興
08/16 00:48, 21F
→
08/16 00:48, 4年前 , 22F
趣的物理量。 如果不大懂我在說什麼，可以站內信，我們可
08/16 00:48, 22F
→
08/16 00:48, 4年前 , 23F
以一起討論。
08/16 00:48, 23F
→
08/16 10:15, 4年前 , 24F
1T site 10階delta也就~10^120 除非在算amr而且系統
08/16 10:15, 24F
→
08/16 10:15, 4年前 , 25F
類似NS不平滑 個人覺得10^307都可以當作delta處理
08/16 10:15, 25F
推
09/07 01:53, , 26F
有沒有可能先行處理一處方程式讓他們更簡化再來跑模擬
09/07 01:53, 26F
‣ 更多 domolcy 的文章
文章代碼(AID): #1THjo2Ag (Mathematica)
Mathematica 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 domolcy 的文章
文章代碼(AID): #1THjo2Ag (Mathematica)