matlab1.zip_matlab_复杂网络 程序

A = load ('D:\matlab\busL空间 无0.txt'); A=sparse(A(:,1),A(:,2),ones(size(A,1),1)); %构造稀疏矩阵存储形式 A=full(A); % 稀疏矩阵转m×n矩阵 disp(A); for i=1:size(A,1) %转为对称矩阵 for j=1:size(A,2) if A(i,j)==1; A(j,i)=1; end end end N=size(A,1); DeD=zeros(1,N); count=zeros(1,N); h=zeros(1,N); for i=1:N DeD(i)=sum(A(i,:)); end H=repmat(DeD,size(A,1),1);%节点度矩阵 kj= A.*H; %邻居节点的度 for i=1:size(A,1) for m=1:size(A,1) for j=1:size(A,1) if kj(i,j)>=m count(i)=count(i)+1; end end if count(i)>=m h(i)=m; end count(i)=0; end end h = h'; xlswrite('D:\matlab\L空间H-index.xlsx',h) figure; M=max(h); for i=1:M+1; N_h(i)=length(find(h==i-1)); end P_h=zeros(1,M+1); P_h(:)=N_h(:)./sum(N_h); bar([0:M],P_h,'r'); xlabel('节点的h指数h'); ylabel('节点h指数为h的概率 P(h)');