メインコンテンツへスキップ
見出し画像

中学でも分かる⑩「微分方程式」

     ここでは微分方程式について解説します。難しい微分方程式を解けるようになるのが目的ではなく、微分積分などの基本をマスターしたレベルで、微分方程式とはどういうものかを理解するのが本記事の目標です。
     読むにあたっては、先に書いた
     中学でも分かる①「微分法」
     中学でも分かる➁「積分法」
    の知識を当たり前のように使うので、分からない場合は最初に読んでおくことをおすすめします。
     また、本シリーズの過去の内容も多く出てくるので、その都度戻るか、最初から本シリーズを読んでおくことをお勧めします。前回までの記事が次の記事で出てくるような構成になっています。

    微分方程式とは

     次のような式を考えます。

    dydx=2x ⋯(∗1)\begin{align*} \dfrac{dy}{dx}=2x \cdots (*1) \end{align*}

     まず、dydx\dfrac{dy}{dx} とは yy を xx で微分するという意味です(中学でも分かる①「微分法」)。
     つまり、この式の意味するところは、
     「なにか y=f(x)y=f(x) という関数があったときに、その関数 yy を xx で微分すると 2x2x になること」
     すなわち、
     「yy の導関数が 2x2x になること」
    を意味しています。 
     するとこの (∗1)(*1) 式から、さかのぼって微分する前の関数  y=f(x)y=f(x) を求めることができます。
     ではどのように求めるのか?
     そのためには、(∗1)(*1) 式の両辺を xx 軸方向に積分します(中学でも分かる➁「積分法」)。

    ∫dydx dx=∫2x dx\begin{align*} \int\dfrac{dy}{dx}\,dx=\int2x\,dx \end{align*}

     ∫□ dx\int\square  dx は xx 軸方向に積分するという記号です。
     ここで、左辺と右辺を別々に見ていきましょう。

    (左辺)について

     左辺は、yy を xx で微分したもの dydx\dfrac{dy}{dx} を xx 軸方向に積分しているので、元に戻って yy となります。

    ddx→ 微分するydydx← 積分する∫□ dx\begin{align*} \\[0pt] &\hspace{24pt}\footnotesize{\dfrac{d}{dx}}\\[4pt] &\hspace{8pt}\xrightarrow[]{  微分する}\\[-18pt] &y\hspace{47pt}\dfrac{dy}{dx}\\[-15pt] &\hspace{8pt}\xleftarrow[  積分する]{}\\[4pt] &\hspace{17pt}\footnotesize{\int□  dx} \end{align*}

     厳密には、積分定数 C1C_1(任意の実数値の定数)が付きます。

    ddx→ 微分するy+C1dydx← 積分する∫□ dx\begin{align*} \\[0pt] &\hspace{47pt}\footnotesize{\dfrac{d}{dx}}\\[4pt] &\hspace{31pt}\xrightarrow[]{  微分する}\\[-18pt] &y+C_1\hspace{47pt}\dfrac{dy}{dx}\\[-15pt] &\hspace{31pt}\xleftarrow[  積分する]{}\\[4pt] &\hspace{39pt}\footnotesize{\int□  dx} \end{align*}

     下の「←\leftarrow」の操作を行うので

    (左辺)=∫dydx dx=y+C1 ⋯(∗2)\begin{align*} (左辺)&=\int\dfrac{dy}{dx}\,dx\\ &=y+C_1 \cdots (*2) \end{align*}

     なお、この計算はあたかも分数の計算のように、形式的に dxdx で約分することが可能です。

    (左辺)=∫dydx dx (dx で約分)=∫dy=∫1 dy=y+C1 ⋯(∗2)\begin{align*} (左辺)&=\int\dfrac{dy}{\cancel{dx}}\,\cancel{dx} \small{(dx  で約分)}\\ &=\int dy\\ &=\int 1  dy\\ &=y+C_1 \cdots (*2) \end{align*}

     最後の ∫dy\int dy は「11 を yy 軸方向に積分する」という意味です。どの方向に積分するかは、dydy の yy という文字に着目します。参考までに、次の違いに注意してください(CC は積分定数)。

     ∫dx=x+C\int dx=x+C(xx 軸方向の積分)
     ∫dy=y+C\int dy=y+C(yy 軸方向の積分)((∗2)(*2) の計算)
     ∫dz=z+C\int dz=z+C(zz 軸方向の積分)

     この(左辺)の計算では、最初は xx 軸方向の積分 ∫□ dx\int□\,dx でしたが、形式的に dxdx で約分することによって、yy 軸方向の積分 ∫□ dy\int□\,dy に切り替わっていることに着目しましょう。

     では次に

    ∫dydx dx=∫2x dx\begin{align*} \int\dfrac{dy}{dx}\,dx=\int2x\,dx \end{align*}

    の右辺を計算します。

    (右辺)について

     右辺は、2x2x を積分をして

    (右辺)=∫2x dx=2∫x2 dx=2⋅12x2+C2=x2+C2 ⋯(∗3)\begin{align*} (右辺)&=\int2x\,dx\\ &=2\int x^2\,dx\\ &=\cancel{2}\cdot\dfrac{1}{\cancel{2}}x^2+C_2\\ &=x^2+C_2 \cdots (*3) \end{align*}

     C2C_2 は積分定数(実数値の定数)です。
     すると (∗2)=(∗3)(*2)=(*3) より

    y+C1=x2+C2\begin{align*} y+C_1=x^2+C_2 \end{align*}

     C1C_1 を移項して

    y=x2+C2−C1\begin{align*} y=x^2+C_2-C_1 \end{align*}

     定数部分をまとめて C2−C1=CC_2-C_1=C とおくと

    y=x2+C\begin{align*} y=x^2+C \end{align*}

     以上により

    dydx=2x ⋯(∗1)\begin{align*} \dfrac{dy}{dx}=2x \cdots (*1) \end{align*}

    を満たす関数 yy は

    y=x2+C\begin{align*} y=x^2+C \end{align*}

    となります。積分定数は1つの文字 CC でまとめておきます。最後に CC を使いために、途中の計算では C1, C2C_1,  C_2 を使っています。

      この (∗1)(*1) のように、未知の関数の導関数を含む等式を微分方程式といいます。
     一般に、ある微分方程式を満たす関数を、その微分方程式の解といい、その方程式のすべての解を求めることを、その微分方程式を解くといいます。
     先の微分方程式

    dydx=2x ⋯(∗1)\begin{align*} \dfrac{dy}{dx}=2x \cdots (*1) \end{align*}

    を解くと、その解は、CC は任意の定数として

    y=x2+C\begin{align*} y=x^2+C \end{align*}

    となります。

    一般解とは

     微分方程式を解くと、任意の定数を含んだ解が得られます。このような任意の定数を含む解を一般解といいます。
     定数部分を、一般的に CC などの文字を用いて表した解が一般解です。

     問題としては次のような形で出題されるので、次のように簡潔に解けばよいでしょう。

    問題 次の微分方程式を解きなさい

    dydx=2x ⋯(∗1)\begin{align*} \dfrac{dy}{dx}=2x \cdots (*1) \end{align*}

    解答

    y=∫2x dx=2⋅12x2+C=x2+C\begin{align*} y&=\int 2x  dx\\ &=\cancel{2}\cdot\dfrac{1}{\cancel{2}}x^2+C\\ &=x^2+C \end{align*}

    より、CC は任意の定数として

    y=x2+C ⋯(答)\begin{align*} y=x^2+C \cdots\small{(答)} \end{align*}

    ■\blacksquare

     結局、(∗1)(*1) の左辺は積分するとダイレクトに yy になるので、右辺を積分するだけで解を求めることができます。

    検算してみると・・・

     ここで、一般解 y=x2+Cy=x^2+C を微分して、それが (∗1)(*1) 式

    dydx=2x ⋯(∗1)\begin{align*} \dfrac{dy}{dx}=2x \cdots (*1) \end{align*}

    を満たすかどうか検算してみましょう。すると

    dydx=(x2+C)′=2x\begin{align*} \dfrac{dy}{dx}=(x^2+C)'=2x \end{align*}

    より、確かに (∗1)(*1) になります。CC は任意の定数なので、どのような定数(任意の実数値)をとってきても (∗1)(*1) は成り立ちます。
     なお、CC は定数なので、微分すると 00 です。微分の復習にもなるので、最初のうちは検算してみるのもよいでしょう。

    特殊解とは

     一方、任意の定数 CC に特定の値を与えて得られる解を特殊解といいます。
     例えば、微分方程式 (∗1)(*1) が
     x=1x=1 のとき y=3y=3
    という条件を満たしているとします。すると、(∗1)(*1) の一般解

    y=x2+C\begin{align*} y=x^2+C \end{align*}

    に x=1, y=3x=1,  y=3 を代入して

    3=12+C3=1+C1+C=3C=3−1=2\begin{align*} &3=1^2+C\\ &3=1+C\\ &1+C=3\\ &C=3-1=2 \end{align*}

    より、定数 CC の値が求まります。この C=2C=2 を、一般解 y=x2+Cy=x^2+C に代入することによって特殊解

    y=x2+2\begin{align*} y=x^2+2 \end{align*}

    が得られます。x=1x=1 のとき y=3y=3 のように、特別な条件を満たしているときの解が特殊解です。

    直接積分形

     先ほどの

    dydx=2x\begin{align*} \dfrac{dy}{dx}=2x \end{align*}

    のように、直接積分すれば一般解が求まる微分方程式を直接積分形といいます。最も基本的な微分方程式です。
     なお、xx で微分していることが明白なときは、この方程式を単に

    y′=2x\begin{align*} y'=2x \end{align*}

    と書く場合があります。「′'」は微分するという記号です。本来 y′y' だけだと見かけ上どの変数で微分しているのか不明ですが、それが自明(普通は xx で微分)であるときは y′y' はよく使われます。同じ微分方程式なので、解き方は先ほどやった解法とまったく同じです。

     では、直接積分形の微分方程式をいくつか解いてみましょう。

    (例題1)基本形

     次の微分方程式

    dydx=4x3−x−2\begin{align*} \dfrac{dy}{dx}=4x^3-x-2 \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=2x=2 のとき y=5y=5 を満たすときの特殊解を求めなさい。

    指針
     左辺を xx 軸方向に積分するとダイレクトに yy になるので、右辺を xx 軸方向に積分するだけで yy を求めることができます。

    解答
    (1) 両辺を xx 軸方向に積分すると

    y=∫(4x3−x−2)dx=4⋅14x4−12x2−2x+C=x4−12x2−2x+C\begin{align*} y&=\int(4x^3-x-2)dx\\ &=\cancel{4}\cdot\dfrac{1}{\cancel{4}}x^4-\dfrac{1}{2}x^2-2x+C\\ &=x^4-\dfrac{1}{2}x^2-2x+C\\ \end{align*}

     よって一般解は、CC は任意の定数として

    y=x4−12x2−2x+C ⋯(答)\begin{align*} y=x^4-\dfrac{1}{2}x^2-2x+C \cdots\small{(答)} \end{align*}

    ■\blacksquare

     では、検算してみましょう(微分の復習にもなります)。
     一般解を微分すると

    dydx=(x4−12x2−2x+C)′=4x3−12⋅2x−2=4x3−x−2\begin{align*} \dfrac{dy}{dx}&=(x^4-\dfrac{1}{2}x^2-2x+C)'\\ &=4x^3-\dfrac{1}{\cancel{2}}\cdot\cancel{2}x-2\\ &=4x^3-x-2 \end{align*}

    となり、問題の微分方程式と同じになります。

    (2) 一般解に x=2, y=5x=2,  y=5 を代入して

    5=24−12⋅22−2⋅2+C5=16−2−4+C5=10+C10+C=5C=5−10=−5\begin{align*} &5=2^4-\dfrac{1}{2}\cdot2^2-2\cdot2+C\\ &5=16-2-4+C\\ &5=10+C\\ &10+C=5\\ &C=5-10=-5 \end{align*}

     これを一般解に代入して

    y=x4−12x2−2x−5 (答)\begin{align*} y=x^4-\dfrac{1}{2}x^2-2x-5 \small{(答)} \end{align*}

    ■\blacksquare

    (例題2)y' の表記

     次の微分方程式

    y′=ex\begin{align*} y'=e^x \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=0x=0 のとき y=2y=2 を満たすときの特殊解を求めなさい。

    指針
     xx で微分していることが明白なときは、yy を xx で微分していると考えて

    y′=dydx\begin{align*} y'=\dfrac{dy}{dx} \end{align*}

    と考えます。つまり、この微分方程式は

    dydx=ex\begin{align*} \dfrac{dy}{dx}=e^x \end{align*}

    と同じ意味です。xx での微分が明白なときは y′y' はよく使われます。
     y′y' を xx 軸方向に積分すると yy になります。

    ddx→ 微分するyy′← 積分する∫□ dx\begin{align*} \\[0pt] &\hspace{24pt}\footnotesize{\dfrac{d}{dx}}\\[4pt] &\hspace{8pt}\xrightarrow[]{  微分する}\\[-13pt] &y\hspace{47pt}y'\\[-12pt] &\hspace{8pt}\xleftarrow[  積分する]{}\\[4pt] &\hspace{17pt}\footnotesize{\int□  dx} \end{align*}

     よって(例題1)と同じように、最初に両辺を xx 軸方向に積分しましょう。左辺は ←\leftarrow の操作で直接 yy になります。

    解答
    (1) 両辺を xx 軸方向に積分すると

    y=∫exdx=ex+C\begin{align*} y&=\int e^xdx\\ &=e^x+C\\ \end{align*}

     よって一般解は、CC は任意の定数として

    y=ex+C ⋯(答)\begin{align*} y=e^x+C \cdots\small{(答)} \end{align*}

    ■\blacksquare

     exe^x は積分しても(微分しても)変わらないことに注意してください。

     明らかですが、検算してみましょう。
     一般解を微分すると

    dydx=(ex+C)′=ex\begin{align*} \dfrac{dy}{dx}&=(e^x+C)'\\ &=e^x \end{align*}

    となり、問題の微分方程式と同じになります。

    (2) 一般解に x=0, y=2x=0,  y=2 を代入して

    2=e0+C2=1+C1+C=2C=2−1=1\begin{align*} &2=e^0+C\\ &2=1+C\\ &1+C=2\\ &C=2-1=1 \end{align*}

     これを一般解に代入して

    y=ex+1 (答)\begin{align*} y=e^x+1 \small{(答)} \end{align*}

    ■\blacksquare

     e0=1e^0=1 であることに注意しましょう。00 乗は 11 と定義されます(中学でも分かる➃「指数関数」)。

    (例題3)両辺を x で割る

     次の微分方程式

    xy′=3\begin{align*} xy'=3 \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=ex=e のとき y=1y=1 を満たすときの特殊解を求めなさい。

    指針
     y′y' は dydx\dfrac{dy}{dx} の意味なので、この微分方程式は

    xdydx=3\begin{align*} x\dfrac{dy}{dx}=3 \end{align*}

    とも表されます。
     xx が左辺にあるので、両辺を xx で割って右辺に持っていくのが最初の一手です。

    解答
    (1) xy′=3xy'=3 の両辺を xx で割って

    y′=3x\begin{align*} y'=\dfrac{3}{x} \end{align*}

     両辺を xx 軸方向に積分して

    y=∫3x dx=3∫1x dx=3log⁡∣x∣+C\begin{align*} y&=\int\dfrac{3}{x}  dx\\ &=3\int\dfrac{1}{x}  dx\\ &=3\log |x|+C \end{align*}

     よって一般解は、CC は任意の定数として

    y=3log⁡∣x∣+C ⋯(答)\begin{align*} y=3\log |x|+C \cdots\small{(答)} \end{align*}

    ■\blacksquare

     途中で積分公式

    ∫1x dx=log⁡∣x∣+C\begin{align*} \int\dfrac{1}{x}\,dx=\log |x|+C \end{align*}

    を用いています。対数の底はネイピア数 ee です。ee が底のときは

    log⁡e∣x∣=log⁡∣x∣\begin{align*} \log_e |x|=\log |x| \end{align*}

    と省略されることに注意しましょう(中学でも分かる➅「ネイピア数 e」)。

     では、検算してみましょう。
     (log⁡∣x∣)′=1x(\log |x|)'=\dfrac{1}{x}
    の微分公式を使います。

    xy′=x(3log⁡∣x∣+C)′=x⋅3⋅1x=3\begin{align*} xy'&=x(3\log |x|+C)'\\ &=\cancel{x}\cdot3\cdot\dfrac{1}{\cancel{x}}\\ &=3 \end{align*}

    となり、問題の微分方程式と同じになります。

    (2) 一般解に x=e, y=1x=e,  y=1 を代入すると

    1=3log⁡∣e∣+C1=3⋅1+C1=3+C3+C=1C=1−3=−2\begin{align*} &1=3\log |e|+C\\ &1=3\cdot1+C\\ &1=3+C\\ &3+C=1\\ &C=1-3=-2 \end{align*}

     これを一般解に代入して

    y=3log⁡∣x∣−2 (答)\begin{align*} y=3\log |x|-2 \small{(答)} \end{align*}

    ■\blacksquare

     log⁡∣e∣=log⁡ee=1\log |e|=\log_e e=1 の計算にも注意しましょう(中学でも分かる➄「対数関数」)。

    変数分離形

     次は変数分離形とよばれる微分方程式です。

    (例題4)右辺に y があるとき

     次の微分方程式

    dydx=2xy\begin{align*} \dfrac{dy}{dx}=2xy \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=0x=0 のとき y=3y=3 を満たすときの特殊解を求めなさい。

    指針
     左辺に xx だけはなく、yy があるときの解法です。
     まずは、yy は xx の関数であるので、y=f(x)y=f(x) となっていることを意識しましょう。

    dydx=2xf(x)\begin{align*} \dfrac{dy}{dx}=2xf(x) \end{align*}

    と同じ意味です。実は右辺は xx だけの式(関数)になっています。この関数 y=f(x)y=f(x) の f(x)f(x) をこれから求めるわけです。
     このときは、最初に両辺を yy で割って、yy を左辺にもっていきます。
     なお、両辺を yy で割る前に
    (case1) y≠0y\ne0 のとき
    (case2) y=0y=0 のとき
    で場合分けをしましょう。00 で割ることはできないので、場合分けが必要です。

    解答
    (1)
    (case1) y≠0y\ne0 のとき
     両辺を y (≠0)y  (\ne0) で割って

    1y dydx=2x\begin{align*} \dfrac{1}{y}  \dfrac{dy}{dx}=2x \end{align*}

     両辺を xx 軸方向に積分して

    ∫1y dydx dx=∫2x dx\begin{align*} \int \dfrac{1}{y}  \dfrac{dy}{\cancel{dx}}  \cancel{dx}&=\int 2x  dx \end{align*}

     形式的に約分すると dydy だけが残って

    ∫1y dy=2∫x dx\begin{align*} \int \dfrac{1}{y}  dy=2\int x  dx \end{align*}

     左辺が yy 軸方向の積分、右辺が xx 軸方向の積分と両辺で分かれることがポイントです。ここで、両辺を積分して

    log⁡∣y∣=2⋅12x2+C1\begin{align*} \log |y|=\cancel{2}\cdot\dfrac{1}{\cancel{2}}x^2+C_1 \end{align*}

    より

    log⁡∣y∣=x2+C1\begin{align*} \log |y|=x^2+C_1 \end{align*}

     左辺は積分公式

    ∫1y dy=log⁡∣y∣+C\begin{align*} \int \dfrac{1}{y}  dy=\log |y|+C \end{align*}

    を用います。右辺は xx を xx 軸方向に積分し、積分定数は C1C_1 と1つにして右辺にまとめます。
     ここで、対数と指数の関係(中学でも分かる➄「対数関数」)

    log⁡aR=r ⟺ R=ar\begin{align*} \log_a R=r  \Longleftrightarrow R=a^r \end{align*}

    より

    log⁡∣y∣=x2+C1 ⟺ ∣y∣=ex2+C1\begin{align*} \log |y|=x^2+C_1  \Longleftrightarrow |y|=e^{x^2+C_1} \end{align*}

    と、∣y∣|y| について求めることができ

    ∣y∣=ex2+C1=ex2eC1=eC1ex2\begin{align*} |y|=e^{x^2+C_1}=e^{x^2}e^{C_1}=e^{C_1}e^{x^2} \end{align*}

     絶対値をはずして

    y=±eC1ex2\begin{align*} y=\pm e^{C_1}e^{x^2} \end{align*}

     定数部分 ±eC1\pm e^{C_1} を CC で置き換えて

    y=Cex2\begin{align*} y=Ce^{x^2} \end{align*}

     定数部分 ±eC1\pm e^{C_1} は、1つの文字 CC でまとめて置き換えます。最後に CC を使うために、途中の計算では C1C_1 を使っています。
     なお、絶対値をはずすと「±\pm」(プラスマイナス)がつくことにも注意。例えば、∣x∣=2|x|=2 の絶対値をはずすと
     x=±2x=\pm2
    となります。

    (case2) y=0y=0 のとき
     y=0y=0 を問題の微分方程式

    dydx=2xy\begin{align*} \dfrac{dy}{dx}=2xy \end{align*}

    に代入すると

    dydx=2x×0\begin{align*} \dfrac{dy}{dx}=2x\times0 \end{align*}

    より

    dydx=0\begin{align*} \dfrac{dy}{dx}=0 \end{align*}

    となります。
     すると、y=0y=0 を xx 軸方向に微分すると dydx=0\dfrac{dy}{dx}=0 より、y=0y=0 はこの微分方程式の解になりますが、y=0y=0 は (case1) で求めた
     y=Cex2y=Ce^{x^2}
    について、C=0C=0 とすれば
     y=0⋅ex2=0y=0\cdot e^{x^2}=0
    より得ることができます。つまり y=0y=0 の場合は、y=Cex2y=Ce^{x^2} に含まれています。
     よって一般解は、(case1)、(case2) より、CC を任意の定数として

    y=Cex2 (答)\begin{align*} y=Ce^{x^2} \small{(答)} \end{align*}

    ■\blacksquare

     では、検算してみましょう。合成関数の微分法の復習になります。
     「外の微分」×\times「中の微分」
    で微分します(中学でも分かる➈「合成関数」)。

    y′=(Cex2)′=C(ex2)′=Cex2外の微分⋅(x2)′中の微分=Cex2‾一般解 y⋅2x=y⋅2x=2xy\begin{align*} y'&=\left(Ce^{x^2}\right)'\\ &=C\left(e^{x^2}\right)'\\ &=C\underset{\tiny 外の微分}{e^{x^2}}\hspace{-4pt}\cdot\hspace{0pt}\underset{\tiny 中の微分}{\left(x^2\right)'}\\ &=\underset{\tiny 一般解  y}{\underline{Ce^{x^2}}}\cdot2x\\ &=y\cdot2x\\ &=2xy \end{align*}

    となり、問題の微分方程式と同じになります。外の微分は変化しないことに注意(中学でも分かる➈「合成関数」)。
     途中で、一般解 y=Cex2y=Ce^{x^2} そのものが現れることにも注意してください(下線部)。その部分を yy で置き換えることによって、問題の微分方程式が得られます。

    (2) 一般解に x=0, y=3x=0,  y=3 を代入して

    3=Ce023=Ce03=C⋅1C=3\begin{align*} &3= Ce^{0^2}\\ &3= Ce^{0}\\ &3=C\cdot1\\ &C=3 \end{align*}

     これを一般解に代入して

    y=3ex (答)\begin{align*} y=3e^x \small{(答)} \end{align*}

    ■\blacksquare

     (例題4)のような形の微分方程式を変数分離形といいます。分離という理由は

    dydx=2x×y↑↑x の関数y の関数\begin{align*} \dfrac{dy}{dx}=&2x \times y\\[-5pt] &\uparrow\hspace{15pt}\uparrow\\[-3pt] \footnotesize{x  の関数}\hspace{-13pt}&\hspace{18pt}\footnotesize{y  の関数} \end{align*}

    と

    (x の関数)×(y の関数)\begin{align*} (x  の関数) \times (y  の関数) \end{align*}

    の積の形に分離されていることからきています。例えば、次のような微分方程式も変数分離形です。

    dydx=(ex−2)y2\begin{align*} \frac{dy}{dx}=(e^x-2)y^2 \end{align*}

     これも

    dydx=(ex−2)×y2↑↑x の関数y の関数\begin{align*} \frac{dy}{dx}&=(e^x-2)\times y^2\\[-5pt] &\hspace{25pt}\uparrow\hspace{27pt}\uparrow\\[-3pt] \footnotesize{x  の関数}\hspace{-42pt}&\hspace{52pt}\footnotesize{y  の関数} \end{align*}

    と

    (x の関数)×(y の関数)\begin{align*} (x  の関数) \times (y  の関数) \end{align*}

    の積の形に分離されています。一般的に書くと

    dydx=P(x)Q(y)\begin{align*} \frac{dy}{dx}=P(x)Q(y) \end{align*}

    という形です。xx の関数 P(x)P(x) と yy の関数 Q(y)Q(y) の積の形になっているのが変数分離形です。
     例えば、次のような場合は P(x)P(x) と Q(y)Q(y) の積に分かれていないので、変数分離形ではありません。

    dydx=(xy−3)xy3\begin{align*} \frac{dy}{dx}=(xy-3)xy^3 \end{align*}

     xx と yy が混じって分離されていません。
     変数分離形の解法は

    dydx=P(x)Q(y)\begin{align*} \frac{dy}{dx}=P(x)Q(y) \end{align*}

    の両辺を Q(y)Q(y) で割って

    1Q(y)dydx=P(x)\begin{align*} \dfrac{1}{Q(y)}\frac{dy}{dx}=P(x) \end{align*}

    としてから、両辺を xx 軸方向で積分するのが基本的な解法となります。
     それでは、変数分離系の微分方程式をいくつか解いていきます。まずは、先ほどあげた具体例を解いてみましょう。

    (例題5)やや複雑

     次の微分方程式

    dydx=(ex−2)y2\begin{align*} \frac{dy}{dx}=(e^x-2)y^2 \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=0x=0 のとき y=1y=1 を満たすときの特殊解を求めなさい。

    指針
     (例題4)のように、両辺を y2y^2 で割って、yy の式を左辺にもっていきます。なお、両辺を y2y^2 で割る前に
    (case1) y≠0y\ne0 のとき
    (case2) y=0y=0 のとき
    で場合分けをしましょう。
     y2=0↔y=0y^2=0\leftrightarrow y=0
    は正しく、および
     y2≠0↔y≠0y^2\ne0\leftrightarrow y\ne0
    も正しいので、y≠0y\ne0 と y=0y=0 で場合分けをします。

    解答
    (1)
    (case1) y≠0y\ne0 のとき
     両辺を y2 (≠0)y^2  (\ne0) で割って

    1y2dydx=ex−2\begin{align*} \frac{1}{y^2}\frac{dy}{dx}=e^x-2 \end{align*}

     両辺を xx 軸方向に積分して

    ∫1y2dydx dx=∫(ex−2) dx∫1y2 dy=∫(ex−2) dx∫y−2 dy=∫(ex−2) dx\begin{align*} \int\frac{1}{y^2}\frac{dy}{\cancel{dx}}  \cancel{dx}&=\int(e^x-2)  dx\\ \int\frac{1}{y^2}  dy&=\int(e^x-2)  dx\\ \int y^{-2}  dy&=\int(e^x-2)  dx \end{align*}

     左辺は yy 軸方向に、右辺は xx 軸方向に積分して

    1−2+1y−2+1=ex−2x+C11−1y−1=ex−2x+C1−y−1=ex−2x+C1−1y=ex−2x+C11y=−ex+2x−C1\begin{align*} \dfrac{1}{-2+1}y^{-2+1}&=e^x-2x+C_1\\ \dfrac{1}{-1}y^{-1}&=e^x-2x+C_1\\ -y^{-1}&=e^x-2x+C_1\\ -\dfrac{1}{y}&=e^x-2x+C_1\\ \dfrac{1}{y}&=-e^x+2x-C_1 \end{align*}

     両辺の逆数をとって

    y=1−ex+2x−C1\begin{align*} y=\dfrac{1}{-e^x+2x-C_1} \end{align*}

     −C1=C-C_1=C とおくと

    y=1−ex+2x+C ⋯(答)\begin{align*} y=\dfrac{1}{-e^x+2x+C} \cdots\small{(答)} \end{align*}

    (case2) y=0y=0 のとき
     y=0y=0 を問題の微分方程式

    dydx=(ex−2)y2\begin{align*} \dfrac{dy}{dx}=(e^x-2)y^2 \end{align*}

    に代入すると

    dydx=0\begin{align*} \dfrac{dy}{dx}=0 \end{align*}

    となります。
     すると、y=0y=0 を xx 軸方向に微分すると dydx=0\dfrac{dy}{dx}=0 より、y=0y=0 はこの微分方程式の解になりますが、(case1) で求めた
     y=1−ex+2x+Cy=\dfrac{1}{-e^x+2x+C}
    に、y=0y=0 は含まれません(CC がどのような実数であっても y=0y=0 にならないため)。
     よって一般解は、(case1)、(case2) より、CC は任意の定数として

    y=1−ex+2x+C または y=0 ⋯(答)\begin{align*} y=\dfrac{1}{-e^x+2x+C} または y=0 \cdots\small{(答)} \end{align*}

    ■\blacksquare

     y=0y=0 が一般解に含まれないときは、上記のように別扱いにして解とします。

     では、検算してみましょう。
     「外の微分」×\times「中の微分」
    で微分します(中学でも分かる➈「合成関数」)。

    y=1−ex+2x+C=(−ex+2x+C)−1\begin{align*} y&=\dfrac{1}{-e^x+2x+C}\\ &=(-e^x+2x+C)^{-1} \end{align*}

    より

    y′=((−ex+2x+C)−1)′=−(−ex+2x+C)−1−1外の微分⋅(−ex+2x+C)′中の微分=−(−ex+2x+C)−2⋅(−ex+2)=−(−ex+2)⋅1(−ex+2x+C)2=(ex−2)(1−ex+2x+C)2この部分が一般解 y=(ex−2)y2\begin{align*} y'&=\left((-e^x+2x+C)^{-1}\right)'\\ &=\underset{\scriptsize 外の微分}{-(-e^x+2x+C)^{-1-1}}\cdot\underset{\scriptsize 中の微分}{\left(-e^x+2x+C\right)'}\\ &=-(-e^x+2x+C)^{-2}\cdot(-e^x+2)\\ &=-(-e^x+2)\cdot\dfrac{1}{(-e^x+2x+C)^2}\\ &=(e^x-2)\underset{\scriptsize この部分が一般解  y}{\left(\dfrac{1}{-e^x+2x+C}\right)^2}\\ &=(e^x-2)y^2 \end{align*}

    となり、問題の微分方程式と同じになります。途中で一般解そのものが現れることに注意してください。後ろのカッコの中身が (case1) で求めた一般解です。その部分を yy で置き換えることによって、問題の微分方程式が得られます。

    (2) x=0,  y=1x=0,   y=1 を代入して

    1=1−e0+2⋅0+C1=1−1+C\begin{align*} 1&=\dfrac{1}{-e^0+2\cdot0+C}\\[6pt] 1&=\dfrac{1}{-1+C}\\ \end{align*}

     両辺に −1+C-1+C をかけて

    −1+C=1\begin{align*} -1+C=1 \end{align*}

    より

    C=1+1=2\begin{align*} C=1+1=2 \end{align*}

     これを一般解に代入して

    y=1−ex+2x+2 ⋯(答)\begin{align*} y=\dfrac{1}{-e^x+2x+2} \cdots\small{(答)} \end{align*}

    ■\blacksquare

    (例題6)右辺が y のみ

     次の微分方程式

    dydx=y\begin{align*} \frac{dy}{dx}=y \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=0x=0 のとき y=1y=1 を満たすときの特殊解を求めなさい。

    指針
     右辺が yy の1変数のみですが、これも変数分離形の微分方程式です。

    dydx=1×y\begin{align*} \frac{dy}{dx}=1\times y \end{align*}

    と分離できます。11 は xx の関数と考えます。どのような xx に対しても定数 11 をとる関数です。
     よって、両辺を yy で割って、yy を左辺にもっていきます。なお、両辺 yy で割る前に
    (case1) y≠0y\ne0 のとき
    (case2) y=0y=0 のとき
    で場合分けをするのは先ほどと同じです。

    解答
    (1)
    (case1) y≠0y\ne0 のとき
     両辺を y (≠0)y  (\ne0) で割って

    1ydydx=1\begin{align*} \frac{1}{y}\frac{dy}{dx}=1 \end{align*}

     両辺を xx 軸方向に積分して

    ∫1ydydx dx=∫1 dx∫1y dy=∫1 dx\begin{align*} \int\frac{1}{y}\frac{dy}{\cancel{dx}}  \cancel{dx}&=\int1  dx\\ \int\frac{1}{y}  dy&=\int1  dx \end{align*}

     左辺は yy 軸方向に、右辺は xx 軸方向に積分して

    log⁡∣y∣=x+C1\begin{align*} \log |y|=x+C_1 \end{align*}

     左辺は積分公式

    ∫1y dy=log⁡∣y∣+C\begin{align*} \int \dfrac{1}{y}  dy=\log |y|+C \end{align*}

    を用います。右辺は 11 を xx 軸方向に積分し、積分定数は CC と1つにして右辺にまとめます。
     ここで、対数と指数の関係(中学でも分かる➄「対数関数」)

    log⁡aR=r ⟺ R=ar\begin{align*} \log_a R=r  \Longleftrightarrow R=a^r \end{align*}

    より

    log⁡e∣y∣=x+C1 ⟺ ∣y∣=ex+C1\begin{align*} \log_e |y|=x+C_1  \Longleftrightarrow |y|=e^{x+C_1} \end{align*}

    と、∣y∣|y| について求めることができ

    ∣y∣=ex+C1=exeC1=eC1ex\begin{align*} |y|=e^{x+C_1}=e^xe^{C_1}=e^{C_1}e^x \end{align*}

     絶対値をはずして

    y=±eC1ex\begin{align*} y=\pm e^{C_1}e^x \end{align*}

     定数部分 ±eC1\pm e^{C_1} を CC で置き換えて

    y=Cex\begin{align*} y=Ce^x \end{align*}

    (case2) y=0y=0 のとき
     y=0y=0 を問題の微分方程式

    dydx=y\begin{align*} \dfrac{dy}{dx}=y \end{align*}

    に代入すると

    dydx=0\begin{align*} \dfrac{dy}{dx}=0 \end{align*}

    となります。
     すると、y=0y=0 を xx 軸方向に微分すると dydx=0\dfrac{dy}{dx}=0 より、y=0y=0 はこの微分方程式の解になりますが、y=0y=0 は (case1) で求めた
     y=Cexy=Ce^x
    について、C=0C=0 とすれば
     y=0⋅ex=0y=0\cdot e^x=0
    より得ることができます。つまり y=0y=0 の場合は、y=Cexy=Ce^x に含まれます。
     よって一般解は、(case1)、(case2) より、CC は任意の定数として

    y=Cex (答)\begin{align*} y=Ce^x \small{(答)} \end{align*}

    ■\blacksquare

     明かですが、検算してみましょう。

    dydx=C(ex)′=Cex‾一般解 y=y\begin{align*} \dfrac{dy}{dx}&=C(e^x)'\\ &=\underset{\tiny 一般解  y}{\underline{Ce^x}}\\ &=y \end{align*}

    となり、問題の微分方程式と同じになります。
     exe^x は微分しても変化しないので、微分しても変わらないを意味する dydx=y\dfrac{dy}{dx}=y の解は、y=Cexy=Ce^x であることは予想がつくところです。

    (2) 一般解に x=0, y=3x=0,  y=3 を代入して

    3=Ce03=C⋅1C=3\begin{align*} &3= Ce^0\\ &3=C\cdot1\\ &C=3 \end{align*}

     これを一般解に代入して

    y=3ex (答)\begin{align*} y=3e^x \small{(答)} \end{align*}

    ■\blacksquare

    変数分離形の解法

     以上のように、変数分離形

    dydx=P(x)Q(y)\begin{align*} \frac{dy}{dx}=P(x)Q(y) \end{align*}

    は、両辺を Q(y)Q(y) で割って

    1Q(y)dydx=P(x)\begin{align*} \dfrac{1}{Q(y)}\dfrac{dy}{dx}=P(x) \end{align*}

    と、yy の関数 Q(y)Q(y) を全部左辺に持っていってから、両辺を xx 軸方向に積分します。

    ∫1Q(y)dydx dx=∫P(x)dx∫1Q(y) dy=∫P(x) dx\begin{align*} \int \dfrac{1}{Q(y)}\dfrac{dy}{\cancel{dx}}  \cancel{dx}&=\int P(x)dx\\[10pt] \int \dfrac{1}{Q(y)}  dy&=\int P(x)  dx \end{align*}

     すると、左辺は形式的に dxdx で約分されて yy 軸方向の積分になり、右辺はそのまま xx 軸方向の積分を行います。

    2階微分方程式

     次は2階微分の入った微分方程式をみてみましょう

    (例題7)基本形

     次の微分方程式

    d2ydx2=3x2+2\begin{align*} \frac{d^2y}{dx^2}=3x^2+2 \end{align*}

    について
    (1) 一般解を求めなさい。
    (2) x=0x=0 のとき y=1y=1、x=2x=2 のとき y=3y=3 を満たすときの特殊解を求めなさい。

    指針
     
    d2ydx2\dfrac{d^2y}{dx^2} は2階微分の意味です(中学で分かる①「微分法」)。つまり、この微分方程式は、
     「y=f(x)y=f(x) を2回連続微分すると 3x2+23x^2+2 になる」
    という意味です。よって、 yy を求めるには 3x2+23x^2+2 を2回連続積分します。

    → 微分する→ 微分するydydxd2ydx2=3x2+2← 積分する← 積分する\begin{align*} &\hspace{8pt}\xrightarrow{  微分する}\hspace{24pt}\xrightarrow{  微分する}\\[-19pt] &y\hspace{48pt}\dfrac{dy}{dx}\hspace{48pt}\dfrac{d^2y}{dx^2}=3x^2+2\\[-15pt] &\hspace{8pt}\xleftarrow[  積分する]{}\hspace{24pt}\xleftarrow[  積分する]{} \end{align*}

     この「←\leftarrow」の操作を行っていきます。
     また、この問題の式は

    y′′=3x2+2\begin{align*} y''=3x^2+2 \end{align*}

    と表されることもあります。xx で微分しているときは、y′=dydxy'=\dfrac{dy}{dx}、y′′=d2ydx2y''=\dfrac{d^2y}{dx^2} という意味になります。

    → 微分する→ 微分するyy′y′′=3x2+2← 積分する← 積分する\begin{align*} &\hspace{8pt}\xrightarrow{  微分する}\hspace{18pt}\xrightarrow{  微分する}\\[-13pt] &y\hspace{48pt}y'\hspace{48pt}y''=3x^2+2\\[-12pt] &\hspace{8pt}\xleftarrow[  積分する]{}\hspace{18pt}\xleftarrow[  積分する]{} \end{align*}

    解答
    (1) 両辺を xx 軸方向に積分すると

    ∫d2ydx2 dx=∫(3x2+2) dx\begin{align*} \int\frac{d^2y}{dx^2}  dx=\int(3x^2+2)  dx \end{align*}

    より

    dydx=3⋅13x3+2x+C1=x3+2x+C1\begin{align*} \frac{dy}{dx}&=\cancel{3}\cdot\dfrac{1}{\cancel{3}}x^3+2x+C_1\\ &=x^3+2x+C_1 \end{align*}

     よって

    dydx=x3+2x+C1\begin{align*} \frac{dy}{dx}=x^3+2x+C_1 \end{align*}

     もう一度 xx 軸方向に積分して

    ∫dydx dx=∫(x3+2x+C1) dx\begin{align*} \int\frac{dy}{dx}  dx=\int(x^3+2x+C_1)  dx \end{align*}

    より

    y=13+1x3+1+2⋅11+1x1+1+C1x+C2=14x4+2⋅12x2+C1x+C2=14x4+x2+C1x+C2\begin{align*} y&=\dfrac{1}{3+1}x^{3+1}+2\cdot\dfrac{1}{1+1}x^{1+1}+C_1x+C_2\\ &=\dfrac{1}{4}x^4+\cancel{2}\cdot\dfrac{1}{\cancel{2}}x^2+C_1x+C_2\\ &=\dfrac{1}{4}x^4+x^2+C_1x+C_2 \end{align*}

     よって一般解は、C1, C2C_1,  C_2 を任意の定数として

    y=14x4+x2+C1x+C2 ⋯(答)\begin{align*} y=\dfrac{1}{4}x^4+x^2+C_1x+C_2 \cdots\small{(答)} \end{align*}

    ■\blacksquare

     任意の定数は、C1, C2C_1,  C_2 と2つあることに注意。2回積分しているので、任意の定数は2つになります。

     では、検算してみましょう。
     一般解を微分して

    dydx=14⋅4x3+2x+C1=x3+2x+C1\begin{align*} \dfrac{dy}{dx}&=\dfrac{1}{4}\cdot4x^3+2x+C_1\\ &=x^3+2x+C_1 \end{align*}

     これをもう一度微分して

    d2ydx2=3x2+2\begin{align*} \dfrac{d^2y}{dx^2}=3x^2+2 \end{align*}

    となり、問題の微分方程式と同じになります。

    (2) 一般解に x=0, y=1x=0,  y=1 を代入すると

    1=14⋅04+02+C1⋅0+C2=C2\begin{align*} 1=\dfrac{1}{4}\cdot0^4+0^2+C_1\cdot0+C_2=C_2\\ \end{align*}

    より

    C2=1\begin{align*} C_2=1 \end{align*}

    一般解に x=2, y=3x=2,  y=3 を代入すると

    3=14⋅24+22+C1⋅2+C23=14⋅16+4+2C1+C23=4+4+2C1+C23=8+2C1+C2\begin{align*} 3&=\dfrac{1}{4}\cdot2^4+2^2+C_1\cdot2+C_2\\ 3&=\dfrac{1}{\cancel{4}}\cdot\cancel{16}+4+2C_1+C_2\\ 3&=4+4+2C_1+C_2\\ 3&=8+2C_1+C_2\\ \end{align*}

     C2=1C_2=1 を代入して

    3=8+2C1+13=9+2C19+2C1=32C1=3−9=−6\begin{align*} &3=8+2C_1+1\\ &3=9+2C_1\\ &9+2C_1=3\\ &2C_1=3-9=-6 \end{align*}

    より

    C1=−3\begin{align*} C_1=-3 \end{align*}

     C1=−3, C2=1C_1=-3,  C_2=1 を一般解に代入して

    y=14x3+x2−3x+1 ⋯(答)\begin{align*} y=\dfrac{1}{4}x^3+x^2-3x+1 \cdots{\small(答)} \end{align*}

    ■\blacksquare

     未知の定数は2つあるので、C1, C2C_1,  C_2 の決定のためには2組の x, yx,  y での f(x, y)f(x,  y) の値が必要です。2組の x, yx,  y を一般解に代入して、その連立方程式を解きます。

    まとめ

     これまでの結果をまとめましょう。C, C1, C2C,  C_1,  C_2 は任意の定数です。

    <直接積分形>

    (冒頭の例題)
     dydx=2x\dfrac{dy}{dx}=2x
    の一般解は
     y=x2+Cy=x^2+C

    (例題1)
     dydx=4x3−x−2\dfrac{dy}{dx}=4x^3-x-2
    の一般解は
     y=x4−12x2−2x+Cy=x^4-\dfrac{1}{2}x^2-2x+C

    (例題2)
     y′=exy'=e^x
    の一般解は
     y=ex+Cy=e^x+C

    (例題3)
     xy′=3xy'=3
    の一般解は
     y=3log⁡∣x∣+Cy=3\log |x|+C

    <変数分離形>

    (例題4)
     dydx=2xy\dfrac{dy}{dx}=2xy
    の一般解は
     y=Cex2y=Ce^{x^2}

    (例題5)
     
    dydx=(ex−2)y2\dfrac{dy}{dx}=(e^x-2)y^2
    の一般解は
     y=1−ex+2x+Cy=\dfrac{1}{-e^x+2x+C} または y=0y=0

    (例題6)
     dydx=y\dfrac{dy}{dx}=y
    の一般解は
     y=Cexy=Ce^x

    <2階微分方程式>
    (例題7)
     
    d2ydx2=3x2+2\dfrac{d^2y}{dx^2}=3x^2+2
    の一般解は
     y=14x4+x2+C1x+C2y=\dfrac{1}{4}x^4+x^2+C_1x+C_2

    ***

     直接積分形は、文字通り直接積分します。
     変数分離形は、両辺を yy の式で割ってから積分します。
     2階微分方程式は、2回積分します。

     ここであげた微分方程式は、最も簡単な微分方程式の一例です。さらに難しい微分方程式はいくらでもありますが、その式を満たす関数 y=f(x)y=f(x) を求めるという精神は同じになります。
     決まった形、決まった解法がそれぞれにあるので、興味があれば、各自専門書などで勉強してみてください。

    小説編


    この記事が参加している募集

    あなたへのおすすめ