计算材料学概 之 蒙特卡洛方法.ppt

计算材料学概述;主要内容;蒙特卡洛法是什么？;MC的基本思想;确定性系统;Monte Carlo方法：;Monte Carlo方法简史;　1777年法国科学家布丰提出的一种计算圆周率的方法——随机投针法，即著名的布丰投针问题。这一方法的步骤是： 　　1) 取一张白纸，在上面画上许多条间距为d的平行线。 　　2) 取一根长度为l（ld） 的针，随机地向画有平行直线的纸上掷n次，观察针与直线相交的次数，记为m 　　3）计算针与直线相交的概率． 　　布丰本人证明了，这个概率是: 　　p=2l/(πd) ，π为圆周率 ：　　利用这个公式可以用概率的方法得到圆周率的近似值。下面是一些资料 ; 实验者 年代 投掷次数 相交次数 圆周率估计值 　　沃尔夫 1850 5000 2531 3.1596 　　史密斯 1855 3204 1219 3.1554 　　德摩根 1880 600 383 3.137 　　福克斯 1884 1030 489 3.1595 　　拉泽里尼 1901 3408 1808 3.1415929 　　赖纳 1925 2520 859 3.1795 　　布丰投针实验是第一个用几何形式表达概率问题的例子，他首次使用随机实验处理确定性数学问题，为概率论和蒙特卡罗方法的发展起到一定的推动作用。 　　;Monte Carlo方法之随机数的产生 许多计算机系统都有随机数生成函数;实例一、计算π值;;用该方法计算π的基本思路是： ?1?、根据圆面积的公式：?s=πR^2?,当R=1时，S=π。?2、由于圆的方程是：x^2+y^2=1(x^2为x的平方的意思),因此1/4圆面积为x轴、y轴和上述方程所包围的部分。?3、如果在1*1的正方形中均匀地落入随机点，则落入1/4圆中的点的概率就是1/4圆的面积。其4倍，就是圆面积。由于半径为1，该面积的值为π的值。; REAL R,R1,R2,PI ISEED=RTC() N0=0 N=300000 DO I=1,N R1=RAN(ISEED) R2=RAN(ISEED) R=SQRT(R1*R1+R2*R2) IF(R1.0)N0=N0+1 END DO PI=4.0*N0/N WRITE(*,*)PI END;面积的计算 ;MC的优点;MC的应用;蒙特-卡洛模拟的意义;Monte Carlo模拟的步骤：;Monte Carlo算法的主要组成部分;实例二 定积分计算;; REAL Y Y=0 N=300000 ISEED=RTC() DO I=1,N X=RAN(ISEED) Y=Y+X**2/N END DO WRITE(*,*)Y END ;Monte Carlo方法另一个重要问题：随机数;伪随机数：;随机数产生的办法;关于随机数的几点注意;FORTRAN 语言产生随机数的实例;Monte Carlo方法之随机数的产生 许多计算机系统都有随机数生成函数;Finite difference approximation of differential equations;FICK 第二定律;Fick 第二定律——稳态扩散解;应用之二 生日问题;INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1,10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,365 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR==1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 END IF END DO END DO DO I=1,10000 IF(M(I).LE.23) SUM=SUM+1 END DO PRINT *,SUM/10000 END