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)有247人回復
哈爾濱理工大學物理系招收物理學考研調(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)有24人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有20人回復

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

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人回復
【求助】先解常微分方程然后畫圖【已完結】已經(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ā)表

[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	6/300	2026-03-17 17:13 by 生物工程調(diào)劑
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +3	share_joy 2026-03-16	3/150	2026-03-17 13:41 by houyaoxu
[考研] 材料與化工304求B區(qū)調(diào)劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學習賺彩禮
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant