根據(jù)鄰接矩陣作圖

1、Matlab中根據(jù)鄰接矩陣做圖function tu_plot(rel,control)%由鄰接矩陣畫圖%輸入為鄰接矩陣，必須為方陣；%第二個(gè)輸入為控制量，0表示無向圖，1表示有向圖。默認(rèn)值為0r_size=size(rel);if nargin2 control=0;endif r_size(1)=r_size(2) disp(Wrong Input! The input must be a square matrix!); return;endlen=r_size(1);rho=10;%限制圖尺寸的大小r=2/1.05len;%點(diǎn)的半徑theta=0:(2*pi/len):2*pi*(1-

2、1/len);pointx,pointy=pol2cart(theta,rho);theta=0:pi/36:2*pi;tempx,tempy=pol2cart(theta,r);point=pointx,pointy;hold onfor 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); %畫點(diǎn)endfor i=1:l

3、en for j=1:len if rel(i,j) link_plot(point(i,:),point(j,:),r,control); %連接有關(guān)系的點(diǎn) end endendset(gca,XLim,-rho-r,rho+r,YLim,-rho-r,rho+r);axis off%function link_plot(point1,point2,r,control)%連接兩點(diǎn)temp=point2-point1;if (temp(1)&(temp(2) return; %不畫子回路；endtheta=cart2pol(temp(1),temp(2);point1_x,point1_y=p

4、ol2cart(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分別為起點(diǎn)和終點(diǎn)%l為箭頭的線長(zhǎng)度，默認(rèn)為主線長(zhǎng)的1/10t=

5、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;endp1=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 onplot(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