| 查看: 241 | 回复: 1 | |||
[交流]
【讨论】大家看看这是怎么了？？？？？？？？？ 已有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 错误出在哪？？？ |
回复此楼» 猜你喜欢
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
听说今天filecode变了
已经有45人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
固定端突然变了，今天
已经有8人回复
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
小木虫(金币+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]; |
2楼2010-06-04 23:07:14
100