hfut-liu 发表于 2009-8-22 18:53

怎么样连续运行蒙特卡罗小程序

各位坛友大家好,本人编程能力不好想请大家帮我看一下下面的程序怎么添加使它自动运行。
运用蒙特卡罗方法能够解决简单的数值积分问题,因为有随机数的原因每次生成的随机数都不一样经过几次重复运行就能得到一个近似解。下面的程序是解y=exp(-x2)在(-1,1)之间的定积分,运行十次取平均值课的近似解,但我不知道怎么加程序只要运行一次就出现十个值(因为运行十次好麻烦),请高手帮忙,谢谢了
xi=unifrnd(-1,1,100000,1);
yi=rand(100000,1);
k=0;
y=exp(-xi.^2);
for i=1:100000
    if yi(i)<=y(i)
      k=k+1;
    end
end
I=k/50000

nevsayno 发表于 2009-8-22 21:38

我写过蒙特卡洛方法的程序 ,只是没看懂你什么意思:)


而且,既然是求近似解,用求高斯积分点和牛顿xx积分方法都可以求出很精确的解啊,那个要简单好多.......

如果仍未解决,可以联系我,qq 595538386    呵呵

sogooda 发表于 2009-8-22 21:56

for jjj=1:10
xi=unifrnd(-1,1,100000,1);
yi=rand(100000,1);
k=0;
y=exp(-xi.^2);
for i=1:100000
    if yi(i)<=y(i)
      k=k+1;
    end
end
I=k/50000
end

hfut-liu 发表于 2009-8-23 09:27

恩谢谢大家的帮助,昨晚回去之后想了一下今天来运行一下对了,依然感谢上面两位的回帖。
for j=1:10
xi=unifrnd(-1,1,100000,1);
yi=rand(100000,1);
k=0;
y=exp(-xi.^2);
for i=1:100000
    if yi(i)<=y(i)
      k=k+1;
    end
end
I(j)=k/50000;
end
I
我是这样写的,跟三楼的差不多
页: [1]
查看完整版本: 怎么样连续运行蒙特卡罗小程序