メインコンテンツへスキップ
見出し画像

今週のフラクタル20 ((-c^4+c^3+c-1)/(z^3-1)+c)

    どうも、108Hassiumです。

    今週は−c4+c3+c−1z3−1+c\frac{-c^4+c^3+c-1}{z^3-1}+cに関するフラクタル図形をお届けします。

    (-c^4+c^3+c-1)/(z^3-1)+c

    画像
    ☝(-c^4+c^3+c-1)/(z^3-1)+cのマンデルブロ集合(z_0=0)

    −c4+c3+c−1z3−1+c\frac{-c^4+c^3+c-1}{z^3-1}+cという関数は、1→∞→c→1というサイクルを持つ3周期発散関数です。

    −cn+1+cn+c−1zn−1+c\frac{-c^{n+1}+c^n+c-1}{z^n-1}+cは「臨界点はz=0z=0のみで、多重度はn−1n-1」「ジュリア集合はnn回回転対称」という性質を持ち、3周期発散関数におけるマルチブロのような存在です。

    画像
    ☝c=0.75+0.65iのジュリア集合
    画像
    ☝c=1.3+0.31iのジュリア集合
    画像
    ☝c=0.05+iのジュリア集合
    画像
    ☝c=0.21+0.85iのジュリア集合
    画像
    ☝c=0.89+0.87iのジュリア集合
    画像
    ☝c=1.13+0.03iのジュリア集合

    三つ葉型の穴があるジュリア集合です。

    画像
    ☝c=-0.64+0.89iのジュリア集合
    画像
    ☝c=-0.14+0.67iのジュリア集合
    画像
    ☝c=-0.15+1.07iのジュリア集合
    画像
    ☝c=-0.67+1.02iのジュリア集合
    画像
    ☝c=-0.1+0.94iのジュリア集合

    三つ葉穴がないジュリア集合です。

    画像
    ☝c=-0.51+0.88iのジュリア集合(x=-16~16,y=-16~16)
    画像
    ☝c=-0.5+0.87iのジュリア集合(x=-32~32,y=-32~32)

    ccの値が−1±3i2\frac{-1\pm\sqrt{3}i}{2}(≒-0.5±0.87i)に近ければ近いほど、−c4+c3+c−1z3−1+c\frac{-c^4+c^3+c-1}{z^3-1}+cのジュリア集合はサイズがデカくなっていくようです。

    多分c=−1±3i2c=\frac{-1\pm\sqrt{3}i}{2}のとき−c4+c3+c−1=0-c^4+c^3+c-1=0になるのと関係があって、理屈としてはc(z22+1z)c(\frac{z^2}{2}+\frac{1}{z})のジュリア集合で起きる現象と同じだと思うのですが、同じく−c4+c3+c−1=0-c^4+c^3+c-1=0の解であるc=1c=1の場合は何故か同じ現象は起きません。

    ※☟c(z22+1z)c(\frac{z^2}{2}+\frac{1}{z})のジュリア集合の大きさについて言及している記事

    画像
    ☝c=0.53+0.88iのジュリア集合(184周期)
    画像
    ☝c=0.47+0.62iのジュリア集合(198周期)
    画像
    ☝c=-0.63+0.89iのジュリア集合(343周期)
    画像
    ☝c=0.04+0.69iのジュリア集合(405周期)
    画像
    ☝c=0.47+0.49iのジュリア集合(432周期)

    いつものやつです。

    (-c^4+c^3+c-1)/(con(z)^3-1)+c

    con(x+iy)=x−iy\text{con}(x+iy)=x-iy

    画像
    ☝(-c^4+c^3+c-1)/(con(z)^3-1)+cのマンデルブロ集合(z_0=0)

    複素共役を混ぜると三つ葉穴どころか発散領域自体が消え(中央の穴とそこから6方向にのびる領域は誤判定によるものです)、全く違う見た目になりました。

    画像
    ☝c=1.2+0.54iのジュリア集合
    画像
    ☝c=1.7+0.1iのジュリア集合
    画像
    ☝c=0.4+0.2iのジュリア集合

    ジュリア集合にも、やはり3周期発散の特徴は現れません。

    画像
    ☝c=-0.8のジュリア集合
    画像
    ☝c=-0.6のジュリア集合
    画像
    ☝c=0.9のジュリア集合

    例外として、ccが実数のとき(c6c^6が実数のとき?)は3周期発散になるようです。

    画像
    ☝c=1.23+0.54iのジュリア集合
    画像
    ☝c=0.49+0.8iのジュリア集合

    3周期発散にならない場合、2種類の吸引的サイクルが存在する場合があるようです。

    画像
    ☝c=0.7+0.4iのジュリア集合
    画像
    ☝c=0.8+0.2iのジュリア集合

    不安定領域と安定領域が混在することもあるようです。

    (-c^4+c^3+c-1)/(B(z)^3-1)+c

    ※B(x+iy)=∣x∣+i∣y∣B(x+iy)=|x|+i|y|

    画像
    ☝(-c^4+c^3+c-1)/(B(z)^3-1)+cのマンデルブロ集合(z_0=0)

    右下の領域(ccの実部も虚部も正)にだけ発散領域が現れました。

    画像
    ☝c=1+0.8iのジュリア集合
    画像
    ☝c=1.2+0.53iのジュリア集合
    画像
    ☝c=1.2+0.67iのジュリア集合
    画像
    ☝c=1.24+0.1iのジュリア集合

    ccが右下にある場合でも、3周期発散っぽい見た目になるときとならないときがあるようです。

    画像
    ☝c=-0.2+0.6iのジュリア集合
    画像
    ☝c=-0.7+0.4iのジュリア集合
    画像
    ☝c=-0.6+0.3iのジュリア集合

    ccが左下にない場合はやはり収束領域自体が現れず、また−c4+c3+c−1con(z)3−1+c\frac{-c^4+c^3+c-1}{\text{con}(z)^3-1}+cと同様に2種類の安定領域が混在することもあるようです。

    画像
    ☝c=1.3+0.6iのジュリア集合

    安定領域と不安定領域と発散領域が混在するジュリア集合です。


     
     
    数学関係の記事を書きます。毎週日曜更新予定。

    あなたへのおすすめ