数学の作り方 How to make Mathematics

第1章 微分方程式とは・初等解法

自然は「変化の法則」で書かれている

物理も生物も経済も、法則はたいてい「変化のしかた」で語られます。速度は位置の変化、加速度は速度の変化、 人口の増え方は今の人口に比例、放射性物質の減り方は今の量に比例——。つまり自然が私たちに教えてくれるのは、 求めたい量 y(t)y(t) そのものではなく、その導関数 yy'yy'' との関係です。

この「導関数を含む方程式」が微分方程式で、それを解いて元の量 y(t)y(t) を復元することが、この分野の営みです。 これは積分の壮大な一般化——「変化の法則から未来を計算する」ことにほかなりません。この章では、 まず微分方程式の言葉を定め、その幾何的な意味(方向場:各点で進む向きを指定する矢印の場)を可視化します。 そして最も基本的な2つの解法——変数分離形一階線形——を手に入れます。

微分方程式 = 未知関数とその導関数の関係式。解く = 変化の法則から関数を復元する。方向場が幾何的な姿。

微分方程式の言葉

定義 常微分方程式・階数・解

未知関数 y(x)y(x) とその導関数の関係式 F(x,y,y,,y(n))=0F(x,y,y',\dots,y^{(n)})=0常微分方程式(ODE)という (最高階の導関数が nn 階ならnn)。これを満たす関数をという。任意定数を含む解の族を 一般解、定数を特定した解を特殊解という。

nn 階の一般解には nn 個の任意定数が現れます(nn 回積分するイメージ)。定数を決めるには、初期条件が要ります。

定義 初期値問題

nn 階方程式に、一点での値 y(x0)=y0, y(x0)=y1,,y(n1)(x0)=yn1y(x_0)=y_0,\ y'(x_0)=y_1,\dots,y^{(n-1)}(x_0)=y_{n-1} を課したものを初期値問題(IVP)という。 (次章以降で見るように、良い条件下で解は一意に定まる。)

方向場:微分方程式の幾何

一階方程式 y=f(x,y)y'=f(x,y) は、幾何的にはこう読めます——平面の各点 (x,y)(x,y) で、進むべき傾き f(x,y)f(x,y) を指定している。 矢印(傾き)の場を作り、その向きに沿って進む曲線が解。これが方向場(傾き場)です。

定義 方向場

一階方程式 y=f(x,y)y'=f(x,y) に対し、各点 (x,y)(x,y) に傾き f(x,y)f(x,y) の短い線分を対応させた図を方向場(傾き場)という。 解曲線は、各点で方向場に接する曲線。

下で体感してください。方程式を選ぶと方向場(灰の線分)が描かれ、キャンバスをクリックすると、その点を通る解曲線が 引かれます。解は「各点の傾きに素直に従って進む道」——微分方程式は道順を指定する場、解はその道を辿った軌跡です。 初期条件(出発点)を決めれば、道が一本に決まる様子が見えます。

キャンバスをクリックで解曲線を追加

方向場は、解を式で求めなくても振る舞いを読める強力な視点です(第7章の定性理論の芽)。 では、実際に式で解く方法へ進みます。

変数分離形

最も基本的な解法。yy の部分と xx の部分に「分けて」両辺を積分します。

定理 変数分離形

dydx=g(x)h(y)\dfrac{dy}{dx}=g(x)h(y) の形(変数分離形)は、h(y)0h(y)\ne0 の範囲で dyh(y)=g(x)dx+C\int\frac{dy}{h(y)}=\int g(x)\,dx+C を計算して解ける。h(y)=0h(y_*)=0 なる yy_* は定数解(平衡解)。

証明

1h(y)dydx=g(x)\frac1{h(y)}\frac{dy}{dx}=g(x) の両辺を xx で積分。左辺は置換積分(微積分第7章)で dyh(y)\int\frac{dy}{h(y)} に等しい。∎

変数分離の例

  • 放射性崩壊・指数増殖 y=kyy'=kydyy=kdx\int\frac{dy}y=\int k\,dx より lny=kx+C\ln|y|=kx+Cy=Cekxy=Ce^{kx}
  • ロジスティック方程式 y=y(1y)y'=y(1-y)dyy(1y)=dx\int\frac{dy}{y(1-y)}=\int dx、部分分数分解して lny1y=x+C\ln\frac{y}{1-y}=x+Cy=11+Cexy=\frac1{1+Ce^{-x}}y=0y=0(不安定)と y=1y=1(安定)が平衡解——上の方向場で確かめられる。
  • 同心円 y=x/yy'=-x/yydy=xdx\int y\,dy=-\int x\,dxy22=x22+C\frac{y^2}2=-\frac{x^2}2+Cx2+y2=x^2+y^2= 一定。

一階線形微分方程式

yy について 11 次の方程式 y+p(x)y=q(x)y'+p(x)y=q(x) は、積分因子という魔法の掛け算で解けます。

定理 一階線形と積分因子

y+p(x)y=q(x)y'+p(x)y=q(x) は、積分因子 μ(x)=ep(x)dx\mu(x)=e^{\int p(x)\,dx} を両辺に掛けると左辺が (μy)(\mu y)' にまとまり、 (μy)=μqy=1μ(x)(μ(x)q(x)dx+C).(\mu y)'=\mu q\quad\Rightarrow\quad y=\frac1{\mu(x)}\Big(\int\mu(x)q(x)\,dx+C\Big).

証明

μ=ep\mu=e^{\int p}μ=pμ\mu'=p\mu を満たす。μy+μpy=μq\mu y'+\mu p y=\mu q の左辺は積の微分則より μy+μy=(μy)\mu y'+\mu' y=(\mu y)'。 両辺を積分して μy=μqdx+C\mu y=\int\mu q\,dx+C。∎

積分因子 μ=ep\mu=e^{\int p} の役割は「左辺 y+pyy'+py を、完全な微分 (μy)(\mu y)' に化けさせる」こと。この一手で、 あとは積分するだけ。なぜ epe^{\int p} なのかは、μ=pμ\mu'=p\muμ\mu 自身が y=pyy'=py の解)を要求すれば自然に出ます。

