力学系とカオス　決定論なのに予測できない世界｜武蔵野個別指導塾

同じ式を繰り返すだけで、収束と不規則な変動が分かれます。ロジスティック写像の計算から、安定性、周期倍分岐、初期値の誤差が広がる仕組みをたどり、予測の限界を数学で捉えます。

力学系とカオス　決定論なのに予測できない世界

停電した部屋で考える、明日の空の決まり方

台風の夜、停電した部屋で、机の上のノートだけを小さな明かりが照らしているとします。窓が鳴るたびに、風は次にどちらから吹くのだろうと気になります。数学の問題なら、条件をそろえて式に代入すれば答えが出るはずです。それなのに、空気の動きに法則があることと、遠い先の天気を細かく言い当てられることは、同じではありません。この違いが、力学系を学ぶ出発点です。

力学系は、状態が時間とともに変わる仕組みを研究する数学です。振り子の位置と速度、ある時点の個体数、繰り返し計算の途中の値などを「状態」としてまとめ、その状態を次へ進める規則を考えます。力学という名前でも、力や物体だけを扱うわけではありません。高校で習う数列も、各項を時刻ごとの状態と読めば、すでに力学系の入口にあります。

決定論的であるとは、同じ初期状態から同じ規則を使えば、その後の状態が一意に定まることです。しかし実際の観測値には幅があり、計算にも丸めが入ります。出発時の小さな違いが時間とともに拡大するなら、法則が正しくても長期の細かな予測は難しくなります。カオスを理解するには、「未来が定まらない」と「未来を必要な精度で計算できない」を分ける必要があります。

ここでは、短い漸化式を実際に計算するところから始めます。落ち着く値、繰り返す値、誤差の増え方を順に調べると、不規則に見える動きにも数学で問えることが残っているとわかります。大学数学が加えるのは、計算量だけでなく、どの条件ならどこまで言えるかを確かめる視点です。

漸化式の一般項から、繰り返した先の姿へ

高校の数列で、x₀ = 3、xₙ₊₁ = xₙ/2 を考えます。順に代入すると、3、3/2、3/4、3/8、3/16 と続きます。等比数列なので一般項は xₙ = 3/2ⁿ です。したがって n → ∞ のとき xₙ → 0 となります。例えば第10回の更新後は x₁₀ = 3/1024 です。この例では、一般項を求める方法と、毎回半分になるから0に近づくという説明が、同じ結論へつながります。

別の初期値 y₀ から同じ規則で進めれば、2つの数列の差は |xₙ-yₙ| = |x₀-y₀|/2ⁿ です。初期値の違いも半分ずつになります。ここには「値が0へ近づく」という性質と、「近い初期値から始めた軌道が近いままでいる」という性質が含まれています。力学系では、この2つを区別して定義します。一般項が見つからない場合にも、後者を調べることで未来の姿を判断できるからです。

次に、次の漸化式へ移りましょう。

xₙ₊₁ = r xₙ(1-xₙ)

これをロジスティック写像と呼びます。r は途中で変えないパラメータで、ここでは 0 ≦ r ≦ 4、0 ≦ x₀ ≦ 1 とします。積 r xₙ は r と xₙ の掛け算を表しています。x(1-x) = 1/4-(x-1/2)² なので、区間 [0,1] では 0 ≦ x(1-x) ≦ 1/4 です。従って 0 ≦ r x(1-x) ≦ 1 となり、更新後も同じ区間にとどまります。この確認なしに反復を始めると、数式が想定する状態の範囲を見失います。

r = 2、x₀ = 0.2 とすれば、x₁ = 2 × 0.2 × 0.8 = 0.32 です。さらに進めると x₂ = 0.4352、x₃ = 0.49160192、x₄ ≈ 0.4998589445 となり、0.5へ近づいているように見えます。ただし、この数個の値だけで「必ず収束する」と証明したことにはなりません。有限回の観察から、無限回の結論へ進むには別の論理が要ります。

