[投稿]經典圖型-物理-Chaos- Marion 5th ed-fig.4-24, p.175

看板Mathematica作者 (超痛恨不求甚解...XXX)時間13年前 (2011/07/07 16:55), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
Legend部分現學現賣用戴老師教給大家的 但是顏色有點怪怪的(還請老師指教),程式部分相當拙劣,是剛學Mathematica時寫的 請大家多多包涵 不過重點是畫完圖我發現書上的兩條線好像應該對調 這一小節主要是在講Chaos的產生與initial conditions有很大的關係(蝴蝶效應) 並且在這個example 4.3中利用一個noliear eqation給予不同的I.C.來驗證 結果得到fig. 4-24。 (分隔線以下直接複製即可...我是用Mathematica 7.0版) =================================================================== (*The two plots in Marion's textbook should be exchanged!!!*) (* The nonlinear equation is: Subscript[x, n+1]=\[Alpha] Subscript[x, n](1-Subscript[x, n]^2) \ *) mylegend[plot_Graphics, legend_List] := Block[ {p = plot, l = legend, color, temp}, color = Cases[p, Hue[a_, b_, c_] :> Hue[a, b, c], Infinity]; temp = {color[[#]], l[[#]]} & /@ Range[Length@color]; Labeled[p, Grid[{Graphics[{#[[1]], Thickness[0.1], Line[{{0, 0}, {1, 0}}]}, ImageSize -> {24, 24}, AspectRatio -> 24/24, ImagePadding -> 0], #[[2]]} & /@ temp], {Top}] ] (*以上為戴老師的Legend*) \[Alpha] = 2.5; recip[x_] := \[Alpha] x (1 - x^2); chaos1 = NestList[recip, 0.700000000, 50];(*blue*) recip[x_] := \[Alpha] x (1 - x^2); chaos2 = NestList[recip, 0.700000001, 50];(*red*) p1 = ListPlot[ {chaos1, chaos2}, PlotRange -> {0, 1}, Joined -> True, PlotStyle -> {{RGBColor[0, 0, 1], Thickness[.003]}, {RGBColor[1, 0, 0], Thickness[.002]}}, BaseStyle -> {FontFamily -> "Helvetica", FontSize -> 15, FontWeight -> "Bold"}]; mylegend[p1, {"0.700000000", "0.700000001"}] -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.116.91.246
文章代碼(AID): #1E5NHr0E (Mathematica)
文章代碼(AID): #1E5NHr0E (Mathematica)