【正文】 續(xù)上表杭州寧波呼和浩特鄂爾多斯附錄13A,B城市配對的MATLAB程序：function [z,ans]=fenpei(marix)%////////////////////////////////////////////////// %輸入效率矩陣 marix 為方陣。 %輸出z為最優(yōu)解，ans為 最優(yōu)分配矩陣；%//////////////////////////////////////////////////a=marix。%確定矩陣維數(shù)s=length(a)。)。end%確定矩陣列最小值，進行列減mr=min(a)。end% start workingnum=0。 index=aamp。 index=~index。 ans = zeros(s)。 l=0。amp。 t=j。 index(:,t)=0。 end end %按列找出“（0）”所在位置，并對“（0）”所在行劃線， %即設置flag,同時修改index,將結果填入ans for j=1:s t=0。 for i=1:s if(flag(i,j)==0amp。index(i,j)==1) r=r+1。 end end if(r==1) flag(t,:)=flag(t,:)+1。 ans(t,j)=1。 % 判斷是否可以終止，若可以則跳出循環(huán) if(s==num) break。 for i=1:s for j=1:s if(flag(i,j)==0) if(a(i,j)m) m=a(i,j)。 end if(flag(i,j)==2) a(i,j)=a(i,j)+m。z=0。程序運行返回的01矩陣：ans = Columns 1 through 19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0Column 20 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0附錄14a=[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0