« シャープのMZシリーズ(8ビットパソコン)で描かれていた3次元グラフをGoogle検索のグラフ表示機能で描く。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2020年5月18日 Microwave Journal の特集はQualcommの通信における歴史、Microwave Magazineの特集は3Dプリンタで作るフィルタ、QorvoのBoostシリーズフィルタ、Rogersのミリ波材料、など。 »

2020年5月18日 (月)

Python+numbaでバーニングシップフラクタルを描いてみる。

マンデルブロ集合、ブッダブロ、の次にメジャーなのはこれかな。

バーニングシップ・フラクタル。

https://ja.wikipedia.org/wiki/%E3%83%90%E3%83%BC%E3%83%8B%E3%83%B3%E3%82%B0%E3%82%B7%E3%83%83%E3%83%97%E3%83%BB%E3%83%95%E3%83%A9%E3%82%AF%E3%82%BF%E3%83%AB

カラーマップを変えて描いてみた。

afmhot

Burningship_afmhot

Greys

Burningship_greys

Spectral

Burningship_spectral

これが一番船が燃えてる感があるな。

 

ソースリストはこちら:

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


@jit
def burning_calc(imag, Nx, Ny, xmin, xmax, ymin, ymax, max_iteration):

    x=np.linspace(xmin,xmax,Nx)
    y=np.linspace(ymin,ymax,Ny)

    for i in range(Nx):
        for j in range(Ny):
            iteration=0
            c=x[i]+1j*y[j]
            z=0+0j
            while z.real**2+z.imag**2<4 and iteration < max_iteration:
                z=(abs(z.real)+1j*abs(z.imag))**2+c
                iteration=iteration+1
            imag[j][i]=np.log(iteration)

Nx=2000
Ny=2000
max_iteration=1000
imag=np.zeros((Nx,Ny))
xmax=-1.7
xmin=-1.8
ymax=0.02
ymin=-0.1
burning_calc(imag, Nx, Ny, xmin, xmax, ymin, ymax,max_iteration)
plt.figure(figsize=(15,15))
plt.imshow(imag,extent=(xmin,xmax,ymin,ymax),cmap=plt.cm.afmhot)
# plt.colorbar()
plt.show()

 

 

 

 

 

« シャープのMZシリーズ(8ビットパソコン)で描かれていた3次元グラフをGoogle検索のグラフ表示機能で描く。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2020年5月18日 Microwave Journal の特集はQualcommの通信における歴史、Microwave Magazineの特集は3Dプリンタで作るフィルタ、QorvoのBoostシリーズフィルタ、Rogersのミリ波材料、など。 »

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

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

« シャープのMZシリーズ(8ビットパソコン)で描かれていた3次元グラフをGoogle検索のグラフ表示機能で描く。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2020年5月18日 Microwave Journal の特集はQualcommの通信における歴史、Microwave Magazineの特集は3Dプリンタで作るフィルタ、QorvoのBoostシリーズフィルタ、Rogersのミリ波材料、など。 »

最近の記事

最近のコメント

2021年9月
      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    
フォト
無料ブログはココログ