数学建模-计算机模拟技术

1模拟基础

2MonteCaHo模拟

3模拟模型案例

1模拟基础

1.0模拟的背景、思路

应用领域：

■第二次世界大战期间，J.V.Neumann等人将进行的“中

子扩散”的科研项目取名为"Monte-Carlo"

・运输系统模拟

■摩天大楼安全疏散系统模拟

•国民经济发展模拟

•人口增长系统模拟

•供水系统模拟

•管理系统模拟

•雷达系统模拟

・战争系统模拟

模拟思路：

•“模拟”一对系统抽象建模

・“试验”一根据模型设计算法，编程进行反复试验

•“估计”一根据试验数据

•“收集”一根据试验结果作出判断

1.1模拟的基本知识

1.1.1模拟的概念及作用

•现实系统的数学或逻辑模型可能十分复杂,

例如大多数具有随机因素的复杂系统，其

中的一些随机性因素很难用准确的数学公

式表述，从而也无法对整个系统采用解析

法求解。模拟是处理这类实际问题的有力

工具。

•模拟通常借助于计算机进行。

计算机模拟：在已经建立的数学、逻辑模

型的基础之上，通过计算机试验，对一个

系统按照一定的决策原则或作业规则，由

一个状态变换为另一个状态的行为进行描

述和分析。

模拟的作用：

•对于很难用解析方法加以处理的问题，模

拟是一种有效的技术；

•对建模过程中的假设进行鉴定，对理论研

究的结论加以检验；

•对不同的实现方案进行多次模拟，按照既

定的目标函数对不同方案进行比较，从中

选择最优方案。

1.1.2模拟的分类

通常，模拟时间是模拟的主要自变量。

设计正确的模拟时间推进机理：模拟过程中

应根据系统的特性正确推进模拟时间，使

系统中各要素与发生的事件保持同步。

推进模拟时间的基本方法：

•下次事件法：将模拟时间由一个事件发生

的时间点推进到紧接着的下一次事件发生

的时间点。

•固定时间步长法：模拟时间每次均以相等

的固定步长向前推进，每到达一个新的模

拟时间点需检查相应时间段内是否发生了

事件。需根据实际问题合理设置模拟时间

发生改变的步长。

根据模拟过程中因变量的变化情况进行分类:

1）离散型模拟：因变量在与事件时间有关的具体模

拟时间点呈离散性变化。大多数系统（如排队服

务系统）可采用离散型模拟。

时间推进方法：一般采用下次事件法

应当重点对系统状态可能发生改变的事件进行描述,

并确定这些事件之间的逻辑关系。

因

变

量

A

事件时间时间

图1离散型模拟

采用离散型模拟模型其中发生系

邦r队系统通常

F。，

T客达二

统状态变化的事件两个一是有顾到；

有.

是服务员完成服务将最近•发生述两种事件之

。上

一的时刻设置为下事件发生点就可将服务过

次，

程描述为图幼扃而疑就讲财比

2）连续型模拟：因变量随时间的改变呈连续性变

化。在大多数计算机模拟过程中，按固定的步长

推进模拟时间。

通常需建立一系列的由系统状态变量组成的

状态方程组，以描述状态变量与模拟时间的关系。

3）混合型模拟：因变量随时间的推移而作连续性

的变化并具有离散性的突变，如库存控制系统。

1.1.3模拟的方式

终态模拟：在规定的时间T内进行模拟运行，时间

达到T时，模拟终止。其性能指标明显取决于系

统的初始状态。

稳态模拟：随着模拟时间的推移，系统的性能逐渐

趋于平稳。其目的是研究非终态系统长期运行条

件下的稳态性能，模拟时间的长短取决于能否获

得系统性能的优良估计（可由模拟输出的精度确

定）。

1.1.4模拟的一般步骤

•明确问题，建立模型。正确描述待研究问题，明

确规定模拟的目的和任务，确定衡量系统性能或

模拟输出结果的目标函数，然后根据系统的结构

及作业规则，分析系统各状态变量之间的关系，

以此为基础建立所研究的系统模型。

•收集和整理数据资料。模拟技术的正确运用，往

往由大量的输入数据作依靠。在随机模拟中，应

认真分析具体收集到的随机性数据资料，确定系

统中随机性因素的概率分布特性，以此为依据产

生模拟过程所必需的抽样数据。

•编制程序，模拟运行。

•分析模拟输出结果：模拟结果的统计特性（样本

均值、方差、置信区间等），灵敏性分析，选择

最优方案。

注：模拟结果的统计分析模拟的输出结果是分布特

征未知的随机变量，每次运行的结果仅仅是对该

随机变量所有观察值总体的一次抽样，对总体的

代表性很差，虽然可以增加模拟运行的时间从而

增加抽样次数，但这些数据总是由一个“种子”

经过一定的算法而获得的伪随机序列，它们是自

相关的，并不能构成统计上独立的随机样本。

1.2随机模拟案例：赶上火车的概率

【问题】如图，一列火车从A站开往B站，某人每天

赶往B站上这趟火车.

工/某人

火车运行方向/

A°--------------------°B

思考：请研究他能否赶上这趟火车。

他已了解到：

1）火车从A站到B站的运行时间是均值为30分

钟，标准差为2分钟的随机变量；

2）火车在下午大约1点离开A站，离开时刻

的频率分布如下：

出发时刻午后1:00午后1:05午后1:10

频率0.70.20.1

他到达B站的时刻的频率分布为

时刻午后1:28午后1:30午后1:32午后1:34

频率0.30.40.20.1

他能否及时赶上火车?

明确问题：他能及时赶上火车的概率是多少?

建模方向（思路）：

i）分析法：用概率统计知识建立分析模型，求解

析解。（思考）

ii）模拟法：用概率统计知识建立模型，通过模拟

求近似解。

即先建立模拟模型，然后通过计算机模拟得到

问题的近似解。在同样条件下多次试验，计算

他能及时赶上火车的频率。

问题分析：能及时赶上火车的充要条件是:

TH

即心―(<T2

其中

彳；一火车从A站出发的时刻;

T2一火车的运行时间；

走什么变

丁3一他到达B站的时刻。重？如何

模拟?

基本假设:

i)假设T1，T,T?都是相互独立的随机变量;

JLJ24T

ii)将午后1时记为t=0,设火车运行时间T?

服从正态分布：T2〜N(30,22)o

建立模型：为了简化计算，将下午L点、记为初始时

刻。得到随机变量「和人的分布律如下：

火车出发时刻！\和人到达B站时刻T3的分布

律分别为：

T1（分）0510

P(t)0.70.20.1

丁(分)

T3（分）2288303033223434

