Youth游乐园客流疏导方案

17Youth游乐园客流疏导方案摘要本文主要研究了游乐园客流疏导方案问题，通过建立TSP模型、分区域疏导游客模型，及时为顾客提供游园线路引导；再通过时间序列分析，在多因素影响情况下对皇冠假日酒店房间预订量进行预测。针对问题1，首先通过游客到达游乐场的时间间隔，建立服从泊松分布的人流到达模型，将游乐场的游客量情况分为高峰期、中低峰期两种状态。然后分别建立TSP模型和M/M/s/K模型，并将这两个模型作为游乐园游客疏导模型。该模型中我们主要考虑的是游客排队等待时间和游玩项目的数量。针对问题2，本文通过对数据的挖掘处理及对影响房间预订量的因素分类，建立时间数列预测模型。并运用二次指数平滑法对下一时期的房间预订量进行预测。最终利用差分公式A=xt-yt，做出差分分析误差条状图，验证出时间序列预测模型的预测结果较为符合实际情况。最后，对模型进行了评价分析与优化，并提出改进的方向。关键字：最优路径区域分块时间序列预测M/M/s/K模型一、问题重述问题背景Youth游乐园即将盛大开园，作为本市建有最多过山车的游乐园，受到了青少年的热捧。预计届时园区将迎来每天1万的大客流。如何根据客流情况，及时分流人群，为顾客提供游园线路引导，保障游客的游园体验显得尤为重要。问题提出（1）附件1为Youth乐园的规划图，共设A-J共10个项目点，游客可沿着图中标出的线路往返下个游乐项目。在保障每位游客体验游乐设施的前提下，建立对每个游乐项目的等候游客进行游览提醒和疏导的模型，以达到游园体验最优。（2）皇冠假日酒店是游乐园内的酒店，目前已开业，为有需要的游客提供住宿便利。请根据该酒店历史预订数据信息,综合考虑影响房间预定量的主要因素（比如季节,工作日/周末,法定假日,暑期等）建立数学模型。并根据酒店2015年全年预定数据（附件2）,预测2016年1月至3月每天预定房间数.二、模型假设1、假设游客到达游乐场的时间间隔服从泊松分布；2、假设每个游客在园内，乐意接受建议并配合相关的疏导工作；3、假设每个游客对每个游乐项目至多体验一次，且在体验完所有项目后一定会选择离开游乐园；4、假设不考虑游乐园内意外情况，如下大雨、设备故障等。5、假设皇冠假日酒店是2015年1月才开业的，前三个月房间预订量相对很低是因为酒店知名度问题。三、问题一问题分析本问要求，在保障每位游客尽量多体验游乐设施的前提下，建立对每个游乐项目的等候游客进行游览提醒和疏导的模型，以达到游园体验最优。主要从时间方面考虑，通过建立相关模型，得出相对用时最短的路径，从而达到游客游园体验最优的目的。根据到游客达游乐场的时间间隔服从泊松分布，分成两种情况：第一种情况，中、低峰期（即10个游乐项目的游客数量都没有超过或刚好等于每场容纳游客数）。在中、低峰期无论游客去哪儿都不用因为排队浪费时间。这种情况下游客只需要走一条最短的路径，就可以达到游客游园体验最优（在不浪费时间的情况下体验完所有项目）的目的。因此，将此情况下的游园体验最优问题转化为TSP经典旅行商问题，再通过建立TSP模型可以求得这条最优路径。第二种情况，高峰期（10个游乐项目的游客数量都超过每场容纳游客数且有一定数量的游客排队等候），此时在每一个游乐项目排队等待的游客都有两个选择：①继续排队等候；②去别的游乐项目。通过建立游客疏导模型，来给游客提供建议，从而保证游客游园体验最优的目的。游髀进入游乐园的排队过程服从泊松分布［茁，艮分布列函数为：瓦屮久为单位吋间内顾客流入屋的期望值。假设该游乐园开放吋间为，9:00-21:00,--天中游客容纳量为10000人,则有七A=JOWO_=13112x603.2建立TSP模型［i］TSP模型是游客从单一起点出发，游玩所有的游乐项目之后，再回到原点，求解通过的最短路径。中、低峰期（在10个游乐项目的游客数量都没有超过或刚好等于每场容纳游客数时），游客可以按照TSP模型求得的这条路径到达每一个游乐项目，已达到游园体验最优（以最短的时间，最少的路程）。游客到达过山车这一类项目，即使不用排队，如果到达的时间合适也需要等到下一场。结合附件1以及题目给出的表1.每个游乐项目的时间安排分析可知，游客遇到过山车一类项目的等待时间均比在路上（最短的距离为250米，按照4000米/每小时计算，至少需要3.85分钟）所用的时间短，所以不考虑因为等待而改变路径的问题。假设A项目如果未达到最多容纳人数，随时去都可以玩。设游乐项目数量为n（n二10），两项目之间的距离为d，x二0或1（1表ijij示有玩过项目i到j的路，0表示没有选择走这条路）。则当满足：每个项目选择当前最短一条路出去，即Yx=1,i=1,2nijj=1每个项目选择当前最短一条路进去，即：n工x=1,j=1,2ni=1注：除起点和中点外，各项目点不构成圈，即：工<|S-1,2<|S<n-1,suh,2,n}s为£,2,n}的真子集i,jes且：xe{o,l},i,j=1,2，n，i丰jij则有最短路径：ijminYdxijijTSP模型的求解i丰j利用lingo（相应的程序见附录1）对以下各式进行求解:minYdxijiji丰jx—1,i—1,2nijj—1刀x—1,j—1,2nijn},s为{1,2,n摘真子集<|s|-1,2<|s|<n},s为{1,2,n摘真子集i,jesr[xg{0,1J,i,j—1,2,n,i丰jijTSP模型的结果分析以项目A为起点，得到最短路程为4350m，路径为：ATJTITHTGTFTDTCTBTETA因为出入口距离与A相距300m，所以最终最短路为4950m，路径为：出入口TATJTITHTGTFTCTBTETAT出入口如果忽略因为到过山车一类项目因等待浪费的时间，从进入游乐园到出游乐园，所需要的最短时间为：minT—4950+兰t—2.1333（小时）4000j（其中，t为每个项目每场所持续的时间）。i故为顾客提供游园线路为出入口TATJTITHTGTFTCTBTETAT出入口，以保障游客达到游园体验最优。建立分区域疏导游客模型由于高峰期时，游客数量众多，排队时间过长会引起游客的不耐烦现象，对游乐园的经营相当不利。对此本题通过参考快速通道模型［2］从分散客流、缩减排队时间、提高游客满意度三个方面考虑，与M/M/s/K模型⑶结合，提出了一种分区域疏导游客的且有多项目可供游客游玩的混合制模型。在高峰期，将游乐园的工作人员分别安排在A、C、E、I、G五个点，当游客到达该点时，游客可以根据工作人员提供的信息进行区域选择游玩，从而达到疏导游客的目的，这样可以防止大量游客在某一项目大量聚集，可以减轻疏导的工作量，增加疏导效率，让游客在游乐园内的分布相对均匀。然后，将每个区域每个项目的相关数据带入M/M/s/K模型进行计算，得到游客的在相应项目的等待时间的数据，根据得到的数据判断游客在该项目是排队等待，还是离开该项目去其他项目。3.3.1区域分块游乐园是一个大的整体，为了提高疏散效率的目的，将游乐园分成联系紧密的几个较小的板块。观察附件1可以将游乐园分成紧密联系的4个部分，具体的分布图如图1

3.3.2建立3.3.2建立M/M/s/K模型M/M/s/K模型是指顾客的相继到达时间在较短一段时间内服从泊松分布。九=严,n=0,1,2n|0,n>Kn卩,0<n<ss卩,s<n<K其中，九：顾客的相继到达时间服从参数的负指数分布；s:项目个数；卩：每个服务台服务时间相互独立的服从参数的负指数分布；K:系统的空间。PnPnpo,0<n<sn!n！Pnp0,s<n<K、s!sn-1其中p=]PnPsp=]PnPs(1—pK-s+1”Is0n！s!(1—p)丿n=0迟巴+匕(K—s+1js!丿-1,Ps丰1n！n=0-1,Ps=1一*一*(n-s)p-snn=0L=L

qn=0由于游乐园的空间是有限的，对于多个区域，顾客的有效到达率九=X(1—pK)e利用Little公式，得到LL1W=冇,W二寸=W-—sXqXsRee经过对每个区域进行合理的分析，得到表1中的参数:

表1：各区域的参数表7、..、参数区域总容纳游客数持续时间SLamda.57041.5450010.375二58041.75450010.4375三26014.552002.9四1801042002.5M/M/s/K模型求解利用Lingo软件（程序见附件2）对M/M/s/K模型求解进行求解得到结果如表2表2：M/M/s/K模型求解结果一区二区三区四区P00.0P00.0P00.0P00.0PLOST0.9PLOST0.9PLOST1.0PLOST1.0LAMDAE41.8LAMDAE14.5LAMDA-E10.0LAMDAE10.0LS579.9LS259.9LS179.9LS179.9LQ575.9LQ254.9LQ175.9LQ175.9WS13.9WS17.9WS18.0WS18.0WQ13.8WQ17.6WQ17.6WQ17.6结果分析对求得的结果进一步分析总结的到表3表3：结果参数分析总结表参数区域Pn（游客能排队游玩该区项目的概率）Lq（该区域中平均排队人数）Ws（在该区域中游客平均滞留的总时间）-一一0.917565.9013.73二0.9165575.9013.89三0.9257254.9017.92四0.95175.9017.90在高峰期时：一区，游客能排队游玩该区项目的概率为0.92，平均排队人数566，游客平均滞留（排队时间加上玩项目的时间）的总时间为13.73min；二区，游客能排队游玩该区项目的概率为0.92，平均排队人数576，游客平均滞留（排队时间加上玩项目的时间）的总时间为13.89min；三区，游客能排队游玩该区项目的概率为0.93，平均排队人数255，游客平均滞留（排队时间加上玩项目的时间）的总时间为17.92min；四区，游客能排队游玩该区项目的概率为0.95，平均排队人数176，游客平均滞留（排队时间加上玩项目的时间）的总时间为17.90min。游客在每个区域可排队游玩的评价概率都在0.9以上，游客在每个区域滞留的时间相对较短。所以在分区域疏导之后，游客可以按顺序游玩每个区域的项目，就可以减少排队时间和因部分项目人员过多而多夸项目在路上浪费的时间。在高峰期，该模型可以根据客流情况，及时分流人群，为顾客提供游园线路引导，保障游客的游园体验。因为游客在每个区域中可以顺利进行，所以游客在每个区域内玩项目的时候，游乐园的相关工作人员可以提升游客在到达一个新项目是进行排队等候。因此该模型可以对每个游乐项目的等候游客进行游览提醒和疏导。四、问题二问题分析本问要求根据皇冠假日酒店历史预订数据信息,综合考虑影响房间预定量的主要因素（比如季节,工作日/周末,法定假日,暑期等）建立数学模型。并根据酒店2015年全年预定数据（附件2）,预测2016年1月至3月每天预定房间数。首先作出了全年的散点图，然后可以很清晰的观测到2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势，所以拟采用在一次指数平滑基础上加以改进的二次指数平滑时间序列来进行预测，并且由于原始数据有90个并且真实可靠，故平滑法采用的初始值以第一天的数据。然后将一次指数平滑值、二次指数平滑值、预测值的结果作在一张excel工作表格中进行对比，利用差分公式A=xt-yt，做出差分分析误差条状图，进而来判断预测的效果。模型准备stepl:对附件2中的数据，我们根据游客入住酒店的时间，按照月份为分类标准进行处理，得到以下结果（如图3）：（单位：间）年月份房间预订量2015年1月142015年2月302015年3月5022015年4月45282015年5月46602015年6月50742015年7月45982015年8月47202015年9月47572015年10月48222015年11月48682015年12月43252016年1月320图3step2：时间数列影响因素分析时间数列的影响因素主要有长期趋势、季节变动、循环变动和不规则变动。长期趋势是指受事物发展的根本因素制约而形成的事物在一段较长时期内基本趋势，可利用二次指数平滑法求解。[3]季节变动是受自然条件（气候）、社会条件（节假日、风俗）影响的。在影响房间预订量的因素中，季节、工作日/周末、法定假日、暑假都属于季节变动因素。循环变动具有周期长、规律性弱且不稳定的特点，因此在建立模型时，我们对循环变动因素不予考虑。不规则变动受偶然因素和意外条件的影响，我们在进行假设时不考虑其对房间预订量的影响。因此，在进行时间数列预测分析时，我们仅考虑长期趋势和季节变动因素。Step3：利用matlab（程序见附件3）画出全年每天预定房间数的趋势图，如图2图2：预定房间走势图可以清晰的观测到2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势。所以可以直接用时间序列预测模型结合前三个月的数据，对2016年前三个月每天的预订房间数量进行预测。4.3建立时间序列预测[4]模型时间序列预测是以时间数列所能反映的社会经济现象的发展过程和规律性，进行引伸外推，预测其发展趋势的方法。从预定房间走势图来看2015年1月至3月每天的预定房间数几乎趋于一个稳定的变化趋势，所以拟采用在一次指数平滑基础上加以改进的二次指数平滑时间序列来进行预测，并且由于原始数据有90个并且真实可靠，故平滑法采用的初始值以第一天的数据。模型求解与结果分析

利用matlab（程序见附录4）进行运算求得2016年前三个月的预测值趋势图见图3时『可-预订房间数40no2040€030100距1月1曰的天数图3：2016年前三个月的预测值趋势图预测的部分具体数据见表4，完整的数据见附录5表4：2016年一月到三月理论上每天预定房间数量的预测表日期房间数日期房间数日期房间数日期房间数1-10.00001-241.34732-16-0.02533-93.66821-22.00001-251.07062-17-0.05353-106.04541-31.70001-260.83582-18-0.07583-117.46041-41.64001-270.63722-19-0.09313-126.45201-52.38501-280.47002-20-0.10623-136.18061-62.01461-290.32982-21-0.11583-145.13901-71.69441-300.21302-22-0.12243-154.44391-81.81811-310.11622-23-0.12653-167.04651-92.12022-10.03672-24-0.12853-178.05411-101.77362-20.17192-250.27113-188.09971-111.87522-30.08972-260.81213-199.32561-121.55872-40.22222-270.87213-208.55541-131.28672-50.13732-280.72213-219.28591-141.85362-60.26722-291.19333-229.89481-151.53422-70.17973-12.19313-2311.99911-161.66022-80.10713-22.84093-2417.37351-174.16562-90.44723-32.78733-2520.32311-183.49332-100.33823-42.73603-2621.79901-192.91392-110.24653-53.28713-2723.0165由于预定的房间数量为整数，但是在预测中出现了小数，所以需要对数据进行处理。原则上只有有0.001个房间也要写为一个房间，但是对于房间数量预测出的数值，一方面反映了预定房间的数量；另一方面反映了预定一个的概率。如

果预定房间的概率小于0.5则说明预定房间的可能性不大，所以对于该数据处理可以采用四舍五入法。2016年一月到三月实际上每天预定房间数量的预测如表5。表5：2016年一月到三月实际每天预定房间数量的预测表日期房间日期房间日期房间日期房间1-201-2512-1703-1041-321-2612-1803-1161-421-2712-1903-1271-521-2812-2003-1361-621-2902-2103-1461-721-3002-2203-1551-821-3102-2303-1641-922-102-2403-1771-1022-202-2503-1881-1122-302-2603-1981-1222-402-2713-2091-1322-502-2813-2191-1412-602-2913-2291-1522-703-113-23101-1622-803-223-24121-1722-903-333-25171-1842-1003-433-26201-1932-1103-533-27221-2032-1203-633-28231-2132-1303-753-29221-2222-1403-843-30211-2322-1503-943-31201-2422-160五、模型检验问题二，时间序列预测模型的检验：对问题二预测的结果进行差分分析（matlab程序见附件4）具体的分析图见图4图4：差分分析图由差分分析误差条状图可以知道，预测值和去年的实际值呈现一阶差分趋势，表明时间序列预测模型的预测结果是符合理想。六、模型评价与推广模型优点问题一，1、考虑了人流到达的不规律性，将人流到达假设为服从泊松分布，进而考虑了低峰期和高峰期两种情况下的疏导模型，低峰期为游客规划了一条最短路径；2、高峰期将游乐场合理的划分为四个区域，进行建议式的疏导，这样既保证了为游客提供了最优游览方案，又让游客拥有自主选择权，互惠互利。问题二，1、简单易行，便于掌握，能够充分运用原时间序列的各项数据；2、计算速度快，对模型参数有动态确定的能力，精度较好。模型缺点问题一，没有考虑一些游客喜欢按照自己的方案游玩的特殊情况，让该疏导模型出现拥堵的情况增加了一定的不确定性。问题二，不能反映事物的内在联系，不能分析两个因素的相关关系，在处理问题时可能存在一定误差。模型改进问题一的方法相当于在为游客提供建议方案，让游客自由选择游玩路径，并不一定能达到我们预期想要的结果。因此需要寻求更优的算法对模型进行求解，例如：利用计算机仿真模拟等方法对模型进行改进。问题二采用时间序列二次指数平滑法来进行预测，只是这种预测的趋势呈现一种平稳、线性的形式，可能会对结果产生较大的误差。需要寻求更精确的方法进行预测。而在预测类方法中有时间序列预测、灰色预测、神经元网络、差分方程等，经过深层次分析，利用时间序列预测模型和灰色预测中的GM（1,1）模型分别进行预测，最后使两种预测进行加权平均法得出一个组合预测模型模型来进行预测，使得出的结果更精确。模型推广TSP模型运用广泛，可以用于最短路径类问题的求解；分区域疏导游客模型可以解决多服务地点高峰期的疏导类问题；时间数列预测模型，可以用于预测公司收入等问题。参考文献[1]司守奎,孙兆亮.数学建模算法与应用.北京.国防工业出版社.2015年4月.58-61.陈治佳，王曦，何苗•大型游乐场快速通道优化模型与仿真模拟J].哈尔滨工业大学学报.哈尔滨.第39卷第7期.2005年09月.101-103.丛国超，朱翼隽•批量到达的多服务台排队模型求解[J].成都信息工程学院学报.成都.第22卷第1期.2007年01月.98-100.司守奎，孙兆亮.数学建模算法与应用.北京.国防工业出版社.2015年4月.170-173.附录附录1运行环境：lingollMODEL:SETS:Entertainment/1..10/:U;!U(I)=sequenceno.ofEntertaiment;LINK(Entertainment,Entertainment):DIST,!DISTmatrix,itneednotbesymmetric;X;!X(I,J)=1ifweuselinkI,J;ENDSETSDATA:DIST=03006001050:350155090010506002503000300750650125012001350900550600300045050095010501400950850105075045009505001150155014001300350650500950012005509004506001550125095050012000650105015001800900120010501150550650040085011501050135014001550900105040004508006009009501400450150085045003502505508501300600180011508003500;ENDDATAN=@SIZE(Entertainment);MIN=@SUM(LINK:DIST*X);@FOR(Entertainment(K):@SUM(Entertainment(I)lI#NE#K:X(I,K))=1;@SUM(Entertainment(J)lJ#NE#K:X(K,J))=1;!Weakformofthesubtourbreakingconstraints;!Thesearenotverypowerfulforlargeproblems;@FOR(Entertainment(J)lJ#GT#1#AND#J#NE#K:U(J)>=U(K)+X(K,J)-(N-2)*(1-X(K,J))+(N-3)*X(J,K)));!MaketheX's0/1;@FOR(LINK:@BIN(X));!Forthefirstandlaststopweknow...;@FOR(Entertainment(K)lK#GT#1:U(K)<=N-1-(N-2)*X(1,K);U(K)>=1+(N-2)*X(K,1));END附录2运行环境：Lingdl附录2sets:state/1..57O/:p;endsetslamda=500;mu=10.375;rho=lamda/mu;s=4;k=570;lamda*pO=mu*p(l);(lamda+mu)*p(l)=lamda*p0+2*mu*p(2);@for(state(i)li#gt#l#and#i#lt#s:(lamda+i*mu)*p(i)=lamda*p(i-l)+(i+l)*mu*p(i+l));@for(state(i)|i#ge#s#and#i#lt#k:(lamda+s*mu)*p(i)=lamda*p(i-l)+s*mu*p(i+l));lamda*p(k_l)=s*mu*p(k);pO+@sum(state:p)=l;P_lost=p(k);lamda_e=lamda*(l-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-l/mu;endmode2:sets:state/1..58O/:p;endsetslamda=500;mu=10.4375;rho=lamda/mu;s=4;k=580;lamda*pO=mu*p(l);(lamda+mu)*p(l)=lamda*p0+2*mu*p(2);@for(state(i)|i#gt#1#and#i#lt#s:(lamda+i*mu)*p(i)=lamda*p(i-l)+(i+l)*mu*p(i+l));@for(state(i)|i#ge#s#and#i#lt#k:(lamda+s*mu)*p(i)=lamda*p(i-l)+s*mu*p(i+l));lamda*p(k_l)=s*mu*p(k);pO+@sum(state:p)=l;P_lost=p(k);lamda_e=lamda*(l-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-l/mu;endmode3:sets:state/1..260/:p;endsetslamda=200;mu=2.9;rho=lamda/mu;s=5;k=260;

lamda*pO=mu*p(l);(lamda+mu)*p(l)=lamda*p0+2*mu*p(2);@for(state(i)li#gt#l#and#i#lt#s:(lamda+i*mu)*p(i)=lamda*p(i-l)+(i+l)*mu*p(i+l));@for(state(i)li#ge#s#and#i#lt#k:(lamda+s*mu)*p(i)=lamda*p(i-l)+s*mu*p(i+l));lamda*p(k_l)=s*mu*p(k);pO+@sum(state:p)=l;P_lost=p(k);lamda_e=lamda*(l-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-l/mu;EndMode4:sets:state/1..180/:p;endsetslamda=200;mu=2.5;rho=lamda/mu;s=4;k=180;lamda*pO=mu*p(l);(lamda+mu)*p(l)=lamda*p0+2*mu*p(2);@for(state(i)|i#gt#1#and#i#lt#s:(lamda+i*mu)*p(i)=lamda*p(i-l)+(i+l)*mu*p(i+l));@for(state(i)|i#ge#s#and#i#lt#k:(lamda+s*mu)*p(i)=lamda*p(i-l)+s*mu*p(i+l));lamda*p(k_l)=s*mu*p(k);pO+@sum(state:p)=l;P_lost=p(k);lamda_e=lamda*(l-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-l/mu;Endmatlabr2007amatlabr2007aa=xlsread（'酒店预定历史数据2015年.xls'）;%将“附件2“更名为“酒店预定历史数据2015年.xls”，并保存在matlab子文件夹下x=a（:,l）;y=a（:,2）;plot（x,y,'k'）;gridon;xlabel（'距1月1日的天数）;ylabel（'预订房间数'）；title（'2015年时间-预订房间数'）；附录4运行环境：matlabr2007a时间序列预测程序：matlabR2007aclcclearyt=load('酒店预定历史数据2015年.txt');%将预订房间数数据以列的形式保存在“酒店预定历史数据2015年.txt”中，并保存在matlab子文件夹下n=length(yt),alpha=0.1;st1(1)=yt(1);st2(1)=yt(1);fori=2:nst1(i)=alpha*yt(i)+(1-alpha)*st1(i-1);st2(i)=alpha*st1(i)+(1-alpha)*st2(i-1);endxlswrite('酒店预定历史数据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');gri