作者lucow (lucow)
看板MATLAB
标题[讨论] 跳跃几何布朗运动模拟
时间Sat May 24 12:40:46 2014
我想模拟几何布朗运动with跳跃
写完後发现跟论文给的图不同，有怀疑过参数问题
不知道有没有人有code可以帮我
或者教学档、note之类的
function Sp =
StockPriceProcess_JumpLogNormalGeometrixBrownianMotion(T,N,S,r,vol,lambda,a,b)
step=T/N;
t=[0:step:T];
W=[0,cumsum(random('Normal',0,sqrt(step),1,N),2)];
J=random('Poisson',(lambda*step),1,N)
loss=ones(1,N)
for i=1:N
if J(1,i)==0
loss(1,i)=loss(1,i);
else
loss(1,i)=exp(sum(random('Normal',a,b,[J(1,i),1])))
end
end
loss=cumprod(loss,2)
Sp=S*exp((r-0.5*vol^2)*t+vol*W);
for i=1:N
if J==0
Sp(1,i+1)=Sp(1,i+1);
else
Sp(1,i+1)=Sp(1,i+1)*loss(1,i);
end
end
plot(t,Sp)
end
这是我的程式码，股价过程是Merton解完随机方程後的形式
最後画图那边我也有问题
股价跳跃上去或下去，我想让有跳跃的地方没有连线，但是我做不出来...
模拟参数 S=20 r=0.06 vol=0.2 a=0.23 b=0.2 T=1 N=1000
lambda=10(跑出来的股价路径很奇怪，後来调成1好像比较正常)
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 27.105.20.232
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1400906448.A.FA1.html
1F：→ mp19990920:你的图？你的程式？你的参考论文？你的模拟参数？ 05/24 16:00
※ 编辑: lucow (219.85.123.44), 05/24/2014 23:44:56
※ 编辑: lucow (219.85.123.44), 05/24/2014 23:56:19
2F：推 jack10313:lambda? 05/25 00:03
※ 编辑: lucow (27.105.11.144), 05/25/2014 10:08:42
3F：推 DIDIMIN: 05/26 13:15