复杂网络简单指标计算_matlab_图_复杂网络_

function [K,meank]=deg(A,x) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% 函数功能：已知图G的邻接矩阵A（不带权重），求节点的度数以及网络的平均度和度分布。 %%% 作者：Liyifan %%% x=1 有向图；x=0 无向图 %%% K为节点对应度数矩阵，meank为网络的平均度 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% N=size(A,1); %%%矩阵行数等于节点数 % K=diag(A*A); %%%计算节点的度（邻接矩阵二次幂的对角） % meank=sum(K)/N; %%%计算网络的平均度 K=zeros(N,1); out_K=zeros(N,1); in_K=zeros(N,1); if x==0 %%%%%无向图求度 for i=1:N K(i)=sum(sum(A(i,:)~=0)); end meank=sum(K)/N; %%%计算网络的平均度 end if x==1 %%%%%有向图求出度 for i=1:N out_K(i)=sum(sum(A(i,:)~=0)); %%%出度 in_K(i)=sum(sum(A(:,i)~=0)); %%%入度 K(i)=out_K(i)+in_K(i); end meank=sum(K)/N; %%%计算网络的平均度 end %%%检验度数矩阵中是否有0元素，如果有说明存在孤立点。 if sum(sum(K==0))>0 disp('网络存在孤立点') end %统计度数为1到kmax的节点分别有多少个 kmax=max(K); p=zeros(1,kmax);%%%用以存放度数为col时，节点的概率（频率代替概率） for col=1:kmax po=find(K==col); p(col)=length(po)/N; end plot(p,'*'); ylabel('P(k)'),xlabel('k') title('网络度分布');