分配问题——匹配.pptVIP

* * 主要内容 基本概念 Kuhn-munkras算法的MATLAB程序 自编函数maxmatch的用法 应用：引例的求解 引例：车辆调度问题 设某车队有8辆车，存放在不同的地点，队长要派出其中5辆到5个工地去运货。各车从存放处调到装货地点所需费用列于下页表，问应选哪5辆车调到何处去运货，才能使各车从车所在地点调到装货地点所需的总费用最少？ 引例：车辆调度问题 18 16 14 16 17 18 20 21 5 21 18 19 25 24 19 30 29 4 26 23 22 19 19 30 29 28 3 19 30 20 21 18 19 31 29 2 26 22 19 27 32 18 25 30 1 8 7 6 5 4 3 2 1 车 地点 引例：车辆调度问题 18 16 14 16 17 18 20 21 5 21 18 19 25 24 19 30 29 4 26 23 22 19 19 30 29 28 3 19 30 20 21 18 19 31 29 2 26 22 19 27 32 18 25 30 1 8 7 6 5 4 3 2 1 车 地点 x5 y1 y2 y3 y4 y5 x1 x2 x3 x4 y6 y7 y8 基本概念 ? 匹配：图G=(V,E)中一些互不相邻的边构成的集合M，称为G的一个匹配；M中边的端点称为M渗透点；V中的其余顶点称为M非渗透点。 ? 理想匹配：若匹配M渗透了图G的所有顶点，则称M为G的理想匹配。 ? 最大匹配：若匹配M是G中边数最多的匹配，则称M为G的最大匹配。 ? M可增长路径：设M是图G=(V,E)的一个匹配，P是G中的一条路径。若P的边在M与E-M中交替出现，则称P是一条M交错路径；若一条M交错路径的起点和终点都是M非渗透点，则称其为M可增长路径。 基本概念 ? 最大权匹配：在加权二部图G中，边权和最大的匹配M称为G的最大权匹配。 ? 二部图的边矩阵A ：A=（aij) mxn ? 二部图的边权矩阵A ：A=（aij)mxn , aij为边（xi，yj)的权，若xi与yj之间无边，则aij取0. MATLAB程序——Kuhn-munkras算法 function sumw=kuhngong(A) n=size(A,1); w=A; l=zeros(n,2); for i=1:n for j=1:n if l(i,1)w(i,j) l(i,1)=w(i,j); end end end FLAG_AGL=zeros(n,n); FLAG_S=zeros(1,n); FLAG_T=zeros(1,n); FLAG_NGLS=zeros(1,n);f=zeros(n,2); for i=1:n for j=1:n if l(i,1)+l(j,2)==w(i,j) FLAG_AGL(i,j)=i; end end end M=zeros(n,2); for i=1:n for j=1:n if (FLAG_AGL(i,j)==i)(~M(j,2)) (~M(i,1)) M(i,1)=i; M(j,2)=i; end end end FLAG3=1; while FLAG3 FLAG3=0; u=0; for i=1:n if ~M(i,1) u=i; break; end end end MATLAB程序——Kuhn-munkras算法（续） if ~u fprintf(1,‘二部图最大权匹配运行结果\n); fprintf(1,‘\n\n求得最大权匹配M={); sumw=0; for i=1:n for j=1:n if M(j,2)==i fprintf(1,x%dy%d,,i,j); sumw=sumw+w(i,j); break; end end end fprintf(1,}\n); fprintf(1,‘最大权W(M)=%g\n,sumw); return else FLAG_S=zeros(1,n); FLAG_T=zeros(1,n); FLAG_S(u)=1;