排队论模型.docx

1、排队论模型2010年东北大学数学建模培训模拟竞赛参赛队：学号 姓名 班级选择题目：A B东北大学大学生创新中心基于排队论的生产系统分析一、摘要本文为B题的生产系统的优化问题做出了自己的解释，通过排队论和蒙特卡洛方法解决了生产系统的效率问题，通过对工具到达时间和服务时间的计算机拟合，将基本模型确定在排队模型，通过对此基本模型的分析和改进，得出了前三问。最后一问在概率论相关理论的基础之上使用计算机模拟仿真（蒙特卡洛法）对生产系统的整个运行过程进行模拟，得出最后的结论。关键词：排队理论、生产系统、蒙特卡洛方法 二、问题的提出 在科学技术突飞猛进的今天，人类的生产力有了很大程度的提高，目前的系统分析已

2、经不仅仅局限于单纯产量的提高，整个系统资源的利用和效率的高低成为衡量生产抑或服务系统的重要指标，而我们的制造业和服务业种系统的效率问题正在成为人们的焦点。本题目正是从实际的生产系统中抽象出来的一个子系统，原题给出了在只有一台工具维修设备的前提下，20名工人的实际操作结果。维护生产正常进行的各项费用已知，题目要求我们对已有问题建立数学模型并改进：(1) 建立数学模型分析该车间是否应该另外添置一台打磨机(2) 现在车间要求达到平均两个小时内只有一人在等候打磨，利用所建模型分析此要求是否可行，车间需要支出费用多少(3) 车间为了确保整个生产安全有序地进行，就要让在打磨间的人数尽量少于2人，利用所建模

3、型分析此要求是否能够实现，车间需要支出费用多少(4) 如果在16:00打磨机停止工作，多少工人可以提前多长时间下班现在该车间想调整作息时间为8:00-12:00，13:00-17:00，打磨机工作时间为8:00-17:00，模型应该如何调整三、条件假设 （1）基于排队系统的假设输入过程：中顾客的总体是有限的，且到了的方式是一个一个的，工具相继到达的间隔时间是随机性的，但是服从一定条件的概率分布；工具的到达也是相互独立的，就是说，以前的到达情况对以后工件的到来没有影响；输入过程是平稳的，也可以认为是时间齐次的，是指描述相继到达的时间间隔分布和所含参数（如期望值、方差等）都是与时间无关的。排队规则

4、：根据题目的已知数据进行分析，工件到达时，如所有的服务台均被占用着，则工具将会等待，服务机制属于先到先服务；从占有的空间上来看，对于等待队伍的长度没有最大限制；从等待队伍的数量上来看，队伍是单列的。服务机构：带有多个服务台的机构中，它们应该是平行并列的（如第一问）服务的方式为每次一台机器对一个工具进行加工；跟本题的输入过程一样，服务时间也是随机性的，但是服从一定条件的概率分布，并且服务时间也是时间齐次性的。（2）对于到达时刻的假设对于工具到达时间：1、根据排队论和概率论的相关理论，我们易知在不相重叠的时间间隔内工具到达数是相互独立的，即为无后效性，2、在充分小的一段时间间隔之内，在区间内有一个

5、顾客到达的概率与时间无关，而约与时间长度成正比，即其中是常数，它表示单位时间有一个顾客到达的概率，称为概率强度。3、另一方面，对于充分小的，在时间间隔内有两个或两个以上顾客到达的概率极小，以至于可以忽略，即（1）泊松分布的概率如下所示，意为在t的时间间隔中到达n个顾客的概率满足以上三个条件的分布被称作泊松流，我们得知该题目的工具到达时间服从泊松分布。通过计算机的相关假设检验，在0.05的置信度基础之上，可以认为工件的到达时间服从该分布。四、符号约定Kendall记号 X表示相继到达间隔时间的分布，Y表示服务时间的分布，Z表示并列的服务台的数目。A系统容量限制，B为顾客源数目，C为填写的服务规则

6、，例如先到先服务（）。 队长，指在系统中的顾客数 指在系统中排队等待服务的顾客数 指一个顾客在系统中的停留时间，期望值 一个顾客在系统中排队等待的时间，期望值 项目总花费 服务设备总投资 平均到达率 服务机构单位时间的费用，为常数6.67 顾客在系统停留单位时间的费用，为常数10五、问题分析1.名词解释排队论：排队论(Queuing Theory) ，是研究系统随机聚散现象和随机服务系统工作过程的数学理论和方法，又称随机服务系统理论，为运筹学的一个分支。【】本题研究的是生产系统的效率问题，可以将磨损的工具认为顾客，将打磨机当做服务系统。：较为经典的一种排队论模式，按照前面的Kendall记号定

