全国研究生数学建模竞赛答案.docx

1、全国研究生数学建模竞赛答案参赛密码 （由组委会填写）全第十二届“中关村青联杯”全国研究生数学建模竞赛学 校上海交通大学参赛队号队员姓名1.2.3.参赛密码 （由组委会填写） 第十二届“中关村青联杯”全国研究生数学建模竞赛题 目 基于聚类分析的Hopfield网络求解旅游路线规划问题摘 要：本文围绕游遍201个5A级景区旅游问题进行了分析，对无费用限制的旅游时间问题、有时限的旅游费用问题利用聚类分析方法和连续的Hopfield网络分别建立了数学模型并设计了每条旅游线路具体的行程表，最后对求解结果进行了分析与验证。问题一在无费用限制情况下，要求用最少的时间游遍所有201个景点。第一步，利用聚类分析

2、方法对201个景点进行聚类。以按省份分类为主，按地理位置分类为辅，考虑实际环境，综合各自的优势为一体，最终划分出20个区域。第二步，根据Hopfield网络的有关方法，以景点间的消耗时间为参考量，建立了适用于问题一的Hopfield网络的计算模型。并用matlab语言编写模型的程序文件，在matlab软件中运行后得出各个区域内的最优旅游路线。第三步，结合题干中所有的旅游限制条件，设计出前往各个区域对应的旅游线路具体行程表。第四步，计算得出游遍201个景点的最短时长为11年。问题二在十年时间限制条件下，要求用最少的费用游遍所有201个景点。第一步，根据题目中的条件，针对十年期间的游玩总费用，建立

3、定价模型。第二步，仍然采用问题一的聚类分析方法的结果，将201个景点聚类成20个区域。第三步，针对问题二的具体情况，以景点间的消耗时间为参考量，对Hopfield网络的计算模型进行改进，得出各个区域内的最优旅游路线。第四步，设计出十年游遍所有景点的最低费用路线，总费用为287486.2元。问题三在前两个问题的基础上，规划出更适合全国旅游爱好者的游玩路线，并以北京市的旅游爱好者为例，给出最佳旅游路线；同时，依据当代旅游爱好者和相关旅游部门的现状，给出合理的建议，以便旅行者获得更好的旅行体验，相关部门提供更好的服务质量。问题四。最后，对整个数学模型进行了总结分析，并作出客观评价。关键词：聚类分析

4、Hopfield网络 matlab 定价模型 最优旅游路线 最佳体验一.问题重述 旅游活动正在成为全球经济发展的重要动力之一，随着我国国民经济的快速发展，人们生活水平得到很大提升，越来越多的人积极参与有益于身心健康的旅游活动。附件给出了全国201个5A级景区的名单，全国高速公路，全国火车、高铁、飞机班次等信息。一位自驾游爱好者拟按这些附件制定旅游计划。根据该旅游爱好者的个人偏好，景点位置及开放时间的实际情况，在旅行中需要达到以下条件：（1）该旅游爱好者每年有不超过30天的外出旅游时间，每年外出旅游的次数不超过4次，每次旅游的时间不超过15天；（2）根据个人偏好，每个5A级景区的游览时间不得小于

5、附件中的要求，最长逗留时间不得超过附件中最少时间的2倍；（3）基于安全考虑，行车时间限定于每天7:00至19:00之间，每天开车时间不超过8小时；（4）若是全天游览，则开车时间控制在3小时内；若是半天游览，开车时间控制在5小时内；（5）在高速公路上的行车平均速度为90公里/小时，在普通公路上的行车平均速度为40公里/小时；（6）该旅游爱好者计划在每一个省会城市至少停留24小时，以安排专门时间去游览城市特色建筑和体验当地风土人情（不安排景区浏览）；（7）选择高铁出行要求当天乘坐高铁的时间不超过6个小时，乘坐高铁或飞机的当天至多安排半天的景区游览;（8）景区开放时间统一为8:00至18:00；（9

6、）旅行中租车费用300元/天，油费和高速过路费另计，租车和还车需在同一城市；（10）住宿费简化为省会城市和旅游景区200元/人天，地级市150元/人天，县城100元/人天；高速公路的油耗加过路费平均为1.00元/公里，普通公路上油耗平均为0.60元/公里； 根据上述条件，需要解决下面问题：（1）该旅行者出行先通过高速公路到达与景区邻近的城市，再自驾到景区。以其常住地在西安市为例，规划设计旅游线路，试确定游遍201个5A级景区至少需要几年？给出每一次旅游的具体行程（每一天的出发地、行车时间、行车里程、游览景区）。（2）若出行方式考虑乘坐高铁或飞机到达与景区相邻的省会城市，而后租车自驾到景区游览。

