Youth游乐园客流疏导方案

Youth游乐园客流疏导方案

摘 要

本文主要研究了游乐园客流疏导方案问题，通过建立TSP模型、分区域疏导 游客模型，及时为顾客提供游园线路引导；再通过时间序列分析，在多因素影响情况下对皇冠假日酒店房间预订量进行预测。

针对问题1，首先通过游客到达游乐场的时间间隔，建立服从泊松分布的人 流到达模型，将游乐场的游客量情况分为高峰期、中低峰期两种状态。然后分别建立TSP模型和M/M/s/K模型，并将这两个模型作为游乐园游客疏导模型。 该模型中我们主要考虑的是游客排队等待时间和游玩项目的数量。

针对问题2，本文通过对数据的挖掘处理及对影响房间预订量的因素分类， 建立时间数列预测模型。并运用二次指数平滑法对下一时期的房间预订量进行预测。最终利用差分公式??xt-yt，做出差分分析误差条状图，验证出时间序列预测模型的预测结果较为符合实际情况。

最后，对模型进行了评价分析与优化，并提出改进的方向。

关键字：最优路径 区域分块 时间序列预测 M/M/s/K模型

1

一、问题重述

1.1问题背景

Youth游乐园即将盛大开园，作为本市建有最多过山车的游乐园，受到了青

少年的热捧。预计届时园区将迎来每天1万的大客流。如何根据客流情况，及时分流人群，为顾客提供游园线路引导，保障游客的游园体验显得尤为重要。

1.2问题提出

（1）附件1为Youth乐园的规划图，共设A-J 共10个项目点，游客可沿着图中标出的线路往返下个游乐项目。在保障每位游客体验游乐设施的前提下，建立对每个游乐项目的等候游客进行游览提醒和疏导的模型，以达到游园体验最优。

（2）皇冠假日酒店是游乐园内的酒店，目前已开业，为有需要的游客提供住宿便利。请根据该酒店历史预订数据信息,综合考虑影响房间预定量的主要因素(比如季节,工作日/周末,法定假日,暑期等)建立数学模型。并根据酒店2015年全年预定数据(附件2),预测2016年1月至3月每天预定房间数.

二、模型假设

1、假设游客到达游乐场的时间间隔服从泊松分布；

2、假设每个游客在园内，乐意接受建议并配合相关的疏导工作；

3、假设每个游客对每个游乐项目至多体验一次，且在体验完所有项目后一定会选择离开游乐园；

4、假设不考虑游乐园内意外情况，如下大雨、设备故障等。

5、假设皇冠假日酒店是2015年1月才开业的，前三个月房间预订量相对很低是因为酒店知名度问题。

三、问题一

3.1问题分析

本问要求，在保障每位游客尽量多体验游乐设施的前提下，建立对每个游乐项目的等候游客进行游览提醒和疏导的模型，以达到游园体验最优。主要从时间方面考虑，通过建立相关模型，得出相对用时最短的路径，从而达到游客游园体验最优的目的。

根据到游客达游乐场的时间间隔服从泊松分布，分成两种情况：

第一种情况，中、低峰期（即10个游乐项目的游客数量都没有超过或刚好等于每场容纳游客数）。在中、低峰期无论游客去哪儿都不用因为排队浪费时间。这种情况下游客只需要走一条最短的路径，就可以达到游客游园体验最优（在不浪费时间的情况下体验完所有项目）的目的。因此，将此情况下的游园体验最优问题转化为TSP经典旅行商问题，再通过建立TSP模型可以求得这条最优路径。

2

第二种情况，高峰期（10个游乐项目的游客数量都超过每场容纳游客数且有一定数量的游客排队等候），此时在每一个游乐项目排队等待的游客都有两个选择：①继续排队等候；②去别的游乐项目。通过建立游客疏导模型，来给游客提供建议，从而保证游客游园体验最优的目的。

3.2建立TSP模型[1]

