 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 2341 | 回復(fù): 1 | ||
ldiny新蟲(chóng) (小有名氣)
|
[求助]
SI模型病毒擴(kuò)散程序問(wèn)題求助 已有1人參與
|
|
小弟最近在做SI模型下的病毒擴(kuò)散，因?yàn)槭褂玫氖菆D的廣度優(yōu)先原則來(lái)進(jìn)行擴(kuò)散，即一次只能是一個(gè)節(jié)點(diǎn)在進(jìn)行擴(kuò)散。實(shí)際上，每個(gè)被感染的節(jié)點(diǎn)都會(huì)同時(shí)向周圍擴(kuò)散，請(qǐng)問(wèn)如下的代碼邏輯該怎么改？為什么出來(lái)的結(jié)果不是我想要的結(jié)果。 %記錄以當(dāng)前節(jié)點(diǎn)作為傳播源的開(kāi)始時(shí)間 t1=cputime; %按照?qǐng)D的廣度優(yōu)先的原則來(lái)進(jìn)行病毒傳播 function [practiceTime,infection]=BFSspreading(A,N,source_infection_node,irate,rrate,susceptible,infection,recover,time_gap,t1,obesever) %初始化隊(duì)列 %隊(duì)列頭 head=1; %隊(duì)列尾，開(kāi)始隊(duì)列為空，tail==head tail=1; %向頭中加入感染源節(jié)點(diǎn) queue(head)=source_infection_node; %隊(duì)列擴(kuò)展 .這樣保證頭在右邊，隊(duì)尾在左邊了 head=head+1; %開(kāi)始按照廣度優(yōu)先搜索順序向鄰居節(jié)點(diǎn)傳播 %判斷隊(duì)列是否為空 while tail~=head %記錄當(dāng)前已經(jīng)經(jīng)過(guò)了多少傳播結(jié)束時(shí)間 t2=cputime; if (t2-t1)>time_gap tail=head; end %取隊(duì)尾節(jié)點(diǎn) i=queue(tail); %如果該節(jié)點(diǎn)不在移除列表之中 if isempty(find(recover==i,1)) for j=1:N %如果節(jié)點(diǎn)j與當(dāng)前節(jié)點(diǎn)i相連并且節(jié)點(diǎn)j不在感染列表和恢復(fù)列表中，即節(jié)點(diǎn)j與當(dāng)前節(jié)點(diǎn)i相連，并且還沒(méi)有被感染過(guò) if A(i,j)==1 && isempty(find(infection==j,1)) && isempty(find(recover==j,1)) %以一定的概率向周圍傳播 infection_random=rand(1); if infection_random < irate %新節(jié)點(diǎn)入列 queue(head)=j; %擴(kuò)展隊(duì)列 head=head+1; %將新節(jié)點(diǎn)j加入感染列表 infection=[infection j]; %從易感染節(jié)點(diǎn)列表中刪除該節(jié)點(diǎn),設(shè)置為-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 |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
非良人
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5
- 帖子: 4
- 在線: 51分鐘
- 蟲(chóng)號(hào): 9929960
- 注冊(cè): 2018-08-18
2樓2018-08-19 21:30:44
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ldiny 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|