原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
第一部分程序：获取不同d直时的网络演化情况并保存成数据文件
clearall;clc;
N_delta_array=4;
Array_delta=[0.0005,0.001,0.0015,0.0025];
fornda=1:N_delta_array
%initialnumberofpoints%totalnumberofstepsneededrunning%numberofnewlinksforeachnewpoints
%initialnumberofpoints
%totalnumberofstepsneededrunning
%numberofnewlinksforeachnewpoints
%probabilityofinitiallinks
%probabilityofisolatepointtobeselected
%probabilityofsaturatedpointtobeselected
%factorofincreasingresource
n0=20;
T=3000;
m2=10;
p_ini=0.5;
p_iso=0.01;
p_sat=0.0002;
delta=Array_delta(nda);
%%%%%%%Getinitialtotalnumber
N=n0; %initialtotalnumberofpoints
ti=zeros(N,1);lij=zeros(N,N);ki=zeros(N,1);
ti=zeros(N,1);
lij=zeros(N,N);
ki=zeros(N,1);
mi=zeros(N,1);
pij=zeros(N,N);
wij=zeros(N,N);
%type
%link
%degree
%resource
%performance
%weight%%%%%%%%%%%%%%%%%
isolate=zeros(N,1);%isolatestate
saturation=zeros(N,1);%saturatedstate
%efficiencyofcooperation%varianceofcooperateefficiency%averagelength%aggregationfactor%averagedegreeofnet%varianceofKt%max(ki)-min(ki)%averageofthedifference%numberofsizesateachstep
%efficiencyofcooperation%varianceofcooperateefficiency%averagelength
%aggregationfactor%averagedegreeofnet%varianceofKt%max(ki)-min(ki)%averageofthedifference%numberofsizesateachstep
Dgt=zeros(T,1);
Lt=zeros(T,1);
Zt=zeros(T,1);
Kt=zeros(T,1);
Vt=zeros(T,1);
lt=zeros(T,1);
ut=zeros(T,1);
St=zeros(T,1);
%%%%%%%Initializationbeforthefirststep
%%%%%%%%%%%%%%%%%%Initializetypesofeachpointsfornti=1:N
ti(nti)=fix(4*rand(1,1))+1;
ifti(nti)==4
ti(nti)=2;
end
end
%%couldorntbeconnected
could_cot=zeros(N,N);
forni=1:N
fornj=1:N
if(ni?=nj)((ti(ni)?=ti(nj))||((ti(ni)==2)(ti(nj)==2)))could_cot(ni,nj)=1;
end
end
end
could_cot_ini=could_cot;
%%%%%%%%%%%%%%%%%%Initializelinks
por_i=rand(1,N);
tempones=ones(1,N);
key_lij=meshgrid(por_i
文档评论（0）