作者:
oddyssey (iuygiuti)
2017-01-18 01:13:42抱歉 因為把問題貼在專版都沒人幫忙回 所以只好來這碰碰運氣了
想請問各位大大及做數理精算的神人 或是為準備期末考而睡不著覺的高手高手高高手
我正在測試first-passage-time 的機率密度 random time 是由標準幾何布朗運動所定義
假設密度函數是 f(t; mu, sigma) mu 是 drift, sigma 是 diffusion 兩者都為常數
按照機率法則 若針對時間 t 的全域(從零到無窮大)做積分 則數值積分應該值會接近一
我在 matlab 下測試這個機率特性沒有問題 但是如果用mathematica的內建 NIntegrate
則結果不會 robust 特別是當 sigma 小於 0.2 時 積分值往往小於一很多 (大約 0.5)
但是只要 sigma 給得值夠大就一切正常 可是類似情況在 matlab 下積分值都是接近一
matlab 的測試結果不受參數值選擇而影響 所以 這是否為 mathematica 的 bug??
喔 對了 小弟用的是 6.0 版本 有沒有其他大大曾經也發現過這種怪異現象 懇請解惑
作者:
gino0717 (gino0717)
2017-01-18 01:14:00我都用matlab
作者: yellowlin (エライザ><) 2017-01-18 01:15:00
只用過matlab啦QQ
作者:
harry901 (harry901)
2017-01-18 01:15:00數值積分如果是很複雜的系統 最好用MATHEMATICA
作者:
ZXCWS (兩分銅幣)
2017-01-18 01:15:00你是物理數學系 還是土木系
作者:
gofee (撞牆的青蛙)
2017-01-18 01:16:00幹你老師不會啦 馬的 救救我的波德圖
作者:
harry901 (harry901)
2017-01-18 01:16:00MATLAB在處理同樣的數值積分 介面是比MATHEMATICA方便
作者: jenkl 2017-01-18 01:16:00
數值積分可作無窮大?
作者:
ZXCWS (兩分銅幣)
2017-01-18 01:16:00跑數值要好運就用Matlab 要效能就用Fortran 誰用MATHEMATIC
作者:
Cybershit (BMI<18der肥宅)
2017-01-18 01:17:00你去 Stack Overflow 問, 比來這問好很多
作者: scatology (scatology) 2017-01-18 01:18:00
Fortran
作者:
harry901 (harry901)
2017-01-18 01:18:00FORTRAN會用的人當然建議用阿 不過現是說兩家數學軟體
沒人在用mathematica啦 那過氣幾十年的東西
作者:
harry901 (harry901)
2017-01-18 01:20:00作者:
oddyssey (iuygiuti)
2017-01-18 01:24:00感謝樓上 ^^
作者:
harry901 (harry901)
2017-01-18 01:24:00再可能存在的問題底下 你可以調整MAXRECURSION通常如果遇到無限大的數值積分 這個MAXRECURSION會出現很多問題 你可以試著調整看看
作者:
harry901 (harry901)
2017-01-18 01:25:00原PO的問題應該是隨機微積分了吧 加油
作者:
oddyssey (iuygiuti)
2017-01-18 01:27:00恩 其實是我做研究遇到的問題 的確和隨機為積分有些關聯
作者:
wen7733 (塔=水怪=阿溫=亡月)
2017-01-18 01:30:00記得在118學Fortran要開還要裝VM跑XP...
6.0多舊了... 而且default設定本來就機會有問題