电梯须停次数

2744557306

%电梯问题随机模拟实验：
%设有r个人在一楼进入电梯，楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同，
%求直到电梯中的人下完为止，电梯须停次数的期望。

1. clear all\" `- t) ~! L0 w8 i1 `
   
2. n=28;                             %电梯层数1 {- x$ `4 V! ]& }/ I% z
   
3. r=14;                             %登梯人数: j. |- ]( o6 l$ ~1 w7 {* K
   
4. x=zeros(1,n);+ s: `0 ]% O5 c' N; H/ s4 S0 E
   
5. N=5000;                           %模拟次数
   & B8 |+ b+ F' T' `
6. s=0;: F; T' {) X. D2 U9 f
   
7. for i=1:N
   + `: [* d) c! s3 P
8.     for k=1:n+ {( l' S- Q9 O\" n
   
9.         x(k)=0;
   & C9 e  b1 L* c2 t( f$ c% a
10.     end
    \" I) P6 l\" Z6 b1 V& b4 r
11.     for j=1:r                     %对每个人进行模拟
    & X! a& q\" `# |& e1 ^
12.     x(floor(rand*n)+1)=1;         %有人下梯6 `4 C  d; T! e/ y: |- j& Y$ J
    
13.     end
    8 }/ e+ l8 w2 Y
14.     s=s+sum(x);                   %统计每次模拟中下的人数和. D8 X7 U- ~% o) q
    
15. end
    & `/ e\" D4 c. W
16. s=s/N;                            %平均
    6 z7 c% X5 i- F5 [
17. fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
    . j! {* A& x& k+ z6 [8 S+ ]
18. fprintf('模拟解为%4.3f\n',s);

复制代码

, [* D6 X( N/ z$ Z