版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

phoenix708

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 9829.8
红花: 1
帖子: 1443
在线: 122.8小时
虫号: 476825
注册: 2007-12-13
性别: GG
专业: 催化化学

引用回帖:

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

我没有matlab的嘛，我没用过，呵呵！麻烦帮忙运算下，把数据也发上来下，我要自己用origin画图，谢谢

赞一下

回复此楼

21楼2012-01-05 21:19:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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发消息送TA红花 TA的回帖

zl030091126

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 9829.8
红花: 1
帖子: 1443
在线: 122.8小时
虫号: 476825
注册: 2007-12-13
性别: GG
专业: 催化化学

送鲜花一朵

引用回帖:

楼: 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发消息送TA红花 TA的回帖

笨笨小熊熊

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jackzhang19

银虫 (小有名气)

应助: 3 (幼儿园)
金币: 636.9
帖子: 71
在线: 38.4小时
虫号: 1443126
注册: 2011-10-14
性别: GG
专业: 功能陶瓷

★ ★
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发消息送TA红花 TA的回帖

jackzhang19

银虫 (小有名气)

应助: 3 (幼儿园)
金币: 636.9
帖子: 71
在线: 38.4小时
虫号: 1443126
注册: 2011-10-14
性别: GG
专业: 功能陶瓷

初值中的对x,y,s在零点的导数值设置的比较粗糙你可以再修改一下结果会更好

赞一下

回复此楼

路漫漫其修远兮吾将上下而求索

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

中和达人

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1443.4
红花: 4
帖子: 128
在线: 36小时
虫号: 1633479
注册: 2012-02-22
性别: MM
专业: 工程热力学

都是牛人啊

回复此楼

务实求新

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 phoenix708 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 307求调剂 +3	73372112 2026-02-28	4/200	2026-02-28 23:49 by 迷糊CCPs
[考研] 272求调剂 +3	材紫有化 2026-02-28	3/150	2026-02-28 22:52 by ms629
[考研] 304求调剂 +3	52hz~~ 2026-02-28	4/200	2026-02-28 21:41 by gaoxiaoniuma
[考研] 290求调剂 +5	材料专硕调剂； 2026-02-28	6/300	2026-02-28 21:40 by gaoxiaoniuma
[考研] 295求调剂 +5	19171856320 2026-02-28	5/250	2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4	想申博！ 2026-02-26	4/200	2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5	提神豆沙包 2026-02-28	5/250	2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3	巴拉巴拉根556 2026-02-28	3/150	2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8	南迦720 2026-02-28	8/400	2026-02-28 21:30 by gaoxiaoniuma
[考研] 材料类求调剂 +6	wana_kiko 2026-02-28	6/300	2026-02-28 21:20 by gaoxiaoniuma
[考研] 284求调剂 +4	天下熯 2026-02-28	4/200	2026-02-28 21:13 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4	好好好1233 2026-02-28	7/350	2026-02-28 20:42 by 好好好1233
[考研] 298求调剂 +8	人间唯你是清欢 2026-02-28	11/550	2026-02-28 20:26 by L135790
[考研] 276求调剂 +3	路lyh123 2026-02-28	4/200	2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10	hyf hyf hyf 2026-02-28	11/550	2026-02-28 18:50 by 无际的草原
[考研] 材料调剂 +3	爱擦汗的可乐冰 2026-02-28	3/150	2026-02-28 18:06 by houyaoxu
[考研] 寻找调剂 +3	LYidhsjabdj 2026-02-28	3/150	2026-02-28 12:59 by miniwendy
[考研] 272求调剂 +3	田智友 2026-02-28	3/150	2026-02-28 12:31 by 王加浩to
[基金申请] 面上可以超过30页吧？ +12	阿拉贡aragon 2026-02-22	13/650	2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4	N1ce_try 2026-02-24	8/400	2026-02-26 08:40 by N1ce_try