数学Iの二次関数は状態の範囲を調べる道具になり、数学Bの数列は時間を表す道具になります。数学II・IIIの微分は、近い入力の違いが次の出力で何倍になるかを調べる道具です。これらの基礎は多くの理系学部では1〜2年で学びますが、力学系を体系的に扱う時期は課程によって異なります。高校から大学への変化は、解く式が難しくなることだけではありません。「一般項は何か」から、「一般項なしでも全体の動きを分類できるか」へと問いが広がるのです。

状態・軌道・安定性を言葉と式で区別する

まず離散力学系を定義します。状態の集合 X と、各状態に次の状態を対応させる写像 f: X → X を与え、非負整数 n に対して xₙ₊₁ = f(xₙ) と定める仕組みです。初期状態 x₀ から得られる列 x₀, x₁, x₂, … を前向き軌道と呼びます。「離散」は、時刻を整数の刻みで数えるという意味です。状態そのものが整数である必要はなく、先ほどの X = [0,1] のように、実数の連続した範囲を取ることもあります。

同じ写像を n 回合成したものを fⁿ と書けば xₙ = fⁿ(x₀) です。この fⁿ は関数値の n 乗ではありません。例えば f²(x) = f(f(x)) です。f⁰ は何も変えない恒等写像とします。力学系で計算が崩れる原因の一つは、この合成と通常の累乗を取り違えることです。まず括弧を残して2回の更新を書けば、区別がつきます。

連続時間の場合は、時刻 t を実数として状態 x(t) を考えます。微分方程式 dx/dt = F(x) は、その代表的な与え方です。解がすべての実数時刻で一意に存在する設定では、時刻 t だけ進める写像 φₜ は φ₀(x) = x、φₜ₊ₛ(x) = φₜ(φₛ(x)) を満たします。これを流れと呼びます。先に s 時間、次に t 時間進むことが、合わせて s+t 時間進むことと同じになる、という条件です。微分方程式なら常に全時刻の流れができるわけではなく、解の存在範囲も調べます。

例として dx/dt = -x の解は x(t) = x(0)e⁻ᵗ です。微分すると右辺の -x(t) に一致し、t = 0 では与えた初期値に戻ります。x(0) = 3 なら x(t) = 3e⁻ᵗ で、連続的に0へ近づきます。離散の等比数列と似ていますが、途中の時刻にも状態があることが違いです。離散系では写像の反復、連続系では微分方程式の時間発展という形で、同じ「軌道」の発想を使えます。

次に不動点を定義します。離散系の点 p が f(p) = p を満たすとき、p を不動点と呼びます。一度そこへ来ると、その後の状態は変わりません。連続系では F(p) = 0 を満たす平衡点が対応します。周期点は、ある正整数 k について fᵏ(p) = p となる点です。そのような k の最小値を最小周期と呼びます。不動点も周期1の点なので、「周期2の点」と言うときには f²(p) = p に加えて f(p) ≠ p を確かめます。

安定性は、まず実数上の離散系で定義しましょう。不動点 p が安定であるとは、任意の ε > 0 に対して、ある δ > 0 が存在し、X 内の初期値 x₀ が |x₀-p| < δ を満たせば、すべての整数 n ≧ 0 について |fⁿ(x₀)-p| < ε となることです。どの程度の近さを要求されても、出発時の近さを十分に選べば、未来の全時刻でその要求を守れるという意味です。δ は ε に応じて選べますが、n ごとに選び直すことはできません。

さらに、ある η > 0 が存在し、|x₀-p| < η を満たす X 内のすべての初期値で fⁿ(x₀) → p となるなら、p は局所的に引きつけます。安定であり、かつこの引きつける性質を持つとき、局所漸近安定と呼びます。恒等写像 f(x) = x の不動点は安定ですが、別の点から出発してもそこへ近づかないので、実数上では局所漸近安定ではありません。「近くに居続ける」と「最後には近づく」の違いが、この例で見えます。

初期値鋭敏性にも量化の順序があります。距離 d を持つ状態空間 X 上で、ある a > 0 が存在し、任意の x ∈ X と任意の δ > 0 に対して、d(x,y) < δ となる y ∈ X と非負整数 n があり、d(fⁿ(x),fⁿ(y)) > a となるとき、初期値に鋭敏に依存するといいます。どれほど小さな近傍にも、将来一定以上離れる出発点があるという条件です。近い2点をどう選んでも必ず離れる、という意味ではありません。実際、ロジスティック写像では x と 1-x の次の値は同じです。

カオスには複数の数学的定式化があります。初期値鋭敏性は重要な性質ですが、それだけをあらゆる文脈でカオスと同一視しないようにします。状態空間をどこに取るか、軌道がどの領域にとどまるか、周期軌道がどう配置されるかも問われます。名称を覚える前に、何を仮定してどの性質を確かめたのかを書けることが大切です。

微分で判定する不動点と、2周期が生まれる計算

1次元の不動点には、微分を使う基本定理があります。f が不動点 p を含む開区間上で連続微分可能であるとします。|f′(p)| < 1 なら p は局所漸近安定であり、|f′(p)| > 1 なら p は不安定です。|f′(p)| = 1 の場合は、この判定だけでは結論が出ません。これらは p の近くでの主張であって、遠くにあるすべての初期値の行き先を決めるものではありません。

安定になる理由を見ましょう。|f′(p)| < q < 1 となる q を選ぶと、導関数の連続性により、p の十分小さな近傍では |f′(x)| ≦ q とできます。平均値の定理を使うと、その近傍の x に対して次が成り立ちます。

|f(x)-p| = |f(x)-f(p)| ≦ q|x-p|

誤差が q 倍以下になるので、その近傍から出ることはありません。繰り返せば |fⁿ(x)-p| ≦ qⁿ|x-p| です。0 < q < 1 だから右辺は0へ向かいます。この「近傍を保ったまま同じ評価を繰り返せる」という一歩が、数回の数値観察を収束の証明へ変えています。

ロジスティック写像 f(x) = r x(1-x) で、r > 0 とします。不動点の式 p = rp(1-p) を整理すると p(r-1-rp) = 0 なので、p = 0 または p = 1-1/r です。後者が [0,1] の内部にあるのは r > 1 のときです。導関数 f′(x) = r(1-2x) を代入すると、0での値は r、もう一つの不動点での値は 2-r です。従って 0 < r < 1 では0が局所漸近安定で、1 < r < 3 では p = 1-1/r が局所漸近安定です。

r = 2 では p = 1/2 となり、f′(p) = 0 です。p からのずれを eₙ = xₙ-1/2 と書くと、近似ではなく厳密に次の式が得られます。

eₙ₊₁ = -2eₙ²

例えば区間 [0.4,0.6] では |f′(x)| ≦ 0.4 です。この区間に入れば、先ほどの収縮の議論を使えます。x₀ = 0.2 の軌道は x₂ = 0.4352 でその区間に入るため、以後は1/2へ収束します。最初の何回かは計算で追い、その後の無限回は定理に任せるという、数値と証明の分担ができます。

では r = 3 を越えるとどうなるでしょうか。r = 3.2 の不動点は p = 0.6875 で、f′(p) = -1.2 です。絶対値が1より大きいので不安定になります。負号は、不動点の左右を入れ替える向きにずれが伝わることを表します。値が増大するだけではなく、左右に振れながら不動点から離れるため、次に周期2の軌道を調べたくなります。

2回更新して戻る条件は f²(x) = x です。展開して因数分解すると、次になります。

f²(x)-x = -x(rx-r+1){r²x²-r(r+1)x+r+1}

前の2つの因子は、すでに求めた不動点に対応します。新しい候補は中括弧を0とした二次方程式から得られ、r > 3 では異なる2点 a、b になります。

a = {r+1-√((r-3)(r+1))}/(2r)

b = {r+1+√((r-3)(r+1))}/(2r)

