加权无标度网络边权演化bbv模型的matlab实现

该代码是生成加权无标度网络边权演化模型的matlab源代码。希望对大家的学习研究有帮助。

tic
%%初始化
n=8;
w0=2;
w=w0*ones（n）;   %%%构造初始网络
for i=1:n
    w（ii）=0;
end
while（n<5000）   
s=sum（w）;
Bj=find（s）; 
bsize=length（Bj）;
k=s（Bj）/sum（s（Bj））;     %%%选择连接：概率公式
k=cumsum（k）;
k=[0k];
%%选择连接点
w=[wzeros（n1）];
w=[w;zeros（1n+1）]; %%增加一个节点（孤立）.
r1=rand（1）;%连接第一个节点
for i=1:bsize
    if（k（i）<=r1 & k（i+1）>r1）
        i1=i;
        w（n+1Bj（i1））=w0;w（Bj（i1）n+1）=w0;  %增加连边
%%%%%原来与新节点相连的边权将发生变化
        s=sum（w）;  
        for j=1:n
                if w（Bj（i1）j）~=0
                    w（Bj（i1）j）=w（Bj（i1）j）+w（Bj（i1）j）*w0/s（Bj（i1））;  %%%  与之相关的边的权发生改变
                    w（jBj（i1））=w（jBj（i1））+w（jBj（i1））*w0/s（Bj（i1））; 
                end
        end
        break
    end
end
r2=rand（1）;%连接第二个节点
while（k（i1）<=r2 & k（i1+1）>r2）
    r2=rand（1）;
end
for i=1:bsize
    if（k（i）<=r2 & k（i+1）>r2）
        i2=i;
        w（n+1Bj（i2））=w0;
        w（Bj（i2）n+1）=w0;  %增加连边
%%%%%原来与新节点相连的边权将发生变化
        s=sum（w）;  
        for j=1:n
                if w（Bj（i2）j）~=0
                    w（Bj（i2）j）=w（Bj（i2）j）+w（Bj（i2）j）*w0/s（Bj（i2））;  %%%  与之相关的边的权发生改变
                    w（jBj（i2））=w（jBj（i2））+w（jBj（i2））*w0/s（Bj（i2））; 
                end
        end
        break
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%
n=n+1;
end
%%%1.strengh的分布图   ---取整再画图    
s2=fix（s）;
pp1=tabulate（s2）;
sw=find（pp1（:3）>0）;

subplot（221）;
loglog（pp1（sw1）pp1（sw3）/100‘k.‘）;
title（‘strengh的分布图‘）
xlabel（‘s‘）;
ylabel（‘P（s）‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%下面画出相应的度分布图
df=sum（w>0）;
pp2=tabulate（df）;
sw=find（pp2（:3）>0）;

subplot（222）;
loglog（pp2（sw1）pp2（sw3）/100‘k*‘）;
lsline     %最小二乘拟合直线
title（‘度分布图‘）
xlabel（‘k‘）;
ylabel（‘P（k）‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% 下面画出:度－势相关图

subplot（223）;
plot（dfs2‘k.‘）;  
lsline     %最小二乘拟合直线
title（‘度－势相关图 ‘）
xlabel（‘k‘）;
ylabel（‘s‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%  下面画出：权分布图
wf=fix（w）;
Aq=zeros（1n*n）;
for i=1:n
    Aq（1（（i-1）*n+1）:i*n）=wf（i:）;
end
wfw=find（wf>0）;
pp3=tabulate（Aq（wfw））;
sw=find（pp3（:3）>0）;

subplot（224）;
plot（pp3（sw1）pp3（sw3）/100‘ko‘）;
title（‘权分布图‘）
xlabel（‘w‘）;
ylabel（‘P（w）‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2429  2008-09-28 22:24  bbv.m

----------- ---------  ---------- -----  ----

                 2429                    1