无标度BA模型matlab源码

clear all;
% 生成标准BA网络 A为邻接矩阵sparse matrix
% m0 = 4;% m0 初始节点，全连接
% m = 3;% m 加一个点的同时加的边数
% N= 1000;% N 总的节点数
m0=4;m=3;N=1000;
A = sparse（NN）;
A（1:m0 1:m0） = 1;
A（1:m0 1:m0） = A（1:m0 1:m0） - eye（m0）;
for i = （m0 + 1）:N
    
     i;

     Degree=sum（A（1:i-1 1:i-1））;
     DegreeSum=sum（Degree）;
     DegreeProability=（Degree）/DegreeSum;
     p = rand（1m）;
     for j = 1:m
         index = 1;
         k = 1;
         s = 0;
         while index == 1
                s = s + DegreeProability（k）;
                if p（j） <= s 
                    if A（i k） == 0
                        index = 0;
                        A（i k） = 1;
                        A（k i） = 1;
                    else