無(wú)標(biāo)度網(wǎng)絡(luò)SIS模型傳播代碼.docx

《無(wú)標(biāo)度網(wǎng)絡(luò)SIS模型傳播代碼.docx》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。

1、%網(wǎng)絡(luò)版functionSpreading_in__Networks_SIS2()clearall;clc;M=300;N=6;beta1=0;delta1=0.4%N個(gè)節(jié)點(diǎn)在M時(shí)刻的感染比例,p_state=zeros(N,M);q_state=zeros(2*N,M);formatlong;%迭代的步長(zhǎng)step=0.025;L1=BA_net(N);%先把laplacian矩陣處理為鄰接矩陣fori=1:NL1(i,i)=0;endL1=-L1;%初始函數(shù)p_state(:,1)=matix(N,0.2);q1=1;fori=1:Nforj=1:Nq1=q1*(1

2、-L1(i,j)*p_state(j,1)*beta1);endq_state(i,1)=q1;endp_11=zeros(1,length(step:step:1));k=1;forbeta1=step:step:1fort=2:Mfori=1:Np_state(i,t)=(1-p_state(i,t-1))*(1-q_state(i,t-1))+(1-delta1)*p_state(i,t-1);endfori=1:Nq1=1;forj=1:Nq1=q1*(1-L1(i,j)*p_state(j,t)*beta1);endq_state(i,t)=q1;enden

3、dp1=0;fori=1:Np1=p1+p_state(i,M);endp_11(1,k)=p1/N;k=k+1;%把數(shù)據(jù)清掉，只留下第一個(gè)初始矩陣的值，這個(gè)很重要。p_state(:,2:M)=0;q_state(:,2:M)=0;endh=plot(step:step:1,p_11,'r-o');set(h,'LineWidth',2);set(gca,'FontSize',12);set(get(gca,'XLabel'),'FontSize',18);set(get(gca,'YLabel'),'FontSize',18);end%生成無(wú)標(biāo)度網(wǎng)絡(luò)functio

4、nA=BA_net(num)%%%從已有的m0個(gè)節(jié)點(diǎn)的網(wǎng)絡(luò)開(kāi)始，采用增長(zhǎng)機(jī)制與優(yōu)先連接的機(jī)制生成BA無(wú)標(biāo)度網(wǎng)絡(luò)%%A——————返回生成網(wǎng)絡(luò)的鄰接矩陣%m0=input('未增長(zhǎng)前的網(wǎng)絡(luò)節(jié)點(diǎn)個(gè)數(shù)m0:');%m=input('每次引入的新節(jié)點(diǎn)時(shí)新生成的邊數(shù)m：');%N=input('增長(zhǎng)后的網(wǎng)絡(luò)規(guī)模N：');%disp('初始網(wǎng)絡(luò)時(shí)m0個(gè)節(jié)點(diǎn)的連接情況：1表示都是孤立；2表示構(gòu)成完全圖；3表示隨機(jī)連接一些邊');%pp=input('初始網(wǎng)絡(luò)情況1，2或3：');m0=2;m=2;N=num;pp=2;ifm>m0disp('輸入?yún)?shù)m不合法');return;en

5、dx=100*rand(1,m0);y=100*rand(1,m0);switchppcase1A=zeros(m0);case2A=-ones(m0);fori=1:m0A(i,i)=0;endcase3fori=1:m0forj=i+1:m0p1=rand(1,1);ifp1>0.5A(i,j)=1;A(j,i)=0;endendendotherwisedisp('輸入?yún)?shù)pp不合法');return;endfork=m0+1:NM=size(A,1);p=zeros(1,M);x0=100*rand(1,1);y0=100*rand(1,1);x(k)=x0;y

6、(k)=y0;ifisempty(find(A==-1))==0p(:)=1/M;elsefori=1:Mp(i)=length(find(A(i,:)==-1))/length(find(A==-1));endendpp=cumsum(p);%求累計(jì)概率fori=1:m%利用賭輪法從已有的節(jié)點(diǎn)中隨機(jī)選擇m個(gè)節(jié)點(diǎn)與新加入的節(jié)點(diǎn)相連random_data=rand(1,1);aa=find(pp>=random_data);jj=aa(1);%節(jié)點(diǎn)jj即為用賭輪法選擇的節(jié)點(diǎn)A(k,jj)=-1;A(jj,k)=-1;endendn=size(A);fori=1:nA(i

7、,i)=-sum(A(i,:));endend%輸入：N行數(shù)，默認(rèn)列數(shù)是：2.a是第3列的和.%輸出:輸出第2列前num/2行和為a,后num/2行和也為a,所有行和都為1的num*2矩陣.functionmatix=matix(num,a)mydata1=rand(num,1);sum1=sum(mydata1);mydata1=mydata1/(sum1/a);matix=[mydata1;];end