| 查看: 202 | 回复: 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 错误出在哪??? |
» 猜你喜欢
一志愿0817化学工程与技术,求调剂
已经有8人回复
271分求调剂学校
已经有3人回复
285求调剂
已经有4人回复
生物学308分求调剂(一志愿华东师大)
已经有7人回复
296材料专硕求调剂
已经有18人回复
材料与化工306分找调剂
已经有18人回复
302求调剂
已经有9人回复
一志愿北京科技大学材料工程085601,求调剂
已经有11人回复
283分材料与化工求调剂
已经有19人回复
326求调剂
已经有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














回复此楼