※小説ではない※専門書 要約資料集 為替(換算)3.9万円でもらう 紐解集生成 専門 初入門 資料

# BOOK-0125 微分方程式 — 振動・相平面・数値解法・偏微分方程式(数学派生 第2巻)

# BOOK-0125 微分方程式 — 振動・相平面・数値解法・偏微分方程式(数学派生 第2巻)

> 学問の宇宙・形式科学の恒星群 数学冠・派生第2巻。ガイド役: Fable 5 監修 / Sonnet 5 執筆(脚本班)
> 読書系統: 前提=第1巻(→BOOK-0078『微分方程式』第1巻)。第1巻で `dy/dx=ky` の指数関数解、単振動 `m(d²x/dt²)=-kx` の三角関数解、オイラー法の入口、三体問題とカオスの入口までを学んだ読者を対象に、本冊(第2巻)ではその続きを扱う。
> 水準: 五〜六(後半水準)。専門用語は初出で必ず水準つき定義を行う。
> 接続方針: ファイルリンクは第1巻(→BOOK-0078)のみに限定する。線形代数・フーリエ解析・確率論など他分野の巻は題名のみを言及し、リンクは張らない(リンク切れ防止のため)。

---

# BOOK-0125 微分方程式 — 振動・相平面・数値解法・偏微分方程式(数学派生 第2巻)

> 学問の宇宙・形式科学の恒星群 数学冠・派生第2巻。ガイド役: Fable 5 監修 / Sonnet 5 執筆(脚本班)

> 読書系統: 前提=第1巻(→BOOK-0078『微分方程式』第1巻)。第1巻で `dy/dx=ky` の指数関数解、単振動 `m(d²x/dt²)=-kx` の三角関数解、オイラー法の入口、三体問題とカオスの入口までを学んだ読者を対象に、本冊(第2巻)ではその続きを扱う。

> 水準: 五〜六(後半水準)。専門用語は初出で必ず水準つき定義を行う。

> 接続方針: ファイルリンクは第1巻(→BOOK-0078)のみに限定する。線形代数・フーリエ解析・確率論など他分野の巻は題名のみを言及し、リンクは張らない(リンク切れ防止のため)。

---

## 入口の物語 — 「行ったり来たり」の奥にある、もう一つの世界

第1巻の終わりで、私たちは奇妙な事実に突き当たった。バネにおもりをつるすと `x(t) = A cos(ωt)` というきれいな三角関数で揺れる。ところが、たった3つの天体が重力で引き合うだけで、その未来はもう「きれいな式」では書き表せなくなる。単振動の式は変数が1つ(位置 x)で、力が単純だったからこそ解けた。しかし現実の振動には摩擦があり、現実の生態系には捕食者と被食者という複数の生き物が絡み合い、現実の熱や波は空間全体に広がって伝わっていく。

本冊(第2巻)では、単振動に「摩擦」を足したときに何が起こるか(減衰と共振)、変数が2つ以上に増えたときに何が起こるか(連立系と相平面)、解析的に解けない式でも精度よく未来を追いかける方法(ルンゲ=クッタ法)、空間にも広がっていく量を記述する偏微分方程式(へんびぶんほうていしき)の入口までを学ぶ。最後には、感染症・振り子・電気回路という一見異なる現象が、驚くほど似た数式の構造を共有していることを見る。

---

## 第五階段(水準五) — 2階線形方程式の一般解と単振動の完全な姿

### 一般解とは何か

第1巻では単振動の解を `x(t) = A cos(ωt)` という形で紹介した。しかしこれは「t=0で最大まで引っ張って手を離した」という特別な初期条件に限った解にすぎない。おもりに初速を与えて手を離すこともできる。この一般の場合を含めた完全な解が、次の**一般解(いっぱんかい、水準五: 微分方程式を満たす関数のうち、任意定数を含み、あらゆる初期条件に対応できる最も広い形の解)**である。

```

x'' = -ω²x の一般解: x(t) = A cos(ωt) + B sin(ωt)

```

A と B は初期条件によって決まる2つの任意定数である。2階の微分方程式には、一般に2つの独立な任意定数が必要になる。これは「積分を2回行う」ことに対応し、1回の積分ごとに1つの積分定数が生まれるという直感と一致する。

### 検算その1 — 初期条件から A と B を決め、解が方程式を満たすことを確認する

角振動数 `ω = 2`(第1巻の m=2kg, k=8N/m の設定を引き継ぐ)、初期条件を「t=0で位置 x=1、初速度 x'=2」とする。一般解を微分すると `x'(t) = -Aω sin(ωt) + Bω cos(ωt)` となるので、t=0 を代入すると次が成り立つ。

```

x(0) = A = 1 → A = 1

x'(0) = Bω = 2 → B = 2/ω = 2/2 = 1

```

したがって解は `x(t) = cos(2t) + sin(2t)` である。この関数が `x'' = -4x` を満たすか、t=0.3 で数値的に検算する。二階微分は第1巻と同じ近似式 `(x(t+h)-2x(t)+x(t-h))/h²`(h=0.0001)を用いる。

```

x(0.3) = cos(0.6) + sin(0.6) ≈ 0.825336 + 0.564642 ≈ 1.389978

数値二階微分 ≈ -5.559912

理論値(-ω²・x(0.3) = -4 × 1.389978) ≈ -5.559912

