数学建模社区-数学中国

标题: 电梯须停次数 [打印本页]

作者: 2744557306    时间: 2023-12-24 17:19
标题: 电梯须停次数
%电梯问题随机模拟实验:
* p! H0 w9 c9 B. i0 n) f; a%设有r个人在一楼进入电梯,楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同,' u# q; A+ a% [! j- ?5 c! b) w
%求直到电梯中的人下完为止,电梯须停次数的期望。
  1. clear all
    * g$ m& S8 v$ V. s: V. a
  2. n=28;                             %电梯层数. w; `8 v& g! E# k2 Q: Y: ^
  3. r=14;                             %登梯人数# f! b/ a0 c+ y; H! o" S
  4. x=zeros(1,n);
    ; J$ b; d; {9 O  N
  5. N=5000;                           %模拟次数
    3 _! V' \7 u( T2 f" \2 s& U$ z- s) W
  6. s=0;- \% @/ ?4 C  M$ J
  7. for i=1:N
    0 h% k3 M  z' j: z: ~& G
  8.     for k=1:n& p! G& x$ T( x( z" k
  9.         x(k)=0;
    3 ~& ]  A2 s9 \9 I: W: e
  10.     end  p, A0 e0 n) J5 \: A. e
  11.     for j=1:r                     %对每个人进行模拟
      `/ h' \' v  \
  12.     x(floor(rand*n)+1)=1;         %有人下梯3 p/ y. r, X$ P% x" j8 M# U
  13.     end
    3 K- o) L- E+ |3 \1 P) F
  14.     s=s+sum(x);                   %统计每次模拟中下的人数和3 K# u+ T: @2 c4 i# T, F) H( X
  15. end
    9 p+ H3 K6 e1 m0 G& W: p+ T
  16. s=s/N;                            %平均# z# A; S* p6 \* j0 H2 i
  17. fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
    ! C$ z$ G" D, O/ f" ^
  18. fprintf('模拟解为%4.3f\n',s);
复制代码
7 H) W9 C. v6 q# U2 {% n, J: N

dianti.m

668 Bytes, 下载次数: 0, 下载积分: 体力 -2 点

售价: 1 点体力  [记录]  [购买]






欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5