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論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有284人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(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é)

我把方程組化簡(jiǎn)了一下，不知對(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ā)表

[考研] 287求調(diào)劑 +9	land xuxu 2026-03-26	9/450	2026-03-27 15:17 by caszguilin
[考研] 328求調(diào)劑 +3	嗯滴的基本都 2026-03-27	3/150	2026-03-27 15:05 by 無際的草原
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4	0857求調(diào)劑. 2026-03-24	4/200	2026-03-27 14:03 by 靜靜靜靜靜靜靜?/a>
[考研] 303求調(diào)劑 +7	安憶靈 2026-03-22	8/400	2026-03-27 11:46 by sanrepian
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7	丹青奶蓋 2026-03-26	8/400	2026-03-27 11:24 by 丹青奶蓋
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4	@taotao 2026-03-26	5/250	2026-03-27 08:10 by hypershenger
[考研] 333求調(diào)劑 +6	wfh030413@ 2026-03-23	6/300	2026-03-26 22:45 by 學(xué)員8dgXkO
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +8	y7czhao 2026-03-26	8/400	2026-03-26 19:51 by 不吃魚的貓
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 化學(xué)調(diào)劑 +6	yzysaa 2026-03-21	6/300	2026-03-25 09:27 by aa331100
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike