[問題] 請問模擬股價走勢加上jump process已回收

看板MATLAB作者 (學海無涯回頭是岸)時間16年前 (2009/12/17 20:41), 編輯推噓0(007)
留言7則, 2人參與, 最新討論串1/1
請問關於模擬股價走勢加上jump process如下 dp(t)= mu(t)*dt + sigma(t)*dW(t) + k(t)*dq(t) 前面沒什麼問題,我用的code如下: Spaths=zeros(Nrepl,1+Nsteps); Spaths(:,1)=S0; dt=T/Nsteps; nudt=(mu-0.5*sigma^2)*dt; sidt=sigma*sqrt(dt); x=randn(Nsteps,1); for i=1:Nrepl for j=1:Nsteps Spaths(i,j+1)=Spaths(i,j)*exp(nudt+sidt*randn); end end 但是最後一個部分的jump process,我想拆成兩個部分, k(t)管jump size,我打算抽常態亂數, q(t)的部分我想用poisson process,讓它有jump就跳1,沒有就跳0,隨機跳, 想請問版友q(t)的poisson process要怎麼寫進去?? 我查了好多地方都找不到,還請大家幫幫忙,謝謝~!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.231.83.233

12/17 21:22, , 1F
股價跳躍部份 你總會有個model吧
12/17 21:22, 1F

12/17 21:35, , 2F
jump的部分原始模型設計就是k(t)*dq(t)
12/17 21:35, 2F

12/17 21:36, , 3F
所以我想要寫出一個只跳0跟1又服從poisson的序列
12/17 21:36, 3F

12/17 21:36, , 4F
不知道這樣有沒有回到到您的問題?
12/17 21:36, 4F

12/18 10:25, , 5F
手上只有以Merton1973提出的jump process model程式碼
12/18 10:25, 5F

12/18 10:26, , 6F
這篇paper跟你的model對股價過程的設定不一樣
12/18 10:26, 6F

12/18 10:27, , 7F
不過相同的是possion跟指數分配的模擬
12/18 10:27, 7F
文章代碼(AID): #1BAYRYgX (MATLAB)