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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3113)
>虫友互识 (233)
>学术会议 (178)
>文献求助 (83)
>导师招生 (55)
>休闲灌水 (51)
>基金申请 (47)
>博后之家 (30)
>考研 (27)
>硕博家园 (26)
>论文投稿 (26)
>教师之家 (22)
>考博 (21)
>论文道贺祈福 (20)
>找工作 (17)
>公派出国 (15)

返回列表

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

[交流] 【讨论】大家看看这是怎么了？？？？？？？？？已有1人参与

function dzdt=jisuan(t,x)
dzdt=zeros(3,1);
dx(1)=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx(2)=-10*x(2)+10*x(3);
dx(3)=-x(1)*x(2)+28*x(2)-x(3);

[t,y]=ode45('jisuan',[0,10],[0.3,0.8,1.3])

t =

0
0.2500
0.5000
0.7500
1.0000
1.2500
1.5000
1.7500
2.0000
2.2500
2.5000
2.7500
3.0000
3.2500
3.5000
3.7500
4.0000
4.2500
4.5000
4.7500
5.0000
5.2500
5.5000
5.7500
6.0000
6.2500
6.5000
6.7500
7.0000
7.2500
7.5000
7.7500
8.0000
8.2500
8.5000
8.7500
9.0000
9.2500
9.5000
9.7500
10.0000

y =

0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
错误出在哪？？？

回复此楼

» 猜你喜欢

1楼 2010-06-04 22:51:54

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

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
robert2020(金币+1):多谢应助！ 2010-06-04 23:20:47
zzuwangshilei(金币+1):多谢应助，还望详细说明下 2010-06-05 10:58:58

function test
clear
clc
x0=eps*ones(1,3);
[t,x]=ode45(@OdeFun,[0:0.02:8],x0);
plot(t,x(:,1),'r.-',t,x(:,2),'b.-',t,x(:,3),'k.-')

function dx = OdeFun(t,x)
dx1=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx2= -10*x(2)+10*x(3);
dx3=-x(1)*x(2)+28*x(2)-x(3);
dx=[dx1;dx2;dx3];

赞一下(3人)