作者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