电梯排队问题的随机模拟

2744557306

这段 MATLAB 代码是在模拟一天和一百天的服务系统中顾客的等待时间情况。
; g( g; u: j& [% 模拟一天
( |6 Z) F2 V9 L) ~8 [- wclear all
c(1) = exprnd(10);  % 第一个顾客到达的时间
b(1) = c(1);
4 }* R, m# {1 hi = 1;
: ^; e1 j( c3 K0 M: t; V6 }w = 0;

while b(i) <= 480  % 经营时间为480分钟（8小时）
" T' l4 d* A2 z- K; w2 y    e(i) = b(i) + unifrnd(4, 15);  % 随机生成每个顾客的服务时间
    w = w + b(i) - c(i);  % 累计等待时间
) H. p! w, D4 @; \/ S- r, y    i = i + 1;
    c(i) = c(i - 1) + exprnd(10);  % 新的顾客到达时间
' o( k4 s3 W& x    b(i) = max([c(i), e(i - 1)]);  % 开始新的服务时间，选择最大值确保服务不会在到达之前开始
end
0 `* v0 Z# u* z; O3 I6 H' x
average_wait_time_one_day = w / (i - 1)  % 一天的平均等待时间
number_of_customers_served_one_day = i - 1  % 一天内服务的顾客数
$ B% b7 S$ O0 a. z
- c6 E6 @) W$ L, |0 f8 M' x2 ~* ?$ V
& o1 o& s4 L/ Z* u% 模拟一百天
9 n6 E- m) |2 ?for j = 1:100
! Z. y0 v$ X! Z9 {4 g3 i8 W, S    c(1) = exprnd(10);  % 第一个顾客到达
    b(1) = c(1);
, u0 k& m9 {+ D- F3 b    i = 1;
    w(j) = 0;
7 m" ~3 t! Q2 C( k5 D! Y1 i    while b(i) <= 480  % 经营时间为480分钟（8小时）
        e(i) = b(i) + unifrnd(4, 15);  % 随机生成每个顾客的服务时间
        w(j) = w(j) + b(i) - c(i);  % 累计等待时间
        i = i + 1;
        c(i) = c(i - 1) + exprnd(10);  % 新的顾客到达时间
        b(i) = max([c(i), e(i - 1)]);  % 开始新的服务时间，选择最大值确保服务不会在到达之前开始
: e0 C: G  N  m9 X; I    end
+ X& ^2 K. g& ~5 }/ b) _# A8 v; a    t(j) = w(j) / (i - 1);  % 记录当天等待平均时间
    m(j) = i - 1;  % 记录当天服务人数
end

3 |  e1 [, P4 p- Hmean_average_wait_time = mean(t)  % 一百天的平均等待时间
mean_number_of_customers_served = mean(m)  % 一百天内每天的平均服务人数

这段代码首先模拟了一天的服务情况，然后在一个循环中模拟了一百天的服务情况。在每个模拟过程中：
  |' d5 F$ N( g' W" q* i  C1 B
1.顾客的到达时间 c(i) 和开始服务时间 b(i) 是根据指数分布和均匀分布随机生成的，分别表示顾客到达和开始服务的时间。
; N+ q2 p$ f; [8 j* ?* t  S7 {4 Q* A2.循环结束条件是超过了一天的营业时间（8小时，即480分钟）。
3.对于每个顾客，他们的服务时间 e(i) 是在均匀分布 [4, 15] 之间随机生成的，模拟了不同顾客的服务时长。
4 H5 w4 J3 O6 H/ q4 ~0 {4.在一天结束后，记录了当天的平均等待时间 t(j) 和服务的顾客数目 m(j)，这样可以在一百次模拟后计算出平均的等待时间和平均服务人数。
( g) d1 P7 G6 |0 T5 Q7 F
整体思路是通过随机模拟顾客到达和服务时长来模拟服务系统的运行情况，并统计不同情况下的等待时间和服务人数，从而评估系统的性能。

6 M+ F' ~: l0 h* W- t  A