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

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

     前回は微分方程式を勉強しました。ここでは、偏微分方程式を勉強します。偏微分方程式とは、文字通り偏微分の入った方程式です。
     まずは、中学でも分かる①「微分法」でやった偏微分を加筆修正して再掲し、復習します。その後、あらたに偏微分方程式を勉強します。難しい内容は避けて、偏微分方程式とは何かについて、その基本のみを紹介します。
     最後に、物理で現れる微分方程式
    ・ニュートンの運動方程式
    ・シュレーディンガー方程式
    を簡単に紹介します。

    偏微分(基本)

     f(x)=x2f(x)=x^2
    という関数は、変数が xx だけの1変数関数です。これを xx で微分すると

    f′(x)=2x\begin{align*} f'(x)=2x \end{align*}

     または

    (x2)′=2xdf(x)dx=2xddxf(x)=2x\begin{align*} (x^2)'&=2x\\ \dfrac{df(x)}{dx}&=2x\\ \dfrac{d}{dx}f(x)&=2x \end{align*}

    と表せることを、中学でも分かる①「微分法」で解説しました。

    2変数関数

     関数には、次のように変数が 22 つの関数もあります。これを2変数関数といいます。

    f(x, y)=3x2y5\begin{align*} f(x,  y)=3x^2y^5 \end{align*}

     これは、x, yx,  y の2つの変数をもつ関数です。x, yx,  y の値が決まると、z=f(x, y)z=f(x,  y) という関数によって zz の値がただ1つに決まります。

    x, y が決まる→ z=f(x, y)  z の値がただ1つに決まる\begin{align*} x,  y  が決まる \xrightarrow[]{  z=f(x,  y)  }  z  の値がただ1つに決まる \end{align*}

    (例)f(x, y)=3x2y5f(x,  y)=3x^2y^5 について、z=f(x, y)z=f(x,  y) とおくと
     x=2, y=1x=2,  y=1 のとき
     z=f(2, 1)=3⋅22⋅15=3⋅4⋅1=12z=f(2,  1)=3\cdot2^2\cdot1^5=3\cdot4\cdot1=12

    2変数関数の微分

     ここで、f(x, y)=3x2y5f(x,  y)=3x^2y^5 の微分を考えます。この場合、xx で微分するか、yy で微分するかで2つの微分があります。

    ➀ yy を定数とみなして xx で微分する。
     これは、∂f(x, y)∂x\dfrac{\partial f(x,  y)}{\partial x} という記号を用いて

    ∂f(x, y)∂x=∂∂x3x2y5=3y5ddxx2=3y5⋅2x=6xy5\begin{align*} \dfrac{\partial f(x,  y)}{\partial x}&=\dfrac{\partial}{\partial x}3x^2y^5\\ &=3y^5\dfrac{d}{dx}x^2\\ &=3y^5\cdot2x\\ &=6xy^5 \end{align*}

     yy を定数とみなし、xx の関数部分を xx で微分します。

    画像

     yy は定数扱いなので、3y53y^5 をまとめて係数(定数)として扱い、偏微分の外に出します。
     このように、2変数関数に関して、1つの変数を定数とみなし、もう1つの変数で微分したものを偏微分といいます。

    偏微分の記号

     偏微分で使う「∂\partial」という記号は、筆記体の「dd」("ディー")が元になっており、”ラウンドディー” や ”デル” と読みます。
     偏微分を表すにあたって
     ∂f(x, y)∂x\dfrac{\partial f(x,  y)}{\partial x} や ∂∂xf(x, y)\dfrac{\partial}{\partial x}f(x,  y)
    のように、分母の xx によって、xx で偏微分することがはっきりと明示されます。これを、分母を yy にして
     ∂f(x, y)∂y\dfrac{\partial f(x,  y)}{\partial y} や ∂∂yf(x, y)\dfrac{\partial}{\partial y}f(x,  y)
    と書けば、次の➁で述べる「yy で偏微分する」という意味になります。

     なお
     df(x)dx\dfrac{df(x)}{dx} や ddxf(x)\dfrac{d}{dx}f(x)
    は高校で習う微分で、この微分を常微分といいます(中学でも分かる①「微分法」)。

     先ほどの①の計算をもう一度みてみましょう。
     xx での偏微分の場合、xx 以外の部分を偏微分の外に出すことによって、最後は xx での常微分 ddxf(x)\dfrac{d}{dx}f(x) になります。

    (再掲)

    ∂f(x, y)∂x=∂∂x3x2y5 ←x での偏微分=3y5ddxx2  ←x での常微分=3y5⋅2x=6xy5\begin{align*} \dfrac{\partial f(x,  y)}{\partial x}&=\dfrac{\partial}{\partial x}3x^2y^5 \leftarrow {\footnotesize x  での偏微分}\\ &=3y^5\dfrac{d}{dx}x^2 \,\leftarrow {\footnotesize x  での常微分}\\ &=3y^5\cdot2x\\ &=6xy^5 \end{align*}

     では次に、yy での偏微分をみていきます。
     f(x, y)=3x2y5f(x,  y)=3x^2y^5 について

    ➁ xx を定数とみなして yy で微分(偏微分)する。

    ∂f(x, y)∂y=∂∂y3x2y5 ←y での偏微分=3x2ddyy5  ←y での常微分=3x2⋅5y4=15x2y4\begin{align*} \dfrac{\partial f(x,  y)}{\partial y}&=\dfrac{\partial}{\partial y}3x^2y^5 \leftarrow {\footnotesize y  での偏微分}\\ &=3x^2\dfrac{d}{dy}y^5 \,\leftarrow {\footnotesize y  での常微分}\\ &=3x^2\cdot5y^4\\ &=15x^2y^4 \end{align*}

     xx を定数とみなし、yy の関数部分を yy で微分(常微分)します。

    画像

     xx は定数扱いなので、3x23x^2 をまとめて係数(定数)として扱い、偏微分の外に出します。最後に yy での常微分が残ります。

     まとめると、f(x, y)=3x2y5f(x,  y)=3x^2y^5 について
    (1) xx で偏微分すると

    ∂f(x, y)∂x=3y5ddxx2=3y5⋅2x=6xy5\begin{align*} \dfrac{\partial f(x,  y)}{\partial x}=3y^5\dfrac{d}{dx}x^2=3y^5\cdot2x=6xy^5 \end{align*}

    (2) yy で偏微分すると

    ∂f(x, y)∂y=3x2ddyy5=3x2⋅5y4=15x2y4\begin{align*} \dfrac{\partial f(x,  y)}{\partial y}=3x^2\dfrac{d}{dy}y^5=3x^2\cdot5y^4=15x^2y^4 \end{align*}

    と、結果が異なることに注意しましょう。

    多変数関数の偏微分

     3変数以上の関数についても偏微分を考えることができます。
     変数が2つ以上の関数を多変数関数といいます。多変数関数の偏微分は、ある1つの変数のみに着目し、それ以外の変数は定数とみなして微分します。

    (例1)f(x, y, z)=5x2y3+z4f(x,  y,  z)=5x^2y^3+z^4 について xx で偏微分すると

    ∂∂xf(x, y, z)=∂∂x(5x2y3+z4)=∂∂x5x2y3+∂∂xz4=5y3ddxx2+0=5y3⋅2x=10xy3\begin{align*} \dfrac{\partial}{\partial x}f(x,  y,  z)&=\dfrac{\partial}{\partial x}(5x^2y^3+z^4)\\ &=\dfrac{\partial}{\partial x}5x^2y^3+\dfrac{\partial}{\partial x}z^4\\ &=5y^3\dfrac{d}{dx}x^2+0\\ &=5y^3\cdot2x\\ &=10xy^3 \end{align*}

     zz は定数とみなしているので

    ∂∂xz4=0\begin{align*} \dfrac{\partial}{\partial x}z^4=0 \end{align*}

    となります。

    (例2)f(x, y, z)=5x2y3+z4f(x,  y,  z)=5x^2y^3+z^4 について yy で偏微分すると

    ∂∂yf(x, y, z)=∂∂y(5x2y3+z4)=∂∂y5x2y3+∂∂yz4=5x2ddyy3+0=5x2⋅3y2=15x2y2\begin{align*} \dfrac{\partial}{\partial y}f(x,  y,  z)&=\dfrac{\partial}{\partial y}(5x^2y^3+z^4)\\ &=\dfrac{\partial}{\partial y}5x^2y^3+\dfrac{\partial}{\partial y}z^4\\ &=5x^2\dfrac{d}{dy}y^3+0\\ &=5x^2\cdot3y^2\\ &=15x^2y^2 \end{align*}

     zz は定数とみなしているので

    ∂∂yz4=0\begin{align*} \dfrac{\partial}{\partial y}z^4=0 \end{align*}

    となります。

    (例3)f(x, y, z)=5x2y3+z4f(x,  y,  z)=5x^2y^3+z^4 について zz で偏微分すると

    ∂∂xf(x, y, z)=∂∂x(5x2y3+z4)=∂∂x5x2y3+ddzz4=0+4z3=4z3\begin{align*} \dfrac{\partial}{\partial x}f(x,  y,  z)&=\dfrac{\partial}{\partial x}(5x^2y^3+z^4)\\ &=\dfrac{\partial}{\partial x}5x^2y^3+\dfrac{d}{dz}z^4\\ &=0+4z^3\\ &=4z^3 \end{align*}

     x, yx,  y は定数とみなしているので

    ∂∂x5x2y3=0\begin{align*} \dfrac{\partial}{\partial x}5x^2y^3=0 \end{align*}

    となります。

    偏微分の意味(引用)

     偏微分について、書籍「ブラックホールと時空の方程式(小林晋平著)」(P82~P83)に非常にわかりやすい優れた解説があるので、(本記事に合うように多少文章を変えて)引用します。

    (引用始め)

     微分とは、何らかの量の微小変化を表すものです。 たとえば、
    時刻 tt での物体の位置(座標) x(t)x(t) の時間変化は dx(t)dt\dfrac{dx(t)}{dt}であり、これは速度 v(t)v(t) を表します。

    dx(t)dt=v(t)\begin{align*} \dfrac{dx(t)}{dt}=v(t) \end{align*}

     位置は x(t)x(t) という書き方からもわかるように時刻のみの関数です。 つまり、時刻 tt が変化すると、それに応じて位置が変化することを表しています。 時刻 tt を動かすとそれにつられて x(t)x(t) も動くという関係があるので、tt を独立変数、xx を従属変数といいます。
     位置 xx がそうであるように、1つの独立変数によって決まる 従属変数のことを1変数関数と言いますが、私たちが日常で見かけるいろいろな物理量にはそれほど単純ではないものもたくさんあります。
     たとえば世界の気温を考えてみましょう。気温は場所や時間に応じてどんどん変化します。今あなたが赤道直下のどこか暑い国にいるとして、赤道に沿って進めば相変わらず気温は高いままでしょうし、北に向かって進んでいけば少しずつ気温は低くなっていくはずです。 途中に山があり、それを登って高いところに行ってもやはり気温は低くなるでしょう。
     赤道方向に沿って x(t)x(t) 軸,、北に向かって yy 軸、空に向かって zz 軸をとれば、気温 TT は x, y, zx,  y,  z のそれぞれによって異なる値をとるので

    T=T(t, x, y, z)\begin{align*} T=T(t,  x,  y,  z) \end{align*}

    と書けます。このように、少なくとも気温は時刻 tt および位置 x, y, zx,  y,  z に依存する 4 変数関数のはずです。このように2つ以上の独立変数に依存して決まる関数を多変数関数と言いますが、 多変数関数の変化は、微分でどのように表されるのでしょうか。
     1変数関数 f(x)f(x) の微分の場合、1つの変数 xx に依存しているので、変化と言えば
     xx を変化させたときに関数 f(x)f(x) がどう変化するか
    を指します。これに対し、気温 T(t, x, y, z)T(t,  x,  y,  z) の場合は、t,x,y,zt, x, y, z という要素のどれを変化させるかによって変化の様子は違ってくるはずです。
     たとえば赤道方向(xx 軸方向)に沿って動いても相変わらず気温は高いままでしょうから、xx 軸方向への 変化はあまりなさそうです。一方、 北上したり(yy 軸の正方向への移動)、山に登っ たり(zz 軸の正方向への移動)すると気温は下がるでしょうから、そちら方向への変化は負の値になりそうです。こうした 「変数ごとの変化」を表すのに役立つのが偏微分です。具体的には、xx 方向に移動したとき気温 TT がどう変化するかを

    ∂∂x T(t, x, y, z)\begin{align*} \dfrac{\partial}{\partial x}\,T(t,  x,  y,  z) \end{align*}

    のように書き表して、 通常の微分 (常微分と言います)と区別します。
     偏微分は「特定の量の変化分だけに注目し、それ以外の量の変化は一切無視する」という計算なので、xx 方向への偏微分であれば xx 以外の変数 t,y,zt, y, z は定数だと思って計算します。
     たとえば f(x, y)=3x+2y+4xyf(x,  y)=3x+2y+4xy という2変数関数があれば、f(x, y)=3x+2y+4xy2f(x,  y)=3x+2y+4xy^2 を xx で偏微分したものは

    ∂∂xf(x, y)=3+0+4y2⋅1=3+4y2\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=3+0+4y^2\cdot1=3+4y^2 \end{align*}

     f(x, y)=3x+2y+4xy2f(x,  y)=3x+2y+4xy^2 を yy で偏微分したものは

    ∂∂xf(x, y)=0+2+4x⋅2y=2+8xy\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=0+2+4x\cdot2y=2+8xy \end{align*}

    となります。

    (引用終わり) 

    偏微分方程式

     偏微分を含んでいる方程式を偏微分方程式といいます。また、常微分を含んでいる方程式を常微分方程式といいます。
     2つの変数をもつ関数として f(x,y)f(x, y) を考えると、例えば

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    のような式が偏微分方程式です。この式の左辺は f(x, y)f(x,  y) の xx での偏微分、右辺は f(x, y)f(x,  y) の yy での偏微分です。
     すると、この式の意味は、
     「f(x, y)f(x,  y) を xx で偏微分したものと、f(x, y)f(x,  y) を yy で偏微分したものは等しい」
    という意味の偏微分方程式です。これは、次のように表記される場合もあります。

    ∂f(x, y)∂x=∂f(x, y)∂y\begin{align*} \dfrac{\partial f(x,  y)}{\partial x}=\dfrac{\partial f(x,  y)}{\partial y} \end{align*}

     では、この偏微分方程式を解いて一般解 f(x, y)f(x,  y) を求めてみましょう。

    問題 次の偏微分方程式

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    について

    (1) f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) として一般解を求めなさい。
    (2) f(0, 0)=1, f(1, 0)=ef(0,  0)=1,  f(1,  0)=e を満たすときの特殊解を求めなさい。

    指針
     f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) して一般解を求める方法を変数分離法といいます。一般解を

    f(x, y)=(x の関数)g(x)×(y の関数)h(y)\begin{align*} f(x,  y)=\underset{g(x)}{(x  の関数)} \times \underset{h(y)}{(y  の関数)} \end{align*}

    と仮定して求める方法です。このように変数分離された解があるかないかを探ることによって、一般解が容易に求まる場合があります。
     なお、変数分離法では解が求まらない場合もありますが、その場合は他の方法を探ることになります。本問題は変数分離法で解ける問題です。

    解答

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    に、f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) を代入して

    ∂∂xg(x)h(y)=∂∂yg(x)h(y)\begin{align*} \dfrac{\partial}{\partial x}g(x)h(y)=\dfrac{\partial}{\partial y}g(x)h(y) \end{align*}

     左辺と右辺に分けて計算します。

    (左辺)について

    (左辺)=∂∂xg(x)h(y) ←x での偏微分=h(y)ddxg(x) ←x での常微分\begin{alignat*}{2} (左辺)&=\dfrac{\partial}{\partial x} g(x)h(y)& &\leftarrow {\footnotesize x  での偏微分}\\ &=h(y)\dfrac{d}{d x}g(x)& &\leftarrow {\footnotesize x  での常微分} \end{alignat*}

     これは xx についての偏微分なので、yy の関数 h(y)h(y) は偏微分の外に出せます。yy は定数扱いです。すると、g(x)g(x) の xx での常微分 ddxg(x)\dfrac{d}{dx}g(x) が最後に残ります。さらに、簡単のため ddxg(x)=g′(x)\dfrac{d}{dx}g(x)=g'(x) とおくことで

    (左辺)=h(y)g′(x) ⋯(∗1)\begin{align*} (左辺)=h(y)g'(x) \cdots(*1) \end{align*}

    となります。

    (右辺)について
     同様に dh(y)dy=h′(y)\dfrac{dh(y)}{dy}=h'(y) とおくと、右辺について

    (左辺)=∂∂yg(x)h(y) ←y での偏微分=g(x)ddyh(y)  ←y での常微分=g(x)h′(y) ⋯(∗2)\begin{align*} (左辺)&=\dfrac{\partial}{\partial y}g(x)h(y) \leftarrow {\footnotesize y  での偏微分}\\ &=g(x)\dfrac{d}{dy}h(y) \,\leftarrow {\footnotesize y  での常微分}\\ &=g(x)h'(y) \cdots(*2) \end{align*}

     (∗1)=(∗2)(*1)=(*2) より

    h(y)g′(x)=g(x)h′(y)\begin{align*} h(y)g'(x)=g(x)h'(y) \end{align*}

     ここで
    (case1) g(x)h(y)≠0g(x)h(y)\ne0 のとき
    (case2) g(x)h(y)=0g(x)h(y)=0 のとき
    で場合分けをします。

    (case1) g(x)h(y)≠0g(x)h(y)\ne0 のとき
     両辺を g(x)h(y)(≠0)g(x)h(y)(\ne0) で割って

    h(y)g′(x)g(x)h(y)=g(x)h′(y)g(x)h(y)g′(x)g(x)=h′(y)h(y)\begin{align*} \dfrac{\cancel{h(y)}g'(x)}{g(x)\cancel{h(y)}}&=\dfrac{\cancel{g(x)}h'(y)}{\cancel{g(x)}h(y)}\\ \dfrac{g'(x)}{g(x)}&=\dfrac{h'(y)}{h(y)} \end{align*}

     すると、左辺は変数 xx だけの関数、右辺は変数 yy だけの関数になります。xx と yy は互いに独立(お互いに関係なく、xx と yy は自由に値をとれる)なのに、左辺と右辺がつねに等しいので、この式を満たすためには、左辺と右辺は

    g′(x)g(x)=Kh′(y)h(y)=K\begin{align*} \dfrac{g'(x)}{g(x)}&=K\\ \dfrac{h'(y)}{h(y)}&=K\\ \end{align*}

    と、同じ定数値 KK をならなければなりません。
     左辺は xx のみの関数であり、右辺は yy のみの関数です。xx と yy は互いに独立に動かせるので、この等式がすべての x, yx,  y で成り立つためには、左辺も右辺もともに一定の値でなければなりません。よって、その共通の定数を KK と置きます。
     抽象的なので、具体例をあげて解説します。

    (追記)直感的には理解できるのですが、下記の解説は、例の取り方も含めてなにか腑に落ちない感じがあるので、もっとうまく言語化できないか調整中です。余計なら外します。 

    (詳説)
     仮に、最も簡単な例として

    g′(x)g(x)=xh′(y)h(y)=y\begin{align*} \dfrac{g'(x)}{g(x)}&=x\\ \dfrac{h'(y)}{h(y)}&=y\\ \end{align*}

    だとします。左辺は xx だけの関数なので、もっとも簡単に xx と置きました。右辺は yy だけの関数なので、もっとも簡単に yy と置きました。つまり、それぞれは一定の値(定数 KK)ではありません。
     すると

    g′(x)g(x)=h′(y)h(y)\begin{align*} \dfrac{g'(x)}{g(x)}=\dfrac{h'(y)}{h(y)} \end{align*}

    のとき

    x=y\begin{align*} x=y \end{align*}

    となりますが、これだと xx と yy は互いに自由な値をとれません。例えば、
     x=1x=1 なら y=1y=1
     y=−3y=-3 なら x=−3x=-3
    と xx と yy は同じ値にならざるを得ず 、互いに無関係に、自由な値をとれないということです。
     一方、KK を定数として

    g′(x)g(x)=Kh′(y)h(y)=K\begin{align*} \dfrac{g'(x)}{g(x)}&=K\\ \dfrac{h'(y)}{h(y)}&=K\\ \end{align*}

    だと、任意の x, yx,  y に対して

    K=K\begin{align*} K=K \end{align*}

    と、恒等式(つねに成り立つ式)になります。例えば、
     x=1, y=3x=1,  y=3 のとき K=KK=K
     x=5, y=−2x=5,  y=-2 のとき K=KK=K
    と、xx と yy が互いに自由に値をとっても、つねに K=KK=K が成り立ちます。
     よって

    g′(x)g(x)=h′(y)h(y)\begin{align*} \dfrac{g'(x)}{g(x)}=\dfrac{h'(y)}{h(y)} \end{align*}

    より、KK を定数として

    g′(x)g(x)=Kh′(y)h(y)=K\begin{align*} \dfrac{g'(x)}{g(x)}&=K\\ \dfrac{h'(y)}{h(y)}&=K\\ \end{align*}

    が得られます。
     後は、前者から g(x)g(x) を求め、 後者から h(y)h(y) を求めれば、求める一般解 f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) が求まります。

     まずは、前者の微分方程式を解いて g(x)g(x) を求めましょう。

    (前者について)

    g′(x)g(x)=K\begin{align*} \dfrac{g'(x)}{g(x)}=K \end{align*}

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

    ∫g′(x)g(x) dx=∫K dx\begin{align*} \int\dfrac{g'(x)}{g(x)}\,dx=\int K\,dx \end{align*}

     C1C_1 を積分定数(任意の定数)として

    log⁡∣g(x)∣=Kx+C1\begin{align*} \log |g(x)|=Kx+C_1 \end{align*}

     左辺は積分公式

    ∫g′(x)g(x) dy=log⁡∣g(x)∣+C\begin{align*} \int \dfrac{g'(x)}{g(x)}  dy=\log |g(x)|+C \end{align*}

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

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

    より

    log⁡∣g(x)∣=Kx+C1 ⟺ ∣g(x)∣=eKx+C1\begin{align*} \log |g(x)|=Kx+C_1  \Longleftrightarrow |g(x)|=e^{Kx+C_1} \end{align*}

    と、∣g(x)∣|g(x)| について求めることができ

    ∣g(x)∣=eKx+C1=eKxeC1=eC1eKx\begin{align*} |g(x)|=e^{Kx+C_1}=e^{Kx}e^{C_1}=e^{C_1}e^{Kx} \end{align*}

     絶対値をはずして

    g(x)=±eC1eKx\begin{align*} g(x)=\pm e^{C_1}e^{Kx} \end{align*}

     定数部分 ±eC1\pm e^{C_1} を C2C_2 とおいて

    g(x)=C2eKx\begin{align*} g(x)=C_2e^{Kx} \end{align*}

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

     これで、前者から g(x)g(x) が求まりました。同様にして、後者から h(y)h(y) が求まります。

    (後者について)

    h′(y)h(y)=K\begin{align*} \dfrac{h'(y)}{h(y)}=K \end{align*}

    について、両辺を yy 軸方向に積分すると

    ∫h′(y)h(y) dy=∫K dy\begin{align*} \int\dfrac{h'(y)}{h(y)}\,dy=\int K\,dy \end{align*}

     C3C_3 を積分定数(任意の定数)として

    log⁡∣h(y)∣=Ky+C3\begin{align*} \log |h(y)|=Ky+C_3 \end{align*}

     対数と指数の関係より

    ∣h(y)∣=eKy+C3=eKyeC3=eC3eKy\begin{align*} |h(y)|=e^{Ky+C_3}=e^{Ky}e^{C_3}=e^{C_3}e^{Ky} \end{align*}

     絶対値をはずして

    h(y)=±eC3eKy\begin{align*} h(y)=\pm e^{C_3}e^{Ky} \end{align*}

     定数部分 ±eC3\pm e^{C_3} を C4C_4 とおいて

    h(y)=C4eKy\begin{align*} h(y)=C_4e^{Ky} \end{align*}

     以上により

    g(x)=C2eKxh(y)=C4eKy\begin{align*} g(x)&=C_2e^{Kx}\\ h(y)&=C_4e^{Ky} \end{align*}

     すると、偏微分方程式

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    の解を f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) と置いたので

    f(x, y)=g(x)h(y)=C2eKx⋅C4eKy=C2C4eKxeKy\begin{align*} f(x,  y)&=g(x)h(y)\\ &=C_2e^{Kx}\cdot C_4e^{Ky}\\ &=C_2C_4e^{Kx}e^{Ky} \end{align*}

     C2C4=CC_2C_4=C とおいて

    f(x, y)=CeKxeKy\begin{align*} f(x,  y)=Ce^{Kx}e^{Ky} \end{align*}

    (case2) g(x)h(y)=0g(x)h(y)=0 のとき

    g(x)h(y)=0↔g(x)=0 または h(y)=0\begin{align*} g(x)h(y)=0 \leftrightarrow g(x)=0  または  h(y)=0 \end{align*}

    となるが、f(x, y)=g(x)h(y)f(x,  y)=g(x)h(y) とおいたので
     g(x)=0g(x)=0 のとき f(x, y)=g(x)h(y)=0×h(y)=0f(x,  y)=g(x)h(y)=0\times h(y)=0
     h(y)=0h(y)=0 のとき f(x, y)=g(x)h(y)=g(x)×0=0f(x,  y)=g(x)h(y)=g(x)\times 0=0
    より

    f(x, y)=0\begin{align*} f(x,  y)=0 \end{align*}

     これは、 (case1) で求めた
     f(x, y)=CeKxeKyf(x,  y)=Ce^{Kx}e^{Ky}
    について、C=0C=0 とすれば
     f(x, y)=0⋅eKxeKy=0f(x,  y)=0\cdot e^{Kx}e^{Ky}=0
    より得ることができます。つまり f(x, y)=0f(x,  y)=0 の場合は、f(x, y)=CeKxeKyf(x,  y)=Ce^{Kx}e^{Ky} に含まれています。
     よって一般解は、(case1)、(case2) より、CC を任意の定数として

    f(x, y)=CeKxeKy (答)\begin{align*} f(x,  y)=Ce^{Kx}e^{Ky} \small{(答)} \end{align*}

    ■\blacksquare

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

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    の左辺を計算すると

    (左辺)=∂∂xf(x, y)=∂∂xCeKxeKy=CeKyddxeKx=CeKy⋅eKx外の微分⋅(Kx)′中の微分=CeKyeKx⋅K=CKeKxeKy\begin{align*} (左辺)&=\dfrac{\partial}{\partial x}f(x,  y)\\ &=\dfrac{\partial}{\partial x}Ce^{Kx}e^{Ky}\\ &=Ce^{Ky}\dfrac{d}{dx}e^{Kx}\\ &=Ce^{Ky}\cdot\underset{\tiny 外の微分}{e^{Kx}}\cdot\underset{\tiny 中の微分}{(Kx)'}\\ &=Ce^{Ky}e^{Kx}\cdot K\\ &=CKe^{Kx}e^{Ky}\\ \end{align*}

     右辺も同様に計算して

    (左辺)=∂∂yf(x, y)=∂∂yCeKxeKy=CeKxddyeKy=CeKx⋅eKy外の微分⋅(Ky)′中の微分=CeKxeKy⋅K=CKeKxeKy\begin{align*} (左辺)&=\dfrac{\partial}{\partial y}f(x,  y)\\ &=\dfrac{\partial}{\partial y}Ce^{Kx}e^{Ky}\\ &=Ce^{Kx}\dfrac{d}{dy}e^{Ky}\\ &=Ce^{Kx}\cdot\underset{\tiny 外の微分}{e^{Ky}}\cdot\underset{\tiny 中の微分}{(Ky)'}\\ &=Ce^{Kx}e^{Ky}\cdot K\\ &=CKe^{Kx}e^{Ky}\\ \end{align*}

    より、等しくなります。

    (2) f(0, 0)=1f(0,  0)=1 のとき
     一般解 f(x, y)=CeKxeKyf(x,  y)=Ce^{Kx}e^{Ky} に x=0, y=0x=0,  y=0 を代入すると、
    e0=1e^0=1 より

    f(0, 0)=Ce0e0=C⋅1⋅1=C=1\begin{align*} f(0,  0)=Ce^{0}e^{0}=C\cdot1\cdot1=C=1 \end{align*}

     f(1, 0)=ef(1,  0)=e のとき
     一般解 f(x, y)=CeKxeKyf(x,  y)=Ce^{Kx}e^{Ky} に x=1, y=0x=1,  y=0 を代入すると

    f(1, 0)=CeKe0=CeK⋅1=CeK=e\begin{align*} f(1,  0)=Ce^{K}e^{0}=Ce^{K}\cdot1=Ce^{K}=e \end{align*}

     これに、C=1C=1 を代入すると

    eK=e\begin{align*} e^{K}=e \end{align*}

    より k=1k=1。よって、C=1, k=1C=1,  k=1 を一般解

    f(x, y)=CeKxeKy\begin{align*} f(x,  y)=Ce^{Kx}e^{Ky} \end{align*}

    に代入して

    f(x, y)=1⋅e1⋅xe1⋅y\begin{align*} f(x,  y)=1\cdot e^{1\cdot x}e^{1\cdot y} \end{align*}

    より

    f(x, y)=exey (答)\begin{align*} f(x,  y)=e^{x}e^{y} \small{(答)} \end{align*}

    ■\blacksquare

     ex, eye^x,  e^y はそれぞれの変数で微分しても変わらないので、この解 f(x, y)=exeyf(x,  y)=e^{x}e^{y} は

    ∂∂xf(x, y)=∂∂yf(x, y)\begin{align*} \dfrac{\partial}{\partial x}f(x,  y)=\dfrac{\partial}{\partial y}f(x,  y) \end{align*}

    の解(の1つ)として予測がつくところです。

     物理の偏微分方程式

     物理の法則を表す方程式には、微分方程式で書かれる場合が多いです。その中で次の2つ
    ・ニュートンの運動方程式
    ・シュレーディンガー方程式
    を紹介します。高校物理の知識を多少使います。

    ニュートンの運動方程式

    F=m d2dt2x(t)\begin{align*} F=m\,\dfrac{d^2}{dt^2}x(t) \end{align*}

     F [N]F  [\text{N}]:物体に働く力
     m [kg]m  [\text{kg}]:物体の質量
     x(t) [m]x(t)  [\text{m}]:物体の位置(座標)
     t [s]t  [\text{s}]:時刻

     イングランドの数学者・物理学者、アイザック・ニュートンの発見した、物体の運動を記述する常微分方程式です。

    画像
    Isaac Newton(1642年-1727年)

     微分とは、何らかの量の微小変化を表します。
     たとえば、時刻 tt での物体の位置 x(t)x(t) の時間変化は速度です。
     つまり、x(t)x(t) を tt で微分したものは速度 v(t)v(t) になります。

    ddtx(t)=v(t)\begin{align*} \dfrac{d}{dt}x(t)=v(t) \end{align*}

     また、時刻 tt での物体の速度 v(t)v(t) の時間変化は加速度です。
     つまり、v(t)v(t) を tt で微分したものは加速度 a(t)a(t) になります。これは、x(t)x(t) を時刻 tt で2回微分したもの(時刻 tt での2階微分)と考えることができます。

    ddtv=ddt(ddtx(t))=d2dt2x(t)2階微分=a(t)\begin{align*} \dfrac{d}{dt}v=\dfrac{d}{dt}\left(\dfrac{d}{dt}x(t)\right)=\underset{2 階微分}{\dfrac{d^2}{dt^2}x(t)}=a(t) \end{align*}

    x(t)→微分v(t)→微分a(t)∣↑‾2階微分\begin{align*} \\[-13pt] &x(t) \xrightarrow[]{微分} v(t) \xrightarrow[]{微分} a(t)\\[-3pt] &\hspace{8pt}|\hspace{78pt}\uparrow\\[-14pt] &\hspace{9pt}\underline{\hspace{85.5pt}}\\[-4pt] &\hspace{35pt}{\footnotesize 2階微分} \end{align*}

     ニュートンの運動方程式は常微分方程式となります。
     高校では微分積分を使わないで物理を勉強するので、運動方程式は、加速度 d2dt2x(t)\dfrac{d^2}{dt^2}x(t) の部分を aa として

    F=ma\begin{align*} F=ma \end{align*}

    という形で習います。

    シュレーディンガー方程式

    iℏ∂∂tψ(x, t)=(−ℏ22m∂2∂x2+V(x))ψ(x,t)\begin{align*} i\hbar\frac{\partial}{\partial t}\psi(x,  t)=\left(-\dfrac{{\hbar}^2}{2m}\frac{{\partial}^2}{\partial x^2}+V(x)\right)\psi(x, t) \end{align*}

     オーストリアの物理学者、エルヴィン・シュレーディンガーが提案した、量子の状態を記述する偏微分方程式です。

    画像
    Erwin Schrödinger(1887年-1961年)

     ψ(x, t)\psi(x,  t) は状態関数とよばれ、電子や陽子のような小さな粒子の状態を表します。位置 xx と時間 tt の2変数関数です。分かりやすさのため、右辺を展開しましょう。

    iℏ∂∂tψ(x, t)=−ℏ22m∂2∂x2ψ(x,t)+V(x)ψ(x,t)\begin{align*} i\hbar\frac{\partial}{\partial t}\psi(x,  t)=-\dfrac{{\hbar}^2}{2m}\frac{{\partial}^2}{\partial x^2}\psi(x, t)+V(x)\psi(x, t) \end{align*}

     この方程式は2階偏微分方程式で、∂∂tψ(x, t)\dfrac{\partial}{\partial t}\psi(x,  t) は ψ(x, t)\psi(x,  t) の tt での偏微分を表し、∂2∂x2ψ(x,t)\dfrac{{\partial}^2}{\partial x^2}\psi(x, t) は ψ(x, t)\psi(x,  t) の xx での2階偏微分を表します。

    iℏ∂∂tψ(x, t)‾t で偏微分=−ℏ22m∂2∂x2ψ(x,t)‾x で2階偏微分+V(x)ψ(x,t)\begin{align*} i\hbar\underset{t  で偏微分}{\underline{\frac{\partial}{\partial t}\psi(x,  t)}}=-\dfrac{{\hbar}^2}{2m}\underset{x  で2階偏微分}{\underline{\frac{{\partial}^2}{\partial x^2}\psi(x, t)}}+V(x)\psi(x, t) \end{align*}

     他の部分については、ii は虚数単位といい、2乗すると −1-1 になる数です。

    i2=−1\begin{align*} i^2=-1 \end{align*}

     ℏ\hbar はディラック定数とよばれる定数で、”エイチバー” と読みます。プランク定数 hh を 2π2\pi で割った定数として定義されます。

    ℏ=h2π\begin{align*} \hbar=\dfrac{h}{2\pi} \end{align*}

     mm は粒子の質量を表します。電子なら電子の質量、陽子なら陽子の質量です。
     V(x)V(x) は粒子のもつポテンシャルエネルギーを表します。高校物理でいうところの位置エネルギーです。高い位置にある物体は重力による位置エネルギーをもち、ばねにつながれた物体は伸びたり縮んだりするばねの力、つまり弾性力による位置エネルギーをもちます。
     ポテンシャルエネルギー V(x)V(x) は、状況に応じたものを自由に与えることができます。その V(x)V(x) に対応する解 ψ(x, t)\psi(x,  t) を求めることによって、原子内の電子状態など、数多くの実験結果を説明することができます。

     偏微分が何なのかが分かると、この難しい数式も、何となくイメージすることができます。と同時に、この偏微分方程式は非常に難しい式であることがわかります。

    (再掲)シュレーディンガー方程式

    iℏ∂∂tψ(x, t)‾t で偏微分=−ℏ22m∂2∂x2ψ(x,t)‾x で2階偏微分+V(x)ψ(x,t)\begin{align*} i\hbar\underset{t  で偏微分}{\underline{\frac{\partial}{\partial t}\psi(x,  t)}}=-\dfrac{{\hbar}^2}{2m}\underset{x  で2階偏微分}{\underline{\frac{{\partial}^2}{\partial x^2}\psi(x, t)}}+V(x)\psi(x, t) \end{align*}

     まず、1階の偏微分と2階の偏微分が1つの式に含まれ、一方は tt での1階偏微分、他方は xx での2階偏微分と異なるものになります。しかも、この式には虚数 ii が含まれます。2乗してマイナスになる虚数は、小さな粒子の世界を説明する上で必要不可欠な数になります。
     そして、右辺を ψ(x, t)\psi(x,  t) でくくったのが、教科書などでよく見られる次の形です。

    iℏ∂∂tψ(x, t)=(−ℏ22m∂2∂x2+V(x))ψ(x,t)\begin{align*} i\hbar\frac{\partial}{\partial t}\psi(x,  t)=\left(-\dfrac{{\hbar}^2}{2m}\frac{{\partial}^2}{\partial x^2}+V(x)\right)\psi(x, t) \end{align*}

     右辺はカッコでくくられていますが、このように書くことによって、xx での2階偏微分、およびポテンシャルエネルギー V(x)V(x) との和を、 ψ(x, t)\psi(x,  t) に作用させるとみることができます。

     なぜこの式が生まれたのか、なぜ虚数が含まれているのか非常に興味深いですが、それを解説するのは本記事のレベルを超えるので、専門書の方をご覧ください。シュレーディンガー方程式の導出自体はそれほど難しいものではなく、高校数学(+α)の知識で理解できると思います。

    微分方程式の簡単な歴史

     微分積分学・力学の創始者として有名なアイザック・ニュートンは、力学の問題を解くために多くの常微分方程式を解いています。常微分方程式は17世紀から研究が始まり、多くの数学者、物理学者によって現在の形に整備されていきました
     その後18世紀には、スイスの数学者・天体物理学者、レオンハルト・オイラーが偏微分方程式の研究を行います。偏微分方程式を最初に研究したのはオイラーとされています。

    画像
    Leonhard Euler(1707年-1783年)

     この10年ほど後に、フランスの数学者・物理学者、ジャン・ル・ロン・ダランベールは、弦振動の運動方程式を偏微分方程式(ダランベールの式)として定式化をし、微分方程式は数学や物理学に限らず、工学、経済学、気象学、生物学など、多くの分野で利用されるようになりました。

    画像
    Jean Le Rond d'Alembert(1717年-1783年)

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

    あなたへのおすすめ