无向图中最大团问题的matlab代码

无向图中最大团问题的matlab代码，利用的算法是回溯法，代码包含MCP函数、测试代码和根据邻接矩阵画无向图的函数；

function graph1（relcontrol）

%由邻接矩阵画图

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

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

r_size=size（rel）;%获取矩阵大小

if nargin<2 %如果参数小于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（2）;

!echo len;

disp（len）;

 

rho=10;%限制图尺寸的大小

r=1/1.05^len;%点的半径

theta=0:（2*pi/len）:2*pi*（1-1/len）;%以0开始，步长2*pi/len，结束2*pi*（1-1/len），相当于有len个点

[pointxpointy]=pol2cart（theta‘rho）;

theta=0:pi/36:2*pi;%73个点

!echo theta:

disp（theta）;

[tempxtempy]=pol2cart（theta‘r）;%73个坐标

point=[pointxpointy];%6个坐标点

!echo point;

disp（point）;

hold on

for i=1:len

    temp=[tempxtempy]+[point（i1）*ones（length（tempx）1）point（i2）*ones（length（tempx）1）];%73行1列的1矩阵与point每个坐标相乘

    plot（temp（:1）temp（:2）‘r‘）;%plot（xy）画出圆圈

     %plot（point（:1）point（:2）‘r‘）;%plot（xy）

    text（point（i1）-0.3point（i2）num2str（i））;

    %画点

end

for i=1:len

    for j=1:len

        if rel（ij）%如果邻接矩阵rel中该点是1，代表相连

            link_plot（point（i:）point（j:）rcontrol）;

            %连接有关系的点

        end

    end

end

set（gca‘XLim‘[-rho-rrho+r]‘YLim‘[-rho-rrho+r]）;

axis off

%%

function link_plot（point1point2rcontrol）

%连接两点

temp=point2-point1;

if （~temp（1））&&（~temp（2））

    return;

    %不画子回路；

end

theta=cart2pol（temp（1）temp（2））;

[point1_xpoint1_y]=pol2cart（thetar）;

point_1=[point1_xpoint1_y]+point1;

[point2_xpoint2_y]=pol2cart（theta+（2*（theta
point_2=[point2_xpoint2_y]+point2;

if control

    arrow（point_1point_2）;

else

    plot（[point_1（1）point_2（1）][point_1（2）point_2（2）]）;

end

%%

function arrow（startstopl）

%startstop分别为起点和终点

%l为箭头的线长度，默认为主线长的1/10

t=0.1;

ang=15/180*pi;

temp=stop（1）-start（1）+1i*（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）]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2422  2018-10-22 16:43  MCP_self\graph1.m

     文件        384  2018-10-22 18:39  MCP_self\main.m

     文件       1074  2018-10-22 16:33  MCP_self\MCP.m

     目录          0  2018-10-22 20:29  MCP_self

----------- ---------  ---------- -----  ----

                 3880                    4