http://blog.sina.com.tw/mny888/
如果能真正創造額外收入 那該有多好
提供您一個可全職也可兼職的好機會
0元創業 在家工作 時間彈性
只有你要, 我們才能伸出援手!!
你我當中，有誰不想增加【額外收入】？ 若大家能【共好】，我們願意免費協助你真正
創造收入!
請上官網了解http://blog.sina.com.tw/mny888/
請上GOOGLE搜尋 樹馬矢星聯盟分公司
如有打擾 萬分抱歉 對不起 感謝您
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 111.184.22.29
老師您好!!
我是前幾個月向您請教過 Mathematica 問題的學生
我們實驗室有另一位學弟也有 Mathematica 的問題想請教您
他作的一些計算模擬同時用到 MATLAB 和 Mathematica
他最近在程式撰寫上遇到瓶頸 想請教您
請問我可以將您的聯絡資料給他嗎? ^^
感謝您 辛苦了 !! 祝一切順利 !!
pxpypz 敬上
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 42.72.114.214
※ 引述《narwhalgary (獨角鯨葛瑞)》之銘言：
: For[i=1,i<201,i++,
: cn[i]= Maximize[f(i,r),r][[2,1,2]]; ← 希望能加速的程式碼 不知如何Compile
: ......
: ]
: 煩請大大們解答 萬分感激
Maximize這些Mathematica的內建函數基本上都已經Compile過了∼
你要的可能是快速的將這200的最佳化問題計算完，不過既然都要計算，
把[[2,1,2]]的部份先保留起來
以下用平行運算
f[i_, r_] := a*i*r^2 - 2 b r + c
cn[i_Integer] := Minimize[f[i, r], r]
DistributeDefinitions[cn, f, Range[200]];
myans=ParallelMap[cn, Range[200]]
myans[[All, 2, 1, 2]]
另外，Mathematica最佳化的指令可以看查一下這個help
ConstrainedOptimizationIntroduction
如果問題沒有要求解析解，NMaximize 或 NMinimize倒是不錯的選擇！
--
養花種魚數月亮賞星星
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.130.134