```

両者は小数第6位まで一致している(python3実測: 数値二階微分-5.5599123527372285、理論値-5.559912353218855。誤差はh=0.0001による近似誤差の範囲内)。単振動の一般解が任意の初期条件のもとで正しく方程式を満たすことが確認できた。

### 位相と振幅への書き換え

`A cos(ωt) + B sin(ωt)` は三角関数の合成公式を使うと、「振幅 R」と「位相のずれ φ(ファイ)」を使った形にも書き換えられる。

```

A cos(ωt) + B sin(ωt) = R cos(ωt - φ)

R = √(A² + B²), φ = atan2(B, A)

```

先ほどの例(A=1, B=1)では `R = √2 ≈ 1.41421`、`φ = atan2(1,1) = π/4 ≈ 0.7854` ラジアンとなる。振幅√2で揺れながら、通常のcos(ωt)よりπ/4だけ位相が進んでいるという意味である。この「振幅と位相のずれ」という見方は、交流回路や波の干渉にも姿を変えて現れる。

---

## 休憩所⑤(ここまでのまとめ箱)

- 2階線形方程式 `x''=-ω²x` の一般解は `x(t)=A cos(ωt)+B sin(ωt)`(A・Bは初期条件から決まる2つの任意定数)。ω=2, 初期条件x(0)=1, x'(0)=2からA=1, B=1が定まり、t=0.3での数値二階微分-5.559912が理論値-4×x(0.3)=-5.559912と一致することを検算済み。

- 一般解は `R cos(ωt-φ)` の形にも書き換えられ(R=√(A²+B²)、φ=atan2(B,A))、振幅と位相のずれという2つの量にまとめられる。

---

## 第六階段(水準五) — 減衰振動と共振、タコマナローズ橋の教訓

### 摩擦を足すとどうなるか — 減衰振動

現実のバネや振り子は摩擦によってやがて揺れが止まる。この摩擦を式に取り込むと、次のような**減衰振動(げんすいしんどう、水準五: 摩擦や抵抗によってエネルギーが失われ、振動の振幅がだんだん小さくなっていく運動)**の方程式が得られる。

```

m x'' + c x' + k x = 0

```

新しく加わった `c x'`(cは正の定数、減衰係数)は速度に比例して働くブレーキのような抵抗力である。`x = e^(rt)` という形を仮定して代入すると(→BOOK-0078第1巻の `y=Ce^(kx)` の自然な拡張)、次の**特性方程式(とくせいほうていしき、水準五: 線形微分方程式に指数関数の解を仮定して代入することで得られる、rについての代数方程式)**が得られる。

```

m r² + c r + k = 0

```

解の公式で解くと `r = (-c ± √(c²-4mk)) / (2m)` となる。ルートの中身、**判別式(はんべつしき、水準五: 2次方程式 ar²+br+c=0 における b²-4ac の値。この符号によって解の種類(実数解が2つ・1つ・複素数解)が決まる)** `D = c² - 4mk` の符号によって、振る舞いが3種類に分かれる。

### 3種類の減衰 — 過減衰・臨界減衰・不足減衰

m=1、k=4(固有角振動数 `ω₀=√(k/m)=2`)に固定し、cだけを変える。

```

c=2 のとき: D = 2²-4×1×4 = 4-16 = -12 (D<0)

c=4 のとき: D = 4²-4×1×4 = 16-16 = 0 (D=0)

c=8 のとき: D = 8²-4×1×4 = 64-16 = 48 (D>0)

```

**過減衰(かげんすい、水準五: D>0のとき。振動せずゆっくりとつり合いの位置へ戻る状態)** (c=8): `r₁≈-0.5359`、`r₂≈-7.4641` という2つの異なる負の実数解が得られる(python3実測)。解は `x(t)=C₁e^(r₁t)+C₂e^(r₂t)` の形になり、指数関数的に0へ収束する。蜂蜜の中でおもりを揺らすイメージである。

**臨界減衰(りんかいげんすい、水準五: D=0のとき。振動せず最も速くつり合いの位置へ戻る、過減衰と不足減衰の境界の状態)** (c=4): 重解 `r=-2.0` が得られ、解は `x(t)=(C₁+C₂t)e^(-2t)` の形になる。振動しない範囲で最も速く静止状態に戻る、ドアクローザーの設計に使われる考え方である。

**不足減衰(ふそくげんすい、水準五: D<0のとき。振幅が小さくなりながら振動を続けつり合いの位置に近づく状態)** (c=2): `r=-1.0±1.7321i` という複素数の解が得られる(iは虚数単位)。実部が減衰の速さ、虚部が振動の速さを表す。解は `x(t)=e^(-t)(C₁cos(1.7321t)+C₂sin(1.7321t))` で、最も日常的な「だんだん揺れが小さくなる」振動がこれである。

### 共振 — 外から力を加えたときの応答

摩擦とは逆に、外部から周期的な力 `F cos(Ωt)`(Ωは外力の角振動数)を加え続けるとどうなるか。

```

m x'' + c x' + k x = F cos(Ωt)

```

外力の角振動数Ωが系の固有角振動数`ω₀=√(k/m)`に近づくと振幅がどんどん大きくなる。これを**共振(きょうしん、水準五: 外部からの周期的な力の振動数が、系そのものが持つ固有の振動数に近いとき、振幅が著しく大きくなる現象)**と呼ぶ。ブランコを揺れの周期に合わせて押すと小さな力で大きく揺らせるのが、身近な共振の例である。

