计算机模拟方法讲述.docVIP

计算机模拟方法 一．四人追逐实验模拟 如图1,在正方形ABCD的四个顶点各有一个人。设在初始时刻时，四人同时出发匀速以沿顺时针走向下一个人。如果他们始终对准下一个人为目标行进，最终结果会如何。作出各自的运动轨迹。 解：该问题可以通过计算机模拟来实现。这需要将时间离散化。设时间间隔为，时刻表示时间 设第个人时刻的位置坐标为： 对前面3个人表达式为： 其中 对第4个人表达式为： 其中 Matlab实现程序run.m如下： %模拟运动 n=240; x=zeros(4,n); y=zeros(4,n); dt=0.05; %时间间隔 v=10; %速度 x(1,1)=100; y(1,1)=0; %第1个人初始坐标 x(2,1)=0; y(2,1)=0; %第2个人初始坐标 x(3,1)=0; y(3,1)=100; %第3个人初始坐标 x(4,1)=100; y(4,1)=100; %第4个人初始坐标 for j=1:n-1 for i=1:3 d=sqrt((x(i+1,j)-x(i,j))^2+(y(i+1,j)-y(i,j))^2); %第i个人和第i+1个人距离 cosx=(x(i+1,j)-x(i,j))/d; %求cos值 sinx=(y(i+1,j)-y(i,j))/d; %求sin值 x(i,j+1)=x(i,j)+v*dt*cosx; %求新x坐标 y(i,j+1)=y(i,j)+v*dt*sinx; %求新y坐标 end %考虑第1,2,3人运动一步 d=sqrt((x(1,j)-x(4,j))^2+(y(1,j)-y(4,j))^2); %第4个人和第1个人距离 cosx=(x(1,j)-x(4,j))/d; %求cos值 sinx=(y(1,j)-y(4,j))/d; %求sin值 x(4,j+1)=x(4,j)+v*dt*cosx; %求第4点新x坐标 y(4,j+1)=y(4,j)+v*dt*sinx; %求第4点新y坐标 end %plot(x,y) for j=1:n plot(x(1,j),y(1,j),x(2,j),y(2,j),x(3,j),y(3,j),x(4,j),y(4,j)) %作点图 hold on %保持每次作图,实现各次图行迭加 end 执行结果见图1 图1 模拟结果图形 二、电梯问题随机模拟 设有个人在一楼进入电梯，楼上共有层。设每个乘客在任何一层楼出电梯的可能性相同，求直到电梯中的人下完为止，电梯需停次数的数学期望。并对进行计算机模拟验证。 计算机模拟算法思想： 楼上层的序号记为，定义整型数组,并将其个元素的初始值置0。计算机随机产生1到的个随机整数，若该数为，则令，表示第层电梯有人要下，从而电梯在该层要停；当有多人在第层要下时，仍用表示。统计数组中为1的元素个数，即为该次模拟实验电梯要停的次数。总共模拟实验次，求平均值就是电梯所停次数的模拟值。 Matlab模拟程序dianti.m N=5000; %模拟次数 n=28; %电梯层数 r=14; %电梯开始进的人数 s=0; x=zeros(n,1); for k=1:N %模拟N次 s1=0; for i=1:n x(i)=0; end for j=1:r %对每个人是否下电梯进行模拟 i=1+floor(rand(1,1)*n); x(i)=1; %第i层有人下 end s1=sum(x); %该次模拟中总共要下的人数 s=s+s1; %统计各次模拟中下的人数 end eq=s/N; %模拟平均值 ei=n*(1-(1-1/n)^r); fprintf(模拟电梯要停次数的理论值%6.2f,模拟值为%6.2f\n,ei,eq); 某两次运行结果如下：11.15,11.18。与理论数值11.17很接近。 三、 矿工选门问题模拟实验 一名矿工陷入一个有三扇门的矿井中，第一扇门通到一个隧道，走2小时可以到达安全区；第二扇门通到另一个隧道，走3小时回到矿井中；第三扇门通到又一个隧道，走5小时回到矿井中。假定矿工总是等可能的选择三扇门的任何一扇门，求到达安全区的平均时间 计算机模拟的Matlab程序：