MATLAB 矩阵生成无向网络图

　　寝室有一个是做网络的，需要把矩阵变为网络图使其形象一点，所以我在这里就写了这样一个函数。

　　功能是将邻接矩阵或关联矩阵变为网络图，不过我这里只能转换为无向图，有向图的箭头我还需要在研究一下，似乎有annotation函数可以调用。

%函数名netplot
%使用方法输入请help netplot
%无返回值%函数只能处理无向图
%作者：tiandsp
%最后修改：2012.12.26
function netplot(A,flag)
    %调用方法输入netplot(A,flag)，无返回值
    %A为邻接矩阵或关联矩阵
    %flag=1时处理邻接矩阵
    %flag=2时处理关联矩阵
    %函数只能处理无向图
    if flag==1      %邻接矩阵表示无向图
        ND_netplot(A);
        return;
    end
    
    if flag==2      %关联矩阵表示无向图
        [m n]=size(A);      %关联矩阵变邻接矩阵
        W=zeros(m,m);
        for i=1:n
            a=find(A(:,i)~=0);
            W(a(1),a(2))=1;
            W(a(2),a(1))=1;
        end
        ND_netplot(W);
        return;
    end
           
    function ND_netplot(A)
        [n n]=size(A);
        w=floor(sqrt(n));       
        h=floor(n/w);        
        x=[];
        y=[];
        for i=1:h           %使产生的随机点有其范围，使显示分布的更广
            for j=1:w
                x=[x 10*rand(1)+(j-1)*10];
                y=[y 10*rand(1)+(i-1)*10];
            end
        end
        ed=n-h*w;
        for i=1:ed
           x=[x 10*rand(1)+(i-1)*10]; 
           y=[y 10*rand(1)+h*10];
        end
        plot(x,y,'r*');    
        
        title('网络拓扑图'); 
        for i=1:n
            for j=i:n
                if A(i,j)~=0
                    c=num2str(A(i,j));                      %将A中的权值转化为字符型              
                    text((x(i)+x(j))/2,(y(i)+y(j))/2,c,'Fontsize',10);  %显示边的权值
                    line([x(i) x(j)],[y(i) y(j)]);      %连线
                end
                text(x(i),y(i),num2str(i),'Fontsize',14,'color','r');   %显示点的序号            
                hold on;
            end
        end  
    end
    
end

调用方法里面都写了，我就不多说了。

下面是一个节点为10的图（搞图像的我果然还是喜欢传图啊）：