版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6591)
>考研 (2333)
>導師招生 (2066)
>文獻求助 (359)
>蟲友互識 (301)
>碩博家園 (169)
>考博 (155)
>博后之家 (146)
>休閑灌水 (117)
>招聘信息布告欄 (112)
>論文投稿 (106)
>基金申請 (87)
>找工作 (46)
>公派出國 (41)
>綠色求助(高懸賞) (40)
>高分子 (39)

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

xwndf250

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時
蟲號: 1943212
注冊: 2012-08-17
專業(yè): 農業(yè)昆蟲學

[求助] matlab表達常微分方程組

流行病模型 sir
dI/dt=a*S*I-b*I
dS/dt=-a*S*I
dR/dt=b*I
S+I+R=1且0 想要做一個橫坐標t=[0,50]，縱坐標S、I、R=[0,1]的3條曲線，求具體程序，本人寫的程序
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條曲線？

回復此樓

» 猜你喜歡

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

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

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

dingd

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

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

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

約束條件不用加，算出來的結果都滿足。用Matlab結果也應該是一樣的吧。

jv1.jpg

贊一下

回復此樓

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

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

xwndf250

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時
蟲號: 1943212
注冊: 2012-08-17
專業(yè): 農業(yè)昆蟲學

★ ★
fegg7502: 金幣+2, 鼓勵交流 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 ...

有個SIR的延伸模型就是 SIRM模型
表達公式如下：
積極免疫
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
求幫助，具體情況如圖中所示得出的結果為第一張圖片。。。

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

流行病模型 4.jpg

贊一下

回復此樓

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

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

相關版塊跳轉我要訂閱樓主 xwndf250 的主題更新

返回列表

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

[考研] 0856，材料與化工321分求調劑 +6	大饞小子 2026-03-27	7/350	2026-03-27 21:52 by cury2062
[考研] 0703化學求調劑，各位老師看看我！�。� +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[有機交流] 高溫高壓反應求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 287求調劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 329求調劑 +4	星野? 2026-03-26	4/200	2026-03-27 12:00 by 不吃魚的貓
[考研] 0856材料化工調劑總分330 +8	zhubinhao 2026-03-27	8/400	2026-03-27 11:46 by Wushiqi17
[考研] 材料與化工考研調劑 +16	孅華 2026-03-22	16/800	2026-03-27 11:23 by 王笑宇一定是研?/a>
[考研] 考研調劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓
[考研] 0703化學一志愿南京師范大學303求調劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 一志愿吉大071010，316分求調劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 材料學碩，求調劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 081200-11408-276學碩求調劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 349求調劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 調劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 286求調劑 +13	Faune 2026-03-21	13/650	2026-03-26 19:52 by peike
[考研] 081700 調劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 289材料與化工（085600）B區(qū)求調劑 +4	這么名字咋樣 2026-03-22	5/250	2026-03-25 08:20 by mx.yue
[考研] 求調劑一志愿武漢理工大學材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 341求調劑(一志愿湖南大學070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

xwndf250

[求助] matlab表達常微分方程組

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

dingd

【答案】應助回帖

xwndf250

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频