xwndf250

應(yīng)助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時(shí)
蟲(chóng)號(hào): 1943212
注冊(cè): 2012-08-17
專(zhuān)業(yè): 農(nóng)業(yè)昆蟲(chóng)學(xué)

[求助] matlab表達(dá)常微分方程組

流行病模型 sir
dI/dt=a*S*I-b*I
dS/dt=-a*S*I
dR/dt=b*I
S+I+R=1且0 想要做一個(gè)橫坐標(biāo)t=[0,50]，縱坐標(biāo)S、I、R=[0,1]的3條曲線，求具體程序，本人寫(xiě)的程序
m文件：function y=SIR(t,x)
a=0.2;b=0.1;
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);
但并不完備怎樣將“S+I+R=1且0 ”加入程序？則樣做出3條曲線？

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有103人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專(zhuān)博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Matlab求解二階偏微分方程組，希望能給出相應(yīng)的求解代碼已經(jīng)有10人回復(fù)
一階常微分方程組，matlab進(jìn)行數(shù)值求解已經(jīng)有3人回復(fù)
matlab求解偏微分方程組已經(jīng)有3人回復(fù)
matlab解一階常微分方程組已經(jīng)有4人回復(fù)
高手進(jìn)來(lái)用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
求助四波混頻耦合波方程數(shù)值求解已經(jīng)有5人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
matlab解微分方程組已經(jīng)有15人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題已經(jīng)有21人回復(fù)
【求助】請(qǐng)教非齊次常微分方程組的解析解法已經(jīng)有4人回復(fù)
【求助】如何求解下列偏微分方程已經(jīng)有13人回復(fù)
【求助】急請(qǐng)微分方程高手看下這個(gè)方程組！已經(jīng)有14人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)

1樓 2013-02-20 21:50:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-02-21 10:18:00
xwndf250: 金幣+20, ★★★很有幫助, 非常感謝您的幫助，謝謝您，能得到您的聯(lián)系方式嗎？對(duì)于這個(gè)模型，有好幾個(gè)延伸版本，我不太懂，希望能得到您的幫助 2013-02-21 13:34:38

用1stOpt試了下：

CODE:

Constant a=0.2, b=0.1;

Variable t=[0,50], I=0.04, S=0.96, R=0;

Plot I, S, R;

ODEFunction I'=a*S*I-b*I;

            S'=-a*S*I;

            R'=b*I;

t i(t) s(t) r(t)
0 0.04 0.96 0
1 0.0438199698527025 0.95199137574788 0.00418865439941787
2 0.0479246796644068 0.943301831986009 0.00877348834958384
3 0.0523191638106323 0.933897582965654 0.0137832532237135
4 0.0570050541663188 0.923747906691054 0.0192470391426273
5 0.0619799165436078 0.912826181758511 0.0251939016978815
6 0.0672365833530442 0.901110994752353 0.0316524218946029
7 0.0727625125962009 0.888587287004862 0.0386502003989369
8 0.0785392095667703 0.875247499940034 0.0462132904931952
9 0.0845417527960957 0.86109266915365 0.0543655780502546
10 0.0907384689488048 0.84613340987996 0.0631281211712348
11 0.0970908016427613 0.830390731714896 0.0725184666423429
12 0.103553415706703 0.813896619568116 0.0825499647251813
13 0.110074570608999 0.796694321796341 0.0932311075946594
14 0.116596784505434 0.778838295959874 0.104564919534692
15 0.123057793921448 0.76039377772209 0.116548428356462
16 0.129391794490216 0.741435958463345 0.129172247046439
17 0.135530927009781 0.722048780849711 0.142420292140508
18 0.141406952437619 0.702323386861487 0.156269660700894
19 0.146953041618962 0.682356277170927 0.170690681210111
20 0.152105592755206 0.662247261658677 0.185647145586117
21 0.15680598361213 0.642097295926693 0.201096720461176
22 0.161002167227083 0.622006306184411 0.216991526588506
23 0.16465002942103 0.602071104047241 0.233278866531729
24 0.167714442786892 0.582383483822453 0.249902073390655
25 0.170169973198862 0.563028578980893 0.266801447820245
26 0.172001218904904 0.544083533701538 0.283915247393558
27 0.173202786382979 0.525616522035906 0.301180691581115
28 0.173778929011366 0.507686123830561 0.318534947158073
29 0.173742892376161 0.490341045239489 0.335916062384349
30 0.173116022554591 0.473620154087765 0.353263823357644
31 0.17192670055368 0.457552787484017 0.370520511962303
32 0.170209167506611 0.442159281224708 0.387631551268682
33 0.168002302009851 0.427451667418802 0.404546030571347
34 0.165348404220617 0.413434487700284 0.421217108079099
35 0.162292032248037 0.400105673460266 0.437602294291697
36 0.158878926139982 0.387457450715047 0.453663623144971
37 0.155155044405282 0.375477234579228 0.46936772101549
38 0.151165728301575 0.364148486021985 0.484685785676441
39 0.146955000595849 0.353451511022385 0.499593488381766
40 0.142564998465521 0.343364188973172 0.514070812561307
41 0.138035534753006 0.333862622956456 0.528101842290538
42 0.133403777866304 0.324921709222201 0.541674512911495
43 0.128704038081417 0.316515626845401 0.554780335073182
44 0.123967646641047 0.308618251200535 0.567414102158418
45 0.11922291362721 0.301203496697017 0.579573589675773
46 0.114495150885485 0.294245595311546 0.591259253802969
47 0.109806747087884 0.287719317979249 0.602473934932867
48 0.105177283161219 0.28160014600198 0.613222570836801
49 0.100623677632725 0.275864399419385 0.62351192294789
50 0.0961603528413486 0.270489328866861 0.633350318291791

