1、数学建模港口问题排队论排队模型之港口系统本文通过排队论和蒙特卡洛方法解决了生产系统的效率问题,通过对工具到 达时间和服务时间的计算机拟合,将基本模型确定在M/M/1排队模型,通过对 此基本模型的分析和改进,在概率论相关理论的基础之上使用计算机模拟仿真 (蒙特卡洛法)对生产系统的整个运行过程进行模拟,得出最后的结论。好。关键词:问题提出:个带有船只卸货设备的小港口,任何时间仅能为一艘船只卸货。船只进港 是为了卸货,响铃两艘船到达的时间间隔在15分钟到145分钟变化。一艘船只 卸货的时间有所卸货物的类型决定,在15分钟到90分钟之间变化。那么,每艘船只在港口的平均时间和最长时间是多少?若一艘船只的
2、等待时间是从到达到开始卸货的时间,每艘船只的平均等待时 间和最长等待时间是多少?卸货设备空闲时间的百分比是多少?船只排队最长的长度是多少?问题分析:排队论:排队论(Queuing Theory),是研究系统随机聚散现象和随机服务系 统工作过程的数学理论和方法,又称随机服务系统理论,为运筹学的一个分支。 本题研究的是生产系统的效率问题,可以将磨损的工具认为顾客.将打磨机当做 服务系统。【1】M/M/1 :较为经典的一种排队论模式,按照前面的Kendall记号定义,前 面的M代表顾客(工具)到达时间服从泊松分布,后面的M则表示服务时间服从 负指数分布,1为仅有一个打磨机。蒙特卡洛方法:蒙特卡洛法蒙
3、特卡洛(Monte Carlo)方法,或称计算机随机模 拟方法,是一种基于“随机数的计算方法。这一方法源于美国在第一次世界大战 进研制原子弹的“曼哈顿计划”。该计划的主持人之一、数学家冯诺伊曼用驰名世 界的赌城一摩纳哥的Monte Carlo-来命名这种方法,为它蒙上了一层神秘色彩。(2)排队论研究的基本问题1排队系统的统计推断:即判断一个给定的排队系统符合于哪种模型,以便根 据排队理论进行研究。2系统性态问题:即研究各种排队系统的概率规律性,主要研究队长分布、等 待时间分布和忙期分布等统计指标,包括了瞬态和稳态两种情形。3最优化问题:即包括最优设计(静态优化),最优运营(动态优化)。3为了得
4、到一些合理的答案,利用计算器或可编程计算器来模拟港口的活动。 假定相邻两艘船到达的时间间隔和每艘船只卸货的时间区间分布,加入两艘船到 达的时间间隔可以是15到145之间的任何数,且这个区间内的任何整数等可能 的出现。再给出模拟这个系统的一般算法之间,考虑有5艘传至的假象情况。对每艘船只有以下数据:船1 船2 船3 船4 船5相邻两艘船到达的时间间隔 20 30 15 120 25卸货时间 55 45 60 75 80因为船1在时钟于t=0分钟计时开始后20分钟到达,所以港口卸货设备在 开始时空空闲了 20分钟。船1立即开始卸货,卸货用时55分,其间,船2在 时钟开始计时后t二20+30二50分
5、中到达在船1与t二20+55二75分钟卸货完毕之 前,船2不能开始卸货,这意味着船2在卸货前必须等待75-50二25分钟。在船2开始卸货之前,船2于t二50+15二65分钟到达,因为船2在t二75分 钟开始卸货,并且卸货需45分钟,所以在船2与t二75+45二120分钟卸货完毕之 前,船3不能开始卸货。这样,船3必须等待120分钟。船4在t二65+120二185分钟之前没有到达,因此船3已经在120+60=180 分钟卸货完毕,港口卸货设备空闲185-180二5分钟,并且,船4到达后立即卸 货。最后,在船4于t二185+75二260分钟卸货完毕之前,船5在t二185+25二210 到达,于是船
6、5在开始卸货前等待260-210=50分钟。模型建立:对于问题中存在的服务系统,建立排队论模型,在仅能为一艘船通过是一个 标准的M/G/1模型:所谓M/G/1模型,就是输入过程为泊松流时,服务时间为任意的条件之下 的,服务机器只有一个得时候。对于M/G/1模型,服务时间T的分布式一般的, (但是要求期望值E(T)和Var(T)方差都存在),其他条件和标准的M/M/1型相 同。为了达到稳态QV1还是必要的,其中有p = E(T)。进入队列 服务台船只到达 o o O . O船只离去接受服务图9-2单服务台单队系统单服务员的排队模型设:(1)船只到来间隔时间服从参数为0.1的指数分布(2)对船只的
7、服务时间服从4, 1 5上的均匀分布.(3)排队按先到先服务规则,队长无限制系统的假设:(1)船只源是无穷的;(2)排队的长度没有限制;(3)到达系统的船只按先后顺序依次进入服务,即“先到先服务”。符号说明w:总等待时间;q :第I个顾客的到达时刻;0 :第I个顾客开始服务时 刻;ej :第i个顾客服务结束时刻;xi :第i-1个顾客与第I个顾客之间到达的间 隔时间;y,:对第i个顾客的服务时间Ci=Ci-l+Xjej 二 bj+yjbj=max(cj1ej-i)模型检验:表1 100艘船港口和系统的模拟结果一艘船呆在港口的平均时间977978818599艘船呆在港口的最长时间17412111
8、1141140159艘船的平均等待时间238591224艘船的最长等待时间994633646893卸货设备空闲时间的百分比0.0670.0790.0930.070.0690.028+系列1线性(系列1)上图为一艘船呆在港口的平均时间+系列1线性(系列1)上图为一艘船的最长等待时间上图为一艘船的最长等待时间以上就是对港口问题的具体分析,其实港口问题还可以从船只的排队角度出 发,我们还可以对多个港口通行做相应的模拟试验,让船主尽量减少等待时间且 港口卸货设备的利用率达到最高,从而是港口的主人获得更大的利润。从排队角 度来解决问题,可以使问题的广度增加,选秘书问题就是一个很典型的例子,可 以从排队角
9、度解决,如果用我在文章中应用的方法来解决也是可以的,这仅仅是一个港口的小问题,甚至可以说是一个非常简单的问题.但是已经 让我感觉到了数学的美,在老师的引导下慢慢接近一种抽象的美,在写论文的这 几天中,数据的整理和分析是最值得享受的时刻,在Excel里输入自己的数据, 是一种忐忑的感觉,因为在那么多的数据面前,我真的不知道将会发生什么,拟 合的过程就更是有意思了,一次一次的尝试,一次一次的比较,在这个过程中, 如果有一点点的进步都会让我兴奋,数学建模在生活中处处存在.如果真的能够 掌握这个本领,生活一定会变得简单而精彩!参考文献:(1)运筹学教材编写组编.运筹学.北京:清华大学出版社,2008(
10、2)Jerry Banks, John S.Carson, Barry L Nelson 等著.离散事件系统仿真.北京: 机械工业出版社,2007vv排队论模型与蒙特卡罗仿真附录一编程如下:clearcs=100;for j=l:csw(j)=0;i=l;x(i)=expr nd(10);c(i)=x(i);b(i)=x(i);while b(i)v二480y(i)=unifrnd(4,15);e(i)=b(i)+y(i);w(j)=w(j)+b(i)-c(i);i=i+l;x(i)=expr nd(10);c(i)=c(i-l)+x(i); b(i)=max(c(i),e(i-l);end
11、i 二 i-1;t(j)二 w(j)/i;m(j)=i;endpt 二 0;pm=0;for j=l:cspt 二pm=pm+m(j);endpt=pt/cspm=pm/cs附录二排队论中一个感兴趣的问题时,当输入过程是Possion流时,顾客相继到达 的间隔时间T服从什么规律,定理设)诃是具有参数2的泊松过程,即PN=畔严”0.2,0.7;.心 是对应的时间间隔序列,则随机变量瑞“ =02.70)是独立同分布的,且服从均 值为的负指数分布,即) = /:; o证明因为是Possion过程中第一个顾客到达的时间,所以时间?;旳等价于0.“内没有顾客到达。故P2 = PN=0=岑_八=宀 进而可得P7;r = l-P7;trs = p在帚+$)内没有顾客到达|7=4=p在u+$)内没有顾客到达 (巾$伽过程的独立性)= PN(+)_N()= O= PAr(f)-Ar(O) = O (Powf期过程的平稳增址性质)= =P7;/即PT2t = l-e ,故7;也是服从均值为A*的负指数分布。2、 对于任意的”21和no有P7;rp=J1=JI,.,7= = (/+-+)-iV(J1+.+JB4)=O = P(f)-Ar(O) = O = 即PTBr = l-e-,所以对任一人(心1),它都服从均值为胆的负指数分布。证毕。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1