有了之前的准备工作,就可以进行接下来的BA无标度网络的生成了:
function Ans = BA(m0,M0,t,m)
%BA(m0,M0,t,m)
%先生成一个m0个节点、共M0条连边的网络,保证每个点都有至少一条连边
%随后生长t个节点,每个节点需要m条连边
%m<m0; 定义m0为偶数
% M0<=1/2*m0*(m0-1)
%----------------生成一个初始网络----------------------
A = regular(m0,M0,t);%规则网络
%--------------开始生长t个节点,每个节点向网络连接m条边------
for i = m0+1:m0+t
x = i;
B = P_degree(A,i-1);
for j = 1:m
p = rand();
y = Choosevertex(p,B);
while A(x,y) == 1 || x-y == 0
p = rand();
y = Choosevertex(p,B);
end
A(x,y) = 1;
A(y,x) = 1;
end
end
Ans = A;
%----------------画出图像------------------------
% figure
% AA = A(1:m0+t-1,1:m0+t-1);
% C = graph(A,'upper');
% CC = graph(AA,'upper');
% D = degree(CC);
% e = max(D);
% f = find(D>=(e-1))';
% F = find(D==