7、义，前面的M代表顾客(工具)到达时间服从泊松分布，后面的M则表示服务时间服从负指数分布，1为仅有一个打磨机。蒙特卡洛方法：蒙特卡洛法蒙特卡洛(Monte Carlo)方法，或称计算机随机模拟方法，是一种基于“随机数”的计算方法。这一方法源于美国在第一次世界大战进研制原子弹的“曼哈顿计划”。该计划的主持人之一、数学家冯诺伊曼用驰名世界的赌城摩纳哥的Monte Carlo来命名这种方法，为它蒙上了一层神秘色彩。（2）在本题的第四问中，我们采用了这种方法。2.问题的背景分析为了研究整个系统的效率问题，提出了排队论的概念，特别是在生产和服务系统中，排队是在日常生活中经常遇到的现象，如顾客到商店购买商品

8、、病人到医院看病常常要排队。现实中的例子还有如电话局的占线问题，车站、码头等交通枢纽的车船堵塞和疏导，故障机器的停机待修（此题目中可以理解为工具的磨损待修），水库水量的库存调节等都是有形或者无形的排队问题。如果增添服务设备（如第一问），就要增加投资或者发生空闲浪费；如果服务设备太少，排队现象就会更加严重，对顾客个人及对社会都会带来不利影响。因此，管理人员必须在这些方面进行平衡，经常检查目前的系统状态是否得当，研究今后改进对策，以期提高服务质量，降低成本。排队论也称为随机系统服务理论，本模型根据最基本的排队理论，对生产系统进行求解，解决其效率问题。3.问题分析根据排队论的相关理论和计算机拟合及系

9、统仿真的相关过程，我们对题目做以下分析：（1）比较及两个系统（此时的G为未知处理时间分布）的总费用问题，可以将费用分为总投资、由于工具的停止运行而产生的生产费用、机器的维修费用三部分组成，得到的结果选取较小的更优。（2）、（3）根据相关排队系统的期望值得到代数方程，并对得到的一元二次方程进行求解，此题在于考察排队系统的相关期望。（4）在计算机拟合的基础上，能够验证在的置信水平下，题目中给出的工具到达时间和服务时间间隔依次符合泊松分布和负指数分布，并且得到相关分布的参数，由蒙特卡洛法对所给出的分布情况进行仿真，根据伯努利大数定理，在样本很大的情况下，总体的概率内涵可以由样本确定，进行几次仿真测试

10、，用期望得到题目的最终解。六、模型建立和问题求解1原始数据的处理根据题目给出的相关数据做出如下的表格：编号中断生产时间开始打磨时间结束打磨时间间隔到达时间服务时间8:00018:078:078:137628:118:138:2841538:198:288:4081248:468:468:58271259:069:069:22201669:499:499:5643779:519:5610:11215810:0610:1110:211510911:2410:2411:307861011:3111:3111:477161111:3511:4712:044171211:4112:0412:196151

11、311:5212:1912:3511161412:0212:3512:4810131512:1212:4813:0710191612:3913:0713:2427171712:4013:2413:401161812:4513:4013:49591912:5913:4913:5914102013:0213:5914:06372模型的建立对于问题中存在的服务系统，建立排队论模型，在仅有一台机器的前提下，这是一个标准的模型：所谓模型，就是输入过程为泊松流时，服务时间为任意的条件之下的，服务机器只有一个得时候。对于模型，服务时间T的分布式一般的，（但是要求期望值和方差都存在），其他条件和标准的型相同。

12、为了达到稳态还是必要的，其中有。在上述条件下，则有上面的七个数只要知道了其中的3个就可以求出其余的数来。这就是着名的公式，下面我们将用这个模型来进行求解。输入过程的确认：对于工具到达时间而言，以分钟作为时间单位，则工具的平均到达率为到达总数/总时间=20/302=0.0662个/分钟，所以得成立。服务过程的确认：由题目中给出的数据得知，=264/20=13.2, =16.9, =0.0662*13.2=0.873841,将得出的常数代入到上面的公式里，可以依次求出我们想要的期望值得=4.19367, 63.34849，50.14849。这与我们假设服务时间为负指数分布变量的情况是基本相同的，由

