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

登录注册

应《网络安全法》要求，自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
错误出在哪？？？

回复此楼

» 猜你喜欢

调剂已经有6人回复
280求调剂已经有6人回复
302求调剂一志愿华中师范大学已经有5人回复
346求调剂已经有4人回复
材料专硕322分已经有9人回复
372分材料与化工（085600）一志愿湖南大学求调剂已经有5人回复
318求调剂，计算材料方向已经有10人回复
310求调剂已经有14人回复
315求调剂已经有3人回复
265求调剂已经有13人回复

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