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

登录注册

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

24小时热门版块排行榜

返回列表

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
错误出在哪？？？

回复此楼

» 猜你喜欢

售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有14人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有13人回复
越来越担心硕士毕业后没有竞争力了已经有15人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有3人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有5人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有7人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有4人回复

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人)