数学建模.docx

1、数学建模海南师范大学大学生数学建模竞赛承 诺 书我们仔细阅读了海南师范大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B中选择一项填写）： B 我们的参赛报名号为（如果设置报名号的

2、话）： 参赛队员 (打印并签名) ：1. 夏林 2. 孔倩 3. 张丽琛 指导教师或指导教师组负责人 (打印并签名)： 日期： 2013 年 4 月 30一、摘要2二、问题的重述3三、问题的分析3四、建模的过程4 1.模型假设4 2.符号定义说明4 3.模型建立5 4.模型的求解与分析9 5.模型检验13五、模型评价13六、参考文献13七、附录14大学食堂窗口开设问题一、 摘要大学校园学生数目日益剧增，食堂开设的打饭窗口供不应求，而开设新窗口的成本愈来愈高，如何实现食堂窗口数目的合理化，从而使大学生排队打饭的效率得以提升，同时为食堂带来最大的经济收益问题，是本文主要探究的对象。对于问题一，我们

3、小组从题中给出数据出发，对表中数据进行分析，根据排队论知识，运用matlab软件绘画出食堂人流的图形，进而验证得出食堂人流模型服从泊松分布。对于问题二，我们小组从学生满意度、食堂运营成本、食堂的服务效率三个方面出发，利用层次分析法，确立食堂服务综合评价模型。进而通过平均等待时间、排队等待人数、服务强度等因素的判断出食堂开设4个窗口是合理的，其综合评价合理度为I=0.82。对于问题三，我们小组不考虑一天中上午、下午开设窗口数可能不同的情况。由数据分析得出各个工作日不同时间段服务的参数值分布情况，结合排队理论知识与问题二的综合评价模型，运用线性规划的方法，得出窗口个数为5时，其合理度最大。关键词：

4、排队论、食堂综合评价模型、泊松分布、层次分析法二、 问题的重述食堂开设新窗口需要一定的成本，快捷的就餐方式是大学生高效率生活的保障（服务每位同学所需时间一般为35s55s），根据提供的资料，回答下列问题：某大学食堂一周内的学生打饭人数统计表：周一周二周三周四周五周六周日中午11:00至11:157354597581845411:15至11:307675806992955911:30至11:458894968998978911:45至12:00120951391567312011612:00至12:1540535666515749晚上17:00至17:152840355462355117:15至

5、17:303345595883575117:30至17:455981719970787017:45至18:007583798065827418:00至18:1568636251586557（1）描述食堂人流的模型；（2）假设饭堂共有八个买饭窗口，如果时常开四个窗口，请分析其合理与否？（3）建立模型描述该食堂开设多少窗口打饭能够满足学生的需要，增加食堂的收益？三、 问题的分析本题的目的在于通过对学校食堂排队系统进行建模分析，使得食堂能够良好的运营。我们通过对食堂某个星期内全部工作日时段排队打饭的学生情况进行统计，运用数学软件MATLAB编程对收集到的数据进行分析，拟合出整个系统中数据呈现的规律或

6、概率以及学生等待时间、学生排队队列长、食堂为每位学生打饭的时间等随机事件的规律或概率。再根据实际数据代入数学模型计算得出相应数值，反映出食堂服务系统的服务效率，得出服务率最高的时候目标函数的值，从而解决问题。四、建模过程4.1模型假设 4.1.1食堂开放时间固定不变，忽略放假，装修等变动； 4.1.2学生来食堂就餐是相互独立的，互不干扰； 4.1.3学生自动选择人少的窗口，每人只打一份饭，且对饭菜胃口无特殊偏好； 4.1.4食堂饭菜供给充足，且每个窗口的服务质量和数量相同； 4.1.5 食堂一天内开设窗口数量是一致的，不考虑上午、下午开设窗口数量不一 致的情况； 4.1.6 不考虑学生到达食堂

7、后的流失人数。4.2符号定义说明：表示系统中的学生数，包括排队等候的和正在接受服务的所学生（称为平均队列长队）；: 表示系统中排队等候的学生数（称为平均队列长）；: 表示学生在系统中的平均逗留时间(包括等待时间和服务时间)； : 表示学生在系统中的平均等待时间（平均排队等待时间）；: 表示排成一大队列时的平均等待时间；: 表示排成一大队列时的平均队列长；: 表示排成k个小队时的平均等待时间；f(k):权重组合函数；: 表示排成k个小队时的平均队列长；：表示排队的平均到达率(称为学生到达速率)；: 表示系统的平均服务率(即食堂服务台的平均服务速率)；:窗口数量 ；：权重（i=1,2,3）;：平均