r = 3.2 なら a ≈ 0.5130445095、b ≈ 0.7994554905 です。実際に f(a) = b、f(b) = a となります。周期軌道の安定性を見るには、f を1回適用したときの傾きだけでなく、f² の不動点としての傾きを調べます。連鎖律により (f²)′(a) = f′(b)f′(a) であり、上の2点ではこれが -r²+2r+4 に等しくなります。

r = 3.2 を代入すると -10.24+6.4+4 = 0.16 で、絶対値は1未満です。従ってこの2周期軌道は局所的に引きつけます。より一般に 3 < r < 1+√6 で同じ判定が使えます。上端は約3.449489743です。不動点が不安定になっても、軌道がすぐ無秩序になるとは限りません。1点へ近づく動きから、2点を交互に訪れる動きへと、落ち着き方そのものが変わります。

パラメータの変化によって軌道の質的な構造が変わることを分岐と呼びます。安定な周期が2倍のものへ移る現象が周期倍分岐です。ただし、導関数の絶対値が1になったという計算だけで、一般の写像すべてに周期倍分岐が起こると断定はできません。ここではロジスティック写像の具体的な因数分解によって、新しい周期軌道が存在することまで確かめました。境界値を機械的に公式へ入れるのでなく、そこで判定が止まる理由を読むことが大学数学の練習になります。

単純な法則と複雑な軌道を結んだ研究

連続時間の研究では、流体の変化を少数の変数にまとめたモデルが重要な役割を果たしました。ローレンツは、決定論的な非線形常微分方程式と相空間の軌道を結びつけ、近い初期状態が大きく異なる状態へ発展しうることを論じています。また、対流を表す簡単な系を数値的に扱っています（Lorenz, 1963）。微分方程式の右辺が決まっていても、軌道が初期状態の微小な違いに対して不安定でありうるという問題が、予測の議論につながります。

離散時間では、メイが、生物・経済・社会科学に現れる1階の差分方程式を取り上げ、安定点から安定な周期の階層、見かけ上ランダムな変動までを示しうると論じました（May, 1976）。ここでいう1階とは、次の状態を直前の状態から決めるという意味です。式の見た目が簡単だから、その先の軌道も簡単だとは限らないのです。先ほどの不動点と2周期の計算は、その幅広い現象の一部分を自分の手で確かめたものです。

次の問いは、異なる非線形写像で似た周期倍分岐が見えるとき、その似方を数量で表せるかというものでした。ファイゲンバウムの1978年の論文は、非線形変換のあるクラスについて「定量的普遍性」を題名に掲げています（Feigenbaum, 1978）。この話題を扱うVulらの総説では、無限の周期倍分岐を考え、写像の反復と尺度の変更を繰り返すことで、変換に対して不変な極限の写像を捉える見方が説明されています。普遍定数 δ ≈ 4.6692 は、その変換を不動点の近くで線形化したときの固有値として現れます（Vul et al., 1984）。

この数の意味を、ロジスティック写像の計算で覗いてみましょう。臨界点 x = 1/2 自身が最小周期2、4、8の軌道に乗るパラメータを順に s₁、s₂、s₃ とします。各周期について f の反復が1/2に戻る方程式を数値的に解くと、s₁ ≈ 3.2360679775、s₂ ≈ 3.4985616993、s₃ ≈ 3.5546408628 となります。このとき (s₂-s₁)/(s₃-s₂) ≈ 4.680771 です。これらは周期が生まれる境界値そのものではなく、各安定周期の中で臨界点が周期軌道に乗る値です。

周期16、32へ計算を延ばすと、同じ形の隣り合う区間幅の比は約4.662960、4.668404になります。初めの数個でも、幅が縮む比率に規則性が見えます。ただし、この数値列だけで極限の存在や普遍性を証明したことにはなりません。普遍性は、どんな写像でも同じになるという主張ではなく、適切なクラスや条件の下で現れる現象です。写像そのものを変換の対象とし、その変換の不動点を調べることで、個々の軌道から一段上の構造へ進むのが、この研究の特徴です。