### タコマナローズ橋 — 「共振」という説明の限界(概略・諸説あり)

共振を語るとき引き合いに出されるのが、1940年11月7日にアメリカのワシントン州で起きた**タコマナローズ橋の崩落**である。多くの入門的な教材では「風が橋の固有振動数と一致した単純な共振の実例」として紹介されることがある。

しかし近年の工学的な調査ではこの説明はやや単純化されすぎているという見方が主流になっている。実際に起きたのは単純な共振ではなく、**自励振動(じれいしんどう、水準五: 周期的でない外力が構造物の運動と相互作用し、構造物自身が振動のエネルギーを外力から絶えず汲み上げて振幅が際限なく増大する現象)**、より専門的には空力フラッター(負の減衰)と呼ばれる現象だったという見方である。橋の断面形状が翼のような働きをし、ねじれ運動と風との相互作用で振幅が増幅されたと理解されている。

この点は工学者の間でも議論が続く**諸説あり**の話題であり、本冊では「1940年11月7日にタコマナローズ橋の主径間が崩落したこと」「教科書的には単純共振とされることがあるが、専門的な調査では自励振動(空力フラッター)によるものとする見方が主流であること」という大枠のみを断定し、詳細な機序には立ち入らない。

---

## 休憩所⑥(ここまでのまとめ箱)

- `mx''+cx'+kx=0` の減衰振動は判別式 `D=c²-4mk` の符号で3種類に分かれる。m=1,k=4固定で、c=8(D=48)は過減衰(r≈-0.536,-7.464)、c=4(D=0)は臨界減衰(重解r=-2.0)、c=2(D=-12)は不足減衰(r=-1.0±1.7321i)であることを検算済み。共振は外力の振動数が系の固有振動数に近づくと振幅が著しく増大する現象。

- タコマナローズ橋(1940年11月7日崩落)は教科書的には単純共振の例とされることがあるが、専門的には自励振動(空力フラッター)による崩落という見方が主流(諸説あり)。

---

## 第七階段(水準五〜六) — 連立線形系と相平面、ロトカ=ヴォルテラモデル

### 変数が2つになると何が変わるか

これまでの式はすべて変数が1つだった。しかし現実には複数の量が互いに影響し合いながら変化する場面が多く、キツネ(捕食者)とウサギ(被食者)の個体数は互いの数に応じて増減する。このような場合を扱うのが**連立微分方程式(れんりつびぶんほうていしき、水準五: 2つ以上の量が、互いの値に依存しながら同時に変化していく様子を表す、複数の微分方程式の組)**である。2変数の連立線形系は一般に次の形で書ける(x, yが変数、a,b,c,dは定数)。

```

dx/dt = a x + b y

dy/dt = c x + d y

```

この系の振る舞いを視覚的に理解する道具が**相平面(そうへいめん、水準五: 系の状態を表す複数の変数(たとえばxとy)を軸にとった平面。時間の経過とともに系がどう変化していくかを、この平面上の1本の軌跡として描き出せる)**である。時間tを軸にする代わりにx,yを軸にとり、(x(t), y(t))という点の軌跡を追うと、その形から系の性質を読み取れる。

### 固有値による分類

連立線形系の振る舞いは、係数行列 `[[a,b],[c,d]]` の**固有値(こゆうち、水準六: ある行列に対して、その行列を掛けても向きが変わらない特別な方向(固有ベクトル)に対する、伸び縮みの倍率を表す数)**によって分類できる。固有値rは単振動の特性方程式と同様、次の2次方程式(**特性方程式**)の解として求まる。

```

r² - (trace)r + (det) = 0

trace = a+d (行列の対角成分の和)

det = ad-bc (行列式)

```

固有値が実数か複素数か、符号がどうかによって、相平面上の軌跡は「安定な渦」「不安定な渦」「鞍点」「中心(閉じた楕円軌道を繰り返す)」などに分類される。この分類は次のロトカ=ヴォルテラモデルの理解に直結する。

### ロトカ=ヴォルテラの捕食者被食者モデル

**ロトカ=ヴォルテラ方程式(水準六: 被食者と捕食者、2種類の生物の個体数が互いに影響し合いながら周期的に増減する様子を表す、非線形の連立微分方程式)**は、アメリカの数理生物学者アルフレッド・ロトカが1925年に(化学反応の理論的研究から)、イタリアの数学者ヴィト・ヴォルテラが1926年に(アドリア海の漁獲データの分析から)、それぞれ独立に導いた。ヴォルテラの研究は第一次世界大戦中の漁獲データがきっかけとされる。

被食者(ウサギなど)の個体数をx、捕食者(キツネなど)の個体数をyとすると、モデルは次のように書ける。

```

dx/dt = a x - b x y (被食者: 自然増加 - 捕食される分)

dy/dt = -c y + d x y (捕食者: 自然減少 + 被食者を食べて増える分)

```

a, b, c, dはすべて正の定数。被食者はaの率で増えるが捕食者と出会う頻度(xyに比例)だけ減り、捕食者はcの率で減るが被食者と出会う頻度だけ増える、という関係を表している。

### 検算その2 — 非自明な平衡点と固有値、周期軌道の確認

この式にはx=y=0(両方とも絶滅)という自明な平衡点のほかに、**非自明な平衡点(ひじめいなへいこうてん、水準六: x=0,y=0以外で、dx/dt=dy/dt=0となる、系が動かなくなる特別な点)**が存在し、次のように求まる。

