[問題] 請問應該算有關迴圈的問題吧?已回收

看板MATLAB作者 (well)時間16年前 (2010/02/19 23:16), 編輯推噓0(004)
留言4則, 2人參與, 最新討論串1/1
d(rt+1)=a(b-rt)dt+σ(r^0.5)dzt (rt dzt 那個t代表的是時間點t) rt+1=rt+drt ................................以下為程式碼................................ a=0.05; b=0.03; r0=0.02; σ=0.01; dt=1; e1=randn(1000,100) r=[r0.*ones(1000,1) r0+a*(b-r0)*dt+σ.*sqrt(r0).*sqrt(dt).*e1] end .............................................................................. 我想請問倒數第二行的程式碼 我覺得寫錯了(應該也錯了吧?) 因為此行程式碼所寫出的矩陣第一行為r0 第二行是算1000個可能的r1(算r1需用到r0) 第三行是算1000個可能的r2(算r2需用到r1) 以此類推 我想請問這要怎麼寫 不知道說明的夠不夠清楚 謝謝大家 下台一鞠躬 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.223.38.91

02/19 23:31, , 1F
for-end r0模擬出1000個dr算出r1再模擬1000個dr算出r2
02/19 23:31, 1F

02/20 16:18, , 2F
可是我要算到r120
02/20 16:18, 2F

02/20 16:19, , 3F
這樣會不會要打很久啊?
02/20 16:19, 3F

02/20 18:00, , 4F
用迴圈
02/20 18:00, 4F
文章代碼(AID): #1BVgjCv2 (MATLAB)
文章代碼(AID): #1BVgjCv2 (MATLAB)