Python+numbaでマンデルブロ集合を描く(今更)、、、colormapsの色を確かめるために!
最近、家にいることが多いのでプログラムでも組んでみようと思った。Pythonなら会社でも使えるしまあいいかなということで
お絵かきから、、、で早速カラーマップで悩む、、、どれがいいんだ、、、
ということでここは実際にお絵かき、Mandelbrot集合あたりがいいだろうということでやってみた。
さすがにnumbaのjitで計算は一瞬です。
Spectral
Spectral_r
amhot
amhot_r
twilight_shifted
coolwarm_r
YlOrRd_r
うーん、afmhotあたりか、、、
プログラムはこちら。cm.????のところを変えれば色が変わる。
import numpy as np
import matplotlib.pyplot as plt
from numba import jit
@jit
def mandel_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=z**2+c
iteration=iteration+1
imag[j][i]=np.log(iteration)
Nx=1200
Ny=1200
max_iteration=1000
imag=np.zeros((Nx,Ny))
xmax=1.0
xmin=-2.5
ymax=1.5
ymin=-1.5
mandel_calc(imag, Nx, Ny, xmin, xmax, ymin, ymax,max_iteration)
plt.figure(figsize=(12,12))
plt.imshow(imag,extent=(xmin,xmax,ymin,ymax),cmap=plt.cm.coolwarm_r)
# plt.colorbar()
plt.show()
« かつやで赤辛カツ丼(竹、ご飯大盛り)+豚汁(小)を食す。 | トップページ | 新型コロナウイルス肺炎の重症患者に使われる人工呼吸器を素粒子・天体物理学者たちが共同で迅速に開発し、FDA(アメリカ食品医薬品局)の緊急使用認可を得たそうだ。Mechanical Ventilator Milano(MVM)。arXivに詳細が出ているのが物理学者らしい。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年2月7日 フジクラがMWCバルセロナでミリ波(28GHz,60GHz)モジュールデモ、Qorvo、NXP、AMDらのヘルスケアウェビナー、Samsung Galaxy S25もやはり5Gミリ波AiPは1個、u-bloxがBluetooth LEモジュール4種発表(2025.02.07)
- 高周波・RFニュース 2025年2月6日 3GPPがRel-20(5G Advancedと6G)の計画を公開、SEMCOが産業用MLCC(0805 inch, X6S, 22uF, 25V)と車載MLCC(1005,2.2uF,10V)を発表、KYOCERA AVXのスーパーキャパシタシミュレーションソフト、I向けフォトニックチップ解説(2025.02.06)
- 高周波・RFニュース 2025年2月5日 Mavenirの5G NTNホワイトペーパー、ミニサーキットの1.5mm角QFNパッケージMMIC、Keysigntがスペインに6G研究ラボ設立、アイルランドとスコットランドの大学がRIS(reconfigurable intelligent surfaces)で協力(2025.02.05)
- 高周波・RFニュース 2025年2月4日 Samsung Galaxy S25plus、S25ultra分解、やっぱり5Gミリ波AiPは1個だけ、スマートワイヤレス聴診器にNordicのICが採用、TDKのTMR角度センサ解説(2025.02.04)
« かつやで赤辛カツ丼(竹、ご飯大盛り)+豚汁(小)を食す。 | トップページ | 新型コロナウイルス肺炎の重症患者に使われる人工呼吸器を素粒子・天体物理学者たちが共同で迅速に開発し、FDA(アメリカ食品医薬品局)の緊急使用認可を得たそうだ。Mechanical Ventilator Milano(MVM)。arXivに詳細が出ているのが物理学者らしい。 »
コメント