```

x* = c/d, y* = a/b

```

a=1.0, b=0.1, c=1.5, d=0.075 とすると(python3実測)、

```

x* = c/d = 1.5/0.075 = 20.0

y* = a/b = 1.0/0.1 = 10.0

```

平衡点でのヤコビ行列(各方程式をx, yでそれぞれ微分して並べた行列)を計算すると、平衡点の定義から対角成分が0になり、次のようになる。

```

J = [ 0 -b・x* ] = [ 0 -2.0 ]

[ d・y* 0 ] [ 0.75 0 ]

```

このJの固有値は、trace=0、det=1.5から特性方程式 `r² - 0・r + 1.5 = 0` を解いて `r = ±√(-1.5) = ±1.2247i` という**純虚数**が得られる(python3実測: 判別式=-6.0、固有値=0.0±1.224745i)。

固有値が実部0の純虚数になるということは、平衡点の近くで系が渦を巻いて収束も発散もせず、**中心(ちゅうしん、水準六: 相平面上で、軌跡が同じ形の閉じた曲線を何度も繰り返す平衡点の分類。エネルギーが増えも減りもしない保存系に典型的に現れる)**型の振る舞いをすることを意味する。ウサギとキツネの個体数は爆発も絶滅もせず、**周期的に増減を繰り返す**という生態学的に観察される現象を、この固有値の計算だけから予言できる。相平面上では(x*,y*)を中心とする閉じた楕円状の軌道を描き、時間変化としては被食者が増える→捕食者が遅れて増える→被食者が減り始める→捕食者も遅れて減る、というサイクルが繰り返される。

---

## 休憩所⑦(ここまでのまとめ箱)

- 連立微分方程式の振る舞いは相平面(状態変数を軸にとった平面上の軌跡)で視覚化でき、係数行列の固有値によって渦・鞍点・中心などに分類される。ロトカ=ヴォルテラ方程式(ロトカ1925年、ヴォルテラ1926年、独立に導出)は捕食者被食者の個体数変動を表す非線形連立系である。

- a=1.0,b=0.1,c=1.5,d=0.075のとき非自明な平衡点は(x*,y*)=(20.0,10.0)、その点でのヤコビ行列の固有値は純虚数±1.224745i(中心型)であることを検算済み。これは個体数が周期的に増減を繰り返すことに対応する。

---

## 第八階段(水準五〜六) — オイラー法を超えて、ルンゲ=クッタ法(RK4)

### オイラー法の弱点を思い出す

第1巻では`dy/dx=y, y(0)=1`をオイラー法で刻み幅h=0.5、2歩分計算し`y(1.0)≈2.25`(真値e≈2.71828との誤差は約17%)を得た。この誤差は、オイラー法が「現在の傾き1つだけ」を使って次の値を予測し、傾きが変化していく途中の様子を無視してしまうことに原因がある。

### ルンゲ=クッタ法(RK4)という改良

**ルンゲ=クッタ法(水準五: 1歩の間に複数の地点での傾きを計算し、それらを重み付けして平均することで、オイラー法よりもはるかに高い精度で次の値を予測する数値解法の一群。特に4次のものは実務で最も広く使われる)**は、ドイツの数学者カール・ルンゲとマルティン・クッタが20世紀初頭に開発した。実務で最もよく使われる**4次のルンゲ=クッタ法(通称RK4)**は、1歩の中で次の4か所の傾きを計算し加重平均する。

```

k1 = f(y) 現在地点での傾き

k2 = f(y + h/2・k1) k1を使って半歩進んだ地点での傾き

k3 = f(y + h/2・k2) k2を使って半歩進んだ地点での傾き

k4 = f(y + h・k3) k3を使って1歩丸ごと進んだ地点での傾き

y(次) = y(今) + (h/6)・(k1 + 2k2 + 2k3 + k4)

```

途中経過の傾きを1点(オイラー法)ではなく4点使い、中間の2点(k2, k3)を重視(重み2倍)して平均することで、傾きの変化をより正確に取り込めるのが要点である。

### 検算その3 — オイラー法とRK4の精度をy'=y, y(0)=1, h=0.1で比較する

`y'=y, y(0)=1, h=0.1`から出発する。オイラー法の1歩は次の通り。

```

オイラー法1歩: y(0.1) ≈ y(0) + h・f(y(0)) = 1 + 0.1×1 = 1.1

```

指定された検算値1.1と一致する(真値`e^0.1≈1.10517092`との誤差は約0.00517、相対誤差約0.47%)。同条件でRK4を1歩計算すると(python3実測)、

```

k1 = f(1) = 1

k2 = f(1+0.05×1) = 1.05

k3 = f(1+0.05×1.05) = 1.0525

k4 = f(1+0.1×1.0525) = 1.10525

y(0.1) = 1 + (0.1/6)×(1 + 2×1.05 + 2×1.0525 + 1.10525)

≈ 1.1051708333...

```

真値`e^0.1≈1.1051709181`と比べると、RK4の誤差はわずか約8.47×10⁻⁸であり、オイラー法の誤差(約0.00517)よりも**桁違いに小さい**。

さらに10歩分計算してt=1.0まで進めると(python3実測)、