7、根据附件材料，建立数学模型设计一个十年游遍所有201个5A景区、费用最优、旅游体验最好的旅游线路，给出每一次旅游的具体线路（含每次具体出行方式；每一天的出发地、费用、路途时间、游览景区、每个景区的游览时间）。 （3）在（2）的基础上加以推广，为全国的自驾游爱好者规划设计类似的旅游线路，进而给出常住地在北京市的自驾游爱好者的十年旅游计划；根据上述三问的结果给旅游爱好者和旅游有关部门提出建议。 （4）根据国家5A级旅游景区评定的相关信息，更合理地规划该旅游爱好者的十年旅游计划。二.模型的建立与求解2.1连续的Hopfield网络概述反馈网络达稳定状态时可以时系统的能量达极小，因而可用于一些最优化问

8、题的计算，如何把实际问题的目标函数表达成下述二次型的能量函数是一个关键问题。或 常用的是连续型Hopfield网络1，如图1所示，每一神经元可由一个（有正反向输出的）放大器模拟，输入端并联的电阻和电容可模拟生物神经元的时间常数，互相连线间的电导Tij则模拟各神经元间突触的特性（相当权系数）。该网络的微分方程为 函数gi常用Sigmoid函数：，u0可控制斜率，时变为阶跃函数。 若g-1()为单调增且连续，Cj0，Tji=Tij，则沿系统轨迹有，当且仅当时，其中，为系统的能量函数。以上表明，随着时间的演变，在状态空间内的网络总是朝着能量函数E减小的方向运动，网络稳定时E取极小值。图1 连续Hop

9、field神经网络的电路形式用Hopfield网络求解优化问题的一般步骤：（1）用罚函数法写出问题的目标函数设优化问题如下： 则目标函数为，其中为足够大的常数，取值可以互不相同。令I与前式中的E相等可定出各连接权Tij的值来。（2）写出网络的动态方程Hopfield网络是一个梯度系统，所以它满足对于常用的连续网络有（3）选择合适的初值，使网络按动态方程演化直到收敛为止。2.2 用Hopfield网络求解旅行商问题（TSP问题）对于N城市的TSP问题，任何一个城市在最终路径上的访问次序可用一个N维向量来表示，就需要N个神经元。如在5-TSP中，设城市1在第3个被访问，则对应的向量为V(1)=00

10、100。N城市TSP问题需用N*N个神经元来实现，而每行每列都只能有一个1，其余为0，该阵称为换位矩阵。换位矩阵中1的和为N，所构造的函数极小值对应于最短路径。构造与TSP相对应的能量函数1-3：A、B、C、D为拉格朗日常数，均为正数。当解合法时前三项为0；当达到最优解时，第四项最小，其值对应于最短路径。在Hopfield网络运行时，采用并行算法：=1.0；u0为符号函数的参量，u0越小，符号函数的离散化程度越高。在进行迭代前，要对uxi赋初值，不妨令，是在（-0.1，0.1）均匀分布的随机数。在迭代时，t为运算步长。因为能量在极小值时变化最慢，所以将能量函数E变化小到一定程度作为结束标志，即

11、。如果超过了一定的迭代次数（如1500次）仍没有收敛，则强行终止。2.3 旅游路线规划问题一的分析与求解假设该旅游爱好者每年外出旅游次数为x次，每次旅游时间为y天。根据条件，该旅游爱好者每年有不超过30天的外出旅游时间，每年外出旅游的次数不超过4次，每次旅游的时间不超过15天，则有： xy 30 x 4 y 15问题一中旅游爱好者需要以最短的时间将201个5A景点全部游览过，则有必要考虑旅行过程中的耗时因素，主要有（1）每次从西安到某个景点的来回时间t1；（2）高速公路及过道的限速t2；（3）到达景点的时间及景点的参观时长t3。针对耗时因素t1，可近似考虑其与从西安出发的次数x成正相关关系，正

12、相关系数为a，则有： t1 = ax 要使得t1时长短，需要减小x的值，同时为了使每年可以去到尽量多的地方，需要该旅游爱好者每年安排30天出行，求得x=2，y=15。 则该旅行者需要每年外出旅游2次，每次15天，从而可以有效地减小t1的时间，缩短整个规划中的时长。2.3.1 聚类分析思想聚类分析指将物理或抽象对象的集合分组为由类似的对象组成的多个类的分析过程4。聚类分析就是通过在相似的基础上收集数据来分类，达到数据简化的目的。聚类分析包括两类方式，（1）层次聚类（Hierarchical Clustering），包括合并法、分解法、树状图；（2）非层次聚类，包括划分聚类、谱聚类。201个5A级

