一階線型微分方程式と積分因子
mathdifferential-equationsfirst-orderintegrating-factorlecture
1導入
この講義では、積分因子を暗記する公式ではなく、左辺を積の微分に作り替えるための因子として導出する。
一階線型の標準形は
y'+p(x)y=q(x)
である。この形では、y' と y が別々に現れるため、そのまま積分しにくい。そこで、ある関数 \mu(x) を掛けて、
\mu y'+\mu p(x)y
を
(\mu y)'
という 1 個の微分にしたい。
2積分因子を導入する理由
\mu=e^{\int p(x)\,dx} という式は、積の微分公式から系統的に導出できる。
まず、直接積分できる形は
y'=q(x)
である。この形なら、左辺が 1 個の微分なので、そのまま積分できる。
一方、一階線型
y'+p(x)y=q(x)
では、左辺に y' と y が別々に現れる。p(x)y が存在するため、左辺をそのまま 1 個の関数の導関数として積分することはできない。
そこで、積の微分公式
(\mu y)'=\mu y'+\mu' y
を逆向きに適用する。もし元の左辺を
\mu y'+\mu' y
の形に作り替えられれば、左辺は 1 個の微分になり、積分できる。
ところが元の方程式には y'+p(x)y しかない。そこで両辺に未知の関数 \mu(x) を掛けて、
\mu y'+\mu p(x)y
を作る。積の微分と比べると、うまくいくためには
\mu'=\mu p(x)
であればよい。この条件から \mu を決めると、積分因子が自然に出る。
したがって、\mu を掛ける目的は、左辺を 1 個の積の導関数に変換し、直接積分できる形へ帰着することである。
3一階線型とは何か
一階線型微分方程式とは、
y'+p(x)y=q(x)
の形に書ける一階微分方程式である。ここで p(x),q(x) は x の関数であり、y や y^2 を含まない。
線型という用語は、L[y]=y'+p(x)y とおいたとき、L[ay_1+by_2]=aL[y_1]+bL[y_2] が成立することを意味する。この性質は、未知関数 y とその導関数 y' が積や非線型関数を介さずに現れることに対応する。たとえば
y'+xy=1
は一階線型である。一方、
y'+y^2=x
は y^2 を含むので、そのままでは一階線型ではない。
4命題1:積分因子の導出
仮定として、p,q が対象区間で連続であり、
y'+p(x)y=q(x)
を考える。\mu が対象区間で 0 にならず、
\mu'=\mu p(x)
を満たすなら、
(\mu y)'=\mu q(x)
が成立する。また、\mu\ne0 であるため、この変換後の方程式は元の方程式と同値である。
4.1証明
元の方程式の両辺に \mu を掛ける。
\mu y'+\mu p(x)y=\mu q(x)
仮定 \mu'=\mu p(x) により、左辺は
\mu y'+\mu' y
である。積の微分公式
(\mu y)'=\mu' y+\mu y'
より、
(\mu y)'=\mu q(x)
を得る。
逆に、(\mu y)'=\mu q(x) を積の微分公式で展開し、\mu'=\mu p(x) を代入すると、\mu y'+\mu p(x)y=\mu q(x) となる。\mu は対象区間で 0 にならないため、\mu で除して元の方程式を得る。
したがって一階線型を解く要点は、\mu'=\mu p(x) を満たす \mu を決定することである。
5積分因子の形
\mu'=\mu p(x)
は \mu\ne 0 の範囲で
\frac{\mu'}{\mu}=p(x)
と書ける。両辺を積分すると、
\log|\mu|=\int p(x)\,dx
である。したがって、積分因子として
\mu(x)=e^{\int p(x)\,dx}
を選べばよい。積分定数を入れると \mu に定数倍が付くだけであり、最終的な解の任意定数に吸収される。
6命題2:一般解の公式
仮定として、p,q が区間 I で連続であるとする。\mu=e^{\int p(x)\,dx} とおけば、一階線型方程式
y'+p(x)y=q(x)
の解は
y=\frac{\int \mu(x)q(x)\,dx+C}{\mu(x)}
で与えられる。
6.1証明
命題1より、
(\mu y)'=\mu q(x)
である。両辺を積分すると、
\mu y=\int \mu q(x)\,dx+C
である。\mu=e^{\int p(x)\,dx} は 0 にならないので、両辺を \mu で割って
y=\frac{\int \mu(x)q(x)\,dx+C}{\mu(x)}
を得る。
この証明は、一般解の公式が、積分因子によって方程式を直接積分できる形へ変換した結果であることを示す。
7例1:y'+y=x
y'+y=x
では p(x)=1,\ q(x)=x である。積分因子は
\mu=e^{\int 1\,dx}=e^x
である。両辺に e^x を掛けると、
e^x y'+e^x y=xe^x
である。左辺は
(e^x y)'
なので、
(e^x y)'=xe^x
である。積分して
e^x y=\int xe^x\,dx+C=(x-1)e^x+C
を得る。したがって
y=x-1+Ce^{-x}
である。
検算する。y'=1-Ce^{-x} だから、
y'+y=(1-Ce^{-x})+(x-1+Ce^{-x})=x
となり、元の方程式を満たす。
8定数変化法との関係
定数変化法は、一階線型を別の角度から解く方法である。
まず非同次項 q(x) を 0 にした同次方程式
y'+p(x)y=0
を解く。y=0 はこの同次方程式の解である。y\ne 0 の範囲では変数分離でき、
\frac{dy}{y}=-p(x)\,dx
から
y=Ce^{-\int p(x)\,dx}
を得る。C=0 とすれば、先に確認した零解もこの解族に含まれる。
非同次方程式では、この定数 C を関数 C(x) に変える。
y=C(x)e^{-\int p(x)\,dx}
とおく。ここで
\mu=e^{\int p(x)\,dx}
とすれば
\mu y=C(x)
である。積分因子法で得た
(\mu y)'=\mu q(x)
は、
C'(x)=\mu q(x)
と同じ意味である。つまり定数変化法は、同次解の定数を動かして非同次項に対応する方法であり、積分因子法と同じ構造を持つ。
9線型作用素による表現
左辺を
L[y]=y'+p(x)y
と書くと、一階線型方程式は
L[y]=q(x)
である。これは線型作用素の方程式である。同次方程式 L[y]=0 の解は核 \ker L にあたり、非同次方程式の一般解は
\text{特殊解}+\text{同次解}
の形になる。
この作用素表現は、後続の二階線型方程式や連立系でも使用する。
data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md
10適用範囲
この方法は、方程式が
y'+p(x)y=q(x)
という線型の形に整理できる場合に適用する。p,q が区間 I で連続なら、I 内の点で指定した初期値に対し、I 全体で一意な解が存在する。
一方、
y'+p(x)y=q(x)y^n
のように y^n が現れる場合は、そのままでは一階線型ではない。ただし Bernoulli 型では、置換によって一階線型へ帰着できる場合がある。
data/lecture/math/differential-equations/bernoulli-equations.lecture.n.md
積分因子は、次の完全微分方程式では、微分形式を全微分へ変換する因子として一般化される。
data/lecture/math/differential-equations/exact-differential-equations.lecture.n.md
11演習リンク
data/exercise/math/differential-equations/first-order-linear-odes-and-integrating-factors.exercise.n.md