24小时热门版块排行榜

返回列表

ldiny

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 3243.8
散金: 10
帖子: 75
在线: 56.7小时
虫号: 1547413
注册: 2011-12-22
专业: 基础物理学

[求助] SI模型病毒扩散程序问题求助已有1人参与

小弟最近在做SI模型下的病毒扩散，因为使用的是图的广度优先原则来进行扩散，即一次只能是一个节点在进行扩散。实际上，每个被感染的节点都会同时向周围扩散，请问如下的代码逻辑该怎么改？为什么出来的结果不是我想要的结果。

  %记录以当前节点作为传播源的开始时间
   t1=cputime;

%按照图的广度优先的原则来进行病毒传播
function [practiceTime,infection]=BFSspreading(A,N,source_infection_node,irate,rrate,susceptible,infection,recover,time_gap,t1,obesever)
%初始化队列
%队列头
head=1;
%队列尾，开始队列为空，tail==head
tail=1;
%向头中加入感染源节点
queue(head)=source_infection_node;
%队列扩展 .这样保证头在右边，队尾在左边了
head=head+1;

%开始按照广度优先搜索顺序向邻居节点传播
%判断队列是否为空
while tail~=head
%记录当前已经经过了多少传播结束时间
t2=cputime;
if (t2-t1)>time_gap
      tail=head;
end
%取队尾节点
i=queue(tail);
%如果该节点不在移除列表之中
if isempty(find(recover==i,1))
   for j=1:N
         %如果节点j与当前节点i相连并且节点j不在感染列表和恢复列表中，即节点j与当前节点i相连，并且还没有被感染过
         if A(i,j)==1 && isempty(find(infection==j,1)) && isempty(find(recover==j,1))

               %以一定的概率向周围传播
               infection_random=rand(1);
               if infection_random < irate
                  %新节点入列
                  queue(head)=j;
                  %扩展队列
                  head=head+1;
                  %将新节点j加入感染列表
                  infection=[infection j];

                  %从易感染节点列表中删除该节点,设置为-1
                  [row,col,v] = find(susceptible==j) ;
                  susceptible(col)=-1;
                  susceptible(find(susceptible==-1))=[];
               end
         end
   end

      tail=tail+1;
end %end if  isempty(find(recover==i,1)
end %end while
end

回复此楼