§C1.6微分方程式入門

最終更新

微分方程式とは、未知の関数とその導関数を含む方程式のことです。 求めるのは「数」ではなく、その関係を満たす関数。自然現象(増減の速さが量に比例する、など)は微分方程式で記述されます。

1 変数分離形

dydx=g(x) h(y)\dfrac{dy}{dx} = g(x)\,h(y)のように、右辺が「xxの式」×「yyの式」に分かれる微分方程式を変数分離形といいます。このとき、h(y)≠0h(y) \ne 0の範囲で∫dyh(y)=∫g(x) dx (+ C)\int \frac{dy}{h(y)} = \int g(x)\,dx \ (+\, C)と解けます。1h(y)dydx=g(x)\dfrac{1}{h(y)}\dfrac{dy}{dx} = g(x)の両辺をxxで積分すると、左辺が置換積分で∫dyh(y)\int \frac{dy}{h(y)}になるからです。dy,dxdy, dxを形式的に移項して別々に積分する、と覚えると速い。積分定数は初期条件(x=x0x = x_0でy=y0y = y_0)で決めます。

2 例:指数的成長・減衰

代表例がdydx=ky\dfrac{dy}{dx} = ky(変化の速さが今の量に比例)です。h(y)=y≠0h(y) = y \ne 0で変数分離すると∫dyy=∫k dx\int \frac{dy}{y} = \int k\,dx、すなわちlog⁡∣y∣=kx+C\log|y| = kx + C、よってy=Aekxy = A e^{kx}(A=±eCA = \pm e^C)。y≡0y \equiv 0(h(y)=0h(y) = 0の定数解)も解で、A=0A = 0として一般解に含まれます。人口増加・放射性崩壊(k<0k < 0)・ニュートンの冷却など、指数関数的な現象がこの一本の式から出てきます。exe^xが「微分しても変わらない」性質がここで効いています。より本格的な存在一意性や解けない方程式の扱いは、大学の常微分方程式で展開されます。

例 2.1 (初期条件で積分定数を決める).y′=yy' = y、y(0)=2y(0) = 2を解きます。k=1k = 1の場合なので、y≠0y \ne 0の範囲でdyy=dx\dfrac{dy}{y} = dxと分けて積分するとlog⁡∣y∣=x+C\log|y| = x + C、すなわちy=Aexy = Ae^{x}です。ここへ初期条件y(0)=2y(0) = 2を入れるとA=2A = 2と決まり、求める解は

y=2exy = 2e^{x}

になります。定数解y≡0y \equiv 0はy(0)=2y(0) = 2を満たさないので、この初期条件のもとでは候補から外れます。

閑話休題:空のバケツは過去を語らない 微分方程式は「変化の法則から未来を予言する」道具ですが、その予言能力には面白い穴があります。底に穴の開いたバケツから水が抜けるとき、水位yyはおよそdydx=−ky\dfrac{dy}{dx} = -k\sqrt{y}に従います(トリチェリの法則)。時間を逆向きに見ればdydx=ky\dfrac{dy}{dx} = k\sqrt{y}で、いまy=0y = 0(空)だとしましょう。このバケツ、いつ空になったのかが方程式から復元できません。ずっと空のまま(y≡0y \equiv 0)も解なら、「11分前まで水があって、ちょうど空になった」も「11時間前に空になった」も、すべて同じ現在(y=0y = 0)に行き着く正当な解だからです。実際、y′=kyy' = k\sqrt{y}はy≡0y \equiv 0のほかに、任意の時刻ccから立ち上がる解y=k24(x−c)2y = \frac{k^2}{4}(x - c)^2(x≥cx \ge c)を持ちます。解が1本に決まらないのです。

犯人はy\sqrt{y}がy=0y = 0の際(きわ)で「急峻すぎる」こと(傾きの変化率が無限大に発散する)です。逆に、右辺がそこそこ穏やかなら解はただ1本に決まる、という解の存在と一意性の定理が大学の常微分方程式で待っています。ニュートン力学の世界は「現在の状態が過去と未来を完全に決める」と語られがちですが、それは方程式の右辺の行儀のよさという数学的な条件あってのこと。バケツ1杯の水が、決定論の境界線を見せてくれます。

例題

条件と何を求めるかを確認してから、式と答えの対応を見比べてください。

次の微分方程式を、初期条件のもとで解け(変数分離形)。

次の微分方程式を、初期条件のもとで解け。

解法の型dy/dx == f(x)g(y) を dy/g(y) == f(x) dx と分けて両辺を積分し、初期条件で積分定数を決める

  1. 例題 1

    dydx=y2,y(0)=14\dfrac{dy}{dx} = y^{2}, \qquad y(0) = \dfrac{1}{4}
  2. 例題 2

    dydx=2xy,y(0)=2\dfrac{dy}{dx} = 2xy, \qquad y(0) = 2
  3. 例題 3

    dydx=−2y,y(0)=1\dfrac{dy}{dx} = -2y, \qquad y(0) = 1
  4. 例題 4

    dydx=y2,y(0)=15\dfrac{dy}{dx} = y^{2}, \qquad y(0) = \dfrac{1}{5}
  5. 例題 5

    dydx=−4xy,y(0)=2\dfrac{dy}{dx} = -4xy, \qquad y(0) = 2
  6. 例題 6

    dydx=−2xy,y(0)=4\dfrac{dy}{dx} = -2xy, \qquad y(0) = 4
  7. 例題 7

    dydx=1+y2,y(0)=0\dfrac{dy}{dx} = 1 + y^{2}, \qquad y(0) = 0
  8. 例題 8

    dydx=y(1−y),y(0)=12\dfrac{dy}{dx} = y(1 - y), \qquad y(0) = \dfrac{1}{2}
  9. 例題 9

    dydx=−4xy,y(0)=1\dfrac{dy}{dx} = -4xy, \qquad y(0) = 1
  10. 例題 10

    dydx=y2,y(0)=13\dfrac{dy}{dx} = y^{2}, \qquad y(0) = \dfrac{1}{3}

演習

問題を解いてから「解答・解説」を開けます。

次の微分方程式を、初期条件のもとで解け(変数分離形)。

演習を読み込み中…

前提記事