13、别的景点覆盖整个中国，如果直接对其进行数据分析和处理，工作量是非常大的。因此我们需要对其进行聚类分析，形成更少数量的点，便于下一步的规划与设计。图2 全国201个5A景点一览图从地图上可以看到，西安市位于中原地区，黄色五角星表示5A级别的景点，这些景点分布在西安市的四周，这在位置上有利于聚类分析的分类，采用合并法，现有两种分类方式，（1）按照省份分类；（2）按照地理位置分类，即相邻的景点划分一类。两类分类方法的优缺点如下：（1）按照省份分类 优点：a. 分类方便，可直接通过省界线确定各区域； b. 同一个省内交通方便，距离较近，可大大缩短旅游时间； c. 符合附件中景点的分类方式，数据处理较为

14、方便； d. 符合现代人游玩的方式，可尽情享受省特色文化。 缺点：相邻省份的某些景点可能更加接近，按照省份分类反而会稍增耗时。（2）按照地理位置分类 优点：相邻景点交通便利，耗时较短。 缺点：a. 数据庞大，分区复杂，工作量大； b. 地图上相邻的某些景点之间交通可能并不发达，反而耗时更大； c. 不能充分感受省特色文化。综上，为了简化数据的处理，为了更好地体验，问题一中采用两种分类结合的方式对中国的5A景点进行划分，以按省份分类为主，按地理位置分类为辅，考虑实际环境，综合各自的优势为一体，最终划分出20个区域，分别是：江苏、浙江、福建、江西、黑龙江、广东、云南、四川、安徽（江苏、陕西）、青海

15、-甘肃（宁夏）、湖北（宁夏）、上海-山西（宁夏）、北京-天津（河北）、山东（河北）、吉林-辽宁-内蒙、湖南(重庆)、河南(陕西)、广西-海南、贵州（重庆）、新疆-西藏，其中符号“-”后的省份被完全包含到该区域中去，括号内的省份仅部分地区合并到该地区。以下所有问题的解决，都是以这20个区域为基础进行规划的。2.3.2 求解的具体过程由于全国的省份众多，这里不再一一说明。下面以旅游景点较多的浙江省为例进行说明。浙江省的旅游景点有：杭州西湖风景区、温州乐清市雁荡山风景区、舟山普陀山风景区、杭州淳安千岛湖风景区、嘉兴桐乡乌镇古镇旅游区、宁波奉化溪口滕头旅游景区、金华东阳横店影视城景区、嘉兴南湖旅游区、

16、杭州西溪湿地旅游区、绍兴市鲁迅故里沈园景区、衢州市开化根宫佛国文化旅游区、湖州市南浔区南浔古镇景区。为了减少计算的复杂性，提高计算结果的准确性，我们把在一个城市的多个景点合并成一个景点来进行计算，并将总游览时间作为该景点的游览时间。比如浙江省的杭州西湖风景区与杭州西溪湿地旅游区两者都在杭州市区，可以被当做一个景点来计算，当最后安排详细行程时再分拆成两个景点进行行程安排。另外由于该旅游爱好者计划在每一个省会城市至少停留24小时，以安排专门时间去游览城市特色建筑和体验当地风土人情，因此我们将每个省的省会当做一个景点来计算，并将游览时间定为8个小时。由于旅游路线规划问题一需要求解的是旅游时间最优，查

17、询各个景点的行车时间（资料来源：XX地图5）与游览时间（资料来源：附件1），建立如下图所示的各景点之间所花费时间的表格。表格中每个花费时间包括行车时间与将要游览景区的游览时间。其中左侧景区为出发地，上侧的景区为目的地。 表1 浙江省各景区间的花费时间表浙江省景区杭州西溪湿地旅游区、西湖杭州淳安千岛湖风景区湖州市南浔区南浔古镇景区嘉兴南湖旅游区嘉兴桐乡乌镇古镇旅游区金华东阳横店影视城景区宁波奉化溪口滕头旅游景区衢州市开化根宫佛国绍兴市鲁迅故里沈园景区温州乐清市雁荡山风景区舟山普陀山风景区杭州西溪湿地旅游区，西湖010.25.55.35.26.16.257.35.4512.611.7杭州淳安千岛湖

