[转载]Matlab中根据邻接矩阵做无向图/有向图

Matlab中根据邻接矩阵做图

 

function tu_plot(rel,control)

%由邻接矩阵画图

%输入为邻接矩阵，必须为方阵；

%第二个输入为控制量，0表示无向图，1表示有向图。默认值为0

 

r_size=size(rel);

if nargin<2

    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*(theta<pi)-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 nargin<3

    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

的有向图：

正在上传…重新上传取消