word word 1 / 15 word 《计算机仿真》期末大作业 Prim算法和Kruskal算法的Matlab实现 05605禹050697〔30〕 连线问题应用举例: 欲铺设连接个城市的高速公路,假如城与城之间的高速公路造价为,试设计一个线路图,使总的造价最低。 连线问题的数学模型就是图论中在连通的赋权图上求权最小的支撑树。试用Matlab分别实现求最小支撑数的Prim算法和Krusal算法〔避圈法〕。 : 画出程序流程图; 对关键算法、变量和步骤进展解释说明; 用如下两图对所写算法的正确性进展验证。即输入图的信息,输出对应图的最小支撑树与其权值。 〔4〕分析两种算法的实现复杂度。 : 〔1〕提供对算法效率〔复杂度〕进展评估的方法,并通过举例验证,与分析得到的算法复杂度结果相对照; 〔2〕从降低存消耗、减少计算时间的角度,对算法进展优化。 三.实验步骤
i.算法分析与需求分析,程序设计 prim算法的根本思想是:设G=〔V,E〕是一个无向连通网,令T=〔U,TE〕是G的最小生成树。T的初始状态为U={v0}〔v0∈V〕TE={},然后重复执行下述操作:在所有u∈U, v∈V-U的边中找一条代价最小的边〔u,v〕并入集合TE,同时v并入U,直至U=V为止。 显然,Prim算法的根本思想是以局部最优化谋求全局的最优化,而且,还涉与到起始结点的问题。 本程序完成的功能是:从图中的任意结点出发,都能够找出最小生成树 实现方案分析: word word 2 / 15 word Prim算法的关键是如何找到连接U和V-U的最短边来扩大生成树T。设当前T中有k个点〔即T的顶点集U中有k个顶点〕如此所有满足u∈U,v∈V-U的边最多有k×(n-k)条,从如此大的边集中选取最短边是不太经济的。利用MST性质,可以用下述方法构造候选最小边集:对应V-U中的每个顶点,保存从该顶点到U中的各顶点的最短边,取候选边最短边集为V-U中的n-k个顶点所关联的n-k条最短边的集合。为表示候选最短边集,需设置两个一维数组lowcost[n]和adjvex[n],其中lowcost用来保存集合V-U中的各顶点与集合U中顶点的最短边的权值,lowcost[v]=0表示顶点v已经参加最小生成树中;adjvex用来保存依附于该边在集合U中的顶点。例如下式明确顶点vi和顶点vk之间的权值为w lowcost[i]=w; adjvex[i]=k; 程序流程图 关键代码说明: 将用于验证算确性的两幅图用邻接矩阵的形式保存,〔注:矩阵的对角线元素设置为零。〕并在主程序finallyprim中直接进展调用。 在输入起点时应该给程序的阅读者就该图的顶点数作出提示,不然的话使用者很可能会输入越界下标。采取的方法如下 len=length(graph_adjacent);%求图中有多少个顶点 k=sprintf('please input the point where you want to start ,do remember it must be between 1 and %d ',len); start_point=input(k);%输入最小生成树产生起点 采取了sprintf格式化字符串的方法,就防止了编程的人每次根据输入图的顶点 数手动对提示作修改。效果如下: 在对第一幅图进展算法验证时在workspace会如下输出: word word 3 / 15 word please input the point where you want to start ,do remember it must be between 1 and 7 在对第二幅图进展算法验证时在workspace会有如下输出: please input the point where you want to start ,do remember it must be between 1 and 8 在输出结果时为了使结果在workspace中输出的清晰,使人一目了然,也使用了sprintf格式化字符串的方法,代码如下 s=0; for ii=1:len-1 k=sprintf('最小生成树第%d条边:〔%d,%d〕,权值为%d',ii,tree(ii,1),tree(ii,2),graph_adjacent(tree(ii,1),tree(ii,2))); disp(k); disp(' '); s=s+graph_adjacent(tree(ii,1),tree(ii,