解析学基礎/微分方程式入門
未知の関数とその導関数を含む方程式を微分方程式と呼びます。本来これは総称であって正確には1変数関数の導関数(常微分)を含む常微分方程式と多変数関数の所謂偏導関数(偏微分)を含む偏微分方程式に大別されるのですが本稿では偏微分方程式は扱わないのでここで単に微分方程式と言った場合は常微分方程式を表すものとします。また理論的な詳細は常微分方程式を参照してもらう事とし本稿では解法の説明を重点的に行う事にします。 {\displaystyle }
1階微分方程式

微分方程式の中に含まれる導関数の最高階数をその微分方程式の階数といいます。この節では1階導関数を含む所謂1階微分方程式について解説します。
1階微分方程式

1階微分方程式に於いて最も単純なのは変数分離形と呼ばれる以下の形の微分方程式です。;
1階微分方程式

この微分方程式は両辺をg(y)で割って積分する事により以下のように解く事ができます。;
1階微分方程式

ここで C = C 2 − C 1 {\displaystyle C=C_{2}-C_{1}} という風にして積分定数(微分方程式の世界では任意定数と呼ぶ事が多いです)を1つにまとめました。(以下この様な任意定数の書き換えを適宜(暗黙裡に)行う場合があります。)この任意定数を含む解を一般解といいます。他方任意定数に特定の値を与えて一意的に決定した解を特殊解と呼びます。任意定数は y 0 = f ( x 0 ) {\displaystyle y_{0}=f(x_{0})} のように独立変数がある値のときの従属変数の値を与えることで決定する事がしばしばあります。この y 0 = f ( x 0 ) {\displaystyle y_{0}=f(x_{0})} の様な関係式 を初期条件と言い、関数値 y 0 {\displaystyle y_{0}} を初期値と言います。初期条件から特殊解を決定する計算は一般解を求めるよりは易しいことが多いですので、一般解が求まったことをもって解が求まったとしてもよさそうにも思われますが、物理現象などを記述するために微分方程式を解いているのならば、特殊解こそが現実の事象を記述しているとも言えますので、初期条件から特殊解を求めることも重要です。
1階微分方程式

未知関数及びその導関数についての1次式で書かれる微分方程式を線型微分方程式と呼びます。ここでは1階のそれについて述べる事にします。微分方程式
1階微分方程式

を1階線型微分方程式といいます。この微分方程式の一般解は以下の公式で与えられます。
1階微分方程式

証明 y = e − ∫ P ( x ) d x { ∫ e ∫ P ( x ) d x Q ( x ) d x + C } {\displaystyle y=e^{-\int P(x)dx}\left\{\int e^{\int P(x)dx}Q(x)dx+C\right\}} とすると、
1階微分方程式

