logo

リッカチ微分方程式の解 📂微分方程式

リッカチ微分方程式の解

定義

下記の1階非線形微分方程式をリカッチ方程式Ricatti equationという。

$$ y^\prime = P(x)y+Q(x)y^2+R(x) $$

説明

$y_{1}$を既に知ってる特別解particular solutionだとすると、一般解は$y=y_{1}+u(x)$の形で表される。この時$u(x)$は任意の定数で、$n=2$の時のベルヌーイ微分方程式を解いて得られる。

解答

リカッチ方程式は見た目があまりにも複雑で解くのが難しい。だから、簡単なトリックを使って、私たちが解きやすい形に変えてから解かないといけない。

  • ステップ 1.

    $y_{1}$を任意の特別解とする。この特別解を得る方法は特に決まっておらず、直感や洞察、または何度もの試みで得る必要がある。$y$に$2$を入れればいいのかな?$3x$を入れればいいのかな?という感じで得る。そして、一般解を$y=y_{1}+u(x)$の形と仮定しよう。$u(x)$は任意の定数だ。

  • ステップ 2.

    与えられた微分方程式に$y=y_{1}+u$を代入すると

    $$ \begin{align*} && y_{1}^\prime + u^\prime = P(y_{1}+u) + Q(y_{1}+u)^2 +R \\ \implies && y_{1}^\prime + u^\prime = P(y_{1}+u) + Q(y_{1}^2 +2y_{1}u +u^2) +R \end{align*} $$

    右辺を整理すると

    $$ y_{1}^\prime + u^\prime = \left( Py_{1} + Qy_{1}^2 + R \right) +\left( P + 2Qy_{1} \right) u +Qu^2 $$

  • ステップ 3.

    この時$y_{1}$が与えられた微分方程式の解であるため$y_{1}^\prime=Py_{1}+Qy_{1}^2+R$を満たす。従って、左右の共通項を消去すると以下のようになる。

    $$ u^\prime = \left( P + 2Qy_{1} \right) u +Qu^2 $$

    これはベルヌーイ微分方程式で$n=2$の時と同じ形だ。だから、ベルヌーイ微分方程式の解法で$u(x)$を求めることができる。そうすれば、最終的に与えられた微分方程式の一般解$y=y_{1}+u(x)$が得られる。

例題

微分方程式$y^\prime = 2- 2xy+y^2$を解け。

$y$に$2x$を代入すると成り立つため、$y_{1}=2x$と置くことができる。そうすると、一般解は

$$ y=y_{1}+u(x)=2x+u \\ y^\prime=2+u^\prime $$

与えられた微分方程式に代入すると

$$ \begin{align*} && 2+u^\prime &= 2 – 2x(2x+u)+(2x+u)^2 \\ \implies && 2+u^\prime &= 2-4x^2-2xu+4x^2+4xu+u^2 \\ \implies && u^\prime &= 2xu+u^2 \\ \implies && u^\prime –2xu &= u^2 \end{align*} $$

つまり、ベルヌーイ方程式で$n=2$の場合だ。

ベルヌーイ方程式

$$ u^\prime + (1-n)pu=q(1-n) $$

両辺を$u^2$で割ると

$$ u^{-2}u^\prime – 2xu^{-1}=1 $$

ここで、$w \equiv u^{1-n}=u^{-1}$と置換してベルヌーイ方程式を解くと$\dfrac{dw}{dx}=-u^{-2}\dfrac{du}{dx}$になるので

$$ \begin{align*} && -w^\prime –2xw &= 1 \\ \implies && w^\prime + 2xw &= -1 \\ \implies && w &= e^{-x^2} \left[ -\displaystyle \int e^{x^2}dx+C \right] \end{align*} $$

しかし、上で$w =u^{-1}$と置換したので

$$ u=\dfrac{e^{x^2}}{C- \displaystyle \int e^{x^2} dx} $$

従って、最終的に一般解は次のようになる。

$$ y=y_{1}+u(x)=2x+\dfrac{e^{x^2}}{C- \displaystyle \int e^{x^2}dx} $$