| 查看: 187 | 回复: 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 错误出在哪？？？ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有11人回复
-
求推荐英文EI期刊
已经有5人回复
-
存款400万可以在学校里躺平吗
已经有27人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
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