13、Matlab的分布拟合结果显示，在置信度为0.05的水平上，可以近似得认为服务时间服从负指数分布，且根据拟合结果 0.07874，拟合的Matlab代码见附录二。此时可以将原模型改写为模型。由于按照打磨技术水平要求，平均每15分钟打磨一件工具车间须花费100元，如果想缩短时间，需要提高打磨机辅助用品质量，每一件工具打磨时间缩短1分钟需要多花费10元。对于整个车间，平均每个工人中断生产1分钟，车间效益将损失15元。总费用为在8:00至16:00八个小时的时间内，一共需要被修理的工具数的期望是8*60/(20/302)=301，在有一个整修机器存在的情况下，总费用由机器的购置费用、修理费和由于工具

14、的不能正常使用而导致的生产损失三部分组成。仅仅有一台机器存在时的总费用为3问题的求解（1） 问题：分析该车间是否应该另外添置一台打磨机另外增加一台打磨机的问题将原来的问题转化为问题，虽然刚开始的投资增加了，但是由于两台打磨机的存在，使整个系统中工具的停留时间得到了一定程度上的减少，单个打磨机的负担得到缓解。模型各种特征的规定与标准的模型的规定相同。另外规定各台服务台工作是相互独立（不搞协作）且平均服务率相同，。于是整个机构的平均服务率为2；经检验得1成立；当系统到达平衡时，设有以下的概率状态分布：，将代入到上面的公式中得， ，根据上面一个机器时的相同做法，得因为，所以选择增加一台机器。（2）

15、问题：现在车间要求达到平均两个小时内只有一人在等候打磨，利用所建模型分析此要求是否可行，车间需要支出费用多少两个小时仅仅有一个人在等待打磨，可以确定，还是在模型的基础之上，即平均每分钟等待的工具数目为0.00833333，即平均等待时间为120分/个，由得，在=0.0662时，解之得=0.1474，得到平均服务时间为1/=6.7843,每个至少应该减少的打磨时间为12.7-6.7843=5.9157=6分钟，所以应该花去的费用为310*6=1860元人民币。（3） 问题：车间为了确保整个生产安全有序地进行，就要让在打磨间的人数尽量少于2人，利用所建模型分析此要求是否能够实现，车间需要支出费用多

16、少同样的道理，这仍然可以按照最基本的模型进行求解，而题目中要求的在打磨间要求的等待队长期望=2，所有有=2成立，由=0.0662，解之得=0.0883，则得到平均服务时间为1/= 11.32503，每个至少应该减少的打磨时间为12.7-11.32503=1.37497=2分钟，可以实现，且需要的花费为310*2=620元人民币。（4） 问题：如果在16:00打磨机停止工作，多少工人可以提前多长时间下班现在该车间想调整作息时间为8:00-12:00，13:00-17:00，打磨机工作时间为8:00-17:00，模型应该如何调整在计算机拟合数据的基础上，我们可以验证在置信度为0.05的水平上，工具

17、到达时间和服务时间间隔分别服从泊松分布和负指数分布，而且得到参数分别为，将其作为参数，在本题目的题设下模拟系统的运行过程，将达到的结果在Excel下进行处理，得到最后的结果。经过Matlab语言的编写和调试，我们对这一过程重复了50次（代码详见附录二），最后得出若16:00停止打磨机，则得到平均有7.9200名工人可以提前下班，车间调整作息时间为8:00-12:00，13:00-17:00，打磨机工作时间为8:00-17:00时得到有4.8400工人可以提前下班。跟事实情况也是符合的，这时我们认为得到的数学模型准确性还是可以的。在16:00停止打磨机器的情况下，由于磨损的工具无法得到及时的打磨

18、，使工人无法继续生产，从而使生产效率低下，下面是用Matlab语言进行仿真的结果，我们将得到的数据导入Excel表格中，举出一个仿真的实例： 中断生产时间开始打磨时间结束打磨时间间隔时间服务时间8:168:168:2916138:208:298:444158:378:448:5817148:468:589:129149:019:129:211599:059:219:364159:059:369:510159:179:5110:0412139:5610:0410:18391410:3110:3110:42351110:4710:4711:00161311:0611:0611:24191811:0

19、711:2411:3611211:1711:3611:48101211:2511:4812:0081212:1012:0012:11451112:3012:3012:40201012:3412:4013:0142112:4313:0113:1891712:4413:1813:2911112:5313:2913:4391413:0613:4313:5013713:0713:5014:1612613:1514:1614:238713:1614:2314:3611313:2014:3614:5241613:5314:5215:05331314:1215:0515:30192514:3515:3015