18、风景区10.207.57.47.27.68.056.77.2513.6513.5湖州市南浔区南浔古镇景区9.511.505.054.677.336.678.66.613.112.05嘉兴南湖旅游区9.311.45.0504.96.656.058.45.512.111.4嘉兴桐乡乌镇古镇旅游区9.211.24.674.907.36.68.756.171311金华东阳横店影视城景区10.111.67.336.657.305.7575.6510.711.85宁波奉化溪口滕头旅游景区10.2512.056.676.056.65.7508.15.751110.04衢州市开化根宫佛国11.310.78.6

19、8.48.7578.107.91314.05绍兴市鲁迅故里沈园景区9.4511.256.65.56.175.655.757.901211.25温州乐清市雁荡山风景区12.613.659.18.196.7798013.05舟山普陀山风景区11.713.58.057.477.856.0410.057.2513.050编写基于Hopfield网络的matlab语言的m文件6（具体程序请见附录一），输入上表数据，运行得出浙江省内的最优路线方案如下矩阵（程序结果vv矩阵的第212列）。杭州西溪湿地旅游区，西湖 0 0 0 1 0 0 0 0 0 0 0杭州淳安千岛湖风景区 0 0 0 0 0 0 0 0

20、 0 0 1湖州市南浔区南浔古镇景区 1 0 0 0 0 0 0 0 0 0 0嘉兴南湖旅游区 0 0 1 0 0 0 0 0 0 0 0嘉兴桐乡乌镇古镇旅游区 0 1 0 0 0 0 0 0 0 0 0金华东阳横店影视城景区 0 0 0 0 0 1 0 0 0 0 0宁波奉化溪口滕头旅游景区 0 0 0 0 0 0 0 1 0 0 0衢州市开化根宫佛国 0 0 0 0 0 0 0 0 0 1 0绍兴市鲁迅故里沈园景区 0 0 0 0 1 0 0 0 0 0 0温州乐清市雁荡山风景区 0 0 0 0 0 0 0 0 1 0 0舟山普陀山风景区 0 0 0 0 0 0 1 0 0 0 0矩阵为一

21、个换位矩阵，矩阵每一列仅有一个1，代表每次游览的城市。从第1列到第11列分别给出了依次游览的城市顺序。TSP问题将会回到起点，因此第11列的游览景区结束后将会回到第1列的游览景区。上表对应的游览景区顺序如下：湖州市南浔区南浔古镇景区嘉兴桐乡乌镇古镇旅游区嘉兴南湖旅游区杭州西溪湿地旅游区，西湖绍兴市鲁迅故里沈园景区金华东阳横店影视城景区舟山普陀山风景区宁波奉化溪口滕头旅游景区温州乐清市雁荡山风景区衢州市开化根宫佛国杭州淳安千岛湖风景区湖州市南浔区南浔古镇景区在地图上标明各景区(黄色五角星)和游览路线(红色直线)后，如下图所示。由图可知整个游览过程在地图上基本成一个环装，说明游览过程是合理的。其中

22、宁波奉化溪口滕头旅游景区这一个景区似乎不太合理，但结合景区间的公路情况我们可以知道，从金华东阳横店影视城景区去往舟山普陀山风景区的公路是途径宁波奉化溪口滕头旅游景区的，而舟山普陀山风景区前往温州乐清市雁荡山风景区也会途径宁波奉化溪口滕头旅游景区，因此舟山普陀山风景区和宁波奉化溪口滕头旅游景区的前后顺序不会增加游览的总时间。图3 浙江省浏览路线图将通过Hopfield算法得出的上述浙江省内游玩的最优时间路线命名为路线a。为了验证路线a是否为时间最少的路线，我们随即抽取两条路线作为比较。路线b： 嘉兴桐乡乌镇古镇旅游区 嘉兴南湖旅游区 杭州西湖风景区 杭州西溪湿地旅游区 湖州市南浔区南浔古镇景区

23、杭州市 绍兴市鲁迅故里沈园景区 舟山普陀山风景区 宁波奉化溪口滕头旅游景区 金华东阳横店影视城景区 温州乐清市雁荡山风景区 衢州市开化根宫佛国文化旅游区 杭州淳安千岛湖风景区 嘉兴桐乡乌镇古镇旅游区。路线c： 杭州淳安千岛湖风景区 嘉兴桐乡乌镇古镇旅游区 嘉兴南湖旅游区 杭州西湖风景区 湖州市南浔区南浔古镇景区 杭州市 绍兴市鲁迅故里沈园景区 舟山普陀山风景区 宁波奉化溪口滕头旅游景区 金华东阳横店影视城景区 温州乐清市雁荡山风景区 杭州西溪湿地旅游区 衢州市开化根宫佛国文化旅游区 杭州淳安千岛湖风景区。参考XX地图查询出来的时间，通过计算，路线a耗时为86.06小时，路线b耗时为105.94

