数学建模社区-数学中国

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

作者: 2744557306    时间: 2023-12-24 17:19
标题: 电梯须停次数
%电梯问题随机模拟实验:4 h. i7 {/ |* W% H# C
%设有r个人在一楼进入电梯,楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同,- P% r2 c' f2 H/ i7 ]
%求直到电梯中的人下完为止,电梯须停次数的期望。
  1. clear all- @( ?5 j/ l: y
  2. n=28;                             %电梯层数
    ) n' V2 T2 }1 W, }" [# w2 W1 A
  3. r=14;                             %登梯人数, @# X# s. Y$ C6 }1 ~* C1 g
  4. x=zeros(1,n);! o0 P6 b. A$ c3 B
  5. N=5000;                           %模拟次数
    ) |5 V# C" v0 t. G/ Q
  6. s=0;
    & x7 ^$ v) P/ x( \
  7. for i=1:N: P% _5 l/ K+ j( Q$ x2 b
  8.     for k=1:n' w$ y8 y7 P, x; I5 B4 I) o
  9.         x(k)=0;
    0 E0 H% y% E! U$ m" \. {/ e
  10.     end# ]$ w0 i! N7 f* [+ A7 r
  11.     for j=1:r                     %对每个人进行模拟
    * ^. j# {& h! [/ B) e9 v0 J
  12.     x(floor(rand*n)+1)=1;         %有人下梯
    0 Q! K- ^0 D2 M* a7 b* y' M
  13.     end. z4 W4 v, P' [2 N! o' ?
  14.     s=s+sum(x);                   %统计每次模拟中下的人数和
    3 J2 E% f! l8 X; `. g
  15. end( A9 e- U7 R/ N7 I; T, u1 g% v
  16. s=s/N;                            %平均- I: P' w: A6 H1 \' C9 Q5 r
  17. fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
    , r" A7 m. D  N$ F% A
  18. fprintf('模拟解为%4.3f\n',s);
复制代码

' x" i7 N' b/ w( u. F  J. g

dianti.m

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

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






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