(圖論)matlab模板程序-免費(fèi)閱讀
【正文】 for (j=1:n) wf=wf+f(1,j)。pd=1。t=n。s(i)=i。 for (i=1:n) a(i,i)=0。0 0 0 6 1。for (j=1:n) wf=wf+f(1,j)。 end dvt=d(n)。 if (d(j)d(i)) d(j)=d(i)。d=zeros(1,n)。243。0 0 0 0 0 0 0 3。inq=find(q==1)。 for i=1:n(1) if all(e(:,2)inp(i))==0 inq=find((e(:,2)inp(i))==0)。p(e(inp,1))=0。 t=tsum(d(ep,:))sum(d(:,eq))+1。 numq=sum(d(:,inp))。t=sum(sum(d)39。 d=(b==0)。)。 W(t1,:)=0。第五講：匹配問題及算法程序一：較大基礎(chǔ)匹配算法function J=matgraf(W)n=size(W,1)。 dot=[str1,str2]。endsubplot(1,2,2)。)。n=size(d,2)。 text(v(i,1)1,v(i,2)2,dot)。subplot(1,2,1)hold on。25 38。 endend程序二：Hamilton改良圈算法（找出比較好的Hamilton路）function [C d1]= hamiltonglf(v)%d表示權(quán)值矩陣%C表示算法最終找到的Hamilton圈。ddd=1。 if sum(m1)==0 dvex1(dd)=edd(kk)。dvex=0。 if f==1 flag=0。 ed(2,eds)~=0 T=ed。 vexs(1,1)=vet。 endendif m~=0 fprintf(39。b(b~=0)=1。 e=e+1。 mina(i,j) min=a(i,j)。程序二：最小生成樹的Prim算法function [T c]=Primf(a)l=length(a)。 c=c+B(3,i)。 [B,i]=sortrows(b39。 k=1。 P=[p1 p2(2:length(p2)) p3(2:length(p3))]。[p2 d2]=n2short(W,t1,t2)。 k=k+1。kk=k2。m=1。end程序二：floyd算法（計(jì)算任意兩點(diǎn)間的最短距離）function [d,r]=floyd(a) n=size(a,1)。程序一：Dijkstra算法（計(jì)算兩點(diǎn)間的最短路）function [l,z]=Dijkstra(W)n = size (W,1)。)。 end end endelseif f==1 %關(guān)聯(lián)矩陣轉(zhuǎn)換為鄰接矩陣 m=size(F,2)。endW。 n=size(F,1)。程序二：無向圖關(guān)聯(lián)矩陣和鄰接矩陣互換算法F表示所給出的圖的相應(yīng)矩陣W表示程序運(yùn)行結(jié)束后的結(jié)果f=0表示把鄰接矩陣轉(zhuǎn)換為關(guān)聯(lián)矩陣f=1表示把關(guān)聯(lián)矩陣轉(zhuǎn)換為鄰接矩陣%無向圖的關(guān)聯(lián)矩陣和鄰接矩陣的相互轉(zhuǎn)換function W=incandadf(F,f)if f==0 %鄰接矩陣轉(zhuǎn)換為關(guān)聯(lián)矩陣 m=sum(sum(F))/2。for i=2:n P=P+A^i。 %給邊的始點(diǎn)賦值為1 W(j,k)=1。 endelse fprint(39。 for i=1:n for j=i:n if F(i,j)~=0 %由i發(fā)出的邊，有向邊的始點(diǎn) W(i,k)=1。 %有向邊由a(1)指向a(2) else W(a(2),a(1))=1。while m=nfor i=1:nfor j=1:n if D(i,j)D(i,m)+D(m,j) D(i,j)+ D(i,m)+D(m,j)。 z(i)=j1。 r(i,j)=r(i,k)。k=1。 end endendk=1。 Pm{i}=P。[p6 d6]=n2short(W,t1,k2)。第三講：最小生成樹程序一：最小生成樹的Kruskal算法function [T c]=krusf(d,flag)if nargin==1 n=size(d,2)。 end end endelse b=d。 k=1。 end end end if k==n break。e=1。 distance(e)=b。第四講：Euler圖和Hamilton圖程序一：Fleury算法（在一個(gè)Euler圖中找出Euler環(huán)游）注：包括三個(gè)文件；, , function [T c]=fleuf1(d)%注：必須保證是Euler環(huán)游，否則輸出T=0,c=0n=length(d)。vexs=zeros(1,eds+1)。 flag=0。 while flagg [flagg ed]=edf(matr,eds,vexs,ed,tem)。 break。 else break。dd1=0。 edd1=[tem。 elseif templength(dvex1) amp。18 54。18 40。
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1