MathLabs

微分方程式と力学系

一階微分方程式

関数とその一階導関数を関係づける方程式で、変数分離法や積分因子法によって解かれる。

直観今、それはどれくらいの速さで変化しているのか?

多くの現象は、量そのものの式によってではなく、その量がどれくらいの速さで変化しているかという規則によって最も自然に記述される:個体群はその時点での大きさに比例して増加し、熱いコーヒーは室温よりどれだけ熱いかに比例した速さで冷め、充電中のコンデンサの電圧は満充電に近づくほどゆっくり上昇する。一階微分方程式とはまさにこのような規則であり——未知関数 y(x)y(x) をその導関数 y′(x)y'(x) 自身と関係づける——それを解くとは、変化率の規則から関数全体を復元することを意味する。

直接積分可能な一階方程式の右辺を示すインタラクティブな3次曲線ウィジェット。
1階常微分方程式 y′=y−cx2y' = y - c x^2 の方向場と、異なる初期値 y(t0)=y0y(t_0) = y_0 から出発する3本の解曲線。

大学変数分離形と線形方程式

定義: 変数分離形方程式

一階方程式が変数分離形であるとは、dydx=f(x)g(y)\frac{dy}{dx}=f(x)g(y) の形に書けることをいう:右辺が xx のみに依存する部分と yy のみに依存する部分の積に分解される。yy を含む項をすべて片側に、xx を含む項をすべてもう片側に移すと、方程式は2つの通常の積分 ∫dyg(y)=∫f(x) dx\int \frac{dy}{g(y)}=\int f(x)\,dx に帰着する。

dydx=f(x)g(y)\frac{dy}{dx}=f(x)g(y)

例えば、自励系方程式 dydx=g(y)\frac{dy}{dx}=g(y)(変化率が xx ではなく現在の値 yy のみに依存する場合)は f(x)=1f(x)=1 の特別な場合である。その平衡解はちょうど g(y)=0g(y)=0 となる箇所に現れ、これらの定数解は不用意に g(y)g(y) で割ると見落としやすい。

dydx+P(x)y=Q(x)\frac{dy}{dx}+P(x)y=Q(x)
一階方程式の種類
種類標準形解法
変数分離形dydx=f(x)g(y)\frac{dy}{dx}=f(x)g(y)両辺を積分する:∫dyg(y)=∫f(x) dx\int \frac{dy}{g(y)}=\int f(x)\,dx
線形一階dydx+P(x)y=Q(x)\frac{dy}{dx}+P(x)y=Q(x)μ(x)=e∫P(x) dx\mu(x)=e^{\int P(x)\,dx} を掛けてから積分する
自励系dydx=g(y)\frac{dy}{dx}=g(y)変数分離形の特別な場合;平衡点は g(y)=0g(y)=0

大学解くこと、そして解の存在を保証すること

P(x)P(x) と Q(x)Q(x) がある区間で連続であり、μ(x)=e∫P(x) dx\mu(x)=e^{\int P(x)\,dx} とすると、線形方程式 dydx+P(x)y=Q(x)\frac{dy}{dx}+P(x)y=Q(x) のすべての解は y(x)=1μ(x)(∫μ(x)Q(x) dx+C)y(x) = \frac{1}{\mu(x)}\left(\int \mu(x) Q(x)\,dx + C\right) で与えられる。

なぜ正しいのか?

この技巧は平方完成に似ている:方程式全体に、うまく選んだ因子 μ(x)\mu(x) を掛けることで、左辺が単一の積 μ(x)y(x)\mu(x)y(x) の導関数にまとまり、微分方程式を直接積分できるものに変える。

証明

μ(x)=e∫P(x) dx\mu(x)=e^{\int P(x)\,dx} とおく;連鎖律で微分すると、指数 ∫P(x) dx\int P(x)\,dx の導関数が P(x)P(x) であることから μ′(x)=P(x)μ(x)\mu'(x)=P(x)\mu(x) が得られる。