導出 まず微分方程式(☆)の右辺を零とおいた y ′ + P ( x ) y = 0 {\displaystyle y'+P(x)y=0} ...(☆☆) という等式を考えます (これを(☆)に付随する同次方程式または斉次方程式といいます)。この同次方程式(☆☆)は変数分離形なので
1階微分方程式

という風に解けます( C 2 = ± e C 1 {\displaystyle C_{2}=\pm e^{C_{1}}} と書き換えました。)。ここでこの解を用いて微分方程式(☆)を満たす関数を作ります。この解の任意定数 C 2 {\displaystyle C_{2}} をある関数 u ( x ) {\displaystyle u(x)} に書き換えた式
1階微分方程式

が(☆)の解になっているとします(この書き換えを定数変化法といいます。任意の関数yは u ( x ) {\displaystyle u(x)} を適切に設定すればこのように表せることに注意)。この(☆☆☆)の両辺を微分すれば
1階微分方程式

が成り立ちます。そしてこの最右辺第2項を移項すれば(☆)より
1階微分方程式

でなければならないことが分かります。従ってこの等式から
1階微分方程式

が得られます。この関数 u ( x ) {\displaystyle u(x)} を(☆☆☆)に代入する事により上述の公式が導かれる事が分かります。(証明終)
1階微分方程式

1階微分方程式
1階微分方程式

をベルヌーイの微分方程式と呼びます。この微分方程式は m = 0 {\displaystyle m=0} のときは1階線型微分方程式、 m = 1 {\displaystyle m=1} のときは変数分離形の微分方程式ですので、既にみたようにして解くことができますが、 m ≠ 0 , 1 {\displaystyle m\neq 0,1} の場合は次のようにして解くことができます。
1階微分方程式

未知関数 y = y ( x ) {\displaystyle y=y(x)} に対して u = u ( x ) = y 1 − m {\displaystyle u=u(x)=y^{1-m}} という関数を考えて、これを微分すると
1階微分方程式

が得られます。これは関数 u = u ( x ) {\displaystyle u=u(x)} に関する1階線型微分方程式に他なりません。
1階微分方程式

1階微分方程式
1階微分方程式

をクレローの微分方程式といいます。この微分方程式に於いて y ′ = p {\displaystyle y'=p} と書いて両辺を微分すると
1階微分方程式

すなわち
1階微分方程式

が成り立ちます。これから p ′ = 0 {\displaystyle p'=0} または x + f ′ ( p ) = 0 {\displaystyle x+f'(p)=0} が得られますが、前者からただちに p = C {\displaystyle p=C} (任意定数)が求まるので、これにより
1階微分方程式

という風に一般解が導かれます。一方の後者の方ですが、この等式と当該微分方程式を連立させて p {\displaystyle p} を消去すると上述の一般解とは一味違う解が導かれます。 この解は上記一般解(✪)の任意定数にいかなる数値を代入しても求める事ができない解であり、特異解と呼ばれています。 (実はこの特異解、直線族(✪)を自身の各点に於ける接線とする曲線(包絡線といいます)を表していたりします。)
2階微分方程式

ここからは2階微分方程式の解法について述べる事にします。
2階微分方程式

以下の微分方程式は置換積分を用いる事により変数分離形の1階微分方程式に帰着させる事ができます。;
2階微分方程式

上式の両辺にdy/dxを掛けて積分すれば
2階微分方程式

となります。ここで u = d y d x {\displaystyle u={\frac {dy}{dx}}} とおけば ∫ d y d x d 2 y d x 2 d x = ∫ u d u {\displaystyle \int {\frac {dy}{dx}}{\frac {d^{2}y}{dx^{2}}}dx=\int udu} ですので、以下のように計算できます。;
2階微分方程式

これは(形は複雑ですが)変数分離形の微分方程式に他なりません。
2階微分方程式

一般の2階の線型微分方程式は、
2階微分方程式

という形で書ける微分方程式です。また微分方程式(*)に付随する同次方程式は次式で表されます。;
2階微分方程式

ここで微分方程式(**)が y 1 = y 1 ( x ) , y 2 = y 2 ( x ) {\displaystyle y_{1}=y_{1}(x),y_{2}=y_{2}(x)} という2つの解を持っていたとします。するとこれらの線型結合
2階微分方程式

も(**)を満たす事が代入により分かります。このような解 y 1 , y 2 {\displaystyle y_{1},y_{2}} を(**)の基本解といいます。一般的な形(関数係数かつ非同次)の2階線型微分方程式(*)を解くのは少し難しいのでそれに関しては後述させてもらう事にしてここでは定数係数かつ同次の2階線型微分方程式の解法について述べてゆく事にします。
2階微分方程式

それでは2階定数係数同次線型微分方程式
2階微分方程式

の解を探してみることにしましょう。やや天下り的ですが y = e λ x {\displaystyle y=e^{\lambda x}} とおくと、 y ′ = λ e λ x , y ′′ = λ 2 e λ x {\displaystyle y'=\lambda e^{\lambda x},y''=\lambda ^{2}e^{\lambda x}} なので、これを方程式に代入すれば
2階微分方程式

となります。指数関数が零になる事はないので上式から二次方程式
2階微分方程式

が得られます(これを特性方程式といいます)。この二次方程式の相違なる二解を λ 1 , λ 2 {\displaystyle \lambda _{1},\lambda _{2}} と書けば二つの指数関数 e λ 1 x , e λ 2 x {\displaystyle e^{\lambda _{1}x},e^{\lambda _{2}x}} が微分方程式(★)の基本解になっている事が分かります。ゆえに(***)より次式で与えられる関数が微分方程式(★)の一般解になっている事が分かりました。;
2階微分方程式

特性方程式(★★)が重解 λ {\displaystyle \lambda } を持つ場合、以下で見るように(★)は y = e λ x {\displaystyle y=e^{\lambda x}} という形以外の解を持ちます。解がある関数 u = u ( x ) {\displaystyle u=u(x)} を用いて y = u e λ x {\displaystyle y=ue^{\lambda x}} と書けるとすると、
2階微分方程式

より
2階微分方程式

が得られます。ここで λ {\displaystyle \lambda } は二次方程式(★★)の解なので、 ( λ 2 + a λ + b ) = 0 {\displaystyle (\lambda ^{2}+a\lambda +b)=0} であり、また解と係数の関係より λ + λ = − a {\displaystyle \lambda +\lambda =-a} すなわち 2 λ + a = 0 {\displaystyle 2\lambda +a=0} が成り立つので、(#)が零となるためには u ′′ = 0 {\displaystyle u''=0} すなわち u = c 1 + c 2 x {\displaystyle u=c_{1}+c_{2}x} ( c 1 , c 2 {\displaystyle c_{1},c_{2}} :積分定数) が成立せねばなりません。従って(★★)が重解を持つとき微分方程式(★)は以下の解を持つ事が分かりました。;
2階微分方程式

特性方程式(★★)の相違なる二解が複素数解 α ± i β {\displaystyle \alpha \pm i\beta } ( i {\displaystyle i} :虚数単位)である場合一般解(♪)は複素数値関数を用いて表される事になりますがこの解は以下のように適当な変形を施す事によって実数値関数を用いた式に書き換えられる事がわかります。 λ 1 = α + i β , λ 2 = α − i β {\displaystyle \lambda _{1}=\alpha +i\beta ,\lambda _{2}=\alpha -i\beta } とおけばオイラーの公式 e i θ = cos θ + i sin θ , ( θ ∈ R ) {\displaystyle e^{i\theta }=\cos {\theta }+i\sin {\theta },(\theta \in \mathbb {R} )} より
2階微分方程式

が成り立ちます。したがって、解(♪)は
2階微分方程式

と書き換えられますので、任意定数を書き換えることで
2階微分方程式

と表されることがわかりました。
2階微分方程式

(♪♪♪)は(♪)とはずいぶん見た目の異なる表現ですので、本当に解になっているか不安になるかもしれません。実際に(★)を満たすか確かめてみましょう。(♪♪♪)を微分すると
2階微分方程式

です。これを(★)の左辺に代入してみます。(★★)の解と係数の関係より a = − 2 α , b = α 2 + β 2 {\displaystyle a=-2\alpha ,b=\alpha ^{2}+\beta ^{2}} であることに注意すると、
2階微分方程式

となり、確かに解であることが確かめられました。
2階微分方程式

(♪♪♪)に於いて c 1 = A sin θ 0 , c 2 = A cos θ 0 {\displaystyle c_{1}=A\sin {\theta _{0}},c_{2}=A\cos {\theta _{0}}} とおけば解は
2階微分方程式

と表せます。(♪♪♪)の表現は解が二次元の線型空間をなすことがわかりやすい表現ですが、(♬)の表現は任意定数を決めたときの関数の様子がよりわかりやすい表現で、どちらも重要な表現です。
2階微分方程式

ここからは定数係数の2階非同次線型微分方程式
2階微分方程式

の解法について考察する事にしましょう。1階線型微分方程式の節で述べた定数変化法がここでも使えます。まず同次方程式(★)の一般解を y = c 1 y 1 + c 2 y 2 {\displaystyle y=c_{1}y_{1}+c_{2}y_{2}} と書く事にします。そしてこの解の任意定数 c 1 , c 2 {\displaystyle c_{1},c_{2}} をそれぞれ二つの関数 u = u ( x ) , v = v ( x ) {\displaystyle u=u(x),v=v(x)} に置き換えたものを
2階微分方程式

で表す事とし、これが(†)を満たすものとします。上式を微分すると
2階微分方程式

が得られます。ここで再び天下り的ですが u ′ y 1 + v ′ y 2 = 0 {\displaystyle u'y_{1}+v'y_{2}=0} を仮定しておく事にします。すると
2階微分方程式

ですが、これを更に微分すれば
2階微分方程式

となります。これらを(†)に代入すれば
2階微分方程式

となりますが、2つの関数 y 1 , y 2 {\displaystyle y_{1},y_{2}} は(★)の基本解だから上式の括弧内は零になります。従って仮定と組み合わせる事により2つの等式
2階微分方程式

が求まります。これを未知数 u ′ , v ′ {\displaystyle u',v'} に関する連立一次方程式と考えて解けば
2階微分方程式

となります。ここで W = W ( y 1 , y 2 ) = y 1 y 2 ′ − y 1 ′ y 2 {\displaystyle W=W(y_{1},y_{2})=y_{1}y'_{2}-y'_{1}y_{2}} とおきました。(この関数Wをロンスキー行列式またはロンスキアンと呼びます。)上述の連立方程式の解を積分すれば次式が得られます。;
2階微分方程式

従ってこれらを y p = u y 1 + v y 2 {\displaystyle y_{p}=uy_{1}+vy_{2}} に代入すれば以下の式が(†)を満たすことがわかります。;
2階微分方程式

この y p {\displaystyle y_{p}} と同次方程式(★)の一般解 y G = c 1 y 1 + c 2 y 2 {\displaystyle y_{G}=c_{1}y_{1}+c_{2}y_{2}} を用いて
2階微分方程式

と表される関数を考えると、
2階微分方程式

ですので、この式が非同次方程式(†)の一般解になっていることがわかります。すなわち、
2階微分方程式

は一般解になっています。
2階微分方程式

(補遺)ここで上述のオイラーの公式について一点補足しておきます。当公式は指数関数及び正弦・余弦のテイラー展開を用いて証明するのが一般的なのですが一応以下のような導出法もあります。;実変数複素数値関数 f ( x ) = g ( x ) + i h ( x ) {\displaystyle f(x)=g(x)+ih(x)} に対してその微積分を
2階微分方程式

で定義しておきます。ここで関数 f ( x ) = cos x + i sin x {\displaystyle f(x)=\cos x+i\sin x} を考えて上式を使って微分すれば
2階微分方程式

であり、すなわち
2階微分方程式

が成り立ちます。これは変数分離形の1階微分方程式であり、 f ( x ) = C e i x {\displaystyle f(x)=Ce^{ix}} という風に解けます。また f ( 0 ) = 1 {\displaystyle f(0)=1} からただちにC=1が求まるので結局
2階微分方程式

が導かれる事が分かります。(略証終)
2階微分方程式

2階冪関数係数線型微分方程式 {\displaystyle }
2階微分方程式

を(2階の)オイラーの微分方程式と呼びます。この線型微分方程式は適当な変数変換を行う事により定数係数のそれに帰着させる事ができます。以下ではその事実を示してゆく事にしましょう。独立変数 x {\displaystyle x} に対して x = e t {\displaystyle x=e^{t}} とおくと、 d t d x = e − t {\displaystyle {\frac {dt}{dx}}=e^{-t}} が得られます。次に関数 y = y ( x ) {\displaystyle y=y(x)} の媒介変数 t {\displaystyle t} に関する微分をドット記号 ❝・❞ を用いて y ̇ = d y d t {\displaystyle {\dot {y}}={\frac {dy}{dt}}} 、 y ̈ = d 2 y d t 2 {\displaystyle {\ddot {y}}={\frac {d^{2}y}{dt^{2}}}} という風に書く事にすれば
2階微分方程式

が成立します。これらの等式を微分方程式 (❤) に代入すれば y ̈ + ( a − 1 ) y ̇ + b y = Q ( e t ) {\displaystyle {\ddot {y}}+(a-1){\dot {y}}+by=Q(e^{t})} ...(♠) が導かれます。これは変数 t {\displaystyle t} を独立変数とする未知関数 y = y ( e t ) {\displaystyle y=y(e^{t})} に関する2階定数係数線型微分方程式に他なりません。
簡単な連立微分方程式

ここでは連立方程式の微分方程式版とも言うべき連立微分方程式の解法について考察する事にします。一般的な解法を述べるのは難しいので以下では極めて単純な1階定数係数同次線型連立微分方程式の解法についてのみ議論してゆく事にします。なお本節では行列の固有値と固有ベクトルに関する知識を仮定しています。
簡単な連立微分方程式

二つの未知関数 y 1 = y 1 ( x ) , y 2 = y 2 ( x ) {\displaystyle y_{1}=y_{1}(x),y_{2}=y_{2}(x)} に対して以下の連立微分方程式が成立しているとします。;
簡単な連立微分方程式

この連立微分方程式は
簡単な連立微分方程式

と表せますので、
簡単な連立微分方程式

とすると
簡単な連立微分方程式

という風に簡潔に表す事ができます。(ベクトル値関数の微分は成分ごとに行うものとします。)
簡単な連立微分方程式

ここで行列Aが異なる2つの固有値 λ 1 , λ 2 {\displaystyle \lambda _{1},\lambda _{2}} を持つとすると、適当な行列Pを用いて
簡単な連立微分方程式

と表せます。ここで z = P − 1 y {\displaystyle \mathbf {z} =P^{-1}\mathbf {y} } とすると、 d y d x = A y {\displaystyle {\frac {d\mathbf {y} }{dx}}=A\mathbf {y} } より
簡単な連立微分方程式

が得られます。
簡単な連立微分方程式

z = ( z 1 z 2 ) {\displaystyle \mathbf {z} ={\begin{pmatrix}z_{1}\\z_{2}\\\end{pmatrix}}} とすると、
簡単な連立微分方程式

となります。この方程式を解くと、
簡単な連立微分方程式

であることはただちにわかります。 P = ( p 11 p 12 p 21 p 22 ) {\displaystyle P={\begin{pmatrix}p_{11}&p_{12}\\p_{21}&p_{22}\\\end{pmatrix}}} とすると y = P z {\displaystyle \mathbf {y} =P\mathbf {z} } より
簡単な連立微分方程式

となります。従って上式に表れる二つのベクトルを
簡単な連立微分方程式

とおく事により以下の連立微分方程式(♡)の解の公式が導かれる事になります。;
簡単な連立微分方程式

前節で述べた2階定数係数同次線型微分方程式
簡単な連立微分方程式

は未知関数及びその導関数を y = y 1 {\displaystyle y=y_{1}} 、 y ′ = y 1 ′ = y 2 {\displaystyle y'=y'_{1}=y_{2}} とおく事により
簡単な連立微分方程式

と連立微分方程式の形で書けます。ゆえに2階線型微分方程式(★)は上述の連立微分方程式(♡)の特別な場合に過ぎない事が分かります。そしてこれは
簡単な連立微分方程式

を用いて d y d x = A † y {\displaystyle {\frac {d\mathbf {y} }{dx}}=A^{\dagger }\mathbf {y} } と書けます。ここでこの行列 A † {\displaystyle A^{\dagger }} の固有多項式を求めてみると
簡単な連立微分方程式

となります(ただし E {\displaystyle E} は単位行列)。即ち前節で天下り的に導入した特性方程式は実は固有方程式 Ψ A † ( λ ) = 0 {\displaystyle \Psi _{A^{\dagger }}(\lambda )=0} に他ならなかったという訳です。