版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

phoenix708

鐵桿木蟲 (著名寫手)

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

引用回帖:

樓: Originally posted by hls85915 at 2012-01-05 21:16:31:
你黏在MATLAB里就OK了

我沒有matlab的嘛，我沒用過，呵呵！麻煩幫忙運算下，把數據也發(fā)上來下，我要自己用origin畫圖，謝謝

贊一下

回復此樓

21樓2012-01-05 21:19:05

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

dbb627

榮譽版主 (著名寫手)

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
phoenix708(金幣+20): ★★★很有幫助答案很完善！ 2012-01-06 09:07:06
cenwanglai(金幣+10): 真是不知道如何感謝你~ 2012-01-11 20:47:04

CODE:

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

    dx(2)-(dx(1)-4.6*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)+1.25*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')

display('t  X   Y  S')

[t‘ x’]

t X Y S
0 0 0 1.0000
0.0000 0.0000 0.0000 1.0000
0.0000 0.0001 0.0001 1.0000
0.0001 0.0002 0.0002 0.9999
0.0002 0.0004 0.0004 0.9998
0.0004 0.0008 0.0008 0.9995
0.0006 0.0012 0.0012 0.9993
0.0008 0.0016 0.0016 0.9991
0.0012 0.0024 0.0023 0.9986
0.0016 0.0032 0.0031 0.9981
0.0024 0.0048 0.0047 0.9971
0.0041 0.0080 0.0078 0.9952
0.0073 0.0144 0.0140 0.9914
0.0139 0.0269 0.0261 0.9839
0.0204 0.0392 0.0378 0.9766
0.0270 0.0512 0.0492 0.9694
0.0335 0.0630 0.0602 0.9623
0.0466 0.0858 0.0812 0.9484
0.0597 0.1078 0.1009 0.9351
0.0858 0.1491 0.1367 0.9098
0.1120 0.1873 0.1681 0.8861
0.1382 0.2226 0.1954 0.8640
0.1643 0.2553 0.2190 0.8432
0.1905 0.2856 0.2394 0.8237
0.2429 0.3399 0.2715 0.7881
0.2952 0.3868 0.2940 0.7564
0.3475 0.4277 0.3085 0.7278
0.3999 0.4634 0.3167 0.7020
0.4522 0.4948 0.3198 0.6785
0.5045 0.5226 0.3187 0.6570
0.5569 0.5473 0.3144 0.6370
0.6615 0.5892 0.2988 0.6011
0.7662 0.6235 0.2773 0.5693
0.8709 0.6521 0.2530 0.5406
0.9755 0.6766 0.2278 0.5142
1.0802 0.6978 0.2031 0.4896
1.1849 0.7164 0.1797 0.4665
1.2895 0.7330 0.1581 0.4446
1.3942 0.7480 0.1384 0.4237
1.4989 0.7616 0.1207 0.4038
1.6035 0.7739 0.1051 0.3848
1.7082 0.7852 0.0914 0.3665
1.8129 0.7956 0.0795 0.3490
1.9175 0.8052 0.0691 0.3321
2.0222 0.8140 0.0601 0.3159
2.1269 0.8221 0.0524 0.3004
2.3362 0.8365 0.0400 0.2712
2.5246 0.8476 0.0316 0.2470
2.7130 0.8572 0.0252 0.2245
2.8826 0.8648 0.0207 0.2058
3.0000 0.8695 0.0182 0.1936

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新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.

22樓2012-01-05 22:56:01

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

zl030091126

23樓2012-01-06 03:22:38

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

phoenix708

鐵桿木蟲 (著名寫手)

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

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2012-01-05 22:56:01:
[code]
odefun=@(t,x,dx)[dx(1)-1.98*x(3)*(1-x(1))+0.73*x(1)*x(2)
dx(2)-(dx(1)-4.6*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)+ ...

Thank you!

回復此樓

24樓2012-01-06 09:07:30

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

笨笨小熊熊

25樓2012-03-06 11:26:07

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

jackzhang19

銀蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 636.9
帖子: 71
在線: 38.4小時
蟲號: 1443126
注冊: 2011-10-14
性別: GG
專業(yè): 功能陶瓷

★ ★
cenwanglai(金幣+2): ~~ 2012-03-10 08:43:18

function dy=jack(t,y)
dy=zeros(6,1);
dy(1)=1.98*y(5)*(1-y(1))-0.73*y(1)*y(3);
dy(3)=y(2)-0.00207*y(3)/(1+y(3)+2*y(5))-0.59*(y(4)*(1+y(3)+2*y(5))-y(3)*(y(4)+2*y(6)))/(1+y(3)+2*y(5))^2;
dy(5)=1/2*(-y(2)-0.00056*y(5)/(1+y(3)+2*y(5))-0.086*(y(6)*(1+y(3)+2*y(5))-y(5)*(y(4)+2*y(6)))/(1+y(3)+2*y(5))^2);
[t,y] = ode23s(@jack,tspan,[0 ,1.98,0 ,1.0067796593,1,-0.9805981300402]);
plot(t,y(:,1),'r',t,y(:,3),'b',t,y(:,5),'y')

贊一下

回復此樓

路漫漫其修遠兮吾將上下而求索

26樓2012-03-07 17:26:52

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

jackzhang19

銀蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 636.9
帖子: 71
在線: 38.4小時
蟲號: 1443126
注冊: 2011-10-14
性別: GG
專業(yè): 功能陶瓷

初值中的對x,y,s在零點的導數值設置的比較粗糙你可以再修改一下結果會更好

贊一下

回復此樓

路漫漫其修遠兮吾將上下而求索

27樓2012-03-07 17:32:11

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

中和達人

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1443.4
紅花: 4
帖子: 128
在線: 36小時
蟲號: 1633479
注冊: 2012-02-22
性別: MM
專業(yè): 工程熱力學

都是牛人啊

回復此樓

務實求新

28樓2012-03-09 11:54:00

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

相關版塊跳轉我要訂閱樓主 phoenix708 的主題更新

返回列表

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

[考研] 275求調劑 +8	Micky11223 2026-03-25	10/500	2026-03-26 15:46 by guosr9609
[考研] 297求調劑 +5	田洪有 2026-03-26	5/250	2026-03-26 15:43 by 無際的草原
[考研] 0703化學/290求調劑/本科經歷豐富/工科也可 +6	丹青奶蓋 2026-03-26	6/300	2026-03-26 15:04 by anny19840123
[考研] 一志愿河工大 081700 276求調劑 +4	地球繞著太陽轉 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 274求調劑 +13	顧九笙要謙虛 2026-03-24	19/950	2026-03-26 13:30 by 小小麥片
[考研] 312求調劑 +5	上岸吧ZJY 2026-03-22	8/400	2026-03-26 13:12 by 公瑾逍遙
[考研] 303求調劑 +6	藍山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 上海電力大學材料防護與新材料重點實驗室招收調劑研究生（材料、化學、電化學，環(huán)境） +4	我愛學電池 2026-03-23	4/200	2026-03-25 00:59 by 1027_324
[考研] 材料調劑 +6	匹克i 2026-03-23	6/300	2026-03-24 21:09 by greychen00
[考研] 一志愿北化315 求調劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 070300化學求調劑 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調劑】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 種大樹
[考研] 一志愿山東大學藥學學碩求調劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 化學308分求調劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 石河子大學（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 324求調劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 材料學碩301分求調劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] A區(qū)線材料學調劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿西南交通專碩材料355 本科雙非求調劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky