微分方程式

解析学において、微分方程式（びぶんほうていしき、（英: differential equation）とは、未知関数とその導関数の関係式として書かれている関数方程式である^[1]。

数学の応用分野においてしばしば、異なる2つの変数の関係を調べることが行われる。2変数を対応付ける関数があらわになっていなくても、その導関数（の満たすべき方程式）を適当な仮定の下で定めることができ、そこから目的とする関数を探し出すことができる。

物理法則を記述する基礎方程式は、多くが時間微分、空間微分を含む微分方程式であり、物理学からの要請もあり微分方程式の解法には多くの関心が注がれてきた。

方程式論は解析学の中心的な分野で、フーリエ変換、ラプラス変換等は元々、微分方程式を解くために開発された手法である。また物理学における微分方程式の主要な問題は境界値問題、固有値問題である^[1]。

微分方程式は大きく線型微分方程式と非線型微分方程式に分類される。線形微分方程式の例として、例えばシュレーディンガー方程式が挙げられる。シュレーディンガー方程式は、量子系の状態の時間発展を記述する方法の一つとして広く用いられている。非線型微分方程式の例として、例えばナビエ–ストークス方程式（NS方程式）が挙げられる。NS方程式は流体の運動を記述する基本方程式であり、物理学の応用としても重要な方程式である。しかし、NS方程式の解の存在性は未解決問題でありミレニアム懸賞問題にも選ばれている。

→その他、有名な微分方程式についてはCategory:微分方程式を参照

概要

微分方程式は方程式に含まれる導関数の階数^{[注釈 1]}によって分類され、最も高い階数が n 次である場合、その微分方程式を n 階微分方程式^{[注釈 2]}と呼ぶ^[1]。

いずれの場合も未知関数は一つとは限らず、また、連立する複数の微分方程式を同時に満たす関数を解とするような連立方程式の形を取る場合もある^[1]。これは連立 n 階微分方程式などと呼ばれる。

常微分方程式と偏微分方程式

→詳細は「常微分方程式」および「偏微分方程式」を参照

一変数関数の導関数の関係式で書かれる常微分方程式と多変数関数の偏導関数を含む関係式で書かれる偏微分方程式に分かれる^[1]。

常微分方程式とは例えば、

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\!f(x)-f(x)=0}$

や、

${\displaystyle {\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}\!f(x)-2{\frac {\mathrm {d} }{\mathrm {d} x}}\!f(x)+f(x)=\sin(x)}$

のような方程式である。

また、偏微分方程式は、

${\displaystyle \left(x{\frac {\partial }{\partial y}}-y{\frac {\partial }{\partial x}}\right)f(x,y)=0}$

や、

${\displaystyle \left({\frac {\partial ^{2}}{\partial x^{2}}}+{\frac {\partial ^{2}}{\partial y^{2}}}\right)f(x,y)=\alpha x+\beta y}$

のような格好をした方程式である。

代数的微分方程式

未知関数とその導関数の関係式が、未知関数や導関数を変数と見たときに解析関数を係数とする多項式である場合、代数的微分方程式と呼ばれる。

線形微分方程式

→詳細は「線形微分方程式」を参照

方程式が未知関数の一次式として書けるような方程式を線形微分方程式と呼ぶ。また、線型でない微分方程式は非線形微分方程式^{[注釈 3]}と呼ばれる。 例えば、g(x) を f(x) を含まない既知の関数とすれば、

${\displaystyle \left({\frac {\mathrm {d} }{\mathrm {d} x}}+\alpha \right)f(x)=g(x)}$

は線形微分方程式であり、

${\displaystyle \left({\frac {\mathrm {d} }{\mathrm {d} x}}f(x){\frac {\mathrm {d} }{\mathrm {d} x}}\right)f(x)=g(x)}$

は非線型形微分方程式である。線形と呼ばれる理由は後述する線形斉次な方程式について、解の線形結合がその方程式の一般解をなすためである。

未知関数が 1 つの場合、高階の線型微分方程式を一階線形微分方程式の形に書き直すことができる。 たとえば、{g_k} を既知関数の組として、以下の線形微分方程式が与えられたとき、

