phoenix708

應(yīng)助: 0 (幼兒園)
金幣: 9829.8
紅花: 1
帖子: 1443
在線: 122.8小時(shí)
蟲號(hào): 476825
注冊(cè): 2007-12-13
性別: GG
專業(yè): 催化化學(xué)

[求助] 會(huì)用matlab的同志幫幫忙哈！

dX/dt=1.98S(1-X)-0.73XY
dY/dt=dX/dt-0.00207Y/(1+Y+2S)-0.059d[Y/(1+Y+2S)]/dt
2dS/dt=-dX/dt-0.00056S/(1+Y+2S)-0.086d[S/(1+Y+2S)]/dt
當(dāng)t=0時(shí),X=0,Y=0,S=1,
其中 t的最大值為3，
求X，Y, S與t的關(guān)系圖
謝謝!

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有159人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收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)價(jià)值貼推薦，對(duì)您同樣有幫助:

需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
幫忙看看這個(gè)平面圖用matlab怎么畫已經(jīng)有19人回復(fù)
速求用matlab寫 armchair型石墨烯納米條帶的能帶圖已經(jīng)有8人回復(fù)
如何用matlab軟件或者用C及C++自己手工編程尋找轉(zhuǎn)折點(diǎn) 已經(jīng)有12人回復(fù)
求助�。∏笠怀绦�，用matlab程序做，用蒙特卡洛方法模擬已經(jīng)有10人回復(fù)
跪求XRD標(biāo)準(zhǔn)卡片，PDF編號(hào)：01-072-7074，希望各位同志幫幫忙吧！！已經(jīng)有15人回復(fù)
【求助】matlab 計(jì)算一個(gè)二重積分，高手們幫幫忙吧已經(jīng)有6人回復(fù)
求助：誰會(huì)用matlab畫圖？已經(jīng)有13人回復(fù)
【求助】會(huì)用Matlab的LMI toolbox的請(qǐng)幫幫忙~~ 已經(jīng)有5人回復(fù)
【求助】有發(fā)過論文的化工前輩過來幫幫忙哈已經(jīng)有3人回復(fù)
【求助】求matlab軟件，要求能用已經(jīng)有5人回復(fù)
【求助】正態(tài)分布如何用MATLAB檢驗(yàn) 已經(jīng)有6人回復(fù)

1樓 2012-01-05 11:47:14

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

scumat

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 50.5
散金: 126
紅花: 53
帖子: 920
在線: 165.7小時(shí)
蟲號(hào): 1191957
注冊(cè): 2011-01-18
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

問題很難，20個(gè)金幣少了。

贊一下

回復(fù)此樓

2樓2012-01-05 16:40:17

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

phoenix708

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 9829.8
紅花: 1
帖子: 1443
在線: 122.8小時(shí)
蟲號(hào): 476825
注冊(cè): 2007-12-13
性別: GG
專業(yè): 催化化學(xué)

引用回帖:

樓: Originally posted by scumat at 2012-01-05 16:40:17:
問題很難，20個(gè)金幣少了。

你能解嗎？能解就可以加金幣

贊一下

回復(fù)此樓

3樓2012-01-05 16:46:35

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

phoenix708

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 9829.8
紅花: 1
帖子: 1443
在線: 122.8小時(shí)
蟲號(hào): 476825
注冊(cè): 2007-12-13
性別: GG
專業(yè): 催化化學(xué)

我把方程組化簡了一下，不知對(duì)否？懂得蟲友幫忙運(yùn)算下，應(yīng)該不難！x'=1.98S*(1-X)-0.73*X*Y;
y'=x'-0.00207Y/(1+Y+2*S)-0.059*(y'/(1+y+2*s)+y*(y'+2*s')/(1+y+2*s)^2);
s'=0.5*(-x'-0.00056S/(1+Y+2*S)-0.086*(s'/(1+y+2*s)+s*(y'+2*s')/(1+y+2*s)^2))
其中y'=dy/dt, s'=ds/dt，x'=dx/dt

贊一下

回復(fù)此樓

4樓2012-01-05 17:49:11

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

dbb627

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

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

感謝參與，應(yīng)助指數(shù) +1

>> syms t
>> Y=sym('Y(t)')

Y =

Y(t)

>> S=sym('S(t)')

S =

S(t)

>> F=Y/(1+Y+2*S)

F =

Y(t)/(2*S(t) + Y(t) + 1)

>> dF=diff(F,t)

dF =

diff(Y(t), t)/(2*S(t) + Y(t) + 1) - (Y(t)*(2*diff(S(t), t) + diff(Y(t), t)))/(2*S(t) + Y(t) + 1)^2

1：先將公式中的項(xiàng)“d[Y/(1+Y+2S)]/dt”求導(dǎo)展開變?yōu)椋?br /> y'/(1+y+2*s)+y*(y'+2*s')/(1+y+2*s)^2，其中y'=dy/dt, s'=ds/dt
2：同理將“d[S/(1+Y+2S)]/dt”變?yōu)椋?br /> s'/(1+y+2*s)+s*(y'+2*s')/(1+y+2*s)^2
用matlab ode15i解隱式微分方程

