数学建模社区-数学中国
标题:
电梯须停次数
[打印本页]
作者:
2744557306
时间:
2023-12-24 17:19
标题:
电梯须停次数
%电梯问题随机模拟实验:
( s1 `+ Z, {6 x1 M( W! Y3 D& {
%设有r个人在一楼进入电梯,楼上共有n层。设每个乘客在任何一层楼出电梯的可能性相同,
0 X, I6 m# C' ^; X' }* [6 r
%求直到电梯中的人下完为止,电梯须停次数的期望。
clear all
# w7 [: n! G* b
n=28; %电梯层数
; `- e' J) G6 ?2 C1 v9 N- B9 |
r=14; %登梯人数
/ n$ L% E3 N0 ^, z9 J2 n8 A
x=zeros(1,n);
8 E8 ^; w0 z( p: f* a3 }
N=5000; %模拟次数
8 }; d; P7 X7 K6 q2 a
s=0;
* d2 B8 I1 k5 L/ ~9 @
for i=1:N
( }7 r4 m# M# [4 e2 t$ y2 [( Y
for k=1:n
! I# u# s& K7 K$ I e( h
x(k)=0;
) l% N" L0 t) R' n8 t0 R! Q
end
+ |& j2 O' L: k, ?1 W( h- J1 F
for j=1:r %对每个人进行模拟
0 d* O3 n& k* C0 L
x(floor(rand*n)+1)=1; %有人下梯
: ?0 r# d, E, {
end
1 b( Y" G4 }! f6 s
s=s+sum(x); %统计每次模拟中下的人数和
7 ~ o5 m& o; x$ X# @+ _8 T( U
end
" Z' e' @7 m# P! P. {
s=s/N; %平均
+ g r5 @' h1 f
fprintf('解析解为%4.3f\n',n*(1-(1-1/n)^r));
. D5 }" V0 v% E/ J: c, x
fprintf('模拟解为%4.3f\n',s);
复制代码
# t6 i7 }5 j* ?% d4 T( ]! q
dianti.m
2023-12-24 17:19 上传
点击文件名下载附件
下载积分: 体力 -2 点
668 Bytes, 下载次数: 0, 下载积分: 体力 -2 点
售价:
1 点体力
[
记录
] [
购买
]
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5