function [DeD,aver_DeD]=Degree_Distribution(A)
%% 求網(wǎng)絡圖中各節(jié)點的度及度的分布曲線
%% 求解算法：求解每個節(jié)點的度，再按發(fā)生頻率即為概率，求P(k) 
%A————————網(wǎng)絡圖的鄰接矩陣
%DeD————————網(wǎng)絡圖各節(jié)點的度分布
%aver_DeD———————網(wǎng)絡圖的平均度
N=size(A,2);
DeD=zeros(1,N);
for i=1:N
   % DeD(i)=length(find((A(i,:)==1)));
   DeD(i)=sum(A(i,:));
end
aver_DeD=mean(DeD);

if sum(DeD)==0
    disp('該網(wǎng)絡圖只是由一些孤立點組成');
    return;
else 
    figure;     
    bar([1:N],DeD);  
    xlabel('節(jié)點編號n');
    ylabel('各節(jié)點的度數(shù)K');
    title('網(wǎng)絡圖中各節(jié)點的度的大小分布圖');
end

figure;
M=max(DeD);
for i=1:M+1;    %網(wǎng)絡圖中節(jié)點的度數(shù)最大為M,但要同時考慮到度為0的節(jié)點的存在性
    N_DeD(i)=length(find(DeD==i-1));
end
P_DeD=zeros(1,M+1);
P_DeD(:)=N_DeD(:)./sum(N_DeD);
bar([0:M],P_DeD,'r');
xlabel('節(jié)點的度 K');
ylabel('節(jié)點度為K的概率 P(K)');
title('網(wǎng)絡圖中節(jié)點度的概率分布圖');