CODE:

odefun=@(t,x,dx)[dx(1)-1.98*x(3)*(1-x(1))+0.73*x(1)*x(2)

    dx(2)-(dx(1)-0.00207*x(2)/(1+x(2)+2*x(3))-0.059*(dx(2)/(1+x(2)+2*x(3))+x(2)*(dx(2)+2*dx(3))/(1+x(2)+2*x(3))^2))

    2*dx(3)+dx(1)+0.00056*x(3)/(1+x(2)+2*x(3))+0.086*(dx(3)/(1+x(2)+2*x(3))+x(3)*(dx(2)+2*dx(3))/(1+x(2)+2*x(3))^2)]; 

 t0=0;

 x0=[0 0 1]';

fix_x0=ones(3,1);

 dx0=[0 1 1]';

fix_dx0=zeros(3,1);

 [x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0); 

solution=ode15i(odefun,[0 3],x02,dx02);

 t=solution.x;

x=solution.y;

plot(t,x)

legend('X','Y','S')

[ Last edited by dbb627 on 2012-1-5 at 18:29 ]

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5樓2012-01-05 18:23:40

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

dbb627

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

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

phoenix708(金幣+30): ★★★很有幫助 2012-01-05 20:02:35

做圖如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

phoenix708

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2012-01-05 18:25:08

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

芝芝花

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 77.6
帖子: 48
在線: 80.7小時(shí)
蟲號(hào): 1042403
注冊(cè): 2010-06-16
性別: MM
專業(yè): 分離過程

引用回帖:

6樓: Originally posted by dbb627 at 2012-01-05 18:25:08:
做圖如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

Thank you!

回復(fù)此樓

7樓2012-01-05 19:23:56

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

scumat

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 50.5
散金: 126
紅花: 53
帖子: 920
在線: 165.7小時(shí)
蟲號(hào): 1191957
注冊(cè): 2011-01-18
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

6樓: Originally posted by dbb627 at 2012-01-05 18:25:08:
做圖如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

佩服！可以得全分。

回復(fù)此樓

8樓2012-01-05 19:40:59

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

phoenix708

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 9829.8
紅花: 1
帖子: 1443
在線: 122.8小時(shí)
蟲號(hào): 476825
注冊(cè): 2007-12-13
性別: GG
專業(yè): 催化化學(xué)

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2012-01-05 18:25:08:
做圖如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

很厲害啊，不過這個(gè)結(jié)果和我想要的結(jié)果有差距，我糾正了第二個(gè)和第三個(gè)方程的系數(shù)，肯定是對(duì)的，勞駕大俠最后幫我計(jì)算一次！多謝多謝！
dX/dt=1.98S(1-X)-0.73XY
dY/dt=dX/dt-4.6Y/(1+Y+2S)-0.059d[Y/(1+Y+2S)]/dt
2dS/dt=-dX/dt-1.25S/(1+Y+2S)-0.086d[S/(1+Y+2S)]/dt
當(dāng)t=0時(shí),X=0,Y=0,S=1,
其中 t的最大值為3，
求X，Y, S與t的關(guān)系圖
謝謝!

贊一下

回復(fù)此樓

9樓2012-01-05 19:52:47

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

kingslcq

新蟲 (正式寫手)

愛科學(xué)的PH.D & Prof.

應(yīng)助: 56 (初中生)
貴賓: 0.028
金幣: 785.2
散金: 450
紅花: 28
帖子: 894
在線: 40.7小時(shí)
蟲號(hào): 585944
注冊(cè): 2008-08-08
專業(yè): 結(jié)構(gòu)工程

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

感謝參與，應(yīng)助指數(shù) +1
phoenix708(金幣+1): 2012-01-05 20:00:10

http://www.gaoyang168.com/bbs/viewthread.php?tid=4017234

贊一下

回復(fù)此樓

10樓2012-01-05 19:54:13

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

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

[考研] 274求調(diào)劑 +15	顧九笙要謙虛 2026-03-24	21/1050	2026-03-27 12:20 by 3Strings
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-27	3/150	2026-03-27 10:17 by guoweigw
[考研] 359求調(diào)劑 +4	王了個(gè)楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 316求調(diào)劑 +5	江辭666 2026-03-26	5/250	2026-03-27 08:08 by hypershenger
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 321求調(diào)劑 +5	材料cailiao 2026-03-21	5/250	2026-03-26 20:41 by fmesaito
[考研] 求調(diào)劑 +8	Auroracx 2026-03-22	8/400	2026-03-26 19:55 by 不吃魚的貓
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6	小小陳小小 2026-03-22	9/450	2026-03-26 15:58 by dick_runner
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 296求調(diào)劑 +5	www_q 2026-03-20	5/250	2026-03-26 12:56 by 3Strings
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO