匈牙利算法的MATLAB代码

精品文档程序文件 fenpei.mfunction z,ans=fenpei(marix)%/ %输入效率矩阵 marix 为方阵; %若效率矩阵中有 M,则用一充分大的数代替; %输出z为最优解，ans为 最优分配矩阵；%/a=marix;b=a;%确定矩阵维数s=length(a);%确定矩阵行最小值，进行行减ml=min(a);for i=1:s a(i,:)=a(i,:)-ml(i);end%确定矩阵列最小值，进行列减mr=min(a);for j=1:s a(:,j)=a(:,j)-mr(j);end% start workingnum=0;while(num=s) %终止条件是“（0）”的个数与矩阵的维数相同 %index用以标记矩阵中的零元素，若a(i,j)=0，则index(i,j)=1,否则index(i,j)=0 index=ones(s); index=a&index; index=index; %flag用以标记划线位，flag=0 表示未被划线， %flag=1 表示有划线过，flag=2 表示为两直线交点 %ans用以记录 a 中“（0）”的位置 %循环后重新初始化flag,ans flag = zeros(s); ans = zeros(s); %一次循环划线全过程，终止条件是所有的零元素均被直线覆盖， %即在flag0位,index=0 while(sum(sum(index) %按行找出“（0）”所在位置，并对“（0）”所在列划线， %即设置flag,同时修改index,将结果填入ans for i=1:s t=0; l=0; for j=1:s if(flag(i,j)=0&index(i,j)=1) l=l+1; t=j; end end if(l=1) flag(:,t)=flag(:,t)+1; index(:,t)=0; ans(i,t)=1; end end %按列找出“（0）”所在位置，并对“（0）”所在行划线， %即设置flag,同时修改index,将结果填入ans for j=1:s t=0; r=0; for i=1:s if(flag(i,j)=0&index(i,j)=1) r=r+1; t=i; end end if(r=1) flag(t,:)=flag(t,:)+1; index(t,:)=0; ans(t,j)=1; end end end %对 while(sum(sum(index) %处理过程 %计数器：计算ans中1的个数，用num表示 num=sum(sum(ans); % 判断是否可以终止，若可以则跳出循环 if(s=num) break; end %否则，进行下一步处理 %确定未被划线的最小元素，用m表示 m=max(max(a); for i=1:s for j=1:s if(flag(i,j)=0) if(a(i,j) a=37.7 32.9 38.8 37 35.443.4 33.1 42.2 34.7 41.833.3 28.5 38.9 30.4 33.629.2 26.4 29.6 28.5 31.10 0 0 0 0; z,ans=fenpei(a)z = 127.8000ans = 0 0 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 0 0 1 0 0 1 1 1 2 2 2 2 0 2 4 4 3 3 1 1 0 7 2 2 4 4 7 7 0 3 5 5 7 7 5 5 0 4 3 3 9 9 6 6 0 8 8 8 1 1 4 4 0 6 6 6 12 12 11 1