A算法的改进课程设计毕业设计.doc

课程设计说明书 no.1a*最短路径算法的改进1.课程设计的目的本课程设计是学习数据通信与通信网技术课程必要的教学环节。由于该课程是专业必修课，需要通过实践巩固基础知识，为使学生取得最现代化的设计技能和研究方法，课程设计训练也就成为了一个重要的教学环节。通过对路由算法的设计和实现，达到进一步完善对通信网基础及应用课程学习的效果。2.设计方案论证算法具体的形式包括：确定起点的最短路径问题：即已知起始结点，求最短路径的问题。确定终点的最短路径问题：与确定起点的问题相反，该问题是已知终结结点，求最短路径的问题。在无向图中该问题与确定起点的问题完全等同，在有向图中该问题等同于把所有路径方向反转的确定起点的问题。确定起点终点的最短路径问题：即已知起点和终点，求两结点之间的最短路径。全局最短路径问题：求图中所有的最短路径。（1）floyd求多源、无负权边的最短路。用矩阵记录图。时效性较差，时间复杂度o(v3)。floyd-warshall算法（floyd-warshall algorithm）是解决任意两点间的最短路径的一种算法，可以正确处理有向图或负权的最短路径问题。floyd-warshall算法的时间复杂度为o(n3)，空间复杂度为o(n2)。floyd-warshall的原理是动态规划：设di,j,k为从i到j的只以(1.k)集合中的节点为中间节点的最短路径的长度。若最短路径经过点k，则di,j,k = di,k,k-1 + dk,j,k-1；若最短路径不经过点k，则di,j,k = di,j,k-1。因此，di,j,k = min(di,k,k-1 + dk,j,k-1 , di,j,k-1)。在实际算法中，为了节约空间，可以直接在原来空间上进行迭代，这样空间可降至二维。 沈 阳 大 学课程设计说明书 no2floyd-warshall算法的描述如下：for k 1 to n do for i 1 to n do for j 1 to n do if (di,k + dk,j o（e*lgv）。当是稀疏图的情况时，此时e=v*v/lgv，所以算法的时间复杂度可为o（v2） 。若是斐波那契堆作优先队列的话，算法时间复杂度，则为o（v*lgv + e）。（3）bellman-ford求单源最短路，可以判断有无负权回路（若有，则不存在最短路），时效性较好，时间复杂度o（ve）。bellman-ford算法是求解单源最短路径问题的一种算法。单源点的最短路径问题是指：给定一个加权有向图g和源点s，对于图g中的任意一点v，求从s到v的最短路径。与dijkstra算法不同的是，在bellman-ford算法中，边的权值可以为负数。设想从我们可以从图中找到一个环路（即从v出发，经过若干个点之后又回到v）且这个环路中所有边的权值之和为负。那么通过这个环路，环路中任意两点的最短路径就可以无穷小下去。如果不处理这个负环路，程序就会永远运行下去。 而bellman-ford算法具有分辨这种负环路的能力。a*（a-star)算法是一种静态路网中求解最短路最有效的直接搜索方法。注意是最有效的直接搜索算法。之后涌现了很多预处理算法（alt，ch，hl等等），在线查询效率是a*算法的数千甚至上万倍。公式表示为： f(n)=g(n)+h(n),其中 f(n) 是从初始点经由节点n到目标点的估价函数，g(n) 是在状态空间中从初始节点到n节点的实际代价，h(n) 是从n到目标节点最佳路径的估计代价。保证找到最短路径（最优解的）条件，关键在于估价函数h(n)的选取： 沈 阳 大 学课程设计说明书 no.3估价值h(n)实际值,搜索的点数少，搜索范围小，效率高，但不能保证得到最优解。 算法是一种启发式搜索算法, 在路径规划中得到广泛的应用, 其中启发函数的设计尤其重要。本文针对路径规划问题, 对a* 算法作了以下改进: 一是在估价函数中考虑以距离和方向两个要素, 通过归一化处理解决了单位不统一的问题; 二是利用 k- d树空间索引结构, 动态加载节点信息, 减小内存使用空间。实验结果表明, 改进后的a* 算法的搜索效率得到了明显的提高。最经典的最短路径搜索算法是 dijkstra算法,属于遍 历搜索, 它简单易用并总能搜索到最短路径。但是当网 络中节点数较多时,该算法搜索的结点数量会很大,效率非常低。因此有人提出了启发式搜索算法,如: 局部择 优搜索法、最好优先搜索法、a*算法等。启发式搜索就是 在状态空间中, 对每一个搜索的位置进行评估, 得到最好的位置,从而在这个位置进行搜索直到搜索到目标为止。 目前在路径优化领域, 最流行的启发式搜索算法当属由 h ar,t nilsson, raphael等人首先提出的 a* 算法。它利用启发函数来估计任意点到目标点的远近程度, 从而减少 搜索空间, 提高搜索效率。许多文献都对 a* 算法进行了研究, 并且都提出在估价函数中引入距离和方向两个要素。但是距离和方向的单位是不统一的,所以在利用时会出现一些问题, 本文针对这一问题进行了研究,并对估价函数进行了改进。另外,为了进一步提高算法的运行效率,本文在算法运行结构上,采用 k- d树空间索引结构,降低内存存储空间。实验结果证明了改进后算法的合理性和可行性。3.设计的过程与分析 a*算法是建立在dijkstra和 bfs(最好优先搜索 )算法基础上的。它的整体框架采用遍历搜索法, 但是它采 用了启发函数来估计地图上任意点到目标点的费用, 从而可以很好地选择搜索方向。a*算法引入了当前节点j的启发函数 f*( j),当前节点 j的启发函数定义为: f*(j)=g(j)+h*(j)(1)其中 g(j)是从起点到当前节点 j的实际费用的量度, h*(j)是从当前节点j到终点的最小费用的估计。 沈 阳 大 学课程设计说明书 no.4注意到若 h *(j) = 0,即没有利用任何全局信息,这时a* 算法就变成了普通的 dijkstra算法。因此普通的dijkstra算法可 看作 a* 算法的特殊情形。h* (j)要满足相容性条件: 即不能高于节点 j到终点的实际最小费用。可以证明,如果估价函数满足相容性条件, 且原问题存在最优解,则 a* 算 法一定能够求出最优路径 。a* 算法的优点在于利用启发函数, 使搜索方向更加智能的趋向于终点,所以其搜索深度较小, 搜索的节点数少,故占用存储空间少,如图 1所示。 图1 a*算法与dijkstra算法搜索区域对比 a* 算法的估价函数 a* 算法的关键在于设计一个合适的启发函数。有文 献对其特点进行了分析, 认为启发函数 f* (j)值是非递减 的,只要能够满足相容性条件即估价函数 h* ( j)小于节点 j到目标点的实际费用,它生成的路径一定是最优的。许多文献都在估价函数的构造中引入了距离和方向两个要素 ,即: h*(j) = w1* l+w2*如图2所示。其中 l 表示当前节点到终点的欧氏距离,表示起点到当前节点的线段与当前节点到终点的线段的夹角即 sje,有文献也用到了中间节点与关联节点的线段和关联节点与终点的线段夹角 nje。w1和w2分别是距离和角度的加权值,w1和w2 的取值范围分别为 055- 0.65和0.35-0.45。但距离和角度的单位不统 一的问题不能忽略,即使角度的单位为弧度、距离的单位 为千米,它们之间也很有可能出现距离值远大于角度值 (即 l )的情况,特别是在大区域路径规划过程中,问题更明显。 沈 阳 大 学课程设计说明书 no.5而当 l时, 方向不再有约束力,而使得估价函数失去意义。 a* 算法的运行结构 当构造合适的估价函数后就要考虑算法的运行,目前大多数方法是将全部数据读入到内存当中, 然后搜索 最短路径。从理论上讲, a* 算法可以通过搜索更少的节 点完成最短路径的搜索。但是算法运行时,必须要考虑两个问题:一是数据读取的速度。即使可以很快地将数 据读入到内存中,我们也还要考虑第二个问题, 即系统内 存的大小。如果系统内存足够大,在算法运行过程中,也将会出现对同一节点进行重复的搜索,从而降低算法的运行效率。针对数据的读取问题,有学者提出了基于限 制区域的a* 算法,减小数据的加载量。但是由于a* 算法本身就是一种有损算法, 这种方法很有可能搜索 不到最短路径, 特别是在考虑道路属性信息和交通限制信息时。图2 估价函数构造示意图 改进的 a* 算法 （1） a* 算法估价函数的改进针对a* 算法估价函数所出现的问题,我们将距离和角度进行归一化处理,即首先计算当前节点所有关联节点相应的距离和角度值,然后求它们的平均值即 l，从而使得估价函数 变为: h* (j) = w1* l + w2* (5) 其中: l= l/l (6)=/(7) 归一化处理以后,只考虑距离和角度对路径的贡献,而不必考虑距离和角度的数值大小。从而避免了距离和 角度单位不统一的问题。虽然算法的运行要增加计算 量, 但是我们可以通过进一步减小算法的搜索空间, 改善 算法的运行结构来提高算法的搜索效率。 沈 阳 大 学课程设计说明书 no.6 针对算法运行效率的问题, 建立k-d树空间索引结构, 动态加载路网数据, 从而提高算法效率不失为一个好方法。 k- d树索引结构是 k( k 1)的二叉检索树,主要用于索引多属性的数据或多维点数据, 它可以通 过坐标快速的访问区域中的路网数据。在算法执行过程中,并非开始就装载所有的路网数据,而是根据算法的需要,读取节点的相关信息,或删除节点信息,虽然会增加运算过程中的 i/o运算,但是这样可以避免无效数据的大量装载,占用大量的内存空间。例如,首先判断当前节点是否在确定的范围内, 如果不在则装载相应区域的数据,如果在确定的范围内,则读取数据的相关信息,并进行节点扩展。然后,在此基础上计算路段的启发值,搜索最短路径。 （2） a* 算法的实现在算法的实现过程中,要构造两个链表。分别存储待扩展的节点和已扩展的节点,分别称为open表和 close 表。算法步骤如下: 初始化设置。将起点信息加载到 open表中, close 表赋值为空。令g(j) = 0。 搜索距离当前节点最近的节点,即求f* (j)的最小值,将节点j从open表中删除并加载到close表中。判断节点 j是否为终点,如果是,终点转向步骤4;否则转向步骤3。判断节点j的节点信息是否在确定的范围内,如果在范围内,则扩展节点j;否则加载节点j的节点信息并进行扩展。转向步骤2。从节点j开始, 利用回溯的方法输出起始节点到目 标节点的最优路径, 以及最短距离,算法终止。在算法的运行过程中,避免了对同一节点的重复访 问, 极大地缩小了搜索空间, 从而缩短了算法的运行时间。对改进后的a*算法进行了实验, 在估价函数归一化处理前后的最短路径搜索结果如图3(a), (b)所示。图3a改进前的搜索路径 沈 阳 大 学课程设计说明书 no.7图3b 改进后的搜索路径（3） 实验采用郑州市地图,节点2606个,路段4127条,在 core i5, 8gb内存的 pc机上运行。与其他算法的实验结果进行了对比, 结果见表1。表 中: t 表示临时标记节点个数,n表示永久标记节点的个数, d表示规划路径长度。表1 算法比较图4 临时标记节点个数比较 沈 阳 大 学课程设计说明书 no.8图5 永久标记节点个数比较 将表1数据中的临时标记节点,永久标记节点比较关 系用图4、图5表示。通过实验由图4可以看出,归一化处理后的a* 算法的搜索区域明显减小,由表1可以看出 a* 算法的搜索效率要比 dijkstra算法的搜索效率高。由图4、图5可知,改进后的a* 算法的搜索效率明显要高,在利用 k- d树建立 空间索引结构以后,使得搜索的点数明显减少,特别是在区域比较大时,改进后的a* 算法的效率提高得更加明显。需要指出的是, 由于a* 算法本身就是有损算法,所以其寻找到的路径长度一般要比 dijkstra算法的规划结果要长,但由于所选的道路更加合理,所以改进后算法的搜索结果更加实用。4设计体会通过这次的课程设计，利用a*算法求解最短路径，加深了对a*算法的了解与认识。课程设计环节中把教材里面的理论知识与实践联系起来，利用理论知识指导实践仿真，取得很多的收获。学习这个算法开始的时候会觉得比较难，花了一天的时间看资料理解算法的原理。在知道原理后的第二天开始编写代码，同样也花了一天时间。最后是程序的显示的优化。总之通过学习这个算法编写程序，可以慢慢的从参考别人的代码转变自己知道原理后编写代码这个过程会比较慢需要不断的联系。a* 算法作为一种启发式搜索算法在路径规划中得到 了非常广泛的应用。利用启发函数减小搜索空间,从而提高搜索效率,因此启发函数在 a* 算法中起着关键的作用。针对a* 算法启发函数中距离和角度两个要素单 沈 阳 大 学课程设计说明书 no.9位不统一的问题做了研究,提出将距离和角度归一化处理,并且利用 k-d树的空间索引结构,减少搜索空间。试验结果表明,改进后的a* 算法的效率明显提高。5参考文献 1t.h.cormen,c.e.leiserson,r.l.rives,teta.lintroductiontoalgorithmsm . mcgraw- h il,l 2001. 2 stevenm. lavatle p lanning a lgorithm sm . cambridge university press, 2006. 3李擎，宋顶立.两种改进的最优路径规划算法j.北京科技大学学报，2011.4周春辉，李诗高.dijkstra算法与a算法研究j.软件导刊，2007.5马进通信网分析m北京：人民交通出版社，2003 沈 阳 大 学课程设计说明书 no.10附录部分程序代码：package com.ubird.astar.core.exception;public class astarpositionexception extends runtimeexception private static final long serialversionuid = 5968574301453821996l; public astarpositionexception(string msg) super(msg); package com.ubird.astar.core;import java.util.arraylist;import java.util.hashmap;import java.util.list;import java.util.map;public class astarmap private list openlist = new arraylist(); private map closemap = new hashmap(); private boolean isfind = false; private list path = new arraylist(); public static final int state_barrier = 2; astarnode target; astarnode source; int astardata; public astarmap(int xgridnum, int ygridnum) this.astardata = new intygridnumxgridnum; this.source = new astarnode(0, 0); this.target = new astarnode(xgridnum - 1, ygridnum - 1); public astarnode gettarget() return this.target; 沈 阳 大 学课程设计说明书 no.11public astarnode getsource() return this.source; public int getastardata() return this.astardata; public list find() init(); astarnode current = null; while (!isend() & (!isfind() current = getminfnodefromopenlist(); if (isachieve(current) buildpath(current); this.isfind = true; else add2closemap(current); for (astarnode neighbor : getneighbor(current) if (neighbor = null) | (isinclosemap(neighbor) | (iscannotgo(current, neighbor) continue; boolean isbetter = true; astarnode nodefromopenlist = getnodefromopenlist(neighbor); if (nodefromopenlist != null) neighbor = nodefromopenlist; int tg = neighbor.distinctg(current); isbetter = tg = 0) & (x = 0) & (y this.astardata.length); private astarnode getminfnodefromopenlist() int index = 0; int minf = (astarnode)this.openlist.get(index).getf(); int length = this.openlist.size(); 沈 阳 大 学课程设计说明书 no.13for (int i = 1; i astarnode.getf() minf = astarnode.getf(); index = i; return (astarnode)this.openlist.remove(index); private boolean isachieve(astarnode current) return current.equals(this.target); private void init() initparameter(); initopenlistandclosemap(); addsource2openlist(); private void initparameter() this.isfind = false; this.source.init(this.target); this.path.removeall(this.path); private void initopenlistandclosemap() clearopenlist(); this.closemap.clear(); public void clearopenlist() this.openlist.removeall(this.openlist); private void addsource2openlist() add2openlist(this.source); 沈 阳 大 学课程设计说明书 no.14 private void add2openlist(astarnode node) this.openlist.add(node); private void add2closemap(astarnode node) this.closemap.put(node.tostring(), node); private boolean isfind() return this.isfind; private boolean isend() return this.openlist.size() = 0; public void loaddata(int data, int barrierval, int clearval) if (data = null) return; if (data.length != this.astardata.length) | (this.astardata0.length != data0.length) throw new runtimeexception(new data should be as large as the old astar map data); for (int i = 0; i this.astardata.length; i+) for (int j = 0; j this.astardatai.length; j+) if (dataij = barrierval) this.astardataij = 2; else if (dataij = clearval) this.astardataij = 0; public void inittargetandsource(int x, int y) this.source.setx(this.target.getx(); this.source.sety(this.target.gety(); this.target.setx(x); this.target.sety(y); package com.ubird.astar.core; 沈 阳 大 学课程设计说明书 no.15import com.ubird.astar.core.exception.astarpositionexception;public class astarnode implements comparable private static final int g_1 = 10; private static final int g_2 = 14; private static final int g_3 = 10; private static final int g_4 = 14; private static final int g_5 = 10; private static final int g_6 = 14; private static final int g_7 = 10; private static final int g_8 = 14; private static final int g_0 = 10; private int g; private int h; private int f; private int x; private int y; private astarnode father; public astarnode(int x, int y) this.x = x; this.y = y; public int getx() return this.x; public void setx(int x) this.x = x; public int gety() return this.y; public void sety(int y) this.y = y; public astarnode getfather() return this.father; 沈 阳 大 学课程设计说明书 no.16public void setfather(astarnode father) this.father = father; public void init(astarnode target) this.g = 0; this.h = heuristiccostestimate(this, target); this.f = (this.g + this.h); public int heuristiccostestimate(astarnode source, astarnode target) return (math.abs(source.x - target.x) + math.abs(source.y - target.y) * 10; public int compareto(astarnode o) return this.f o.f ? -1 : 1; public boolean equals(object obj) if (obj = null) | (!(obj instanceof astarnode) return false; astarnode node = (astarnode)obj; return (node.x = this.x) & (node.y = this.y); public string tostring() return this.x + , + this.y; public void recalculatorgandh(astarnode father, astarnode target) this.g = distinctg(father); this.h = heuristiccostestimate(this, target); this.f = (this.g + this.h); this.father = father; 沈 阳 大 学课程设计说明书 no.17public int distinctg(astarnode father) int offsetx = this.x - father.x; int offsety = this.y - father.y; int distinct = 0; if (offsetx = 0) & (offsety = -1) distinct = 10; else if (offsetx = 1) & (offsety = -1) distinct = 14; else if (offsetx = 1) & (offsety = 0) distinct = 10; else if (offsetx = 1) & (offsety = 1) distinct = 14; else if (offsetx = 0) & (offsety = 1) distinct = 10; else if (offsetx = -1) & (offsety = 1) distinct = 14; else if (offsetx = -1) & (offsety = 0) distinct = 10; else if (offsetx = -1) & (offsety = -1) distinct = 14; else throw new astarpositionexception(unvalid relation between current node( + this + ) and fater node( + father + ); return distinct + father.g; public int getg() return this.g; public boolean isbetter(astarnode node) return isgbetter(node); public boolean isgbetter(astarnode node) return this.g + distinctg(node) node.g; 沈 阳 大 学课程设计说明书 no.18public boolean isfbetter(astarnode node) return this.f node.f; public int getf() return this.f; package com.ubird.astar.demo;import com.ubird.astar.ui.astarmenubar;import com.ubird.astar.ui.astarpanel;import com.ubird.astar.ui.controlpanel;import com.ubird.astar.ui.statuspanel;import com.ubird.astar.ui.uframe;import java.awt.container;import javax.swing.swingutilities;public class astardemo public static void main(string args) swingutilities.invokelater(new runnable() public void run() uframe frame = new uframe(); astarpanel astarpanel = new astarpanel(15, 15, 60, 40); statuspanel statuspanel = new statuspanel(); astarpanel.setstatuspanel(statuspanel); frame.getcontentpane().add(astarpanel, center); frame.getcontentpane().add(new controlpanel(astarpanel), north); frame.getcontentpane().add(statuspanel, south); frame.setjmenubar(new astarmenubar(astarpanel); frame.pack(); frame.setvisible(true); frame.setdefaultcloseoperation(3); astarpanel.requestfocus(); ); 沈 阳 大 学课程设计说明书 no.19在算法的实现过程中, 要构造两个链表。分别存储 待扩展的节点和已扩展的节点, 分别称为 open表和 c lose 表。算法步骤如下: 1)初始化设置。将起点信息加载到 open表中, c lose 表赋值为空。令 g( j) = 0。 2)搜索距离当前节点最近的节点, 即求 f* ( j)的最小 值, 将节点 j从 open表中删除并加载到 c lose表中。判断 节点 j是否为终点, 如果是, 终点转向步骤 4; 否则转向步 骤 3。 3)判断节点 j的节点信息是否在确定的范围内, 如果 在范围内, 则扩展节点 j ; 否则加载节点 j的节点信息并进 行扩展。