Python+NumbaでWikipediaのList of chaotic mapsを１つずつ描く（その3） Bogdanov map

さて今日は

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

の中の、

Bogdanov mapだ。

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

x(n+1)=x(n)+y(n+1)

y(n+1)=y(n)+ε*y(n)+k*x(n)*(x(n)-1)+μ *x(n)*y(n)

という写像。っても何やってもWikipediaのような図にならない、、、

ということで

https://arxiv.org/ftp/chao-dyn/papers/9402/9402006.pdf

も参考にして初期値をいろいろ変えてプロットしてみる。

おお、それっぽい図になった。パラメータはε=0,k=1.2,μ=0.0。

別のパラメータはε=0.0001,k=1.44,μ=-0.1。

リストはこちら：

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

ep=0.0
kp=1.2
up=0.0

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

    for i in range(N-1):
        y[i+1]=(1+ep)*y[i]+kp*x[i]*(1-x[i])+up*x[i]*y[i]
        x[i+1]=x[i]+y[i+1]

N=300
xmax=2
xmin=-2
ymax=2
ymin=-2
Nx=100
Ny=100

fig=plt.figure(figsize=(10,10))
plt.xlim(-0.5,2)
plt.ylim(-1,1)

for i in range(Nx):
    for j in range(Ny):
        x=np.zeros(N)
        y=np.zeros(N)
        x[0],y[0]=xmin+(xmax-xmin)*i/Nx, ymin+(ymax-ymin)*j/Ny
        Bogdanov_calc(x,y,N)
        plt.scatter(x,y,s=1,c='black',alpha=0.05)

plt.show()