${\displaystyle \left(\sum _{k=0}^{n}g_{k+1}(x){\frac {\mathrm {d} ^{k}}{{\mathrm {d} x}^{k}}}\right)f(x)=g_{0}(x),\quad g_{n+1}(x)=1}$

未知関数 f(x) のk 階の導関数を y_k(x) として (k = 0,..., n − 1)、以下の一組の微分方程式を得る。

${\displaystyle {\begin{cases}\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}y_{k-1}(x)=y_{k}(x),\quad k=1,\dots ,n-1\\\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}y_{n-1}(x)=g_{0}(x)-\sum _{k=0}^{n-1}g_{k+1}(x)y_{k}(x)\end{cases}}}$

この微分方程式は、より一般的に、ベクトルと行列の記法を用いて

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\mathbf {y} (x)=A(x)\mathbf {y} (x)+\mathbf {b} (x)}$

と書くことができる。ここで y は未知関数 y₀,..., y_n−1 を成分に持つベクトル、A は既知関数 {a_ij}_{i,j=0,...,n−1} を成分に持つ n × n の行列、b は既知関数 b₀,..., b_n−1 を成分に持つベクトルである。

斉次方程式と非斉次方程式

すべての項が未知関数を含むか 0 であるような線形微分方程式を線形斉次微分方程式^{[注釈 4]}と呼び、斉次でない線形微分方程式は線型非斉次微分方程式^{[注釈 5]}と呼ばれる。同じ意味の言葉として斉次方程式をしばしば同次方程式と呼ぶことがある。 例えば、

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\!f(x)+f(x)=0}$

は斉次な方程式であり、右辺に α を加えた、

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}\!f(x)+f(x)=\alpha }$

は非斉次な方程式である。

より一般の線形常微分方程式について、

${\displaystyle \left(\sum _{k=0}^{n}g_{k+1}(x){\frac {\mathrm {d} ^{k}}{{\mathrm {d} x}^{k}}}\right)f(x)=g_{0}(x),\quad g_{n+1}(x)=1}$

右辺の関数 g₀(x) がゼロならこの方程式は斉次である。 斉次方程式の特徴として、方程式の解 s(x) が得られたとき、その定数倍 cs(x) も方程式の解となる。また、斉次方程式の解の線形結合もその斉次方程式の解になる。

また、非斉次な方程式の解 s_in(x) が得られたとき、元の方程式を斉次な形にしたときの解 s_hom(x) を用いて、非斉次方程式の新たな解 s_in(x) + s_hom(x) を作ることができる。実際、

${\displaystyle \left(\sum _{k=0}^{n}g_{k+1}(x){\frac {\mathrm {d} ^{k}}{\mathrm {d} x^{k}}}\right)\left(s_{\mathrm {in} }(x)+s_{\mathrm {hom} }(x)\right)=g_{0}(x)}$

としたとき、s_in(x), s_hom(x) はそれぞれ

${\displaystyle {\begin{aligned}\left(\sum _{k=0}^{n}g_{k+1}(x){\frac {\mathrm {d} ^{k}}{{\mathrm {d} x}^{k}}}\right)s_{\mathrm {hom} }(x)&=0\\\left(\sum _{k=0}^{n}g_{k+1}(x){\frac {\mathrm {d} ^{k}}{{\mathrm {d} x}^{k}}}\right)s_{\mathrm {in} }(x)&=g_{0}(x)\end{aligned}}}$

を満たすので、s_in(x) + s_hom(x) は元の方程式の解になっている。

確率微分方程式

→詳細は「確率微分方程式」を参照

方程式に含まれる既知関数が確率変数によって記述されるような微分方程式を確率微分方程式^{[注釈 6]}と呼ぶ。確率常微分方程式や確率偏微分方程式はしばしば英語の頭文字を取って“SODE”, “SPDE”と略記される。代表的な例は物理学におけるランジュバン方程式や金融工学におけるブラック-ショールズ方程式がある。確率微分方程式の既知関数は、自身の期待値や相関関数によって特徴付けられる。

時間遅れをもつ微分方程式

常微分方程式　${\displaystyle x'(t)=-ax(t-r)}$（ただしｒは正の実数）のように、時間遅れをもつ微分方程式が、数理生物学のモデルなどに現われる^[2]。

解法

微分方程式に限らず一般の方程式は必ずしも厳密解が得られるとは限らない。従って多く場合は摂動などの手法を用いて近似的な評価を与えるか、ルンゲ＝クッタ法やSOR法、有限要素法のような数値解法によって具体的な解を得ることになる。しかしながらいくつかの基本的な微分方程式については、厳密解が得られたり、形式的に解を書き表せる。

微分方程式の具体的な解法としては代表的なものに、斉次方程式の解を利用して解く定数変化法、グリーン関数を用いた解法、差分方程式を用いた解法、ラプラス変換や逆ラプラス変換を用いた解法などが知られている。

指数関数と微分方程式

一階の線型斉次常微分方程式の中で最も基本的な方程式として次のものがある。

${\displaystyle {\frac {\mathrm {d} y(x)}{\mathrm {d} x}}=y(x).}$

一般の線形微分方程式を解く際も、まずこの種の斉次微分方程式に帰着させるため、この方程式は微分方程式の解法を調べる上で基本的な役割を果たす。 この方程式の解はよく知られているように指数関数となる^{[注釈 7]}。

${\displaystyle y(x)=C\mathrm {e} ^{x}.}$

ここで C は任意定数である。解法は脚注にて紹介する^{[注釈 10]}。

指数関数の有用な性質として、微分作用素を別の定数や関数に置き換えられることが挙げられる。係数が定数の斉次方程式

${\displaystyle \left(\sum _{k=0}^{n}c_{k+1}{\frac {\mathrm {d} ^{k}}{{\mathrm {d} x}^{k}}}\right)f(x)=0,\quad c_{n+1}=1}$

の解として指数関数で書けるものを探すと、f(x) = Cexp(λx) と置き換えて、

${\displaystyle \sum _{k=0}^{n}c_{k+1}\lambda ^{k}=0}$

と書くことができる^{[注釈 11]}。これは λ に対する n 次の代数方程式になっている。 重根がなければ方程式の解が n 個求まることになり、斉次方程式の一般解はそれらの線型結合として表される。

この形の方程式の一般解を求める方法としては定数変化法がある^{[注釈 12]}。

一階線型常微分方程式

一つの未知関数に対する、一般の一階線型常微分方程式は、既知関数を P(x)、Q(x) として、次のように書かれる。

${\displaystyle {\frac {\mathrm {d} y(x)}{\mathrm {d} x}}+P(x)y(x)=Q(x).}$

この一階線型常微分方程式は、一般解が求積法で解ける。 まず、斉次方程式

${\displaystyle {\frac {\mathrm {d} y(x)}{\mathrm {d} x}}+P(x)y(x)=0}$

の一般解は、積分定数を A ≠ 0 として、

${\displaystyle y(x)=A\exp \left(-\int P(x')\,\mathrm {d} x'\right)}$

となる。一階線型常微分方程式の一般解は、斉次方程式の解を利用し A を x の関数とみなす定数変化法によって求められる。

${\displaystyle y(x)=\left\{\int Q(x')\exp \left(\int P(x'')\,\mathrm {d} x''\right)\mathrm {d} x'+C\right\}\exp \left(-\int P(x')\,\mathrm {d} x'\right).}$

ここで C ≠ 0 は積分定数である。

二階線型常微分方程式

二階線型常微分方程式の一般形は、既知関数を P(x), Q(x), R(x) として、次のように書かれる。

${\displaystyle {\frac {\mathrm {d} ^{2}y(x)}{\mathrm {d} x^{2}}}+P(x){\frac {\mathrm {d} y(x)}{\mathrm {d} x}}+Q(x)y=R(x).}$

この二階線型常微分方程式は、このままの形では求積法を用いて一般解を表示することはできない。 もし、右辺を 0 とした斉次方程式の特殊解として、y = y₁ が存在すれば、

${\displaystyle {\frac {\mathrm {d} ^{2}y_{1}(x)}{\mathrm {d} x^{2}}}+P(x){\frac {\mathrm {d} y_{1}(x)}{\mathrm {d} x}}+Q(x)y_{1}(x)=0}$

が成り立つので、z なる未知関数を導入して、

${\displaystyle y(x)=y_{1}(x)z(x)}$

とすれば、二階線型常微分方程式が、z に関する常微分方程式、

${\displaystyle {\begin{aligned}&y_{1}{\frac {\mathrm {d} ^{2}z(x)}{\mathrm {d} x^{2}}}+{\Bigl (}2y'_{1}(x)+P(x)y_{1}(x){\Bigr )}{\frac {\mathrm {d} z(x)}{\mathrm {d} x}}=R(x),\\&y'_{1}(x)={\frac {\mathrm {d} y_{1}(x)}{\mathrm {d} x}},\end{aligned}}}$

に変換される。この常微分方程式は、導関数 dz/dx に関して一階線型常微分方程式なので、求積法で解ける。その一般解を

${\displaystyle z=\psi (C_{1},C_{2},x)~~}$

とすると、二階線型常微分方程式の一般解は、

${\displaystyle y=y_{1}\psi (C_{1},C_{2},x)~~}$

で与えられる。なお、C₁, C₂ は積分定数である。 x の既知関数を含む二階線型常微分方程式で、求積法で解ける微分方程式は少ないが、 次の微分方程式などが知られている^[3][4][5]。

求積法で解ける方程式の例^{[注釈 13]}

方程式	一般解[3]
${\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}-xP(x){\frac {\mathrm {d} y}{\mathrm {d} x}}+P(x)y=0}$	${\displaystyle y=x{\Bigl \{}C_{1}+C_{2}\int {\frac {1}{\,x^{2}\,}}\exp {\Bigl (}\int \!xP(x)\,dx{\Bigr )}\,dx{\Bigr \}}}$
${\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+P(x){\frac {\mathrm {d} y}{\mathrm {d} x}}-a(a+P(x))y=0}$	${\displaystyle y=e^{ax}{\Bigl \{}C_{1}+C_{2}\!\int \exp {\Bigl (}\!-2ax-\!\!\int \!P(x)\,dx{\Bigr )}\,dx{\Bigr \}}}$
${\displaystyle P(x){\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+(a+bx){\frac {\mathrm {d} y}{\mathrm {d} x}}-by=0}$	${\displaystyle y=C_{1}\!\!\int \!\!\int \!\!{\frac {1}{\,P(x)\,}}\exp {\Bigl (}\!-\!\!\int \!{\frac {\,a+bx\,}{P(x)}}\,dx{\Bigr )}\,dx\,dx+C_{2}{\Bigl (}x+{\frac {a}{\,b\,}}{\Bigr )}}$
${\displaystyle {\frac {\mathrm {d} ^{2}y\,}{\mathrm {d} x^{2}}}-\left({\frac {1}{2P(x)}}\cdot {\frac {\mathrm {d} P(x)}{\mathrm {d} x}}\right){\frac {\mathrm {d} y}{\mathrm {d} x}}+P(x)y=0}$	${\displaystyle y=C_{1}\sin \left(\int {\sqrt {P(x)}}\,\mathrm {d} x\right)+C_{2}\cos \left(\int {\sqrt {P(x)}}\,\mathrm {d} x\right)}$
${\displaystyle {\frac {\mathrm {d} ^{2}y\,}{\mathrm {d} x^{2}}}-\left({\frac {1}{P(x)}}\cdot {\frac {\mathrm {d} P(x)}{\mathrm {d} x}}\right){\frac {\mathrm {d} y}{\mathrm {d} x}}-\left(P(x)\right)^{2}y=0}$	${\displaystyle y=C_{1}\exp \left(\int P(x)\,\mathrm {d} x\right)+C_{2}\exp \left(-\int P(x)\,\mathrm {d} x\right)}$
${\displaystyle x{\frac {\mathrm {d} ^{2}y\,}{\mathrm {d} x^{2}}}+(\alpha +\beta x){\frac {\mathrm {d} y}{\mathrm {d} x}}+\beta y=0.}$	${\displaystyle y=x^{1-\alpha }e^{-\beta x}\left(C_{1}\int {}x^{\alpha -2}e^{\beta x}\,\mathrm {d} x+C_{2}\right)}$