```

真値: e ≈ 2.718281828

オイラー法(10歩,h=0.1): y(1.0) ≈ 2.593742460 誤差 ≈ 0.124539(約4.58%)

RK4(10歩,h=0.1): y(1.0) ≈ 2.718279744 誤差 ≈ 0.00000208(約0.0000766%)

```

同じ刻み幅・歩数でありながら、RK4の誤差はオイラー法のおよそ6万分の1に抑えられている。傾きの計算回数を1回から4回に増やす代わりに得られる大きな精度向上であり、天気予報やロケットの軌道計算など実務ではRK4かその発展形が標準的に使われている。

---

## 休憩所⑧(ここまでのまとめ箱)

- ルンゲ=クッタ法(RK4)は1歩の間に4か所(k1〜k4)の傾きを計算し重み付き平均することで、オイラー法より格段に高精度な数値解法である。`y'=y, y(0)=1, h=0.1`でオイラー法1歩は1.1(誤差約0.00517)、RK4は1歩で約1.10517083(誤差約8.47×10⁻⁸)であることを検算済み。

- 同条件で10歩進めてt=1.0にすると、オイラー法誤差≈0.1245に対しRK4誤差≈0.0000021と、RK4の誤差はオイラー法のおよそ6万分の1。

---

## 第九階段(水準六) — 偏微分方程式の入口、熱と波

### 変化するのは時間だけではない

ここまでの微分方程式は「時間tに対する変化」だけを記述していた(常微分方程式)。しかし熱いフライパンの温度は「どの場所か」によっても違い、ギターの弦の振動も位置によって変位が異なる。このように複数の変数にわたる変化を扱うのが**偏微分方程式(へんびぶんほうていしき、水準六: 2つ以上の独立変数(たとえば時間と空間の位置)を持つ関数の、それぞれの変数についての変化率(偏微分)の関係を表す方程式)**である。「偏微分」とは複数ある変数のうち1つだけに注目し他を固定して微分する意味であり、温度`u(x,t)`で「xを固定してtだけで微分する」ことを`∂u/∂t`と書く。

### 熱伝導方程式

もっとも基本的な偏微分方程式の一つが、**熱伝導方程式(ねつでんどうほうていしき、水準六: 物体の中で熱がどのように時間とともに空間へ広がっていくかを表す偏微分方程式)**である。

```

∂u/∂t = α (∂²u/∂x²)

```

u(x,t)は位置x・時刻tにおける温度、αは熱の伝わりやすさを表す正の定数(熱拡散率)である。「ある地点の温度の時間変化は、その地点の温度分布の“曲がり具合”に比例する」という意味で、周りより温度が低い地点は上がり、高い地点は下がる、「でこぼこをならす」働きを表している。

### 検算その4 — 差分法による1ステップの数値計算

数値的に確かめるため、位置を5点に離散化し時間もΔtで進める**差分法(さぶんほう、水準六: 微分を、ごく短い間隔での値の差(差分)で近似することで、偏微分方程式を数値的に解く手法)**を使う。空間の2階微分と更新式は次の通り(rはα・Δt/Δx²)。

```

∂²u/∂x² ≈ (u[i+1] - 2u[i] + u[i-1]) / Δx²

u[i]^{次} = u[i] + r・(u[i+1] - 2u[i] + u[i-1]), r = α・Δt/Δx²

```

5点(両端は温度0固定=境界条件)、中央の1点だけ最初に温度1、他は0、α=1、Δx=1、Δt=0.4(r=0.4)として1ステップ計算すると(python3実測)、初期状態`[0,0,1,0,0]`が1ステップ後`[0,0.4,0.2,0.4,0]`になる。中央の熱の一部が両隣に伝わり(0.4ずつ)、中央自身は減った(1→0.2)。「でこぼこをならす」働きの数値的な現れである。この種の差分計算には**安定性条件**(この場合`r≦0.5`)があり、超えると計算が破綻することが知られている(今回のr=0.4は条件を満たす)。詳しい導出は本冊の範囲を超えるため「上限がある」という事実のみを断定する。

### 波動方程式とフーリエの貢献

もう一つの基本的な偏微分方程式が、弦の振動や音・光の伝わり方を表す**波動方程式(はどうほうていしき、水準六: ある量が空間を一定の速さで伝わっていく波の運動を表す偏微分方程式)**である。

```

∂²u/∂t² = c² (∂²u/∂x²)

```

熱伝導方程式と似ているが左辺が時間について2階微分(加速度に相当)になっている点が異なり、熱伝導方程式の解が「でこぼこを滑らかにならしながら広がる」拡散的な振る舞いを見せるのに対し、波動方程式の解は「形を保ったまま伝わる」波としての振る舞いを見せる。cは波が伝わる速さを表す。

これらの偏微分方程式、とりわけ熱伝導方程式の解法に大きく貢献したのが、フランスの数学者・物理学者**ジョゼフ・フーリエ(1768年〜1830年)**である。フーリエは1807年頃から熱伝導の研究を進め、1822年に代表作『熱の解析的理論(Théorie analytique de la chaleur)』を出版した。この著作で、複雑な温度分布や波形を単純なサイン・コサイン(正弦・余弦)の波の足し合わせで表現できる**フーリエ級数**を体系化した。音声や画像処理など現代技術の土台にもなっている(詳しい仕組みは別巻に譲る)。

---

## 休憩所⑨(ここまでのまとめ箱)

