BA模型拓扑复杂网络

% A=[0 1 1;1 0 0;1 0 0];
% find（A==1）
% length（find（A==1））
% p=zeros（1size（A1））;
% if length（find（A==1））==0
% size（A1）
% p（:）=1/size（A1）
% else
% for i=1:size（A1）
% p（i）=length（find（A（i:）==1））/length（find（A==1））;
% end
% end
% pw=cumsum（p）
m0=input（‘未增长前的网络节点个数m0:  ‘）;
m=input（‘ 每次引入的新节点时新生成的边数m： ‘）;
N=input（‘增长后的网络规模N： ‘）;
disp（‘初始网络时m0个节点的连接情况：1表示都是孤立；2表示构成完全图；3表示随机连接一些边‘）;
pp=input（‘初始网络情况1，2或3： ‘）;
if m>m0
    disp（‘输入参数m不合法‘）;
    return;
end
x=100*rand（1m0）;
y=100*rand（1m0）;

switch  pp
    case 1
        A=zeros（m0）;
    case 2
        A=ones（m0）;
        for i=1:m0
            A（ii）=0;%%除了自己和自己不连接外其余都是任意两两连接
        end
    case 3
        for i=1:m0
            for j=i+1:m0
                p1=rand（11）;
                if p1>0.5