全微分の意味と計算方法

こんにちはコーヤです。

このページでは、多変数関数に含まれる変数すべてに注目して微分する、全微分の意味と計算方法を勉強します。

偏微分と全微分の違い

偏微分と比べながら全微分の意味を確認します。

左側が偏微分のイメージ、右側が全微分のイメージです。

偏微分の変化量

xΔxだけ変化した時のzの変化量Δzは以下の式で表されます。

Δz=f(x+Δx,y)f(x,y)

全微分の変化量

xΔxだけ、yΔyだけ変化した時のzの変化量Δzは以下の式で表されます。

Δz=f(x+Δx,y+Δy)f(x,y)

偏微分は1つの変数に注目して微分しましたが、全微分では全ての変数に注目して微分します。

全微分の条件と公式

全微分可能となる条件は以下の式が成り立つことです。

lim(ΔxΔy)(0,0)f(x+Δx,y+Δy){f(x,y)+fx(x,y)Δx+fy(x,y)Δy}(Δx)2+(Δy)2=0

このとき全微分は以下のように表されます。

df=fx(x,y)dx+fy(x,y)dy

条件と公式だけ見ても意味が分からないので、この条件と公式の導出を行います。

全微分の導出

曲面z=f(x,y)を平面で近似できるか計算します。

図の曲面がf(x,y)、平面αが点Aにおける近似平面を表しています。

A0,B0,C0,D0xy平面上の点で、それぞれ座標は以下のようにします。

A0=(x,y,0)B0=(x+Δx,y,0)C0=(x,y+Δy,0)D0=(x+Δx,y+Δy,0)

A0B0間の距離はΔxで、点A0C0間の距離はΔyです。

曲面f(x,y)は点ABCDで表し、近似平面αは点ABCDで表します。

ここで近似平面αの点Bと点Cの位置は、曲線ABの点Aでの接線が直線ABとなるように、曲線ACの点Aでの接線が直線ACとなるようにします。

次に、この3次元の図を3方向から見ていきます。ABA0B0の面、ACA0C0の面、ADA0D0の面の3種類を以下に表します。

図左側のABA0B0の面について考えます。

新たな点B1を図のように設定します。

A0B0間の距離はΔxなので、点AB1間の距離もΔxです。

直線ABの傾きは曲線ABの接線なので傾きfx(x,y)です。

したがって点BB1間の距離bは、傾きfx(x,y)の直線がΔxだけ変化したときの変化量なので

b=fx(x,y)Δx

と表すことができます。

図中央のACA0C0の面についても同様の計算です。点CC1間の距離cは以下のように表されます。

c=fy(x,y)Δy

最後に図右側のADA0D0の面です。

DD1間の距離dABCDが平面であることを利用して

d=b+c=fx(x,y)Δx+fy(x,y)Δy

となります。

次に点DD間の距離をε(x,y)とします。

Dz座標はf(x+Δx,y+Δy)で、点Dz座標は点Az座標f(x,y)に距離dを足したものです。

したがってε(x,y)

ε(x,y)=f(x+Δx,y+Δy){f(x,y)+fxΔx+fyΔy}

となります。

これで曲面f(x,y)と平面αの位置関係が分かりました。

曲面f(x,y)を平面αで近似する計算をしましたが、点DD間の距離ε(x,y)が曲面と近似平面の誤差となっています。

この誤差ε(x,y)が無視できるほど小さければ近似成功です。

ここで(Δx,Δy)(0,0)のとき、点A0D0間も0になります。点A0D0間が0になるより速くε(x,y)が0になってくれれば、誤差は無視できると考えられます。

この条件を式で表すと

lim(ΔxΔy)(0,0)ε(x,y)(Δx)2+(Δy)2=0

となります。これが全微分可能となる条件です。

この条件を満たす場合の点AD間のz座標の差Δzを計算しましょう。今Δz

Δz=d+ε(x,y)=fx(x,y)Δx+fy(x,y)Δy+ε(x,y)

と表すことができます。両辺を(Δx)2+(Δy)2で割って

Δz(Δx)2+(Δy)2=fx(x,y)Δx(Δx)2+(Δy)2+fy(x,y)Δy(Δx)2+(Δy)2+ε(x,y)(Δx)2+(Δy)2