- 偏微分方程式は時間と空間など複数の変数にわたる変化を記述する。熱伝導方程式 `∂u/∂t=α∂²u/∂x²` は「でこぼこをならす」拡散的な振る舞いを表し、差分法(α=1,Δx=1,Δt=0.4,r=0.4)による1ステップ計算で中央温度1の分布が[0,0.4,0.2,0.4,0]へ変化することを検算済み。

- 波動方程式 `∂²u/∂t²=c²∂²u/∂x²` は形を保ったまま伝わる波の運動を表す。フーリエは1822年『熱の解析的理論』でフーリエ級数を体系化し、熱伝導方程式の解法に貢献した。

---

## 第十階段(水準六) — ロジスティック方程式と個体群成長、そしてカオスへの扉

### 指数関数的増加の限界

第1巻で見たマルサスの人口論`dy/dt=ky`は人口が際限なく増え続けることを前提にしていたが、現実には食料や生息地に限りがあり、個体数がある上限に近づくと増加は頭打ちになる。これを表すのが**ロジスティック方程式(水準六: 個体数が増えるほど増加率が鈍り、上限(環境収容力)に近づくにつれ増加がほぼ止まる、個体群成長を表す非線形微分方程式)**である。

```

dP/dt = r P (1 - P/K)

```

Pは個体数、rは本来の増加率、Kは**環境収容力(かんきょうしゅうようりょく、水準六: その環境が持続的に養うことのできる、個体数の上限)**と呼ばれる定数である。`(1-P/K)`という項がミソで、Pがゼロに近いときはこの項がほぼ1(マルサス的な指数関数的増加)、PがKに近づくとほぼ0になり増加率も0に近づく。

### 検算その5 — S字曲線を数値で追いかける

解析解は次の式で表される。

```

P(t) = K / (1 + ((K-P₀)/P₀)・e^(-rt))

```

K=1000、P₀=10、r=0.3として計算すると(python3実測)、

```

t=0: P=10.000

t=10: P=168.665

t=20: P=802.957

t=30: P=987.930

t=40: P=999.392

t=50: P=999.970

```

最初はゆっくり増え(t=0→10)、中盤で急激に増え(t=10→20で634以上の増加)、終盤は再びゆっくり1000へ収束する。この「ゆっくり→急激→再びゆっくり」の軌跡はグラフに描くとSの字になることから**S字曲線(ロジスティック曲線)**と呼ばれる。増加が最も急激になる**変曲点(へんきょくてん、水準六: 曲線の曲がる向きが凸から凹へ切り替わる点)**はP=K/2(この例ではP=500)のときに訪れ、`t=ln((K-P₀)/P₀)/r≈15.317`で`P(15.317)≈500.000`と理論値に一致する(python3実測)。

### ロジスティック写像 — 離散版という別物への変身

ここまでは「時間が連続的に流れる」モデルだったが、1年に1度世代交代する昆虫のように「時間が飛び飛び(離散的)」な現象を扱いたい場合もある。この場合に使われるのが**ロジスティック写像(水準六: ロジスティック方程式の考え方を、連続時間ではなく1世代ごとの離散的なステップに書き直した漸化式。同じ名前の由来を持つが、連続版とは全く異なる複雑な振る舞い(カオス)を示すことで知られる)**である。

```

x_{n+1} = r・x_n・(1 - x_n)

```

`x_n`は0から1に規格化されたn世代目の個体数(の割合)、rは成長率パラメータ。連続版とよく似た形だが**離散的な漸化式という点で数学的には別物**であり、rを変化させると連続版には見られない劇的な振る舞いを見せる。

### 検算その6 — rを変えると振る舞いが激変する

初期値`x₀=0.2`から、rの値を変えて十分に世代を重ねた後の様子を確認する(python3実測)。

```

r=2.8: 20世代後、値は0.64付近に落ち着く(1つの値への収束)

r=3.2: 50世代後、値は0.51304と0.79946を交互に繰り返す(周期2)

r=3.5: 100世代後、値は0.50088,0.875,0.38282,0.82694を周期的に繰り返す(周期4)

```

rを2.8→3.2→3.5と増やすと「1つの値に収束」→「2つの値を往復」→「4つの値を巡回」と周期が2倍2倍に増えていく。これは**分岐(ぶんき、水準六: パラメータの値がある閾値を超えるたびに、系の周期的な振る舞いが2倍・4倍・8倍…と分裂していく現象)**と呼ばれる。

### カオスへ — 初期値のごくわずかな違いが運命を分ける

r=3.9まで増やすと周期的な繰り返しは崩れ、一見ランダムな振る舞いが現れる。第1巻の三体問題で見た「初期条件への鋭敏性」と同じ現象を、この単純な1行の式で直接確かめられる。

初期値をわずか0.0001だけ変えた2つの計算(`x₀=0.5`と`x₀=0.5001`)を、r=3.9で30世代進めると(python3実測)、

```

世代0: 差 = 0.0001000 (最初はほぼ同じ)

世代10: 差 = 0.0000054 (一時的に縮まる)

世代20: 差 = 0.0036582

世代30: 差 = 0.6892401 (30世代でほぼ無関係な値に)

```

最初はごくわずかだった違いが、30世代後には0.6892という「ほぼ無関係な2つの値」と言えるほどの差に育っている。これは**決定論的カオス(けっていろんてきかおす、水準六: サイコロのような偶然の要素を一切含まない、完全に決まったルール(漸化式)に従っているにもかかわらず、初期条件のごくわずかな違いが時間とともに指数関数的に拡大し、長期の予測を事実上不可能にする現象)**の、最も単純で有名な実例の一つである。

