第29卷第9期系统工程理论与实践
SystemsVd.29.No.9Sept.,2009200睥9月En西neering—Theory&Practice
文章编号:1000-6788(2009)09一0078一06
基于排队论的一个物流模型
吴锦标t,刘再明t,尹小玲2,付延冰3
(1.中南大学概率统计研究所,长沙410075;2.中山大学数学与计算科学学院,广州510275;3.中南大学交通运输工程学院,长沙410075)
摘要利用排队论研究了一个物流模型.货物到达货运站形成一复合泊松过程.每辆货车每次装载货物数量必须刚好为Ⅳ.求得了系统稳态存在的充分必要条件,利用补充变量法求得了任意时刻货运站内的平均货物量以及货车刚到达货运站时货运站内待运的平均货物量.最后给出了一个数值实例.
关键词随机运筹学;排队论;物流;批到达;批服务
中图分类号0122文献标志码A
Logistic8modelbased
LIUonqueueingtheorywuJnbia01,
andz罅nlin91,ⅥNxi瓣lin92,FuYambin矿namc蛐dn衄8portation(1.In8tituteofProbabili锣andStati8tics,CentralSouthUniVe鹅i锄Chan即ha410075,China;2・DepartmentofMathemati馏Scientinccomputation,Sun、,at-蹲nUniver8ity'Gu锄gzhou510275,China;3.Departmentof
En舀neering,centra卫southuniver8ity,ch蚰gBha4l0075,china)
Abst阳ctWb
accordingto
FbrthisacoIl8ideralogistic8modelba8edonqueueingtheo眄.The肺ight8甜rivetothe舭ightteminalcompoundPoi880nproce晒.A行eightvehiclem11stc8r巧jllstNpiece8of如ights
wi乞hatatime.model,thenece硝aryand8umcientconditionforthe8ystemstabiHtyisobtained.W宅cIeriVetheatamean砌ountsof&e逗hts
Keywordsrandom印ochaasweUa8a七a矗eightvehiclearriv出印ochtheⅢethodof8upplement踟了variabl馏.Atlastwe百venumericalex锄ple.the0盯;10百8tic8;bukarrive.bulk一∞r、ricestoch嬲ticoperationsre∞aurch;queueing
1引言
排队论被广泛应用于交通物流领域,特别是批到达、批服务排队模型可以很好地模拟一些交通、货运、物流等现象.对于批到达排队模型,文献【l_6】进行了详细的研究.本文基于排队论,对一个货运站建立排队模型,研究了货物的到达形成一复合泊松过程,每辆货车每次限载货物数量为Ⅳ.不允许货车超载,也不允许不满载就发车,且具有一般分布装车时间的一个物流模型.
2模型描述
1)到达货运站的货物形成—个复合齐次泊松过程,即货物的到达形成一参数为入的泊松过程,在同一时收稿日期:2008-07-18
资助项目:国家自然科学基金(60574002);湖南省研究生创新基金(3340-7423600000I);中南大学研究生学位论文创新基金(3960-71131100003)作者简介:吴锦标(1982-),男,湖南株洲人,博士研究生,主要从事排队系统与库存控制方面的研究-
第9期吴锦标,等:基于排队论的—个物流模型
刻到达的货物数量(批量)是一正整数值随机变量x,其分布是P伍=i)=q,l=l,2,…,且Ex2<oo.各批到达的货物数量相互独立.
2)每辆货车每次装载的货物数量限为Ⅳ.不允许超载也不允许不满载就发车.若货车到达货运站时,货运站内待运的货物数量大于或等于Ⅳ,则该货车装运前Ⅳ件货物离开货运站;若货运站内待运的货物数量小于Ⅳ,则该货车在货运站等待,直到货运站内的货物数萤大于等于Ⅳ时才装运前Ⅳ件货物离开货运站.
3)货车到达货运站时且货运站内待运的货物数量大于等于Ⅳ时,开始装车,装车时间S服从一般分布,其分布函数、密度函数、风险率函数分别记为A@),口(茁),叩@)=T墨杰。s有有限的一阶矩和二阶矩。
4)前一辆货车离开货运站到后—辆货车到达货运站的时间间隔B服从一般分布,其分布函数、密度函数、风险率函数分别记为B(z),6(z),p(z)=T兰警裔.B有有限的一阶矩和二阶矩.
5)以上各随机变量相互独立.
以下我们记F(z)=1一F@),F+(8)=.f≯e“2dF扛).
3系统稳态存在的充分必要条件
将前一辆货车离开货运站时刻看做服务开始时刻,后一辆货车到达货运站时刻看做服务结束时刻,则本模型相当于一个M/G/1排队系统,而M/G/1排队系统的再生点为顾客服务完离开系统时刻,故考虑货车到达货运站时刻嵌入马氏链.令{‰,铭∈Ⅳ}表示第佗辆货车到达货运站时刻,则{M;=Ⅳ(t。+O),铭∈Ⅳ}为状态空间Ⅳ上的嵌入马氏链(见文献[7】).
定理l以{Ⅳk,他∈Ⅳ}表示第礼辆货车到达货运站时货运站内待运的货物数量,则{Ⅳt;,竹∈Ⅳ)是遍历的,当且仅当不等式p=AEx(ES+EB)<Ⅳ成立.
证明由文献【7】可知{Ⅳt。,佗∈Ⅳ)为一个不可约、非周期的马氏链,由For8ter【8】准则:—个不可约、非周期的马氏链是遍历的,当且仅当存在一个非负函数,U),J∈Ⅳ及£>o,使得:
%=E【,(%+1)一,(帆)1%=J】
对所有的自然数歹∈Ⅳ有限,除有限个自然数外几乎所有的歹∈Ⅳ都有巧≤一£.
令,◇)=』,则有:
一JAEx(Es+EB)一Ⅳ,J=1,2,…
4’一1AEx(Es+EB),J=o
故若不等式入Ex(Es+EB)<Ⅳ成立,由For8ter准则知,嵌入马氏链{^k,n∈Ⅳ)为遍历链.
若嵌入马氏链满足Kaplan【8l条件,且对任意J∈Ⅳ,有%<oo,存在Jo∈Ⅳ,当J≥Jo时有吻≥o,则嵌入马氏链{Ⅳr;,礼∈Ⅳ}为非遍历链.由于存在七=1使得当歹<i一而,l>o时,有%=o,其中R=(nJ)为{M;,佗∈Ⅳ)的一步转移概率矩阵,故满足K印lan条件.若有不等式AEx(ES+EB)≥Ⅳ成立,则显然对任意的J≥1,有巧<∞且存在—个如∈Ⅳ对任意J≥如有巧≥0,故马氏链{Ⅳt¨佗∈Ⅳ)为非遍历链,因此若马氏链{^k,竹∈Ⅳ)为遍历链,则不等式AEx(Es+E日)<Ⅳ成立.证毕.
注即要使系统稳态存在,必须相邻两辆货车离开货运站的时间间隔内新产生的平均待运货物数量小于运输车辆最大装载量.
4模型求解
令Ⅳ(t)为时刻t货运站的货物总数量(包括装车时车上的货物量).令C(t)为系统在时刻t所处的状态.C(t)=o表示时刻t货车在货运站等待货物的到来;C(t)=1表示时刻t在装货;G(t)=2表示时刻t上—辆货车已发车,货物在等待下—辆货车的到来.
定义∈(£)如下:
1)若c∽=1,f(£)表示已逝去的装车时间;
2)若C(t)=2,毒(t)表示已逝去的运输时间,则{c(t),Ⅳ(t),∈(t),t≥o)是—个向量马氏过程(见文献【9】).
系统工程理论与实践第29卷
记
K@)=P{C@)=O,Ⅳ(t)=竹,,t≥0,0≤n<Ⅳ,
t≥0,ⅣSn,&(z,£)dt=P{C(t)=1,Ⅳ(£)=扎,z<∈(t)≤z+dz),
jk(z,t)dt=P{G0)=2,Ⅳ(t)=扎,z<∈(t)Sz+dz),t≥0,0≤佗.
由于货物的到达过程为Poi8son流,由B1lrke定理【10】知,马尔可夫过程{C(t),Ⅳ(t),专(如亡≥o}的稳态概率分布存在当且仅当AEx(ES+EB)<Ⅳ成立.
当AEx(ES十EB)<Ⅳ时,我们记
%2熙K(t),&(z)2。魄&(z,£),R(z)2。骢R(z,t)
根据过程的转移规律,可得到系统的稳态方程组:
A%=&。≥1,∑蛔K—t+五R(。)p(z)血,o≤佗<Ⅳ
掣=_(A州圳坼)+如Ⅳ,葛AQ鼬(巩扎≥Ⅳ
掣-_(A州圳砷)+沁。莲A色%㈤,n≥o
&(o)=正R(z)p(z)妇+.,o∑AQK“礼≥Ⅳt:o焉+1
j名(o)=/i%+Ⅳ(z),7(z)dz,n≥o
Ⅳ一1∞^ooo。,∞
三K+,三Z&@)如+至上R(甸血=ln=0n=Ⅳo”n=0。”
Ⅳ一l∞∞oo
y(z)=∑严K,s(z,名)=∑扩&(z),F(邓)=∑z”R(z),c(名)=∑q∥n=On=Ⅳfl=0t=l
望鱼菩≥三!=一(A+叼(z))s(z,石)+久c(z)s(z,z)
皇墨羞三叠皇=-一(入+p(z))F(z,名)+入c(名)F(z,z)
F(o,z)zⅣ=/s(z,z)t7(z)dz
s(o,z)=队c如)一刈y0)+/F0,z)p0)(k
S@,z)=S(o,z)唧“入c(z)一Nz,A(z)
F0,z)=F(o,z)唧“入c(z)一刈z)百(z)
F(O,z)zⅣ=S(0,z)A’Q—Ac0))
s(0,z)=(AC0)一刈y(z)+F(o,2)口+Q—AG0))(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)为了解上述方程,我们定义下面的母函数:由(2)式得到:由(3)式得到:由(5)式得到:由(1)式和(4)式得到:解(7)式和(8)式分别得到:将(11)式代入(9)式得到:将(12)式代入(10)式得到:
第9期吴锦标,等:基于排队论的—个物流模型81联立(13)式和(14)式可解得:
脚,=拦‰哉焉箫当‰
由(13)式和(15)式得到:
将(15)式代入(12)式得到:∞,㈣,
将(16)式代入(11)式得到:脚,=避学拦蒜裂湍舞骘铲
%,=唑若觜等黜等篙秽
J‰=∑仇‰一七,n=o,1,…,Ⅳ一1嘶)=正高瓮恭篙‰硒∞,∞,下面我们求y(z):记在货车等待货物期间,某批货物到达货运站时,货运站已有n件货物的概率为‰,则‰满足:
I丌0::爿
所以%,n=0,1,…,Ⅳ一1,满足:
%=凰‰,佗=O,1,…,Ⅳ一1.其中‰是归一化常数.为了确定梳,我们利用归一化条件(6)得到:‰=—T{芝}1,其中p=入Ex(Es+EB).
Ⅳ【点。叫
从而
(Ⅳ一p)∑扩‰
y(z)=(19)
Ⅳf蔓1‰1Ln=oJ
将(19)式代入(17)式和(18)式可分别得到:
s(z)=/s(z,名)血
F(z)=/F(z,z)血
5相关指标==茄一・箫—----—-———————--———————=——-—-一・一==搿一・箫———————————————————’・一A+(A一入c(z))B+(A一入c(z))一zⅣ。,f≮二1—1㈣,IZI,-、_叫A+(A一入c@))B+(入一入c(z))一名Ⅳ。,『。芒1.1㈣,、。叫●Z_l
=茄一・箫2万丙五瓦瓦两再可不脬‘—了辱丁∞,‘””J
82系统工程理论与实践第29卷
定理3在稳态下,货车刚到达货运站时,货运站内待运的货物数量分布的母函数
证明记Ⅱj是货车刚到达货运站时,货运站内待运的货物数鼍恰好等于歹的概率,则
’码=三。铲乃(z)p(z)dz,其中‰是归一化常数.由(17)式得到:弘,=氅帮一・脊∞,
瓢力3篆‰Z杩@江(州一‰Z脚∥m(功如
利用归一化条件三(1)=1,得到:Lo=菇r_,将此表达式代入上式,立即得到(23)式.
推论l令L是稳态下,任意时刻货运站内待运的平均货物数量,则
L=掣L幽(Ⅳ)+脚肼
Ⅳ一l咄学粉一・箫入EⅨ(x一1)】(ES+EB)+入2(Ex)2(ES2+EB2+2E。SEB)+p—pⅣ
∑Jq
∑口n其中Lo(Ⅳ)=专}.n=0
推论2令L+是稳态下,货车刚到达货运站时货运站内待运的平均货物数量,则
肚掣L幽(Ⅳ)+酬x-1)/2跚。舛
AE[x(x一1)】(Es+EB)+入2(Ex)2(Es2+EB2+2EsEB)一Ⅳ(Ⅳ一1)
6数值实例
首先假设货物批昔服从几何分布,即c知=(1一p)p知~,后≥1.则可得到它的母函数c(z)=车碧.再假设s和B都服从指数分布,其均值Es=击,E口=击,有有限二阶矩Es2=寿,EB2=磊.令pl=10,p2=1,p:o.5.则可得到Ex:2,Ex2=6,Es=o.1,Es2=o。02,刀B=1,EB2=2,%=(≥)‘.
下面设入=5,由于p=AEX(ES+EB)=11,所以为了满足p<Ⅳ,取Ⅳ>11.得到图1、图2.
图1图2
第9期吴锦标,等:基于排队论的—个物流模型
由图1可以看出,当1l<Ⅳ≤22.5967时,任意时刻货运站内待运的平均货物数量厶随着Ⅳ的增长而急速减少;当Ⅳ>22.5967时,任意时刻货运站内待运的平均货物数量L随着Ⅳ的增长而增加;当Ⅳ=22.5967时,任意时刻货运站内待运的平均货物数量L最少陋≈11.618).
由图2可以看出,当Ⅳ比较小时,货车刚到达货运站时货运站内待运的平均货物数量L+随着Ⅳ的增长而急剧减少;当Ⅳ比较大时,货车刚到达货运站时货运站内待运的平均货物数量L+下降很平缓,且有Ⅳ_∞Um三+=6.
下面设Ⅳ=50,由于p=AEx(Es+EB)=2.2A,所以为了满足p<Ⅳ,我们取A≤22.得到图3,图4.
^(Ⅳ=50)A(Ⅳ;5回
图3圈4
由图3可以看出,当A≤6.0061时,任意时刻货运站内待运的平均货物数量L随着A的增长而缓慢减少;当A>6.0061时,任意时刻货运站内待运的平均货物数量L随着Ⅳ的增长而增加;当A=6.0061时,任意时刻货运站内待运的平均货物数量三最少(L≈21.20335).
由图4可以看出货车刚到达货运站时货运站内待运的平均货物数量工+随着A的增长而增加.
参考文献
阻
忙choudln町G,PaulM.Abatch舭iⅦlqueue丽th锄additionalaseIvicechannellmderN-policy【J】.AppHedtime8withMathematic8andComputation,2004,156:115—130.Arumugana七hanR,Je)rakumarS.Steadysta土eanaly8i8of
N—poHcyandbull【queuewithnndtiplevacation8,8etupclo∞downtim鹤【J】.AppliedMathematicalModemng,2005,29:972-986.
降Kri8hna毗ddyGV,NadarajanR,Ar砌uganathanR.Analysisofabulkqueue而thN—poucynlultipleVacatio璐
and8etuptime8【J】.Computer8&Operatio邶R鹤e缸ch,1998,25:957_967.
阻Leess,L∞Hw,YoonsH,eta1.Batch哪ivaLlqueuewithN-policyandsinglevacation【J】.Computer8&
operatio姗Re∞缸ch,1995,22:173—189.
瞄
陋MadanKC,Al-Ra_聊ashM.OntheM£/G/1queue谢th妣dback缸d叩tionm∞rverVacatio璐b鹪edonasin酉ev8ucationJau—chuanpolicy【J】.AppHedMathematicsandcomputation,2005,160:909_919.Ke.AnM…/G/1system诵th8tartup舱rverand‘,additionaloptio璐for∞Ⅳice[J】.AppuedMath-
ematicalModelung,2008,32:443_458.
p邓永录.运筹学中的随机模型【M】.广州:中山大学出版社,1996.
DengYL.Thestocha8ticModelsin0peratioIl8m税盯ch【M】.Gu卸留hou:sunYat-跎nunive粥ityPre蹈,1996.降sennotLI,HumbletPA,Twe砌ieRL.Me趾dri托sandthenon-ergodicityofMarkovchaLi瑚【J】.operation8
m娉earch,1983,31:783_789.
p史定华.随机模型的密度演化方法【M].北京:科学出版社,1999.
shiDH.TheDen8ityEvolvementMethodofstoch鹪ticMod幽【M】.Beijing:sciencePubushingPre鹄,1999.坶c00perRB.IntroductiontoQueueingTheo盯[M】.NewYbrk:North—Holla∞止1981.
基于排队论的一个物流模型
作者:
作者单位:吴锦标, 刘再明, 尹小玲, 付延冰, WU Jin-biao, LIU Zai-ming, YIN Xiao-ling, FU Yan-bing吴锦标,刘再明,WU Jin-biao,LIU Zai-ming(中南大学,概率统计研究所,长沙,410075), 尹
小玲,YIN Xiao-ling(中山大学,数学与计算科学学院,广州,510275), 付延冰,FU Yan-
bing(中南大学,交通运输工程学院,长沙,410075)
系统工程理论与实践
SYSTEMS ENGINEERING—THEORY & PRACTICE
2009,29(9)
2次刊名:英文刊名:年,卷(期):被引用次数:
参考文献(10条)
1.Krishna Reddy G V;Nadarajan R;Arumuganathan R Analysis of a bulk queue with N-policy multiplevacations and setup times 1998
2.Arumuganathan R;Jeyakumar S Steady state analysis of a bulk queue with multiple vacations,setuptimes with N-policy and closedown times 2005
3.Choudhury G;Paul M A batch arrival queue with an additional service channel under N-policy[外文期刊] 2004(1)
4.Cooper R B Introduction to Queueing Theory 1981
5.史定华 随机模型的密度演化方法 1999
6.Seunot L I;Humblet P A;Tweedie R L Mean drifts and the non-ergodicity of Markov chains 1983
7.邓永录 运筹学中的随机模型 1996
8.Jau-Chuan Ke An M[X]/G/1 system with startup server and J additional options for service 2008
9.Madan K C;Al-Rawwash M On the MX/G/1 queue with feedback and optional server vacations based on asingle vacation policy 2005
10.Lee S S;Lee H W;Yoon S H Batch arrival queue with N-policy and single vacation 1995
本文读者也读过(2条)
1. 周智勇.陈峻.王炜.Zhou Zhiyong.Chen Jun.Wang Wei 基于排队论的出行车辆停放接受条件[期刊论文]-东南大学学报(自然科学版)2006,36(4)
2. 彭绪亚.刘长玮.刘国涛.伍翔.贾传兴.邓镓佳.陈志剑.PENG Xu-ya.LIU Chang-wei.LIU Guo-tao.WU Xiang.JIAChuan-xin.DENG Jia-jia.CHEN Zhi-jian 排队论在垃圾转运站设备优化配置中的应用[期刊论文]-重庆大学学报(自然科学版)2008,31(3)
引证文献(4条)
1.李庆印.谭章禄.于少伟 随机需求下农产品供应链利润分析[期刊论文]-农机化研究 2012(2)
2.孙亮.王晓原.周涛 随机需求下供应链持续盈利能力分析[期刊论文]-技术经济与管理研究 2011(11)
3.邓小瑜.李引珍.赵亚玲 港口航道通过能力研究综述[期刊论文]-水运工程 2011(3)
4.高文杰 基于灰概率的M/M/1/∞排队系统及应用研究[期刊论文]-天津师范大学学报(自然科学版) 2010(3)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xtgcllysj200909011.aspx
第29卷第9期系统工程理论与实践
SystemsVd.29.No.9Sept.,2009200睥9月En西neering—Theory&Practice
文章编号:1000-6788(2009)09一0078一06
基于排队论的一个物流模型
吴锦标t,刘再明t,尹小玲2,付延冰3
(1.中南大学概率统计研究所,长沙410075;2.中山大学数学与计算科学学院,广州510275;3.中南大学交通运输工程学院,长沙410075)
摘要利用排队论研究了一个物流模型.货物到达货运站形成一复合泊松过程.每辆货车每次装载货物数量必须刚好为Ⅳ.求得了系统稳态存在的充分必要条件,利用补充变量法求得了任意时刻货运站内的平均货物量以及货车刚到达货运站时货运站内待运的平均货物量.最后给出了一个数值实例.
关键词随机运筹学;排队论;物流;批到达;批服务
中图分类号0122文献标志码A
Logistic8modelbased
LIUonqueueingtheorywuJnbia01,
andz罅nlin91,ⅥNxi瓣lin92,FuYambin矿namc蛐dn衄8portation(1.In8tituteofProbabili锣andStati8tics,CentralSouthUniVe鹅i锄Chan即ha410075,China;2・DepartmentofMathemati馏Scientinccomputation,Sun、,at-蹲nUniver8ity'Gu锄gzhou510275,China;3.Departmentof
En舀neering,centra卫southuniver8ity,ch蚰gBha4l0075,china)
Abst阳ctWb
accordingto
FbrthisacoIl8ideralogistic8modelba8edonqueueingtheo眄.The肺ight8甜rivetothe舭ightteminalcompoundPoi880nproce晒.A行eightvehiclem11stc8r巧jllstNpiece8of如ights
wi乞hatatime.model,thenece硝aryand8umcientconditionforthe8ystemstabiHtyisobtained.W宅cIeriVetheatamean砌ountsof&e逗hts
Keywordsrandom印ochaasweUa8a七a矗eightvehiclearriv出印ochtheⅢethodof8upplement踟了variabl馏.Atlastwe百venumericalex锄ple.the0盯;10百8tic8;bukarrive.bulk一∞r、ricestoch嬲ticoperationsre∞aurch;queueing
1引言
排队论被广泛应用于交通物流领域,特别是批到达、批服务排队模型可以很好地模拟一些交通、货运、物流等现象.对于批到达排队模型,文献【l_6】进行了详细的研究.本文基于排队论,对一个货运站建立排队模型,研究了货物的到达形成一复合泊松过程,每辆货车每次限载货物数量为Ⅳ.不允许货车超载,也不允许不满载就发车,且具有一般分布装车时间的一个物流模型.
2模型描述
1)到达货运站的货物形成—个复合齐次泊松过程,即货物的到达形成一参数为入的泊松过程,在同一时收稿日期:2008-07-18
资助项目:国家自然科学基金(60574002);湖南省研究生创新基金(3340-7423600000I);中南大学研究生学位论文创新基金(3960-71131100003)作者简介:吴锦标(1982-),男,湖南株洲人,博士研究生,主要从事排队系统与库存控制方面的研究-
第9期吴锦标,等:基于排队论的—个物流模型
刻到达的货物数量(批量)是一正整数值随机变量x,其分布是P伍=i)=q,l=l,2,…,且Ex2<oo.各批到达的货物数量相互独立.
2)每辆货车每次装载的货物数量限为Ⅳ.不允许超载也不允许不满载就发车.若货车到达货运站时,货运站内待运的货物数量大于或等于Ⅳ,则该货车装运前Ⅳ件货物离开货运站;若货运站内待运的货物数量小于Ⅳ,则该货车在货运站等待,直到货运站内的货物数萤大于等于Ⅳ时才装运前Ⅳ件货物离开货运站.
3)货车到达货运站时且货运站内待运的货物数量大于等于Ⅳ时,开始装车,装车时间S服从一般分布,其分布函数、密度函数、风险率函数分别记为A@),口(茁),叩@)=T墨杰。s有有限的一阶矩和二阶矩。
4)前一辆货车离开货运站到后—辆货车到达货运站的时间间隔B服从一般分布,其分布函数、密度函数、风险率函数分别记为B(z),6(z),p(z)=T兰警裔.B有有限的一阶矩和二阶矩.
5)以上各随机变量相互独立.
以下我们记F(z)=1一F@),F+(8)=.f≯e“2dF扛).
3系统稳态存在的充分必要条件
将前一辆货车离开货运站时刻看做服务开始时刻,后一辆货车到达货运站时刻看做服务结束时刻,则本模型相当于一个M/G/1排队系统,而M/G/1排队系统的再生点为顾客服务完离开系统时刻,故考虑货车到达货运站时刻嵌入马氏链.令{‰,铭∈Ⅳ}表示第佗辆货车到达货运站时刻,则{M;=Ⅳ(t。+O),铭∈Ⅳ}为状态空间Ⅳ上的嵌入马氏链(见文献[7】).
定理l以{Ⅳk,他∈Ⅳ}表示第礼辆货车到达货运站时货运站内待运的货物数量,则{Ⅳt;,竹∈Ⅳ)是遍历的,当且仅当不等式p=AEx(ES+EB)<Ⅳ成立.
证明由文献【7】可知{Ⅳt。,佗∈Ⅳ)为一个不可约、非周期的马氏链,由For8ter【8】准则:—个不可约、非周期的马氏链是遍历的,当且仅当存在一个非负函数,U),J∈Ⅳ及£>o,使得:
%=E【,(%+1)一,(帆)1%=J】
对所有的自然数歹∈Ⅳ有限,除有限个自然数外几乎所有的歹∈Ⅳ都有巧≤一£.
令,◇)=』,则有:
一JAEx(Es+EB)一Ⅳ,J=1,2,…
4’一1AEx(Es+EB),J=o
故若不等式入Ex(Es+EB)<Ⅳ成立,由For8ter准则知,嵌入马氏链{^k,n∈Ⅳ)为遍历链.
若嵌入马氏链满足Kaplan【8l条件,且对任意J∈Ⅳ,有%<oo,存在Jo∈Ⅳ,当J≥Jo时有吻≥o,则嵌入马氏链{Ⅳr;,礼∈Ⅳ}为非遍历链.由于存在七=1使得当歹<i一而,l>o时,有%=o,其中R=(nJ)为{M;,佗∈Ⅳ)的一步转移概率矩阵,故满足K印lan条件.若有不等式AEx(ES+EB)≥Ⅳ成立,则显然对任意的J≥1,有巧<∞且存在—个如∈Ⅳ对任意J≥如有巧≥0,故马氏链{Ⅳt¨佗∈Ⅳ)为非遍历链,因此若马氏链{^k,竹∈Ⅳ)为遍历链,则不等式AEx(Es+E日)<Ⅳ成立.证毕.
注即要使系统稳态存在,必须相邻两辆货车离开货运站的时间间隔内新产生的平均待运货物数量小于运输车辆最大装载量.
4模型求解
令Ⅳ(t)为时刻t货运站的货物总数量(包括装车时车上的货物量).令C(t)为系统在时刻t所处的状态.C(t)=o表示时刻t货车在货运站等待货物的到来;C(t)=1表示时刻t在装货;G(t)=2表示时刻t上—辆货车已发车,货物在等待下—辆货车的到来.
定义∈(£)如下:
1)若c∽=1,f(£)表示已逝去的装车时间;
2)若C(t)=2,毒(t)表示已逝去的运输时间,则{c(t),Ⅳ(t),∈(t),t≥o)是—个向量马氏过程(见文献【9】).
系统工程理论与实践第29卷
记
K@)=P{C@)=O,Ⅳ(t)=竹,,t≥0,0≤n<Ⅳ,
t≥0,ⅣSn,&(z,£)dt=P{C(t)=1,Ⅳ(£)=扎,z<∈(t)≤z+dz),
jk(z,t)dt=P{G0)=2,Ⅳ(t)=扎,z<∈(t)Sz+dz),t≥0,0≤佗.
由于货物的到达过程为Poi8son流,由B1lrke定理【10】知,马尔可夫过程{C(t),Ⅳ(t),专(如亡≥o}的稳态概率分布存在当且仅当AEx(ES+EB)<Ⅳ成立.
当AEx(ES十EB)<Ⅳ时,我们记
%2熙K(t),&(z)2。魄&(z,£),R(z)2。骢R(z,t)
根据过程的转移规律,可得到系统的稳态方程组:
A%=&。≥1,∑蛔K—t+五R(。)p(z)血,o≤佗<Ⅳ
掣=_(A州圳坼)+如Ⅳ,葛AQ鼬(巩扎≥Ⅳ
掣-_(A州圳砷)+沁。莲A色%㈤,n≥o
&(o)=正R(z)p(z)妇+.,o∑AQK“礼≥Ⅳt:o焉+1
j名(o)=/i%+Ⅳ(z),7(z)dz,n≥o
Ⅳ一1∞^ooo。,∞
三K+,三Z&@)如+至上R(甸血=ln=0n=Ⅳo”n=0。”
Ⅳ一l∞∞oo
y(z)=∑严K,s(z,名)=∑扩&(z),F(邓)=∑z”R(z),c(名)=∑q∥n=On=Ⅳfl=0t=l
望鱼菩≥三!=一(A+叼(z))s(z,石)+久c(z)s(z,z)
皇墨羞三叠皇=-一(入+p(z))F(z,名)+入c(名)F(z,z)
F(o,z)zⅣ=/s(z,z)t7(z)dz
s(o,z)=队c如)一刈y0)+/F0,z)p0)(k
S@,z)=S(o,z)唧“入c(z)一Nz,A(z)
F0,z)=F(o,z)唧“入c(z)一刈z)百(z)
F(O,z)zⅣ=S(0,z)A’Q—Ac0))
s(0,z)=(AC0)一刈y(z)+F(o,2)口+Q—AG0))(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)为了解上述方程,我们定义下面的母函数:由(2)式得到:由(3)式得到:由(5)式得到:由(1)式和(4)式得到:解(7)式和(8)式分别得到:将(11)式代入(9)式得到:将(12)式代入(10)式得到:
第9期吴锦标,等:基于排队论的—个物流模型81联立(13)式和(14)式可解得:
脚,=拦‰哉焉箫当‰
由(13)式和(15)式得到:
将(15)式代入(12)式得到:∞,㈣,
将(16)式代入(11)式得到:脚,=避学拦蒜裂湍舞骘铲
%,=唑若觜等黜等篙秽
J‰=∑仇‰一七,n=o,1,…,Ⅳ一1嘶)=正高瓮恭篙‰硒∞,∞,下面我们求y(z):记在货车等待货物期间,某批货物到达货运站时,货运站已有n件货物的概率为‰,则‰满足:
I丌0::爿
所以%,n=0,1,…,Ⅳ一1,满足:
%=凰‰,佗=O,1,…,Ⅳ一1.其中‰是归一化常数.为了确定梳,我们利用归一化条件(6)得到:‰=—T{芝}1,其中p=入Ex(Es+EB).
Ⅳ【点。叫
从而
(Ⅳ一p)∑扩‰
y(z)=(19)
Ⅳf蔓1‰1Ln=oJ
将(19)式代入(17)式和(18)式可分别得到:
s(z)=/s(z,名)血
F(z)=/F(z,z)血
5相关指标==茄一・箫—----—-———————--———————=——-—-一・一==搿一・箫———————————————————’・一A+(A一入c(z))B+(A一入c(z))一zⅣ。,f≮二1—1㈣,IZI,-、_叫A+(A一入c@))B+(入一入c(z))一名Ⅳ。,『。芒1.1㈣,、。叫●Z_l
=茄一・箫2万丙五瓦瓦两再可不脬‘—了辱丁∞,‘””J
82系统工程理论与实践第29卷
定理3在稳态下,货车刚到达货运站时,货运站内待运的货物数量分布的母函数
证明记Ⅱj是货车刚到达货运站时,货运站内待运的货物数鼍恰好等于歹的概率,则
’码=三。铲乃(z)p(z)dz,其中‰是归一化常数.由(17)式得到:弘,=氅帮一・脊∞,
瓢力3篆‰Z杩@江(州一‰Z脚∥m(功如
利用归一化条件三(1)=1,得到:Lo=菇r_,将此表达式代入上式,立即得到(23)式.
推论l令L是稳态下,任意时刻货运站内待运的平均货物数量,则
L=掣L幽(Ⅳ)+脚肼
Ⅳ一l咄学粉一・箫入EⅨ(x一1)】(ES+EB)+入2(Ex)2(ES2+EB2+2E。SEB)+p—pⅣ
∑Jq
∑口n其中Lo(Ⅳ)=专}.n=0
推论2令L+是稳态下,货车刚到达货运站时货运站内待运的平均货物数量,则
肚掣L幽(Ⅳ)+酬x-1)/2跚。舛
AE[x(x一1)】(Es+EB)+入2(Ex)2(Es2+EB2+2EsEB)一Ⅳ(Ⅳ一1)
6数值实例
首先假设货物批昔服从几何分布,即c知=(1一p)p知~,后≥1.则可得到它的母函数c(z)=车碧.再假设s和B都服从指数分布,其均值Es=击,E口=击,有有限二阶矩Es2=寿,EB2=磊.令pl=10,p2=1,p:o.5.则可得到Ex:2,Ex2=6,Es=o.1,Es2=o。02,刀B=1,EB2=2,%=(≥)‘.
下面设入=5,由于p=AEX(ES+EB)=11,所以为了满足p<Ⅳ,取Ⅳ>11.得到图1、图2.
图1图2
第9期吴锦标,等:基于排队论的—个物流模型
由图1可以看出,当1l<Ⅳ≤22.5967时,任意时刻货运站内待运的平均货物数量厶随着Ⅳ的增长而急速减少;当Ⅳ>22.5967时,任意时刻货运站内待运的平均货物数量L随着Ⅳ的增长而增加;当Ⅳ=22.5967时,任意时刻货运站内待运的平均货物数量L最少陋≈11.618).
由图2可以看出,当Ⅳ比较小时,货车刚到达货运站时货运站内待运的平均货物数量L+随着Ⅳ的增长而急剧减少;当Ⅳ比较大时,货车刚到达货运站时货运站内待运的平均货物数量L+下降很平缓,且有Ⅳ_∞Um三+=6.
下面设Ⅳ=50,由于p=AEx(Es+EB)=2.2A,所以为了满足p<Ⅳ,我们取A≤22.得到图3,图4.
^(Ⅳ=50)A(Ⅳ;5回
图3圈4
由图3可以看出,当A≤6.0061时,任意时刻货运站内待运的平均货物数量L随着A的增长而缓慢减少;当A>6.0061时,任意时刻货运站内待运的平均货物数量L随着Ⅳ的增长而增加;当A=6.0061时,任意时刻货运站内待运的平均货物数量三最少(L≈21.20335).
由图4可以看出货车刚到达货运站时货运站内待运的平均货物数量工+随着A的增长而增加.
参考文献
阻
忙choudln町G,PaulM.Abatch舭iⅦlqueue丽th锄additionalaseIvicechannellmderN-policy【J】.AppHedtime8withMathematic8andComputation,2004,156:115—130.Arumugana七hanR,Je)rakumarS.Steadysta土eanaly8i8of
N—poHcyandbull【queuewithnndtiplevacation8,8etupclo∞downtim鹤【J】.AppliedMathematicalModemng,2005,29:972-986.
降Kri8hna毗ddyGV,NadarajanR,Ar砌uganathanR.Analysisofabulkqueue而thN—poucynlultipleVacatio璐
and8etuptime8【J】.Computer8&Operatio邶R鹤e缸ch,1998,25:957_967.
阻Leess,L∞Hw,YoonsH,eta1.Batch哪ivaLlqueuewithN-policyandsinglevacation【J】.Computer8&
operatio姗Re∞缸ch,1995,22:173—189.
瞄
陋MadanKC,Al-Ra_聊ashM.OntheM£/G/1queue谢th妣dback缸d叩tionm∞rverVacatio璐b鹪edonasin酉ev8ucationJau—chuanpolicy【J】.AppHedMathematicsandcomputation,2005,160:909_919.Ke.AnM…/G/1system诵th8tartup舱rverand‘,additionaloptio璐for∞Ⅳice[J】.AppuedMath-
ematicalModelung,2008,32:443_458.
p邓永录.运筹学中的随机模型【M】.广州:中山大学出版社,1996.
DengYL.Thestocha8ticModelsin0peratioIl8m税盯ch【M】.Gu卸留hou:sunYat-跎nunive粥ityPre蹈,1996.降sennotLI,HumbletPA,Twe砌ieRL.Me趾dri托sandthenon-ergodicityofMarkovchaLi瑚【J】.operation8
m娉earch,1983,31:783_789.
p史定华.随机模型的密度演化方法【M].北京:科学出版社,1999.
shiDH.TheDen8ityEvolvementMethodofstoch鹪ticMod幽【M】.Beijing:sciencePubushingPre鹄,1999.坶c00perRB.IntroductiontoQueueingTheo盯[M】.NewYbrk:North—Holla∞止1981.
基于排队论的一个物流模型
作者:
作者单位:吴锦标, 刘再明, 尹小玲, 付延冰, WU Jin-biao, LIU Zai-ming, YIN Xiao-ling, FU Yan-bing吴锦标,刘再明,WU Jin-biao,LIU Zai-ming(中南大学,概率统计研究所,长沙,410075), 尹
小玲,YIN Xiao-ling(中山大学,数学与计算科学学院,广州,510275), 付延冰,FU Yan-
bing(中南大学,交通运输工程学院,长沙,410075)
系统工程理论与实践
SYSTEMS ENGINEERING—THEORY & PRACTICE
2009,29(9)
2次刊名:英文刊名:年,卷(期):被引用次数:
参考文献(10条)
1.Krishna Reddy G V;Nadarajan R;Arumuganathan R Analysis of a bulk queue with N-policy multiplevacations and setup times 1998
2.Arumuganathan R;Jeyakumar S Steady state analysis of a bulk queue with multiple vacations,setuptimes with N-policy and closedown times 2005
3.Choudhury G;Paul M A batch arrival queue with an additional service channel under N-policy[外文期刊] 2004(1)
4.Cooper R B Introduction to Queueing Theory 1981
5.史定华 随机模型的密度演化方法 1999
6.Seunot L I;Humblet P A;Tweedie R L Mean drifts and the non-ergodicity of Markov chains 1983
7.邓永录 运筹学中的随机模型 1996
8.Jau-Chuan Ke An M[X]/G/1 system with startup server and J additional options for service 2008
9.Madan K C;Al-Rawwash M On the MX/G/1 queue with feedback and optional server vacations based on asingle vacation policy 2005
10.Lee S S;Lee H W;Yoon S H Batch arrival queue with N-policy and single vacation 1995
本文读者也读过(2条)
1. 周智勇.陈峻.王炜.Zhou Zhiyong.Chen Jun.Wang Wei 基于排队论的出行车辆停放接受条件[期刊论文]-东南大学学报(自然科学版)2006,36(4)
2. 彭绪亚.刘长玮.刘国涛.伍翔.贾传兴.邓镓佳.陈志剑.PENG Xu-ya.LIU Chang-wei.LIU Guo-tao.WU Xiang.JIAChuan-xin.DENG Jia-jia.CHEN Zhi-jian 排队论在垃圾转运站设备优化配置中的应用[期刊论文]-重庆大学学报(自然科学版)2008,31(3)
引证文献(4条)
1.李庆印.谭章禄.于少伟 随机需求下农产品供应链利润分析[期刊论文]-农机化研究 2012(2)
2.孙亮.王晓原.周涛 随机需求下供应链持续盈利能力分析[期刊论文]-技术经济与管理研究 2011(11)
3.邓小瑜.李引珍.赵亚玲 港口航道通过能力研究综述[期刊论文]-水运工程 2011(3)
4.高文杰 基于灰概率的M/M/1/∞排队系统及应用研究[期刊论文]-天津师范大学学报(自然科学版) 2010(3)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xtgcllysj200909011.aspx