匈牙利解法的matlab算法,匈牙利算法的matlab编程求助

function z=Hungarian(A,m,n);

%求最大匹配

M(m,n)=0;

for(i=1:m)

for(j=1:n)

if(A(i,j)) M(i,j)=1;break;

end;

end %求初始匹配M

if(M(i,j))break;

end;

end %获得仅含一条边的初始匹配M

while(1)

for(i=1:m)

x(i)=0;

end %将记录X中点的标号和标记*

for(i=1:n)

y(i)=0;

end %将记录Y中点的标号和标记*

for(i=1:m)

pd=1; %寻找X中M的所有非饱和点

for(j=1:n)

if(M(i,j))pd=0;

end;

end

if(pd)x(i)=-n-1;

end;

end %将X中M的所有非饱和点都给以标号0 和标记*, 程序中用n+1 表

%示0 标号, 标号为负数时表示标记*

pd=0;

while(1)xi=0;

for(i=1:m)

if(x(i)<0)xi=i;break;

end;

end %假如X 中存在一个既有标号又有标记*的点, 则任取X中一个既有标号又有标记*的点xi

if(xi==0)pd=1;break;

end %假如X中所有有标号的点都已去掉了标记*, 算法终止

x(xi)=x(xi)*(-1); %去掉xi 的标记*

k=1;

for(j=1:n)

if(A(xi,j)&y(j)==0)y(j)=xi;yy(k)=j;k=k+1;

end;

end %对与xi 邻接且尚未给标号的yj 都给以标号i

if(k>1)k=k-1;

for(j=1:k)

pdd=1;

for(i=1:m)

if(M(i,yy(j)))x(i)=-yy(j);pdd=0;break;

end;

end %将yj 在M中与之邻接的点xk (即xkyj∈M), 给以标号j 和标记*

if(pdd)break;

end;

end

if(pdd)k=1;j=yy(j); %yj 不是M的饱和点

while(1)P(k,2)=j;P(k,1)=y(j);j=abs(x(y(j))); %任取M的一个非饱和点yj, 逆向返回

if(j==n+1)break;

end %找到X中标号为0 的点时结束, 获得M-增广路P

k=k+1;

end

for(i=1:k)

if(M(P(i,1),P(i,2)))M(P(i,1),P(i,2))=0; %将匹配M 在增广路P 中出现的边去掉

else M(P(i,1),P(i,2))=1;end;end %将增广路P 中没有在匹配M中出现的边加入到匹配M中

break;end;end;end

if(pd)break;end;end %假如X中所有有标号的点都已去掉了标记*, 算法终止M %显示最大匹配M, 程序结束

z=M