**注意点**: このカオス的振る舞いは**離散的な漸化式**の性質である。冒頭の**連続的なロジスティック方程式**`dP/dt=rP(1-P/K)`は1変数の常微分方程式としては滑らかなS字曲線に収束するだけでカオスは生じない(カオスには最低でも3変数以上の連続系、あるいは離散系が必要)。同じ「ロジスティック」の名を持ちながら連続版と離散版で運命が異なるこの事実は、連続(微分方程式)と離散(漸化式)という2つの言語の違いを際立たせる好例である。

---

## 休憩所⑩(ここまでのまとめ箱)

- ロジスティック方程式(連続)`dP/dt=rP(1-P/K)`はK=1000,P₀=10,r=0.3のときt=0→10→20→30→40→50でP≈10→168.665→802.957→987.930→999.392→999.970とS字曲線を描いて収束することを検算済み(変曲点t≈15.317でP=500)。

- ロジスティック写像(離散)`x_{n+1}=r x_n(1-x_n)`はr=2.8で周期1、r=3.2で周期2(0.51304/0.79946)、r=3.5で周期4に分岐することを検算済み。r=3.9では初期値0.0001の差が30世代で0.6892の差にまで拡大する決定論的カオスが生じる(検算済み)が、連続版のロジスティック方程式そのものにはカオスは生じない(離散と連続の違い)。

---

## 第十一階段(水準六) — 微分方程式で世界を写す、モデルの同型性

### SIRモデル — 感染症の広がりを写す(概略)

感染症が集団の中でどう広がるかを記述する基本的な数理モデルが**SIRモデル(水準六: 集団を「感受性者(Susceptible、まだ感染していない人)」「感染者(Infectious、感染していて他人にうつしうる人)」「回復者(Recovered、回復して免疫を得た、あるいは離脱した人)」の3グループに分け、その人数の移り変わりを連立微分方程式で表す感染症モデル)**である。

```

dS/dt = -β S I / N

dI/dt = β S I / N - γ I

dR/dt = γ I

```

Nは総人口(S+I+R=Nで一定)、βは感染のしやすさ、γは回復のしやすさを表す定数。比`R₀=β/γ`は**基本再生産数(きほんさいせいさんすう、水準六: 1人の感染者が、免疫を持たない集団の中で平均して何人にうつすかを表す数)**と呼ばれ、感染拡大か収束かの分かれ目になる指標として知られている。このモデルは先ほどのロトカ=ヴォルテラモデルと同じ「掛け算の項によって2グループが影響し合う」骨格を共有している。β=0.3, γ=0.1(R₀=3.0)、N=1000、最初の感染者1人という設定で計算すると(python3実測)、S+I+Rの合計が常に1000に保たれたまま(保存則)、感染者数Iはday0の1人からday30には161.4人、day60には114.1人へと山型の推移を描く。詳細は本冊では構造の概略のみを扱う。

### 振り子とバネ、そして電気回路 — 同じ式の3つの衣装

第1巻の「バネの振動」と「振り子の揺れ」が同じ数式構造(`ω²=k/m`と`ω²=g/L`)を共有していたのと同様に、電気回路にも同じ構造を持つ現象がある。**RLC回路(水準六: 抵抗R・コイルL・コンデンサCを直列につないだ、電気の振動を扱う基本的な回路)**を流れる電荷`q(t)`は次の式に従う。

```

L (d²q/dt²) + R (dq/dt) + (1/C) q = 0

```

第六階段の減衰振動`m x'' + c x' + k x = 0`と見比べると、

```

機械振動: m x'' + c x' + k x = 0

電気回路: L q'' + R q' + 1/C q = 0

対応関係: m(質量) ↔ L(インダクタンス、コイルの性質)

c(減衰係数) ↔ R(電気抵抗)

k(バネ定数) ↔ 1/C(コンデンサの容量の逆数)

```

数式が完全に同一ということは、片方で成り立つ性質がもう片方にもそのまま当てはまることを意味する。これを**モデルの同型性(どうけいせい、水準六: 見た目や分野がまったく異なる2つの現象が、記述する微分方程式の数式としての構造が完全に一致していること)**と呼ぶ。L=1H(ヘンリー)、R=2Ω(オーム)、C=0.25F(ファラド)とすると`1/C=4.0`(機械振動のk=4に相当)。固有角振動数は`ω₀=√(1/(LC))=√4=2.0`(python3実測、対応関係経由の計算とも一致)。判別式に相当する`R²-4L/C=4-16=-12`は負であり、第六階段のc=2の不足減衰と完全に同じ状況である。

つまり抵抗2Ω・コイル1H・コンデンサ0.25FのRLC回路は、m=1kg、c=2、k=4のバネ-おもり系と数式上まったく同じ振る舞いをする。電気振動の波形を、バネにおもりをつるした実験で目に見える形で再現できる、と言い換えてもよい。バネとおもり・振り子・RLC回路という3つの世界が、係数の呼び名こそ違え、すべて同じ`x''+2x'+4x=0`型の数式構造に支配されているのである。

---

## 休憩所⑪(ここまでのまとめ箱)

