看板 NTUNL 關於我們 聯絡資訊
(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
LittleSnow19:這篇貼太快了 05/14 21:29
LittleSnow19:也許該多想一下的 05/14 21:30
LittleSnow19:我想是因為 sol=NDSolve[...] 長這樣子的關係 : 05/14 21:30
LittleSnow19:{{x -> InterpolatingFunction[{{0., 100.}}, <>], 05/14 21:32
LittleSnow19: y -> InterpolatingFunction[{{0., 100.}}, <>]} 05/14 21:33
ANUBISANKH:應該是的,不過用習慣也沒有想這麼多 XD 05/14 22:23