一階線形の例

yy=xy'-y=xp=1,q=xp=-1,q=x):μ=ex\mu=e^{-x}(exy)=xex(e^{-x}y)'=xe^{-x}、積分して exy=(x+1)ex+Ce^{-x}y=-(x+1)e^{-x}+Cy=(x+1)+Cexy=-(x+1)+Ce^x。解は特殊解 (x+1)-(x+1) と同次解 CexCe^x の和——この「特殊解+同次解」の構造は 線形方程式一般で本質的になる。

つまずきポイント

注意 よくある誤解

  • 一般解は任意定数を含む。 nn 階なら nn 個。初期条件で定数を決めて特殊解を得る。定数を忘れると解の族の一部しか得られない。
  • 変数分離で h(y)=0h(y)=0 の定数解を見落とさない。 y=y(1y)y'=y(1-y)y=0,y=1y=0,y=1dyh\int\frac{dy}{h} の外にある平衡解。割り算で消える。
  • 積分因子は epe^{\int p} y+py=qy'+py=qppyy の係数、yy' の係数を 11 に正規化してから)。符号・正規化に注意。
  • 方向場は「解を求めずに振る舞いを読む」道具。 式で解けなくても、平衡・漸近・安定性が見える。

この章のまとめ

  • 微分方程式は「未知関数とその導関数の関係式」。一般解は任意定数を含み、初期条件で特殊解に定まる。一階 y=f(x,y)y'=f(x,y) の幾何は方向場(各点の傾き)で、解は方向場に沿う曲線。
  • 変数分離形 y=g(x)h(y)y'=g(x)h(y)dyh=gdx\int\frac{dy}{h}=\int g\,dxh=0h=0 の定数解に注意。指数増殖・ロジスティック・同心円。
  • 一階線形 y+py=qy'+py=q積分因子 μ=ep\mu=e^{\int p} を掛けて (μy)=μq(\mu y)'=\mu q、積分して解く。解は「特殊解+同次解」。

次章は、初等解法をさらに広げます——完全形・積分因子、ベルヌーイ・リッカチ方程式を扱います。