1、%網(wǎng)絡版function Spreading_in_Networks_SIS2()clear all;clc;M=300;N=6;beta1=0; delta1=0.4%N個節(jié)點在M時刻的感染比例,p_state=zeros(N,M);q_state=zeros(2*N,M);format long;%迭代的步長step=0.025;L1=BA_net(N);%先把laplacian矩陣處理為鄰接矩陣for i=1:N L1(i,i)=0;endL1=-L1;%初始函數(shù)p_state(:,1)=matix(N,0.2);q1=1;for i=1:N for j=1:N q1=q1*(1-L1(

2、i,j)*p_state(j,1)*beta1); end q_state(i,1)=q1;endp_11=zeros(1,length(step:step:1); k=1;for beta1=step:step:1 for t=2:M for i=1:N p_state(i,t)= (1-p_state(i,t-1)*(1-q_state(i,t-1)+(1-delta1)*p_state(i,t-1); end for i=1:N q1=1; for j=1:N q1=q1*(1-L1(i,j)*p_state(j,t)*beta1); end q_state(i,t)=q1; end e

3、nd p1=0; for i=1:N p1=p1+p_state(i,M); end p_11(1,k)=p1/N; k=k+1; %把數(shù)據(jù)清掉，只留下第一個初始矩陣的值，這個很重要。 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ǎng)絡function A=BA_net(num

4、) % 從已有的m0個節(jié)點的網(wǎng)絡開始，采用增長機制與優(yōu)先連接的機制生成BA無標度網(wǎng)絡 % A 返回生成網(wǎng)絡的鄰接矩陣 % m0=input(未增長前的網(wǎng)絡節(jié)點個數(shù)m0: ); % m=input( 每次引入的新節(jié)點時新生成的邊數(shù)m： ); % N=input(增長后的網(wǎng)絡規(guī)模N： ); % disp(初始網(wǎng)絡時m0個節(jié)點的連接情況：1表示都是孤立；2表示構成完全圖；3表示隨機連接一些邊); % pp=input(初始網(wǎng)絡情況1，2或3： ); m0=2; m=2; N=num; pp=2; if mm0 disp(輸入?yún)?shù)m不合法); return; end x=100*rand(1,m0);

5、 y=100*rand(1,m0); switch pp case 1 A=zeros(m0); case 2 A=-ones(m0); for i=1:m0 A(i,i)=0; end case 3 for i=1:m0 for j=i+1:m0 p1=rand(1,1); if p10.5 A(i,j)=1;A(j,i)=0; end end end otherwise disp(輸入?yún)?shù)pp不合法); return; end for k=m0+1:N M=size(A,1); p=zeros(1,M); x0=100*rand(1,1);y0=100*rand(1,1); x(k)=x0

6、;y(k)=y0; if isempty(find(A=-1)=0 p(:)=1/M; else for i=1:M p(i)=length(find(A(i,:)=-1)/length(find(A=-1); end end pp=cumsum(p); %求累計概率 for i=1:m %利用賭輪法從已有的節(jié)點中隨機選擇m個節(jié)點與新加入的節(jié)點相連 random_data=rand(1,1); aa=find(pp=random_data);jj=aa(1); % 節(jié)點jj即為用賭輪法選擇的節(jié)點 A(k,jj)=-1;A(jj,k)=-1; end end n=size(A); for i=1:n A(i,i)=-sum(A(i,:); endend%輸入：N行數(shù)，默認列數(shù)是：2. a是第3列的和.%輸出:輸出第2列前num/2行和為a,后num/2行和也為a,所有行和都為1的num*2矩陣.function ma