function?[?ss?ii?rr?]?=?CPSIR（?DegreeBetaSourse_IDTppvv??）
%??刻畫SIR模型
%???返回[ssiirr]
%???ss表示每一時刻傳播節點的密度
%???ii表示每一時刻未感染節點的密度
%???rr表示每一時刻免疫節點的密度
%???輸入變量的含義：
%??????Degree---------度分布矩陣
%??????Beta------------傳播概率矩陣
%??????links------------全邊矩陣
%??????Sourse_ID---------源傳播節點

N?=?length（Degree）;
A?=?zeros（N3）;%產生一個N行3列的零矩陣存儲節點狀態第一列為傳播狀態，第二列為未感染狀態，第三列為免疫狀態
A（Sourse_ID1）?=?1;%假設第一個節點為傳播節點
for?e?=?1:N
????if?e?~=Sourse_ID
????????A（e2）=1;
????end
end
%?p1?=?0.09;
p2=pp;%免疫概率
v=vv;%傳播自然轉移為免疫的速度
ss?=?zeros（1T）;%產生一個1行T列的零矩陣，存儲傳播節點密度
ii?=?zeros（1T）;
rr?=?zeros（1T）;

for?t?=?1:T
????disp（num2str（t））;
????S?=?0;%傳播節點數量
????I?=?0;
????R?=?0;
????for?i?=?1:N
????????S?=?S?+?A（i1）;
????????I?=?I?+?A（i2）;
????????R?=?R?+?A（i3）;
????end
????ss（1t）?=?S/N;
????ii（1t）?=?I/N;
????rr（1t）?=?R/N;
????
????%%
????%對節點遍歷，根據傳播規則進行交互
????for?i?=?1:N
????????%如果節點i為傳播節點，且它的度大于0
????????if?Degree（i）>0?&&?A（i1）?==?1
????????????%確定節點i在Beta矩陣中的最開始那行的id（i_uplimit）
????????????i_uplimit?=?1;
????????????if?i?==?1???????????????%將序號為1的特殊處理
????????????????i_uplimit?=?1;
????????????else
????????????????for?j?=?2:i
????????????????????i_uplimit?=?i_uplimit?+?Degree（j-1）;
????????????????end
????????????end
????????????
????????????%contact?process?CP接觸過程，每一時刻只與一個鄰居交互
????????????rand_nei?=?randi（Degree（i）?1）?+?i_uplimit?-1;%找出該鄰居的編號
????????????%判斷該鄰居是未感染節點，則進行傳播
????????????if?A（Beta（rand_nei2）2）==1
????????????????a?=?rand（）;
????????????????
????????????????if?a????????????????????A（Beta（rand_nei2）1）?=?1;
????????????????????A（Beta（rand_nei2）2）?=?0;
????????????????end
????????????end
????????????
????????????if?A（Beta（rand_nei2）3）==1?||?A（Beta（rand_nei2）1）==1%Beta（k1）表示的是Beta矩陣中第k行第一列的節點id
????????????????b?=?rand（）;
????????????????if?b????????????????????A（i3）?=?1;
????????????????????A（i1）?=?0;
????????????????end
????????????end????????????
????????end
????end
????%%
????%傳播節點可有較小概率自然衰減成免疫節點
????for?i?=?1:N
????????if?A（i1）?==1
????????????c?=rand（）;
????????????if?c????????????????A（i3）?=?1;
????????????????A（i1）?=?0;
????????????end
????????end
????end
????
end
%?plot（ss‘+-‘）;
%?hold?on;
%?plot（ii‘*‘）;
%?plot（rr‘.‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2603??2012-09-29?08:55??CPSIR.m