电梯须停次数

2744557306

%电梯问题随机模拟实验：
%设有r个人在一楼进入电梯，楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同，
0 d9 V: }) y, ~# d9 F! I- x4 ^( u0 ?. l%求直到电梯中的人下完为止，电梯须停次数的期望。

1. clear all
   7 @4 a$ o: b2 y: [* {; g! u
2. n=28;                             %电梯层数) x% ?) M& o9 L1 R% k* m- r
   
3. r=14;                             %登梯人数2 p1 ?) u+ `! z* C$ c
   
4. x=zeros(1,n);
   $ B$ e6 d* n! G) z, I4 l) ^$ O& o
5. N=5000;                           %模拟次数3 |9 x; z2 n! k  ]. z
   
6. s=0;
   9 g- u! N- c& a7 V# @2 N, J+ S
7. for i=1:N1 m' o5 z- R0 g0 K/ I
   
8.     for k=1:n: n; c5 P/ ^7 e' b% H
   
9.         x(k)=0;8 V, @* f7 Z5 @  N6 z
   
10.     end9 I* R5 w\" g  t- y# r
    
11.     for j=1:r                     %对每个人进行模拟
    4 Q\" z\" Q# {* O
12.     x(floor(rand*n)+1)=1;         %有人下梯0 P& W# R1 \, B3 J
    
13.     end6 C; B, ]- z( m
    
14.     s=s+sum(x);                   %统计每次模拟中下的人数和
    \" y% t+ J, i0 o9 H3 E
15. end
    ) S# V1 X, G' ^$ E3 c9 f( }
16. s=s/N;                            %平均9 s' t2 j4 d9 X+ p
    
17. fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));; I+ E0 a5 [5 W  w
    
18. fprintf('模拟解为%4.3f\n',s);

复制代码