約束條件不用加，算出來(lái)的結(jié)果都滿足。用Matlab結(jié)果也應(yīng)該是一樣的吧。

jv1.jpg

贊一下

回復(fù)此樓

2樓2013-02-21 09:50:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xwndf250

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時(shí)
蟲(chóng)號(hào): 1943212
注冊(cè): 2012-08-17
專(zhuān)業(yè): 農(nóng)業(yè)昆蟲(chóng)學(xué)

★ ★
fegg7502: 金幣+2, 鼓勵(lì)交流 2013-02-22 12:03:49

引用回帖:

2樓: Originally posted by dingd at 2013-02-21 09:50:27
用1stOpt試了下：

Constant a=0.2, b=0.1;
Variable t=, I=0.04, S=0.96, R=0;
Plot I, S, R;
ODEFunction I'=a*S*I-b*I;
S'=-a*S*I;
R'=b*I;

t i(t) s(t) r(t)
0 0.04 0.96 ...

有個(gè)SIR的延伸模型就是 SIRM模型
表達(dá)公式如下：
積極免疫
dS/dt=γβI^2+γ(-2α+β)SI+γβ(M_i+M_l )I+k_is M_l
dI/dt=-γβI^2-γ(-2α+β)SI-γβ(M_i+M_l )I-(γ-k_a)I
dR/dt=γI+k_ir M_i
(dM_i)/dt=k_a I-k_il M_i-k_ir M_i
(dM_l)/dt=k_il M_i-k_ls M_l
和
消極免疫
dS/dt=γβI^2+γ(-2α+β)SI+γβ(M_i+M_l )I-k_b SI+k_is M_l
dI/dt=-γβI^2-γ(-2α+β)SI-γβ(M_i+M_l )I-γI
dR/dt=γI+k_ir M_i
(dM_i)/dt=k_b SI-k_il M_i-k_ir M_i
(dM_l)/dt=k_il M_i-k_ls M_l
α=0.1、β=0.5、γ=0.1、r=0.1、ka=0.5、kb=0.5、kir=0.01、kil=0.01、kls=0.1。
S=5/6、I=1/6 其他初始值為0
求幫助，具體情況如圖中所示得出的結(jié)果為第一張圖片。。。

2012 Social Transfer of Pathogunisation in Ant Colonies PLOS 9.jpg

流行病模型 4.jpg

贊一下

回復(fù)此樓

3樓2013-02-21 13:47:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 295材料工程專(zhuān)碩求調(diào)劑 +3	1428151015 2026-03-27	3/150	2026-03-27 19:49 by JourneyLucky
[考研] 289求調(diào)劑 +7	新時(shí)代材料 2026-03-27	7/350	2026-03-27 18:08 by arrow8852
[考研] 286求調(diào)劑 +8	PolarBear11 2026-03-26	8/400	2026-03-27 18:05 by yu221
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 085602 化工專(zhuān)碩 338分求調(diào)劑 +8	路癡小琪 2026-03-27	8/400	2026-03-27 16:50 by lenovolsw
[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	12/600	2026-03-27 09:41 by caszguilin
[考研] 一志愿北化求調(diào)劑 +3	Jsman 2026-03-22	3/150	2026-03-26 21:06 by ajpv風(fēng)雷
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-26	6/300	2026-03-26 18:03 by 邱gl
[考研] 22 350 本科985求調(diào)劑，求老登收留 +4	李軼男003 2026-03-20	4/200	2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 調(diào)劑310 +3	溫柔的晚安 2026-03-25	4/200	2026-03-25 23:16 by peike
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 26考研-291分-廈門(mén)大學(xué)（085601）-柔性電子學(xué)院材料工程專(zhuān)業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] 0703化學(xué)求調(diào)劑 +6	奶油草莓. 2026-03-22	7/350	2026-03-25 10:00 by shangxh
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4	我愛(ài)學(xué)電池 2026-03-23	4/200	2026-03-25 00:59 by 1027_324
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00