24、小时，路线c耗时为125.84小时。路线时间对比如图所示，路线a耗时最短，从而验证了模型算法的正确性。图4 三种路线耗时对比图 由于各个省内的游览过程是一个巡回的过程，在地图上呈现为一个环形，当我们从西安前往该省游玩时，该省存在一个景点离西安在距离上最近，选择这个景点作为西安去这个省的目的地可以减少整体的行车时间。如浙江省内的湖州市南浔区南浔古镇景区离西安最近，选择这个城市为西安前往浙江的目的地。这样最终前往浙江的这一次旅游的具体行程如下表。表2 浙江旅游的具体行程天数出发时间出发地游览景区行车时间行车里程到达时间游览时间游览完毕时间17:00西安湖州市南浔区南浔古镇景区8:00015:002

25、7:00西安湖州市南浔区南浔古镇景区6:34128313:3437:00湖州市南浔区南浔古镇景区嘉兴桐乡乌镇古镇旅游区0:4319.47:434:0011:43311:43嘉兴桐乡乌镇古镇旅游区嘉兴南湖旅游区1:0237.312:454:0016:45316:45嘉兴南湖旅游区杭州西溪湿地旅游区1:279718:124:0012:00412:00杭州西溪湿地旅游区杭州西湖风景区0:2812.812:284:0016:28416:28杭州西湖风景区杭州市0016:2824:0016:28516:28杭州市绍兴市鲁迅故里沈园景区1:2069.917:484:0011:48611:48绍兴市鲁迅故里

26、沈园景区金华东阳横店影视城景区1:39118.713:274:0011:44711:44金华东阳横店影视城景区舟山普陀山风景区3:52286.515:368:0013:36813:36舟山普陀山风景区宁波奉化溪口滕头旅游景区2:2315215:594:009:5999:59宁波奉化溪口滕头旅游景区温州乐清市雁荡山风景区2:53188.412:528:0010:521010:52温州乐清市雁荡山风景区衢州市开化根宫佛国文化旅游区4:5639515:484:009:48119:48衢州市开化根宫佛国文化旅游区杭州淳安千岛湖风景区2:29121.312:178:0010:171210:17杭州淳安千

27、岛湖风景区湖州市南浔区南浔古镇景区3:27260.213:444:0017:441217:44湖州市南浔区南浔古镇景区西安1:16019:00137:00湖州市南浔区南浔古镇景区西安8:00015:00147:00湖州市南浔区南浔古镇景区西安5:18012:18前面已经论述了将每条线路尽可能的安排成接近15天可以减少浪费在来回路途上的时间。下面以湖南省的线路举例说明，这条线路同时又包括了重庆省的部分景点。因为这些景点没有被游览完成，如果单独前往一次重庆需要浪费4天来回的路程，这样加上游玩的时间将达到7天左右。而湖南省离重庆省非常近，在湖南省游玩完成后还有多余的时间即可游玩重庆省的景点，从而减少

28、了一次出行的时间，大大减少了总的旅游行程时间。前往湖南(重庆)的这一次旅游的具体行程如下表。表3 湖南(重庆)旅游的具体行程表天数出发时间出发地游览景区行车时间行车里程到达时间游览时间游览完毕时间17:00西安张家界武陵源天门山旅游区8:00015:0027:00西安张家界武陵源天门山旅游区4:4297411:4216:0017:4247:00张家界武陵源天门山旅游区湘潭韶山旅游区4:4912911:494:0015:49415:49湘潭韶山旅游区长沙市宁乡县花明楼景区0:3931.816:284:0010:28510:28长沙市宁乡县花明楼景区郴州市东江湖旅游区4:09335.814:374:008:3768:37郴州市东江湖旅游区衡阳南岳衡山旅游区2:43190.511:204:0015:20615:20衡阳南岳衡山旅游区长沙岳麓山橘子洲旅游区2:55134.918:154:0012:00712:00长沙岳麓山橘子洲旅游区长沙市0:00012:0024:0012:00812:00长沙市岳阳岳阳楼君山岛景区3:01182.415:014:009:0199:01岳阳岳阳楼君山岛景区万盛黑山谷-龙鳞石海风景区5:00014:01107:00岳阳岳阳楼君山岛景区万盛黑山谷-龙鳞石海风景区6:05872.813