拟蒙特卡罗方法

✍ dations ◷ 2025-04-04 11:20:35 #蒙地卡罗方法,拟随机数

数值分析中,拟蒙特卡罗方法(Quasi-Monte Carlo method)是使用低差异列(一种确定生成的超均匀分布列,也称为拟随机列、次随机列)来进行数值积分和研究其它一些数值问题的方法。而普通的蒙特卡罗方法或蒙地卡罗积分方法使用的是伪随机数。MATLAB中提供了生成如哈尔顿列、索博尔列等超均匀分布列的函数。

拟蒙特卡罗方法和蒙特卡罗方法的具体内容相似,要解决的问题都是通过测量某个可测函数 在某些点上的取值,而在数值上求它的积分的近似值。例如要求在单位体积 s {\displaystyle ^{s}} 1, ..., ,那么:

其中的都是s维向量。拟蒙特卡罗方法和普通蒙特卡罗方法的区别在于的具体选取方式。蒙特卡罗方法用的是伪随机列,而拟蒙特卡罗方法用到的是哈尔顿列、索博尔列等低差异列。使用低差异列的优点是收敛速率较快。拟蒙特卡罗方法可以达到O(1/N)的收敛速率,而普通蒙特卡罗方法的收敛速率则是 O(N-0.5)。

近年来,拟蒙特卡罗方法在金融数学和计算机数学领域里得到了越来越多的应用,因为其中常常会需要计算高维积分的数值近似。蒙特卡罗方法和拟蒙特卡罗方法可以快捷简单地得到较好的结果。

拟蒙特卡罗方法的近似误差可以用取点1, ..., 的差异度作为上限。具体来说,Koksma-Hlawka不等式表明,误差项

限制,其中V(f)为函数的Hardy-Krause变差,DN是(x1,...,xN)的差异度,定义为

其中Q是任何s中边界与坐标轴平行的方形“块”。 | ϵ | V ( f ) D N {\displaystyle |\epsilon |\leq V(f)D_{N}} 表明拟蒙特卡罗方法的近似误差大约是 O ( 1 N ) {\displaystyle O({\frac {1}{N}})} 的量级,于此相对的是普通蒙特卡罗方法的近似误差为 O ( 1 N ) {\displaystyle O({\frac {1}{\sqrt {N}}})} 量级。注意这里的不等式给出的是误差上限,事实上拟蒙特卡罗方法的收敛速率要比其上限所示的速率快得多。因此,一般来说拟蒙特卡罗方法比起普通的蒙特卡罗方法来说大大加快了收敛的速率。

相关

  • Ir4f14 5d7 6s22, 8, 18, 32, 15, 2蒸气压第一:880 kJ·mol−1 第二:1600 kJ·mol主条目:铱的同位素铱是原子序77号的元素,其元素符号为Ir,是一种十分刚脆,银白色的铂系过渡金属
  • HER2/neu1MFG, 1MFL, 1MW4, 1N8Z, 1QR1, 1S78, 2A91, 2JWA, 2KS1, 2L4K, 3BE1, 3H3B, 3MZW, 3N85, 3PP0, 3RCD· protein tyrosine kinase activity · transmembrane receptor pro
  • Psub2/subOsub3/sub三氧化二磷(化学式:P4O6)是一种由氧和磷组成的化合物,为无色蜡状有大蒜气味的极毒晶体,虽然它的正确名称应该是六氧化四磷,但以前一直以为分子结构是P2O3,因此三氧化二磷的名称一直
  • 武夷黑猪武夷黑猪,又称闽北黑猪、赣东黑猪,产于武夷山脉两侧的家猪品种。目前已处于极度濒危状况。因武夷山分属福建省和江西省,武夷黑猪在江西境内称赣东黑猪,福建境内称闽北黑猪。主产
  • 并发计算并发计算(英语:Concurrent computing,或译为并发处理、共时计算),是一种程序计算的形式,在系统中,至少有两个以上的计算在同时运作,计算结果可能同时发生。用来实现并发系统(Concurre
  • 泡利方程泡利方程或称薛定谔-泡利方程,为描述带有自旋1/2的粒子在与电磁场相互作用下的修正方程(自旋1/2粒子例如电子)。在此之前,用以描述粒子行为的薛定谔方程则未考虑到粒子自旋的性
  • 按人口排列的世界城市列表
  • 东泇河东泇河,又称柴沟河,位于山东省南部和江苏省北部,为邳苍分洪道右岸支流。发源于山东省费县南部,源头分东西两支:冬支源于马庄镇梧桐峪村刘家岗东麓;西支发源于新庄镇官庄东北的箕山
  • 瓦莱丽·霍布森瓦莱丽·霍布森(1917年4月4日-1998年11月13日)是一名英国女演员。她于1935年在《科学怪人的新娘》中扮演弗兰肯斯坦男爵夫人一角。同年她在《伦敦狼人》中与亨利·赫尔演对手戏
  • 弗兰克·鲍沙其弗兰克·鲍沙其(英语:Frank Borzage,/bɔːrˈzeɪɡiː/,1894年4月23日-1962年6月19日)是一位美国导演。弗兰克·鲍沙其在整个1920年代是一个相当成功的导演,达到在20世纪默片时代