zh.wikipedia.org/wiki/%E8%92%99%E5%9C%B0%E5%8D%A1%E7%BE%85%E6%96%B9%E6%B3%95
3 corrections found
非权重蒙特卡罗积分,也称确定性抽样
这里把普通蒙特卡罗积分说成“确定性抽样”不对。普通(plain/crude)蒙特卡罗积分本质上是随机抽样;“确定性”的对应物通常指拟蒙特卡罗(quasi-Monte Carlo)。
Full reasoning
这句话把普通蒙特卡罗积分和确定性采样/拟蒙特卡罗混为一谈了。
- GNU Scientific Library 的官方文档明确写道:plain Monte Carlo algorithm samples points randomly,并且积分器需要用户提供random number generator。这说明普通蒙特卡罗积分是随机抽样方法,不是“确定性抽样”。
- 相反,文献通常把quasi-Monte Carlo (QMC)称为 Monte Carlo 的deterministic counterpart/version。也就是说,若要说“确定性抽样”,对应的应是拟蒙特卡罗,而不是这里的普通非加权蒙特卡罗积分。
因此,把“非权重蒙特卡罗积分”直接称为“确定性抽样”会误导读者,把随机 Monte Carlo 和确定性的 quasi-Monte Carlo 混为一类。
2 sources
- Monte Carlo Integration - GSL 2.8 documentation
The plain Monte Carlo algorithm samples points randomly from the integration region to estimate the integral and its error... The integration uses a fixed number of function calls ... and obtains random sampling points using the random number generator r.
- Quasi-Monte Carlo integration using digital nets with antithetics
... quasi-Monte Carlo (QMC) integration, a deterministic counterpart of Monte Carlo, of functions defined over the s-dimensional unit cube.
此种方法的正确性是基于概率论的中心极限定理。
把蒙特卡罗积分的“正确性”归因于中心极限定理不准确。其收敛正确性来自大数定律;中心极限定理主要用于描述误差的近似分布和置信区间。
Full reasoning
这句话把两个不同的概率论结论混在了一起。
- 对于基本蒙特卡罗积分,样本平均会收敛到真实积分值这一“正确性/一致性”结论,来自大数定律(Law of Large Numbers)。
- **中心极限定理(Central Limit Theorem)**的作用则是进一步描述:当样本数足够大时,误差经过适当缩放后近似服从正态分布,因此可以据此构造标准误差和置信区间。
卡内基梅隆大学的课程讲义在同一页上把这两点区分得很清楚:先给出 Monte Carlo integration 定理,并明确说“Theorem 3 follows immediately from the Law of Large Numbers”;随后又单独写“By the central limit theorem ...”来推导误差的概率界。
因此,若说“此种方法的正确性是基于中心极限定理”,就把收敛依据(大数定律)和误差近似分布依据(中心极限定理)混淆了。
2 sources
- Monte Carlo Integration
Theorem 3 follows immediately from the Law of Large Numbers... By the central limit theorem and the three sigma rule ...
- Monte Carlo Integration - GSL 2.8 documentation
The routines also provide a statistical estimate of the error on the result... For large N this variance decreases asymptotically as Var(f)/N ... The error estimate itself should decrease as 1/sqrt(N).
使用此种方法所得近似解的统计误差只与m有关
把普通蒙特卡罗积分的误差说成“只与样本数 m 有关”过于绝对。误差的衰减率是 1/√m,但误差大小还取决于被积函数的方差(以及由此带来的维度/问题结构影响)。
Full reasoning
这里把误差的阶数和误差的实际大小混为一谈了。
- 对普通蒙特卡罗积分,常见结论是误差按 O(1/√m) 衰减;这说的是随样本数增长的速率。
- 但误差并不是“只与 m 有关”。GSL 官方文档给出的公式是:估计量的方差对大样本满足 Var(f)/N,也就是说除了样本数 N 之外,还显式依赖于被积函数在积分区域上的方差 Var(f)。
- CMU 的讲义也给出:
[
\mathrm{Var}\left(\frac1N\sum_{n=1}^N f(X_n)\right)=\frac1N\int |f(x)-\mu|^2 p(x),dx.
]
这再次说明误差不仅取决于样本数,还取决于函数本身和采样分布;这些量在不同问题、不同维度下都可能改变。
所以,更准确的说法应当是:普通蒙特卡罗积分的误差率通常是 1/√m 量级,但误差常数还依赖于被积函数/分布,并非只由 m 决定。
2 sources
- Monte Carlo Integration - GSL 2.8 documentation
For large N this variance decreases asymptotically as Var(f)/N, where Var(f) is the true variance of the function over the integration region. The error estimate itself should decrease as 1/sqrt(N).
- Monte Carlo Integration
If further ∫ |f(x)|^2 p(x) dx < ∞, then Var((1/N) Σ f(X_n)) = (1/N) ∫ |f(x)-μ|^2 p(x) dx ... By the central limit theorem ... N = O(1/ε^2).