作者LittleSnow19 (这是什麽??)
看板NTUNL
标题[问题] 18.7
时间Sun May 14 21:19:17 2006
(Mathematica code)
Duffing = {x'[t] == y[t], y'[t] == -0.25 y[t] + x[t] - x[t]^3 + 0.3 Cos[t]}
//定义 Duffing 这个微分方程组系统。并设定一些基本参数。
sol = NDSolve[{Duffing, x[0] == 1, y[0] == 0.1}, {x, y}, {t, 0, 100},
MaxSteps -> 10000]
//利用 numerical solve 微分方程组，并设定最大运作数
ParametricPlot[Evaluate[{x[t], y[t]} /. sol], {t, 0, 100}]
//利用参数式 plot 相图。
//然後在设定其他参数就可以观察是否出现非周期解的现象了!!
有没有人可以告诉我
为什麽一定要写 Evaluate[{x[t], y[t]} /. sol]
而不是 Evaluate[sol]
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.240.139
1F：→ LittleSnow19:这篇贴太快了 05/14 21:29
2F：→ LittleSnow19:也许该多想一下的 05/14 21:30
3F：→ LittleSnow19:我想是因为 sol=NDSolve[...] 长这样子的关系 : 05/14 21:30
4F：→ LittleSnow19:{{x -> InterpolatingFunction[{{0., 100.}}, <>], 05/14 21:32
5F：→ LittleSnow19: y -> InterpolatingFunction[{{0., 100.}}, <>]} 05/14 21:33
6F：推 ANUBISANKH:应该是的，不过用习惯也没有想这麽多 XD 05/14 22:23