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

今週のフラクタル27 (-(c^2-3)/(c^2z^3-3z)+1 他)

    どうも、108Hassiumです。

    今回は−c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1に関するフラクタル図形をお届けします。

    ※2024/01/01追記

    記事タイトルを含むほぼすべての数式が間違っていたので修正しました。

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

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

    −c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1という関数は、1→0→∞→1…という発散サイクルを持つ3周期発散関数です。

    臨界点は±1c\pm\frac{1}{c}ですが、今回はz0=1cz_0=\frac{1}{c}をメインとします。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.54+0.38iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.52+0.44iのジュリア集合

    3周期発散関数の特徴である三つ葉穴があるジュリア集合です。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.55+0.33iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-1.71+0.34iのジュリア集合

    三つ葉穴の無いジュリア集合です。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-0.44+0.19iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.95+0.91iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.95+0.91iのジュリア集合

    以前の記事で「周期発散関数のジュリア集合の見た目が極の位数と関係しているかもしれない」という話をしましたが、発散周期が3以上だったり臨界点が複数あったりすると関係がはっきりしなくなるようです。(−c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1は1位の極しか持ちませんが、上の3つのジュリア集合は位数2以上の極があるときの特徴が出ているように見えます)

    ※☟極の位数とジュリア集合の見た目の関係について書いた記事

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=2.1のジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.15+0.09iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-0.16+0.14iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.3+0.2iのジュリア集合

    白領域のあるジュリア集合です。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.88iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-0.02+0.77iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-0.15+0.3iのジュリア集合

    3周期発散関数なのに何故か2周期発散関数のような双葉型の穴があるジュリア集合です。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.44+0.24iのジュリア集合(140周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.47+0.46iのジュリア集合(159周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.82+0.1iのジュリア集合(160周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-1.82+0.1iのジュリア集合(160周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=1.66+0.34iのジュリア集合(174周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-1.66+0.34iのジュリア集合(174周期)
    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=-0.12+0.23iのジュリア集合(223周期)

    いつものやつです。

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

    ※con(z)\text{con}(z)はzzの複素共役

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

    以前の記事で「f(z)f(z)が周期発散関数でも、f(con(z))f(\text{con}(z))は周期発散関数にならないことがある」という話をしましたが、f(z)=−c2−3c2z3−3z+1f(z)=-\frac{c^2-3}{c^2z^3-3z}+1のときはちゃんと周期発散性を維持できるようです。

    ※☟f(con(z))f(\text{con}(z))の周期発散性の話をした記事

    どうやら発散サイクルの∞以外の値の中に実数でない値が混ざっていると、con関数によって発散サイクルが破壊されてしまうようです。

    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=-0.96+0.29iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=0.39+0.33iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=0.45+0.13iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=-0.41+0.1iのジュリア集合

    三つ葉穴はcの値によって現れたり現れなかったりしますが、現れた場合には−c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1では見られなかった「三つ葉型の接点と双葉型の接点を両方持つ穴」が出現することがあります。(c=-0.41+0.1iの右側がわかりやすいです)

    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=1.55+0.12iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=1.74+0.1iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=1.84+0.24iのジュリア集合

    収束領域の配置にcon系関数っぽさが強く出ているジュリア集合です。

    ※☟con系関数のジュリア集合の特徴の説明がある記事

    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=-0.26+0.25iのジュリア集合(126周期)
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=1.73+0.09iのジュリア集合(133周期)
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=-1.73+0.09iのジュリア集合(133周期)
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=-0.36+0.3iのジュリア集合(186周期)
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=1.86+0.23iのジュリア集合(324周期)
    画像
    ☝-(c^2-3)/(c^2con(z)^3-3con(z))+1のc=0.09+0.78iのジュリア集合(736周期)

    いつものです。

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

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

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

    B(z)B(z)を含む関数のマンデルブロ集合としては珍しく、点対称な形になりました。

    どうもf(z)=cg(z)f(z)=cg(z)(g(z)g(z)はccに依存しない関数)のときはf(B(z))f(B(z))のマンデルブロ集合が点対称になるっぽいのですが、今回はこのパターンに当てはまらないので点対称になるのは予想外でした。

    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=0.45のジュリア集合
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-0.3+0.1iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=0.1のジュリア集合

    一応「1点に穴が3つずつ集まるジュリア集合」はあるようですが、−c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1のものとはかけ離れた見た目になります。

    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=0.1+0.71iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=1.82+0.2iのジュリア集合
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=1.79+0.03iのジュリア集合

    いかにも「絶対値を含む関数のジュリア集合」っぽい見た目のジュリア集合です。

    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=1.86+0.09iのジュリア集合(25周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-0.84+0.68iのジュリア集合(26周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-0.85+0.66iのジュリア集合(26周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-0.87+0.78iのジュリア集合(28周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-1.11+0.54iのジュリア集合(28周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-0.17+1.04iのジュリア集合(36周期)
    画像
    ☝-(c^2-3)/(c^2B(z)^3-3B(z))+1のc=-1.63+0.35iのジュリア集合(88周期)

    アレです。

    周期発散と穴の形状

    3周期発散関数のジュリア集合の穴の構造は、

    • 三つ葉型:3つの穴が1点でつながっている箇所がある

    • 飛び地型:どの2つの穴の間にも別の穴か収束領域があり、点でつながっている穴はない

    ・・・の2種類しかないというのが今までの記事での定説でした。

    ところが、前述の通り−c2−3c2z3−3z+1-\frac{c^2-3}{c^2z^3-3z}+1ではどちらでもない双葉型の穴のあるジュリア集合が見つかってしまいました。

    画像
    ☝-(c^2-3)/(c^2z^3-3z)+1のc=0.88iのジュリア集合(再掲)

    そもそもなぜ3周期発散関数のジュリア集合があの2種類だけだと考えていたのかというと、有理関数のジュリア集合の収束領域の形状と関連付けていたからでした。

    画像
    ☝三つ葉型と飛び地型の収束領域を持つジュリア集合

    3周期のジュリア集合の収束領域が3つ葉型と飛び地型の2種類だから3周期発散関数の穴もそれと同じで2種類なんだろう、と考えていたのでした。

    ならあの双葉穴は何なんだと考えていたところ、あるジュリア集合の存在を思い出しました。

    画像
    ☝z^3/(z+0.1i)+0.38-0.22iのジュリア集合(11周期)

    ※☟この記事の最後の方に出てきます

    白領域を取り囲んでいる収束領域をよく見ると、2個の領域が1点でつながった形になっています。

    ということは3周期でも2個の収束領域が1点でつながった形のジュリア集合が存在するのではないかと気付き、少し探してみたところ以下のようなジュリア集合が見つかりました。

    画像
    ☝z^3/(z+0.1-0.2i)-0.1+0.6iのジュリア集合

    ということは、3周期発散関数のジュリア集合に双葉型の穴を持つものが存在するのは、有理関数のジュリア集合の収束領域の形状として双葉型というパターンが存在することと関係しているだろう、と考えることができます。

    三つ葉型と飛び地型はz2+cz^2+cのような簡単な関数のジュリア集合でも見られますが、双葉型は遭遇する機会が少なく「ごく一部の関数でしか見られないレアケース」という風に認識していたため、周期発散関数の話題とリンクするのはかなり意外に感じました。

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

    あなたへのおすすめ