电梯模拟问题.doc

电梯问题 有个人在一楼进入电梯，楼上共有层。设每个乘客在任何一层楼出电梯的概率相同，试建立一个概率模型，求直到电梯中的人下完为止，电梯需停次数的数学期望。并对进行计算机模拟验证。 分析 我们把电梯作为考虑对象，电梯每层要么停要么不停，只有这两种情况。而停与不停是随机的，因此可把电梯每层停与不停用一个随机变量表示，从而可定义一个随机变量序列如下： =0时当且仅当第层电梯没有一个人出电梯，每个人在第?层不出电梯的概率为，而每个人出与不出电梯又是独立的，因此所有个人都不出电梯的概率为 则 的数学期望 记，则为电梯停的次数。 的数学期望 计算机模拟 计算机模拟算法思想： 为C语言编程方便，考虑的楼上层的序号记为，该记法不影响结果。 定义整型数组,并将其个元素的初始值置0。计算机随机产生0到的个随机整数，若该数为，则令，表示第层电梯有人要下，从而电梯在该层要停；当有多人在第层要下时，仍用表示。统计数组中为1的元素个数，即为该次模拟实验电梯要停的次数。总共模拟实验次，求其平均值就是电梯所停次数的模拟值。 下面是C语言程序，采用标准C编制，可在Turbo C2.0、Borland C++3.0或Visual C++6.0下运行。 #include stdio.h #include stdlib.h #include time.h #include math.h /* 电梯问题 */ #define N 5000 /* 模拟总次数 */ #define m 10 /* 总共人数 */ #define n 7 /* 总共层数 */ void main() { int i,j,k; int x[n]; double p,q,total,num,err; p=pow(1.0-1.0/n,m); /* 指定层不停的概率 */ q=1-p; /* 指定层要停的概率 */ total=n*q; /* 需听层数的理论值 */ srand(time(NULL)); /* 初始化随机种子数 */ num=0; for(k=0;kN;k++) { for(i=0;in;i++) x[i]=0; for(j=0;jm;j++) { i=rand()%n;/* 模拟第j个人在哪一层出电梯 */ x[i]=1; /* 第j个人若在第i层停则令该层为1 */ } for(i=0;in;i++) if(x[i]==1) num++; /* 统计一次试验中停的层数 */ } /* 完成总共模拟次数 */ num=num*1.0/N; /* 需停次数的模拟计算值 */ err=fabs(total-num)/total; /* 计算相对误差 */ printf(电梯需停次数的期望值=%7.4f 模拟值=%7.4f 相对误差= %8.4f%%\n, total,num, err*100); } 某两次运行结果如下： 第一次结果： 电梯需停次数的期望值=5.5016 模拟值=5.5262 相对误差=0.4473% 第二次结果： 电梯需停次数的期望值=5.5016 模拟值=5.5128 相对误差=0.2037% 从模拟结果来看，期望值与模拟值很接近，相对误差在0.5%以内，可以验证该该理论模型与计算的正确性。