数学建模提高班第六讲-网络优化模型及案例分析.pptVIP

网络优化模型及案例分析 本专题学习目的 掌握把实际问题转化为图或网络问题的方法 了解图的基本概念和矩阵表示方法 掌握最短路问题、最小生成树问题的算法 了解旅行商问题 最短路问题及算法 1) 赋权图中从给定点到其余顶点的最短路 2) 求赋权图中任意两顶点间的最短路 算法的基本思想 算法的基本思想 （I）求距离矩阵的方法. （II）求路径矩阵的方法. （III）查找最短路路径的方法. （IV）Floyd算法：求任意两顶点间的最短路. 最小生成树及算法 许多应用问题都是一个求无向连通图的最小生成树问题。例如：要在n个城市之间铺设光缆，主要目标是要使这 n 个城市的任意两个之间都可以通信，但铺设光缆的费用很高，且各个城市之间铺设光缆的费用不同；另一个目标是要使铺设光缆的总费用最低。这就需要找到带权的最小生成树。 1) 树的定义与树的特征 定理2 设G是具有n个顶点的图，则下述命题等价： 2）图的生成树 （II）找图中生成树的方法 A 避圈法 B 破圈法 B 破圈法 3) 最小生成树与算法 A Kruskal算法（或避圈法） B破圈法 旅行售货员问题 旅行售货员问题或货郎担问题 一个可行的办法 : 例9 对下图,用二边逐次修正法求较优H圈. 附录：图论算法的matlab实现 最短路问题 例 某公司在六个城市中有分公司，从ci到cj的直接航程票价记在下述矩阵的位置上。（表示无直接航路），请帮助该公司设计一张城市c1到其它城市间的票价最便宜的路线图。 用矩阵 （为顶点个数）存放各边权的邻接矩阵，行向量pb、index1、index2、d分别用来存放标号信息、标号顶点顺序、标号顶点索引、最短通路的值。其中分量 index2(i)存放始点到第点最短通路中第顶点前一顶点的序号； d(i)存放由始点到第点最短通路的值。 Dijkstra 的Matlab程序如下： M=10000; a(1,:)=[0,50,M,40,25,10]; a(2,:)=[zeros(1,2),15,20,M,25]; a(3,:)=[zeros(1,3),10,20,M]; a(4,:)=[zeros(1,4),10,25]; a(5,:)=[zeros(1,5),55]; a(6,:)=zeros(1,6); a=a+a; pb(1:length(a))=0; pb(1)=1; index1=1; index2=ones(1,length(a)); d(1:length(a))=M;d(1)=0;temp=1; while sum(pb)length(a) tb=find(pb==0); d(tb)=min(d(tb),d(temp)+a(temp,tb)); tmpb=find(d(tb)==min(d(tb))); temp=tb(tmpb(1)); pb(temp)=1; index1=[index1,temp]; index=index1(find(d(index1)==d(temp)-a(temp,index1))); if length(index)=2 index=index(1); end index2(temp)=index; end d, index1, index2 Floyd算法的Matlab程序如下： clear; clc; M=10000; a(1,:)=[0,50,M,40,25,10]; a(2,:)=[zeros(1,2),15,20,M,25]; a(3,:)=[zeros(1,3),10,20,M]; a(4,:)=[zeros(1,4),10,25]; a(5,:)=[zeros(1,5),55]; a(6,:)=zeros(1,6); b=a+a;path=zeros(length(b)); for k=1:6 for i=1:6 for j=1:6 if b(i,j)b(i,k)+b(k,j) b(i,j)=b(i,k)+b(k,j); path(i,j)=k; end end end end b, path 最小生成树问题 prim算法如下 ： clc;clear; M=1000; a(1,2)=50; a(1,3)=60; a(2,4)=65; a(2,5)=40; a(3,4)=52;a(3,7)=45; a(4,5)=50; a(4,6)=30;a(4,7)=42; a(5,6)=70; a=[a;zeros(2,7)]; a=a+a;a(find(a==0))=M; result=[];p=1;tb=2:leng