计算物理蒙特卡罗方法基础.ppt

1,计算机模拟方法,(1)蒙特卡洛方法：随机性模拟方法或统计试验方法，又称蒙特卡洛(Monte Carlo)方法。它是通过不断产生随机数序列来模拟过程。自然界中有的过程本身就是随机的过程，物理现象中如粒子的衰变过程、粒子在介质中的输运过程.等。当然蒙特卡洛方法也可以借助概率模型来解决不直接具有随机性的确定性问题。(2)分子动力学方法：确定性模拟方法。它是通过数值求解一个个的粒子运动方程来模拟整个系统的行为。在统计物理中称为分子动力学（Molecular Dynamics）方法。(3)离散型模拟方法-元胞自动机等,2,What is a Monte Carlo method?,2-1 蒙特卡罗方法的基础知识,the Comte de Buffon needle experiment,AD 1777,3,Stanislaw Ulam(1909-1984),Nicholas Metropolis(1915-1999),蒙特卡洛方法的起源,4,The Name of the Game,Metropolis coined the name“Monte Carlo”,from its gambling casino.,Monte-Carlo,Monaco,5,从蒙特卡洛模拟的应用来看，该类型的应用可以分为三种形式：（1）直接蒙特卡洛模拟。它采用随机数序列来模拟复杂随机过程的效应。（2）蒙特卡洛积分。这是利用随机数序列计算积分的方法。积分维数越高，该方法的积分效率就越高。（3）Metropolis蒙特卡洛模拟 这种模拟是以所谓“马尔科夫”（Markov）链的形式产生系统的分布序列。该方法可以使我们能够研究经典和量子多粒子系统的问题。,6,一 基本思想,直接蒙特卡洛模拟法：对求解问题本身就具有概率和统计性的情况。如：中子在介质中的传播，核衰变过程等，思想是按照实际问题所遵循的概率统计规律，用计算机进行直接的抽样试验，然后计算其统计参数。该方法也就是通常所说的“计算机实验”。,间接蒙特卡洛方法：蒙特卡洛方法也可以人为地构造出一个合适的概率模型，依照该模型进行大量的统计实验，使它的某些统计参量正好是待求问题的解。,7,代表了该运动员的成绩。换言之，为积分的估计值，或近似值。,现假设该运动员进行了N次射击，每次射击的弹着点依次为r1，r2，rN，则N次得分g(r1)，g(r2)，g(rN)的算术平均值,例1 射击问题（打靶游戏）-直接蒙特卡洛方法,假设射击100次，,平均成绩,8,设r表示射击运动员的弹着点到靶心的距离，g(r)表示击中r处相应的得分数（环数），f(r)为该运动员的弹着点的分布密度函数，它反映运动员的射击水平。该运动员的射击成绩为,用概率语言来说，g(r)是随机变量，的数学期望，即,在该例中，用N次试验所得成绩的算术平均值作为数学期望的估计值（积分近似值）。,9,(1)巴夫昂(Buffon)投针实验,实验方案：在平滑桌面上划一组相距为s的平行线，向此桌面随意地投掷长度l的细针，那末从针与平行线相交的概率就可以得到的数值。,例2 圆周率的数值计算-间接蒙特卡洛方法,10,数学统计理论计算：,针的投影长度,确定的，相交的概率,的平均值,假如在N次投针中，有M次和平行线相交。当N充分大时，相交的频数 M/N就近似为细针与平行线相交的概率。,11,经过n次投针后得到值的精度,针与平行线相交的概率,针与平行线相交的次数应满足二项式分布,其期望值为,的方差,的标准误差,的标准误差,相交和不相交,12,这意味着试验所得的值的不确定性的范围如下：对100次投针为，0.2374 对10,000次投针为，0.0237 对1,000,000次投针为，0.0024 可见，增加模拟的次数可以减小误差，但不可消除误差。,的标准误差,13,前人进行了实验，其结果列于下表：,14,(2)投点法实验,实验方案：在平滑桌面上划正方形，同时划一内切圆，向此正方形随意地投点，那末投点落在圆内的概率就可以得到的数值。,任意投点落在圆内的概率,15,的标准误差,的标准误差,的标准误差,的标准误差,投针实验的误差分析,投点实验的误差分析,对100次投针为，0.1642 对10,000次投针为，0.0164 对1,000,000次投针为，0.0016,16,投点法实验程序流程图,17,program mainuse freconstantuse randomnameimplicit noneinteger nmax,minteger i,ncountreal*8 lenr,lens,lenr2real*8 x,y,dxy2open(10,file=Pi.dat)call randomval()lenr=1.0d0lens=2.0d0*lenrlenr2=lenr*lenrm=0ncount=0write(*,*)Input nmax:read(*,*)nmax,do i=1,nmax call randomnum()x=lenr*(rand-0.5d0)*2.0d0 call randomnum()y=lenr*(rand-0.5d0)*2.0d0 dxy2=x*x+y*y if(dxy2.le.lenr2)then m=m+1 end if ncount=ncount+1 if(mod(ncount,100).eq.0)then write(10,(I10,F15.6)ncount,4.0d0*dble(m)/dble(ncount)end if end do end,投点法实验源程序,18,结果和分析,(1)总计投点1.0105次(2)该算法收敛，计算值平均值为3.1392,19,例3 定积分计算,这时我们可以随机地向正方形内投点，最后统计落在曲线下的点数M，当总的掷点数N充分大时，M/N就近似等于积分值I。,20,间接蒙特卡罗方法的思想s,当问题可以抽象为某个确定的数学问题时，(1)首先建立一个恰当的概率模型，即确定某个随机事件A或 随机变量X，使得待求的解等于随机事件出现的概率或随 机变量的数学期望值。(2)然后进行模拟实验，即重复多次地模拟随机事件A或随机 变量X。(3)最后对随机实验结果进行统计平均，求出A出现的频数或 X的平均值作为问题的近似解。,该方法是按照实际问题所遵循的概率统计规律，用计算机进行直接的抽样试验，然后计算其统计参数。,直接蒙特卡罗方法的思想,21,“Buffon投针法”计算圆周率。,作业,22,二 随机变量和随机变量的分布,随机变量：是一个不止是一个值的变量（通常是连续的），并且人们可能无法事先预言某一个特定的值。,但是：其分布是可以了解的，假设我们研究某一连续性的变量，由随机变量的分布我可以得到它取某值的概率：,称为u的概率分布密度函数，它表示随机变量u在u到u+du之间值的概率。,称为g(u)的分布函数。,G(u)在区间取值的单调递增函数,23,三 随机变量的独立性,假如我们考虑两个随机变量u和v的分布，则必须引进这两个变量的联合分布密度函数h(u,v)，此时带来的数学问题就更为复杂。,若h(u,v)=p(u)q(v)，则两个随机变量u和v彼此独立。,对如下三个变量,(x,y,z)相互关联。,24,四 期望值、方差和协方差,一个函数f(u)的数学期望值定义为该函数的平均值,称为u的分布函数。,通常u是在a,b区间均匀分布的随机变量，有,f 的数学期望值：,类似地，自由变量u的期望值为u的平均值，有,25,一个函数或变量的方差:,标准误差或均方根误差：方差的平方根。,由于标准误差与其真值有相同的量纲，因而它比方差更具有物理意义。,假如x和y是随机变量，c是一个常数，则：,(1)数学期望是线性算符(2)方差是非线性算符,x,y间的协方差,26,协方差0，正关联协方差0，负关联,注意：,(1)协方差=0,x,y 为独立变量,(2)x,y 为独立变量,协方差=0,27,五 大数法则和中心极限定理,概率论中的大数法则和中心极限定理是蒙特卡洛方法的基础。,1 大数法则反映了大量随机数之和的性质。,如果函数在a,b区间，以均匀的概率分布密度随机地取n个数ui，对每个计算出函数值hui。大数法则告诉我们这些函数值之和除以n所得的值将收敛于函数h在区间a,b的期望值，即,大数法则保证了在抽取足够多的随机样本后，计算得到的积分的蒙特卡洛估计值将收敛于该积分的正确结果。,28,2 中心极限定理,中心极限定理告诉我们：在有足够大，但又有限的抽样数n的情况下，蒙特卡洛估计值是如何分布的。,该定理指出：无论随机变量的分布如何，它的若干个独立随机变量抽样值之和总是满足正则分布(即高斯分布)。,例如：有一个随机变量，它满足分布密度函数f(x)。如果我们将n个满足分布密度函数f(x)的独立随机数相加：,则Rn满足高斯分布。高斯分布可以由给定的期望值和方差完全确定下来。,29,当n充分大时，对任意的，由列维定理知：,这说明，该积分的期望值与蒙特卡罗估计值之差在范围,内的概率为1-。,30,积分的期望值与蒙特卡罗估计值之差在范围,内的概率为1-。,显著水平：，置信水平：1-。,减小蒙特卡罗估计值标准误差的办法：,(1)适当选取最优的随机变量，使其方差变小；(2)提高实验次数；,这里蒙特卡罗方法精度的概念不是通常意义下收敛真值，而是在某一置信度，或者说是某一概率下收敛于真值。,2为蒙特卡罗估计值标准误差，,31,蒙特卡罗方法的特点,优点能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。受几何条件限制小。收敛速度与问题的维数无关。具有同时计算多个方案与多个未知量的能力。误差容易确定。程序结构简单，易于实现。,缺点收敛速度慢。误差具有概率性。在粒子输运问题中，计算结果与系统大小有关。,