matlab 根据邻接矩阵作图.docVIP

Matlab中根据邻接矩阵做图 function tu_plot(rel,control) %由邻接矩阵画图 %输入为邻接矩阵，必须为方阵； %第二个输入为控制量，0表示无向图，1表示有向图。默认值为0 r_size=size(rel); if nargin2 control=0; end if r_size(1)~=r_size(2) disp(Wrong Input! The input must be a square matrix!); return; end len=r_size(1); rho=10;%限制图尺寸的大小 r=2/1.05^len;%点的半径 theta=0:(2*pi/len):2*pi*(1-1/len); [pointx,pointy]=pol2cart(theta,rho); theta=0:pi/36:2*pi; [tempx,tempy]=pol2cart(theta,r); point=[pointx,pointy]; hold on for i=1:len temp=[tempx,tempy]+[point(i,1)*ones(length(tempx),1),point(i,2)*ones(length(tempx),1)]; plot(temp(:,1),temp(:,2),r); text(point(i,1)-0.3,point(i,2),num2str(i)); %画点 end for i=1:len for j=1:len if rel(i,j) link_plot(point(i,:),point(j,:),r,control); %连接有关系的点 end end end set(gca,XLim,[-rho-r,rho+r],YLim,[-rho-r,rho+r]); axis off %% function link_plot(point1,point2,r,control) %连接两点 temp=point2-point1; if (~temp(1))(~temp(2)) return; %不画子回路； end theta=cart2pol(temp(1),temp(2)); [point1_x,point1_y]=pol2cart(theta,r); point_1=[point1_x,point1_y]+point1; [point2_x,point2_y]=pol2cart(theta+(2*(thetapi)-1)*pi,r); point_2=[point2_x,point2_y]+point2; if control arrow(point_1,point_2); else plot([point_1(1),point_2(1)],[point_1(2),point_2(2)]); end %% function arrow(start,stop,l) %start,stop分别为起点和终点 %l为箭头的线长度，默认为主线长的1/10 t=0.1; ang=15/180*pi; temp=stop(1)-start(1)+j*(stop(2)-start(2)); L=abs(temp);P=angle(temp); if nargin3 l=t*L; end p1=P-ang;p2=P+ang; a=[stop(1)-l*cos(p1) stop(2)-l*sin(p1)]; b=[stop(1)-l*cos(p2) stop(2)-l*sin(p2)]; hold on plot([start(1) stop(1)],[start(2) stop(2)]); plot([a(1) stop(1)],[a(2) stop(2)]); plot([b(1) stop(1)],[b(2) stop(2)]); end 效果图如下： 邻接矩阵为 0 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 0 0 1 1 0 1 0 1 0 0 0 1 1 0 0 0 1 0 的有向图： 小程序，小兴趣。欢迎加入编程爱好者群一起成长。