怎么计算蒙特卡洛模拟结果的期望值？R语言做的，代码放在下面，谢谢各位了！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

S0=2986.6;      #initial price S0
mu=0.019623;    #drift mu
sigma=0.2893;  # volatility sigma
T=1;        #1 year
NSteps=48;
NReps=10000;
SPaths = matrix(NA,NReps,NSteps+1);
SPaths[,1] = S0;
dt = T/NSteps
nudt = (mu-0.5*sigma^2)*dt
sidt = sigma*sqrt(dt)
x=seq(1:(NSteps+1))
colno=c(1,2,3,4,5,6,7,8,9,10)
set.seed(123)
for (i in 1:NReps){
  for (j in 1:NSteps){
    SPaths[i,j+1] = SPaths[i,j]*exp(nudt + sidt*rnorm(1))
  }
  if(i == 1) plot(x,SPaths[i,],"l", xlab="时间t",ylim=c(1500,5000),ylab="上证50指数",col =sample(colno,1))
  if (i != 1) lines(x,SPaths[i,],xlab="时间t",ylim=c(1500,5000),ylab="上证50指数",col =sample(colno,1))

}

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：蒙特卡洛模拟 蒙特卡洛 期望值 蒙特卡 R语言