[問題] 雙重積分速度問題已回收
請教大家...現在的目的是要解下面這個方程式
T 10 10 ∞ 0
Σ Σ Σ ∫ ∫ normpdf(Rt(t),U+D,3).*gampdf(U,m,5).*gampdf(-D,n,2) dDdU
t m n 0 -∞
其中t,m,n都是從1開始,T約為2700
三個加總是用for迴圈處理,雙重積分是用dblquad處理
程式碼如下
p3=0
for t=1:T
for m=1:10
for n=1:10
p3=p3+dblquad(@(D,U) normpdf(Rt(t),U+D,3).*...
gampdf(U,m,5).*gampdf(-D,n,2) ,-1000,0,0,1000);
end
end
end
但是計算的時間大約要40分鐘左右
想請問大家有沒有讓處理速度比較快的指令或是方法
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.119.143.110
※ 編輯: qpoiuy 來自: 140.119.143.110 (09/26 12:23)
推
09/26 15:06, , 1F
09/26 15:06, 1F