« Python+NumbaでWikipediaのList of chaotic mapsを1つずつ描く(その1) Gingerbreadman map。クッキーよりゼットンかキングジョーに見える。 | トップページ | 円内部(2次元)の一様乱数でアホなミスをする、、、球内部(3次元)のときだけ角度に気を付けると思っていた、、、半径もか、、、Pythonで図示して確認する。 »

2020年6月18日 (木)

Python+NumbaでWikipediaのList of chaotic mapsを1つずつ描く(その2) Tinkerbell Map。あのティンカーベルに似ている?

ネタがないのでWikipediaのList of Chaotic Mapsをいろいろ見ていくシリーズ。

https://en.wikipedia.org/wiki/List_of_chaotic_maps

 

今日はTinkerbell Mapだ。

https://en.wikipedia.org/wiki/Tinkerbell_map

x[i+1]=x[i]**2-y[i]**2+a*x[i]+b*y[i]
y[i+1]=2*x[i]*y[i]+c*x[i]+d*y[i]

こんな写像。これはとても簡単に

Tinkerbell

 

綺麗に描ける。うーん、ティンカーベルかー。

前にも書いたけど阪神高速のロゴに見える、、、

Hanshinkosoku

 

Pythonのコードはこちら:

import numpy as np
import matplotlib.pyplot as plt
from numba import jit

a=0.9
b=-0.6013
c=2.0
d=0.50

@jit
def tinker_calc(x,y,N):

    for i in range(N-1):
        x[i+1]=x[i]**2-y[i]**2+a*x[i]+b*y[i]
        y[i+1]=2*x[i]*y[i]+c*x[i]+d*y[i]

N=20000


fig=plt.figure(figsize=(10,10))
plt.xlim(-1.6,0.6)
plt.ylim(-1.6,0.6)



x=np.zeros(N)
y=np.zeros(N)
x[0],y[0]=-0.72,-0.64
tinker_calc(x,y,N)
plt.scatter(x,y,s=1)


plt.show()

« Python+NumbaでWikipediaのList of chaotic mapsを1つずつ描く(その1) Gingerbreadman map。クッキーよりゼットンかキングジョーに見える。 | トップページ | 円内部(2次元)の一様乱数でアホなミスをする、、、球内部(3次元)のときだけ角度に気を付けると思っていた、、、半径もか、、、Pythonで図示して確認する。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« Python+NumbaでWikipediaのList of chaotic mapsを1つずつ描く(その1) Gingerbreadman map。クッキーよりゼットンかキングジョーに見える。 | トップページ | 円内部(2次元)の一様乱数でアホなミスをする、、、球内部(3次元)のときだけ角度に気を付けると思っていた、、、半径もか、、、Pythonで図示して確認する。 »

最近の記事

2021年3月
  1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30 31      
フォト
無料ブログはココログ