20、:45231514:4815:4515:58131314:5815:5816:02101015:023915:0641515:1481715:2410815:42181715:4861215:557615:59413从上面的表格中可以较为清晰得知道，如果机器在16:00停止，8个工人可以平均提前28.75分钟下班，效率低下。如果以此数据为例，车间调整作息时间为8:00-12:00，13:00-17:00，打磨机工作时间为8:00-17:00，得到的表格如下：中断生产时间开始打磨时间结束打磨时间间隔时间服务时间8:168:168:2916138:208:298:444158:378:448:58

21、17148:468:589:129149:019:129:211599:059:219:364159:059:369:510159:179:5110:0412139:5610:0410:18391410:3110:3110:42351110:4710:4711:00161311:0611:0611:24191811:0711:2411:3611211:1711:3611:48101211:2511:4812:0081213:0613:0613:1313713:0713:1313:3912613:1513:3913:468713:1613:4613:5911313:2013:5914:15416

22、13:5314:1514:28331314:1214:2814:53192514:3514:5315:08231514:4815:0815:21131314:5815:2115:31101015:0215:3115:403915:0615:4015:5541515:1415:5516:1281715:2416:1216:2010815:4216:2016:37181715:4816:3716:4961215:5516:4916:557615:5916:5517:0841316:0121116:2120816:276716:39121016:511217由上面可知，一共有6名员工可以提前下班，平

23、均提前时间为37分钟，较上面有了提高。七、附录1参考文献（1）运筹学教材编写组编. 运筹学. 北京：清华大学出版社，2008（2）Jerry Banks，John S.Carson，Barry L Nelson 等着. 离散事件系统仿真.北京：机械工业出版社，2007（3）薛定宇、陈阳泉着.高等应用数学问题的Matlab求解（第二版）.北京：清华大学出版社，20072 计算机拟合工具到达时间和服务时间间隔的Matlab代码：（1）验证服务时间服从泊松分布X = 6 15 12 12 16 7 15 10 6 16 17 15 16 13 19 17 16 9 10 7;size(X)%可以知道

24、自己输入是20个数字然后借助alam = poissfit (X)，求出alam作为泊松分布的参数r=sort(X);H,p=kstest(r,r poisscdf(r,alam),0.05)结果是：H = 0p = 0.0560当H=0时，根据Kolmogorov-Smirnov检验方法可以知道：X服从泊松分布。（2）验证间隔时间是服从指数分布Y=7 4 8 27 20 43 2 15 78 7 4 6 11 10 10 27 1 5 14 3;size(Y)%可以知道自己输入是20个数字根据排队论知识可以知道，参数可以是：0.0662r=sort(Y);H,p=kstest(r,r exp

25、cdf(r,mean(Y),0.05)结果是：H = 0p = 0.8285当H=0时，根据Kolmogorov-Smirnov检验方法可以知道：Y服从负指数分布。3、对题目的第四问进行仿真首先利用服务时间服从泊松分布，可以利用MATLAB里面现有的函数去产生随机数用于模拟。A=poissrnd(12.7,40,1);在利用时间间隔是服从指数分布，B=round(exprnd(mean(Y),40,1);然后利用编写好的MATLAB的代码进行模拟：final.m为在16:00打磨机停止工作，function count=final(A,B) for i=1:40if(i=1)D(1)=A(1)

26、;else D(i)=A(i)+D(i-1);end end D=D;%C=A,B,D; for i=1:40if(i=1) G(1)=A(1); H(1)=G(1)+B(1);elseif(D(i)H(i-1) G(i)=H(i-1); H(i)=G(i)+B(i);else G(i)=D(i); H(i)=G(i)+B(i);end end G=G; H=H; C=A,B,D,G,H; j=0;for i=1:40 if(C(i,3)=480) j=j+1; endendcount=j; endfinal_new.m为在16:00打磨机停止工作，function count=final_n

27、ew(A,B) for i=1:40if(i=1)D(1)=A(1);else D(i)=A(i)+D(i-1);end end for i=1:40 if(D(i)240) D(i)=D(i)+60; end end D=D;%C=A,B,D; for i=1:40if(i=1) G(1)=A(1); H(1)=G(1)+B(1);elseif(D(i)H(i-1) G(i)=H(i-1); H(i)=G(i)+B(i);else G(i)=D(i); H(i)=G(i)+B(i);end end G=G; H=H; C=A,B,D,G,H; j=0;for i=1:40 if(C(i,3)=480) j=j+1; endendcount=j;end这两个函数之中的count为最后计算的可以提前下班人数。然后编写下面这样一段代码：for i=1:50A=poissrnd(12.7,40,1);B=round(exprnd(15.1000,40,1);X(i)=final(A,B);Y(i)=