終於。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

睡覺去。

No related posts.