8、每日排队打饭学生到达人数；：周一至周五平均每日各时段排队学生到达人数；：周六周日平均每日各时段排队学生到达人数；：每小时到达人数；：排队流失学生人数；：窗口完全空闲的概率；：系统中有n个学生的概率；: 表示服务强度，其值为有效的平均到达率与平均服务率之比，即其中主要性指标是，。主要性指标其值越小，说明系统排队越少，等待时间越少，因而系统性能越好。显然，它们是打饭排队学生与食堂服务部门都很关注的因素。排队学生希望等待时间和队列长越短越好，当然对食堂服务人员来说，服务强度越小越好。4.3模型的建立4.3.1 问题一 模型的建立排队理论系统说明：所谓M/M/k的排队系统是指这样的一种服务：当天打饭排

9、队学生的到达服从参数为的泊松分布；排队学生接受服务时间服从参数为的指数分布；有k个服务台（食堂服务窗口），打饭排队学生按到达的先后次序接受服务。泊松分布： (为常数，k=0,1,2,)即在时间T内有k位排队学生的到达的概率为： 其中是在时间内学生到达的平均学生数，平均到达率。4.3.2问题二模型建立服务时间服从负指数分布：其中为大于0的常数，代表单位时间内的平均服务率。设在任意时刻t系统中有n个排队学生的概率为。当系统达到稳定状态后，趋于稳定状态概率此时，与t无关，称系统处于统计平衡状态，并称为统计平衡状态下的稳态概率，它表示系统在稳定状态下有n个学生的概率，此时，特别，表示稳态系统所有服务台

10、全部空闲的概率。其中： ；服务强度：；平均队长：；平均队列长：；平均逗留时间：；平均等待时间：图1如图1所示，此时问题归结为一个M/M/k/排队系统，即排成一个大队对k个窗口的情况。根据排队理论，当服务强度1；(2) 当开设窗口数k=2时: =1.81131；(3) 当开设窗口数k=3时:=1.20751； 综合（1）（2）（3）可知，当k=1或k=2或k=3时，服务强度大于1，即系统内排队学生的到达率大于系统的平均服务率，可见系统不存在平衡状态，且排队的人会越来越多，排队等候的时间也会越来越长，按照题目所给每天的工作时间为2.5个小时，k=1时，会有个学生无法完成排队打饭，k=2时会有个学生

11、无法完成排队打饭，k=3会有个学生无法完成排队打饭，这样不仅会造成大量的排队打饭学生流失，而且可能对学生吃饭有很大影响，因此学校食堂开设1个2个或者3个窗口均不满足要求，需要增开窗口才能满足排队打饭学生的需求。(4) 当开设窗口数k=4时: =0.90561, 服务强度小于1,系统可以达到平衡状态。此时：平均排队长度：=26.34平均等待时间：=5.45分钟其食堂优化合理度为I=0.82，故食堂开设4个窗口是合理的。4.4.3问题三的求解 以问题二所建立的综合评价食堂优化合理度模型为基础，经分析可知：当食堂窗口开放数量k3时，服务强度，即显示学生到达食堂率大于食堂服务学生的速率，也就表明食堂的

12、服务存在着不平衡状态，学生前来排队的人数会越来越多，排队等待时间则会延长。这时需要增加窗口数量。当k=4时，由问题二可知，此时窗口的服务强度为0.911, T=5.45, =26.34, 此时，需要排5.45分钟，其食堂优化合理度为I=0.82当开设窗口数k=5时:=0.7245f(5)食堂开设窗口数为5个时，学生平均等待时间较短，不足3分钟，排队人数是k=4时的一半，所以基本不会存在排队现象。K=5时对于排队学生来说，满意率更为提高， 其合理度I=0.833.而当k=6,7,8时，学生平均等待的时间更短，排队人数更少。学生的满意度就更高了，但对于食堂而言，随着窗口数的增加，运营成本增加，而服

13、务效率和学生满意度的增加速度小于运营成本的增加数，其合理度I A=68.57 78.33 93.80 116.25 52.33 50.50 55.14 75.43 76.00 58.60alpha=0.05;lamda=poissfit(A,alpha);p3=poisscdf(A,lamda);H3,s3=kstest(A,A,p3,alpha)n=length(A);if H3=0disp(该数据服从泊松分布) elsedisp(该数据不服从泊松分布) endA = 68.5700 78.3300 93.8000 116.2500 52.3300 50.5000 55.1400 75.4300 76.0000 58.6000H3 = 0s3 = 0.1267该数据服从泊松分布附录三：人流到达模型图 x=0:16;y1=poisspdf(x,4.83);hold onplot(x,y1,:r*); hold off;title(食堂人流分布)附录四：学生满意度、运营成本、服务效率权重代码 A=1 1/2 1/2;2 1 1;2 1 1 x,y=eig(A)%求得x为特征向量矩阵 m m=find(y=max(max(y)%找到y中对应最大的特征值所在列m w=x(:,m)/sum(x(:,m) %w即为矩阵A的权重！