誤差の成長、対流のモデル、折り重なる図形

予測の限界を考える道具として、リアプノフ指数があります。微分可能な1次元写像について、軌道上で導関数が0でなく、次の極限が存在するとき、その軌道に沿う指数を次のように表せます。log は自然対数です。

λ = lim［n → ∞］(1/n)∑［j = 0,…,n-1］log|f′(xⱼ)|

連鎖律によって、n 回の更新の導関数は各時刻の導関数の積になります。積の対数は和になるため、この式は微小なずれの倍率を対数で平均したものです。λ > 0 なら、線形化が有効な間、ずれはおおよそ e の λn 乗の割合で成長します。λ < 0 なら微小なずれを縮める方向に働きます。極限が存在しない軌道や導関数が0になる軌道では、そのままこの有限値の式を使うことはできず、扱いを調整します。

計算しやすい例として、r = 4 の不動点 p = 3/4 を取ります。f′(p) = -2 なので、この不動点に沿う指数は λ = log 2 ≈ 0.693147 です。初期誤差が 0.000001 なら、線形化による20回後の見積もりは 2²⁰ × 0.000001 = 1.048576 になります。しかし状態が [0,1] 内にある以上、2つの状態の差は1を越えません。この数は実際の差の予言ではなく、途中で微小誤差の近似が破れると読むべき値です。

この例はもう一つの注意点も教えます。p 自身の軌道は、ずっと3/4にとどまる周期1の軌道です。指数が正であることは、その軌道が自ら不規則に動くことと同義ではありません。周囲からのずれが増幅されることと、軌道全体のカオス的な性質を区別する必要があります。長い数値列から指数を推定する場合も、過渡的な動き、観測の長さ、丸めの影響を分けて考えます。

連続時間の具体例として、ローレンツ方程式は次の形で書けます。

dx/dt = σ(y-x)

dy/dt = x(ρ-z)-y

dz/dt = xy-βz

σ、ρ、β は固定したパラメータです。状態は1個の実数でなく3成分の組 (x,y,z) であり、時間とともに動く点の軌跡を3次元の相空間で調べます。非線形項 xz と xy があるため、それぞれの成分を独立の等比的な増減として扱うことはできません。対流のモデルと初期状態への不安定性という接点が、ローレンツの研究の中心にあります（Lorenz, 1963）。

具体的に σ = 10、ρ = 28、β = 8/3 として、まず平衡点だけを計算しましょう。最初の式から y = x、2番目から x(27-z) = 0 です。x = 0 なら3番目から z = 0 となります。x ≠ 0 なら z = 27 で、3番目から x² = 72 です。従って平衡点は (0,0,0)、(√72,√72,27)、(-√72,-√72,27) となります。平衡点を求めたことと、その周囲の軌道がどう動くかを証明したことは別です。まず基準となる状態を確定し、その後に安定性や長期の軌道を調べます。

この系では、ベクトル場の発散は -σ-1-β で、今の値なら -41/3 です。これは流れに沿う微小体積を縮める性質を表します。ただし体積が縮むことと、すべての方向の距離が縮むことは違います。ある方向には伸び、別の方向にはそれ以上に縮むことがありえます。伸長と収縮に折り返しが組み合わさると、軌道が有限の領域にありながら複雑な配置を作る仕組みを考えられます。負の発散だけでカオスの存在を証明できるわけではありません。

そのような幾何を学ぶ入口がフラクタルです。代表例として、線分 [0,1] から中央の開いた3分の1を除き、残る各線分に同じ操作を繰り返す集合を考えます。n 段階後の線分は2ⁿ本で、各長さは1/3ⁿ、長さの合計は (2/3)ⁿ です。合計は0に近づきますが、極限の集合には端点などが残ります。全体が縮尺1/3の2つのコピーからなる自己相似性を持ちます。

