早知道有向圖和無向圖差別沒有想象中的大我就寫到一起了。
函數中使用的arrow畫箭頭函數是在這個網站下的。
1 %函數名netplot 2 %使用方法輸入請help netplot 3 %無返回值 4 %函數只能處理有向圖 5 %作者：tiandsp 6 %最后修改：2012.12.26
7 function netplot(A,flag) 8 %調用方法輸入netplot(A,flag)，無返回值 9 %A為鄰接矩陣或關聯矩陣 10 %flag=1時處理鄰接矩陣 11 %flag=2時處理關聯矩陣 12 %函數只能處理有向圖 13 if flag==1 %鄰接矩陣表示有向圖 14 D_netplot(A); 15 return; 16 end 17
18 if flag==2 %關聯矩陣表示有向圖 19 [m n]=size(A); %關聯矩陣變鄰接矩陣 20 W=zeros(m,m); 21 for i=1:n 22 a=find(A(:,i)~=0); 23 if A(a(1),i)==1
24 W(a(1),a(2))=1; 25 else
26 W(a(2),a(1))=1; 27 end 28 end 29 D_netplot(W); 30 return; 31 end 32
33 function D_netplot(A) 34 [n n]=size(A); 35 w=floor(sqrt(n)); 36 h=floor(n/w); 37 x=[]; 38 y=[]; 39 for i=1:h %使產生的隨機點有其范圍，使顯示分布的更廣 40 for j=1:w 41 x=[x 10*rand(1)+(j-1)*10]; 42 y=[y 10*rand(1)+(i-1)*10]; 43 end 44 end 45 ed=n-h*w; 46 for i=1:ed 47 x=[x 10*rand(1)+(i-1)*10]; 48 y=[y 10*rand(1)+h*10]; 49 end 50 plot(x,y,'r*'); 51
52 title('網絡拓撲圖'); 53 for i=1:n 54 for j=1:n 55 if A(i,j)~=0
56 c=num2str(A(i,j)); %將A中的權值轉化為字符型 57 text((x(i)+x(j))/2,(y(i)+y(j))/2,c,'Fontsize',10); %顯示邊的權值 58 arrow([x(i) y(i)],[x(j) y(j)]); %帶箭頭的連線 59 end 60 text(x(i),y(i),num2str(i),'Fontsize',14,'color','r'); %顯示點的序號 61 hold on; 62 end 63 end 64 end 65
66 end
運行結果：