TSP模型是游客从单一起点出发，游玩所有的游乐项目之后，再回到原点，求解通过的最短路径。中、低峰期（在10个游乐项目的游客数量都没有超过或刚好等于每场容纳游客数时），游客可以按照TSP模型求得的这条路径到达每一个游乐项目，已达到游园体验最优（以最短的时间，最少的路程）。

游客到达过山车这一类项目，即使不用排队，如果到达的时间合适也需要等到下一场。结合附件1以及题目给出的表1. 每个游乐项目的时间安排分析可知，游客遇到过山车一类项目的等待时间均比在路上（最短的距离为250米，按照4000米/每小时计算，至少需要3.85分钟）所用的时间短，所以不考虑因为等待而改变路径的问题。假设A项目如果未达到最多容纳人数，随时去都可以玩。

设游乐项目数量为n（n?10），两项目之间的距离为dij，xij?0或1（1表

示有玩过项目i到j的路，0表示没有选择走这条路）。

则当满足：每个项目选择当前最短一条路出去，即

?xj?1nij?1,i?1,2??n

每个项目选择当前最短一条路进去，即：

?xi?1i,j?snij?1,j?1,2??n

注：除起点和中点外，各项目点不构成圈，即：

??s?1,2?s?n?1,s??1,2,??n?，s为?1,2,??n?的真子集 且：

xij??0,1?,i,j?1,2??，n，i?j

则有最短路径：

min?dijxij

i?jTSP模型的求解

利用lingo（相应的程序见附录1）对以下各式进行求解：

min?dijxij

i?j 3

?n??xij?1,i?1,2??n?j?1?n?xij?1,j?1,2??n s.t??i?1?1,2,??n?，s为?1,2,??n?的真子集???s?1,2?s?n?1,s???i,j?s?xij??0,1?,i,j?1,2??，n，i?j?TSP模型的结果分析

以项目A为起点，得到最短路程为 4350m，路径为：

A?J?I?H?G?F?D?C?B?E?A

因为出入口距离与A相距300m，所以最终最短路为4950m，路径为： 出入口?A?J?I?H?G?F?C?B?E?A?出入口 如果忽略因为到过山车一类项目因等待浪费的时间，从进入游乐园到出游乐园，所需要的最短时间为：

