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