用MATLAB结合四种方法搜寻罗马尼亚度假问题
选修了cs的AI课,开始有点不适应,只能用matlab硬着头皮上了,不过matlab代码全网仅此一份,倒有点小自豪。
一、练习题目
分别用宽度优先、深度优先、贪婪算法和 A*算法求解“罗马利亚度假问题”。具体地图我这里不给出了,有兴趣的可以去搜索。即找到从初始地点 Arad到 目的地点 Bucharest 的一条路径。
要求:分别用文件存储地图和启发函数表,用生成节点数比较以上四种算法在同一问题求解时的效率,列表给出结果。
附:罗马尼亚度假问题图(图1.1)
图1.1 罗马尼亚度假问题
1.2 题目分析
本题要求分别用宽度优先、深度优先、贪婪算法和 A*算法求解“罗马利亚度假问题”。罗马尼亚度假问题本质上属于“图类”问题,该地图上共有20个地点,要求从Arad出发,到达Bucharest,从图中搜索到达路径,并比较四种方法的优缺点。因此主要的数据结构可以采用图存储的方法,搜索方法题目已经给出。
二、数据结构
2.1 图结构
图:由有穷、非空点集和边集合组成,简写成G(V,E),其中,G表示一个图,V表示图中的顶点,E表示图中的边。在本题,顶点为20个罗马尼亚城市,边则为相邻城市之间的距离。边之间有方向,图为有向图,无方向的图为无向图。本题所用的图为无向图。
尽管二维数组比较占用内存,但是由于MATLAB对矩阵运算非常方便,运算速度也很快,我采用二维数组的方法存储邻接矩阵。对20个地点编号1-20,其中Arad为3号地点,对边采用数值的方法,例如3号到4号距离为75,则令矩阵中点(3,4)的值为75。并令自身距离为0,不相邻的点之间也设为0。
2.2 队列结构
队列(Queue):是只允许在一端进行插入操作,在另一端进行删除操作的线性表。队列也是一种特殊的线性表,是一种先进先出的线性表。允许插入的一端称为表尾,允许删除的一端称为表头。我们将在广度优先搜索中使用到这个结构存储已搜索过的节点。其结构如图2.2所示。
图2.2 队列结构图
2.3 栈结构
栈(Stack):也是一种线性存储结构,栈中的数据元素遵守“先进后出”(First In Last Out)的原则,简称FILO结构。只能在栈顶进行插入和删除操作。我们将在深度优先搜索中使用到这个结构存储已搜索过的节点。其结构如图2.3所示。
图2.3 栈结构图
三、算法思想
3.1 宽度优先
宽度优先搜索算法(又称广度优先搜索)其别名又叫BFS( Breadth First Search)。属于一种盲目搜寻法,目的是系统地展开并检查图中的所有节点,以找寻结果。算法采用队列的数据结构,所有因为展开节点而得到的子节点都会被加进一个先进先出的队列中。其邻居节点尚未被检验过的节点会被放置在一个被称为 open 的队列,而被检验过的节点则被放置在被称为 closed 的容器中(open-closed表)算法自始至终一直通过已找到和未找到顶点之间的边界向外扩展,首先搜索和s距离为k的所有顶点,然后再去搜索和S距离为k+l的其他顶点。算法流程图如图3.1所示。
图3.1 DFS算法流程图
3.2 深度优先
深度优先搜索方法,又称DFS(Depth First Search),和树的先序遍历比较类似。假设初始状态是图中所有顶点均未被访问,则从某个顶点v出发,首先访问该顶点,然后依次从它的各个未被访问的邻接点出发深度优先搜索遍历图,直至图中所有和v有路径相通的顶点都被访问到。 若此时尚有其他顶点未被访问到,则另选一个未被访问的顶点作起始点,重复上述过程,直至图中所有顶点都被访问到为止。算法流程图如图3.2所示。
图3.2 BFS算法流程图
3.3 贪婪方法
贪婪算法(又称贪心算法)是指,在对问题求解时,总是做出在当前看来是最好的选择。为了解决问题,需要寻找一个构成解的候选对象集合,起初,算法选出的候选对象的集合为空。接下来的每一步中,根据选择函数,算法从剩余候选对象中选出最有希望构成解的对象。如果集合中加上该对象后不可行,那么该对象就被丢弃并不再考虑;否则就加到集合里。每一次都扩充集合,并检查该集合是否构成解。
本题中具体实现方法为,先进行深度搜索,但是不急进入堆栈操作,而是存储当前所有搜索到的点的距离,选择距离最短的点,并放弃搜索其他同一深度的点,进入堆栈操作。算法流程图如图3.3所示。
图3.3 贪婪算法流程图
3.4 A*方法
A*搜寻算法俗称A星算法。A*算法是比较流行的启发式搜索算法之一,被广泛应用于路径优化领域。它的独特之处是检查最短路径中每个可能的节点时引入了全局信息,对当前节点距终点的距离做出估计,并作为评价该节点处于最短路线上的可能性的量度。
本题中的实现方法为,同贪婪类似,A*就相当于有一个智慧的老人为搜寻的对象打分,搜索过程中将距离和打分值相加,作为新的距离即可。其算法流程图如3.4所示。
图3.3 A*算法流程图
四、关键代码
4.1 BFS方法
while tail~=head %判断 i=queue(tail); %取点for j=1:20 %搜索所有适合的节点if A(i,j)>=1 && isempty(find(flag==j,1))queue(head)=j; head=head+1; %数数 flag=[flag j]; %扩容 result=[result;i,j,A(i,j)]; %记录endendtail=tail+1; %队列增加 end
4.2 DFS方法
while top~=0 %判断 pre_len=length(stack); %记录栈长度i=stack(top); %取栈顶for j=1:20if A(i,j)>=1 && isempty(find(flag==j,1)) top=top+1; stack(top)=j; flag=[flag j]; re=[re;i,j,A(i,j)]; %记录 break; endend if length(stack)==pre_len %如果栈未增加,则出栈stack(top)=[];top=top-1;end end
4.3 贪婪方法
while top~=0 pre_len=length(stack); i=stack(top); record=[]; for j=1:20if A(i,j)>=1 && isempty(find(flag==j,1)) %记录所有相邻节点record=[record;i,j,A(i,j)] endendif isempty(record)breakend[s,k]= min(record(:,3,:)) %取距离最小节点i=record(k,1,:);j=record(k,2,:); if isempty(find(flag==j,1))top=top+1; stack(top)=j; flag=[flag j]; re=[re;i,j,A(i,j)]; endif length(stack)==pre_len stack(top)=[];top=top-1;end end
4.4 A*方法
绝大部分与贪婪算法相同,只是更新了距离值。
for j=1:20if A(i,j)>=1 && isempty(find(flag==j,1)) F(i,j)=A(i,j)+H(j);record=[record;i,j,F(i,j)]; %启发值endendif isempty(record)break end
4.5 反向寻址
所有的算法均采用相同的反向寻址方法。
while (1)x=find(re(:,2,:)==m) %找到到达目的地所有的经过地m=re(x,1,:) %迭代法反向寻找来的路径if 1-isempty(x)lujin=[city{re(x,1,:)},lujin];juli=juli+re(x,3,:)elsebreak end end
4.6 读取EXCEL
city={'Oradea','Zerind','Arad','Timisonra','Lugoj','Mehadia','Dobreta','Craiova','Rimmicu Vikea','Sibiu',...'Fagaras','Pitesti','Bucharest','Giurglu','Uiziceni','Hirsova','Eforie','Vaslui','Lasi','Neamt'}; %存储城市名 filename = 'mymap.xlsx'; sheet = 1; xlRange = 'C3:V22'; map = xlsread(filename,sheet,xlRange); %读取excel map(isnan(map)) = 0; %将不相邻的点之间也设为0
五、运行结果
BFS方法的运行结果显示路径为:{'Arad' 'Sibiu' 'Fagaras' 'Bucharest'}
DFS方法的运行结果显示路径为:{ 'Arad' 'Zerind' 'Oradea' 'Sibiu' 'Rimmicu Vikea' 'Craiova' 'Pitesti' 'Bucharest'}
贪婪方法的运行结果显示路径为:{ 'Arad' 'Zerind' 'Oradea' 'Sibiu' 'Rimmicu Vikea' 'Pitesti' 'Bucharest'}
A*方法的运行结果显示路径为:{ 'Arad' 'Sibiu' 'Rimmicu Vikea' 'Pitesti' 'Bucharest'}
比较见表5.1
表5.1 四种算法的运行结果
算法 |
生成节点数 |
求解时间 |
距离 |
BFS方法 |
11 |
3.725s |
450 |
DFS方法 |
12 |
3.057s |
762 |
贪婪方法 |
7 |
3.606s |
575 |
A*方法 |
7 |
3.049s |
418 |
注:求解时间包括计时函数自用时间
六、比较结论
得益于MATLAB高速的矩阵运算能力,四种方法均在3-4秒之间完成,速度相差不大,但是在生成节点数上,DFS方法搜索了12个节点最多,贪婪方法和A*方法均为7最少。比较四种搜索方法得到的搜索路径,有启发值的A*方法搜索到的路径距离最短,为418,其次是宽度优先搜索,距离为450,距离最长的路径是由DFS方法产生,为762,贪婪方法为575。通过比较我们可以得出如下结论:
- 四种搜索方法在处理小型网络的搜索问题时,速度相差不大。
- 贪婪方法和A*方法生成节点数较少,理论上能够更快搜索出到达路径,在处理大型图的问题时,会表现得比较明显。
- 贪婪方法每一步都是选择当前状态下的最优解进行搜索,很容易陷入局部最优,从而使得搜索时间延长。
- 尽管BFS方法和DFS方法都一定可以找到路径,但是BFS方法搜索到的路径距离要明显优于DFS方法。
---恢复内容结束---
转载于:https://www.cnblogs.com/Hangingter/p/7784042.html
用MATLAB结合四种方法搜寻罗马尼亚度假问题相关推荐
- matlab图像拼接融合(四种方法)
matlab图像拼接的四种方法 1.直接拼接, 2.亮度调整后拼接, 3.按距离比例融合, 4.亮度调整后按距离比例融合 流程: 1.读入左,右图,并取出重合部分,并转化为亮度图 2.分别把每 ...
- matlab中单独存图_Matlab中图片保存的四种方法
Matlab 中图片保存的四种方法 matlab 的绘图和可视化能力是不用多说的, 可以说在业内是家喻户晓的. Matlab 提供了丰富 的绘图函数,比如 ez** 系类的简易绘图函数, surf . ...
- Matlab中将数据保存为txt或dat格式四种方法
转载于 Tsingke 老师的文章--Matlab中将数据保存为txt或dat格式四种方案 - Tsingke - 博客园. ...
- Simulink中构造时变传递函数的四种方法
1. 原由 传递函数是针对线性时不变(Linear Time-Invariant,LTI)系统定义的,严格来说,时变传递函数的说法并不十分严谨.但 "时变"和"非线性&q ...
- java 中lock,java中lock获取锁的四种方法
在java接口中会存放着许多方法,方便线程使用时的直接调用.对于lock接口大家都不陌生,我们已经初步对概念进行了理解.那么在获取锁的方法上想必还不是很清楚.下面我们就lock获取锁的四种方法分别进行 ...
- 简单介绍C语言使用四种方法初始化结构体
这篇文章说明了什么是结构体,介绍了结构体的概念和使用优点,在C语言中如何使用和初始化结构体方法,通过详细的代码展开进行说明,希望该篇文章对你有所帮助 什么是结构体 在实际问题中,一组数据往往有很多种不 ...
- svg鼠标响应事件的四种方法(其中两种可支持火狐)
svg鼠标响应事件的四种方法 鼠标响应事件的四种方法,以click事件为例. Mouse Events - SMIL <?xml version="1.0" encoding ...
- oracle调整字段精度的四种方法
oracle调整字段精度的四种方法: 01_执行用户_ddl/dml_表名_注释(建表/授权/同义词) 调整方式一:(精度只能调大不能调小) alter table table_name_a MODI ...
- (转).NET导出Excel的四种方法及评测
.NET导出Excel的四种方法及评测 导出Excel是.NET的常见需求,开源社区.市场上,都提供了不少各式各样的Excel操作相关包.本文,我将使用NPOI.EPPlus.OpenXML.Aspo ...
最新文章
- Python-内建函数
- 微软程序员利用测试账户套现千万美元,或面临 20 年监禁
- c语言程序算一元二次方程,如何用C语言来计算一元二次方程
- linux发布tomcat项目项目流程,linux下tomcat的配置及项目的部署流程
- linux的RPM支持图片格式,Linux下rpm、yum和源码三种安装方式详细介绍
- Unity3d UGUI 通用Confirm确认对话框实现(Inventory Pro学习总结)
- 【数据库原理实验(openGauss)】认识DBMS系统
- Lightroom Classic 教程,如何将照片从 Lightroom 移至Ps,在 Ps 中为照片应用滤镜?
- linux下启动管理,rEFInd启动管理器配置文件详解
- 怎样获得淘宝商品类目API接口
- 黄杏元《地理信息系统概论》考研复习考点精讲(五)
- 家藏四库图书系列书目(经典)
- 计算机如何调整对比度,怎么调整电脑显示器亮度和对比度
- 如何构建稳的商期货cary合
- 美通企业日报 | 拜耳健康消费品中国研发中心开幕;TCL上半年电视机销量稳居全球第二...
- 【EtherCAT从站(SPI)扩展模块电路分析】
- 启用vsftpd日志及其解读
- 2017年全国大学生电子设计竞赛 单相用电器分析监测装置(k题)
- Win10 Linux 子系统(WSL)监听端口报错Error `IN6_IS_ADDR_V4MAPPED (sin6-sin6_addr.s6_addr32)` Failed的处理
- 【Mo 人工智能技术博客】图卷积网络概述及其在论文分类上的应用
热门文章
- 子进程会继承父进程的哪些内容_【学习教程】Node.js创建子进程方法
- 表单的默认提交方式_对于PHP表单提交有哪集中方式讲解
- foreach判断最后一个_ArrayList集合为什么不能使用foreach增删改?
- Android app 标签,android 获取APP的唯一标识applicationId的实例
- php mysql电商网站设计与制作_赢在电子商务——php+mysql电商网站设计与制作
- c语言的c语句ab2,当初我在重庆交通大学读书的时候C语言成绩也相当不错的
- java对xml解析_Java中对xml的解析
- Java JDK 自带排序(Arrays.sort(数组名))与自行编写的降序
- java 添加等待时间_Java中线程等待特定时间的最有效方法 - java
- 从网络访问计算机不能添加用户名,windows7 和 XP 能互相用\\ip访问,但是不能通过计算机名访问...