采用A算法解决八数码问题.docx

采用A算法解决八数码问题.docx

《采用A算法解决八数码问题.docx》由会员分享，可在线阅读，更多相关《采用A算法解决八数码问题.docx（29页珍藏版）》请在冰豆网上搜索。

采用A算法解决八数码问题

人工智能实验一报告

题目：

采用A*算法解决八数码问题

**********

学号：

*********

专业：

信息与计算科学

提交日期：

2014-05-22

1问题描述

所谓八数码问题是指这样一种游戏：

将分别标有数字1，2，3，…，8的八块正方形数码牌任意地放在一块3×3的数码盘上。

放牌时要求不能重叠。

于是，在3×3的数码盘上出现了一个空格。

现在要求按照每次只能将与空格相邻的数码牌与空格交换的原则，不断移动该空格方块以使其和相邻的方块互换，直至达到所定义的目标状态。

空格方块在中间位置时有上、下、左、右4个方向可移动，在四个角落上有2个方向可移动，在其他位置上有3个方向可移动，问题描述如下图1-1所示：

初始状态中间状态目标状态

图1-1八数码问题的求解过程

1.1待解决问题的解释

首先，八数码问题包括一个初始状态（START）和目标状态（END），所谓解决八数码问题就是在两个状态间寻找一系列可过渡状态：

（START>STATE1>STATE2>...>END）

这个状态是否存在就是我们要解决的第一个问题：

Q1：

每一个状态及每一次操作的表示方法？

有许多表示方法，比如一个3*3的八数码盘可以压缩成一个int值表示，但不适用于15puzzle或大于8的puzzle问题。

如果对空间要求很高，应该还可以再压缩。

本文采用一个int表示的方法。

表示方法如下：

由于int的表示范围大于1e9，所以我们取一个int的低9位，为了方便寻找空格的位置，int的个位我们用来放空格的位置（19）。

而前8位，按照行从上到下，列从左到右的顺序依次记录对应位置上的数字。

1.2问题的搜索形式描述

八数码问题形式化描述：

初始状态：

初始状态向量：

规定向量中各分量对应的位置，各位置上的数字。

把3×3的棋盘按从左到右，从上到下的顺序写成一个一维向量。

我们可以设定初始状态：

<1,5,2,4,0,3,6,7,8>

后继函数：

按照某种规则移动数字得到的新向量。

例如：

<1,5,2,4,0,3,6,7,8><1,0,2,4,5,3,6,7,8>

目标测试：

新向量是都是目标状态。

即<1,2,3,4,5,6,7,8,0>是目标状态？

路径耗散函数：

每次移动代价为1，每执行一条规则后总代价加1。

1.3解决方案介绍（原理）

该问题是一个搜索问题。

它是一种状态到另一种状态的变换。

要解决这个问题，必须先把问题转化为数字描述。

由于八数码是一个3*3的矩阵，但在算法中不实用矩阵，而是将这个矩阵转化为一个一维数组，使用这个一维数组来表示八数码，但是移动时要遵守相关规则。

（1）可用如下形式的规则来表示数字通过空格进行移动：

→

（2）共24条移动规则，对应与每个位置的移动规则。

（3）搜索顺序举例：

<1>优先移动行数小的棋子（数字）

<2>同一行中优先移动列数大的棋子

（4）约束规则：

不使离开既定位置的数字数增加

八数码的节点扩展应当遵循棋子的移动规则。

按规则，每一次可以将一个与空格相邻的棋子移动到空格中，实际上也可以看做空格的相反方向移动。

空格的移动方向可以是上下左右，当然不能出边界。

棋子的位置，也就是保存状态的数组元素的下标，空格移动后，相应位置发生变化，在不移出边界的条件下，空格向右，下，左，上移动后，新位置是原位置分别加上1,3，-1，-3。

在这里，空格可以用任意数字表示。

操作本文用urdl分别表示空格的向上向右向下向左四个操作。

图的搜索策略：

经分析，8数码问题的搜索策略共有：

1.广度优先搜索、2.深度优先搜索、3.有界深度优先搜索、4.最好优先搜索、5.局部择优搜索，等等。

其中广度优先搜索法是可采纳的，有界深度优先搜索法是不完备的，最好优先和局部择优搜索法是启发式搜索法。

