 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 基于WS小世界下的SIR模型
| 1 | 1/1 | 返回列表 |
| 查看: 2764 | 回復: 0 | ||
xwndf250銀蟲 (小有名氣)
|
[求助]
基于WS小世界下的SIR模型
|
|
流行病模型 sir dI/dt=a*S*I-b*I dS/dt=-a*S*I dR/dt=b*I m文件 function y=SIR(t,x) a=10;b=0.4; y=[a*x(1)*x(2)-b*x(1); -a*x(1)*x(2); b*x(3)]; 窗口命令 ts=0:50; x0=[0.04,0.96,0]; [t,x]=ode45('SIR',ts,x0); figure,plot(t,x(:,1),t,x(:,2),t,x( :,3)) 最近聽說基于WS小世界網(wǎng)絡系統(tǒng)可以解決SIR模型問題，而且更加真實，我想請問下怎樣基于WS小世界網(wǎng)絡系統(tǒng)構建SIR模型，求matlab代碼。 這是找到的一個其他小世界模型，可能對大神們有幫助 function ws_net() disp('WS小世界網(wǎng)絡模型') N=input('請輸入網(wǎng)絡節(jié)點數(shù)'); K=input('請輸入與節(jié)點左右相鄰的K/2的節(jié)點數(shù)'); p=input('請輸入隨機重連的概率'); angle=0:2*pi/N:2*pi-2*pi/N; x=100*cos(angle); y=100*sin(angle); plot(x,y,'r.','Markersize',30); hold on; %生成最近鄰耦合網(wǎng)絡； A=zeros(N); disp(A); for i=1:N if i+K<=N for j=i+1:i+K A(i,j)=1; end else for j=i+1:N A(i,j)=1; end for j=1  (i+K)-N)A(i,j)=1; end end if K A(i,j)=1; end else for j=1:i-1 A(i,j)=1; end for j=N-K+i:N A(i,j)=1; end end end disp(A); %隨機化重連 for i=1:N for j=i+1:N if A(i,j)==1 pp=unifrnd(0,1); if pp<=p A(i,j)=0; A(j,i)=0; b=unidrnd(N); while i==b b=unidrnd(N); end A(i,b)=1; A(b,i)=1; end end end end %根據(jù)鄰接矩陣連線 for i=1:N for j=1:N if A(i,j)==1 plot([x(i),x(j)],[y(i),y(j)],'linewidth',1); hold on; end end end hold off aver_path=aver_pathlength(A); disp(aver_path); |
回復此樓» 猜你喜歡
329求調(diào)劑
已經(jīng)有3人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有11人回復
-
一志愿蘇州大學材料求調(diào)劑，總分315（英一）
已經(jīng)有3人回復
-
0703化學調(diào)劑
已經(jīng)有14人回復
-
320求調(diào)劑0856
已經(jīng)有3人回復
-
307求調(diào)劑
已經(jīng)有9人回復
-
288求調(diào)劑
已經(jīng)有15人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有5人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有14人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有18人回復
相關版塊跳轉
我要訂閱樓主 xwndf250 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學0703化學318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 考研化學學碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|