495010 minT???ti?2.1333(小时）4000i?1（其中，ti为每个项目每场所持续的时间）。 故为顾客提供游园线路为

出入口?A?J?I?H?G?F?C?B?E?A?出入口，以保障游客达到游园体验最优。 3.3建立分区域疏导游客模型

由于高峰期时，游客数量众多，排队时间过长会引起游客的不耐烦现象，对游乐园的经营相当不利。

对此本题通过参考快速通道模型[2]从分散客流、缩减排队时间、提高游客满意度三个方面考虑，与M/M/s/K模型[3]结合，提出了一种分区域疏导游客的且有多项目可供游客游玩的混合制模型。

在高峰期，将游乐园的工作人员分别安排在A、C、E、I、G五个点，当游客到达该点时，游客可以根据工作人员提供的信息进行区域选择游玩，从而达到疏导游客的目的，这样可以防止大量游客在某一项目大量聚集，可以减轻疏导的工作量，增加疏导效率，让游客在游乐园内的分布相对均匀。然后，将每个区域每个项目的相关数据带入M/M/s/K模型进行计算，得到游客的在相应项目的等待时间的数据，根据得到的数据判断游客在该项目是排队等待，还是离开该项目去其他项目。 3.3.1区域分块

游乐园是一个大的整体，为了提高疏散效率的目的，将游乐园分成联系紧密的几个较小的板块。观察附件1可以将游乐园分成紧密联系的4个部分，具体的分布图如图1

4

H 四 区 G F 三 区 D I 一 区 E 二 区 C J A 图1：游乐园区域分块图

B

3.3.2建立M/M/s/K模型

M/M/s/K模型是指顾客的相继到达时间在较短一段时间内服从泊松分布。

??,n?0,1,2??，K?1 ?n???0，n?K?n?,0?n?s

?s?,s?n?K其中，?：顾客的相继到达时间服从参数的负指数分布；s：项目个数；?：每个服务台服务时间相互独立的服从参数的负指数分布；K:系统的空间。

于是

??np0,0?n?s??n！ pn??n0??p,s?n?Kn?1??s!s 其中

??s?1?n?s(1??K?s?1??1ss????,??1???n！s!(1??)??n?0s?p0?? ?1ns??s?1???s?????K?s?1,??1???????n?0n！s!?n??该区域中平均滞留的总人数

Lq?Ls???n?s?pn?s

n?0s?1由于游乐园的空间是有限的，对于多个区域，顾客的有效到达率

?e???1?pK?

利用Little公式，得到

LL1Ws?s,Wq?s?Ws?

?e?e?经过对每个区域进行合理的分析，得到表1中的参数：

5

表1：各区域的参数表

参数 区域 一 二 三 四 总容纳游客数 持续时间 S Lamda ? 570 580 260 180 41.5 41.75 14.5 10 4 4 5 4 500 500 200 200 10.375 10.4375 2.9 2.5 M/M/s/K模型求解

利用Lingo软件（程序见附件2）对M/M/s/K模型求解进行求解得到结果如表2

表2：M/M/s/K模型求解结果

一区 P0 P_LOST LAMDA_E L_S L_Q W_S W_Q 0.0 0.9 41.8 579.9 575.9 13.9 13.8 二区 P0 P_LOST LAMDA_E L_S L_Q W_S W_Q 0.0 0.9 14.5 259.9 254.9 17.9 17.6 P0 P_LOST LAMDA-E L_S L_Q W_S W_Q 三区 0.0 1.0 10.0 179.9 175.9 18.0 17.6 P0 P_LOST LAMDA_E L_S L_Q W_S W_Q 四区 0.0 1.0 10.0 179.9 175.9 18.0 17.6 结果分析

对求得的结果进一步分析总结的到表3

表3：结果参数分析总结表

参数 Pn (游客能排队游区域 玩该区项目的概率) 一 二 三 四 0.917 0.9165 0.9257 0.95 Lq（该区域中平均排队人数） 565.90 575.90 254.90 175.90 Ws (在该区域中游客平均滞留的总时间) 13.73 13.89 17.92 17.90 在高峰期时：

一区，游客能排队游玩该区项目的概率为0.92，平均排队人数566，游客平均滞留（排队时间加上玩项目的时间）的总时间为13.73min；

二区，游客能排队游玩该区项目的概率为0.92，平均排队人数576，游客平均滞留（排队时间加上玩项目的时间）的总时间为13.89min；

6

三区，游客能排队游玩该区项目的概率为0.93，平均排队人数255，游客平均滞留（排队时间加上玩项目的时间）的总时间为17.92min；

四区，游客能排队游玩该区项目的概率为0.95，平均排队人数176，游客平均滞留（排队时间加上玩项目的时间）的总时间为17.90min。

游客在每个区域可排队游玩的评价概率都在0.9以上，游客在每个区域滞留的时间相对较短。所以在分区域疏导之后，游客可以按顺序游玩每个区域的项目，就可以减少排队时间和因部分项目人员过多而多夸项目在路上浪费的时间。

在高峰期，该模型可以根据客流情况，及时分流人群，为顾客提供游园线路引导，保障游客的游园体验。

因为游客在每个区域中可以顺利进行，所以游客在每个区域内玩项目的时候，游乐园的相关工作人员可以提升游客在到达一个新项目是进行排队等候。

因此该模型可以对每个游乐项目的等候游客进行游览提醒和疏导。

四、问题二

4.1问题分析

本问要求根据皇冠假日酒店历史预订数据信息,综合考虑影响房间预定量的主要因素(比如季节,工作日/周末,法定假日,暑期等)建立数学模型。并根据酒店2015年全年预定数据(附件2),预测2016年1月至3月每天预定房间数。

首先作出了全年的散点图，然后可以很清晰的观测到2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势，所以拟采用在一次指数平滑基础上加以改进的二次指数平滑时间序列来进行预测，并且由于原始数据有90个并且真实可靠，故平滑法采用的初始值以第一天的数据。然后将一次指数平滑值、二次指数平滑值、预测值的结果作在一张excel工作表格中进行对比，利用差分公式??xt-yt，做出差分分析误差条状图，进而来判断预测的效果。 4.3模型准备

step1：对附件2 中的数据，我们根据游客入住酒店的时间，按照月份为分类标准进行处理，得到以下结果（如图3）： （单位：间） 年月份房间预订量 2015 年1 月14 2015 年2 月30 2015 年3 月502 2015 年4 月4528 2015 年5 月4660 2015 年6 月5074 2015 年7 月4598 2015 年8 月4720 2015 年9 月4757 2015 年10 月4822 2015 年11 月4868 2015 年12 月4325

7

2016 年1 月320

图3

step2：时间数列影响因素分析时间数列的影响因素主要有长期趋势、季节变动、循环变动和不规则变动。

长期趋势是指受事物发展的根本因素制约而形成的事物在一段较长时期内基本趋势，可利用二次指数平滑法求解。[3]

季节变动是受自然条件（气候）、社会条件（节假日、风俗）影响的。在影响房间预订量的因素中，季节、工作日/周末、法定假日、暑假都属于季节变动因素。 循环变动具有周期长、规律性弱且不稳定的特点，因此在建立模型时，我们对循环变动因素不予考虑。

不规则变动受偶然因素和意外条件的影响，我们在进行假设时不考虑其对房间预订量的影响。

因此，在进行时间数列预测分析时，我们仅考虑长期趋势和季节变动因素。 Step3：利用matlab（程序见附件3）画出全年每天预定房间数的趋势图，如图2

图2：预定房间走势图

可以清晰的观测到2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势。所以可以直接用时间序列预测模型结合前三个月的数据，对2016年前三个月每天的预订房间数量进行预测。 4.3建立时间序列预测[4]模型

时间序列预测是以时间数列所能反映的社会经济现象的发展过程和规律性，进行引伸外推，预测其发展趋势的方法。

从预定房间走势图来看2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势，所以拟采用在一次指数平滑基础上加以改进的二次指数平滑时间序列来进行预测，并且由于原始数据有90个并且真实可靠，故平滑法采用的初始值以第一天的数据。 4.4模型求解与结果分析

8

利用matlab（程序见附录4）进行运算求得2016年前三个月的预测值趋势图见图3

图3：2016年前三个月的预测值趋势图

预测的部分具体数据见表4，完整的数据见附录5

表4：2016年一月到三月理论上每天预定房间数量的预测表

日期 1-1 1-2 1-3 1-4 1-5 1-6 1-7 1-8 1-9 1-10 1-11 1-12 1-13 1-14 1-15 1-16 1-17 1-18 1-19 房间数 0.0000 2.0000 1.7000 1.6400 2.3850 2.0146 1.6944 1.8181 2.1202 1.7736 1.8752 1.5587 1.2867 1.8536 1.5342 1.6602 4.1656 3.4933 2.9139 日期 1-24 1-25 1-26 1-27 1-28 1-29 1-30 1-31 2-1 2-2 2-3 2-4 2-5 2-6 2-7 2-8 2-9 2-10 2-11 房间数 1.3473 1.0706 0.8358 0.6372 0.4700 0.3298 0.2130 0.1162 0.0367 0.1719 0.0897 0.2222 0.1373 0.2672 0.1797 0.1071 0.4472 0.3382 0.2465 日期 2-16 2-17 2-18 2-19 2-20 2-21 2-22 2-23 2-24 2-25 2-26 2-27 2-28 2-29 3-1 3-2 3-3 3-4 3-5 房间数 -0.0253 -0.0535 -0.0758 -0.0931 -0.1062 -0.1158 -0.1224 -0.1265 -0.1285 0.2711 0.8121 0.8721 0.7221 1.1933 2.1931 2.8409 2.7873 2.7360 3.2871 日期 3-9 3-10 3-11 3-12 3-13 3-14 3-15 3-16 3-17 3-18 3-19 3-20 3-21 3-22 3-23 3-24 3-25 3-26 3-27 房间数 3.6682 6.0454 7.4604 6.4520 6.1806 5.1390 4.4439 7.0465 8.0541 8.0997 9.3256 8.5554 9.2859 9.8948 11.9991 17.3735 20.3231 21.7990 23.0165 由于预定的房间数量为整数，但是在预测中出现了小数，所以需要对数据进行处理。原则上只有有0.001个房间也要写为一个房间，但是对于房间数量预测出的数值，一方面反映了预定房间的数量；另一方面反映了预定一个的概率。如

9

果预定房间的概率小于0.5则说明预定房间的可能性不大，所以对于该数据处理可以采用四舍五入法。2016年一月到三月实际上每天预定房间数量的预测如表5。

表5：2016年一月到三月实际每天预定房间数量的预测表

日期 1-2 1-3 1-4 1-5 1-6 1-7 1-8 1-9 1-10 1-11 1-12 1-13 1-14 1-15 1-16 1-17 1-18 1-19 1-20 1-21 1-22 1-23 1-24 房间 0 2 2 2 2 2 2 2 2 2 2 2 1 2 2 2 4 3 3 3 2 2 2 日期 1-25 1-26 1-27 1-28 1-29 1-30 1-31 2-1 2-2 2-3 2-4 2-5 2-6 2-7 2-8 2-9 2-10 2-11 2-12 2-13 2-14 2-15 2-16 房间 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 日期 2-17 2-18 2-19 2-20 2-21 2-22 2-23 2-24 2-25 2-26 2-27 2-28 2-29 3-1 3-2 3-3 3-4 3-5 3-6 3-7 3-8 3-9 房间 0 0 0 0 0 0 0 0 0 0 1 1 1 1 2 3 3 3 3 5 4 4 日期 3-10 3-11 3-12 3-13 3-14 3-15 3-16 3-17 3-18 3-19 3-20 3-21 3-22 3-23 3-24 3-25 3-26 3-27 3-28 3-29 3-30 3-31 房间 4 6 7 6 6 5 4 7 8 8 9 9 9 10 12 17 20 22 23 22 21 20 五、模型检验

10

问题二，时间序列预测模型的检验：

对问题二预测的结果进行差分分析（matlab程序见附件4）具体的分析图见图4

图4：差分分析图

由差分分析误差条状图可以知道，预测值和去年的实际值呈现一阶差分趋势，表明时间序列预测模型的预测结果是符合理想。

六、模型评价与推广

6.1模型优点

问题一，1、考虑了人流到达的不规律性，将人流到达假设为服从泊松分布，进而考虑了低峰期和高峰期两种情况下的疏导模型，低峰期为游客规划了一条最短路径；2、高峰期将游乐场合理的划分为四个区域，进行建议式的疏导，这样既保证了为游客提供了最优游览方案，又让游客拥有自主选择权，互惠互利。 问题二，1、简单易行，便于掌握，能够充分运用原时间序列的各项数据；2、计算速度快，对模型参数有动态确定的能力，精度较好。 6.2模型缺点

问题一，没有考虑一些游客喜欢按照自己的方案游玩的特殊情况，让该疏导模型出现拥堵的情况增加了一定的不确定性。

问题二，不能反映事物的内在联系，不能分析两个因素的相关关系，在处理问题时可能存在一定误差。 6.3模型改进

问题一的方法相当于在为游客提供建议方案，让游客自由选择游玩路径，并不一定能达到我们预期想要的结果。因此需要寻求更优的算法对模型进行求解，例如：利用计算机仿真模拟等方法对模型进行改进。

问题二采用时间序列二次指数平滑法来进行预测，只是这种预测的趋势呈现一种平稳、线性的形式，可能会对结果产生较大的误差。需要寻求更精确的方法进行预测。而在预测类方法中有时间序列预测、灰色预测、神经元网络、差分方程等，经过深层次分析，利用时间序列预测模型和灰色预测中的GM?1,1?模型分

11

别进行预测，最后使两种预测进行加权平均法得出一个组合预测模型模型来进行预测，使得出的结果更精确。 6.4模型推广

TSP模型运用广泛，可以用于最短路径类问题的求解；分区域疏导游客模型可以解决多服务地点高峰期的疏导类问题；时间数列预测模型，可以用于预测公司收入等问题。

参考文献

[1] 司守奎,孙兆亮.数学建模算法与应用.北京.国防工业出版社.2015年4月.58-61.

[2] 陈治佳,王曦,何苗.大型游乐场快速通道优化模型与仿真模拟[J].哈尔滨工业大学学报.哈尔滨.第39卷第7期.2005年09月.101-103.

[3] 丛国超,朱翼隽.批量到达的多服务台排队模型求解[J].成都信息工程学院学报.成都.第22卷第1期.2007年01月.98-100.

[4] 司守奎,孙兆亮.数学建模算法与应用.北京.国防工业出版社.2015年4月.170-173.

12

附 录

附录1 运行环境：lingo11 MODEL: SETS: Entertainment / 1.. 10/: U;! U( I) = sequence no. of Entertaiment; LINK( Entertainment , Entertainment ): DIST, ! DIST matrix, it need not be symmetric; X; ! X( I, J) = 1 if we use link I, J; ENDSETS DATA: DIST =0 300 600 1050 350 1550 900 1050 600 250 300 0 300 750 650 1250 1200 1350 900 550 600 300 0 450 500 950 1050 1400 950 850 1050 750 450 0 950 500 1150 1550 1400 1300 350 650 500 950 0 1200 550 900 450 600 1550 1250 950 500 1200 0 650 1050 1500 1800 900 1200 1050 1150 550 650 0 400 850 1150 1050 1350 1400 1550 900 1050 400 0 450 800 600 900 950 1400 450 1500 850 450 0 350 250 550 850 1300 600 1800 1150 800 350 0 ; ENDDATA N = @SIZE( Entertainment ); MIN = @SUM( LINK: DIST * X); @FOR( Entertainment ( K): @SUM( Entertainment ( I)| I #NE# K: X( I, K)) = 1; @SUM( Entertainment ( J)| J #NE# K: X( K, J)) = 1; ! Weak form of the subtour breaking constraints; ! These are not very powerful for large problems; @FOR( Entertainment ( J)| J #GT# 1 #AND# J #NE# K: U( J) >= U( K) + X ( K, J) - ( N - 2) * ( 1 - X( K, J)) + ( N - 3) * X( J, K)) ); ! Make the X's 0/1; @FOR( LINK: @BIN( X)); ! For the first and last stop we know...; @FOR( Entertainment ( K)| K #GT# 1: U( K) <= N - 1 - ( N - 2) * X( 1, K); U( K) >= 1 + ( N - 2) * X( K, 1)); END

Lingo11 附录2 运行环境：

13

model: sets: state/1..570/:p; endsets lamda=500;mu=10.375;rho=lamda/mu;s=4;k=570; lamda*p0=mu*p(1); (lamda+mu)*p(1)=lamda*p0+2*mu*p(2); @for(state(i)|i #gt#1 #and# i #lt# s: (lamda+i*mu)*p(i)=lamda*p(i-1)+(i+1)*mu*p(i+1)); @for(state(i)|i #ge# s #and# i #lt# k: (lamda+s*mu)*p(i)=lamda*p(i-1)+s*mu*p(i+1)); lamda*p(k-1)=s*mu*p(k); p0+@sum(state:p)=1; P_lost=p(k);lamda_e=lamda*(1-P_lost); L_s=@sum(state(i):i*p(i)); L_q=L_s-lamda_e/mu; W_s=L_s/lamda_e; W_q=W_s-1/mu; end mode2: sets: state/1..580/:p; endsets lamda=500;mu=10.4375;rho=lamda/mu;s=4;k=580; lamda*p0=mu*p(1); (lamda+mu)*p(1)=lamda*p0+2*mu*p(2); @for(state(i)|i #gt#1 #and# i #lt# s: (lamda+i*mu)*p(i)=lamda*p(i-1)+(i+1)*mu*p(i+1)); @for(state(i)|i #ge# s #and# i #lt# k: (lamda+s*mu)*p(i)=lamda*p(i-1)+s*mu*p(i+1)); lamda*p(k-1)=s*mu*p(k); p0+@sum(state:p)=1; P_lost=p(k);lamda_e=lamda*(1-P_lost); L_s=@sum(state(i):i*p(i)); L_q=L_s-lamda_e/mu; W_s=L_s/lamda_e; W_q=W_s-1/mu; end mode3: sets: state/1..260/:p; endsets lamda=200;mu=2.9;rho=lamda/mu;s=5;k=260; 14

lamda*p0=mu*p(1); (lamda+mu)*p(1)=lamda*p0+2*mu*p(2); @for(state(i)|i #gt#1 #and# i #lt# s: (lamda+i*mu)*p(i)=lamda*p(i-1)+(i+1)*mu*p(i+1)); @for(state(i)|i #ge# s #and# i #lt# k: (lamda+s*mu)*p(i)=lamda*p(i-1)+s*mu*p(i+1)); lamda*p(k-1)=s*mu*p(k); p0+@sum(state:p)=1; P_lost=p(k);lamda_e=lamda*(1-P_lost); L_s=@sum(state(i):i*p(i)); L_q=L_s-lamda_e/mu; W_s=L_s/lamda_e; W_q=W_s-1/mu; End Mode4: sets: state/1..180/:p; endsets lamda=200;mu=2.5;rho=lamda/mu;s=4;k=180; lamda*p0=mu*p(1); (lamda+mu)*p(1)=lamda*p0+2*mu*p(2); @for(state(i)|i #gt#1 #and# i #lt# s: (lamda+i*mu)*p(i)=lamda*p(i-1)+(i+1)*mu*p(i+1)); @for(state(i)|i #ge# s #and# i #lt# k: (lamda+s*mu)*p(i)=lamda*p(i-1)+s*mu*p(i+1)); lamda*p(k-1)=s*mu*p(k); p0+@sum(state:p)=1; P_lost=p(k);lamda_e=lamda*(1-P_lost); L_s=@sum(state(i):i*p(i)); L_q=L_s-lamda_e/mu; W_s=L_s/lamda_e; W_q=W_s-1/mu; End

附录3 运行环境：matlaba=xlsread('酒店预定历史数据2015年.xls'); %将“附件2“更名为“酒店预定历史数据2015年.xls”，并保存在matlab子文件夹下 x=a(:,1);y=a(:,2); plot(x,y,'k');grid on; xlabel('距1月1日的天数');ylabel('预订房间数'); title('2015年时间-预订房间数'); r2007a

15

附录4 运行环境：matlab时间序列预测程序：matlab R2007a clc clear yt=load('酒店预定历史数据2015年.txt'); r2007a %将预订房间数数据以列的形式保存在“酒店预定历史数据2015年.txt”中，并保存在matlab子文件夹下 n=length(yt),alpha=0.1; st1(1)=yt(1);st2(1)=yt(1); for i=2:n st1(i)=alpha*yt(i)+(1-alpha)*st1(i-1); st2(i)=alpha*st1(i)+(1-alpha)*st2(i-1); end xlswrite('酒店预定历史数据2015年.xls ',[st1',st2'],'Sheet1','C2:D366') %将预订房间数一次指数平滑值、二次指数平滑值写在在“酒店预定历史数据2015年.xls”中的C、D列 at=2*st1-st2; bt=alpha/(1-alpha)*(st1-st2); yhat=at+bt; xlswrite('酒店预定历史数据2015年.xls',yhat','Sheet1','E2') str=['E',int2str(n+2)]; xlswrite('酒店预定历史数据2015年.xls',at(n)+2*bt(n),'Sheet1',str) %将预订房间数的预测值写在在“酒店预定历史数据2015年.xls”中的E列 xt=xlsread('酒店预定历史数据2015年.xls',1,'A2:E91'); %将“酒店预定历史数据2015年.xls”中[A2:E91]的数据读取入xt矩阵中 t=xt(:,1);y1=xt(:,2);y2=xt(:,3);y3=xt(:,4); y4=xt(:,5);y5=xt(:,5)-xt(:,2); plot(t,y1,'o');grid on;hold on; plot(t,y2,'b-');grid on;hold on; plot(t,y3,'r--');grid on;hold on; plot(t,y4,'k*');grid on; legend('原始数据','一次指数平滑值','二次指数平滑值','预测值',0); xlabel('距1月1日的天数');ylabel('预订房间数'); title('时间-预订房间数');% 并将作出的图形进行灰度处理 plot(t,y5,'o');grid on;hold on; xlabel('当日预测值-当日实际值');ylabel('误差房间数量'); title('差分分析');%并将作出的图形进行灰度处理和线性改为条状

附录5 日期 1-1 1-2

房间数 0.0000 2.0000 日期 1-24 1-25 房间数 1.3473 1.0706 日期 2-16 2-17 16

房间数 -0.0253 -0.0535 日期 3-9 3-10 房间数 3.6682 6.0454

1-3 1-4 1-5 1-6 1-7 1-8 1-9 1-10 1-11 1-12 1-13 1-14 1-15 1-16 1-17 1-18 1-19 1-20 1-21 1-22 1-23 1.7000 1.6400 2.3850 2.0146 1.6944 1.8181 2.1202 1.7736 1.8752 1.5587 1.2867 1.8536 1.5342 1.6602 4.1656 3.4933 2.9139 3.0154 2.4974 2.0529 1.6723 1-26 1-27 1-28 1-29 1-30 1-31 2-1 2-2 2-3 2-4 2-5 2-6 2-7 2-8 2-9 2-10 2-11 2-12 2-13 2-14 2-15 0.8358 0.6372 0.4700 0.3298 0.2130 0.1162 0.0367 0.1719 0.0897 0.2222 0.1373 0.2672 0.1797 0.1071 0.4472 0.3382 0.2465 0.1698 0.1059 0.0532 0.0099 2-18 2-19 2-20 2-21 2-22 2-23 2-24 2-25 2-26 2-27 2-28 2-29 3-1 3-2 3-3 3-4 3-5 3-6 3-7 3-8 -0.0758 -0.0931 -0.1062 -0.1158 -0.1224 -0.1265 -0.1285 0.2711 0.8121 0.8721 0.7221 1.1933 2.1931 2.8409 2.7873 2.7360 3.2871 4.5506 4.0185 4.1573 3-11 3-12 3-13 3-14 3-15 3-16 3-17 3-18 3-19 3-20 3-21 3-22 3-23 3-24 3-25 3-26 3-27 3-28 3-29 3-30 7.4604 6.4520 6.1806 5.1390 4.4439 7.0465 8.0541 8.0997 9.3256 8.5554 9.2859 9.8948 11.9991 17.3735 20.3231 21.7990 23.0165 22.2125 21.4892 20.2384

17

本文来源：https://www.bwwdw.com/article/b0tt.html