本实验采用启发式A*搜索算法来实现。

2算法介绍

问题的求解实际上就是在这个图中找到一条路径可以从开始到结果。

这个寻找的过程就是状态空间搜索。

常用的状态空间搜索有深度优先和广度优先。

广度优先是从初始状态一层一层向下找，直到找到目标为止。

深度优先是按照一定的顺序前查找完一个分支，再查找另一个分支，以至找到目标为止。

启发式搜索就是在状态空间中的搜索对每一个搜索的位置进行评估，得到最好的位置，再从这个位置进行搜索直到目标。

这样可以省略大量无畏的搜索路径，提高了效率。

2.1A*搜索算法一般介绍

A*算法实际是一种启发式搜索，所谓启发式搜索，就是利用一个估价函数评估每次的的决策的价值，决定先尝试哪一种方案，这样可以极大的优化普通的广度优先搜索。

一般来说，从出发点（A）到目的地（B）的最短距离是固定的，我们可以写一个函数judge（）估计A到B的最短距离，如果程序已经尝试着从出发点A沿着某条路线移动到了C点,那么我们认为这个方案的AB间的估计距离为A到C实际已经行走了的距离H加上用judge（）估计出的C到B的距离。

如此，无论我们的程序搜索展开到哪一步，都会算出一个评估值，每一次决策后，将评估值和等待处理的方案一起排序，然后挑出待处理的各个方案中最有可能是最短路线的一部分的方案展开到下一步，一直循环到对象移动到目的地，或所有方案都尝试过却没有找到一条通向目的地的路径则结束。

A*算法是一个可采纳的最好优先算法。

A*算法的估价函数可表示为：

f'（n）=g'（n）+h'（n）

这里，f'（n）是估价函数，g'（n）是起点到终点的最短路径值，h'（n）是n到目标的最断路经的启发值。

由于这个f'（n）其实是无法预先知道的，所以我们用前面的估价函数f（n）做近似。

g（n）代替g'（n），但g（n）>=g'（n）才可（大多数情况下都是满足的，可以不用考虑），h（n）代替h'（n），但h（n）<=h'（n）才可。

可以证明应用这样的估价函数是可以找到最短路径的，也就是可采纳的。

2.2算法伪代码

首先定义两个表，open表用于存放已经生成，且已用启发式函数进行过估计或评价，但尚未产生它们的后继节点的那些结点，这些结点也称未考察结点；

closed表用于存放已经生成，且已考察过的结点。

设S0为初态，Sg为目标状态。

具体过程如下：

（1）把S0放入open表，记为f=h，令closed为空表；

（2）重复下列过程，直至找到目标结点为止。

若open为空表，则失败；

（3）选取open表中未设置过的具有最小f值的结点为最佳节点，并放入closed表中

（4）若最佳节点不是目标节点，则扩展之，产生后继节点。

（5）对每个后继结点进行下列过程：

Ø建立从该后继结点返回最佳节点的指针；

Ø计算g（后继结点）=g（最佳节点）+k（最佳节点，后继结点）；

ØSs：

如果该后继节点∈open，则称此节点为old，并把它添加至最佳节点的后继节点中

Ø比较新旧路径代价，如果个g（后继节点）

Ø若至old节点的代价比较低或一样，则停止扩展节点

Ø若后继节点不再open表中，则看其是否在closed中

Ø若后继节点在open表中，则转向Ss；

Ø若后继节点既不在open表中，又不在closed表中，则把它放入open表中，并填入最佳节点的后裔表，然后走下一步

（6）计算f值

（7）GOLOOP

3算法实现

3.1实验环境与问题规模

（1）实验环境：

WindowsXP

（2）实验编程工具：

VC++6.0

3.2数据结构

1）定义状态图中的结点数据结构

typedefstructNode

{

statusdata;//结点所存储的状态

structNode*parent;//指向结点的父亲结点

structSpringLink*child;//指向结点的后继结点

structNode*next;//指向open或者closed表中的后一个结点

intfvalue;//结点的总的路径

intgvalue;//结点的实际路径

inthvalue;//结点的到达目标的苦难程度

}NNode,*PNode;

2）定义存储指向结点后继结点的指针的地址

