テトレーションの列
z,zz,zzz,…
の極限を zzz⋯ と書くことにします.これが ∞ に発散するような z はどのようなものであろうかと考えてみました.例えば
222⋯=∞
は当たり前ですが,
√2√2√2⋯=2
となります. z が正の実数 a∈R>0 の場合の ∞ への発散領域は初等解析的な方法によって求めることができます.正の実数 a に対して
1,a,aa,aaa,aaaa,…
という値の列は, y=ax と y=x のグラフを用いて以下のように描かれます.
このとき aaa⋯ は発散するかどうかのしきい値は, y=ax と y=x のグラフがただ1点で接しているときの a であると考えることができます.それは次のように計算します.まず y=ax の微分 y′=axloga が1となる点を求めると, axloga=1x=−loga(loga) となります. y=ax と y=x はただ1点で接しており,傾きが 1 となる点 (x,y) は y=x のグラフの上にあるはずです.このような点 x において x=y=ax が成り立つので x=ax−loga(loga)=a−loga(loga)=(loga)−1=logaeloga=a−logae=e−1a=ee−1 というように a の値を求めることができます.つまり a>ee−1 となるような正の実数 a に対して aaa⋯ は ∞ に発散することになります.
ここで,この a を複素数を拡張したときに, aaa⋯ が発散する領域がどうなるのかを考えてみました.しかし計算によって何かしらの a の条件を求めることができなかったので,プログラムを組んで aaa⋯ の発散点をプロットしてみました.すると結構エグいカブトムシが現れました.
この図は zzz⋯ が発散するような点を黒くプロットしたものです.想像以上にやばい図が出来上がってしまいました.ちなみに zz を計算するときの log の偏角は −π から π までとしています.

この図の一部を拡大するとものすごいグルグルしています.プロット制度を上げれば「無限の回転」になると思われます.この図を持っておけばいつでも黄金の回転を放つことができます.以下に中心部のカブトムシはズームしていきます.
関数の反復により得られる図は面白いものを秘めていますね.有名なところではマンデルブロ集合があります.他の関数の反復についてもいろいろと気になったのでやってみました.以下の図はすべて範囲 −4<x<4,−4<y<4 の範囲で発散点を黒くプロットしたものです.
z(−z)z(−z)z(−z)⋯ が発散する点の集合 ↓
z¯zz¯zz¯z⋯ が発散する点の集合 ↓
z(Rez)z(Rez)z(Rez)⋯ が発散する点の集合 ↓
z(Imz)z(Imz)z(Imz)⋯ が発散する点の集合 ↓
⋯cos(cos(cos(cosz))))⋯ が発散する点の集合 ↓
⋯sin(sin(sin(sinz))))⋯ が発散する点の集合 ↓
⋯tan(cos(tan(cos(tan(cosz)))))⋯ が発散する点の集合 ↓
⋯cos(tan(cos(tan(cos(tanz))))))⋯ が発散する点の集合 ↓
⋯z+sin(z+sin(z+sin(z+sinz))))⋯ が発散する点の集合 ↓
⋯z+cos(z+cos(z+cos(z+cosz))))⋯ が発散する点の集合 ↓
このとき aaa⋯ は発散するかどうかのしきい値は, y=ax と y=x のグラフがただ1点で接しているときの a であると考えることができます.それは次のように計算します.まず y=ax の微分 y′=axloga が1となる点を求めると, axloga=1x=−loga(loga) となります. y=ax と y=x はただ1点で接しており,傾きが 1 となる点 (x,y) は y=x のグラフの上にあるはずです.このような点 x において x=y=ax が成り立つので x=ax−loga(loga)=a−loga(loga)=(loga)−1=logaeloga=a−logae=e−1a=ee−1 というように a の値を求めることができます.つまり a>ee−1 となるような正の実数 a に対して aaa⋯ は ∞ に発散することになります.
ここで,この a を複素数を拡張したときに, aaa⋯ が発散する領域がどうなるのかを考えてみました.しかし計算によって何かしらの a の条件を求めることができなかったので,プログラムを組んで aaa⋯ の発散点をプロットしてみました.すると結構エグいカブトムシが現れました.
真ん中あたりにカブトムシ( −3<x<3,−3<y<3 の範囲でプロット)
ここでいう「発散」は無限大への発散を表し,震動する点は除外している
ここでいう「発散」は無限大への発散を表し,震動する点は除外している
この図は zzz⋯ が発散するような点を黒くプロットしたものです.想像以上にやばい図が出来上がってしまいました.ちなみに zz を計算するときの log の偏角は −π から π までとしています.

この図の一部を拡大するとものすごいグルグルしています.プロット制度を上げれば「無限の回転」になると思われます.この図を持っておけばいつでも黄金の回転を放つことができます.以下に中心部のカブトムシはズームしていきます.

z(−z)z(−z)z(−z)⋯ が発散する点の集合 ↓
z¯zz¯zz¯z⋯ が発散する点の集合 ↓
z(Rez)z(Rez)z(Rez)⋯ が発散する点の集合 ↓
z(Imz)z(Imz)z(Imz)⋯ が発散する点の集合 ↓
⋯cos(cos(cos(cosz))))⋯ が発散する点の集合 ↓
⋯sin(sin(sin(sinz))))⋯ が発散する点の集合 ↓
⋯tan(cos(tan(cos(tan(cosz)))))⋯ が発散する点の集合 ↓
⋯cos(tan(cos(tan(cos(tanz))))))⋯ が発散する点の集合 ↓
⋯z+sin(z+sin(z+sin(z+sinz))))⋯ が発散する点の集合 ↓
⋯z+cos(z+cos(z+cos(z+cosz))))⋯ が発散する点の集合 ↓
「チルノの数学ノート」からここにたどり着きました。
返信削除他に類を見ない、おもしろい本ですね! 目からウロコでした。ぜひ続刊希望です。
さて、この「カブトムシ」、以前私も思いついて拙いプログラムを組んだでみたことがあります。
http://d.hatena.ne.jp/rikunora/20100504/p1
どうもこれは出来が良くなかったのですが、itaさんという方が美しく描き直してくれました。
http://d.hatena.ne.jp/ita/20100505/p1
それが、こんなところで出会えるとは、何とも嬉しくなりました。
なるほど、・・・sin(sin(sin z)・・・なども複雑な集合を形作るのですね。
http://d.hatena.ne.jp/rikunora/20170304/p1
ひょっとして、実数の世界でいわゆるカオスを生み出すような写像は、
複素数に拡大すると、皆このような、えも言われぬ図になるのでしょうか。
素人質問ですみませんが、冪乗の冪乗は結合法則が成り立たないため3の3乗の3乗を27の3乗と考えるか3の27乗と考えるかによって話が変わってくるかと思います。
返信削除漸化式的に書けば
1. I(n+1) = x^I(n)
2. I(n+1) = I(n)^x
とするかによって変わります。
本記事では自然と1を使っているように見えますが、前置きなしに1を使うのは四則演算の計算順序のようにどこかで決められているものなのでしょうか?
出典等あれば教えていただけると幸いです。