この集合の自己相似次元 D は、2 × (1/3)ᴰ = 1 を満たすので D = log 2/log 3 ≈ 0.630930 です。この例では、自己相似性と部分の分離があり、ハウスドルフ次元もこの値になります。フラクタルという語には広い使い方がありますが、少なくとも「形が複雑だから次元が必ず整数でない」とは推論できません。力学系の極限集合を調べるときも、どの次元をどの条件で求めたかを明確にします。

生物の個体数のモデルでは、r を変えたときの安定性の違いが、定常状態を想定してよいかを検討する入口になります。メイの論考が扱うのも、簡単な差分方程式が多様な分野で複雑な振る舞いを生むという問題です（May, 1976）。ただし、観測された不規則な変動を見ただけでは、決定論的カオスなのか、外からの変動なのか、モデルの不足なのかは決まりません。モデルを作ること、数学的な性質を調べること、データとの対応を確かめることを、別々の問いとして進めます。

表・グラフ・証明を行き来する学び方

力学系でつまずきやすいのは、数値が合っていることを、その先の無限時間の保証と受け取る場面です。もう一つは、見た目が乱れたグラフを、そのままカオスの証明だと思うことです。ここまでの計算では、軌道の列、導関数、不動点の式がそれぞれ違う情報を与えました。これらを1枚のノートで対応させると、「何を見れば答えられる問いか」が整理できます。

ロジスティック写像を視覚化する教材として、SernaとJoshiは、区間 [0,1] を10等分し、反復で得た値の属する区間に対応するLEDを点灯させる仕組みを扱っています（Serna and Joshi, 2012）。研究の対象はマイクロコントローラによる可視化教材であり、特定の学年の学習者の成績向上を比較して示したものではありませんが、数値の移り変わりを位置の移動として捉える授業の発想につながります。機器がなくても、数直線に値を打ち、更新ごとに印をつける方法を考えられます。

McCartneyは、のこぎり波写像を基にする結合写像格子について、ネットワークのリアプノフ指数と行列の結果を結ぶ例を提示しています（McCartney, 2024）。対象はカオスや力学系の授業に用いる高次元写像の数学的教材であり、特定の国や学年で学習効果を検証した結果ではありませんが、数列の反復と線形代数をつなぐ学習課題として位置づけられます。計算の道具が増えるほど、以前の単元を別の役割で使えるようになります。

学習を進める手順は、次のように具体化できます。

1. 状態空間、更新式、初期値、パラメータを最初に書きます。x₀ を変える操作と r を変える操作を分け、何を比較しているかを固定します。特に [0,1] から出ない理由を、二次関数の最大値で説明してみます。
2. 数回の値を表にし、横軸を時刻、縦軸を状態として描きます。さらに y = f(x) と y = x を同じ平面に描き、出力を次の入力へ渡す反復を追います。前者は時間変化、後者は更新の仕組みを示す図です。
3. f(p) = p を解き、f′(p) を計算します。表から予想した収束先と、不動点の判定を照合します。導関数の絶対値が1の場合には、結論の欄に「この判定では未決定」と書き、その先の検討を残します。
4. 2周期を疑ったら、2回の合成を作り、不動点の解を取り除きます。最後に f(a) = b、f(b) = a を代入で確かめ、f′(a)f′(b) を調べます。候補を求める計算と、周期・安定性の確認を分けます。

高校生のうちには、一般項を求める練習に加えて、漸化式の値を手で追い、二次関数のグラフと結ぶ練習ができます。微分を学んだら、接線の傾きの大きさが入力誤差をどう変えるかを言葉にしてみてください。計算機を使う場合も、表示桁を変える前に、何の誤差を小さくしたいのかを決めます。初期値を正確に知れないことと、計算の丸めを減らすことは、別の問題だからです。