typedefstructSpringLink

{

structNode*pointData;//指向结点的指针

structSpringLink*next;//指向兄第结点

}SPLink,*PSPLink;

3）初始化一个空链表

voidinitLink（PNode&Head）

3.3实验结果

经测试，程序运行良好，结果正确。

输入测试数据，初试状态<0,1,2,3,4,5,6,7，8>，目标状态<,1,4,2,3,5,8,6,7,0>，运行结果有解，共经过四步。

3.4系统中间及最终输出结果

图3-1程序运行中间结果

图3-2程序运行最终结果

4参考文献

1．《Artificialintelligence:

;amodernapproach人工智能:

一种现代方法》作者：

Russell,StuartJ.出版社：

清华大学出版社

2．CSDN博客

5附录—源代码及其注释

#include"iostream"

#include"stdlib.h"

#include"conio.h"

#definesize3

usingnamespacestd;

//定义二维数组来存储数据表示某一个特定状态

typedefintstatus[size][size];

structSpringLink;

//定义状态图中的结点数据结构

typedefstructNode

{

statusdata;//结点所存储的状态

structNode*parent;//指向结点的父亲结点

structSpringLink*child;//指向结点的后继结点

structNode*next;//指向open或者closed表中的后一个结点

intfvalue;//结点的总的路径

intgvalue;//结点的实际路径

inthvalue;//结点的到达目标的苦难程度

}NNode,*PNode;

//定义存储指向结点后继结点的指针的地址

typedefstructSpringLink

{

structNode*pointData;//指向结点的指针

structSpringLink*next;//指向兄第结点

}SPLink,*PSPLink;

PNodeopen;

PNodeclosed;

//开始状态与目标状态

statusstartt={0,1,2,3,4,5,6,7,8};

statustarget={1,4,2,3,5,8,6,7,0};

//初始化一个空链表

voidinitLink（PNode&Head）

{

Head=（PNode）malloc（sizeof（NNode））;

Head->next=NULL;

}

//判断链表是否为空

boolisEmpty（PNodeHead）

{

if（Head->next==NULL）

returntrue;

else

returnfalse;

}

//从链表中拿出一个数据

voidpopNode（PNode&Head,PNode&FNode）

{

if（isEmpty（Head））

{

FNode=NULL;

return;

}

FNode=Head->next;

Head->next=Head->next->next;

FNode->next=NULL;

}

//向结点的最终后继结点链表中添加新的子结点

voidaddSpringNode（PNode&Head,PNodenewData）

{

PSPLinknewNode=（PSPLink）malloc（sizeof（SPLink））;

newNode->pointData=newData;

newNode->next=Head->child;

Head->child=newNode;

}

//释放状态图中存放结点后继结点地址的空间

voidfreeSpringLink（PSPLink&Head）

{

PSPLinktmm;

while（Head!

=NULL）

{

tmm=Head;

Head=Head->next;

free（tmm）;

}

}

//释放open表与closed表中的资源

voidfreeLink（PNode&Head）

{

PNodetmn;

tmn=Head;

Head=Head->next;

free（tmn）;

while（Head!

=NULL）

{

//首先释放存放结点后继结点地址的空间

freeSpringLink（Head->child）;

tmn=Head;

Head=Head->next;

free（tmn）;

}

}

//向普通链表中添加一个结点

voidaddNode（PNode&Head,PNode&newNode）

{

newNode->next=Head->next;

Head->next=newNode;

}

//向非递减排列的链表中添加一个结点

voidaddAscNode（PNode&Head,PNode&newNode）

{

PNodeP;

PNodeQ;

P=Head->next;

Q=Head;

while（P!

=NULL&&P->fvaluefvalue）

{

Q=P;

P=P->next;

}

//上面判断好位置之后，下面就是简单的插入了

newNode->next=Q->next;

Q->next=newNode;

}

//计算结点额h值

intcomputeHValue（PNodetheNode）

{

intnum=0;

for（inti=0;i<3;i++）

{

for（intj=0;j<3;j++）

{

if（theNode->data[i][j]!

=target[i][j]）

num++;

}

}

returnnum;

}

//计算结点的f，g，h值

voidcomputeAllValue（PNode&theNode,PNodeparentNode）

{

if（parentNode==NULL）

theNode->gvalue=0;

else

theNode->gvalue=parentNode->gvalue+1;

theNode->hvalue=computeHValue（theNode）;

theNode->fvalue=theNode->gvalue+theNode->hvalue;

}

//初始化函数，进行算法初始条件的设置

voidinitial（）

{

//初始化open以及closed表

initLink（open）;

initLink（closed）;

//初始化起始结点，令初始结点的父节点为空结点

PNodeNULLNode=NULL;

PNodeStart=（PNode）malloc（sizeof（NNode））;

for（inti=0;i<3;i++）

{

for（intj=0;j<3;j++）

{

Start->data[i][j]=startt[i][j];

}

}

Start->parent=NULL;

Start->child=NULL;

Start->next=NULL;

computeAllValue（Start,NULLNode）;

//起始结点进入open表

addAscNode（open,Start）;

}

//将B节点的状态赋值给A结点

voidstatusAEB（PNode&ANode,PNodeBNode）

{

for（inti=0;i<3;i++）

{

for（intj=0;j<3;j++）

{

ANode->data[i][j]=BNode->data[i][j];

}

}

}

//两个结点是否有相同的状态

boolhasSameStatus（PNodeANode,PNodeBNode）

{

for（inti=0;i<3;i++）

{

for（intj=0;j<3;j++）

{

if（ANode->data[i][j]!

=BNode->data[i][j]）

returnfalse;

}

}

returntrue;

}

//结点与其祖先结点是否有相同的状态

boolhasAnceSameStatus（PNodeOrigiNode,PNodeAnceNode）

{

while（AnceNode!

=NULL）

{

if（hasSameStatus（OrigiNode,AnceNode））

returntrue;

AnceNode=AnceNode->parent;

}

returnfalse;

}

//取得方格中空的格子的位置

voidgetPosition（PNodetheNode,int&row,int&col）

{

for（inti=0;i<3;i++）

{

for（intj=0;j<3;j++）

{

if（theNode->data[i][j]==0）

{

row=i;

col=j;

return;

}

}

}

}

//交换两个数字的值

voidchangeAB（int&A,int&B）

{

intC;

C=B;

B=A;

A=C;

}

//检查相应的状态是否在某一个链表中

boolinLink（PNodespciNode,PNodetheLink,PNode&theNodeLink,PNode&preNode）

{

preNode=theLink;

theLink=theLink->next;

while（theLink!

=NULL）

{

if（hasSameStatus（spciNode,theLink））

{

theNodeLink=theLink;

returntrue;

}

preNode=theLink;

theLink=theLink->next;

}

returnfalse;

}

//产生结点的后继结点（与祖先状态不同）链表

voidSpringLink（PNodetheNode,PNode&spring）

{

introw;

intcol;

getPosition（theNode,row,col）;

//空的格子右边的格子向左移动

if（col!

=2）

{

PNoderlNewNode=（PNode）malloc（sizeof（NNode））;

statusAEB（rlNewNode,theNode）;

changeAB（rlNewNode->data[row][col],rlNewNode->data[row][col+1]）;

if（hasAnceSameStatus（rlNewNode,theNode->parent））

{

free（rlNewNode）;//与父辈相同，丢弃本结点

}

else

{

rlNewNode->parent=theNode;

rlNewNode->child=NULL;

rlNewNode->next=NULL;

computeAllValue（rlNewNode,theNode）;

//将本结点加入后继结点链表

addNode（spring,rlNewNode）;

}

}

//空的格子左边的格子向右移动

if（col!

=0）

{

PNodelrNewNode=（PNode）malloc（sizeof（NNode））;

statusAEB（lrNewNode,theNode）;

changeAB（lrNewNode->data[row][col],lrNewNode->data[row][col-1]）;

if（hasAnceSameStatus（lrNewNode,theNode->parent））

{

free（lrNewNode）;//与父辈相同，丢弃本结点

}

else

{

lrNewNode->parent=theNode;

lrNewNode->child=NULL;

lrNewNode->next=NULL;

computeAllValue（lrNewNode,theNode）;

//将本结点加入后继结点链表

addNode（spring,lrNewNode）;

}

}

//空的格子上边的格子向下移动

if（row!

=0）