MC积分
MC积分
MC积分。蒙特卡洛方法最重要的应用在于对高维定积分的计算。采用蒙卡模拟的方法来求解这种确定性的数学物理问题时,必须先选择一个合适的概率模型,利用此概率模型试验所得的随机事件的统计结果等价于待求问题的解。。设想一个物理系统由相互作用的多个粒子组成,如凝聚态物质中的原子、原子中的电子,而每个粒子都可以有几个自由度,要描述这个系统就要涉及高维的积分。●如对m个原子组成的气体,经典配分函数是:=dr...drmeap-β2>U(rii)i<i这里的β=(kBT)-1代表温度,U(rii)是两体相互作用势。除非m相当小时,这个式子要用其他任意一种数值积分计算方法都是不可能算出的
MC积分
。假设我们要计算地图上某个区域的面积,可以均匀的在包围此区域的一投点法个方框里撒上小石子,待计算面积与正方形面积之比即为该区域中的石子数目与总数之比。显然,石子越小,撒的数自越多越均匀,则求得的面积越准确。这就是用MonteCarlo方法计算定积分的原理。因为f(x)的定积分值即为曲线下的面积值Sf(r)dr = So(n/N)这里的S.为方形区域的面积,N是总点数,n是郑入f(x)下面积区域中的点数,该方法从原理上与舍选抽样法是一致的。例如,用MonteCarlo方法求元值。产生一对在一1,1区间中均匀分布的随机数作为点的坐标x,y)值,判断条件2+y2<1是否成立,成立则计数n值,当总点数N足够大时,有:Tlim22=N8
投点法
平均值法平均值法是建立在大数定理的基础之上:。如果函数f(x在[a.bl区间,以均匀的概率分布密度随机地取N个数xi,对每个x;计算出函数值f(x)。大数法则告诉我们这些函数值之和除以N所得的值将收敛于函数的期望值,即:N1<f)=(a)f(a)dc(b-ai=1Var(<f>)●故有:Var(f)=v2Var(f)N(b-aN2Nf(ci)f()drNi=1。在简单抽样中,我们由均匀分布中选取随机数,并不考虑被积函数的具体情况。因此,被积函数的极大处和极小处有相同的抽样权重。。但是对积分贡献较大的更多在函数值较大处,故直观上就可以看出,非权重的简单抽样效率较低。即为了获得较高计算结果的精度,需要大量的抽样
平均值法 𝑉𝑎𝑟(< 𝑓 >) = 𝑉 2 𝑁2 𝑉𝑎𝑟(𝑓) = 𝑁 1 𝑉 2 𝑉𝑎𝑟(𝑓) 𝑁
中心极限定理与大数法则一样,概率论中的中心极限定理同样是MonteCarlo方法应用于统计计算的基础。它是对积分的蒙特卡洛估计值收敛于该积分的正确结果的收敛程度的研究中心极限定理告诉我们:在有足够大,但又有限的抽样数N的情况下蒙特卡洛估计值的误差分布。即:→ d(β)af/vA。其中μ=一aJf(α)dc,f是函数f(x)的标准偏差,Φ(β)是Gauss正态分布,那么对于积分值得标准偏差。s,就有:ag=f)-μlα,从上式可以看到蒙特卡洛积分的误差和与N有关。在方差固定时,增加模拟次数可以有效地减小误差
中心极限定理