| 2 | 1/1 | 返回列表 |
| 查看: 2342 | 回復(fù): 1 | ||
ldiny新蟲 (小有名氣)
|
[求助]
SI模型病毒擴(kuò)散程序問題求助 已有1人參與
|
|
小弟最近在做SI模型下的病毒擴(kuò)散,因?yàn)槭褂玫氖菆D的廣度優(yōu)先原則來進(jìn)行擴(kuò)散,即一次只能是一個(gè)節(jié)點(diǎn)在進(jìn)行擴(kuò)散。實(shí)際上,每個(gè)被感染的節(jié)點(diǎn)都會(huì)同時(shí)向周圍擴(kuò)散,請(qǐng)問如下的代碼邏輯該怎么改?為什么出來的結(jié)果不是我想要的結(jié)果。 %記錄以當(dāng)前節(jié)點(diǎn)作為傳播源的開始時(shí)間 t1=cputime; %按照?qǐng)D的廣度優(yōu)先的原則來進(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ì)列尾,開始隊(duì)列為空,tail==head tail=1; %向頭中加入感染源節(jié)點(diǎn) queue(head)=source_infection_node; %隊(duì)列擴(kuò)展 .這樣保證頭在右邊,隊(duì)尾在左邊了 head=head+1; %開始按照廣度優(yōu)先搜索順序向鄰居節(jié)點(diǎn)傳播 %判斷隊(duì)列是否為空 while tail~=head %記錄當(dāng)前已經(jīng)經(jīng)過了多少傳播結(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相連,并且還沒有被感染過 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 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 梁梁校校 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|