Bernoulli 方程式
mathdifferential-equationsbernoullifirst-orderlecture
1導入
この講義では、Bernoulli 方程式を適切な置換によって一階線型方程式へ帰着する方法を説明する。
一階線型
y'+p(x)y=q(x)
では、y は 1 次で現れる。Bernoulli 方程式では、右辺に y^n が現れるため、そのままでは線型ではない。しかし y^n の形が規則的なので、u=y^{1-n} という置換で線型化できる。
2標準形
Bernoulli 方程式とは、
y'+p(x)y=q(x)y^n
の形に書ける一階微分方程式である。本講義ではまず n を整数として扱う。n を一般の実数として扱う場合は、y^{1-n} が実数値で定義され微分可能な区間、たとえば y>0 の範囲で考える。
ただし、n=0 と n=1 は特別である。n=0 なら
y'+p(x)y=q(x)
であり、最初から一階線型である。n=1 なら
y'+p(x)y=q(x)y
すなわち
y'+(p(x)-q(x))y=0
であり、これも一階線型である。したがって Bernoulli 置換が本質的に必要なのは n\ne 0,1 の場合である。
3命題:Bernoulli 方程式は一階線型へ帰着する
仮定として、n\ne 0,1 であり、u=y^{1-n} が定義され微分可能な範囲で
y'+p(x)y=q(x)y^n
を考える。たとえば n が整数なら y\ne 0 の範囲でよい。u=y^{1-n} とおくと、u は
u'+(1-n)p(x)u=(1-n)q(x)
を満たす。これは u に関する一階線型微分方程式である。
3.1証明
u=y^{1-n} とおく。合成関数の微分より、
u'=(1-n)y^{-n}y'
である。元の方程式を y^n で割ると、
y^{-n}y'+p(x)y^{1-n}=q(x)
である。ここで y^{1-n}=u、また y^{-n}y'=u'/(1-n) だから、
\frac{1}{1-n}u'+p(x)u=q(x)
である。両辺に 1-n を掛けて、
u'+(1-n)p(x)u=(1-n)q(x)
を得る。
したがって Bernoulli 方程式は、置換によって一階線型へ帰着する。
4解法の手順
- y'+p(x)y=q(x)y^n に整理する。
- n=0 または n=1 でないかを確認する。
- y=0 が元の式を満たすかを先に確認する。
- u=y^{1-n} が定義され微分可能な範囲で u=y^{1-n} と置く。
- u の一階線型方程式を積分因子で解く。
- u=y^{1-n} の実数値で微分可能な逆変換の各枝を、その成立区間で選択し、元の方程式へ代入して検証する。
y^n で割る操作を行うため、y=0 の解を落とさないように注意する。
5例:y'+y=xy^2
y'+y=xy^2
は Bernoulli 方程式であり、n=2、p(x)=1、q(x)=x である。まず y=0 は元の式を満たすので、定数解である。
つぎに y\ne 0 として、u=y^{1-2}=y^{-1} とおく。命題より、
u'+(1-2)\cdot 1\cdot u=(1-2)x
だから
u'-u=-x
である。これは一階線型である。積分因子は
\mu=e^{\int -1\,dx}=e^{-x}
である。両辺に e^{-x} を掛けると、
(e^{-x}u)'=-xe^{-x}
である。積分して
e^{-x}u=(x+1)e^{-x}+C
を得る。したがって
u=x+1+Ce^x
である。u=1/y なので、
y=\frac{1}{x+1+Ce^x}
を得る。これに加えて、最初に確認した y=0 も解である。
この式は、分母 x+1+Ce^x が 0 にならない連結区間ごとに解を定める。分母の零点をまたいで解を延長することはできない。
6Logistic 方程式との関係
r>0, K>0 とする。Logistic 方程式
y'=ry\left(1-\frac{y}{K}\right)
は
y'-ry=-\frac{r}{K}y^2
と整理できるため、Bernoulli 方程式の特殊例である。
ただし Logistic 方程式は自律方程式でもある。そのため、明示解を求める前に、右辺の符号から解の増減と平衡解の近傍での向きを判定できる。
data/lecture/math/differential-equations/logistic-equation.lecture.n.md
7適用範囲
Bernoulli 型は、非線型方程式のうち特別に線型化できる型である。非線型というだけで Bernoulli 型とは限らない。y^n が標準形に従って現れているかを確認する必要がある。
また、置換は y\ne 0 の範囲で行う。y=0 で元の方程式が定義され、かつ y=0 が実際に解となる場合は、その定数解を別に扱う。
8演習リンク
data/exercise/math/differential-equations/classifying-and-solving-first-order-odes.exercise.n.md