塾で担いたいのは、答えに付いた数字だけでなく、推論の切れ目を一緒に見つけることです。停電した部屋で開いたノートに「値が0.5に近い」と書いたなら、その横へ「この先も近づく理由」を置けるようにします。家庭では、お子さまに難しい公式を再現させるより、「どこまでは計算で、どこからは証明で確かめたの」と尋ねると、考えた筋道を話すきっかけになります。計算結果だけを急いで求めず、条件や途中式を説明する時間を持つことが支えになります。

ここで手が止まるなら、止まった単元まで一緒に戻ります。武蔵境の数学塾（1 対 1）／お問い合わせ

カオスを学ぶ前に残りやすい3つの疑問

決定論なら、もっと桁数を増やせば無限に先まで予測できますか

初期値と計算の精度を上げれば、予測できる時間を延ばせる場合があります。しかし初期値の不確かさが残り、それが増幅されるなら、同じ精度で際限なく先まで予測できるとはいえません。正のリアプノフ指数を用いる近似では、許される誤差を固定したとき、予測時間の伸びは初期誤差を減らした比の対数に対応します。精度を上げることは有効でも、予測時間が同じ比で伸びるとは限らないのです。

カオスと乱数は同じものですか

同じではありません。ここでの決定論的な写像は、同じ初期値から同じ状態列を作ります。一方、確率モデルは起こりうる結果とその確率を扱います。決定論的な系を統計的に調べることもありますが、それは両者が同じという意味ではありません。観測された短い数列の見た目だけで、生成の仕組みまで断定することは避けます。

フラクタルの図が描ければ、カオスを理解したことになりますか

図は構造を考える助けになりますが、それだけでは十分ではありません。反復回数を増やしたときに何が残るのか、どの集合に属する点を描いたのか、次元は何を意味するのかを説明する必要があります。中3分の1を除く集合なら、線分の本数と長さを追う計算が図を支えます。力学系でも、軌道の図と安定性の条件を行き来して理解を深めます。

力学系を学ぶと、数学が答えを1つ出すためだけの道具ではないことが見えてきます。初期値を少し変えても保たれるもの、パラメータを動かすと変わるもの、有限の計算からは決められないものを、それぞれの言葉で扱えます。出発点は、漸化式に値を代入する高校の計算でした。そこへ微分、極限、証明を重ねることで、安定な点から周期軌道、予測の限界までをつなげられます。決まった法則の先にある複雑さも、条件を明らかにする一歩から数学の問いになります。

• 前の回: 圏論　数学の構造どうしの関係を扱う言葉
• 特集「大学数学」全30回の一覧へ
• 数学の個別指導のご案内：武蔵境の数学塾（1 対 1）
• 高校生の指導のご案内：学年別個別指導コース
• 塾への相談・お問い合わせ：お問い合わせ

あわせて読む（数学の地図・高校数学の単元解説）

• 数学の地図 第19回 割合と比から指数的な増え方、微分方程式へ
• 数学の地図 第21回 文章題から数学的モデリング、応用数学へ

参考文献

1. Edward N. Lorenz（1963）「Deterministic Nonperiodic Flow」Journal of the Atmospheric Sciences. https://doi.org/10.1175/1520-0469(1963)020<0130:dnf>2.0.co;2
2. Robert M. May（1976）「Simple mathematical models with very complicated dynamics」Nature. https://doi.org/10.1038/261459a0
3. Mitchell J. Feigenbaum（1978）「Quantitative universality for a class of nonlinear transformations」Journal of Statistical Physics. https://doi.org/10.1007/bf01020332
4. E B Vul, Yakov G Sinai, K M Khanin（1984）「Feigenbaum universality and the thermodynamic formalism」Russian Mathematical Surveys. https://doi.org/10.1070/rm1984v039n03abeh003162
5. Juan D. Serna, Amitabh Joshi（2012）「Visualizing the Logistic Map with a Microcontroller」Physics Education. https://doi.org/10.1088/0031-9120/47/6/736
6. Mark McCartney（2024）「Chaos on the Hypercube and Other Places」International Journal of Mathematical Education in Science and Technology. https://doi.org/10.1080/0020739x.2022.2113466