版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xwndf250

銀蟲 (小有名氣)

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

[求助] matlab處理常微分方程作圖問題

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

回復(fù)此樓

» 猜你喜歡

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

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

Matlab常微分方程組求解問題已經(jīng)有4人回復(fù)
一階常微分方程組，matlab進(jìn)行數(shù)值求解已經(jīng)有3人回復(fù)
matlab作圖怎樣使坐標(biāo)刻度線明顯已經(jīng)有6人回復(fù)
matlab解一階常微分方程組已經(jīng)有4人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
matlab中二維圖作圖已經(jīng)有8人回復(fù)
matlab 線性常微分方程求解已經(jīng)有5人回復(fù)
關(guān)于Matlab作圖的問題已經(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)的問題求助已經(jīng)有3人回復(fù)
【求助】matlab解常微分方程，暫停和繼續(xù) 已經(jīng)有5人回復(fù)
【求助】求Matlab循環(huán)作圖已經(jīng)有12人回復(fù)
科研作圖時(shí)候matlab中使用latex 已經(jīng)有3人回復(fù)
【求助】先解常微分方程然后畫圖【已完結(jié)】已經(jīng)有11人回復(fù)
【求助】請教Matlab高手積分作圖已經(jīng)有3人回復(fù)

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

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

csgt0

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

彩色掛圖

專家經(jīng)驗(yàn): +2
計(jì)算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號(hào): 449046
注冊: 2007-11-02
專業(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沒有看懂，我記得求導(dǎo)貌似用diff函數(shù)的，能不能詳細(xì)講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

);”的 ...

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

贊一下

回復(fù)此樓

showmethemoney

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

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

查看全部 8 個(gè)回答

csgt0

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

彩色掛圖

專家經(jīng)驗(yàn): +2
計(jì)算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號(hào): 449046
注冊: 2007-11-02
專業(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

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時(shí)
蟲號(hào): 1943212
注冊: 2012-08-17
專業(yè): 農(nóng)業(yè)昆蟲學(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沒有看懂，我記得求導(dǎo)貌似用diff函數(shù)的，能不能詳細(xì)講解下odeset（）函數(shù)，intvalue作用，“for i=1:length(T) 和vdy(:,i)=rigid(T(i),Y(i,

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

贊一下

回復(fù)此樓

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

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

csgt0

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

彩色掛圖

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xwndf250: 金幣+15, ★★★很有幫助, 非常感謝您的幫助，給您金幣，順便問下，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的回帖

查看全部 8 個(gè)回答

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

[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao

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

24小時(shí)熱門版塊排行榜

xwndf250

[求助] matlab處理常微分方程作圖問題

» 猜你喜歡

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

csgt0

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

csgt0

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

xwndf250

csgt0

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

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

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