Matlab最小生成树.pdfVIP

数学建模——Matlab 中求解最小生成树 关键词：Matlab 、最小生成树、破圈法、避圈法、运筹学 要求：选择一道编程题自己独立完成，必须自己编写源代码，不能从网上 下载。先编写算法的通用程序，然后以例子运行，论文内容包括程序代码、程 序说明、例子运行结果，最终程序文件连同论文一起发至e-mail，便于老师运行 程序是否正确。编程使用C 或MATLAB 。 选择题目：编写实现生成树、最小生成树的程序（包括避圈法、破圈法）。 题目分析：本题要求编写实现生成树、最小生成树的程序，首先来了解一下关于 关于生成树的概念： 1、树的概念： 树是无向图的特殊情况，即对于一个N 个节点的无向图，其中只有N-1 条边， 且图中任意两点间有且只有一条路径，即图中不存在环，这样的图称为树。 2 、生成树的概念： 对于一个无向连通图G=(V,E)，其中V 代表顶点，E 代表边，对它做一次遍历， 每个节点经过一次，那么图中的N 个节点再加上遍历过程中经过的N-1 条边所 构成的子图就是图G 的一个生成树。 3、最小生成树的概念： 对于一个无向连通图 G=(V,E),给它的每条边 （u,v ）赋一个权值w(u,v) 。若图 G 的生成树不止一个，那么其中包含的N-1 条边的权值之和的最小的生成树就是图 G 的最小生成树。 4 、关于运筹学中最小生成树有2 种不错的算法，即避圈法和破圈法，下面来看 一下求解最小生成树的算法： 先看一下图示法表示的最小生成树： 1 1避圈法求解上面无向带权连通图的基本步骤是： 每步从未选的边中选取边e，使它与已选边不构成圈，且e 是未选边中的最小权 边，直到选够n-1 条边为止。 2 2破圈法基本思想如下 : (1) 每次从图中选取任意一个圈, 然后去掉该圈中权值最大的边(如果存在多条 相同权值的最大边，可以任意选择一条去掉即可) 使之不构成圈. (2) 重复上述过程. 直到图中不再含圈且所有顶点均包含在图中为止, 就构成最 小生成树. 它的算法图解如下： 5、我们在用matlab 程序求解最小生成树的时候需要用到无向图的邻接矩阵，首 先来了解一下邻接矩阵的概念。 邻接矩阵（Adjacency Matrix ）：是表示顶点之间相邻关系的矩阵。设G=(V,E)是 3 一个图，其中V={v1,v2, …,vn} 。G 的邻接矩阵是一个具有下列性质的n 阶方阵： ①对无向图而言，邻接矩阵一定是对称的，而且对角线一定为零，有向图则不 一定如此。 ②在无向图中，任一顶点i 的度为第i 列所有元素的和，在有向图中顶点i 的出 度为第i 行所有元素的和，而入度为第i 列所有元素的和。 ③用邻接矩阵法表示图共需要n^2 个空间，由于无向图的邻接矩阵一定具有对 称关系，所以扣除对角线为零外，仅需要存储上三角形或下三角形的数据即可， 因此仅需要n （n-1）/2 个空间。 看一个例子： 上面的邻接矩阵是带权邻接矩阵，有两种表示方法，我们往往只使用第一 种，它的正对角线上全部是 0，方便在程序中求解。 编程内容： 有了以上的数学基础，就可以在Matlab 中编写适当的程序来求解无向带权 图的最小生成树。下面先用避圈法实现求解，代码中有详细的解析说明： 4 1避圈法Matlab 程序代码： %编程工具Matlab； %这是一个通过避圈法求解连通带权图的最小生成树的程序. n=input(请输入图的顶点数目:n= ) W=input( 请 输 入 图 的 加 权 邻 接 矩 阵:[W(1,1),..,W(1,n);..;W(n,1),..,W(n,n)]=) %用W(i,i)=inf 代替 =0 %准备工作 T=zeros(n); %最小生成树的加权邻接矩阵 WW=W; for i=1:n for j=1:n if W(i,j)==inf WW(i,j)=0; end