ここで(Δx,Δy)(0,0)とすると、条件よりε(x,y)の項は無視できるので

dz(dx)2+(dy)2=fx(x,y)dx(dx)2+(dy)2+fy(x,y)dy(dx)2+(dy)2

両辺(dx)2+(dy)2をかけて

dz=fx(x,y)dx+fy(x,y)dy

となります。これで条件と公式の導出が終わりました。

全微分計算の例題

それでは例題2つで全微分の計算をしてみましょう。

全微分可能かどうか判定し、その後全微分を求めます。

例題1

f(x,y)=x2+y2

全微分可能かどうか確かめます。

fx(x,y)=2xfy(x,y)=2y

より誤差ε(x,y)

ε(x,y)=f(x+Δx,y+Δy){f(x,y)+fx(x,y)Δx+fy(x,y)Δy}=(x+Δx)2+(y+Δy)2{x2+y2+2xΔx+2yΔy}=(Δx)2+(Δy)2

となります。これを全微分の条件式に代入して

lim(ΔxΔy)(0,0)ε(x,y)(Δx)2+(Δy)2=lim(ΔxΔy)(0,0)(Δx)2+(Δy)2(Δx)2+(Δy)2=lim(ΔxΔy)(0,0)(Δx)2+(Δy)2=0

となるのでf(x,y)は全微分可能です。よって全微分は

df=fx(x,y)dx+fy(x,y)dy=2xdx+2ydy

となります。

例題2

f(x,y)=xy

全微分可能かどうか確かめます。

fx(x,y)=yfy(x,y)=x

より誤差ε(x,y)

ε(x,y)=f(x+Δx,y+Δy){f(x,y)+fx(x,y)Δx+fy(x,y)Δy}=(x+Δx)(y+Δy){xy+yΔx+xΔy}=ΔxΔy

となります。これを全微分の条件式に代入して

lim(ΔxΔy)(0,0)ε(x,y)(Δx)2+(Δy)2=lim(ΔxΔy)(0,0)ΔxΔy(Δx)2+(Δy)2

このままだと計算できないので

Δx=rcosθΔy=rsinθ

としてr0の極限に変換します。

lim(ΔxΔy)(0,0)ΔxΔy(Δx)2+(Δy)2=limr0r2cosθsinθr=limr0rcosθsinθ=0

となるのでf(x,y)は全微分可能です。よって全微分は

df=fx(x,y)dx+fy(x,y)dy=ydx+xdy

となります。

全微分の変数変換

関数f(x,y)xyが、他の関数で表されている場合の微分を考えます。

1変数で置換する場合

変数tを用いてx,y

x=x(t)y=y(t)

と表される場合、関数f(x,y)tで微分すると

dfdt=fx(x,y)dx+fy(x,y)dydt=fx(x,y)dxdt+fy(x,y)dydt=fxdxdt+fydydt

となります。

例えば

f(x,y)=xyx=t3y=sint

この関数をtで微分すると

dfdt=fxdxdt+fydydt=y3t2+xcost=sint3t2+t3cost=3t2sint+t3cost

となります。

2変数で置換する場合

変数u,vを用いてx,y

x=x(u,v)y=y(u,v)

と表される場合、関数f(x,y)uvでそれぞれ微分すると

fu=fx(x,y)dx+fy(x,y)dyu=fx(x,y)xu+fy(x,y)yu=fxxu+fyyu

fv=fx(x,y)dx+fy(x,y)dyv=fx(x,y)xv+fy(x,y)yv=fxxv+fyyv

となります。

例えば

f(x,y)=xyx=u+vy=uv

この関数をuvでそれぞれ微分すると

fu=fxxu+fyyu=y1+xv=uv1+(u+v)v=2uv+v2

fv=fxxv+fyyv=y1+xu=uv1+(u+v)u=2uv+u2

となります。

まとめ

全微分は多変数関数に含まれる全ての変数に注目して微分します。

全微分可能となる条件は以下の式が成り立つことです。

lim(ΔxΔy)(0,0)f(x+Δx,y+Δy){f(x,y)+fxΔx+fyΔy}(Δx)2+(Δy)2=0

このとき全微分は以下のように表されます。

df=fxdx+fydy

コメント欄