matlab图论程序算法大全推荐文档.docxVIP

精心整理 图论算法matlab实现 求最小费用最大流算法的 MATLAB程序代码如下： n=5;C=[0 15 16 0 0 0 0 0 13 14 011 017 0 0 0 0 0 8 0 0 0 0 0]; %弧容量 b=[0 4 1 0 0 0 0 0 6 1 》 0 2 0 3 0 0 0 0 0 2 0 0 0 0 0]； %弧上单位流量的费用 wf=0;wf0=lnf; %wf表示最大流量，wf0 表示预定的流量值 for (i=1:n) for (j=1:n)f(i,j)=0; end; end %取初始可行流 f 为零流 while (1) for (i=1:n) for (j=1:n) if (j~=i)a(i,j)=lnf; end; end; end%勾造有向赋权图 for (i=1:n) for (j=1:n) if (C(i,j)0f(i,j)==0)a(i,j)=b(i,j); elseif (C(i,j)0f(i,j)==C(i,j))a(j,i)二-b(i,j); elseif (C(i,j)0)a(i,j)=b(i,j);a(j,i)=-b(i,j); end; end; end for (i=2:n)p(i)=lnf;s(i)二i; end %用 Ford 算法求最短路，赋初值 for (k=1:n)pd=1; %求有向赋权图中vs到vt的最短路 2019年9月 精心整理 endfor (i=2:n) for (j=1:n) if (p(i)p(j)+a(j,i))p(i)=p(j)+a(j,i);s(i)二j;pd=O; end ;end; end if (pd) break; end;end %求最短路的Ford算法结束 if (p(n)==lnf) break; end %不存在vs到vt的最短路，算法终止.注意在求最小费 用最大流时构造有 向赋权图中不会含负权回路，所以不会出现k=n dvt=Inf;t=n; %进入调整过程，dvt 表示调整量 while (1) %计算调整量 if (a(s(t),t)0)dvtt=C(s(t),t)-f(s(t),t); %前向弧调整量 elseif (a(s(t),t)0)dvtt=f(t,s(t)); end %后向弧调整量 if (dvtdvtt)dvt=dvtt; end ■i ? I x 1 i if (s(t)==1) break; end %当t的标号为vs时，终止计算调整量 j-J— — _____ t=s(t); end %继续调整前一段弧上的流f pd=0; if (wf+dvt=wf0)dvt=wf0-wf;pd=1; end%如果最大流量大于或等于预定的流量 值 ,--j I f t=n; while (1) %调整过程 if (a(s(t),t)O)f(s(t),t)=f(s(t),t)+dvt; %前向弧调整 elseif (a(s(t),t)0)f(t,s(t))=f(t,s(t))-dvt; end %后向弧调整 if (s(t)==1) break; end %当t的标号为vs时，终止调整过程 t=s(t); end if (pd) break; end%如果最大流量达到预定的流量值 wf=0; for (j=1: n)wf=wf+f(1,j); end; end %计算最大流量 end; end%计算最小费用zwf=0; for (i=1:n) for (j=1:n)zwf二zwf+b(i,j)*f(i,j); 2019年9 end; end%计算最小费用 精心整理 f %显示最小费用最大流 图 6-22 wf %显示最小费用最大流量 zwf %显示最小费用，程序结束_ Kruskal避圈法： Kruskal避圈法的MATLAB程序代码如下： n=8;A=[0 2 8 1 0 0 0 0 2 0 6 0 1 0 0 0 8 6 0 7 5 1 2 0 1 0 7 0 0 0 9 0 0 1 5 0 0 3 0 8 ■i ? I x 1 i 0 0 1 0 3 0 4 6 0 0 2 9 0 4 0 3 0 0 0 0 8 6 3 0]; k=1; %己录A中不同正数的个数 ,--j I f for (i=1:n-1) for (j=i+1:n) %此循环是查找A中所有不同的正数 if (A(i,j)0)x(k)=A(i,j); %数组x 记录A中不同的正数 kk=1; %缶时变量 for (s=1:k-1) if (x(k)==x(s))kk=0; break; end; end %排除相同的正数 k=k+kk; end; end; end k=k-1 %显示A中所有不同正数的个数 for (i=1: