sas 看涨期权蒙特卡洛模拟方法定价代码编制问题 高手请入

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

此处 t 代表行权时间
r代表连续复利收益率    sig 波动率      m代表模拟m次      n代表每次模拟n步
div代表连续股利收益率
我下面提供的是部分欧式看涨期权的蒙特卡洛模拟的核心代码  如有问题请纠正
但是我本人sas初学 但知道其有强大的分析和作图功能
请高手将我的代码丰富下 能够实现图像形式的股价模拟路径 非常感谢



dt = t/n;                                                                                                                           
nudt = (r-div-0.5*sig^2)*dt;                                                                                                         
sigsdt = sig*sqrt(dt);                                                                                                               
lns = ln(s);                                                                                                                                                                                                                                      
sum_CT = 0;                                                                                                                          
sum_CT2 = 0;                                                                                                                           
do  j = 1 to M;     *for each simulation;                                                                                            
lnst = lns;                                                                                                                           
  do i = 1 to  N ;     *for each time step;                                                                                       
  X=rannor(0);                                                                                                                     
  lnst=lnst+nudt+sigsdt*x;   *evolve the stock price;                                                                              
  end;                                                                                                         

ST = exp(lnst);                                                                                                                     
CT = max{0,ST-K};                                                                                                                    
sum_CT = sum_CT + CT;                                                                                                                 
sum_CT2 = sum_CT2 + CT*CT;                                                                                                            
end;     

call_value = sum_CT/M*exp(-r*t);                                                                                                  
SD = sqrt((sum_CT2 - sum_CT*sum_CT/M)*exp(-2*r*T)/(M-1));                                                                             
SE = sd/sqrt(M);

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：蒙特卡洛模拟 蒙特卡洛 蒙特卡 Simulation 连续复利收益率 蒙特卡洛 收益率

学习了。真的是知识的海洋

个人观点：
你上面想实现M次模拟的结果然后对这模拟值求平均在得到看涨期权价格，但是我觉得你选择的是同一个种子产生随机正态分布数，这样实际上是一个程序，还有你想得到这个价格的路径值，你可以在程序里面加上OUTPUT来产生，但是下面的程序就不是这样了，以下是一个看涨看跌期权的蒙特卡洛模拟：
proc iml;
s0=j(10000,1,100);
x=100;
t=1;
n=250;
sigma=0.1;
rf=0.05;
m=10;
yeild=j(nrow(s0),n,.);
s=j(nrow(s0),n,.);
call=j(nrow(s0),1,.);
put=j(nrow(s0),1,.);
result=j(m,2,.);
do k=1 to m;
do i=1 to nrow(yeild);
do j=1 to ncol(yeild);
yeild[i,j]=(rf-0.5*sigma**2)*t/n+sigma*sqrt(t/n)*rannor(round(k+exp(k),1));
s[i,j]=s0[i]*exp(yeild[i,1:j][,+]);
end;
end;
do i=1 to nrow(s0);
call[i]=max(s[i,n]-x,0);
put[i]=max(x-s[i,n],0);
end;
result[k,1]=call[:,]*exp(-rf*t);
result[k,2]=put[:,]*exp(-rf*t);
end;
print result;
quit; 可以具体研究下！