Re: [投稿]經典圖型-物理-Chaos- Marion 5th ed-fi …

看板Mathematica作者chungyuandye (養花種魚數月亮賞星星)時間13年前 (2011/07/07 18:21)推噓2(2推 0噓 3→)

留言5則, 1人參與討論串1/1

※ 引述《kichigop (超痛恨不求甚解...XXX)》之銘言： : 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"}] 稍微做一下修改，有些事情用動畫比較能夠讓人了解 Chaos[inita_, initb_, iter_] := Block[{ita = inita, itb = initb, n = iter, recip, \[Alpha]}, \[Alpha] = 2.5; recip[x_] := \[Alpha] x (1 - x^2); ListLinePlot[ Transpose@ NestList[{recip[#[[1]]], recip[#[[2]]]} &, {ita, itb}, iter], PlotStyle -> {{Blue, Thickness[.003]}, {Red, Thickness[.002]}}, Frame -> True, BaseStyle -> {FontFamily -> "Helvetica", FontSize -> 15, FontWeight -> "Bold"}, PlotRange -> {{0, 100}, {0, 1.1}}] ] Manipulate[ mylegend[Chaos[0.700000000, 0.7000000001, x], {"0.700000000", "0.7000000001"}], {{x, 2, "Time"}, 2, 100, 1}] -- 養花種魚數月亮賞星星 http://chungyuandye.twbbs.org -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 120.119.96.8