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

今週のフラクタル25 (c/(z^3-3z^2))

    どうも、108Hassiumです。

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

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

    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(z_0=2,x=-10~-2)
    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(z_0=2,x=-16~4)

    cz3−3z2\frac{c}{z^3-3z^2}は0→∞→0というサイクルを持つ2周期発散関数ですが、czn−1+1\frac{c}{z^n-1}+1等で見られた網目状の模様が現れません。

    ※☟「周期発散」という概念の説明がある記事

    ※☟czn−1+1\frac{c}{z^n-1}+1に関する記事

    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(3.2+3.5i付近)
    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(-8.901001+3.537224i付近)
    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(-5.0607+0.100201i付近)
    画像
    ☝c/(z^3-3z^2)のマンデルブロ集合(-8.5951+0.0089i付近)

    マンデルブロ集合の細部を拡大すると、z2+cz^2+cのマンデルブロ集合やジュリア集合に似た模様が見られます。

    やはりczn−1+1\frac{c}{z^n-1}+1のときのような細かい網目模様はなく、z2+cz^2+cのマンデルブロ集合でも見られるような模様に謎の輪っかがついていたり細い線でつながっていたりといった感じです。

    画像
    ☝(-3.47+0.61i)/(z^3-3z^2)のジュリア集合
    画像
    ☝(-4.88+0.55i)/(z^3-3z^2)のジュリア集合
    画像
    ☝(-6+1.93i)/(z^3-3z^2)のジュリア集合
    画像
    ☝(-6.85+2.63i)/(z^3-3z^2)のジュリア集合

    ジュリア集合もマンデルブロ集合と同じ感じで、網目模様が現れなかったり、網状になっても線が細かったりと目が非常に粗かったりして今までに紹介した周期発散関数とはかなり様子が異なります。

    画像
    ☝(-8.66+0.04i)/(z^3-3z^2)のジュリア集合(495周期)
    画像
    ☝(-8.11+2.16i)/(z^3-3z^2)のジュリア集合(860周期)
    画像
    ☝(-8.07+2.18i)/(z^3-3z^2)のジュリア集合(924周期)
    画像
    ☝(-4.82+0.65i)/(z^3-3z^2)のジュリア集合(1496周期)
    画像
    ☝(-8.28+2.58i)/(z^3-3z^2)のジュリア集合(3240周期)

    いつものです。

    マンデルブロ集合のサイズが大きいと、長い周期のジュリア集合が多く現れがちです。

    周期発散と極の位数

    cz3−3z2\frac{c}{z^3-3z^2}は2周期発散関数でありながら周期発散関数によく見られる網状の模様が現れないのが特徴ですが、この性質はcz3−3z2\frac{c}{z^3-3z^2}が持つ「極きょく」の「位数」というものと関係している可能性があります。

    まず極というのは関数の特異点の一種で、有理関数の場合は分母が0になって関数の値が無限大に発散するような点のことです。

    cz3−3z2\frac{c}{z^3-3z^2}の場合はz=0z=0とz=3z=3のときにz3−3z2z^3-3z^2が0になってcz3−3z2\frac{c}{z^3-3z^2}が発散するので、この2点が極になります。

    そして極の位数というのは極の性質を表す値で、aaがf(z)f(z)の極であるとき(z−a)nf(z)(z-a)^nf(z)がz=az=aで極を持たなくなるような最小のnnが「aaにおける極の位数」だそうです。

    f(z)=cz3−3z2f(z)=\frac{c}{z^3-3z^2}の場合、z2f(z)z^2f(z)でz=0z=0の極が無くなるのでz=0z=0での位数は2、z=3z=3の方は(z−3)1f(z)(z-3)^1f(z)で消えるので位数は1となります。

    画像
    画像

    周期発散関数は∞を含むサイクルを持つ関数ですが、このサイクルには極の値も含まれています。

    2周期発散関数だと極と無限遠を往復するサイクルになり、3周期以上の場合でも極から無限遠に飛ばされてからまた極の場所に戻ってくる、という構造になっています。

    cz3−3z2\frac{c}{z^3-3z^2}での0→∞→0というサイクルも極と無限遠を行き来するサイクルですが、前述の通りz=0z=0における極の位数は2です。

    どうも「発散サイクルに含まれる極の位数」がジュリア集合やマンデルブロ集合の見た目と関係しているようで、位数が1のときはczn−1+1\frac{c}{z^n-1}+1のような網状の見た目に、2以上だと今回のように糸状になるっぽいです。

    例えばcz4−2z3\frac{c}{z^4-2z^3}のz=0z=0の極の位数は3ですが、この関数でもcz3−3z2\frac{c}{z^3-3z^2}と同様に糸状の模様が現れます。

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

    また、cz3−3z2+3\frac{c}{z^3-3z^2}+3という関数は、定数項があることによって位数1の方の極から発散していく(=3→∞→3というサイクルを持つ)ようになっており、czn−1+1\frac{c}{z^n-1}+1と同じような網状の模様が現れます。

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

    f(z)=az3−3z2+bf(z)=\frac{a}{z^3-3z^2}+bという形をした3周期発散関数は、2つの極のうちどちらを発散サイクルに含めるかで2種類に分かれます。

    画像
    ☝(-c^4+3c^3)/(z^3-3z^2)+cのマンデルブロ集合(0→∞→c→0)
    画像
    ☝(-c^4+3c^3)/(z^3-3z^2)+cのマンデルブロ集合(3→∞→c→3)

    ご覧の通り、位数2の極をサイクルに含む方に糸状の模様が現れています。

    ちなみに過去の記事で紹介した周期発散関数はほとんどが位数1の極しか持ちませんが、1z2−1\frac{1}{z^2}-1は位数2の極を持つ3周期発散関数であり、ジュリア集合は糸状(鎖状?)になっています。

    画像
    ☝1/z^2-1のジュリア集合

    こういった現象が発生する仕組みは不明(そもそも感覚に基づく非数学的現象なので解明のしようがない)ですが、何か数学的な背景があるとすれば臨界点の多重度の話と関係あるんじゃないか、と思っています。

    ※☟臨界点の多重度の話が載っている記事


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

    あなたへのおすすめ