dydx+P(x)y=Q(x)\frac{dy}{dx}+P(x)y=Q(x) の両辺に μ(x)\mu(x) を掛けると、μ(x)y′+μ(x)P(x)y=μ(x)Q(x)\mu(x)y' + \mu(x)P(x)y = \mu(x)Q(x) となる。

積の法則により ddx[μ(x)y(x)]=μ(x)y′(x)+μ′(x)y(x)\frac{d}{dx}[\mu(x)y(x)] = \mu(x)y'(x) + \mu'(x)y(x) であり、μ′(x)=P(x)μ(x)\mu'(x)=P(x)\mu(x) なので、これは上式の左辺 μ(x)y′+μ(x)P(x)y\mu(x)y' + \mu(x)P(x)y にちょうど等しい。したがって方程式は ddx[μ(x)y(x)]=μ(x)Q(x)\frac{d}{dx}[\mu(x)y(x)] = \mu(x)Q(x) となる。

両辺を xx について積分すると、任意定数 CC に対して μ(x)y(x)=∫μ(x)Q(x) dx+C\mu(x)y(x) = \int \mu(x)Q(x)\,dx + C が得られる。μ(x)\mu(x)(指数関数なので決してゼロにならない)で両辺を割ると y(x)=1μ(x)(∫μ(x)Q(x) dx+C)y(x) = \frac{1}{\mu(x)}\left(\int \mu(x) Q(x)\,dx + C\right) が得られる。

f(x,y)f(x,y) が (x0,y0)(x_0,y_0) の周りの長方形上で連続であり、yy に関してリプシッツ条件 ∣f(x,y1)−f(x,y2)∣≤L∣y1−y2∣|f(x,y_1)-f(x,y_2)|\le L|y_1-y_2|(ある定数 LL について)を満たすとする。このとき初期値問題 y′=f(x,y), y(x0)=y0y'=f(x,y),\ y(x_0)=y_0 は、x0x_0 を含むある区間上で一意な解を持つ。

なぜ正しいのか?

リプシッツ条件は ff が yy 方向にどれだけ急峻に変化しうるかを制限し、これにより近くの解の軌道が分裂したり交差したりするのを防ぐ。これはまさに、各点をちょうど一つの軌道が通る、明確に定義された流れを保証するために必要な条件である。

証明

初期値問題を積分方程式に書き換えると、連続関数 y(x)y(x) が y′=f(x,y), y(x0)=y0y'=f(x,y),\ y(x_0)=y_0 を解くことは y(x)=y0+∫x0xf(t,y(t)) dty(x) = y_0 + \int_{x_0}^x f(t,y(t))\,dt を解くことと同値である。この同値性は微積分学の基本定理から従う。

連続関数に対するピカール反復作用素 TT を T[y](x)=y0+∫x0xf(t,y(t)) dtT[y](x) = y_0 + \int_{x_0}^x f(t,y(t))\,dt で定義し、仮定が成り立つ長方形内に ff がとどまるように選んだ小さな区間 I=[x0−h,x0+h]I=[x_0-h,x_0+h] に制限する。

II 上の2つの連続関数 y1,y2y_1,y_2 に対して、リプシッツ条件より ∣T[y1](x)−T[y2](x)∣≤∫x0x∣f(t,y1(t))−f(t,y2(t))∣ dt≤L h sup⁡t∈I∣y1(t)−y2(t)∣|T[y_1](x)-T[y_2](x)| \le \int_{x_0}^x |f(t,y_1(t))-f(t,y_2(t))|\,dt \le L\,h\,\sup_{t\in I}|y_1(t)-y_2(t)| となる。Lh<1Lh<1 となるほど hh を小さく選べば、TT は上限ノルムにおける縮小写像となる。

バナッハの不動点定理により、完備距離空間(上限ノルムを備えた II 上の連続関数の空間)上の縮小写像はちょうど一つの不動点 yy を持ち、y=T[y]y=T[y] を満たす。この不動点こそが、II 上での元の初期値問題の一意な解である。

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

一階方程式は、限られた資源に制限された個体群成長(生態学・疫学におけるロジスティック成長)、電子工学における充放電するRC回路、熱力学におけるニュートンの冷却法則、物理学における放射性崩壊や炭素年代測定、金融における時間変化する利率での連続複利などをモデル化する。以下の2つの具体例は、最も一般的な2つのパターンを最後まで解く。

