(問題)

ある鉛直面(t,x)内において、定点Oから定点A(a,b)まで、質量Mの質点が、ある経路に沿って重力のみにより運動するとき、質点が点Oから点Aまで最も早く到達するときの経路x(t)はどのような曲線になるだろうか?


P1001249.jpg



質点がx落下した時の速度は√2gxであるから、質点がOから位置xを通過するまでの時間をτ、点Oから位置xまでの経路長をsとすると、

ds/dτ=√(2gx)

よって

dτ=ds/2gx=(dx^2+dt^2)/2gx=(1+x’^2)/2gxdt

よって、τは以下のtについての0からaまでの定積分

τ=∫√(1+x’^2)/2gxdt

となる。

よって求めるx(t)の関数形は

評価関数(汎関数)I(x(t))=∫√(1+x’^2)/xdt


を最小にする関数x(t)を求める変分問題に帰着する。

ここで、Iの右辺の被積分関数をf(t,x,x’)=(1+x’^2)/xとおく。

オイラーの変分の必要条件は

(1) f-∂f/x’x’=(f/t)・dt

ここで、(以下計算)

f/t=0

f/x’=x’/(x(1+x’^2))

これを(1)式に代入して整理すると

x(1+x’^2)=C (Cは定数)

上式において、x’=cotθとおいて

x=C/(1+cotθ^2)=Csinθ^2=C(1cos2θ)/2

また、

dt/dθ=(dx/dθ)(dt/dx)=Csin2θ/cotθ=2Csinθ^2=C(1cos2θ)

積分して

t=C(θ-(sin2θ)/2)D

ここで初期条件t=0の時x=0を考慮してD=0

以上より、求める関数はθをパラメータとして

t=C(2θ-sin2θ)

x=C(1cos2θ)


となる。これはサイクロイド曲線である。

以下の比較的簡単な人口モデルの漸化式を考えます。(Ricker方程式といいます)



  x(k+1)=α(k)x(k)exp(-x(k))



kは第k

x(k)k年の人口密度(変数)、

α(k)は人口増加率(パラメータ)

exp(-x(k))は高人口密度での競争による人口減少を考慮した項


ここで、k、αがどんな値であろうと、漸化式①は自明な不動点x(k)=0を持ちます。

また、①はαに依存する不動点として、x(k,α)=lnαを持ちます。

ここで、不動点がx(k,α)=lnαのときの①のαに依存する固有値をλとすると、 λ(α)=1-lnαとなります。

よって、x(1)1<α<α(1)で安定。α(1)<αのとき不安定となります。

すなわち、システムはα=α(1)=e^2のときに分岐します。

同様にkの値を増やしていくと、以下のようなx(k)とα(k)の関係図(分岐ダイヤグラム)が得られます。



P1001243.jpg

上図でx(1)は上でのべたように、α(1)を境に安定性が失われ、x(2)なるリミットサイクルを生じます。分岐点α(2)でも同様の分岐が二つ生じることになります。分岐点α(3)、α(4)・・・でも同様の分岐が逐次生じます。

ここで、次にパラメータα(k)の比β


②  β={α(k)-α(k-1)}/{α(k+1)-α(k)}

を考えます。この比βはkが増加するに従ってある一定の値に漸近します。この比βの極限値をFeigenbaum(ファイゲンバウム)定数といい、

β=4.669201....となります。

ここで、上記のようなflip分岐のカスケードを示す全てのシステムにおける比βは普遍値となりFeigenbaum定数をとります。この普遍値は、円周率やオイラー数、ひいては微細構造定数などに負けず劣らず、非常に神秘的といえると思います。



























最果ての荒野を歩いていた
ひとつのキャラバンが通りすがった

美しい瞳の少女がいた
僕の存在に気づいたのだろうか
少し戸惑いながら僕は
近づき少女に挨拶した

黒い瞳に僕は吸い込まれた

僕は迷うことなく彼女の生涯を唄った

少女は不思議そうに微笑んだ
僕は唄い続けた
言葉が溢れて止どまらなかった

隊商は流れに任せて進んでいく
群衆のなか僕は彼女を見失った

それでも僕は彼女のために唄い続けた

夕闇に捕らわれ
やがて声も枯れた
それでも僕は止めることができなかった




最後に僕はひとつの祈りを少女に捧げて
その場を通りすぎた