匈牙利算法.docVIP

在下将《运筹学》上第五章第五节的四步求解最优指派 的过程用Matlab实现。共编写了三个子程序。请各位高 手指正。 fc01是主子程序，在其中调用其余两个子程序 输入参数flag=1时,函数求解极小化问题, flag=0或不输入时,函数求解极大化问题 fc02是寻找最大对集的子程序，实现书上第二步。 上次贴出来过，但其中有个隐蔽小错误，现 已改正，上次贴出来的地方也改正了 fc03实现第三步和第四步 function [e,s]=fc01(a,flag) if nargin==1 flag=0; end b=a; if flag==0 cmax=max(max(b)); b=cmax-b; end m=size(b); for i=1:m(1) b(i,:)=b(i,:)-min(b(i,:)); end for j=1:m(1) b(:,j)=b(:,j)-min(b(:,j)); end d=(b==0); [e,total]=fc02(d); while total~=m(1) b=fc03(b,e); d=(b==0); [e,total]=fc02(d); end inx=sub2ind(size(a),e(:,1),e(:,2)); e=[e,a(inx)]; s=sum(a(inx)); function [e,total]=fc02(d) total=0; m=size(d); e=zeros(m(1),2); t=sum(sum(d)); nump=sum(d); while t~=0 [s,inp]=sort(nump); inq=find(s); ep=inp(inq(1)); inp=find(d(ep,:)); numq=sum(d(:,inp)); [s,inq]=sort(numq); eq=inp(inq(1)); total=total+1; e(total,:)=[ep,eq]; inp=find(d(:,eq)); nump(inp)=nump(inp)-1; nump(ep)=0; t=t-sum(d(ep,:))-sum(d(:,eq))+1; d(ep,:)=0*d(ep,:); d(:,eq)=0*d(:,eq); end function b=fc03(b,e) m=size(b); t=1; p=ones (m(1),1); q=zeros(m(1),1); inp=find(e(:,1)~=0); p(e(inp,1))=0; while t~=0 tp=sum(p+q); inp=find(p==1); n=size(inp); for i=1:n(1) inq=find(b(inp(i),:)==0); q(inq)=1; end inp=find(q==1); n=size(inp); for i=1:n(1) if all(e(:,2)-inp(i))==0 inq=find( ( e(:,2)-inp(i) )==0 ); p(e(inq))=1; end end tq=sum(p+q); t=tq-tp; end inp=find(p==1); inq=find(q==0); cmin=min(min(b(inp,inq))); inq=find(q==1); b(inp,:)=b(inp,:)-cmin; b(:,inq)=b(:,inq)+cmin; 程序较短，不用加注释。而且我编的是“黑箱”程序， 给出输入输出参数的意义即可使用。 主程序只需调用fc01即可完成任务 fc01的输入参数 a 表示你所求问题的加权矩阵 flag控制求解极大化还是极小化。flag=1时, 函数 求解极小化问题,flag=0或不输入时,函数求解 极大化问题 fc01的输出参数 e是所求最优对集的表示，第一列表示行号，第二列表示列号 第三列表示其对应元素的值 s是所求的最优指派的最优值 举个例子，说明使用方法 例如《运筹学》课本P131例8的求解， 因为是求解极小化问题，flag=1