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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

xwndf250

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 474.7
帖子: 117
在線: 46.1小時
蟲號: 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 想要做一個橫坐標(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的變化率隨時間的變化），注意，不是t與S、I、R的圖像！最好用matlab語言。

回復(fù)此樓

» 猜你喜歡

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

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

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微分方程作圖-報錯報錯已經(jīng)有4人回復(fù)
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復(fù)
Matlab中如何調(diào)用M文件作圖已經(jīng)有9人回復(fù)
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
【求助】matlab解常微分方程，暫停和繼續(xù) 已經(jīng)有5人回復(fù)
【求助】求Matlab循環(huán)作圖已經(jīng)有12人回復(fù)
科研作圖時候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的回帖

xwndf250

銀蟲 (小有名氣)

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

引用回帖:

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

好的謝謝你

回復(fù)此樓

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

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

查看全部 8 個回答

csgt0

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

彩色掛圖

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

【答案】應(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時的初值

[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小時
蟲號: 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時的初值
= 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)驗: +2
計算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

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

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

查看全部 8 個回答

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

[考研] 086000調(diào)劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 295求調(diào)劑 +4	1428151015 2026-03-27	5/250	2026-03-27 19:53 by JourneyLucky
[考研] 283求調(diào)劑（080500） +4	A child 2026-03-27	4/200	2026-03-27 15:34 by XPU李慶
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4	0857求調(diào)劑. 2026-03-24	4/200	2026-03-27 14:03 by 靜靜靜靜靜靜靜?/a>
[考研] 321求調(diào)劑 +4	璞玉~~ 2026-03-25	4/200	2026-03-27 11:41 by sanrepian
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-27	3/150	2026-03-27 10:17 by guoweigw
[考研] 求調(diào)劑 +5	蘆lty 2026-03-25	6/300	2026-03-26 20:49 by 不吃魚的貓
[考研] 279求調(diào)劑 +6	紅衣隱官 2026-03-21	6/300	2026-03-26 18:32 by 不吃魚的貓
[考研] 289求調(diào)劑 +17	碩星赴 2026-03-23	17/850	2026-03-26 16:18 by 不吃魚的貓
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 284求調(diào)劑 +15	Zhao anqi 2026-03-22	15/750	2026-03-25 12:51 by wht0531
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00

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

24小時熱門版塊排行榜

xwndf250

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

» 猜你喜歡

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

xwndf250

csgt0

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

xwndf250

csgt0

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

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

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