《最短路径算法在物流运输中的应用.docx》由会员分享,可在线阅读,更多相关《最短路径算法在物流运输中的应用.docx(33页珍藏版)》请在三一办公上搜索。
1、本科生毕业设计(论文)题 目:学生姓名:学 号:院 系:线性表的设计和实现张三专业年级:指导教师:201107011153基础科学学院信息技术系2012级信息与计算科学专业李四注:1论文封面单独打印一张纸;中英文摘要正反 打印一张纸;目录、正文、参考文献、致谢、附录 均独立正反打印!2.部分专业对格式有特殊要求的,教学院(系)可 自行商定。摘 要随着现代物流业的发展,如何优化和配置物流的运输路径成为了一个热点的问题。 其中,最具代表性的问题就是如何在一个道路网络中选择两点之间的合适路径,使其 距离最短。为了解决这个问题,本文介绍了两种最常用的最短路径求解方法一一 DIJKSTRA算法与FLOY
2、D算法,分析了它们的适用范围以及时间复杂度。最后,对一 个具体的航空公司物流配送问题进行了求解,得到了理论最优路径。关键词:最短路径问题;DIJKSTRA算法;物流运输ABSTRACTWith the development of modern logistics industry, how to optimize and configure the transport path of logistics has become a hot issue. Among them, the most representative problem is how to select the approp
3、riate path between two points in a road network to minimize the distance. In order to solve this problem, this paper introduces two most common shortest path solutions Dijkstra algorithm and Floyd algorithm, and analyzes their application range and time complexity. Finally, a specific airline logist
4、ics distribution problem is solved, and the theoretical optimal path is obtained.Keywords: Minimum path problem; Dijkstra algorithm; Logistics transportation目 录第一章引言11.1研究背景11.2研究现状11.2.1最短路径算法研究现状11.2.2最短路径算法分类2第二章最短路径问题的基本理论知识32.1最短路问题的定义32.2 最短路问题的 Dijkstra 算法32.2.1 Dijkstra 算法的局限性32.2.2 Dijkstra
5、 算法求解步骤32.2.3 Dijkstra算法的时间复杂度42.2.4简单案例分析42.3 最短路问题的 Floyd 算法52.3.1算法定义52.3.2算法思想原理52.3.3算法过程描述62.3.4算法适用范围62.3.5算法简单实例6第三章 实际案例分析73.1问题描述73.1.1问题的背景及假设73.1.2符号说明73.2模型的建立与求解83.2.1 模型一83.2.2 模型二10第四章总结154.1优点154.2缺点15参考文献16致谢17附录18附录A实际案例背景数据18第一章引言1.1研究背景在现实生活中中,我们经常会遇到图类问题,图是一种有顶点和边组成,顶点代表 对象,在示意
6、图中我们经常使用点或者原来表示,边表示的是两个对象之间的连接关系, 在示意图中,我们使用连接两点G点直接按的下端来表示。顶点的集合是V,边的集合 是E的图记为GV,E,连接两点u和v的边用e(u,v)表示。最短问题是图论中的基础问 题,也是解决图类问题的有效办法之一,在数学建模中会经常遇到,通常会把一个实际 问题抽象成一个图,然后来进行求的接任意两点之间的最短距离。因此掌握最短路问题 具有很重要的意义。1.2研究现状本节主要讨论两个方面的问题,首先简要回顾最短路径算法研究现状,然后概要总 结最短路径算法分类。1.2.1最短路径算法研究现状最短路径问题一直是计算机科学、运筹学、地理信息科学等学科
7、领域的研究热点。 国内外大量专家学者对此问题进行了深入研究。经典的图论与不断发展完善的计算机数据结构及算法的有效结合使得新的最短路 径算法不断涌现。常用的路径规划方法有:平行最短路径搜索算法,蚁群算法,基于矩阵 负载平衡的启发算法,EBSP*算法和Dijkstra算法等。创门在空间复杂度、时间复杂度、 易实现性及应用范围等方面各具特色但是因为Dijkstra算法可以给出最可靠的最短路 径,并且容易实现,所以备受青睐和并被广泛应用。经典的Dijkstra算法的时间复杂度为,直接应用到大规模城市路网时,最短路 径查询时间难以令人接受,专家学者纷纷开展Dij kstra优化算法研究,概括起来,以往
8、研究者主要是从5个方面对最短路径算法进行性能优化:(1)基于数据存储结构的优化, 以空间换取时间;(2 )基于路网规模控制的优化;(3)基于搜索策略的优化;(4 )优先级 队列结构的优化;(5 )基于双向搜索的并行计算优化。1.2.2最短路径算法分类由于问题类型、网络特性的不同,最短路径算法也表现出多样性。(1) 按照最短路径问题分类,最短路径问题通常可分为2个基本类型:一是单源最 短路径问题,即查找某一源点到其余各点的最短路径;另一类是查找某个节点对之间的 最短路径。最短路径问题具体可细分为以下几种,单源最短路径问题,单对节点间最短路径、 所有节点间最短路径、k则最短路径、实时最短路径、指定
9、必经节点的最短路径以及前 N条最短路径问题等,本文的研究范畴属于单对节点间最短路径问题。(2) 按照网络类型和表示方法分类,网络可以分为稀疏网络和非稀疏网络,常用的 表示方法有邻接矩阵和邻接表。邻接矩阵方法能够在门时间内查询到任意两个节点之间是否有一条边,它的空 间复杂度为夕。现实生活中网络节点往往很多,动辄上万,而且是稀疏网络居 多,比如城市路网,所以用邻接矩阵表示既不现实,又浪费空间。邻接表是另一种存储网络拓扑的数据结构,它是一种链式存储结构,对于交通网络 等稀疏图 ,采用邻接表数据结构存储网络拓扑数据空间复杂度仅 为6:我+兴,不存在存储空间的浪费。邻接表数据结构已被证明是网络表达中最有
10、 效率的数据结构,在最短路径算法中得到了广泛应用。第二章最短路径问题的基本理论知识2.1最短路问题的定义最短路问题(short-path problem):若网络中的每条边都有一个数值(长度、成本、 时间等),则找出两节点,(通常是源节点和目标节点)之间总权和最小的路径就是最 短路问题。最短路问题是网络理论解决的典型问题之一,可用来解决管道铺设,线路安 装,厂区布局和设备更新等实际问题。2.2最短路问题的Dijkstra算法2.2.1 Dijkstra算法的局限性在了解和使用某种算法之前,我们先要明白这种算法有怎样的局限性。只有深入理 解来每一种算法的局限性,才能根据具体的问题选择合适的算法来
11、求解。Dijkstra算法最大的局限性在于不能够处理带有负边的图,即图中任意两点之间的 权值必须非负。如果某张图中存在长度为负数的边,那么Dijkstra算法将不再适用,需 要寻找其他算法求解。2.2.2 Dijkstra算法求解步骤(1) 先给图中的点进行编号,确定起点的编号。(2) 得到图的构成,写出图的矩阵(u , u ) (u , u )000 n G =. (u , u ) (u , u )n 0n n(3) 根据要求求出发点S到终点E的最短距离,那么需要从当前没被访问过的结点集合unvist=u I u e U23.n中找到一个距离已经标记的点的集合中vist=u | u e 1,
12、2,3. n的最短距离,得到这个顶点;(4) 利用这个顶点来松弛其它和它相连的顶点距离S的值(5) 重复步骤(2)和(3),直到再也没有点可以用来松弛其它点,这样我们就得到了 由起点S到其它任意点的最短距离。2.2.3 Dijkstra算法的时间复杂度我们可以用大符号将Dijkstra算法的时间复杂度表示成边数 m与顶点数n 的函数。Dijkstra算法最简单的实现方法是用一个链表或者数组来存储所有顶点的集合 Q,因此搜索Q中最小元素的运算(Extract-Min(Q)只需要线性搜索Q中的所有元 素,据此我们可以知道算法的时间复杂度是 U:。对于边数少于n2稀疏图而言,我们可以用邻接表来更有效
13、的实现Dijkstra算法。为了达到这一目的,需要将一个二叉堆或者斐波纳契堆用作优先队列来寻找最小 的顶点(Extract-Min)。当用到二叉堆的时候,算法所需的时间为 心,.w + i沽二, 斐波纳契堆能稍微提高一些性能,使得算法运行时间达到9 : ,” +。2.2.4简单案例分析给出对应的结点之间的关系(表2-1为对应的结点之间的关系)长度ABCDEA02151010B201115C15111207D10102003E105730表2-1需要注意的是,其中(A,B)= 2表示结点A到B的长度为2 第一步:进行编号,假定A点即为起点。第二步:得到图02151010201115G151102
14、071012003105730第三步:首先从起点A开始找到距离A最近的点,那就是A点了;第四步:把A点标记到已经用过的的集合IM /I 用A来更新其它点W, D.扩到起点的距离得到的集合A B C D E dist =0 2 15 10 10表示起点到B,C,D,E的距离分别为2,15,10,10第五步:重复上述步骤:得到A B C D E dist =vist = A, B unvist = C, D, E0 2 13 3 7继续重复上述步骤,最后的到vist = A,B,C,D,E, unvist =0,得到的A BC D Edist = 0 2 1336即最短路求解完毕。2.3最短路问题
15、的Floyd算法2.3.1算法定义除了 Dijkstra算法,另外一种简单的最短路算法是floyd算法,它也经常被用于解决含有有向图或者是负权的最短路径问题,并且能够用于计算有向图的传递闭包。该算 法的时间复杂度为日Vi,空间复杂度为(.矿!。2.3.2算法思想原理Floyd算法是非常常见的使用动态规划来寻找最优路径的算法。如果我们用简单的 语言解释,总体要实现的目标是找到从点i到点j的最短路径。如果我们换一个角度, 从动态规划的角度观察,那就必须得要对这个目标进行一个重新的诠释。从任意节点i到任意节点j的最短路径不外乎2种可能,1是直接从i到j,2是从i 经过若干个节点k到j。所以,我们假设
16、Dis(i,j)为节点u到节点v的最短路径的距离, 对于每一个节点k,我们检查.万鬲;一,局=。*;温:)+庭(一j)是否成立,如果成立,证 明从i到k再到j的路径比i直接到j的路径短,我们便设置 Ss(们项)=以田(挪)十功;M&J),这样一来,当我们遍历完所有节点k,Dis(i,j)中记录 的便是i到j的最短路径的距离。2.3.3算法过程描述(1)从任意一条单边路径开始。所有两点之间的距离是边的权,如果两点之间没 有边相连,则权为无穷大。(2)对于每一对顶点u和v,看看是否存在一个顶点w使得从u到w再到v比己 知的路径更短,如果是的话需要更新它。2.3.4算法适用范围无向图最短路问题; 稠
17、密图效果最佳;边权可正可负。2.3.5算法简单实例图3-2无向图根据图3-2,用Floyd算法找出任意两点的最短路径步骤如下表3-2:distk1distk2distk3MINA-B1371A-C135*1A-D3353B-C2262B-D*44*4C-D2462表3-2 Floyd算法步骤流程第三章实际案例分析3.1问题描述3.1.1问题的背景及假设网上购物一直是常见的消费方式,其依托于物流业逐渐蓬勃发展,每个送货人员都 需要以最快的速度送货,而且往往会发送多个地方,所以有必要设计耗时最小的路线。现在考虑一个快递公司,总部在地图上的O点,派送人员需要将货物发往城市很多, 如何设计交货方案,以
18、便花费最少的时间。物流地图如图1所示,下表中显示了每个点 的信息,假设托运人只能沿着这些连接的线行进,而不采取任何其他路线。(1)最大承载能力为50公斤,货量最大为1立方米。(2)调度员的平均速度为24公里/小时。假设每件货物要交出需要3分钟,为了简单起 见,在同一地点有几件商品,这些货物只需每3分钟一次交割即可。现在派送者将向50个地点发送100件货物。问题要求如下:1. 假设货物从1到30到指定地点并返回。设计最快的路线和方式来求出结果。要求标 记送货线路。2. 假设送货人员从上午8点开始交货,1到30天货物交货时间不能超过规定的时间, 请设计最快的路线和方式。要求标记送货线。具体数据请参
19、见附录。3.1.2符号说明Mi所载货物的质量总和;V所载货物的体积总和;imi第i件货物的质量;匕第i件货物的体积;d从i点到j点的距离权值;B任意两节点之间最短路径距离矩阵;3.2模型的建立与求解3.2.1模型一我们首先对题中所给的数据进行汇总分析,得出M = 0mi =48.5公斤, 30厂1V30=X、=088立方米,所以均未超出送货员的载重,所以送货员可以一次性将货物送 r 1完。而题中数据显示送货员需到达的节点数位22个(包括出发点O)如下表0131416171821232426273132343638394042434549表3-1节点数利用程序用Floyd算法我们可以得出任意两点
20、之间最短路径的距离矩阵3寸.其中 (i,j=122,(1) 先根据题目数据给初始矩阵B(i, j)赋值,其中没有给出距离的赋inf,以便于更新。(2) 进行迭代计算。对任意两点(i, j),若存在k,使B(i,k) + B(k, j) v B(i, j),则更新B(i, j) = B(i, k) + B(k, j)。(3) 直到所有点的距离不再更新停止计算。则得到最短路距离矩阵B(i, j)。由旅行售货员问题(TSP)建立矩阵d( ) m = 22 ;其中d 表示点i到点j的距离权值。d为对称 矩阵,令d=0。现求节点0到各节点再到节点0的最短距离,要求各线路上的权值最 小。设立变量匕&.),
21、其关系如下:当节点,和节点j连通,广L当节点和节点,不连通,气舟目标函数为寻找 一条从起点0到各节点再到节点0的最短距离,要求各线路上的权值和最小,故目标函(3-1)数为:最短路径min yMji1 厂 1(3-2)(1)对起始点0至少有一条路径出去,故有j j 1(2)对其余各节点,恰有一条路径进去,故有=广1ki(3)所有节点不出现闭合回路,约束为u(3-3)总的线性规划模型为:min Z = d x munIJ IJ IJ ,T J=1(3-4)(2):ki= 1,i = 2,3,., nk=11ki约束条件s.t. (3)u u + nx n 1,i,顶 1,2,., n利用lingo
22、软件编程算出在各约束条件下的最短路径距离、最短路径所经过的各节 点的顺序得:最短距离 min/ =沦1 7S 一.最短时间iiiirir =.比.各节点行进路线为:0一26一27一39一36一38一43一42一49一45一40一34一24一13一18一14一16一32一23一17一21一0图3-1节点路线表3.2.2模型二问题2题目增加了时间约束,所以我们需要在模型一的基础上进行改进。送货员从 早上8:00出发,需要分别在9:00、9: 30、10: 15、12:00之前件货物送到各指定点。 根据“时间要求越早,先送的原则”,将22个节点按时间限制划分为四个区域,然后分 阶段依次解决各区域的最
23、短路径,得出各出发点和各终点。从而得出总距离最短路径。首先我们在图中描出各节点坐标,找到各节点位置。如下图:图3-2节点位置表阶段1:要求9:00送到:限制在此时间段的节点为三个:13、18、24,送货员8:00从O点出发,需选择最 短路径在9:00之前将货物送达。根据各节点之间的距离和上图,我们很容易得出此段最 短路径出发点为18,终点为24,从而路线为18-13-24,再根据示以及问题1所得数 据,确定最优线路为18131924。最后验证结果:根据路径我们算得此路径距离:2182.0289+3113.4627+5704.3372=10999.83 m.从而得出此段用去的时间=10999.8
24、3*3/20=27.50min60min,从而可以知道按此路径 送货员能按时将货物送到。阶段2:要求9:30送到:根据题中信息知,限制在此时间的点为:31,34,40,45。同上阶段相同,结合数据和上图分析的出发点为31,终点为45,从而我们可以得到两 种行程路线:31344045或31403445。需要对两条路线进行对比优化。两种路线的行程总路程如下:路线12431313434404045路程(m)1780.14752324.74731630.7823217.0056路线22431314040343445路程(m)1780.14753954.92931630.7824847.7876表3-2
25、行程路线表对比两组数据,可以选定线路1为最佳方案。按此路径行进距离=1780.1475+3954.9293+1630.782+4847.7876=12213.6464 米。得出耗时=12213.6464*3/20=30.53min30min+60min-27.50min。即得出满足时间要求。阶段3:要求10:15送到:此时间要求共有四个指定地点:49, 42, 43,38。分析可得起点为42,终点为38,从而得到两种送货路线:42494338或42434938。两种路线的总路程如下:路线14542424949434338各段路程(m)2351.72281971.37642889.0501261
26、8.4442路线24542424343494938各段路程(m)2351.7228917.67372889.05015507.4943表3-3总路程表分析比较两组数据,可以选定线路1为最佳方案。同上对数据进行验证:按此路径行进距离=2351.7228+1971.3764+2889.0501+2618.4442=9830.5935 米得出耗时=9830.5935*3/20=24.58min45min.即能够按时件货物送到。阶段4:要求12: 00到达:此时间段共有十个指定地点:26,21,14,17,23,32,39,36,27,16。分析可确定36为起点。起点确定终点不确定。对此我们进行迭代算
27、法:即从出发点开始,找出到出发点的最近距离的一个点,然 后依次迭代计算,再以找出的点为出发点,找出到该点的最短距离的点,如此进行下去, 则可以找出一条较短的行进路线。首先以36为起点,具体计算过程如下:点36到其他点的距离(单位:m)如下:36272139321726231416距离2203.9172880.1783983.844061.1444704.0894808.6965373.0136176.9117470.654表3-4点36距离表所以确定27为第二次所到地点。点27到其他点的距离(单位:m)如下:273926213217231416距离1779. 9232604.7814796.5
28、216265.066620.4327576.938093.2549674.571表3-5点27距离表所以确定39为第三次所到地点。由线路图可知,离开39后需要回到27,才能送货到其它地点,则再根据上述表格选取 除39外距离27最近的地点的,即是第四次所到地点,易知26为第四次所到地点(途 经 31)。点26到其他点的距离(单位:m)如下:26211714233216距离2191.744015.6515488.4735790.1657102.0347887.806表3-6点26距离表可得21为第五次所到地点。点21到其他点的Euclid距离(单位:m)如下:211714233216距离1823.
29、911_3296.7333598.4254910.294 5696.066表3-7点21距离表可得17为第五次所到地点。点17到其他点的Euclid距离(单位:m)如下:1723143216距离1774.5149215.7233086.3833872.156表3-8点17距离表理论上讲应该选取点23,但根据线路图以及剩余送货的地点,综合考虑后选取次优解即 14为第六次送货地点。点14到其他点的Euclid距离(单位:m)如下:14162332距离2607.6813970.2375282.106表3-9点14距离表可得16为第七次所到地点。点16到其他点的距离(单位:m)如下:162332距离2
30、097.63409.5表3-10点17距离表可得23为第八次所到地点,32为终点。由以上结果可得最佳送货路线如下:362739(27)(31)2621 1714162332在图中标出路线:图3-3路线表所以综上考虑将各阶段的路径连接起来形成的最终最短路径为:0一18一13一19一24一31一34一40一45一42一49一43一38一36一27一39一27一31一26一21一17一14一18一23一32。得出最短路径距离minZ=54629.65m时间 minT=2.276h第四章总结4.1优点对于题中各数据的处理,采用的工具、方法比较先进,各种计算方法精确,误差 较小。在解决问题时,我们把原图
31、变为完全图,利用全局线性规划法充分发挥lingo软 件包求最优解功能。并且成功地对01变量进行了使用和约束,简化了模型建立难度, 同时给出了在各种约束条件下的最短路径的计算方法,具有较强的实用性和通用性,在 日上生活中经常可以用到。4.2缺点由于数据较多,难以对模型结果进行验证,只能一步一步的对模型进行优化。参考文献I 朱道立.运筹学M.高等教育出版社,2006. 周建兴.MATLAB从入门到精通.第2版M.人民邮电出版社,2012.3 徐立华.求解最短路问题的一个计算机算法J.系统工程,1989(5):46-51.4 林澜,闫春钢,蒋昌俊,等.动态网络最短路问题的复杂性与近似算法J.计算机学
32、 报,2007, 30(4):608-614.5 牛学勤,王炜.基于最短路搜索的多路径公交客流分配模型研究靠J.东南大学学 报(自然科学版),2002, 32(6):917-919.6 马良河,刘信斌,廖大庆.城市公交线路网络图的最短路与乘车路线问题J.数学的 实践与认识,2004, 34(6):38-44.7 魏航,李军,刘凝子.一种求解时变网络下多式联运最短路的算法J.中国管理科学, 2006, 14(4):56-63.8 张运河,林柏梁,梁栋,等.优化多式联运问题的一种广义最短路方法研究J.铁道 学报, 2006, 28(4):22-26.9 李帮义,姚恩瑜.最短路网络及应用J.系统工程
33、理论与实践,2000, 20(6):104-107.10 龙光正,杨建军.改进的最短路算法J.系统工程与电子技术,2002, 24(6):106-108.II 周经伦,吴唤群.受顶点数限制的最短路问题及其算法J.系统工程, 1996(5):37-44.12贺国先.集装箱公铁联运的费用加权最短路计算机算法J.铁道学报,2006,28(1):1-5.致谓寸 一级标题,中间空4格。空1行大学四年的学习生活即将结束,在此,我要感谢所有曾经教导过我的老师和关心过 我的同学,他们在我成长过程中给予了我很大的帮助。本文能够成功的完成,要特别感 谢我的导师XXX教授的关怀和教导。附录A实际案例背景数据1600
34、014000120001000080006000400020000z0200040006000800010000120001400016000图1快递公司送货地点示意图O点为快递公司地点,O点坐标(11000,8250),单位:米表1各货物号信息表货物号送达地点重量(公斤)体积(立方米)不超过时间1132.500.03169: 002180.500.03549: 003311.180.02409: 304261.560.035012: 005212.150.030512: 006141.720.010012: 007171.380.010912: 008231.400.042612: 0093
35、20.700.048112: 0010381.330.021910: 1511451.100.02879: 3012430.950.022810: 1513392.560.059512: 0014452.280.03019: 3015422.850.019010: 1516431.700.078210: 1517320.250.041212: 0018361.790.018412: 0019272.450.044512: 0020242.930.04209: 0021310.800.01089: 3022272.250.001812: 0023261.570.021012: 0024342.8
36、00.01039: 3025401.140.01559: 3026450.680.03829: 3027491.350.014410: 1528320.520.002012: 0029232.910.048712: 0030161.200.042912: 003111.260.02503221.150.05013331.630.04833441.230.00063551.410.03873660.540.00673770.700.01293880.760.03463992.140.008740101.070.012441111.370.051042122.390.042843130.990.0
37、04844141.660.049145150.450.020946162.040.009847171.950.032448182.120.055449193.870.026250202.010.032451211.380.041952220.390.000153231.660.050254241.240.053455252.410.001256261.260.005957270.420.022458281.720.058059291.340.037260300.060.040261310.600.027462322.190.050363331.890.049464341.810.0325653
38、51.000.005566361.240.017767372.510.036168382.040.011069391.070.044070400.490.032971410.510.009472421.380.045573431.310.012174441.260.000575450.980.041376461.350.024177472.120.023078480.540.054279491.010.056680501.120.028481250.790.001182462.120.049283322.770.003484232.290.005485200.210.049086251.290
39、.008887191.120.024988410.900.003889462.380.043490371.420.002091321.010.030092332.510.013393361.170.002094381.820.030895170.330.034596110.300.017297154.430.053698120.240.005699101.380.017510071.980.0493表250个位置点的坐标位置点X坐标(米)Y坐标(米)1918550021445560372705704373567052620995610080143571002522808716025259138
40、452680101193530501178503545126585418513763052001413405532515212559751615365704517141657385188825807519585581652078083552112770856022220088352314765905524779093302544359525261086096352710385105002856597652925809865301565995531939510100321483510365331250109003472801106535153051137536123901141537641011510381391511610399510120504083451230041493013650421326514145431418014215443