匈牙利算法 matlab程序出现问题 运行时出现 Undefined function or variable 'a'.

调用 Edmonds 时，忘了个 实际参虚判数 a（有值得）或数值了。

也可以 用默认参轮誉做数值。腊衡

function [Matching,Cost] = Edmonds(a)

if nargin==0

a=.......

end

Matching = zeros(size(a))

.........

是不是你的程团李序掉进死循环了

你检查下

输入的新的数据不符合程序允许,或者迟滑超过极码或腊限值

寻匈牙利算法:

function [result,msum]=sbppp(cost,m)

if nargin==1

dd=cost

else

dd=max(max(cost))-cost

end

[nop,nop]=size(cost)msum=0

for i=1:nop

dd(i,:)=dd(i,:)-min(dd(i,:))

end

for j=1:nop

dd(:,j)=dd(:,j)-min(dd(:,j))

end

backup=dd

for z=1:nop

bh=nopbl=nopresult=[]

for k=1:nop

for i=1:nop

h=find(dd(i,:)==0)

if length(h)~=0&length(h)<bh

bh=length(h)

ch=i

end

end

L=find(dd(ch,:)==0)

for j=1:length(L)

l=find(dd(:,L(j))==0)

if length(l)<bl

bl=length(l)

cl=L(j)

end

end

result(1,k)=chresult(2,k)=cl

dd(ch,:)=1dd(:,cl)=1

bl=nopbh=nop

if length(find(dd==0))==0

break

end

end

if length(result(1,:))==nop

break

end

dd=backupDD=ddd=zeros(nop)

for i=1:length(result(1,:))

d(result(1,i),result(2,i))=1

end

D=~(d+dd)

p=[]q=[]k=1zx=inf

for i=1:nop

if sum(d(i,:))==0

p(k)=i

k=k+1

end

end

for j=1:length(p)

q=find(D(p(j),:)==1)

for e=1:length(q)

pp=find(d(:,q(e))==1)

if pp~=0

p(k)=pp

k=k+1

end

end

end

for l=1:length(p)

q=find(D(p(l),:)==1)

for u=1:length(q)

DD(:,q(u))=inf

end

end

for l=1:length(p)

if min(DD(p(l),:))<zx

zx=min(DD(p(l),:))

end

end

for l=1:length(p)

q=find(D(p(l),:)==1)

for u=1:length(q)

dd(:,q(u))=dd(:,q(u))+zx

end

end

for l=1:length(p)

dd(p(l),:)=dd(p(l),:)-zx

end

backup=dd

end

for i=1:length(result(1,:))

msum=msum+cost(result(1,i),result(2,i))

end

匈牙利算法的MatLab实现 收藏

程序文件 fenpei.m

//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

function [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 working

num=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)

%一次循环划线全过程，终止条件是所有的零元素均被直线覆盖，

%即在flag>0位,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)<m)

m=a(i,j)

end

end

end

end

%未被划线，即flag=0处减去m；线交点，即flag=2处加上m

for i=1:s

for j=1:s

if(flag(i,j)==0)

a(i,j)=a(i,j)-m

end

if(flag(i,j)==2)

a(i,j)=a(i,j)+m

end

end

end

end %对while(num~=s)

%计算最优（min）值

zm=ans.*b

z=0

z=sum(sum(zm))

/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

运行实例：

>>a=[37.732.938.83735.4

43.433.142.234.741.8

33.328.538.930.433.6

29.226.429.628.531.1

00000]

>>[z,ans]=fenpei(a)

z =

127.8000

ans =

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

---