======================== 合成関数の微分への利用 ======================== 任意の関数を偏微分しようと思ったら、関数の形だけ微分の公式を見つけなければなりません。 その公式を見つける際に、もとの関数を合成関数として見る方法があります。 そうすることで関数の偏微分するという問題は、合成関数の偏微分の問題に変わります。 この方法の利点は、ごく限られた微分の公式を覚えているだけで、多くの偏微分の公式を見つけることが可能になるのです。 ここでは、掲示板でよせられた質問と電磁気学への応用についてのみ考えていきます。一般的な議論についてはここでは扱いませんので、 あらかじめご了承下さい。 質問と解説 -------------------------------- 掲示板で次のような問題について質問がよせられてきました。 質問 ^^^^^^^^^^^ 位置ベクトル ${\bm{r} = x\bm{e}_{x} + y \bm{e}_{y} + z \bm{e}_{z}}$ の動径距離 ${|\bm{r}| = r = \sqrt{x^2 +y^2 +z^2}}$ を直交成分 ${x,y,z}$ で偏微分すると、どうして次のようになるのですか。 \nabla r &= \left( \frac{\partial }{\partial x} \bm{e}_{x} + \frac{\partial }{\partial y} \bm{e}_{y} + \frac{\partial }{\partial z} \bm{e}_{z} \right) r\\ &= \frac{\partial r}{\partial x} \bm{e}_{x} + \frac{\partial r}{\partial y} \bm{e}_{y} + \frac{\partial r}{\partial z} \bm{e}_{z}\\ &= \frac{x}{r} \bm{e}_{x} + \frac{y}{r} \bm{e}_{y} + \frac{z}{r} \bm{e}_{z}\\ &= \frac{x\bm{e}_{x} + y \bm{e}_{y} + z \bm{e}_{z}}{r}\\ &= \frac{\bm{r}}{r} \tag{1} 解説 ^^^^^^^^^^^ この問題の場合、関数 $r$ を次のような合成関数として見るのが便利です。 &f = f(x,y,z) \tag{2}\\ &r = r(x,y,z) = r(f) \tag{3} すると偏微分は合成関数の偏微分の公式 \frac{\partial r}{\partial x} = \frac{\partial r}{\partial f} \frac{\partial f}{\partial x} \tag{4} が使えます。 $\tag{2}$ , $\tag{3}$ に具体的な関数を書き込むと、次のようになります。 &f = x^2 +y^2 +z^2 \tag{5}\\ &r = \sqrt{x^2 +y^2 +z^2 } = \sqrt{f} \tag{6} これでもう計算の準備が整いました。ここでひとこと説明を加えておきまと、 結局ここで説明している計算の方法としては ${\tag{5}}$ , $\tag{6}$ のように、適当に微分の公式が使えるような 関数のかたちとして見て、それらを $\tag{4}$ 式に代入して計算をおしすすめていこうという方針です。数式で示すと、次の通りです。 &\frac{\partial f}{\partial x} = 2x \tag{7}\\ &\frac{\partial r}{\partial f} = \frac{1}{2} f^{\frac{1}{2} - 1} =\frac{1}{2} f^{-\frac{1}{2}} =\frac{1}{2 \sqrt{x^2 + y^2 +z^2}}\tag{8}\\ &\frac{\partial r}{\partial x} = 2x \frac{1}{2 \sqrt{x^2 + y^2 +z^2}} = \frac{x}{\sqrt{x^2 + y^2 +z^2}} = \frac{x}{r} \tag{9} ただし $\tag{7}$ と $\tag{8}$ で、次の公式の $n=2$ の場合と $n=1/2$ の場合を使っています。( $g:$ 関数) \frac{\partial g^{n}}{\partial g} = n g^{n-1} (n \neq 0 ) \tag{10} 後の $y,z$ についての偏微分についても同じ方法でもとまります。結果だけ書いておくと &\frac{\partial r}{\partial y} = \frac{y}{r} \tag{11}\\ &\frac{\partial r}{\partial y} = \frac{y}{r} \tag{12} です。こうして $\tag{9}$ , $\tag{11}$ , $\tag{12}$ から $\tag{1}$ が得られる事が示されました。 補足説明 ^^^^^^^^^^ ここでは与えた式 $\tag{4}$ の導出を説明しておきます。 実はこの式は、全微分 $df,dr$ から自然に出てくる式です。これらの全微分の具体的なかたちは $\tag{2}$ , $\tag{3}$ より &df = \frac{\partial f}{\partial x}dx + \frac{\partial f}{\partial y}dy + \frac{\partial f}{\partial z}dz \tag{13}\\ &dr = \frac{\partial r}{\partial x}dx + \frac{\partial r}{\partial y}dy + \frac{\partial r}{\partial z}dz = \frac{\partial r}{\partial f}df \tag{14} になります。 $\tag{13}$ を $\tag{14}$ の最後のところに代入することによって \frac{\partial r}{\partial x}dx + \frac{\partial r}{\partial y}dy + \frac{\partial r}{\partial z}dz = \frac{\partial r}{\partial f}\frac{\partial f}{\partial x}dx + \frac{\partial r}{\partial f}\frac{\partial f}{\partial y}dy + \frac{\partial r}{\partial f}\frac{\partial f}{\partial z}dz が成り立つことが分かります。両辺を比較すると公式として与えた式 $\tag{4}$ が出てくるというわけです。 &\frac{\partial r}{\partial x} = \frac{\partial r}{\partial f} \frac{\partial f}{\partial x} \\ &\frac{\partial r}{\partial y} = \frac{\partial r}{\partial f} \frac{\partial f}{\partial y} \tag{16}\\ &\frac{\partial r}{\partial z} = \frac{\partial r}{\partial f} \frac{\partial f}{\partial z} \tag{17} 重要事項(まとめ) ^^^^^^^^^^^^^^^^^ 関数を何の関数として見るかが、重要。 練習問題 ---------- ここに載せている問題は主に電磁気学でする計算のものです。答えもできるだけ丁寧に書いたので、ぜひ挑戦してみてください。 練習問題1. ^^^^^^^^^^^ 次の式を $\tag{1}$ , $\tag{4}$ , $\tag{10}$ を使って ${\bm{e}_x}$ , ${\bm{e}_y}$ , ${\bm{e}_z}$ の各成分を求めてください。 &\nabla \bm{r} \tag{18}\\ &\nabla \frac{1}{r} \tag{19} (答え) \nabla \bm{r} &= \left( \frac{\partial }{\partial x} \bm{e}_{x} + \frac{\partial }{\partial y} \bm{e}_{y} + \frac{\partial }{\partial z} \bm{e}_{z} \right) \left(x \bm{e}_{x} + y \bm{e}_{y} + z \bm{e}_{z} \right)\\ &= \left( \frac{\partial x}{\partial x} \bm{e}_{x}\bm{e}_{x} +\frac{\partial y}{\partial x} \bm{e}_{x}\bm{e}_{y} +\frac{\partial z}{\partial x} \bm{e}_{x}\bm{e}_{z} +\frac{\partial x}{\partial y} \bm{e}_{y}\bm{e}_{x} +\frac{\partial y}{\partial y} \bm{e}_{y}\bm{e}_{y} +\frac{\partial z}{\partial y} \bm{e}_{y}\bm{e}_{z} +\frac{\partial x}{\partial z} \bm{e}_{z}\bm{e}_{x} +\frac{\partial y}{\partial z} \bm{e}_{z}\bm{e}_{y} +\frac{\partial z}{\partial z} \bm{e}_{z}\bm{e}_{z} \right)\\ &= \bm{e}_{x}\bm{e}_{x} + \bm{e}_{y}\bm{e}_{y} + \bm{e}_{z}\bm{e}_{z} (=1) \tag{20} \\ \\ \nabla \frac{1}{r} &= \nabla r^{-1}\\ &= \left( \frac{\partial }{\partial x} \bm{e}_{x} + \frac{\partial }{\partial y} \bm{e}_{y} + \frac{\partial }{\partial z} \bm{e}_{z} \right) \frac{1}{r} \\ &= \frac{\partial r^{-1}}{\partial x}\bm{e}_{x} + \frac{\partial r^{-1}}{\partial y }\bm{e}_y + \frac{\partial r^{-1}}{\partial z} \bm{e}_{z} \\ &= \frac{\partial r}{\partial x} \frac{\partial r^{-1}}{\partial r}\bm{e}_{x} + \frac{\partial r}{\partial y} \frac{\partial r^{-1}}{\partial r}\bm{e}_{y} + \frac{\partial r}{\partial z} \frac{\partial r^{-1}}{\partial r}\bm{e}_{z} [ \because (4)] \\ &= \frac{\partial r^{-1}}{\partial r} \left( \frac{\partial r}{\partial x}\bm{e}_x + \frac{\partial r}{\partial y}\bm{e}_y + \frac{\partial r}{\partial z}\bm{e}_z \right) \\ &= \frac{\partial r^{-1}}{\partial r} \left( \frac{\partial }{\partial x}\bm{e}_x + \frac{\partial }{\partial y}\bm{e}_y + \frac{\partial }{\partial z} \bm{e}_z \right) r \\ &= \frac{\partial r^{-1}}{\partial r} \nabla r \\ &= (-1) r^{(-1-1)} \nabla r [\because (10)] \\ &= - \frac{1}{r^2} \frac{\bm{r}}{r} [\because (1)] \\ &= - \frac{\bm{r}}{r^{3}} \\ &= \frac{-x}{r^{3}}\bm{e}_{x} + \frac{-y}{r^{3}}\bm{e}_{y} + \frac{-z}{r^{3}}\bm{e}_{z} \tag{21} 式 $\tag{20}$ の最後のカッコ書きの部分については演算規則が単位元を満たしている事を確認してみてください。 あと、式 $\tag{19}$ のように微分される関数がある1つの変数 $s$ の関数のときは、次のように $\nabla$ を書き換えると複雑な計算をまとめるときに便利です。 \nabla &= \frac{\partial}{\partial x} \bm{e}_{x} + \frac{\partial}{\partial x} \bm{e}_{x} + \frac{\partial}{\partial x} \bm{e}_{x}\\ &= \left( \frac{\partial s}{\partial x} \bm{e}_{x} + \frac{\partial s}{\partial y} \bm{e}_{y} + \frac{\partial s}{\partial z} \bm{e}_{z} \right) \frac{\partial}{\partial s} \tag{22} 例として $s=r$ のときの結果を書いておきます。これも $\tag{1}$ を使えば、簡単に分かるので確認してください。 \nabla = \frac{\bm{r}}{r} \frac{\partial}{\partial r} \tag{23} 練習問題2. ^^^^^^^^^^^ 真空中で位置 ${\bm{r}_1}$ に電荷が ${-q}$ 、位置 $\bm{r}$ に電荷が ${q}$ があるとします。 ${\left( \bm{r}_2 - \bm{r}_1 = \bm{l} \right)}$ で一定 のとき、位置 ${\bm{r}}$ における電場を ${\left( r \gg l \right)}$ の範囲での近似解をもとめてください(電場 $\bm{E}$ と $V$ の関係は ${\bm{E} = - \nabla V}$ の関係で結ばれます)。 このときの電位 ${V(\bm{r})}$ は近似的に次のように書けます。 V(\bm{r}) &\simeq \frac{q \bm{l}\cdot \bm{r}}{4 \pi \epsilon_{0} r^{3}} \tag{24-1}\\ &= - \frac{q}{4 \pi \epsilon_{0}}\bm{l} \cdot \nabla \frac{1}{r} [\because (19)] \tag{24-2} (答え) \bm{E}(\bm{r}) &= -\nabla V(\bm{r})\\ &\simeq \frac{1}{4 \pi \epsilon_{0}} q\bm{l} \cdot \nabla \nabla \frac{1}{r} [\because (24-2)]\\ &= - \nabla \frac{q \bm{l}\cdot \bm{r}}{4 \pi \epsilon_{0} r^{3}} [\because (24-1)]\\ &= - \frac{1}{4 \pi \epsilon_{0} } q\bm{l} \cdot \nabla (\bm{r}r^{-3}) \\ &= \frac{1}{4 \pi \epsilon_{0}} q\bm{l} \cdot ( - \bm{r} \nabla r^{-3} - r^{-3} \nabla \bm{r} ) [\because (4),(20)]\\ &= \frac{1}{4 \pi \epsilon_{0}} q \bm{l} \cdot \left( - \bm{r} \frac{\partial r^{-3}}{\partial r}\nabla r - r^{-3} \right) \\ &= \frac{1}{4 \pi \epsilon_{0}} q \bm{l} \cdot \left( - \bm{r} (-3)r^{-3-1} \nabla r - r^{-3} \right) [\because (10)] \\ &= \frac{1}{4 \pi \epsilon_{0}} q \bm{l} \cdot \left( 3 \bm{r}r^{-4} \bm{r}r^{-1} - r^{-3} \right) [\because (1)]\\ &= \frac{1}{4 \pi \epsilon_{0}} q \bm{l} \cdot \left( 3 \frac{\bm{r}\bm{r}}{r^{5}} - \frac{1}{r^{3}} \right) \\ &= \frac{1}{4 \pi \epsilon_{0}} \left( 3 \frac{q\bm{l}\cdot\bm{r}\bm{r}}{r^{5}} - \frac{q\bm{l}}{r^{3}} \right) \tag{25} また、この式の上から2段目と最後から2番目の段を見比べると次の事も分かります。 \nabla \nabla \frac{1}{r} = 3 \frac{\bm{r}\bm{r}}{r^{5}} - \frac{1}{r^{3}} \tag{26} 補足 ^^^^^^^ 練習問題2.で与えた式 $\tag{21}$ について説明を加えておきます。 まずは、クーロンの法則より電荷 ${-q}$ がつくる電位 ${V_{1}}$ , 電荷 $q$ がつくる電位 ${V_{2}}$ は &V_{1}(\bm{r},\bm{r}_{1}) = \frac{-q}{4 \pi \epsilon_{0}} \frac{1}{|\bm{r} - \bm{r}_{1}|}\\ &V_{2}(\bm{r},\bm{r}_{2}) = \frac{q}{4 \pi \epsilon_{0}} \frac{1}{|\bm{r} - \bm{r}_{2}|} と書けます。そこから全電位は V(\bm{r},\bm{r}_{1},\bm{r}_{2}) &= V_{1}(\bm{r},\bm{r}_{1}) + V_{2}(\bm{r},\bm{r}_{2})\\ &= \frac{1}{4 \pi \epsilon_{0}} \left( \frac{-q}{|\bm{r} - \bm{r}_{1}|} + \frac{q}{|\bm{r} - \bm{r}_{2}|} \right) \\ &= \frac{q}{4 \pi \epsilon_{0}} \left( -|\bm{r} - \bm{r}_{1}|^{-1} + |\bm{r} - \bm{r}_{2}|^{-1} \right) \tag{27} と書けることが分かります。ここで次の式をそれぞれ ${\frac{r_{1}}{r} = X}$ , ${\frac{r_{2}}{r} = Y}$ の関数として見ます。 &\left( 1 - 2 \frac{{r}_{1}}{r} \cos {\theta}_{1} + \left( \frac{r_{1}}{r} \right)^2 \right)^{-\frac{1}{2}} = F \left( \frac{r_{1}}{r} \right) = F(X) \tag{28-1}\\ &\left( 1 - 2 \frac{{r}_{2}}{r} \cos {\theta}_{2} + \left( \frac{r_{2}}{r} \right)^2 \right)^{-\frac{1}{2}} = G \left( \frac{r_{2}}{r} \right) = G(Y) \tag{28-2} 次にこの二つの式を原点を中心に テイラー級数_ に展開します。すると F(X) &= \sum_{n=0}^{\infty} \frac{1}{n!} \left. \frac{\partial^{n}F}{\partial X^{n}} \right|_{X=0} X^{n} \\ &= \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} \cos^{(n-2k)} \theta_{1} \left( \frac{r_{1}}{r} \right)^{n}\\ &= \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} \left( \frac{\bm{r}_{1} \cdot \bm{r}}{rr_{1}} \right)^{n-2k} \left( \frac{r_{1}}{r} \right)^{n}\\ &= \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} \left( \bm{r}_{1} \cdot \frac{\bm{r}}{r} \right)^{n-2k} \left( \frac{1}{r_{1}} \right)^{n-2k} \left( \frac{r_{1}}{r} \right)^{n} \\ &= \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} \frac{r^{n-2k+1}}{(n-2k)! (-1)^{n-2k}}\left( \frac{1}{r_{1}} \right)^{n-2k} \left( \frac{r_{1}}{r} \right)^{n} \left( \bm{r}_{1} \cdot \frac{\bm{r}}{r} \right)^{n-2k} \frac{\partial^{n-2k}}{\partial r^{n-2k}} \frac{1}{r} \\ &= r \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!}\frac{1}{(n-2k)! (-1)^{n-2k}} \left( \frac{r}{r_{1}} \right)^{n-2k} \left( \frac{r_{1}}{r} \right)^{n} \left( \bm{r}_{1} \cdot \frac{\bm{r}}{r} \frac{\partial}{\partial r} \right)^{n-2k} \frac{1}{r} \\ &= r \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]} \frac{(-1)^{n-k}(2n-2k)!}{2^{n} k! (n-k)! (n-2k)!(n-2k)!} \left( \frac{r_{1}}{r} \right)^{2k} (\bm{r}_{1} \cdot \nabla)^{n-2k}\frac{1}{r} [\because (23) ] \\ &= r \left( \frac{1}{r} - (\bm{r}_{1} \cdot \nabla) \frac{1}{r} +\frac{1}{2} (\bm{r}_{1} \cdot \nabla)^{2} \frac{1}{r} \cdot\cdot\cdot \right) \tag{29-1} \\ \\ G(Y) &= \sum_{n=0}^{\infty} \frac{1}{n!} \left. \frac{\partial^{n}G}{\partial Y^{n}} \right|_{Y=0} Y^{n} \\ &= \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} \cos^{(n-2k)} \theta_{2} \left( \frac{r_{2}}{r} \right)^{n}\\ &= r \sum_{n=0}^{\infty} \sum_{k=0}^{[n/2]} \frac{(-1)^{n-k}(2n-2k)!}{2^{n} k! (n-k)! (n-2k)!(n-2k)!} \left( \frac{r_{2}}{r} \right)^{2k} (\bm{r}_{2} \cdot \nabla)^{n-2k}\frac{1}{r} \\ &= r \left( \frac{1}{r} - (\bm{r}_{2} \cdot \nabla) \frac{1}{r} +\frac{1}{2} (\bm{r}_{2} \cdot \nabla)^{2} \frac{1}{r} \cdot\cdot\cdot \right) \tag{29-2} となります( $[*]$ とは $*$ を超えない最大の整数のことです)。すると以下の式は近似的に次のように書けます。 (| \bm{r} - \bm{r}_{2} |^{-1} - |\bm{r} - \bm{r}_{1} |^{-1} ) &= \frac{1}{r} \left( 1 - 2 \frac{{r}_{2}}{r} \cos {\theta}_{2} + \left( \frac{r_{2}}{r} \right)^2 \right)^{-\frac{1}{2}} - \frac{1}{r} \left( 1 - 2 \frac{{r}_{1}}{r} \cos {\theta}_{1} + \left( \frac{r_{1}}{r} \right)^2 \right)^{-\frac{1}{2}} \\ &= \frac{G(Y)}{r} - \frac{F(X)}{r} \\ &= (\bm{r}_{2} \cdot \nabla) \frac{1}{r} +\frac{1}{2} (\bm{r}_{2} \cdot \nabla)^{2} \frac{1}{r} \cdot\cdot\cdot \\ &- (\bm{r}_{1} \cdot \nabla) \frac{1}{r} +\frac{1}{2} (\bm{r}_{1} \cdot \nabla)^{2} \frac{1}{r} \cdot\cdot\cdot \\ &\simeq (\bm{r}_{2} - \bm{r}_{1}) \cdot \nabla \frac{1}{r} [\because r \gg l]\\ &= \bm{l} \cdot \nabla \frac{1}{r}\\ &= \frac{\bm{l} \cdot \bm{r}}{r^{3}} \tag{30} ただし、式 $\tag{30}$ の5段目で ${(r \gg l)}$ から $\frac{l}{r}$ 2次以降の項は非常に小さいので切り捨ています。後は $\tag{30}$ を $\tag{27}$ に代入して電位 $V$ が $\tag{24}$ になる事が分かります。ところで $\tag{29}$ で現れた複雑なかたちをした 多項式 P_n (x) = \sum_{k=0}^{[n/2]}\frac{(-1)^{k} (2n - 2k)!}{2^{n} k! (n-k)!(n-2k)!} &x^{n-2k} \\ & (-1 \le x \le 1) \tag{31} ですが、これにはルジャンドル多項式という名前がついています。これについてはまた別の機会に説明したいと思います。 .. _テイラー級数: http://www12.plala.or.jp/ksp/mathInPhys/taylor/index.html @@author: おこめ@@ @@accept: YYYY-MM-DD@@