下面给出几个常用的α与的数值: ? 关于蒙特卡罗方法的误差需说明两点:第一,蒙特卡罗方法的误差为概率误差,这与其他数值计算方法是有区别的。第二,误差中的均方差σ是未知的,必须使用其估计值 来代替,在计算所求量的同时,可计算出 。 α 0.5 0.05 0.003 0.6745 1.96 3 减小方差的各种技巧 显然,当给定置信度α后,误差ε由σ和N决定。要减小ε,或者是增大N,或者是减小方差σ2。在σ固定的情况下,要把精度提高一个数量级,试验次数N需增加两个数量级。因此,单纯增大N不是一个有效的办法。 另一方面,如能减小估计的均方差σ,比如降低一半,那误差就减小一半,这相当于N增大四倍的效果。因此降低方差的各种技巧,引起了人们的普遍注意。后面课程将会介绍一些降低方差的技巧。 效率 一般来说,降低方差的技巧,往往会使观察一个子样的时间增加。在固定时间内,使观察的样本数减少。所以,一种方法的优劣,需要由方差和观察一个子样的费用(使用计算机的时间)两者来衡量。这就 是蒙特卡罗方法中效率的概念。它定义为 ,其中c 是观察一个子样的平均费用。显然 越小,方法越有效。 蒙特卡罗方法的特点 优点 能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。 受几何条件限制小。 收敛速度与问题的维数无关。 具有同时计算多个方案与多个未知量的能力。 误差容易确定。 程序结构简单,易于实现。 缺点 收敛速度慢。 误差具有概率性。 在粒子输运问题中,计算结果与系统大小有关。 能够比较逼真地描述具有随机性质的事物的特点及物理实验过程 从这个意义上讲,蒙特卡罗方法可以部分代替物理实验,甚至可以得到物理实验难以得到的结果。用蒙特卡罗方法解决实际问题,可以直接从实际问题本身出发,而不从方程或数学表达式出发。它有直观、形象的特点。 受几何条件限制小 在计算s维空间中的任一区域Ds上的积分 时,无论区域Ds的形状多么特殊,只要能给出描述Ds的几何特征的条件,就可以从Ds中均匀产生N个点 ,得到积分的近似值。 其中Ds为区域Ds的体积。这是数值方法难以作到的。 另外,在具有随机性质的问题中,如考虑的系统形状很复杂,难以用一般数值方法求解,而使用蒙特卡罗方法,不会有原则上的困难。 收敛速度与问题的维数无关 由误差定义可知,在给定置信水平情况下,蒙特卡罗方法的收敛速度为    ,与问题本身的维数无关。维数的变化,只引起抽样时间及估计量计算时间的变化,不影响误差。也就是说,使用蒙特卡罗方法时,抽取的子样总数N与维数s无关。维数的增加,除了增加相应的计算量外,不影响问题的误差。这一特点,决定了蒙特卡罗方法对多维问题的适应性。而一般数值方法,比如计算定积分时,计算时间随维数的幂次方而增加,而且,由于分点数与维数的幂次方成正比,需占用相当数量的计算机内存,这些都是一般数值方法计算高维积分时难以克服的问题。 具有同时计算多个方案与多个未知量的能力 对于那些需要计算多个方案的问题,使用蒙特卡罗方法有时不需要像常规方法那样逐个计算,而可以同时计算所有的方案,其全部计算量几乎与计算一个方案的计算量相当。例如,对于屏蔽层为均匀介质的平板几何,要计算若干种厚度的穿透概率时,只需计算最厚的一种情况,其他厚度的穿透概率在计算最厚一种情况时稍加处理便可同时得到。 另外,使用蒙特卡罗方法还可以同时得到若干个所求量。例如,在模拟粒子过程中,可以同时得到不同区域的通量、能谱、角分布等,而不像常规方法那样,需要逐一计算所求量。 误差容易确定 对于一般计算方法,要给出计算结果与真值的误差并不是一件容易的事情,而蒙特卡罗方法则不然。根据蒙特卡罗方法的误差公式,可以在计算所求量的同时计算出误差。对干很复杂的蒙特卡罗方法计算问题,也是容易确定的。 一般计算方法常存在着有效位数损失问题,而要解决这一问题有时相当困难,蒙特卡罗方法则不存在这一问题。 程序结构简单,易于实现 在计算机上进行蒙特卡罗方法计算时,程序结构简单,分块性强,易于实现。 收敛速度慢 如前所述,蒙特卡罗方法的收敛速度为 ,一般不容易得到精确度较高的近似结果。对于

