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
Bernoulli Equations
1Introduction
This lecture explains how an appropriate substitution reduces a Bernoulli equation to a first-order linear equation.
In the first-order linear equation
y'+p(x)y=q(x),
the unknown y occurs to the first power. A Bernoulli equation also contains a term y^n and is therefore generally nonlinear. Its power-law structure nevertheless permits linearization through the substitution u=y^{1-n}.
2Standard Form
A Bernoulli equation is a first-order differential equation that can be written as
y'+p(x)y=q(x)y^n.
This lecture first treats integer values of n. For a general real exponent n, the analysis must be restricted to an interval on which y^{1-n} is real-valued and differentiable; the region y>0 is a standard choice.
The cases n=0 and n=1 are exceptional. If n=0, the equation is
y'+p(x)y=q(x),
which is already first-order linear. If n=1, then
y'+p(x)y=q(x)y,
or equivalently
y'+(p(x)-q(x))y=0,
which is also first-order linear. Thus the Bernoulli substitution is essential only when n\ne0,1.
3Proposition: Reduction of a Bernoulli Equation to First-Order Linear Form
Assume that n\ne0,1 and consider
y'+p(x)y=q(x)y^n
on a region where u=y^{1-n} is defined and differentiable. For integer n, it is sufficient to work where y\ne0. Under the substitution u=y^{1-n}, the function u satisfies
u'+(1-n)p(x)u=(1-n)q(x).
This is a first-order linear differential equation in u.
3.1Proof
Set u=y^{1-n}. The chain rule gives
u'=(1-n)y^{-n}y'.
Dividing the original equation by y^n gives
y^{-n}y'+p(x)y^{1-n}=q(x).
Since y^{1-n}=u and y^{-n}y'=u'/(1-n), this equation becomes
\frac{1}{1-n}u'+p(x)u=q(x).
Multiplication by 1-n yields
u'+(1-n)p(x)u=(1-n)q(x).
Therefore, the substitution reduces the Bernoulli equation to a first-order linear equation.
4Solution Procedure
- Rewrite the equation as y'+p(x)y=q(x)y^n.
- Check whether n=0 or n=1.
- Before division, determine whether y=0 satisfies the original equation.
- On a region where u=y^{1-n} is defined and differentiable, set u=y^{1-n}.
- Solve the resulting first-order linear equation in u by an integrating factor.
- On each interval of validity, select every real-valued differentiable branch of the inverse transformation from u=y^{1-n} to y, and verify it in the original equation.
Because the derivation divides by y^n, the zero solution must be checked beforehand so that it is not discarded.
5Example: y'+y=xy^2
The equation
y'+y=xy^2
is a Bernoulli equation with n=2, p(x)=1, and q(x)=x. The function y=0 satisfies the original equation and is therefore a constant solution.
For y\ne0, set u=y^{1-2}=y^{-1}. The proposition gives
u'+(1-2)u=(1-2)x,
and hence
u'-u=-x.
This is first-order linear. Its integrating factor is
\mu=e^{\int -1\,dx}=e^{-x}.
Multiplication by e^{-x} gives
(e^{-x}u)'=-xe^{-x}.
Integration yields
e^{-x}u=(x+1)e^{-x}+C,
so
u=x+1+Ce^x.
Since u=1/y,
y=\frac{1}{x+1+Ce^x}.
Together with this family, the previously identified function y=0 is also a solution.
For each value of C, the displayed formula defines a solution on every connected interval on which the denominator x+1+Ce^x is nonzero. The solution cannot be extended across a zero of the denominator.
6Connection with the Logistic Equation
Assume r>0 and K>0. The logistic equation
y'=ry\left(1-\frac{y}{K}\right)
can be rewritten as
y'-ry=-\frac{r}{K}y^2,
so it is a special Bernoulli equation.
The logistic equation is also autonomous. Before deriving an explicit solution, the sign of its right-hand side can therefore be used to determine monotonicity of solutions and their direction near equilibrium solutions.
data/lecture/math/differential-equations/logistic-equation.lecture.n.md
7Scope
Bernoulli equations form a special class of nonlinear equations that can be linearized. A differential equation is not of Bernoulli type merely because it is nonlinear; the power y^n must occur in the standard form specified above.
Moreover, the substitution is performed where y\ne0. If the original equation is defined at y=0 and y=0 is actually a solution, that constant solution must be treated separately.
8Exercises
data/exercise/math/differential-equations/classifying-and-solving-first-order-odes.exercise.n.md