終於。Orz 。睡覺

今天終於把力學take home exam 寫完交出去了。
在與周函毅通力合作下,我們成功的寫出了Poincaré!

sys = Sequence[
x''[t] == (-3Sin[x[t]] – Sin[x[t] – 2y[t]] – 2Sin[x[t] – y[t]](y'[t]^2 + 0.5x'[t]^2Cos[x[t] – y[t]]))/(1.5 – 0.5Cos[2x[t] – 2y[t]]),
y''[t] == 2Sin[x[t] – y[t]](x'[t]^2 + 2Cos[x[t]] + y'[t]^2Cos[x[t] – y[t]])/(3 – Cos[2x[t] – 2y[t]])];
sol1 = NDSolve[{sys, x[0] == 0.5π, x'[0] == 0, y[0] == 0, y'[0] == 0}, {x[t], x'[t], y[t], y'[t]}, {t, 0, 200}, MaxSteps -> 300000]

ParametricPlot[Evaluate[{y[t], y'[t]} /. sol1], {t, 0, 200}, PlotPoints -> 1000]

Module[{temp1, temp2},
sys = Sequence[
x''[t] == (-3Sin[x[t]] – Sin[x[t] – 2y[t]] – 2Sin[x[t] – y[t]](y'[t]^2 + 0.5x'[t]^2Cos[x[t] – y[t]]))/(1.5 – 0.5Cos[2x[t] – 2y[t]]),
y''[t] == 2Sin[x[t] – y[t]](x'[t]^2 + 2Cos[x[t]] + y'[t]^2Cos[x[t] – y[t]])/(3 – Cos[2x[t] – 2y[t]])];
temp1 = NDSolve[{sys, x[0] == 0.5π, x'[0] == 0, y[0] == 0, y'[0] == 0}, {x[t], x'[t], y[t], y'[t]}, {t, 0, 18000},
MaxSteps -> 3000000] // Flatten ;
temp2 = Table[ { x[t] , x'[t]} /. temp1, {t, 0, 18000, 2Pi}]; |ListPlot[temp2, PlotStyle -> PointSize[0.01]]]
<-Poincaré

五題寫了四題。他說,數學家就是用咖啡來換取答案的。
我想這話是可信的,是值得完全接受的。
這周末兩天我只睡了三小時@@,沒有咖啡,我只能對著書本一直點頭。

若大家對Mathematica有興趣,我有找到一個範例nb網站

大家可以玩玩看雙擺的java。
我這張畫的是它垂直高度(y)與角度(θ)的關係,它會呈現混沌的現象。

這一週都一直坐在座位前面,都沒有運動,
晚上熬夜又一定要吃點東西。
量了體脂率,上升了 3%(有沒有那麼誇張?)

喔,要我用以最肥胖的樣子來面對我的生日,
真叫人情何以堪。

然後今天熬夜完,想說補補身子,中午吃了麥當勞?!XD
睡覺去。

發表迴響

你的電子郵件位址並不會被公開。 必要欄位標記為 *