数学建模社区-数学中国
标题:
电梯须停次数
[打印本页]
作者:
2744557306
时间:
2023-12-24 17:19
标题:
电梯须停次数
%电梯问题随机模拟实验:
4 h. i7 {/ |* W% H# C
%设有r个人在一楼进入电梯,楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同,
- P% r2 c' f2 H/ i7 ]
%求直到电梯中的人下完为止,电梯须停次数的期望。
clear all
- @( ?5 j/ l: y
n=28; %电梯层数
) n' V2 T2 }1 W, }" [# w2 W1 A
r=14; %登梯人数
, @# X# s. Y$ C6 }1 ~* C1 g
x=zeros(1,n);
! o0 P6 b. A$ c3 B
N=5000; %模拟次数
) |5 V# C" v0 t. G/ Q
s=0;
& x7 ^$ v) P/ x( \
for i=1:N
: P% _5 l/ K+ j( Q$ x2 b
for k=1:n
' w$ y8 y7 P, x; I5 B4 I) o
x(k)=0;
0 E0 H% y% E! U$ m" \. {/ e
end
# ]$ w0 i! N7 f* [+ A7 r
for j=1:r %对每个人进行模拟
* ^. j# {& h! [/ B) e9 v0 J
x(floor(rand*n)+1)=1; %有人下梯
0 Q! K- ^0 D2 M* a7 b* y' M
end
. z4 W4 v, P' [2 N! o' ?
s=s+sum(x); %统计每次模拟中下的人数和
3 J2 E% f! l8 X; `. g
end
( A9 e- U7 R/ N7 I; T, u1 g% v
s=s/N; %平均
- I: P' w: A6 H1 \' C9 Q5 r
fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
, r" A7 m. D N$ F% A
fprintf('模拟解为%4.3f\n',s);
复制代码
' x" i7 N' b/ w( u. F J. g
dianti.m
2023-12-24 17:19 上传
点击文件名下载附件
下载积分: 体力 -2 点
668 Bytes, 下载次数: 0, 下载积分: 体力 -2 点
售价:
1 点体力
[
记录
] [
购买
]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5