xwndf250

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

[求助] matlab處理常微分方程作圖問(wèn)題

流行病模型 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)dI/dt、dS/dt、dR/dt的3條曲線(xiàn)，求具體程序，本人寫(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)];
end
命令:應(yīng)該怎樣寫(xiě)，才是t與dI/dt、dS/dt、dR/dt的圖像（就是S、I、R的變化率隨時(shí)間的變化），注意，不是t與S、I、R的圖像！最好用matlab語(yǔ)言。

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有259人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)

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

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

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

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

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

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+2, 感謝應(yīng)助 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時(shí)的初值

[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

贊一下

回復(fù)此樓

showmethemoney

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

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

xwndf250

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

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

引用回帖:

2樓: Originally posted by csgt0 at 2013-02-26 11:06:11
function xwn
options = odeset('RelTol',1e-4,'AbsTol',);
intvalue=; %t=0時(shí)的初值
= 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沒(méi)有看懂，我記得求導(dǎo)貌似用diff函數(shù)的，能不能詳細(xì)講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

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

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

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

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

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

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

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

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 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沒(méi)有看懂，我記得求導(dǎo)貌似用diff函數(shù)的，能不能詳細(xì)講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

);”的 ...

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

贊一下

回復(fù)此樓

showmethemoney

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

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

xwndf250

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

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

引用回帖:

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

好的謝謝你

回復(fù)此樓

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

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

657801288

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 319.8
帖子: 39
在線(xiàn): 16.6小時(shí)
蟲(chóng)號(hào): 2098895
注冊(cè): 2012-10-28
專(zhuān)業(yè): 控制理論與方法

引用回帖:

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

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

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

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

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

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

[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +3	WW.' 2026-03-23	4/200	2026-03-23 15:31 by cocolv
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 307求調(diào)劑 +11	冷笙123 2026-03-17	11/550	2026-03-22 20:16 by edmund7
[考研] 315分，誠(chéng)求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶(hù) 2026-03-17	9/450	2026-03-20 14:24 by 無(wú)懈可擊111
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438