function  [C,aver_C]=Clustering_Coefficient(A)
%% 求網絡圖中各節點的聚類系數及整個網絡的聚類系數
%% 求解算法：求解每個節點的聚類系數，找某節點的所有鄰居，這些鄰居節點構成一個子圖
%% 從A中抽出該子圖的鄰接矩陣，計算子圖的邊數，再根據聚類系數的定義，即可算出該節點的聚類系數
%A————————網絡圖的鄰接矩陣
%C————————網絡圖各節點的聚類系數
%aver———————整個網絡圖的聚類系數
N=size(A,2);
C=zeros(1,N);
for i=1:N
    aa=find(A(i,:)==1);  %尋找子圖的鄰居節點
    if isempty(aa)
        disp(['節點',int2str(i),'為孤立節點，其聚類系數賦值為0']);
        C(i)=0;
    else        
        m=length(aa); 
        if m==1
            disp(['節點',int2str(i),'只有一個鄰居節點，其聚類系數賦值為0']);
            C(i)=0;
        else
           B=A(aa,aa);          % 抽取子圖的鄰接矩陣
           C(i)=length(find(B==1))/(m*(m-1));
        end
    end
end
aver_C=mean(C);