- SIRモデルは感受性者・感染者・回復者の3区分を連立微分方程式で表す感染症モデル(概略)。基本再生産数R₀=β/γが感染拡大の分かれ目となる指標として知られる。β=0.3,γ=0.1,N=1000の設定でS+I+R=1000の保存則が常に成立することを検算済み。

- RLC回路の式`Lq''+Rq'+(1/C)q=0`は機械振動の式`mx''+cx'+kx=0`と完全に同型(L↔m, R↔c, 1/C↔k)。L=1H,R=2Ω,C=0.25Fの場合、ω₀=√(1/(LC))=2.0、判別式相当R²-4L/C=-12<0(不足減衰)であり、これはm=1,c=2,k=4のバネ系と数式上完全に同じ振る舞いをすることを検算済み。

---

## ノウハウ集 — 手で解けない方程式への3つの実践解

微分方程式の多くは指数関数や三角関数のようなきれいな式で解けない。ここでは「手で解けない」方程式への3つの実践的アプローチを、材料・工程・なぜ効くか・限界の形でまとめる。

### 実践解1: 定性分析(相平面と平衡点の分類)

- **材料**: 微分方程式の右辺、係数行列(連立系の場合)。

- **工程**: ①平衡点を求める。②ヤコビ行列(特性方程式)を計算する。③固有値の符号・実部虚部から渦・鞍点・中心などに分類する。

- **なぜ効くか**: 方程式を解かなくても長期的な振る舞いは固有値の計算だけで分かることが多い。本冊のロトカ=ヴォルテラの検算(固有値が純虚数→中心→周期軌道)がその実例である。

- **限界**: 平衡点の「近く」での局所的な振る舞いしか分からず、非線形性が強い系では大域的な描像とずれることもある。

### 実践解2: 数値解法(オイラー法・ルンゲ=クッタ法)

- **材料**: 初期条件、変化率を計算する関数、刻み幅h。

- **工程**: ①刻み幅hを決める。②オイラー法かRK4で1歩ずつ進める。③精度を上げたい場合はhを小さくするか高次手法に切り替える。

- **なぜ効くか**: 「今の値と傾きから次の値を予測する」操作は常に実行できる。本冊の検算(RK4がオイラー法より桁違いに高精度)の通り、コストと精度のバランスを取りながら数値の推移を得られる。

- **限界**: あくまで近似値であり誤差が常にある。カオス的な系では初期条件や丸め誤差が指数関数的に拡大し、長期予測はできない。

### 実践解3: 近似(線形化・摂動的な考え方)

- **材料**: 複雑な非線形の方程式、注目したい平衡点や特定の状況。

- **工程**: ①注目する点の近くだけを考え、非線形の項を直線的な近似に置き換える。②線形化された簡単な方程式を解く。③元の方程式の「近くの振る舞い」をどれだけ再現しているか検討する。

- **なぜ効くか**: 本冊のヤコビ行列の計算も「平衡点の近くだけを線形近似する」考え方の一種であり、範囲を十分絞れば単純な線形方程式として扱える場面が多い。

- **限界**: 近似が成り立つのは注目した点の「近く」に限られ、前提が崩れる状況では的外れな結論を導く。

---

## なぜ面白いのか — 摩擦・複数・空間・偶然、4つの拡張がもたらす豊かさ

第1巻では「変化の規則を書く」という微分方程式の基本文法を学んだ。本冊ではその文法に4つの拡張を加えた。摩擦を1つ足すだけで単振動は3つの運命に枝分かれし、外力を加えれば共振という劇的な振る舞いが現れた。変数を2つに増やせばウサギとキツネの数が周期的に増減するロトカ=ヴォルテラの渦が姿を現し、空間という次元を追加すれば熱の拡散と波の伝播が生まれ、1行の漸化式にパラメータの変化を1つ許しただけで周期運動からカオスへの扉が開いた。そして感染症の広がりも電気回路の振動も、根っこの数式はバネや振り子と同じ骨格を共有している。少数のパターンの組み合わせで自然界の多くの現象が説明できるという事実こそが、微分方程式という言語の懐の深さである。摩擦・複数・空間・偶然という拡張を経てもなお、→BOOK-0078第1巻の`dy/dx=ky`の精神は生き続けている。

---

## 次への一歩 — この先に待っている風景

本冊では水準五〜六にあたる後半の基礎──2階線形方程式の一般解、減衰振動と共振、相平面・固有値による分類(ロトカ=ヴォルテラ)、数値解法の高度化(RK4)、偏微分方程式の入口、ロジスティック方程式とカオス、SIRモデルやRLC回路との同型性──を扱った。この先には厳密なカオス理論(ストレンジアトラクタ、リアプノフ指数)、高次元の偏微分方程式、確率微分方程式など、まだ見ぬ風景が広がる(後続の巻に委ねる)。→BOOK-0078第1巻の「未来を予言する式」という入口から始まった旅は、本冊でその適用範囲を大きく広げた。

---

## 参照リンク一覧(整合チェック用)

- →BOOK-0078『微分方程式』第1巻(単振動`x(t)=A cos(ωt)`の基礎・オイラー法・三体問題とカオスの入口。読書系統の前提として全篇で参照)

（線形代数・フーリエ解析・数値解析などの他分野の巻は、リンク切れ防止のため題名のみの言及とし、ファイルリンクは張らない。）

---

# BOOK-0125 微分方程式 — 振動・相平面・数値解法・偏微分方程式(数学派生 第2巻)