数学建模社区-数学中国
标题:
电梯须停次数
[打印本页]
作者:
2744557306
时间:
2023-12-24 17:19
标题:
电梯须停次数
%电梯问题随机模拟实验:
& U$ P: O; o( q* ^) O
%设有r个人在一楼进入电梯,楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同,
5 f; K4 w6 ?1 F, R
%求直到电梯中的人下完为止,电梯须停次数的期望。
clear all
" @) g) C& {! \( }- k. ?
n=28; %电梯层数
- o6 Y: W6 z$ S3 H
r=14; %登梯人数
4 {1 s' T6 ?' i& o3 d) t- R
x=zeros(1,n);
2 {3 Y$ I1 r) s6 r, E2 b3 w( X# X
N=5000; %模拟次数
& u' ~" I; f; X3 U
s=0;
0 [0 F* R Y8 m/ s9 g' t
for i=1:N
; k- S; |& ~/ g
for k=1:n
3 V4 f. {! A* j1 W& j7 Q
x(k)=0;
9 S. s' e6 [7 b: p9 v
end
4 m2 n0 W O$ d$ Z# ]5 ]
for j=1:r %对每个人进行模拟
8 t# i( s- s4 `1 D0 i
x(floor(rand*n)+1)=1; %有人下梯
7 ^' x+ ]( O6 a" b
end
' l5 M/ y5 U( v
s=s+sum(x); %统计每次模拟中下的人数和
! m( a* p+ O( h& X$ b8 |
end
$ s4 Y7 b! O) W2 @" ]
s=s/N; %平均
: W D, o, H4 j; N$ ^; r
fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
# y. B2 `9 x& o; o* Z
fprintf('模拟解为%4.3f\n',s);
复制代码
C9 R8 R- b: {
dianti.m
2023-12-24 17:19 上传
点击文件名下载附件
下载积分: 体力 -2 点
668 Bytes, 下载次数: 0, 下载积分: 体力 -2 点
售价:
1 点体力
[
记录
] [
购买
]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5