matlab蒙特卡罗方法求体积_蒙特卡罗方法详细讲解与MATLAB实现.ppt相关推荐

  1. 灰度重采样的方法分为_重采样方法比较~

    展开全部 最邻近法(Nearest Neighbor):最邻近法直接将与某像元位置最邻近的像元值作为该像32313133353236313431303231363533e58685e5aeb93133 ...

  2. matlab相关性分析频谱_信号与系统:用matlab分析wav音频的频谱

    信号与系统:用matlab分析wav音频的频谱 2018-11-25 [天水浪客一门课程的平时作业,matlab确实强大,用它做这个有点杀鸡用牛刀的意味...软件太大了,直接靠安装了此软件的同学搞定, ...

  3. matlab蒙特卡罗方法求体积_实验二-蒙特卡罗方法计算三维体积

    班级: 学号: 姓名: 实验时间: 2014 年 月 日 实验 项目 实验二 蒙特卡罗方法计算三维体积 所属 课程 数学实验 实 验 目 的 了解蒙特卡罗方法的原理,掌握随机数使用技术. 实 验 内 ...

  4. python蒙特卡洛方法圆周率_蒙特卡罗方法 python 实现

    蒙特卡罗(Monte Carlo)方法的精髓:用统计结果去计算频率,从而得到真实值的近似值. 一.求圆周率的近似值,采用 投点法 import numpy as np import matplotli ...

  5. 抽象工厂和工厂方法示例_工厂方法设计模式示例

    抽象工厂和工厂方法示例 本文是我们名为" Java设计模式 "的学院课程的一部分. 在本课程中,您将深入研究大量的设计模式,并了解如何在Java中实现和利用它们. 您将了解模式如此 ...

  6. 工厂方法模式_工厂方法模式

    工厂方法模式是简单工厂模式的升级版,简单工厂模式不符合设计模式的原则(即:单一职责,开闭原则) 优点: 职责明确,扩展方便 缺点:需要创建多个工厂 实现步骤: 1.将工厂通用方法抽取接口 (例如:IF ...

  7. 构造函数也可以进行方法重载_防蛀方法不止涂氟、窝沟封闭,还有一种方法在家就可以进行。...

    我们都对牙膏很熟悉,而且天天用它清洁牙齿.如果不这样做,就会有一层牙菌斑长时间附着在牙齿表面,从而损害牙釉质和牙龈. 每天清理牙菌斑是保持良好的口腔卫生,维持牙齿健康的首要因素.然而,除了使用牙膏,现 ...

  8. python 魔法方法常用_魔法方法使用指南(python)

    1 什么是魔法方法? 它们在面向对象的Python的处处皆是 它们是一些可以让你对类添加"魔法"的特殊方法 它们经常是两个下划线包围来命名的(比如 __init__,__lt__) ...

  9. split()方法的用法,超详细讲解

    @java 认识split()方法 在java中常常会有对字符串进行的操作, 而有时候我们想把一个字符串拆分成为多个字符串, 这个时候我们就用到了**split()**方法. 本文就是解释split( ...

最新文章

  1. 分割catalina.out 每天生成一个文件
  2. Linux 帧缓存 数据,linux--帧缓存技术
  3. 平面设计师如何训练自己创意思维
  4. 利用Flask来构建项目的大概步骤
  5. python傅里叶函数图像_python实现傅里叶级数展开的实现
  6. 1-4:学习shell之操作文件与目录
  7. python3 try except or_Python基础10:try except异常处理详解
  8. Backbone 学习笔记
  9. 小米MIX4最新渲染图曝光;苹果首款OLED iPad将于2023年问世;格力将引进鸿蒙操作系统|极客头条...
  10. 汉字Unicode编码表
  11. 基于python的网络聊天室论文_Python基于Socket实现简单聊天室
  12. xp mysql 5.6安装版下载地址_xp系统安装版下载,xp系统安装版下载地址
  13. Pythonic写法
  14. bilibili有的视频没声音解决办法
  15. ug电子产品计算机建模,UG10网格曲面咖啡壶
  16. 四阶魔方邻角互换公式
  17. EOJ - 我决不会TLE (一个智障的题目)
  18. 关于grep,sed,awk的几个用法
  19. java导出word复选框变成圆圈的解决办法
  20. 2016年8月19日 星期五 --出埃及记 Exodus 16:20

热门文章

  1. 高通骁龙410e/APQ8016E嵌入式物联网模块芯片
  2. SAP EWM TCODE list
  3. 如何反编译 cocos creator 生成 的jsc文件/反编译jsc文件(一)
  4. C# linq的学习及使用
  5. 违反唯一约束条件 ORA-00001
  6. 智能优化算法(源码)-食肉植物算法(Carnivorous Plant Algorithm ,CPA)
  7. 4g运行内存手机还能用多久_手机4G和6G运行内存有多大区别?看完秒懂
  8. html 实时统计字数,记一次前端 input、textarea输入框实时 统计字数(真实字数)...
  9. 一文详解工业视觉中的光源及应用
  10. PHP 蒙太奇马赛克拼图,AndreaMosaic制作一幅马赛克拼图