《计算物理学》课件 第五章 蒙特卡罗模拟方法.pptx

第五章蒙特卡罗模拟方法;2;3;1、蒙特卡罗方法概述;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24;蒙特卡罗方法的基本思想;实验结果;27;28;29;30;31;32;减小方差的各种技巧;34;35;蒙特卡罗方法的优缺点;描述具有随机性质的事物;受几何条件限制小;收敛速度和问题的维数无关;同时计算多个方案与多个未知量;误差容易确定;不足之处：;不足：计算结果与系统大小有关;主要应用范围;2、随机数与伪随机数;46;47;随机数的定义及性质;49;50;51;由于随机数在蒙特卡罗方法中所处的特殊地位，它们虽然也属于由具有已知分布的总体中产生简单子样的问题，但就产生方法而言，却有着本质上的差别。

为了产生随机数，可以使用随机数表。随机数表是由0，1，…，9十个数字组成，每个数字以0.1的等概率出现，数字之间相互独立。这些数字序列叫作随机数字序列。

如果要得到n位有效数字的随机数，只需将表中每n个相邻的随机数字合并在一起，且在最高位的前边加上小数点即可。例如，某随机数表的第一行数字为86397584102…，要想得到三位有效数字的随机数依次为0.863，0.975，0.841。

因为随机数表需在计算机中占有很大内存，而且也难以满足蒙特卡罗方法对随机数需要量非常大的要求，因此，该方法不适于在计算机上使用。

;随机数的产生方法——物理方法（1）;随机数的产生方法——物理方法（2）;55;存在的问题;第一个问题，不能从本质上加以改变，但只要递推公式选得比较好，随机数间的相互独立性是可以近似满足的。

第二个问题，则不是本质的。因为用蒙特卡罗方法解任何具体问题时，所使用的随机数的个数总是有限的，只要所用随机数的个数不超过伪随机数序列出现循环现象时的长度就可以了。

用数学方法产生的伪随机数容易在计算机上得到，可以进行复算，而且不受计算机型号的限制。因此，这种方法虽然存在着一些问题，但仍然被广泛地在计算机上使用，是在计算机上产生伪随机数的主要方法。;计算程序产生的随机数并不是真正的随机数，它们是确定的，但看上去是随机的，且能通过一些随机性的检验，故常称为伪随机数。

发生周期性循环现象的伪随机数的个数称为伪随机数的??期。

从伪随机数序列的初始值开始，到出现循环现象为止，所产生的伪随机数的个数称为伪随机数的最大容量。;59;60;线性同余方法在计算机上的使用;62;63;64;65;66;67;对指数分布的直接抽样;舍选抽样法（acceptance-rejectionsampling);简单舍选抽样法;;抽样效率:;改进的舍选抽样法;改进的舍选抽样法;改进的舍选抽样法;例1:标准正态分布的抽样,x?[-a,a];77;78;79;80;81;82;83;84;85;86;87;88;89;MC模拟;;排队系统：排队系统模拟：M/G/1排队系统;连续系统模拟实例:追逐问题;1.建立平面直角坐标系:A(x1,y1),B(x2,y2),C(x3,y3),D(x4,y4).;v=1;

dt=0.05;

x=[001010];

y=[010100];

fori=1:4

plot(x(i),y(i),.),holdon

end

d=20;

while(d0.1)

x(5)=x(1);y(5)=y(1);

fori=1:4

d=sqrt((x(i+1)-x(i))^2+(y(i+1)-y(i))^2);

x(i)=x(i)+v*dt*(x(i+1)-x(i))/d;

y(i)=y(i)+v*dt*(y(i+1)-y(i))/d;

plot(x(i),y(i),.),holdon

end

end;;某设备上安装有四只型号规格完全相同的电子管，已知电子管寿命为1000--2000小时之间的均匀分布。当电子管损坏时有两种维修方案，一是每次更换损坏的那一只；二是当其中一只损坏时四只同时更换。已知更换时间为换一只时需1小时，4只同时换为2小时。更换时机器因停止运转每小时的损失为20元，又每只电子管价格10元，试用模拟方法决定哪一个方案经济合理？;98;99;100;101;102;103;104;105;106;107;108;109;110;111;112;113;114;115;116;117;118;119;120;121;122;123;124;IsingModel;126;IsingModel;Mean-FieldTheoryofIsingModel;Mean-FieldTheoryofIsingModel;Mean-FieldTheoryof