xwndf250

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時
蟲號: 1943212
注冊: 2012-08-17
專業(yè): 農(nóng)業(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]，縱坐標dI/dt、dS/dt、dR/dt的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)];
end
命令:應該怎樣寫，才是t與dI/dt、dS/dt、dR/dt的圖像（就是S、I、R的變化率隨時間的變化），注意，不是t與S、I、R的圖像！最好用matlab語言。

回復此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有156人回復
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有5人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有26人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有31人回復

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

Matlab常微分方程組求解問題已經(jīng)有4人回復
一階常微分方程組，matlab進行數(shù)值求解已經(jīng)有3人回復
matlab作圖怎樣使坐標刻度線明顯已經(jīng)有6人回復
matlab解一階常微分方程組已經(jīng)有4人回復
高手進來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復
matlab中二維圖作圖已經(jīng)有8人回復
matlab 線性常微分方程求解已經(jīng)有5人回復
關于Matlab作圖的問題已經(jīng)有14人回復
【【吐血求助】】Matlab微分方程作圖-報錯報錯已經(jīng)有4人回復
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復
Matlab中如何調(diào)用M文件作圖已經(jīng)有9人回復
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復
【求助】matlab解常微分方程，暫停和繼續(xù) 已經(jīng)有5人回復
【求助】求Matlab循環(huán)作圖已經(jīng)有12人回復
科研作圖時候matlab中使用latex 已經(jīng)有3人回復
【求助】先解常微分方程然后畫圖【已完結(jié)】已經(jīng)有11人回復
【求助】請教Matlab高手積分作圖已經(jīng)有3人回復

1樓 2013-02-25 13:02:03

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★
感謝參與，應助指數(shù) +1
dbb627: 金幣+2, 感謝應助 2013-02-26 14:07:58

CODE:

function xwn

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);

intvalue=[0.1 0.4 0.5];       %t=0時的初值

[T,Y] = ode45(@rigid,[0 50],intvalue,options);

plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')

title('T-Y圖')

for i=1:length(T)     

vdy(:,i)=rigid(T(i),Y(i,:));

end

figure

plot(T,vdy(1,:),'-',T,vdy(2,:),'-.',T,vdy(3,:),'.')

title('T-dY圖')

end

function dy = rigid(t,y)

dy = zeros(3,1);    % a column vector

a=0.2;

b=0.1;

dy(1) = a*y(2)*y(1)-b * y(1);

dy(2) = -a*y(2) * y(1);

dy(3) = b * y(1);

end

贊一下

回復此樓

showmethemoney

2樓2013-02-26 11:06:11

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

xwndf250

銀蟲 (小有名氣)

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

引用回帖:

2樓: Originally posted by csgt0 at 2013-02-26 11:06:11
function xwn
options = odeset('RelTol',1e-4,'AbsTol',);
intvalue=; %t=0時的初值
= ode45(@rigid,,intvalue,options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
title('T-Y圖')
for i=1: ...

function dy = rigid(t,y)后面的我看懂了，但是前面function xwn沒有看懂，我記得求導貌似用diff函數(shù)的，能不能詳細講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

);”的意思。我是小白，求指教。好的話我多給金幣。

贊一下

回復此樓

3樓2013-02-26 21:06:43

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xwndf250: 金幣+15, ★★★很有幫助, 非常感謝您的幫助，給您金幣，順便問下，t0初值怎么確定的？也就是0.1，0.4，0.5是怎么確定的？ 2013-02-27 10:00:25
fegg7502: 金幣+1, 應助指數(shù)+1, 鼓勵交流 2013-04-02 09:30:14

odeset（）函數(shù)設置數(shù)值解的計算精度，其實一般不用也可以
intvalue是初值啊，t=0時的初值
“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,);”的為了循環(huán)計算每個t下的3個導數(shù)值

贊一下

回復此樓

showmethemoney

4樓2013-02-27 09:54:46

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★
fegg7502: 金幣+1, 應助指數(shù)+1, 鼓勵交流 2013-04-02 09:30:24

引用回帖:

3樓: Originally posted by xwndf250 at 2013-02-26 21:06:43
function dy = rigid(t,y)后面的我看懂了，但是前面function xwn沒有看懂，我記得求導貌似用diff函數(shù)的，能不能詳細講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

);”的 ...

就是t=0是的ISR啊，具體多少得看你的實際情況，我只是隨意寫的個數(shù)。

贊一下

回復此樓

showmethemoney

5樓2013-02-27 10:37:52

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

xwndf250

銀蟲 (小有名氣)

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

引用回帖:

5樓: Originally posted by csgt0 at 2013-02-27 10:37:52
就是t=0是的ISR啊，具體多少得看你的實際情況，我只是隨意寫的個數(shù)。...

好的謝謝你

回復此樓

6樓2013-02-27 11:09:25

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

657801288

鐵蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 319.8
帖子: 39
在線: 16.6小時
蟲號: 2098895
注冊: 2012-10-28
專業(yè): 控制理論與方法

引用回帖:

你好，兩種方法做出的圖怎么不一樣呢？求指教。

贊一下

回復此樓

7樓2013-03-30 08:49:39

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準晶和納米晶
管轄: 計算模擬

【答案】應助回帖

★
fegg7502: 金幣+1, 應助指數(shù)+1, 鼓勵交流 2013-04-02 09:30:38

不是兩種方法啊，前面畫的是T-Y圖，后面是你要的T-dY圖。

贊一下

回復此樓

showmethemoney

8樓2013-04-01 10:48:30

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

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

[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[考研] 0703化學336分求調(diào)劑 +4	zbzihdhd 2026-03-15	5/250	2026-03-17 17:33 by ruiyingmiao
[考研] 311求調(diào)劑 +8	冬十三 2026-03-15	8/400	2026-03-17 16:59 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[文學芳草園] 伙伴們，祝我生日快樂吧 +17	myrtle 2026-03-10	26/1300	2026-03-16 18:32 by 青橙Ln
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[基金申請] 有必要更換申報口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky