PCが描く奇妙な画像集(数学的万華鏡と生物形態等の世界)

・インタープリタBASICによるフラクタルとカオスの奇妙な画集。

226 ヴェルンハルスト方程式画像(その3)

2014-07-20 13:36:32 | ヴェルンハルスト方程式
ヴェルンハルスト方程式は下記で示される漸化式であった。(記事274参照)
  X(t+1)=A{X(t)(1-X(t)}

即ち、初期値を t=0 の時のX(0)とすれば、上式より次の数列が得られる。

X(0),X(1),X(2),X(3),・・・・・,X(t),・・・・・

この数列がある値に収束するのか、発散するのか、振動するのか、あるいはカオス状態になるのかは記事m275で見てきた。 その画像を再び下に示す。



上図は t が充分大きいときの漸化値であるが、
X(0),X(1),X(2),X(3),・・・・・,X(t),・・・・・
全てを表示したらどうなるだろうか。それは、A が与えられた時の漸化の様子を示す画像となる。

下図は、X(0),X(1),X(2),X(3),・・・・・,X(t),・・・・・,X(500)-----------------(1)
ときの図である。



上図で、t の値は色で区別している。t=0,1,2,3,・・・,500 としている。
BASIC/98 では下図に示すように16色しか使えず、色にはパレット番号が付けられている。
そこで、t の区別は、パレットNo.C=t MOD 16 として区別した。



従って、色から、上記数列のアトラクタへの接近の様子が分かる。

この漸化の様子も、図の横軸:Aに対して、ある種の規則性がありそうだ。
例えば、図中のP点のように、上記(1)の全ての数列の値が一点で交わるようなA点が存在している。これは別の意味でのアトラクタと言えるのではないか? 詳細に調べると、この図には、いろいろな規則性があるようだ。

下図は、t が、ある値以下のX(t)は表示させない場合の図である。
上の図から、t を 50,100,400 以下のX(t)を表示させない図である。
順次、数列{X(t)}の過渡値が削除されめため漸近値(アトラクタ)が明確になっていく。但し、カオス状態は変わらない。








--------------------------------------------------------
以下は、A=3.5~4の場合の画像である。
下図は記事274で示したアトラクタ及びカオスの画像である。



下図は、t=0,1,2,3,・・・,500 の全てのX(t)を表示した図である。
上のA=1~4図同様に、上記(1)の( ほとんど )全ての数列の値が或る一点 P に交わるような A 点が存在する。



下図は、t が、ある値以下のX(t)は表示させない場合の図である。数列{X(t)}の過渡値が削除されめため漸近値(アトラクタ)が明確になっていく。










275 ヴエルンハルスト方程式画像(その2)

2014-07-20 11:07:19 | ヴェルンハルスト方程式
記事274にヴエルンハルスト方程式:X(t+1)=AX(t){1-X(t)} ・・・・・(1)
(ここで、A は定数、X(t)は変数)
画像を掲載した。
これはtの漸化式であるが、X(t)の漸化値が定数Aによって様相が一変することを見た。
A の増加に従い、X(t)の漸化値は、ある固定値から始まり分岐する。
その分岐は更に分岐し続け、あるA値からX(t)の値はカオス状態に突入したのであった。
下図は記事274に掲載した、その様子を示す画像である。



では、ヴエルンハルスト方程式(1)式において、Aを固定して、横軸にtをとり、
縦軸にX(t)にとって、t の変化に対してX(t)はどのように変化するのだろうか?
それを調べたのが下図である。A の代表的に四つの値で見てみよう。
A=2, A=3.3, A=3.5とA=4の場合である。
各場合の画像を以下に示す。









上図から分かるように一番上の画像と、A=2, A=3.3, A=3.5とA=4を固定した t 対 X(t)画像は
一致している。当然のことだが、A=4の場合のX(t)のカオス画像の様子は面白い。
なるほど、X(t)はtに対して予想不可能に見える。

---------------------------------------
追記:


274 ヴェルンハルスト方程式画像(その1)カオス画像の神秘さ

2014-07-20 10:35:38 | ヴェルンハルスト方程式
カオス関連の一般読者向けの本には恐らく必ず出てくるものに、ヴェルンハルスト方程式というものがある。これは中学生でも分かる以下の式だ。

X(t+1)=AX(t){1-X(t)} (ここで、A は定数、X(t)は変数)

これは漸化式で、X(t)が与えられた時、上式の右辺よりX(t+1)が求められるが、今度そのX(t+1)を上式の右辺のX(t)に置き換えて新しいX(t+1)を求める。
この操作を繰り返すことによって下式の数列が得られる。

X(0),X(1),X(2),・・・・,X(t),X(t+1),・・・・

この数列は際限なく単純に増加していったり、あるいは際限なく減少していったり、あるいは上がったり下がったりして振動して、ともかく、ある特定の値に漸近したり、振動したりする・・・と、フツーの人(勿論、私も含めて)思うだろう。

こんな単純な式には何の変哲もない退屈なコトしか無いだろう、と思うのが人情というものだ。

ところが、だ。知っている人には釈迦に説法だが、この単純極まりない式には、トンデモナイことが実は潜んでいるのだ。

確かに、この式は何の変哲もないコトから始まるのだが、ある時、突然、奇っ怪なことが突然始まるのだ。その奇っ怪なこととは、今や日常語にもなった『カオス』のことだが、つまり、『次に何が起こるか予想が付かない』が発生するのだ。

この現象を見事に目の当りに見せてくれる画像がある。
実は、この画像は今やよく知られた画像だが、私は自分で(プログラムを作ってみて)体験してみた。PC の画面に徐々に現れる、その画像を見ていると、その奇っ怪さが身にしみてくる。( 何事も自身で体験してみなければワカラナイことは世の中には多い。
この体験もその一つだろう。)

その奇っ怪な画像とは下に添付した画像だが、この画像の縦軸は数列{X(t)}の漸化値で(但し1に規格化してある)、横軸は定数 A だ。

A がある値までは何事もなく漸化値(アトラクタと呼ばれる)は単調に増加していくのだが、A=3あたりになると、アトラクタは二つに分岐する。つまり、数列{X(t)}の漸化値は、二つの値を行ったり来たりする。そしてAが増加すると更に分岐が始まる。更に又この分岐が急速に進み、ついにカオス状態へ突入していく。

カオス状態では『数列{X(t)}の漸化値』は予測不能な状態となる。

しかし不思議なことに、このカオス状態の画像をよく見ると、ある種の秩序がある
ことが分かる。カオスの画像の濃淡の中にある複雑な放射線状のモノがそれだ。
そして更によく見ると、このカオス画像はフラクタル性(自己同一性)もある。

こんな簡単な式(ヴェルンハルスト方程式)に、一体、何ゆえに、このような複雑さ・豊穣さが
潜んでいるのだろう!!!??? 不思議とは、こういうことだ。

実は、この画像のカオス性は数学の単なるお遊びではなくて、この世の森羅万象のには、よく見られる現象なのだ。

ここらの話の本はたくさんあるだろうが、私のお気に入りの本を紹介しておこう。
これも古い本だ。『鏡の伝説--カオス-フラクタル理論が自然を見る目を変えた』(J.ブリッグス、F.D.ピート著、ダイヤモンド社)。










273 Z^2マンデル集合の中の部分の拡大画像の画像処理(その2)

2014-07-20 10:11:49 | Z^2マンデルブロー集合
上図の1-14画像の中の画像(1-14-1~1-14-5画像)を拡大して,
それらの画像の、LOG(No)画像及び赤黒縞画像を見よう。
先ず、1-14-1~1-14-5画像の場所を以下に示す。







------------------------------------

以下に各拡大画像を示すが、比較しやすいように、各拡大画像において、

・オリジナル画像(C=No MOD 16)と、
・LOG(No)画像
・赤黒縞画像( Noが偶数の時→赤、奇数の時→黒で、かつ、No>Naの場合のみ表示 )

とを並列して示す。

--------------------------------------

・1-14-1 画像







--------------------------------------

・1-14-2 画像







-------------------------------------

・1-14-3 画像







------------------------------------

・1-14-4 画像







-----------------------------------

・1-14-5 画像










272 Z^2マンデル集合の中の部分の拡大画像の画像処理(その1)

2014-07-20 08:52:07 | Z^2マンデルブロー集合
下図の1-13画像の中の画像を拡大して,それらの画像の、LOG(No)画像及び赤黒縞画像を見よう。先ず、1-13画像の場所を以下に示す。





------------------------------------





以下に各拡大画像を示すが、比較しやすいように、各拡大画像において、

・オリジナル画像(C=No MOD 16)
・LOG(No)画像
・赤黒縞画像( Noが偶数の時→赤、奇数の時→黒で、かつ、No>Naの場合のみ表示 )

とを並列して示す。

--------------------------------------

・1-13-1 画像







--------------------------------------

・1-13-2 画像







-------------------------------------

・1-13-3 画像






------------------------------------

・1-13-4 画像







-----------------------------------

・1-13-5 画像








上図において、特に興味深いのは( 他の記事の画像もそうだが )、LOG(No)画像である。

この図の色は、Noの分布( 但し、LOG(No )での分布)を表しているのだが、その分布は決して一様ではなく、いろいろな形状で分布している。この形状には何か構造があるのだろうか?

全くデタラメには見えないし、かといって、種々の似たような形状が随所にある(それは、この自然界で見かける形状にも似ている!!。)また、ミニ・マンデル集合も随所に散乱して見える。

これらには、何か構造があるのだろうか・・・?


271 f(Z)マンデルブロ集合画像:f(z)=sinZ,cosZ,tanZ,e^Z,sinhZ,sinsinZ

2014-07-20 07:52:44 | その他の複素関数のマンデルブロ集合
f(z)=sinZ, cosZ, tanZ, e^Z, sinhZ ,sinsinZ のマンデルブロ集合画像はどのような画像になるのだろうか?

下図はその画像である。図の黄色の部分がマンデルブロ集合部分である。

画像の上段左から、sinZ, cosZ, tanZ
画像の下段左から、e^Z, sinhZ, sinsinZ である。

CXS, CXE は下記のとおり。
sinZの場合:CXS=0 CXE=2π
cosZの場合:CXS=0 CXE=2π
tanZの場合:CXS=0 CXE=π
e^Zの場合:CXS=-π/2 CXE=π
sinhZの場合:CXS=0 CXE=π
sinsinZの場合:CXS=-π CXE=π

なお、N-loop脱出条件は、X^2+Y^2>2500である。




-----------------------------------------------------------
sinsin Z マンデルブロ画像を以下に示す。

画像作成条件は以下のとおり。
・複素関数:sinsin Z
・CX,CY範囲:|CX|<1.5π,|CY|<1.1π
・N-loop脱出条件:Q=X^2+Y^2>2500
if Q>2500 then pset (K,J,),C
  C=No mod 16 但し、C=6→5,C=7→8
・N-loop貫通時(収束時)はC=6(黄)
・Nmax=500


270 Z^nマンデルブロ集合画像(n=3,4,5,6,7)

2014-07-20 07:43:05 | その他の複素関数のマンデルブロ集合
Z^nマンデルブロ集合(n=3,4,5,6,7)は、どのような画像になるのだろうか?

下図はその画像である。黄色の部分がマンデルブロ集合部分である。

画像の上段左から、Z^2, Z^3, Z^4
画像の下段左から、Z^5, Z^6, Z^7 である。

CXS, CXE は下記のとおり。
Z^2の場合:CXS=-1.5 CXE=0.5
Z^3の場合:CXS=-1.2 CXE=1.2
Z^4の場合:CXS=-1.4 CXE=1.2
Z^5の場合:CXS=-1.2 CXE=1.2
Z^6の場合:CXS=-1.2 CXE=1.2
Z^7の場合:CXS=-1.2 CXE=1.2

なお、N-loop脱出条件は、x^2+Y^2>4である。


269 (Z^2+sinZ)マンデルブロ画像

2014-07-20 07:26:22 | その他の複素関数のマンデルブロ集合


上図は、(Z^2+sinZ)マンデルブロ画像の全体図である。独特な形をしていて面白い。
上図の5個所の部分を下図に示すように選び拡大してみる。





下図は各拡大部分の2以降の画像である。

各拡大部分の図は上から順に下記の画像作成条件の画像である。

(1) 色C=No MOD 16 の画像。但し、収束部分は黄色にしている。
(2) 色C=INT(LOG(10^6*No) MOD 16 の画像。
  このLOG画像については、アルバム500の11を参照。
(3) Noが偶数の時→赤、奇数の時→黒とした画像。かつ、NoがNa以下の場合のみ表示。
(4) NoがNa以上の場合のみ表示した画像。

結論として、これらの画像より、No の構造(分布)の概要が分かるが特別に特記すべき事柄は見出せない。 
また、ミニ・マンデルブロ画像も見つからない。

---------------------------------------

・1 の拡大部分の画像



この部分に何か面白いモノが潜んでいると思って拡大してみたが、格別なモノはない。

---------------------------------------

・2 の拡大部分の画像









---------------------------------------

・3 の拡大部分の画像









二番目のLOG(No)画像で、画像の途中で黄色→白色となっているが、これはプログラム・ミスで特別な理由はない。

---------------------------------------

・4 の拡大部分の画像









--------------------------------------

・5 の拡大部分の画像