Revised at 22:13, July 4, 2016
解析学A 第12
回http://my.reset.jp/˜gok/math/ 1
12
合成関数の偏導関数12.1
関数の合成とその微分2変数関数
f (x, y)
と2つの1変数関数g(t), h(t)
があった時に、x = g(t), y = h(t)
と云う関係式を接着剤にして新しい関数k(t)
を次の様に作る事が出来ます:
f (x, y) g(t) h(t)
+
のり
x = g(t)
y = h(t) −→ k(t) = f (g(t), h(t)) .
もちろん、
(g(t), h(t))
がf
の定義域に入る様なt
の範囲に限られますが、まあ、合成出 来るかどうか気にするよりも、合成出来たとして、更にf, g, h
が(偏)微分可能であっ たならk
も微分可能になる様に思われますので導関数がどうなるか調べてみましょう。k(t + p) − k(t) p
= f (g(t + p), h(t + p)) − f (g(t), h(t)) p
分子にダミー項を挿入して変形すると
= f (g(t + p), h(t + p)) − f (g(t), h(t + p))
p + f (g(t), h(t + p)) − f (g(t), h(t)) p
= f (g(t + p), h(t + p)) − f (g(t), h(t + p))
g(t + p) − g(t) · g(t + p) − g(t) p + f (g(t), h(t + p)) − f (g(t), h(t))
h(t + p) − h(t) · h(t + p) − h(t) p
となりますが、平均値の定理によれば= f x (c, h(t + p)) g(t + p) − g(t)
p + f y (g(t), d) h(t + p) − h(t) p
となる様な
c
がg(t)
とg(t + p)
の間に、d
がh(t)
とh(t + p)
の間にそれぞれ存在しま す。従って、g, h
が微分可能、f
が偏微分可能であって、更にf x , f y
が連続関数である ならばここでp → 0
の極限をとれば→ f x (g(t), h(t))g 0 (t) + f y (g(t), h(t))h 0 (t)
が成り立つ事が分かります。従って
k 0 (t) = f x (g(t), h(t))g 0 (t) + f y (g(t), h(t))h 0 (t)
となりますが、直感的に言ってf = k, x = g, y = h
なわけですからf 0 (t) = f x (x(t), y(t))x 0 (t) + f y (x(t), y(t))y 0 (t)
と書いても悪くないでしょう。ただし、
f
は本来2変数関数ですからやっぱりf 0 (t)
って えのは正確にはまずいですよ。そこら辺きちんと理解出来てるならこうした書き方をし ても悪くはないと言うことです。あるいは、別の記法ではdk
dt (t) = @f
@x (g(t), h(t)) dg
dt (t) + @f
@y (g(t), h(t)) dh dt (t)
としたり、これを簡略化してdf dt = @f
@x dx dt + @f
@y dy dt
と書く事もあります。12.2
具体的な計算例今見た一般的な結果を利用して
f (x, y) = x 2 − y 2 , g(t) = e t + e − t , h(t) = e t − e − t
のときに合成関数F (t) = f (g(t), h(t))
の微分を計算してみましょう。f x (x, y) = 2x, f y (x, y) = − 2y, g 0 (t) = e t − e − t , h 0 (t) = e t + e − t
ですからF 0 (t) = 2x(t)(e t − e − t ) − 2y(t)(e t + e − t ) = 2(e 2t − e − 2t ) − 2(e 2t − e − 2t ) = 0
と云う風に計算する事が出来ますが、もちろんこのような具体的な場合にはそんな事を しなくても実際に合成した関数がどうなるのかを書き下してやって:F (t) = (e t + e − t ) 2 − (e t − e − t ) 2 = 4
その1変数関数を普通に微分すればオッケーです(そりゃあ
0
になりますわね)。じゃあ、さっきやった様な一般的な計算は意味ないじゃん、と云う事になるのかと言 うと、そうではありません。あとのセクションで考える様な事柄(方向微分など)を理 解するには是非必要です。
Revised at 22:13, July 4, 2016
解析学A 第12
回http://my.reset.jp/˜gok/math/ 2
12.3
2変数の合成関数2変数関数
f (x, y)
と別の2つの2変数関数g(v, w), h(v, w)
があった時に、x = g(v, w), y = h(v, w)
と云う関係式を使って新しい関数k(v, w)
を次の様に作る事が出来 ます:
f (x, y) g(v, w) h(v, w)
+
のり
x = g(v, w)
y = h(v, w) −→ k(v, w) = f(g(v, w), h(v, w)) .
この合成関数のv, w
に関する偏導関数はどうやって計算されるでしょうか。さっきやった合成は中に入れる方の関数が1変数でしたが今回は2変数です。しか し、どうせ偏微分するわけですし、その際もう一方の変数は定数だと思ってあたかも1 変数関数であるかの様に扱って微分すれば良いのだと云う事から考えれば、中に入れる ものが2変数になったところで特に何も変わりません。
k(v, w)
をv
で偏微分する時は、2変数関数f (x, y)
と 1変数関数x = g(v, w), y =
h(v, w)
(w
は定数だと思っています)を合成した関数のv
に関する 普通の 微分ですから、さっきと同様に
k v (v, w) = f x (g(v, w), h(v, w))g v (v, w) + f y (g(v, w), h(v, w))h v (v, w)
とすれば良い事になります(もちろん各偏導関数が連続になっていないとダメです)。
これを象徴的に書けば(
k = f, g = x, h = y
などと同一視します)f v = f x x v + f y y v , @f
@v = @f
@x
@x
@v + @f
@y
@y
@v
であり、全く同様にw
に関する偏微分はk w (v, w) = f x (g(v, w), h(v, w))g w (v, w) + f y (g(v, w), h(v, w))h w (v, w)
あるいは象徴的に
f w = f x x w + f y y w , @f
@w = @f
@x
@x
@w + @f
@y
@y
@w
となります。もちろん、具体的な計算においては、まず合成を実行して
v, w
の関数として書き下 してからそれを実際に偏微分すれば良いでしょう。12.3.1
極座標による計算例2変数関数
f (x, y)
を、x = r cos θ, y = r sin θ
によってr, θ
の関数と見るとf r = f x x r + f y y r = f x cos θ + f y sin θ
f θ = f x x θ + f y y θ = − f x r sin θ + f y r cos θ
f rr = (f x ) r cos θ + (f y ) r sin θ
= (f xx x r + f xy y r ) cos θ + (f yx x r + f yy y r ) sin θ
= f xx cos 2 θ + f xy sin θ cos θ + f yx cos θ sin θ + f yy sin 2 θ f θθ = − (f x ) θ r sin θ − f x r cos θ + (f y ) θ r cos θ − f y r sin θ
= − (f xx x θ + f xy y θ )r sin θ − f x r cos θ + (f yx x θ + f yy y θ )r cos θ + f y r sin θ
= −{ f xx ( − r sin θ) + f xy r cos θ } r sin θ − f x r cos θ
+ { f yx ( − r sin θ) + f yy r cos θ } r cos θ + f y r sin θ
= f xx r 2 sin 2 θ − f xy r 2 cos θ sin θ − f x r cos θ
− f yx r 2 sin θ cos θ + f yy r 2 cos 2 θ + f y r sin θ 1
r 2 f θθ = f xx sin 2 θ − f xy cos θ sin θ − f x
1 r cos θ
− f yx sin θ cos θ + f yy cos 2 θ + f y
1 r sin θ
から、f xx + f yy = f rr + 1 r f r + 1
r 2 f θθ
が成り立つことが分かります(
f (x, y)
の偏導関数の連続性は仮定します)。12.4
方向微分2変数関数
f (x, y)
のグラフを考えます。これはつまりz = f (x, y)
が表す曲面ですね。そして
xy-
平面の(もちろん関数の定義域内の)点(a, b)
とこの点を通る直線L
を考 えます。直線L
は方程式で表現しても良いのですがここでは方向ヴェクターを√ v w
!
と してパラメータ表示しておきます。L :
√ x y
!
=
√ a b
! + t
√ v w
!
(
t
はパラメータ)Revised at 22:13, July 4, 2016
解析学A 第12
回http://my.reset.jp/˜gok/math/ 3
このような準備をした上で、2変数関数f (x, y)
をこの直線
L
の上だけで見てみるとどうなっているか 調べてみましょう。これは即ち直線
L
を含みxy-
平面に垂直な平面P
でz = f (x, y)
のグラフ(曲面です)を切った時の 切り口の曲線がグラフになっていますね(左図)。しかし、直線
L
はただの直線ですから座標軸では ありません。従って 目盛 も付いていません。そ こでこの直線に目盛を打って(座標軸と考えて)本 当にこの直線上の関数と考えられる様にするにはど うしたら良いでしょうか。要するにどんな目盛の付け方をすれば良いのかと云う事ですが、どんな風に目盛を 打ってもグラフの高さを変えているわけではありませんから、やはり目盛の間隔は
x-
軸の目盛の間隔と同じ1にしなければなりません。もし間隔2で目盛を打ってしまったら、底辺と高さの縮尺が変化する事になってしま い、おかしな事になってしまいます(直線
L
がたまたまx-
軸だった場合を考えてみれ ば良いでしょう)。それを実行するにはどうしたら良いかと云えば、パラメータ
t
が1増えるごとに直線 上では方向ヴェクター√ v w
!
の大きさ分だけ点が移動するわけですから、元々方向ヴェ クター
√ v w
!
を単位ヴェクターで取っておけばパラメータの変化1が丁度点の距離1の 変化に一致する事になります。
これを数式で表現すると、合成関数
F (t) = f (x(t), y(t)) = f (a + tv, b + tw)
をt
の 関数として考えると云う事になる筈です。例えばこの合成関数を
t
に関して(点t = 0
で)微分すれば、曲面z = f(x, y)
を平面P
で切った切り口の曲線の接戦の傾きが出る筈です。これを、関数f (x, y)
の点(a, b)
におけるヴェクター√ v w
!
方向の方向微分係数と言います。
もし
f (x, y)
が偏微分可能で更に偏導関数が連続であれば、先に見た合成関数の微分の様子から、この方向微分係数は
F 0 (0) = f x (a, b)v + f y (a, b)w =
√ f x (a, b) f y (a, b)
!
·
√ v w
!
となっています。当然ですが
√ 1 0
!
方向(
x-
軸方向)の方向微分はf x (a, b)
ですし、√ 0 1
!
方向(
y-
軸方向)の方向微分はf y (a, b)
になっていますね。従って
f x (a, b)
とf y (a, b)
と云うたった2方向での方向微分さえ求めておけば、後は 内積を取る事によって全ての方向での方向微分が求まってしまうと見る事が出来るで しょう。そう云う意味でたった2方向しかやっていない印象のある偏微分ではあります が、実はこれだけで十分多くの事が分かってしまうわけなんです。このヴェクター(関数)
√ f x (a, b) f y (a, b)
!
は今後もよく利用しますが、特別に名前がつけ られています。これを関数
f (x, y)
の(点(a, b)
に於ける)勾配(gradient
)と言い、記 号:gradf (a, b)
で表す事にします。12.5
参考:曲線に沿った微分この様に、直線に沿って
(x, y)
を(a, b)
に近づけた場合の関数f (x, y)
の値の変化率 は方向微分係数(これは特殊な場合として偏微分係数を含んでいます)で計算される事 が分かりましたが、以前2変数関数の極限値を計算した時に見た様に平面内での点の近 づけ方は他にも沢山あった筈です。そこで曲線に沿って(x, y)
を(a, b)
に近づけた場合 の変化率はどうやって計算されるか考えてみましょう。下図の様に
xy-
平面内の点(a, b)
を通る曲線C : (t) =
√ x(t) y(t)
!
を考えて、この曲線 上での関数
f (x, y)
の値f (x(t), y(t))
の変化を見ます。曲線
C
をz-
軸方向に伸ばして得られる曲面Q
と2変数関数z = f (x, y)
のグラフの 交線が曲線C
上でのf (x, y)
のグラフになっています。本来、点
(x, y)
の変化量に対する値f (x, y)
の変化量の比を見ようようとしているわけで すが、ここで1つ問題があります。同じ(軌跡 をもつ)曲線でも、内部のパラメータの取り方 は様々です(同じ道を走っていても車の速度は 色々です)。Revised at 22:13, July 4, 2016
解析学A 第12
回http://my.reset.jp/˜gok/math/ 4
だからパラメーターt
の変化に対する値f(x(t), y(t))
の変化量の比で計算しようとするならば違うパラメーター(=違う速度)で計算すると当然
f (x(t), y(t))
の値の変化率 も変わって来てしまいます。そう云う意味で、パラメーターの変化に対する点の位置の 変化が常に1である様な状況で考えなければならない事が分かります。これは即ち速度 ヴェクターの大きさが常に1である様な等速運動を考える事ですが、パラメータの言葉 で云えば弧長パラメータと云う事になります。実際、点
(a, b)
(この点はパラメータt = t 0
とします)からパラメーターの正の向きの弧長を
s(t)
とすれば、s(t) = Z t
t
0| 0 (p) | dp, s 0 (t) = | 0 (t) |
ですので、この逆関数
t = t(s)
を考えればs = s(t(s))
の両辺を微分して1 = s 0 (t(s))t 0 (s), t 0 (s) = 1
s 0 (t(s)) = 1
| 0 (t(s)) |
ですから0 (s) = 0 (t(s))t 0 (s) = 0 (t(s))
| 0 (t(s)) |
となって、弧長をパラメータに取ると速度ヴェクターは常に単位ヴェクターである事が 分かります。
そこで最初からパラメータ
t
は弧長で取ってあると仮定すれば、F (t) = f (x(t), y(t))
に対して、F 0 (t) = f x (x(t), y(t))x 0 (t) + f y (x(t), y(t))y 0 (t) = gradf(x(t), y(t)) · 0 (t)
が曲線C
に沿った変化率を表していますが、0 (t)
は単位ヴェクターでしたので結局こ れは先に見た0 (t)
方向の方向微分係数に一致しています。この様に、曲線に沿って近づけた時の変化率を見ても、それは結局接線方向の方向微 分に他ならない事が分かります。しかも方向微分と云うもの自体、偏微分さえ計算して やれば後は内積でオッケーでしたね。
その意味で、偏微分さえやってしまえば、全ての近づけ方を制覇した事になると言え るでしょうか。
12.6
参考:全微分と偏微分の関係ここで『全微分可能性』の定義を振り返ると、適当な定数
p, q
が存在して(v,w) lim → (a,b)
f (v, w) − f (a, b) − p(v − a) + q(w − b) p (v − a) 2 + (w − b) 2 = 0
となるとき、
f (x, y)
は(a, b)
において全微分可能であると定義したわけでしたが、もし もf (x, y)
が点(a, b)
で全微分可能であるならば、この極限計算においてw = b, v → a
の場合を考えればv lim → a
f (v, b) − f (a, b) − p(v − a)
| v − a | = 0
であり、極限値が0
なので特にv lim → a
f (v, b) − f (a, b) − p(v − a)
v − a = 0
であることも分かりますから、
f (v, b) − f (a, b)
v − a = f (v, b) − f (a, b) − p(v − a) + p(v − a) v − a
= f (v, b) − f (a, b) − p(v − a)
v − a + p
→ p (as v → a)
となり、実は
f (x, y)
は点(a, b)
において変数x
に関して偏微分可能であって、p = @f @x (a, b)
である事が分かります。全く同様にq = @f @y (a, b)
でもありますから、全微分の定義に現 れた 微分係数p, q
は実は偏微分係数だった事が分かります。この様に偏微分は関数を局所的に1次関数で近似する考え方とも通じており、いずれ 明らかになるように、関数の局所的な性質を解析する目的において十分大きな役割を果 たす事が出来ます。
Exercise
基本演習