~~P⑴(tj00^.33004.400^.22O0H.1

能及时赶上火车的概率〃=P{T3VTi+T2}

如果r为在(0J)均匀分布的随机数，为了模拟随

机变量%和T3,可以通过如下方法：

0,0<r<0.7,

4=<5,0.7<r<0.95

10,0.9<r<1.贝和，3可分另U

用来模拟随机变

〃28,0<r<0,3

9量％和T3。

30,0.3<r<0,7,

/3=<

32，0.7<r<0,9?

34,0.9<r<1.0.

模拟算法设计

输入：赶火车次数（天数）

输出：赶上火车的频率

两种不同风格的算法描述

主要变量说明：

n模拟次数

k临时变量，存储当前累计模拟次数

count存储赶上火车的次数

算法1（分步骤描述）

第1步输入模拟次数n

第2步k=l,count=0

第3步当k〈=n,执行Step4-9,否则执行第11步

第4步生成均匀分布随机数赋给r

第5步由r及公式确定T1模拟火车出发时刻

争6步生成均匀分布随机数赋给r;.

第7步由r及公式确定T3模拟人达到时刻

第8步生成正态分布随机数T2模拟火车运行时间

第9步IFT1+T2>T3,count=count+l,END

第10步执行第3步

第H步输出赶上火车频率p=count/n

算法2（伪代码描述）

i）初始化：

输入模拟次数n;

count=0;

道）模拟11次

fori=lton,

模拟随机变量£,丁2，丁3,

分别赋给tl,t2,t3;

iftl+t2>t3,

count=count+l

endif

endfor

app_prb=count/n;

n=input「输入模拟次数：1）;

count=0;

fori=l:n,

模拟rtl=rand;当模拟随机变量tl（火车从A站出发的时刻）

程序ifrtl<0.7

Tl=0;

elseifrtl>=0.7&rtl<0.9

Tl=5;

else•exid

T2=30+randn*2;%模取随机菱量t2（火车的运行时间）

%模拟随机变量t3（他到达B站的时亥lj）

rt3=rand;

ifrt3<0.3

T3=28;

elseifrt3>=0.3&rt3<0.7

T3=30;

elseifrt3>=0.7&rt3<0.9

T3=32;

else

T3=34;

end

ifT3<T1+T2zcount=count+l;end

end%for

prob=count/n

模拟结果：每次模拟1000次或5000次

序号模拟次数近似概率p

110000.6280

210000.6920

310000.6530

450000.6490

550000.6260

650000.6288

系统模拟注意事项：

•一次模拟结果毫无意义！

•模拟是试验性的,是思维结果的验证。

•必须进行足够多次的模拟，并对结果进行统计分

析。

系统模拟特点：系统模拟是研究系统，特别是动态

系统的重要方法，对于：

•结构复杂的系统；

•很难用解析方法求出变量关系的系统；

•内部机理不明的“黑箱”系统；

•为验证用其他方法建立的模型及结果,应是较好

的选择。

1.3随机变量的建模

•利用理论分布，基于对问题的实际的、合理的假

设，选择适当的理论分布模拟随机变量

优点是给出了各种理论结果出现的概率，便于进行数学分

析和处理。但此方法仅限于十分简单的情况，当问题越复

杂，数学处理变得越困难，并且丢失了试验数据的信息。

•基于实际数据的频率做近似模拟

优点在于完全与观察数据相符，并且随实际问题的复杂程

度增大不会产生更大的困难，仅增大工作量而已。缺点是

不便于进行数学分析，不得不依赖于模拟得到的统计结果。

・应用常将两种模拟方法结合起来使用

1.4均匀分布随机变量模拟

1.4.1平方取中法

。1.从一个4位数对开始，称为种子；

。2.将它平方得到一个8位数(必要时前面加0)；

。3.取中间的4位数作为下一个随机数。

。按照如上方式，就能得到0—9999随机出现的整

数。因而可以用来模拟(0,1)上均匀分布的随机

数。

。缺点：会退化为0或几个数的循环，所以现在一般不

采用这种方法

。思考：如何皆断随机数生成算法的好坏？

■functionsimrandl

%平方取中法:会退化为0

x0=6532161%2040,2041不好

%653217出现循环

n=6;%n位整数

count=0;disp(sprintf(n'))

whilexO〜=0%fori=l:100,

t=x02;t=mod(t,10(n+n/2));

t=fix(t/10(n/2));disp(sprintf('%10d\t))

x0=t;

count=count+1;

end

1.4.2线性同余法

。给定4个正整数包仇gio

。产生规则：

一叶1=（。/九+6）modc,c可以取很大，如c=231。

•特点（问题）：

。1.序列不断重复（循环周期至多「个数）：

❷2.序列各数之间缺少统计独立性

。3,序列是个确定性

O两组在实践中效果较好的参数设置：

七7计1=（5。/）mod-31）

工片+531

1=(8<rn)mod(2—1)

线性同余法程序：产生o—c整数

。产生n个随机整数

•functionr=randint(x,aFb,c,n)

%输入x为种子

%线性同余法

fori=l:n,

x=mod((a*x+b),c);

®=x;

end

。调用例子

。x=7;a=l,b=7,c=10;n=20;

r—randint(x,a,b.c,n)

。x=7;a=l,b=7,c=231;n=20;

r=randint(x,a,b,c,n)

。思考题：利用线性同余法实现如下功能

。产生(0,1)上的随机数；

。产生U(a,b)均匀分布的随机数；

。产生0,1，一n的随机整数；

。产生从j到k(jWk)上的随机整数。

。思考解答

(设randint为线性同余法随机数生成函数)

❶(1)randint/c;

■(2)a+(b-a)*randint/c;

。(3)INT((n+l)*randint/c)

•(4)j+INT((k-j+l)*randint/c)

1.5Matlab随机模拟函数:

1.5.1常见分布随机变量的模拟

♦1.产生mxn阶［0,1］均匀分布的随机数矩

阵：rand(m,n)

。思考题：请问如何用rand产生均匀分布U(a,b)的

随机数？

♦2.产生mxn阶［a,b］均匀分布U(a,b)的随

机数矩阵：unifrnd(a,b,m,n)

。思考题：在Matlab中用rand和unifrnd分别产生10个均

匀分布U(0,5)的随机数。

。3.产生mxn阶均值为〃，标准差为。的正态分布的

随机数矩阵：normrndn)

/⑺=$expT(号丹/wR

V27T(7/°

❷randn(m,n)产生标准正态分布的随机数

。当研究对象视为大量相互独立的随机变量之和，且其

中每一种变量对总和的影响都很小时，可以认为该对

象服从正态分布。

。机械加工得到的零件尺寸的偏差、射击命中点与目标

的偏差、各种测量误差、人的身高、体重等，都可近

似看成服从正态分布。

❷4.产生mxn阶期望值为〃的指数分布的随机数矩

阵：exprnd（出m,n）

。若连续型随机变量X的概率密度函数为

。其中X>0为常数，则称X服从参数为人的指数分

布。

0指数分布的期望值为1/八

。排队服务系统中顾客到达率为常数时的到达间

隔、故障率为常数时零件的寿命都服从小指数分

布。

。指数分布在排队论、可靠性分析中有广泛应用。

•注意：Matlab中，产生参数为八的指数分布的命

令为exprnd(l/A)

顾客到达某商店的间隔时间服从参数为0」的指数分布

指数分布的均值为1/0.1=10。指两个顾客到达商店的

平均间隔时间是10个单位时间.即平均10个单位时间到

达1个顾客.顾客到达的间隔时间可用exprnd(lO)模

拟。

。5.产生mxn阶参数为人的泊松分布的随机数矩

阵：poissrnd(AFmFn)

e设离散型随机变量X的所有可能取值为0J,2,...

，且取各个值的概率为

Xke~x

P(X=A-)="=0J2・・・，

A-!

•其中》＞0为常数，则称X服从参数为人的泊松分

布。

9泊松分布的期望值为人

•泊松分布在排队系统、产品检验、天文、物理等

领域有广泛应用。

。指数分布与泊松分布的关系：

如相继两个事件出现的间隔时间服从参数为A的指数

分布，则在单位时间间隔内事件出现的次数服从参数

为人的泊松分布.即单位时间内该事件出现k次的概率

为：

.Xke~x

p.X=卜）=———.卜=0.12…，

卜！

反之亦然。

(1)顾客到达某商店的间隔时间服从参数为0.1的指数

分布㈡(2)该商店在单位时间内到达的顾客数服从参

数为0」的泊松分布

。说明：

。(1)指两个顾客到达商店的平均间隔时间是io个单位时

间.即平均10个单位时间到达1个顾客.

。(2)指一个单位时间内平均到达0.1个顾客

敌坦克分队对我方阵地实施突袭,其到达规律服从泊

松分布，平均每分钟到达4辆.

(1)模拟敌坦克在3分钟内到达目标区的数量,以

及在第1、2、3分钟内各到达几辆坦克.

(2)模拟在3分钟内每辆敌坦克的到达时刻。

❷求解思路：先明确随机变量及其分布

。(1)每分钟到达车辆数：4用poissmd(4)进行模

拟。入=4

。(2)两辆车到达的平均间隔时间0.25分钟.

一坦克到达的间隔时间应服从参数为4的负指数分

布，用exprnd(1/4)模拟。

。练习：如何用Matlab编程实现上述问题的模拟？

1.5.2其它随机变量的模拟

X01

概率0.50.5

。模拟方法:

()</•<().5.

r〜U(0,l).

0.5<r<1.

o实现：

。u=rand

ifu<0.5,

x=o；

else

X0123

概率p0.30.20.40.1

9模拟方法：

'0.0<7-<0.3,

1.0.3<r<0,5,

・，'=j2.0,5<r<0,9,r〜U(0,l).

3,0.9<r<1,

•u=rand

ifu<0.3,

X=0;

elseifu<0.5%0.3+0.2

X=l;

elseifu<0.9%0.3+0.2+0.4

X=2;

else

X=3;

end

2Monte-Carlo模拟

2.1Monte-Carlo原理

。蒙特卡罗(MonteCarlo)方法，也称为随机模拟

(randomsimulation)

♦基本思想：为了解决数学、物理、工程技术等方

面的问题

。1.首先建立一个概率模型或随机过程，使它的参数等

于问题的解；

•2.然后通过对模型或过程的观察或抽样试验来计算所

求参数的统计特征，最后给出所求解的近似值。

2.2蒙特卡罗法应用

2.2.1求解非线性规划

。对于非线性规划问题:

min/(/).£eE"

((//(X)>0J=1,2,­••,m.

s.t.

\cij<Xj<%,j=1,2,・・・

・用蒙特卡罗法求解的基本思想是：在估计的区域

{(“1"：2,,・・6

.，。几升勺[ajj}j]j1,2,...,n

誓骋若f实验点，然后从试验点中找出可行

再从可行点中选择最小点.

9基本假设

。试验点的第/个分量叼服从对.%］内的均匀分布.

。符号假设

。P：试验点总数：MAXP：最大试验点总数；

。K:可行点总数：MAXK：最大可行点数；

。X*：迭代产生的最优占・

。Q：迭代产生的最小宿i('x*),其初始值为计算机所能

表示的最大数.

。算法思路

。先产生一个随机数作为初始试验点，以后则将上一个试

验点的第j个分量随机产生，其它分量不变而产生一新

的试验点.这样，每产生一个新试验点只需一个新的

随机数分量.当K>MAXK或P>MAXP时停止迭代.

例子1

maxz———+&ii+

(3J：I+12=10

S.t<N1>0

X2>0

9在Matlab软件包中编程，一般需要设计三个M-文

件：

。主程序；

。目标函数程序；

。约束条件测试程序.

例子2（飞行管理模型的蒙特卡罗法求解程

序）

2.2.2估算圆周率

•(1)采用确定性模型估算圆周率

。।可题:请建立确定性模型估算国周率.

。利用定积分的几何意义：估算曲边梯形面积

。(2)采用蒙特卡罗法估算圆周率

平行四边形ABCD的面积为22=4

圆的面积6=7F12=7T

。现在模拟产生在正方形ABCD中均匀分布的点n个，如

果这n个点中有m个点在该圆内，则圆的面积与正方

形ABCD的面积之比可近似为m/n；

°即工仁如今兀仁坦

'4nn

。模拟程序：估算圆周率：

9n=input('请输入产生点的个数：)m=0;

fori=l:n,

x=-l+2*rand;y=-l+2*rand;

ifx~2+y-2<=Lm=m+l;end

end

mypi=4*m/n

a更简洁实现：

9n=10000;xy=-l+2*rand(2,n);

x=xy(l,:);y=xy(2,:);

mypi=4*sum(x「2+y「2j=l)/n

。下面是一组运行结果：

模拟次数”7T的近似

100003.1028

100003.1395

1000003.1418

1000003.1403

。小结

。运用蒙特卡罗法求解问题要先建立概率模型或者随机

过程。

。思考

•请建立更精确的估计圆周率的数学模型。

2.2.3估算定积分

例子3

求定积分J32dl=宗

图3.1:直线〃=2,直线1=1与曲线〃=x2

1)频率法

。分析：，

•条件:

。考虑〃上)20的情形

•构造一个矩形包含了曲边梯形，矩形区域高"满足：

(1>max/(z)

a<i<b

。利用大数定理，事件发生频率依概率收敛与事件发生

的概率。

。具体思路：

。在矩形区域内产生〃(足够大)个点

。落在曲边梯形内的点为m个，则

小⑺一〜竺

1—57

其中矩形区域面积S=(6-a)cl

。则所求定积分

rb

If(x)dx«—(fc—d)d

•n=106;%产生随机点的个数

a=0;%积分下限

b=l;%积分上限

maxfx=l;%函数最大值(估计)

d=maxfx+l;%确定一个d(作为大矩形的高)

m=0;

fori=l:n,

x=a+rand*(b-a);%随机产生x坐标

y=d*rand;%随机产生y坐标

if丫〈=*八2,%在曲线下方

m=m+l;%计算增力口

end

end

s=m/n*d*(b-a)

。运行程序:

/fd.r70.3332160()0()()0()0

Jo

。思考题

。(1)模型中的参数(如矩形区域的高d)对结果有影响

吗？

。(2)这种方法对于函数值可正可负的情形该怎么应用

呢？

2)平均值法

o平均值法是利用随机变量的平均值(数学期望)

来计算定积分1

I=f3dH(3.1)

Ja

。考虑一系列在(0,1)上均匀分布的独立随机变

量：△,&「・・&，则｛/(&)｝，”1,2,…，7?是

一系列相互独立同分布的随机变量，有

1rbI

同/(&)=----/，3心=-----

I)—aLb—a

9由强大数定律知以概率为1成立

n

山】】！£/(&)=£(3.2)

72T8,=]u«

9因此当〃足够大时，可得近似公式

an

1=//(1)"七他一°叶£/出)(3.3)

Jai=l

。由此得到计算/的平均值法；

。(1)产生在(0.1)上均匀分布随机数4.孙…，九；

9(2)令Uq—a+(6—1限2：=1,2，7，；

9(3)计算宁作为/的估计值。

。分析式(3.1)〜缶1(3.3),可发现本质上平均值法是用样

学期望的估计工

9%program:simJni.mean

%蒙特卡罗模拟法求定积分(平均值法)

「input('输入产生随机数个数:')；

ifisempty(n)|(n<10000),

n=10000;

end

s=0;%存储函数值之和

a=0;%积分下限

b=l;%积分上限

fori=l:n,

x=a+(b-a)*rand;%rand(l)

funvalue=x2;%计奠函数宿

S=s+funvalue;%结合函数计算函数值，并累计到变

量s

end

I=(b-a)*s/n

。运行结果：

sim-int_inean

。输入产生随机数个数：iuo

I=

0.33795894941753

练习1

有30个电子部件D「D?,…,。30，它们的使用情况

如下:一开始D、工作:若Di损坏则D2立即开始工

作;若D2损坏则D3立即开始工作,其余类推。

没D的寿命X*i=l,2,…,30）服从参数

为入=0.1（小时）的指数分布且相互独立。

令丁为30个器件的总使用时间，试用M—C（蒙特卡

罗）模拟法求解P{T230}和石（T尸。

3模拟模型案例

例1：某港口提供有足够的泊位供船舶停靠，但现

在仅有一个可供装卸的泊位，船舶先到则先进行

装卸，如果船舶得不到及时装卸而造成的滞期费

为每小时100元。现要弄清该系统的性能，重点

考察船舶进入该港后等待装卸的滞留时间以及卸

位（即装卸用的泊位）的利用率，从而进行经济

效益分析。

首先，对进入该港口的100艘船进行实际调查，记录其活动情况，得到

这100艘船到达港口的时间间隔和装卸时间的分布情况的频数和累积频

率分布。

表1100艘船到达港口的时间间隔频数表

到达间隔56789101112131415161718

（小时）

到达船舶136791011111197654

数（艘）

表2100艘船装卸时间频数表

装卸时间（小时）910111213141516

装卸船舶数（艘）2022191610832

表3船舶到达港口的时间间隔累积分布表

到达间隔56789101112131415161718

累积频率0.010.040.100.170.260.360.470.580.690.780.850.910.961.00

表4船舶装卸时间累积分布表

装卸时间910111213141516

累积频率0.200.420.610.770.870.950.981.00

为了比较准确地反映系统的性能，我们采用稳态

模拟方式，可以考察系统运行360天的情况。假

定港口每天24小时连续工作，因此模拟时间T设

置为8640小时。考虑到船舶到达港口的事件是影

响整个系统状态的主要因素，可以将模拟事件设

置为该事件的发生时刻。

所用变量说明：

t:一艘船到达港口的时间；

tl:前一艘船驶离港口的时间；

td:两艘船到达港口的时间间隔；

ts：当前船舶装卸所需时间；

tw：t时刻所有已到达港口的船舶等待装卸时

间^总、和；

tf:t时刻装卸位已空闲时间总和；

表5某港口对船舶服务的模拟结果

模拟到达船到港船舶滞留船舶平均滞港口支付滞装卸泊位空装卸泊位

次数舶数总时间（小时）留时间（小时）留费（元）闲时间（小时）利用率

1736748010.1674800035396%

272243696.0543690051594%

3737910912.3690190042695%

473657587.8257580038196%

5741750110.1275010034196%

673968409.2668400042295%

经过6次模拟计算，装卸泊位的平均利用率为

95.3%,但到达的船舶的却因得不到及时的服务

而造成平均每艘滞留9.295小时，每年到港船舶总

滞留6843小时，因此港口每年约需支付68万元的

滞留费，这显然是一笔不小的开支。为此，可修

改模拟模型，考虑增设一个装卸泊位，并重新进

行模拟运行，将这些数据提供给决策者一确定是

否需要再新建一个装卸泊位。

练习：某公共汽车站每隔30分钟到达一辆汽车，但

可能有［0,3］分钟的误差，此误差大小与前一辆汽

车的运行无关。汽车最多容纳50名旅客，到达该

汽车站时车内旅客人数服从［20,50］的均匀分布。

到站下车的旅客人数服从［3,7］的均匀分布，每

名旅客下车的时间服从［1,7］秒的均匀分布。旅

客按每30分钟到达12个人的泊松分布到达汽车站,

单列排队等车，先到先上，如果某位旅客未能上