早知道有向图和无向图差别没有想象中的大我就写到一起了。
函数中使用的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
运行结果: