MathLabs

微分方程式と力学系

カオス理論

初期条件に極めて敏感な決定論的系であり、長期的には予測不可能な挙動を生み出す。

直観なぜ天気を二週間先まで予測できないのか?

天気予報は明日を予測しても来月を予測しても同じ物理法則——ニュートンの方程式、熱力学——を使うが、予報は約十日を過ぎると役に立たなくなる。その理由は方程式が間違っているからでも、コンピュータが遅すぎるからでもない。大気がカオス的だからである。ブラジルで一羽の蝶が羽ばたき、気温をわずかに変えるだけで、原理的には数週間後に大西洋でハリケーンが発生するかどうかを変えうる。これは神秘主義ではなく、初期条件への鋭敏な依存性と呼ばれる正確な数学的性質であり、1890年代に三体問題を研究していたアンリ・ポアンカレが独立に発見し、1963年に気象学者エドワード・ローレンツがコンピュータの印字結果を小数点以下六桁から三桁に丸めた際に再発見した。

カオス的アトラクタの蝶形の幾何学を示すインタラクティブな3Dパラメトリック曲面。
ロジスティック写像 xn+1=rxn(1−xn)x_{n+1} = r x_n(1 - x_n) のクモの巣図:パラメータ rr を 2.52.5(安定不動点)から 3.23.2(22 周期軌道)、そして 3.83.8(カオス)へと上げてみよう。

大学鋭敏な依存性とローレンツ系

定義: 初期条件への鋭敏な依存性

力学系が初期条件への鋭敏な依存性を持つとは、互いに距離 δ0\delta_0 だけ離れた近接点から出発した二つの軌道が指数関数的な速さ ∣δn∣≈∣δ0∣ eλn|\delta_n| \approx |\delta_0|\,e^{\lambda n} で分離していくことをいい、ここで λ\lambda はリアプノフ指数である。λ>0\lambda>0 のとき、ごくわずかな測定誤差 δ0\delta_0 でさえ有限の時間の後に予報を台無しにするほどの巨視的な大きさへと成長する——これこそ天気予報が約十日で役に立たなくなる理由である。

x˙=σ(y−x)y˙=x(ρ−z)−yz˙=xy−βz\begin{aligned} \dot x &= \sigma(y-x) \\ \dot y &= x(\rho - z) - y \\ \dot z &= xy - \beta z \end{aligned}

これはローレンツ系であり、大気対流を大胆に単純化したものである。xx は対流ロールの速度に、yy は上昇流と下降流の間の温度差に、zz は鉛直温度分布の線形からのずれに比例する。パラメータ σ\sigma はプラントル数、ρ\rho はレイリー数(下から流体がどれだけ強く熱せられているか)に比例し、β\beta は幾何学的な係数である。ローレンツの古典的なカオス的選択は σ=10\sigma=10、ρ=28\rho=28、β=8/3\beta=8/3 であり、これらの値ではあらゆる軌道が最終的に同じ蝶形のフラクタル集合(ローレンツ・アトラクタ)に落ち着くが、ある時刻にどちらの翅の上にいるかは事前にはほぼ予測不可能である。

xn+1=r xn(1−xn)x_{n+1} = r\,x_n(1-x_n)

はるかに単純な系でも同じ現象が見られる。ロジスティック写像 xn+1=r xn(1−xn)x_{n+1} = r\,x_n(1-x_n) は、資源に限りのある人口増加の一次元モデルである(ここで xn∈[0,1]x_n\in[0,1] は最大値に対する人口の割合、rr は増加率を制御する)。rr が小さいときはあらゆる軌道が単一の平衡点に収束する。rr が増加すると平衡点は安定性を失い、安定な2周期軌道に置き換わり、次いで4周期、8周期へと続く——周期倍増分岐のカスケードが r≈3.5699r\approx3.5699(ファイゲンバウム点)で蓄積し、それを超えると軌道は通常もはや繰り返さない。注目すべきことに、分岐点が集積していく速さ δ≈4.6692\delta\approx4.6692(ファイゲンバウム定数)は、無関係な膨大な種類の一次元写像に対して同一である——非線形力学において真に普遍的な定数のまれな例である。

rの増加に伴うロジスティック写像の長期的挙動
パラメータ長期的挙動
r<3r<3単一の安定な平衡点
r=3.2r=3.2安定な周期2軌道
r=3.5r=3.5安定な周期4軌道
r=3.57r=3.57カオスの始まり(周期倍増の集積)
r=4r=4完全にカオス的;リアプノフ指数 ln⁡2\ln 2 > 0

大学カオスの定量化:リアプノフ指数

λ=lim⁡n→∞1n∑i=0n−1ln⁡∣f′(xi)∣\lambda = \lim_{n\to\infty}\frac{1}{n}\sum_{i=0}^{n-1}\ln\left|f'(x_i)\right|

リアプノフ指数 λ\lambda は、典型的な軌道 x0,x1,x2,…x_0,x_1,x_2,\dots に沿った、各ステップでの局所的な引き伸ばし係数 ∣f′(xi)∣|f'(x_i)| の対数の平均であり、近接点が平均としてどれほど速く分離するかを測る。これは利率が平均的な指数関数的成長を測るのと同じである。λ>0\lambda>0 は近接点が指数関数的に発散すること——カオスの特徴——を意味し、λ<0\lambda<0 は安定な平衡点のように収束することを、λ=0\lambda=0 は境界的な(周期的または準周期的な)場合を意味する。

ロジスティック写像 f(x)=rx(1−x)f(x)=rx(1-x) について、非自明な平衡点 x∗=1−1rx^* = 1-\tfrac{1}{r} は 1<r<31<r<3 で安定であり、ちょうど r=3r=3 で安定性を失い、安定な周期2軌道へと周期倍増分岐する。

なぜ正しいのか?

平衡点における写像の導関数は、そこでの微小な摂動が一段階後にどれだけ成長または縮小するかを測る。その大きさが1を超えると、その点は引き込み型から反発型へと反転し、写像の2回合成写像の傾きがちょうどその閾値で (−1)2=1(-1)^2=1 となるため、まさに新たな周期2軌道が分岐しうる。

証明

ステップ1(平衡点を求める)。rx(1−x)=xrx(1-x)=x を解く、すなわち x[r(1−x)−1]=0x[r(1-x)-1]=0。これにより x=0x=0 または非自明な平衡点 x∗=1−1rx^* = 1-\tfrac{1}{r} が得られ、これは r>1r>1 のときちょうど (0,1)(0,1) に属する。

ステップ2(線形化)。微分すると f′(x)=r−2rxf'(x)=r-2rx。非自明な平衡点で評価すると f′(x∗)=r−2r(1−1r)=r−2r+2=2−rf'(x^*) = r - 2r\left(1-\tfrac1r\right) = r - 2r + 2 = 2-r となり、f′(x∗)=2−rf'(x^*) = 2-r。

ステップ3(安定性の判定基準)。平衡点が局所的に安定であるのはちょうど ∣f′(x∗)∣<1|f'(x^*)|<1 のときであり、すなわち ∣2−r∣<1|2-r|<1 で、これは 1<r<31<r<3 と書き換えられる。したがって x∗x^* はこの範囲全体で近くの軌道を引き寄せる。

ステップ4(r = 3 での分岐)。ちょうど r=3r=3 で f′(x∗)=−1f'(x^*)=-1 となり、線形化は境界的である。rr が3をわずかに超えると ∣f′(x∗)∣>1|f'(x^*)|>1 となり x∗x^* は反発点になる。一方、二回合成写像 f∘ff\circ f は連鎖律により、r=3r=3 ちょうどで x∗x^* における導関数が f′(x∗)2=1f'(x^*)^2=1 となる。f∘ff\circ f をより高次に展開すると、この退化点が rr が3を超えて増加するにつれて f∘ff\circ f の二つの真に新しい平衡点(すなわち ff の安定な周期2軌道)へと分裂することが示される——これが周期倍増分岐である。

「r=4r=4」において、ロジスティック写像 xn+1=r xn(1−xn)x_{n+1} = r\,x_n(1-x_n) は、ルベーグ測度に関してほとんど全ての初期条件 x0∈(0,1)x_0\in(0,1) に対してリアプノフ指数 ln⁡2\ln 2 を持つ。

なぜ正しいのか?

r = 4 での写像は一見単純な倍加写像とはまったく異なって見えるが、巧妙な変数変換(滑らかな共役)によって、あらゆる点で引き伸ばし率が明らかに2である倍加写像そのものに変換される——変数変換によって生じる余分な歪みは、長時間にわたって平均するとゼロに帰着する。

証明

ステップ1(共役)。x=sin⁡2(πy)x=\sin^2(\pi y) とおく。倍角の公式を用いると 4x(1−x)=4sin⁡2(πy)cos⁡2(πy)=sin⁡2(2πy)4x(1-x)=4\sin^2(\pi y)\cos^2(\pi y)=\sin^2(2\pi y)。倍加写像 T(y)=2y mod 1T(y)=2y \bmod 1 について y′=T(y)y' = T(y) とおくと、sin⁡2(πy′)=sin⁡2(2πy)\sin^2(\pi y') = \sin^2(2\pi y) は上の右辺とちょうど一致するので、h(y)=sin⁡2(πy)h(y)=\sin^2(\pi y) として f(h(y))=h(T(y))f(h(y)) = h(T(y)) が成り立つ:r = 4 におけるロジスティック写像は倍加写像と滑らかに共役である。

ステップ2(倍加写像のリアプノフ指数)。倍加写像は微分可能な箇所では至るところ T′(y)=2T'(y)=2 である区分線形写像なので、あらゆる軌道に沿って、任意の nn に対して 1n∑i=0n−1ln⁡∣T′(yi)∣=ln⁡2\tfrac1n\sum_{i=0}^{n-1}\ln|T'(y_i)| = \ln 2 がちょうど成り立つ——極限は自明に ln⁡2\ln 2 である。

ステップ3(共役を通して指数を運ぶ)。連鎖律により f(h(y))=h(T(y))f(h(y))=h(T(y)) を微分すると f′(h(y))h′(y)=h′(T(y))T′(y)f'(h(y))h'(y)=h'(T(y))T'(y)、すなわち x=h(y)x=h(y) において ln⁡∣f′(x)∣=ln⁡∣h′(T(y))∣+ln⁡∣T′(y)∣−ln⁡∣h′(y)∣\ln|f'(x)| = \ln|h'(T(y))| + \ln|T'(y)| - \ln|h'(y)| となる。この等式を軌道 y0,y1,…,yn−1y_0,y_1,\dots,y_{n-1} に沿って足し合わせ nn で割ると、中間の項がドミノ式に相殺し 1n∑i=0n−1ln⁡∣f′(xi)∣=ln⁡2+1n[ln⁡∣h′(yn)∣−ln⁡∣h′(y0)∣]\tfrac1n\sum_{i=0}^{n-1}\ln|f'(x_i)| = \ln2+\tfrac1n\big[\ln|h'(y_n)|-\ln|h'(y_0)|\big] となる。

ステップ4(境界項が消える)。h′(y)=πsin⁡(2πy)h'(y)=\pi\sin(2\pi y) は有界(その大きさは決して π\pi を超えない)なので、括弧内の項 ln⁡∣h′(yn)∣−ln⁡∣h′(y0)∣\ln|h'(y_n)|-\ln|h'(y_0)| は、h′h' の可算個の零点——ルベーグ測度ゼロの集合——を避けるあらゆる軌道について n→∞n\to\infty としても有界のままである。有界な量を nn で割ると 00 に収束するので、r = 4 におけるロジスティック写像のリアプノフ指数は、ほとんど全ての初期条件について ln⁡2\ln 2 に等しい。

大学実世界での応用と具体例

カオス理論は天気だけの話ではない。心拍モニターがカオス的な細動の発生を監視する理由、安全なカオス通信システムがメッセージをカオス的な搬送信号に混ぜ込み、まったく同じパラメータを持つ受信機だけが復号できる理由、繁殖率 rr の高い個体群生態モデルが外部からのランダム性なしに激しく振動しうる理由、そして技術者が乱数生成や暗号化のためにわざとカオス信号を生み出す電子回路(チュア回路など)を設計する理由を説明する。いずれの場合でも同じ二つの計算が繰り返される。平衡点や周期軌道を求めて線形化することと、小さな誤差がどれほど速く爆発的に拡大するかを示すリアプノフ指数を見積もることである。

例: 永遠に振動する個体群:r = 3.2 における周期2軌道

ある漁業では、来年の資源量の割合をロジスティック写像 xn+1=r xn(1−xn)x_{n+1} = r\,x_n(1-x_n)(成長パラメータ r=3.2r=3.2、周期倍増の閾値 r=3r=3 を超える)でモデル化している。過渡的な振る舞いが消えた後、資源量が交互に取る二つの個体数水準を求めよ。

解答

周期2軌道の点は、それ自体は平衡点ではないが f(f(x))=xf(f(x))=x を満たすので、f(f(x))−xf(f(x))-x はゼロにならなければならない。f(f(x))−xf(f(x))-x を展開し、すでに f(x)=xf(x)=x を満たす通常の平衡点 x=0x=0 と x∗=1−1/rx^*=1-1/r(したがって自明に f(f(x))=xf(f(x))=x)である二つの根を除くと、真に新しい因子 r2x2−r(r+1)x+(r+1)=0r^2x^2-r(r+1)x+(r+1)=0 が残る。

これは xx に関する二次方程式である。二次方程式の解の公式により x=(r+1)±(r−3)(r+1)2rx=\dfrac{(r+1)\pm\sqrt{(r-3)(r+1)}}{2r}。

r=3.2r=3.2 を代入すると、分子の項は r+1=4.2r+1=4.2 および (r−3)(r+1)=0.2×4.2=0.84(r-3)(r+1)=0.2\times4.2=0.84 となり、0.84≈0.9165\sqrt{0.84}\approx0.9165、2r=6.42r=6.4 である。これにより x≈{0.513, 0.799}x\approx\{0.513,\ 0.799\} が得られる。

したがって漁業資源は、収容力の約 51.3%51.3\% と 80.0%80.0\% の間を毎年交互に取る振動へと落ち着く——これは系がカオス的になるずっと前に、完全に決定論的な規則から生じる、真に周期的な個体群サイクルである。

例: なぜ十日なのか?天気の予報可能限界を見積もる

大気モデルのリアプノフ指数は一日あたり約 λ≈0.9\lambda\approx0.9 と見積もられている。二つの予報計算が(正規化された単位で)δ0=10−6\delta_0=10^{-6} だけ離れた状態で始まるとし、これはほぼ完璧な初期測定を表す。n=10n=10 日後にこのごくわずかな差がどれほど大きくなっているかを見積もり、それが予報にとって何を意味するか述べよ。

解答

リアプノフ指数の定義により、小さな誤差は(平均として)δn≈δ0 eλn\delta_n \approx \delta_0\,e^{\lambda n} のように成長する。

λ≈0.9\lambda\approx0.9(1日あたり)と n=10n=10 日から、指数は λn≈0.9×10=9\lambda n \approx 0.9\times10=9 となり、成長率は e9≈8103e^{9}\approx8103 である。

δ0=10−6\delta_0=10^{-6} から出発すると、十日後には誤差は δ10≈10−6×8103≈8.1×10−3\delta_{10}\approx10^{-6}\times8103\approx8.1\times10^{-3} にまで成長する——当初はごくわずかで実質的に測定不可能だった差が、約八千倍に増幅されたことになる。

実際の測定・モデル化の誤差は相対的にすでに 10−610^{-6} よりはるかに大きいため、この指数関数的な急増こそが、物理モデルがどれほど優れていても、コンピュータがどれほど強力であっても、業務用の天気予報が約一〜二週間で予測能力を完全に失う理由である。それは λ\lambda によって定まる数学的な限界であり、工学的な制約ではない。

研究研究最前線:厳密なカオスと高次元乱流

非自明な平衡点が x∗=1−1rx^* = 1-\tfrac{1}{r} であるロジスティック写像 f(x)=rx(1−x)f(x)=rx(1-x) について、導関数 f′(x∗)f'(x^*) を rr で表すとどうなるか?

二つの天気予報の実行が δ0=10−6\delta_0=10^{-6}(正規化単位)離れて始まる。δn≈δ0eλn\delta_n\approx\delta_0 e^{\lambda n} と λ≈0.9\lambda\approx0.9(1日あたり)を用いると、n=10n=10 日後に隔たりはおよそ何倍に成長しているか?

ロジスティック写像の非自明な平衡点 x∗=1−1rx^* = 1-\tfrac{1}{r} が周期倍増によって最初に安定性を失うのはパラメータ rr がどの値のときか?

ある力学系のリアプノフ指数が λ>0\lambda>0 である。これは何を示すか?

参考文献

  1. Steven H. Strogatz (2015). Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering
  2. Edward N. Lorenz (1963). Deterministic Nonperiodic Flow
  3. Robert M. May (1976). Simple mathematical models with very complicated dynamics
  4. Warwick Tucker (2002). A Rigorous ODE Solver and Smale's 14th Problem