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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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