例: ロジスティック人口成長

個体群 P(t)P(t) がロジスティックモデル dPdt=rP(1−PK)\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right) に従って成長し、初期個体数 P(0)=P0P(0)=P_0 であるとする。P(t)P(t) を求めよ。

解答

この方程式は変数分離形である:dPP(1−P/K)=r dt\frac{dP}{P(1-P/K)}=r\,dt と書き直せる。

部分分数分解 1P(1−P/K)=1P+1/K1−P/K\frac{1}{P(1-P/K)}=\frac1P+\frac{1/K}{1-P/K} により左辺を項ごとに積分でき、ln⁡∣P∣−ln⁡∣1−P/K∣=rt+C1\ln|P|-\ln|1-P/K|=rt+C_1、すなわち定数をまとめて KK を掛けた後 ln⁡∣PK−P∣=rt+C2\ln\left|\frac{P}{K-P}\right|=rt+C_2 が得られる。

指数をとると、ある定数 AA に対して PK−P=Aert\frac{P}{K-P}=Ae^{rt} となる。t=0t=0 で P(0)=P0P(0)=P_0 を用いると A=P0K−P0A=\frac{P_0}{K-P_0} が得られる。

PP について解いて整理すると閉じた形 P(t)=KP0ertK+P0(ert−1)P(t) = \frac{K P_0 e^{rt}}{K + P_0(e^{rt}-1)} が得られ、これは P0P_0 付近から始まり t→∞t\to\infty で環境収容力 KK に近づく、実際の個体群データで見られるまさにS字曲線である。

例: RC回路の充電

定電圧源 VV から抵抗を通して充電されるコンデンサは dQdt+QRC=VR\frac{dQ}{dt} + \frac{Q}{RC} = \frac{V}{R}(Q(0)=0Q(0)=0)を満たす。ここで Q(t)Q(t) は時刻 tt における電荷である。Q(t)Q(t) を求めよ。

解答

これは P(t)=1RCP(t)=\frac{1}{RC}、右辺項 VR\frac{V}{R} を持つ線形方程式であり、積分因子は μ(t)=et/(RC)\mu(t)=e^{t/(RC)} である。

両辺に掛けると ddt[et/(RC)Q]=VRet/(RC)\frac{d}{dt}\left[e^{t/(RC)}Q\right] = \frac{V}{R}e^{t/(RC)} となり、両辺を積分するとある定数 C1C_1 に対して et/(RC)Q=VC et/(RC)+C1e^{t/(RC)}Q = VC\,e^{t/(RC)} + C_1 が得られる。

et/(RC)e^{t/(RC)} で割ると Q(t)=VC+C1e−t/(RC)Q(t) = VC + C_1 e^{-t/(RC)} となる。Q(0)=0Q(0)=0 を適用すると C1=−VCC_1=-VC が得られる。

よって Q(t)=CV(1−e−t/(RC))Q(t)=CV\left(1-e^{-t/(RC)}\right):電荷は 00 から最終値 CVCV へと上昇し、その近づき方は時定数 RCRC に支配される——抵抗やコンデンサが大きいほど遅くなる。

線形一階方程式 dydx+P(x)y=Q(x)\frac{dy}{dx}+P(x)y=Q(x) の積分因子は何か?

ロジスティック成長モデル dPdt=rP(1−PK)\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right) において、定数 KK は何を表すか?

変数分離形方程式 dydx=xy\frac{dy}{dx}=xy を y(0)=1y(0)=1 のもとで解け。

ピカール・リンデレフの定理が y′=f(x,y), y(x0)=y0y'=f(x,y),\ y(x_0)=y_0 に一意な解を与えることを保証する条件はどれか?

参考文献

  1. William E. Boyce, Richard C. DiPrima, Douglas B. Meade (2017). Elementary Differential Equations and Boundary Value Problems
  2. Morris Tenenbaum, Harry Pollard (1985). Ordinary Differential Equations