« 扇町公園をぶらぶら。カンテレの社屋を眺めながら、、、 | トップページ | ローレンツ方程式3種類(通常のLorenzモデル、 Lorenz84モデル、Lorenz96モデル)をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算(その18) »

2021年10月13日 (水)

Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使う(その17) 4つの羽をもつFour-wing attractorを描く。ところでMatplotlib3.4.0からfig.gca(projection=”3d”)が警告が出てfig.add_subplot(projection=”3d”)にしないとだめっぽい。

久しぶりのシリーズ。

たまたまこのサイトを見た。

https://www.dynamicmath.xyz/strange-attractors/

知らないアトラクタが結構ある。。。ということでいくつか描いてみよう。

式はこちら。

Fourwing_eqjpg

描くと確かに4つ羽が(わかりにくいけど)ある!

Four_wing

わかりにくいのでGIFアニメにしてみた。クリックで始まります。

Fourwing

ところで最初、3次元のグラフを描くのに

fig = plt.figure(figsize=(12,12))
ax = fig.gca(projection='3d')

とすると、

anaconda3\lib\site-packages\ipykernel_launcher.py:37: MatplotlibDeprecationWarning: Calling gca() with keyword arguments was deprecated in Matplotlib 3.4. Starting two minor releases later, gca() will take no keyword arguments. The gca() function should only be used to get the current axes, or if no axes exist, create new axes with default keyword arguments. To create a new axes with non-default arguments, use plt.axes() or plt.subplot().


のようなWarningがでた。

https://matplotlib.org/stable/api/prev_api_changes/api_changes_3.4.0.html

から

fig = plt.figure(figsize=(12,12))
ax = fig.add_subplot(projection='3d')

こんな感じに変えないといけないっぽい。これで警告亡くなった。

« 扇町公園をぶらぶら。カンテレの社屋を眺めながら、、、 | トップページ | ローレンツ方程式3種類(通常のLorenzモデル、 Lorenz84モデル、Lorenz96モデル)をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算(その18) »

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

コメント

コメントを書く

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

« 扇町公園をぶらぶら。カンテレの社屋を眺めながら、、、 | トップページ | ローレンツ方程式3種類(通常のLorenzモデル、 Lorenz84モデル、Lorenz96モデル)をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算(その18) »

最近の記事

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