POINT
- ベクトル解析の公式と,その導出方法の一覧.
- 行列計算も統一的に理解できる.
ベクトル解析の公式と,その導出方法を一覧にまとめました.
力学・電磁気学・流体力学などを学ぶ上で,これらの計算はとても重要です.計算練習をして,すぐに公式を導出できるようになりましょう!
- 前提知識と記法
- 公式と導出
- 積の微分:$\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right) $$= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A} +\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)$
- スカラー三重積:$\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}$
- ベクトル三重積:$\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\cdot \boldsymbol{c})\boldsymbol{b} -(\boldsymbol{a}\cdot \boldsymbol{b})\boldsymbol{c}$
- $\mathrm{rot\,}\mathrm{grad\,}=0$, $\mathrm{div\,}\mathrm{rot\,}=0$
- $\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A}$$=\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}$
- $\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})$$=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}$
- $\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B})$$=(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$+\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}$
- $\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})$$=\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}-\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}-(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}$
- $(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\dfrac{1}{2}\bigl[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) +\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B}) $$-\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}+\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}-\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}\bigr]$
- $(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\boldsymbol{\nabla} (|\boldsymbol{A}|^2/2)$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{A}$
- 積分公式
- ダイアド積
前提知識と記法
微分演算子の略号を使って計算すると,負担を減らすことができます.微分演算子
この記事では簡単のため,微分演算子を\begin{aligned}
\partial_i:=\frac{\partial}{\partial x^i}
\end{aligned}
と書くことにします($x_1=x$, $x_2=y$, $x_3=z$).\partial_i:=\frac{\partial}{\partial x^i}
\end{aligned}
※ 一般的に使われる記法です.分数を使わず,一行で書くことが出来るので便利です.
例えば,
- $\partial_i\phi$は,「$\boldsymbol{\nabla\,}\phi$の$i$成分」という意味
- $\partial_2 f$は,$\dfrac{\partial f}{\partial y}$という意味 ($\partial_y f$とかくこともあります)
です.
ベクトル演算の成分表記
以下では,ダブっている添字については和を取ることにします(Einsteinの規約).すると,ベクトル演算の成分表記が次のように簡単に表せます.名称 | ベクトル表記 | 成分表記 |
---|---|---|
内積 | $\boldsymbol{A}\cdot\boldsymbol{B}$ | $A_i B_i$ |
外積 | $\boldsymbol{A}\times\boldsymbol{B}$ | $\epsilon_{ijk}A_j B_k$ |
勾配 / Gradient | $\mathrm{grad\,}\phi=\boldsymbol{\nabla}\phi$ | $\partial_i\phi$ |
回転 / Rotation / Curl | \begin{aligned}\mathrm{rot\,}\boldsymbol{A}&=\mathrm{curl\,}\boldsymbol{A}\\&=\boldsymbol{\nabla}\times\boldsymbol{A}\end{aligned} |
$\epsilon_{ijk}\partial_jA_k$ |
発散 / Divergence | $\mathrm{div\,}\boldsymbol{A}=\boldsymbol{\nabla}\cdot\boldsymbol{A}$ | $\partial_iA_i$ |
ラプラシアン / Lapracian | \begin{aligned} \Delta \phi&=\mathrm{div\,grad\,}\phi\\&=\boldsymbol{\nabla\,}\cdot\boldsymbol{\nabla\,}\phi \end{aligned} |
$\partial_i\partial_i \phi$ |
行列の成分表記
一般的には,ベクトル演算は行列やテンソルの演算に含まれます.そのような場合には,計算結果が行列式の形にまとめられる場合があります(例:積分の座標変換でヤコビアンを計算する場合).そこで,行列式がどのように表されるのか,知っておく必要があります:$n\times n$行列$A=(a_{i,j})$の行列式は,
\begin{aligned}
\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{1i_1}a_{2i_2}\cdots a_{ni_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_11}a_{i_22}\cdots a_{i_nn}
\end{aligned}
と書くことができます.\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{1i_1}a_{2i_2}\cdots a_{ni_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_11}a_{i_22}\cdots a_{i_nn}
\end{aligned}
この式は,$\epsilon_{i_1i_2\cdots i_n}$の反対称性から
行列式
\begin{aligned}
\epsilon_{j_1j_2\cdots j_n}\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{j_1 i_1}a_{j_2 i_2}\cdots a_{j_n i_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_1 j_1}a_{i_2 j_2}\cdots a_{i_n j_n}
\end{aligned}
\epsilon_{j_1j_2\cdots j_n}\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{j_1 i_1}a_{j_2 i_2}\cdots a_{j_n i_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_1 j_1}a_{i_2 j_2}\cdots a_{i_n j_n}
\end{aligned}
公式の導出手順
ベクトル解析のどんな公式も,以下の手順で導けます!手順
- Step. 1
- 成分表記で書き下す.
- Step. 2
- Step. 3
- 綺麗な形にまとめる.
公式と導出
上の手順に従い,公式を導出します.積の微分:$\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right) $$= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A} +\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)$
積の微分のベクトル解析版です.成分表記すれば,通常の積の微分の公式が使えます.\begin{aligned}
\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right)
&=\partial_i \left( \psi A_i\right) \\
&=\left(\partial_i \psi \right) A_i + \psi \left(\partial_i A_i \right) \\
&= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A}
+\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)
\end{aligned}
\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right)
&=\partial_i \left( \psi A_i\right) \\
&=\left(\partial_i \psi \right) A_i + \psi \left(\partial_i A_i \right) \\
&= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A}
+\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)
\end{aligned}
スカラー三重積:$\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}$
ベクトルの外積が\begin{aligned}
(\boldsymbol{a}\times\boldsymbol{b})_i=\epsilon_{ijk}a_j b_k
\end{aligned}
で表わされることから,(\boldsymbol{a}\times\boldsymbol{b})_i=\epsilon_{ijk}a_j b_k
\end{aligned}
\begin{aligned}
\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})
&=a_i (\epsilon_{ijk} b_j c_k)\\
&=\color{red}{\epsilon_{ijk} a_i b_j c_k}\\
&=(\epsilon_{kij} a_i b_j ) c_k\\
&=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}
\end{aligned}
\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})
&=a_i (\epsilon_{ijk} b_j c_k)\\
&=\color{red}{\epsilon_{ijk} a_i b_j c_k}\\
&=(\epsilon_{kij} a_i b_j ) c_k\\
&=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}
\end{aligned}
ベクトル三重積:$\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\cdot \boldsymbol{c})\boldsymbol{b} -(\boldsymbol{a}\cdot \boldsymbol{b})\boldsymbol{c}$
縮約公式を使うと,\begin{aligned}
[\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})]_i
&=\epsilon_{ijk}a_j(\epsilon_{klm}b_l c_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) a_j b_l c_m \\
&= (\boldsymbol{a}\cdot \boldsymbol{c})b_i -(\boldsymbol{a}\cdot \boldsymbol{b})c_i
\end{aligned}
[\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})]_i
&=\epsilon_{ijk}a_j(\epsilon_{klm}b_l c_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) a_j b_l c_m \\
&= (\boldsymbol{a}\cdot \boldsymbol{c})b_i -(\boldsymbol{a}\cdot \boldsymbol{b})c_i
\end{aligned}
$\mathrm{rot\,}\mathrm{grad\,}=0$, $\mathrm{div\,}\mathrm{rot\,}=0$
微分の順序を交換できるとき*1 以下の重要な性質が成り立ちます:\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k=0
\end{aligned}
\epsilon_{ijk}\partial_j\partial_k=0
\end{aligned}
この式は,反対称性から簡単に導けます:
【証明】
\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k=\epsilon_{ijk}\partial_k\partial_j.
\end{aligned}
ここで$j\leftrightarrow k$とし,$\epsilon_{ijk}$の反対称性を用いると\epsilon_{ijk}\partial_j\partial_k=\epsilon_{ijk}\partial_k\partial_j.
\end{aligned}
\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k
=\epsilon_{ikj}\partial_j\partial_k
=-\epsilon_{ijk}\partial_j\partial_k.
\end{aligned}
であるから,上の式が導かれた.//\epsilon_{ijk}\partial_j\partial_k
=\epsilon_{ikj}\partial_j\partial_k
=-\epsilon_{ijk}\partial_j\partial_k.
\end{aligned}
これを利用すると,以下の公式が導けます:
\begin{aligned}
&(\mathrm{rot\,}\mathrm{grad\,}\phi)_i
=\epsilon_{ijk}\partial_j\partial_k\phi
=0\\
&\mathrm{div\,}\mathrm{rot\,}\boldsymbol{A}
=\partial_i\epsilon_{ijk}\partial_j A_k
=0
\end{aligned}
&(\mathrm{rot\,}\mathrm{grad\,}\phi)_i
=\epsilon_{ijk}\partial_j\partial_k\phi
=0\\
&\mathrm{div\,}\mathrm{rot\,}\boldsymbol{A}
=\partial_i\epsilon_{ijk}\partial_j A_k
=0
\end{aligned}
$\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A}$$=\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}$
\begin{aligned}
(\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A})_i
&=\epsilon_{ijk}\partial_j(\epsilon_{klm}\partial_l A_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl})\partial_j\partial_l A_m \\
&=\partial_i\partial_j A_j-\partial_j\partial_jA_i \\
&=[\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}]_i
\end{aligned}
(\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A})_i
&=\epsilon_{ijk}\partial_j(\epsilon_{klm}\partial_l A_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl})\partial_j\partial_l A_m \\
&=\partial_i\partial_j A_j-\partial_j\partial_jA_i \\
&=[\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}]_i
\end{aligned}
ここで,右辺の$\boldsymbol{\Delta\,A}$は「ベクトルラプラシアン」です.
$\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})$$=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}$
\begin{aligned}
\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})
&=\partial_i(\epsilon_{ijk}A_jB_k) \\
&=\epsilon_{ijk}(\partial_iA_j)B_k+\epsilon_{ijk}A_j(\partial_iB_k) \\
&=B_k\epsilon_{kij}\partial_iA_j - A_j\epsilon_{jik}\partial_iB_k \\
&=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}
\end{aligned}
\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})
&=\partial_i(\epsilon_{ijk}A_jB_k) \\
&=\epsilon_{ijk}(\partial_iA_j)B_k+\epsilon_{ijk}A_j(\partial_iB_k) \\
&=B_k\epsilon_{kij}\partial_iA_j - A_j\epsilon_{jik}\partial_iB_k \\
&=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}
\end{aligned}
$\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B})$$=(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$+\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}$
\begin{aligned}
(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B})_i
&=\epsilon_{ijk} A_j (\mathrm{rot\,}\boldsymbol{B})_k \\
&=\epsilon_{ijk} A_j \cdot \epsilon_{klm}\partial_l B_m \\
&=\epsilon_{kij}\epsilon_{klm} A_j \partial_l B_m \\
&=(\delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}) A_j \partial_l B_m \\
&=A_j \partial_i B_j - A_j \partial_j B_i
\end{aligned}
に着目する.$\boldsymbol{A}$と$\boldsymbol{B}$を入れ替えたものを足し合わせると(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B})_i
&=\epsilon_{ijk} A_j (\mathrm{rot\,}\boldsymbol{B})_k \\
&=\epsilon_{ijk} A_j \cdot \epsilon_{klm}\partial_l B_m \\
&=\epsilon_{kij}\epsilon_{klm} A_j \partial_l B_m \\
&=(\delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}) A_j \partial_l B_m \\
&=A_j \partial_i B_j - A_j \partial_j B_i
\end{aligned}
\begin{aligned}
&(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A})_i \\
&=(A_j \partial_i B_j+B_j \partial_i A_j) - (A_j \partial_j B_i+B_j \partial_j A_i) \\
&=\partial_i(A_j B_j) - (A_j \partial_j B_i) - (B_j \partial_j A_i) \\
&=[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) - (\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B} - (\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A} ]_i
\end{aligned}
&(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A})_i \\
&=(A_j \partial_i B_j+B_j \partial_i A_j) - (A_j \partial_j B_i+B_j \partial_j A_i) \\
&=\partial_i(A_j B_j) - (A_j \partial_j B_i) - (B_j \partial_j A_i) \\
&=[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) - (\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B} - (\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A} ]_i
\end{aligned}
$\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})$$=\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}-\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}-(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}$
\begin{aligned}
&[\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})]_i \\
&=\epsilon_{ijk} \partial_j (\boldsymbol{A}\times \boldsymbol{B})_k \\
&=\epsilon_{ijk} \partial_j (\epsilon_{klm} A_l B_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) [(\partial_j A_l) B_m + A_l (\partial_j B_m)]\\
&=[(\partial_j A_i) B_j + A_i (\partial_j B_j)] - [(\partial_j A_j) B_i + A_j (\partial_j B_i)] \\
&=[(\boldsymbol{B}\cdot \boldsymbol{\nabla})\boldsymbol{A} + \boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}
- \boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A} - (\boldsymbol{A}\cdot \boldsymbol{\nabla})\boldsymbol{B}]_i
\end{aligned}
&[\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})]_i \\
&=\epsilon_{ijk} \partial_j (\boldsymbol{A}\times \boldsymbol{B})_k \\
&=\epsilon_{ijk} \partial_j (\epsilon_{klm} A_l B_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) [(\partial_j A_l) B_m + A_l (\partial_j B_m)]\\
&=[(\partial_j A_i) B_j + A_i (\partial_j B_j)] - [(\partial_j A_j) B_i + A_j (\partial_j B_i)] \\
&=[(\boldsymbol{B}\cdot \boldsymbol{\nabla})\boldsymbol{A} + \boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}
- \boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A} - (\boldsymbol{A}\cdot \boldsymbol{\nabla})\boldsymbol{B}]_i
\end{aligned}
$(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\dfrac{1}{2}\bigl[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) +\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B}) $$-\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}+\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}-\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}\bigr]$
上で求めた2つの関係式を辺々加えることで示される.$(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\boldsymbol{\nabla} (|\boldsymbol{A}|^2/2)$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{A}$
上で$\boldsymbol{A}=\boldsymbol{B}$とすれば導かれる.この関係式は,流体力学でよく使われます.
➡物質微分の意味と関係式(流体力学) - Notes_JP
積分公式
積分公式については以下を参照してください:➡ベクトル解析の公式(積分編) - Notes_JP
ダイアド積
ダイアド積の公式は以下を参照してください:➡ディアディック(ダイアド積)の計算 - Notes_JP
*1:フツーできるときしか考えないので,あまり気にしないように.