版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zyj8119

木虫 (著名写手)

模拟EPI: 10
应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【转帖】Monte Carlo 方法计算数值积分已有1人参与

一位朋友说要贴出Monte Carlo计算积分的源程序，我就随便做一个简单的吧，复杂的程序完全可以从这个来演化，我想Monte Carlo积分的最大优势就在于高维积分，以及不规则区域，可以节约很多计算机时。

下面只是演示一个2重积分，可以扩展到20维的只要添加相应的loop项。

被积函数： exp(sqrt(x^2+y^2));
x 上下限：x^2 < sin(y)
y 上下限： y^2CODE:

% MONTE CARLO INT

% by caoer

clear all

N =100000;

x = 2*rand(N,1)-1;

y = rand(N,1);

f = 0;

fsq = 0;

n = 0;

for i=1:N

        if x(i)^2 < sin(y(i)) & y(i)^2
        n=n+1;

        f=f+exp(sqrt(x(i)^2+y(i)^2)); 

        x_plot(n) = x(i);

        y_plot(n) = y(i);

    end

end

f = f/N;

p=n/N;

a = 2*1;

I1 = f*a;  %

I2 = a*p;  %interesting area

I = I1/I2

plot(x_plot,y_plot,'o')

蒙特卡洛法用于求积分时，与积分重数无关，这点非常重要。虽然四维以下的积分用蒙特卡洛法效率可能不如传统的一些数值积分方法，但是维数高的时候，蒙特卡洛法比传统方法要有效的多，而且实现起来也非常容易。可以说，计算高维积分是蒙特卡洛方法最成功和典型的应用。

基本的蒙特卡洛法具有计算不可重复性的缺点。这里共享采用等序列分布的蒙特卡洛法，具有计算可重复性，误差阶比采用基本蒙特卡洛法好的优点。就基本的蒙特卡洛法求积分来说，不管积分重数多少，基本上是计算规模增加100倍，精度提高10倍。

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

Monte carlo模拟已经有5人回复
【书籍】《应用随机过程教程及在算法和智能计算中的随机模型》龚光鲁、钱敏平著已经有46人回复
【转帖】一道面试题的monte carlo实现已经有8人回复

好好学习，天天向上。

1楼 2010-09-17 09:46:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

luomingqi

木虫 (正式写手)

应助: 47 (小学生)
金币: 4356.8
散金: 20
红花: 1
帖子: 870
在线: 66.5小时
虫号: 896596
注册: 2009-11-07
性别: GG
专业: 控制理论与方法

支持一下的