← 戻る ライブラリ › 微分方程式と力学系 › 常微分方程式 微分方程式と力学系
一階微分方程式 関数とその一階導関数を関係づける方程式で、変数分離法や積分因子法によって解かれる。
直観 今、それはどれくらいの速さで変化しているのか? 多くの現象は、量そのものの式によってではなく、その量がどれくらいの速さで変化しているかという規則によって最も自然に記述される:個体群はその時点での大きさに比例して増加し、熱いコーヒーは室温よりどれだけ熱いかに比例した速さで冷め、充電中のコンデンサの電圧は満充電に近づくほどゆっくり上昇する。一階微分方程式とはまさにこのような規則であり——未知関数 y ( x ) y(x) y ( x ) をその導関数 y ′ ( x ) y'(x) y ′ ( x ) 自身と関係づける——それを解くとは、変化率の規則から関数全体を復元することを意味する。
1階常微分方程式 y ′ = y − c x 2 y' = y - c x^2 y ′ = y − c x 2 の方向場と、異なる初期値 y ( t 0 ) = y 0 y(t_0) = y_0 y ( t 0 ) = y 0 から出発する3本の解曲線。 大学 変数分離形と線形方程式 定義: 変数分離形方程式
一階方程式が変数分離形であるとは、d y d x = f ( x ) g ( y ) \frac{dy}{dx}=f(x)g(y) d x d y = f ( x ) g ( y ) の形に書けることをいう:右辺が x x x のみに依存する部分と y y y のみに依存する部分の積に分解される。y y y を含む項をすべて片側に、x x x を含む項をすべてもう片側に移すと、方程式は2つの通常の積分 ∫ d y g ( y ) = ∫ f ( x ) d x \int \frac{dy}{g(y)}=\int f(x)\,dx ∫ g ( y ) d y = ∫ f ( x ) d x に帰着する。
d y d x = f ( x ) g ( y ) \frac{dy}{dx}=f(x)g(y) d x d y = f ( x ) g ( y ) 例えば、自励系方程式 d y d x = g ( y ) \frac{dy}{dx}=g(y) d x d y = g ( y ) (変化率が x x x ではなく現在の値 y y y のみに依存する場合)は f ( x ) = 1 f(x)=1 f ( x ) = 1 の特別な場合である。その平衡解はちょうど g ( y ) = 0 g(y)=0 g ( y ) = 0 となる箇所に現れ、これらの定数解は不用意に g ( y ) g(y) g ( y ) で割ると見落としやすい。
d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) 一階方程式の種類 種類 標準形 解法 変数分離形 d y d x = f ( x ) g ( y ) \frac{dy}{dx}=f(x)g(y) d x d y = f ( x ) g ( y ) 両辺を積分する:∫ d y g ( y ) = ∫ f ( x ) d x \int \frac{dy}{g(y)}=\int f(x)\,dx ∫ g ( y ) d y = ∫ f ( x ) d x 線形一階 d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) μ ( x ) = e ∫ P ( x ) d x \mu(x)=e^{\int P(x)\,dx} μ ( x ) = e ∫ P ( x ) d x を掛けてから積分する自励系 d y d x = g ( y ) \frac{dy}{dx}=g(y) d x d y = g ( y ) 変数分離形の特別な場合;平衡点は g ( y ) = 0 g(y)=0 g ( y ) = 0
大学 解くこと、そして解の存在を保証すること P ( x ) P(x) P ( x ) と Q ( x ) Q(x) Q ( x ) がある区間で連続であり、μ ( x ) = e ∫ P ( x ) d x \mu(x)=e^{\int P(x)\,dx} μ ( x ) = e ∫ P ( x ) d x とすると、線形方程式 d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) のすべての解は y ( x ) = 1 μ ( x ) ( ∫ μ ( x ) Q ( x ) d x + C ) y(x) = \frac{1}{\mu(x)}\left(\int \mu(x) Q(x)\,dx + C\right) y ( x ) = μ ( x ) 1 ( ∫ μ ( x ) Q ( x ) d x + C ) で与えられる。
なぜ正しいのか? この技巧は平方完成に似ている:方程式全体に、うまく選んだ因子 μ ( x ) \mu(x) μ ( x ) を掛けることで、左辺が単一の積 μ ( x ) y ( x ) \mu(x)y(x) μ ( x ) y ( x ) の導関数にまとまり、微分方程式を直接積分できるものに変える。
証明 μ ( x ) = e ∫ P ( x ) d x \mu(x)=e^{\int P(x)\,dx} μ ( x ) = e ∫ P ( x ) d x とおく;連鎖律で微分すると、指数 ∫ P ( x ) d x \int P(x)\,dx ∫ P ( x ) d x の導関数が P ( x ) P(x) P ( x ) であることから μ ′ ( x ) = P ( x ) μ ( x ) \mu'(x)=P(x)\mu(x) μ ′ ( x ) = P ( x ) μ ( x ) が得られる。
d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) の両辺に μ ( x ) \mu(x) μ ( x ) を掛けると、μ ( x ) y ′ + μ ( x ) P ( x ) y = μ ( x ) Q ( x ) \mu(x)y' + \mu(x)P(x)y = \mu(x)Q(x) μ ( x ) y ′ + μ ( x ) P ( x ) y = μ ( x ) Q ( x ) となる。
積の法則により d d x [ μ ( x ) y ( x ) ] = μ ( x ) y ′ ( x ) + μ ′ ( x ) y ( x ) \frac{d}{dx}[\mu(x)y(x)] = \mu(x)y'(x) + \mu'(x)y(x) d x d [ μ ( x ) y ( x )] = μ ( x ) y ′ ( x ) + μ ′ ( x ) y ( x ) であり、μ ′ ( x ) = P ( x ) μ ( x ) \mu'(x)=P(x)\mu(x) μ ′ ( x ) = P ( x ) μ ( x ) なので、これは上式の左辺 μ ( x ) y ′ + μ ( x ) P ( x ) y \mu(x)y' + \mu(x)P(x)y μ ( x ) y ′ + μ ( x ) P ( x ) y にちょうど等しい。したがって方程式は d d x [ μ ( x ) y ( x ) ] = μ ( x ) Q ( x ) \frac{d}{dx}[\mu(x)y(x)] = \mu(x)Q(x) d x d [ μ ( x ) y ( x )] = μ ( x ) Q ( x ) となる。
両辺を x x x について積分すると、任意定数 C C C に対して μ ( x ) y ( x ) = ∫ μ ( x ) Q ( x ) d x + C \mu(x)y(x) = \int \mu(x)Q(x)\,dx + C μ ( x ) y ( x ) = ∫ μ ( x ) Q ( x ) d x + C が得られる。μ ( x ) \mu(x) μ ( x ) (指数関数なので決してゼロにならない)で両辺を割ると y ( x ) = 1 μ ( x ) ( ∫ μ ( x ) Q ( x ) d x + C ) y(x) = \frac{1}{\mu(x)}\left(\int \mu(x) Q(x)\,dx + C\right) y ( x ) = μ ( x ) 1 ( ∫ μ ( x ) Q ( x ) d x + C ) が得られる。
f ( x , y ) f(x,y) f ( x , y ) が ( x 0 , y 0 ) (x_0,y_0) ( x 0 , y 0 ) の周りの長方形上で連続であり、y y y に関してリプシッツ条件 ∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ L ∣ y 1 − y 2 ∣ |f(x,y_1)-f(x,y_2)|\le L|y_1-y_2| ∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ L ∣ y 1 − y 2 ∣ (ある定数 L L L について)を満たすとする。このとき初期値問題 y ′ = f ( x , y ) , y ( x 0 ) = y 0 y'=f(x,y),\ y(x_0)=y_0 y ′ = f ( x , y ) , y ( x 0 ) = y 0 は、x 0 x_0 x 0 を含むある区間上で一意な解を持つ。
なぜ正しいのか? リプシッツ条件は f f f が y y y 方向にどれだけ急峻に変化しうるかを制限し、これにより近くの解の軌道が分裂したり交差したりするのを防ぐ。これはまさに、各点をちょうど一つの軌道が通る、明確に定義された流れを保証するために必要な条件である。
証明 初期値問題を積分方程式に書き換えると、連続関数 y ( x ) y(x) y ( x ) が y ′ = f ( x , y ) , y ( x 0 ) = y 0 y'=f(x,y),\ y(x_0)=y_0 y ′ = f ( x , y ) , y ( x 0 ) = y 0 を解くことは y ( x ) = y 0 + ∫ x 0 x f ( t , y ( t ) ) d t y(x) = y_0 + \int_{x_0}^x f(t,y(t))\,dt y ( x ) = y 0 + ∫ x 0 x f ( t , y ( t )) d t を解くことと同値である。この同値性は微積分学の基本定理から従う。
連続関数に対するピカール反復作用素 T T T を T [ y ] ( x ) = y 0 + ∫ x 0 x f ( t , y ( t ) ) d t T[y](x) = y_0 + \int_{x_0}^x f(t,y(t))\,dt T [ y ] ( x ) = y 0 + ∫ x 0 x f ( t , y ( t )) d t で定義し、仮定が成り立つ長方形内に f f f がとどまるように選んだ小さな区間 I = [ x 0 − h , x 0 + h ] I=[x_0-h,x_0+h] I = [ x 0 − h , x 0 + h ] に制限する。
I I I 上の2つの連続関数 y 1 , y 2 y_1,y_2 y 1 , y 2 に対して、リプシッツ条件より ∣ T [ y 1 ] ( x ) − T [ y 2 ] ( x ) ∣ ≤ ∫ x 0 x ∣ f ( t , y 1 ( t ) ) − f ( t , y 2 ( t ) ) ∣ d t ≤ L h sup t ∈ I ∣ y 1 ( t ) − y 2 ( 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)| ∣ T [ y 1 ] ( x ) − T [ y 2 ] ( x ) ∣ ≤ ∫ x 0 x ∣ f ( t , y 1 ( t )) − f ( t , y 2 ( t )) ∣ d t ≤ L h sup t ∈ I ∣ y 1 ( t ) − y 2 ( t ) ∣ となる。L h < 1 Lh<1 L h < 1 となるほど h h h を小さく選べば、T T T は上限ノルムにおける縮小写像となる。
バナッハの不動点定理により、完備距離空間(上限ノルムを備えた I I I 上の連続関数の空間)上の縮小写像はちょうど一つの不動点 y y y を持ち、y = T [ y ] y=T[y] y = T [ y ] を満たす。この不動点こそが、I I I 上での元の初期値問題の一意な解である。
大学 実世界での応用と具体例 一階方程式は、限られた資源に制限された個体群成長(生態学・疫学におけるロジスティック成長)、電子工学における充放電するRC回路、熱力学におけるニュートンの冷却法則、物理学における放射性崩壊や炭素年代測定、金融における時間変化する利率での連続複利などをモデル化する。以下の2つの具体例は、最も一般的な2つのパターンを最後まで解く。
例: ロジスティック人口成長
個体群 P ( t ) P(t) P ( t ) がロジスティックモデル d P d t = r P ( 1 − P K ) \frac{dP}{dt}=rP\left(1-\frac{P}{K}\right) d t d P = r P ( 1 − K P ) に従って成長し、初期個体数 P ( 0 ) = P 0 P(0)=P_0 P ( 0 ) = P 0 であるとする。P ( t ) P(t) P ( t ) を求めよ。
解答 この方程式は変数分離形である:d P P ( 1 − P / K ) = r d t \frac{dP}{P(1-P/K)}=r\,dt P ( 1 − P / K ) d P = r d t と書き直せる。
部分分数分解 1 P ( 1 − P / K ) = 1 P + 1 / K 1 − P / K \frac{1}{P(1-P/K)}=\frac1P+\frac{1/K}{1-P/K} P ( 1 − P / K ) 1 = P 1 + 1 − P / K 1/ K により左辺を項ごとに積分でき、ln ∣ P ∣ − ln ∣ 1 − P / K ∣ = r t + C 1 \ln|P|-\ln|1-P/K|=rt+C_1 ln ∣ P ∣ − ln ∣1 − P / K ∣ = r t + C 1 、すなわち定数をまとめて K K K を掛けた後 ln ∣ P K − P ∣ = r t + C 2 \ln\left|\frac{P}{K-P}\right|=rt+C_2 ln K − P P = r t + C 2 が得られる。
指数をとると、ある定数 A A A に対して P K − P = A e r t \frac{P}{K-P}=Ae^{rt} K − P P = A e r t となる。t = 0 t=0 t = 0 で P ( 0 ) = P 0 P(0)=P_0 P ( 0 ) = P 0 を用いると A = P 0 K − P 0 A=\frac{P_0}{K-P_0} A = K − P 0 P 0 が得られる。
P P P について解いて整理すると閉じた形 P ( t ) = K P 0 e r t K + P 0 ( e r t − 1 ) P(t) = \frac{K P_0 e^{rt}}{K + P_0(e^{rt}-1)} P ( t ) = K + P 0 ( e r t − 1 ) K P 0 e r t が得られ、これは P 0 P_0 P 0 付近から始まり t → ∞ t\to\infty t → ∞ で環境収容力 K K K に近づく、実際の個体群データで見られるまさにS字曲線である。
例: RC回路の充電
定電圧源 V V V から抵抗を通して充電されるコンデンサは d Q d t + Q R C = V R \frac{dQ}{dt} + \frac{Q}{RC} = \frac{V}{R} d t d Q + R C Q = R V (Q ( 0 ) = 0 Q(0)=0 Q ( 0 ) = 0 )を満たす。ここで Q ( t ) Q(t) Q ( t ) は時刻 t t t における電荷である。Q ( t ) Q(t) Q ( t ) を求めよ。
解答 これは P ( t ) = 1 R C P(t)=\frac{1}{RC} P ( t ) = R C 1 、右辺項 V R \frac{V}{R} R V を持つ線形方程式であり、積分因子は μ ( t ) = e t / ( R C ) \mu(t)=e^{t/(RC)} μ ( t ) = e t / ( R C ) である。
両辺に掛けると d d t [ e t / ( R C ) Q ] = V R e t / ( R C ) \frac{d}{dt}\left[e^{t/(RC)}Q\right] = \frac{V}{R}e^{t/(RC)} d t d [ e t / ( R C ) Q ] = R V e t / ( R C ) となり、両辺を積分するとある定数 C 1 C_1 C 1 に対して e t / ( R C ) Q = V C e t / ( R C ) + C 1 e^{t/(RC)}Q = VC\,e^{t/(RC)} + C_1 e t / ( R C ) Q = V C e t / ( R C ) + C 1 が得られる。
e t / ( R C ) e^{t/(RC)} e t / ( R C ) で割ると Q ( t ) = V C + C 1 e − t / ( R C ) Q(t) = VC + C_1 e^{-t/(RC)} Q ( t ) = V C + C 1 e − t / ( R C ) となる。Q ( 0 ) = 0 Q(0)=0 Q ( 0 ) = 0 を適用すると C 1 = − V C C_1=-VC C 1 = − V C が得られる。
よって Q ( t ) = C V ( 1 − e − t / ( R C ) ) Q(t)=CV\left(1-e^{-t/(RC)}\right) Q ( t ) = C V ( 1 − e − t / ( R C ) ) :電荷は 0 0 0 から最終値 C V CV C V へと上昇し、その近づき方は時定数 R C RC R C に支配される——抵抗やコンデンサが大きいほど遅くなる。
よくある誤り. 変数分離をする際、g ( y ) g(y) g ( y ) で割ると g ( y ) = 0 g(y)=0 g ( y ) = 0 となる定数解が知らぬ間に失われる——ロジスティック方程式では、P ≡ 0 P\equiv0 P ≡ 0 と P ≡ K P\equiv K P ≡ K は完全に有効な平衡解であるが、一般公式 P ( t ) = K P 0 e r t K + P 0 ( e r t − 1 ) P(t) = \frac{K P_0 e^{rt}}{K + P_0(e^{rt}-1)} P ( t ) = K + P 0 ( e r t − 1 ) K P 0 e r t は有限の t t t に対してこれらを決して生成しないため、別途確認して明記しなければならない。もう一つのよくある誤りは、まだ標準形 d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) になっていない方程式に積分因子の公式を適用してしまうことである——y ′ y' y ′ の係数が 1 1 1 でない場合、先に両辺を割らなければ P ( x ) P(x) P ( x ) 、ひいては μ ( x ) \mu(x) μ ( x ) が誤ったものになる。 歴史的ノート
アイザック・ニュートンの1671年の流率法に関する論文は、既に級数展開によって単純な微分方程式を解いており、変化率を研究の主要な対象として扱っていた。独立に、ゴットフリート・ヴィルヘルム・ライプニッツの1690年代の論文は、変数分離の技法と、今日でも使われている簡潔な記法 d y / d x dy/dx d y / d x を導入した。流率法と微分という二つのアプローチが合わさって、現代の常微分方程式の講義が今なお最初に教える体系的な分類と解法(変数分離形、線形、完全形)を生み出した。
アイザック・ニュートン ゴットフリート・ヴィルヘルム・ライプニッツ
線形一階方程式 d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx}+P(x)y=Q(x) d x d y + P ( x ) y = Q ( x ) の積分因子は何か?
μ ( x ) = e ∫ P ( x ) d x \mu(x)=e^{\int P(x)\,dx} μ ( x ) = e ∫ P ( x ) d x μ ( x ) = e − ∫ P ( x ) d x \mu(x)=e^{-\int P(x)\,dx} μ ( x ) = e − ∫ P ( x ) d x μ ( x ) = ∫ P ( x ) d x \mu(x)=\int P(x)\,dx μ ( x ) = ∫ P ( x ) d x μ ( x ) = P ( x ) \mu(x)=P(x) μ ( x ) = P ( x ) ロジスティック成長モデル d P d t = r P ( 1 − P K ) \frac{dP}{dt}=rP\left(1-\frac{P}{K}\right) d t d P = r P ( 1 − K P ) において、定数 K K K は何を表すか?
初期個体数 内的成長率 環境収容力 平衡に達するまでの時間
変数分離形方程式 d y d x = x y \frac{dy}{dx}=xy d x d y = x y を y ( 0 ) = 1 y(0)=1 y ( 0 ) = 1 のもとで解け。
y = e x 2 / 2 y=e^{x^2/2} y = e x 2 /2 y = e x 2 y=e^{x^2} y = e x 2 y = 1 + x 2 2 y=1+\frac{x^2}{2} y = 1 + 2 x 2 y = x 2 2 y=\frac{x^2}{2} y = 2 x 2 ピカール・リンデレフの定理が y ′ = f ( x , y ) , y ( x 0 ) = y 0 y'=f(x,y),\ y(x_0)=y_0 y ′ = f ( x , y ) , y ( x 0 ) = y 0 に一意な解を与えることを保証する条件はどれか?
f f f が連続であるだけf f f が ( x 0 , y 0 ) (x_0,y_0) ( x 0 , y 0 ) の近くで y y y に関してリプシッツ連続であることf f f が x x x について線形であることf f f がいたるところ有界であること