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

今週のフラクタル2 (2z^3/(z^2+5iz+0.68)+c)

    どうも、108Hassiumです。

    先週から手抜き投稿用の企画を始めて見たはいいものの、加減がわからないせいで結構なボリュームの記事になってしまい、結果的に全然楽になってない気がします。

    今週も結構長めになりましたが、次回あたりから画像もテキストも減っていくと思います。

    2z^3/(z^2+5iz+0.68)+c

    画像
    ☝2z^3/(z^2+5iz+0.68)+cのマンデルブロ集合(x=-2~2,y=-2.9~1.1,z_0=0)

    2z3z2+5iz+0.68+c\frac{2z^3}{z^2+5iz+0.68}+cの臨界点は0、0.2i、-10.2iの3つですが、今回は主にz0=0z_0=0から生成されるフラクタル図形のみを扱います。

    画像
    ☝2z^3/(z^2+5iz+0.68)+cのマンデルブロ集合(x=-0.2~0.2,y=-2.2~-1.8,z_0=0)
    画像
    ☝2z^3/(z^2+5iz+0.68)+cのマンデルブロ集合(x=1.1~1.5,y=-1.1~-0.7,z_0=0)
    画像
    ☝2z^3/(z^2+5iz+0.68)+cのマンデルブロ集合(x=1.1~1.5,y=-1.1~-0.7,z_0=0)

    前回の2z3z2+4iz+0.27+c\frac{2z^3}{z^2+4iz+0.27}+cと同じように細部が崩壊していますが、こっちの方がより激しく崩れているような感じがします。

    ※☟前回

    画像
    ☝c=0.32+1.03i

    収束領域内に奇妙な模様のあるジュリア集合です。

    陰影が付いているので立体的にも見えます。

    画像
    ☝c=1.66+0.12i
    画像
    ☝c=1.39-0.86i

    z=0z=0という臨界点は2重解なので、同じく2重臨界点を持つz3+cz^3+cのジュリア集合と似たような三角形状の収束領域が現れることがあります。

    画像
    ☝c=-2.8i

    環状に繋がった飛び地型ジュリア集合です。

    周期発散関数のジュリア集合っぽい見た目ですが、2z3z2+5iz+0.68+c\frac{2z^3}{z^2+5iz+0.68}+cの次数は1なので周期発散関数ではありません。

    画像
    ☝c=0.14-2.06i
    画像
    ☝c=0.1-2.08i

    収束領域が環状に並び、その中に別サイクル(z0=0.2iz_0=0.2i)の収束領域があるパターンのジュリア集合です。

    画像
    ☝c=0.13-2.07i
    画像
    ☝c=1.42-0.9i
    画像
    ☝c=1.28-0.92i

    環の内側の白い領域が花びら状に枝分かれしているパターンです。

    c=1.28-0.92iの方は、なんと1440周期の480分岐です。

    画像
    ☝c=0.15-2.06i
    画像
    ☝c=0.21-2.02i
    画像
    ☝c=1.09+0.76i

    収束領域が環状に並んでいるものの、環の中に白い領域が無いパターンです。

    画像
    ☝c=0.16-2.05i

    花弁状の収束領域の周りに小さい収束領域が散らばっています。

    「環の中の白領域が枝分かれするやつ」と似ていますが、花もその外側もどちらも同じサイクルの収束領域になっています。

    画像
    ☝c=1.24-0.9i

    2種類の収束領域が分岐点を共有し、なおかつ一方の領域(白い方)がもう一方の領域の中に埋まっている形のジュリア集合です。

    実は前回の2z3z2+4iz+0.27+c\frac{2z^3}{z^2+4iz+0.27}+cでも似たようなジュリア集合が見つかっている(未掲載)のですが、第3回高周期ジュリア集合選手権のために実施したz3z+0.1i+c\frac{z^3}{z+0.1i}+cの調査では1つも見つかりませんでした。

    他の関数でもあまり見かけないので、個人的にこの形のジュリア集合は結構レアだと認識しています。

    画像
    ☝c=0.89+0.92i
    画像
    ☝c=1.33+0.5i

    よくわからないけど何となくレアっぽい感じがするものです。

    画像
    ☝c=0.16-1.41i(297周期)
    画像
    ☝c=-2.27i(360周期)
    画像
    ☝c=0.5+1.03i(404周期)
    画像
    ☝c=0.98-1.04i(672周期)
    画像
    ☝c=1.1+0.7i(1683周期)

    cの実部と虚部の小数点以下が2桁までで周期が長いものTOP5です。

    1位の1.1+0.7iは小数点以下が1桁で、2桁レギュレーションで1桁が1位になるのは(未公開の記録も含めて)史上初でした。

    摂動

    前回の2z3z2+4iz+0.27+c\frac{2z^3}{z^2+4iz+0.27}+cと2z3z2+5iz+0.68+c\frac{2z^3}{z^2+5iz+0.68}+cは、どちらも2z3z2+az+b+c\frac{2z^3}{z^2+az+b}+cという形になっています。

    aaとbbを実数にすると、y=2x3x2+ax+b+cy=\frac{2x^3}{x^2+ax+b}+cのグラフは以下のようになります。

    画像
    ☝y=2x^3/(x^2+ax+b)+cのグラフ

    このグラフ、どうやら∣a∣|a|を大きくするか∣b∣|b|を小さくしていくとy=2x+2a2x+a−2a+cy=2x+\frac{2a^2}{x+a}-2a+cのグラフと似た形になっていくようです。

    画像
    ☝a=6
    画像
    ☝b=0.2

    ただし、似た形になるといってもx=0x=0の近くに特異点(極)が存在することはかわらず、関数としての性質が2x+2a2x+a−2a+c2x+\frac{2a^2}{x+a}-2a+cと大きく異なるような範囲が残り続けます。

    また、前回と今回で「似たような形のジュリア集合を生成する関数」として度々触れているz3z+0.1i+c\frac{z^3}{z+0.1i}+cにも似たような性質があり、aaを実数とするとy=x3x+a+cy=\frac{x^3}{x+a}+cのグラフは∣a∣|a|が小さいほどy=x2+cy=x^2+cに近づいていきます。

    画像
    ☝y=x^3/(x+a)+cとy=x^2+cのグラフ

    似た関数どうしは大抵似たような形のジュリア集合を生成するのですが、特異点の有無というのはジュリア集合に少なからず影響するようです。

    なのでz3z+a+c\frac{z^3}{z+a}+cのジュリア集合は、全体的なシルエットはz2+cz^2+cのものと似ていながら穴が開いていたりバッキバキに割れたり、別サイクルの収束領域が食い込んでいたりといった感じで細部は大きく異なります。

    2z3z2+az+b+c\frac{2z^3}{z^2+az+b}+cも同様に、y=2x+2a2x+a−2a+cy=2x+\frac{2a^2}{x+a}-2a+cのジュリア集合を崩壊させたようなジュリア集合を生成します。(y=2x+2a2x+a−2a+cy=2x+\frac{2a^2}{x+a}-2a+cのジュリア集合はz2+cz^2+cのものとも似た形です)

    こんな感じで普通(?)のジュリア集合の関数に細工をして、ジュリア集合の細部を変形させたものを"perturbated julia set"と呼ぶらしいです。

    ちなみに"perturbation"は日本語では「摂動」といい、大きな作用による運動(例えば太陽からの引力による地球の公転)が小さな作用(例えば月の引力)によって変化させられる現象を指す物理学・数学用語だそうです。

    ジュリア集合とは別に、マンデルブロ集合の描画法にも摂動論を使った手法があるようです。(私は試したことないですが)

    話をジュリア集合に戻します。

    面白いジュリア集合を生成したい場合、摂動の影響は小さすぎても大きすぎても良くないようです。

    影響が小さいと変化も小さくなりますが、大きすぎる場合も何故か普通な感じの見た目になってしまいます。

    画像
    ☝摂動が大きすぎる例:z^3/(z+i)+0.2+0.96i

    2x3x2+ax+b+c\frac{2x^3}{x^2+ax+b}+cについて、∣b∣|b|が小さいほど2x+2a2x+a−2a+c2x+\frac{2a^2}{x+a}-2a+cに近くなっていく(=摂動が小さくなる)ことはすでに述べました。

    2z3z2+4iz+0.27+c\frac{2z^3}{z^2+4iz+0.27}+cは2z3z2+5iz+0.68+c\frac{2z^3}{z^2+5iz+0.68}+cと比べると面白いジュリア集合が少なかったのですが、これはaaの値(4iと5i)が近いのに対してbbの値(0.27と0.68)が小さいのが原因であると考えることができます。


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

    あなたへのおすすめ