ベクトル解析

この項目では、数学理論としてのベクトル解析について説明しています。1901年にエドウィン・ビッドウェル・ウィルソンとウィラード・ギブスによって出版されたベクトル解析に関する著作『Vector Analysis』については「ベクトル解析 (著書)」をご覧ください。

ベクトル解析（ベクトルかいせき、英語：vector calculus）は空間上のベクトル場やテンソル場に関する微積分に関する数学の分野である。

多くの物理現象はベクトル場やテンソル場として記述されるため、ベクトル解析は物理学の様々な分野に応用を持つ。

物理学では3次元ユークリッド空間上のベクトル解析が特によく用いられるが、ベクトル解析は一般のn次元多様体上で展開できる。

3次元ユークリッド空間におけるベクトル解析

ベクトル場とスカラー場

定義

3次元ユークリッド空間${\displaystyle \mathbb {R} ^{3}}$上のベクトル場Xとは、${\displaystyle \mathbb {R} ^{3}}$上の各点Pに対し、Pを始点とする3次元ベクトルX(P)を対応させる写像のことである。

ベクトル場の例（2次元の場合）

本項では特に断りのない限り、この写像がPに関して滑らかな場合を考える。すなわち、${\displaystyle \mathbb {R} ^{3}}$の座標を使って

${\displaystyle \mathbf {X} (P)=(X_{1}(x_{1},x_{2},x_{3}),X_{2}(x_{1},x_{2},x_{3}),X_{3}(x_{1},x_{2},x_{3}))}$

と表したとき、各X_i(x₁,x₂,x₃)が任意回微分可能である場合を考える。

なお、ベクトル場の記法としてX(P)の代わりに

${\displaystyle \mathbf {X} _{P}}$

のように Pを下付きに書くことも多い。しかしこの下付きの記法だと、成分表記したときに煩雑になるので、本項ではX(P)とX_Pの両方の記法を混用する。

同様に3次元ユークリッド空間${\displaystyle \mathbb {R} ^{3}}$上のスカラー場Fとは、${\displaystyle \mathbb {R} ^{3}}$上の各点Pに対し、実数F(P)を対応させる写像のことである。

具体例

ベクトル場の例としては、電場、磁場、重力場などがある。また流体上の各点にその点での粒子の速度ベクトルを対応させることで速度場を定義する事もできる。

スカラー場の線積分

定義

Fをスカラー場とし、 ユークリッド空間上の曲線Cが${\displaystyle \mathbf {x} (u)=(x_{1}(u),x_{2}(u),x_{3}(u))}$、u∈[a,b]とパラメーター表示されているとする。積分

${\displaystyle I=\int _{a}^{b}F(\mathbf {x} (u))\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} u}(u)\right\|\operatorname {d} u}$

を曲線Cに沿ったスカラー場Fの線積分という。

線積分の定義はCとその向き付けには依存するが、同じ向き付けを与える限りパラメーターの取り方に依存しない。実際曲線Cを別のパラメーターwに ${\displaystyle \mathbf {x} (w)=(x_{1}(u(w)),x_{2}(u(w),x_{3}(u(w)))}$、w∈[c,d]と変数変換して同様の積分

${\displaystyle J=\int _{c}^{d}F(\mathbf {x} (w(u)))\left\|{\operatorname {d} \mathbf {x} (u(w)) \over \operatorname {d} w}\right\|\operatorname {d} w}$

を考えると、この変数変換が曲線の向きを変えないとき、すなわち

${\displaystyle {\operatorname {d} u \over \operatorname {d} w}>0}$

が恒等的に言えるときには

${\displaystyle {\begin{array}{ll}J&=\displaystyle \int _{c}^{d}F(\mathbf {x} (u(w)))\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} w}\right|\operatorname {d} w\\&=\displaystyle \int _{a}^{b}F(\mathbf {x} (u(w)))\left|{\operatorname {d} \mathbf {x} \over \operatorname {d} u}\right\|{\operatorname {d} u \over \operatorname {d} w}\operatorname {d} w\\&=\displaystyle \int _{a}^{b}F(\mathbf {x} (u))\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} u}\right\|\operatorname {d} u\\&=I\end{array}}}$

が成立する。

そこで線素dsを

${\displaystyle \operatorname {d} s:=\|\operatorname {d} \mathbf {x} \|={\sqrt {\operatorname {d} x_{1}{}^{2}+\operatorname {d} x_{2}{}^{2}+\operatorname {d} x_{3}{}^{2}}}}$

と定義し、スカラー場の線積分を

${\displaystyle \int _{C}F\operatorname {d} s}$

と表記する。

Cの始点と終点が一致するとき（すなわちCが閉曲線のとき）はそのことを強調して

${\displaystyle \oint _{C}F\operatorname {d} s}$

とも表記する。

弧長

線積分の特殊なケースとして

${\displaystyle \int _{C}\operatorname {d} s}$

を考えると、曲線Cの長さ（弧長）に一致する事が知られている。

厳密な証明は弧長の項目にゆずるが、直観的には以下の理由による。Cを${\displaystyle \mathbf {x} (u)=(x_{1}(u),x_{2}(u),x_{3}(u))}$、u∈[a,b] と向きをはじめとする保つようにパラメトライズし、[a,b]を長さΔuの微小区間に分けると、Cの長さはおよそ

${\displaystyle \sum _{i}\|\mathbf {x} ((i+1)\Delta u)-\mathbf {x} (i\Delta u)\|}$ ${\displaystyle \sum _{i}\|\Delta \mathbf {x} /\Delta u\|\Delta u}$

なので、Δuを0に近づけると、線積分

${\displaystyle \int \left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} u}\right\|\operatorname {d} u}$

に一致する。従って上述の線積分で弧長を求める事ができる。

弧長パラメーター

曲線Cを、

x(s)=(Cの始点からs離れた位置)

とパラメトライズできる^{[注 1]}。このようなsをCの弧長パラメーターという。 Cを弧長パラメーターで${\displaystyle \mathbf {x} (s)=(x_{1}(s),x_{2}(s),x_{3}(s))}$と表したとき、定義より

${\displaystyle s=\int _{0}^{s}\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} s}\right\|\operatorname {d} s}$

なので、両辺を微分すると、

${\displaystyle 1=\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} s}\right\|}$

が恒等的に成り立つ。

従って線積分とは、

${\displaystyle \int _{C}F(x(s))\operatorname {d} s=\int _{0}^{1}F(x(s))\left\|{\operatorname {d} \mathbf {x} \over \operatorname {d} s}\right\|\operatorname {d} s=\int _{0}^{1}F(x(s))\operatorname {d} s}$

より、弧長でパラメトライズされた場合のF(x(s))の積分である。

ベクトル場の線積分

定義

Xをベクトル場とし、 ユークリッド空間上の曲線Cが${\displaystyle \mathbf {x} (u)=(x_{1}(u),x_{2}(u),x_{3}(u))}$、u∈[a,b]とパラメーター表示されているとし、積分

${\displaystyle \int _{a}^{b}\mathbf {X} (\mathbf {x} (u))\cdot {\operatorname {d} \mathbf {x} \over \operatorname {d} u}(u)\operatorname {d} u}$

を考える。ここで「・」は内積である。

スカラー場に対する線積分と同様の議論で、上述の積分はベクトル場Xと曲線Cのみに依存し、Cのパラメトライズの方法によらない。そこで上述の積分を

${\displaystyle \int _{C}\mathbf {X} \cdot \operatorname {d} \mathbf {x} }$

と表記し、ベクトル場Xの曲線Cに沿った線積分という。ここで

${\displaystyle \operatorname {d} \mathbf {x} =(\operatorname {d} x_{1},\operatorname {d} x_{2},\operatorname {d} x_{3})}$

である。成分で書けば、線積分は

${\displaystyle \int _{C}X_{1}\operatorname {d} x_{1}+X_{2}\operatorname {d} x_{2}+X_{3}\operatorname {d} x_{3}}$

とも表示できる。

Cの始点と終点が一致するとき（すなわちCが閉曲線のとき）はそのことを強調して

${\displaystyle \oint _{C}\mathbf {X} \cdot \operatorname {d} \mathbf {x} }$

とも表記する。

弧長パラメーターによる表示

Cを弧長パラメーターで${\displaystyle \mathbf {x} (s)=(x_{1}(s),x_{2}(s),x_{3}(s))}$と表すと、 Cに沿った線積分は、

${\displaystyle \int \mathbf {X} (\mathbf {x} (s))\cdot {\operatorname {d} \mathbf {x} \over \operatorname {d} s}(s)\operatorname {d} s}$

と表記できる。　すでに示したように

${\displaystyle \left|{\operatorname {d} \mathbf {x} \over \operatorname {d} s}\right|=1}$

が恒等的に成り立つので、内積

${\displaystyle \mathbf {X} (\mathbf {x} (s))\cdot {\operatorname {d} \mathbf {x} \over \operatorname {d} s}(s)}$

は ${\displaystyle \mathbf {X} (\mathbf {x} (s))}$をCのx(s)での接線方向の射映である。

すなわち線積分は、ベクトル場Xの、 Cの接線方向成分を積分したものである。

スカラー場の面積分

定義

3次元ユークリッド空間内の曲面Sが

${\displaystyle \mathbf {x} =\mathbf {x} (u_{1},u_{2})}$

とパラメトライズされていたとする。このとき、スカラー場FのS上での面積分を

${\displaystyle \int F(\mathbf {x} (u_{1},u_{2}))\left\|{\partial \mathbf {x} \over \partial u_{1}}\times {\partial \mathbf {x} \over \partial u_{2}}\right\|\operatorname {d} u_{1}\operatorname {d} u_{2}}$

により定義する。

Sのパラメーターを

${\displaystyle (u_{1},u_{2})=(u_{1}(w_{1},w_{2}),u_{2}(w_{1},w_{2}))}$

と変数変換しても、この変数変換がSの向き付けを変えないなら、すなわちヤコビアン

${\displaystyle {\partial (u_{1},u_{2}) \over \partial (w_{1},w_{2})}>0}$

が恒等的に成り立つなら、面積分の値は替わらないことを容易に示せる。

そこでFのS上での面積分を

${\displaystyle \int _{S}F\operatorname {d} S}$

とSのパラメトライズの方法によらない形で表記する。

1の面積分

${\displaystyle \int _{S}\operatorname {d} S}$

はSの面積に等しい事が知られており、従って ${\displaystyle \operatorname {d} S}$は面積の微小量を表していると考えられる。この${\displaystyle \operatorname {d} S}$の事を面素という。

Sが閉曲面のときはそのことを強調して、面積分の事を

${\displaystyle \oint _{S}\operatorname {d} S}$

とも表記する。

ベクトル場の面積分

向き付けられた曲面S上の点PにおけるSの流さ1の法線（単位法線）をn_Pとする。なお、 PにおけるSの単位法線は2本あるが、そのうちSの向きとn_Pが右手系になるものをn_Pとする。

このとき、ベクトル場XのS上での面積分を

${\displaystyle \int _{S}\mathbf {X} \operatorname {d} \mathbf {S}$ :=\int _{S}\mathbf {X} \cdot \mathbf {n} \operatorname {d} S}

により定義する。

Sが

${\displaystyle \mathbf {x} =\mathbf {x} (u_{1},u_{2})}$

とパラメトライズされている場合、面積分の定義から、

${\displaystyle \int _{S}\mathbf {X} \operatorname {d} \mathbf {S} =\int \mathbf {X} \cdot \left({\partial \mathbf {x} \over \partial u_{1}}\times {\partial \mathbf {x} \over \partial u_{2}}\right)\operatorname {d} u_{1}\operatorname {d} u_{2}}$

である。積分内はベクトル3重積であるので、

${\displaystyle \mathbf {X} \cdot \left({\partial \mathbf {x} \over \partial u_{1}}\times {\partial \mathbf {x} \over \partial u_{2}}\right)=\operatorname {det} \left(\mathbf {X} ,{\partial \mathbf {x} \over \partial u_{1}},{\partial \mathbf {x} \over \partial u_{2}}\right)}$

でもある。

勾配、回転、発散

定義

Fをスカラー場とするとき、Fの勾配 gradFをベクトル場

${\displaystyle \operatorname {grad} F={\partial F \over \partial \mathbf {x} }=\left({\partial F \over \partial x_{1}},{\partial F \over \partial x_{2}},{\partial F \over \partial x_{3}}\right)}$

により定義する。

さらにベクトル場Xの回転 rot X、発散 div Xをそれぞれベクトル場

${\displaystyle \operatorname {rot} \mathbf {X} =\left({\frac {\partial X_{3}}{\partial x_{2}}}-{\frac {\partial X_{2}}{\partial x_{3}}},{\frac {\partial X_{1}}{\partial x_{3}}}-{\frac {\partial X_{3}}{\partial x_{1}}},{\frac {\partial X_{2}}{\partial x_{1}}}-{\frac {\partial X_{1}}{\partial x_{2}}}\right)}$

${\displaystyle \operatorname {div} \mathbf {X} ={\partial X_{1} \over \partial x_{1}}+{\partial X_{2} \over \partial x_{2}}+{\partial X_{3} \over \partial x_{3}}=\operatorname {tr} \left({\partial X_{i} \over \partial x_{j}}\right)_{ij}}$

により定義する。

微分演算子ナブラ${\displaystyle \nabla }$を

${\displaystyle \nabla =\left({\partial \over \partial x_{1}},{\partial \over \partial x_{2}},{\partial \over \partial x_{3}}\right)}$

と定義すると、勾配、回転、発散は

${\displaystyle \operatorname {grad} F=\nabla F}$

${\displaystyle \operatorname {rot} \mathbf {X} =\nabla \times \mathbf {X} }$

${\displaystyle \operatorname {div} \mathbf {X} =\nabla \cdot \mathbf {X} }$

と表記できる。

ストークスの定理とガウスの定理

勾配、回転、発散と線積分、面積分は以下の関係を満たす。ここでF、Xはそれぞれ3次元ユークリッド空間上のスカラー場とベクトル場、C、S、Vは3次元ユークリッド空間内の有界な曲線、曲面、および3次元領域で、「∂」は境界を表し、P、QはそれぞれCの始点と終点を表す。

${\displaystyle \int _{C}\operatorname {grad} F\operatorname {d} s=F(Q)-F(P)}$

${\displaystyle \int _{S}\operatorname {rot} \mathbf {X} \operatorname {d} \mathbf {S} =\oint _{\partial S}\mathbf {X} \operatorname {d} s}$　（（ケルビン・）ストークスの定理）

${\displaystyle \int _{V}\operatorname {div} \mathbf {X} \operatorname {d} V=\oint _{\partial V}\mathbf {X} \operatorname {d} \mathbf {S} }$　（ガウスの発散定理）

1パラメーター変換

発散divの幾何学的意味を見るため、ベクトル場の1パラメーター変換という概念を導入する。

Xを3次元ユークリッド空間${\displaystyle \mathbb {R} ^{3}}$上のベクトル場とし、xを${\displaystyle \mathbb {R} ^{3}}$の点とする。

Φ_u(x)を以下のように定義する：

${\displaystyle \Phi _{u}(\mathbf {x} )}$をベクトル場Xに沿ってuだけ進んだ点、すなわち

${\displaystyle \left.{\partial \over \partial u}\Phi _{u}(\mathbf {x} )\right|_{u=u'}=\mathbf {X} _{\Phi _{u'}(\mathbf {x} )}}$

が全てのu' ∈[0,u]に対して成り立つ点とする。

このようなΦ_u(x)は全ての(x,u)に対して定義できるとは限らないが^{[注 2]}、xの近傍Uとε>0を十分小さく選べば、任意のx' ∈Uと任意のu' ∈[0,ε]に対してこのようなΦ_u(x' )を定義できることが知られている。このような写像Φ_uをベクトル場Xの1パラメーター変換という。

発散divの幾何学的意味

1パラメーター変換をもちいると、発散divを幾何学的に意味づける事ができる。 Φ_u(x' )を成分で y₁(x' ), y₂(x' ), y₃(x' ), と書くことにすると、体積要素はヤコビアンを用いて

${\displaystyle \operatorname {d} y_{1}\operatorname {d} y_{2}\operatorname {d} y_{3}=\left|{\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }\right|\operatorname {d} x_{1}\operatorname {d} x_{2}\operatorname {d} x_{3}}$

という関係式を満たす。すなわち、 Φ_uは点xにおいて微小体積を体積比

${\displaystyle \left|{\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }\right|}$

で変換する写像である。

ヤコビの公式より、

${\displaystyle {\partial \over \partial u}\left|{\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }\right|=\operatorname {tr} \left({\widetilde {\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }}{\partial \over \partial u}{\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }\right)}$

ここで、${\displaystyle {\tilde {A}}}$はAの余因子行列である。 Φ₀は恒等写像なので、Iを単位行列とすると、1パラメーター変換の定義より、

${\displaystyle {\partial \over \partial u}\left|{\operatorname {d} \Phi _{u} \over \operatorname {d} \mathbf {x} }\right|_{u=0}=\operatorname {tr} \left({\tilde {I}}{\operatorname {d} \mathbf {X} \over \operatorname {d} \mathbf {x} }\right)=\operatorname {div} \mathbf {X} }$

すなわち、div Xは 微小体積の1パラメーター変換による変化率を表している。

ポアンカレの補題とポテンシャル

簡単な計算により、任意のスカラー場Fと任意のベクトル場Xに対し

${\displaystyle \operatorname {rot} (\operatorname {grad} F)=0}$

${\displaystyle \operatorname {div} (\operatorname {rot} \mathbf {X} )=0}$

が恒等的に成立する事が簡単な計算により確認できる。

また3次元ユークリッド空間上では次が成立する（ポアンカレの補題）：

${\displaystyle \operatorname {rot} \mathbf {X} }$が恒等的に0 ${\displaystyle \iff \operatorname {grad} \phi =\mathbf {X} }$となるφが存在する

${\displaystyle \operatorname {div} \mathbf {X} }$が恒等的に0 ${\displaystyle \iff \operatorname {rot} \mathbf {A} =\mathbf {X} }$となるAが存在する

このようなφ、Aが存在するとき、φ、AをそれぞれXのスカラー・ポテンシャル、ベクトル・ポテンシャルという。

なお、ポアンカレの補題が成り立つのはユークリッド空間では1次以上のコホモロジー（ド・ラームコホモロジー）が消えている事と関係しており、一般の多様体では必ずしもこの補題は成り立たない。

スカラー・ポテンシャル、ベクトル・ポテンシャルとも、存在する場合には一意ではない。しかし、φ₁、φ₂を同一のベクトル場Xのスカラー・ポテンシャルとするとき、

${\displaystyle \phi _{2}(\mathbf {x} )=\phi _{1}(\mathbf {x} )+{\text{const.}}}$

である事が容易に示せる。

またA₁、A₂を同一のベクトル場Xのスカラー・ポテンシャルとするとき、

${\displaystyle \mathbf {A} _{2}=\mathbf {A} _{1}+\operatorname {grad} F}$

を満たす Fが必ず存在する。 実際、ベクトル・ポテンシャルの定義より、

${\displaystyle \operatorname {rot} (\mathbf {A} _{2}-\mathbf {A} _{1})=\operatorname {rot} \mathbf {A} _{2}-\operatorname {rot} \mathbf {A} _{1}=\mathbf {X} -\mathbf {X} =0}$

なので、ポアンカレの補題より

${\displaystyle \mathbf {A} _{2}-\mathbf {A} _{1}=\operatorname {grad} F}$

となるFが存在する。

歴史

現代の学校教育では古典力学の導入からベクトルを用いた物理教育が行われ、数学でも幾何ベクトル・線型代数学・ベクトル解析といったベクトルの概念が普通に教えられている。しかし古典力学の登場と同時にベクトルも誕生したのではなく、物理法則などを表記するために19世紀に生まれ^[1]、20世紀になり高次元ベクトル場にまで一般化された。

ベクトルが誕生するまでは直交座標系を用いた解析幾何学やウィリアム・ローワン・ハミルトンが考案した四元数を用いた記法が主流であり、力学・電磁気学の教育・研究でも解析幾何学的な多変数微積分学を用いた力学や四元数表記の電磁気学が普通であった^[1]。余談だが、同じようにベクトルを扱う数学理論である線型代数も登場時期はほぼ同じであり、こちらは完成が遅れたため教育に本格的に導入されるのは20世紀後半、数学教育の現代化が言われ出した頃である。20世紀前半は教えられている物理数学が現代とは違っていたのであり、ベクトルは数学ではなく物理学の授業で導入され、行列式が先に教えられていたし^[2]、行列を用いて量子力学を定式化したヴェルナー・ハイゼンベルクも線型代数を習っていなかった。日本でも明治初期の物理教育では、四元数に基づく電磁気学が教えられていたことは有名である。

ベクトルを初めて教育に導入したのはウィラード・ギブスとされ、1880年代のイェール大学の講義で記号こそ現代とは違うものの、外積・内積やベクトル解析の概念などが当時使われていたが、イギリスの四元数の著書もある物理学者ピーター・ガスリー・テイトの評判も大変不評であったという^[1]。今日用いられている記号や専門用語の大半は1901年に出版されたギブスとエドウィン・ウィルソン（英語版）の共著『ベクトル解析』によって確立された。

しかし、ギブス以降の物理学の教育ではベクトルは四元数を推進していたハミルトンやテイトのいたイギリスにおいて寧ろ盛んに用いられるようになり、物理学における常識的な概念となった^[1]。（イギリスのオリヴァー・ヘヴィサイドの存在が影響していると考えられる。）しかしながら20世紀に入ってからはむしろスピン角運動量などの概念も四元数に非常に類似しており、ハミルトンには先見性があったのではないかとされる^[1]。

関連概念

• 場の微分
  • 曲率
  • ナブラ
    • 勾配・発散・回転
    • 偏微分

• 線積分・面積分（微分形式の積分）
  • グリーンの定理
  • 発散定理（ジョージ・グリーン、カール・フリードリヒ・ガウス）

脚注

[脚注の使い方]

注釈

1. ここでは話を簡単にするためパラメーターsの原点をCの始点にしたが、C上の任意のん一点を原点としてよい。この場合、原点より始点側は負の値でパラメトライズする。
2. 例えば、 u→1のときΦ_u(x)→∞となれば、1以上のuに対して Φ_u(x)を定義できない。

出典

1. 湯川秀樹　『物理講義』、1975年、講談社、58-62頁
2. 銀林浩、『線型代数学序説』、現代数学社、2002年 まえがきより

参考文献

• H・K・ニッカーソン、D・C・スペンサー、N・E・スティーンロッド著、原田重春・佐藤正次訳『現代ベクトル解析』、岩波書店、1965年、ISBN 4-00-005293-4 C3041
• 岩堀長慶、『ベクトル解析　―力学の理解のために―』、裳華房、〈数学選書2〉、1960年。 ISBN 978-4-7853-1302-9

関連項目

• ベクトル解析の公式の一覧
• 物理学関係
  • マクスウェルの方程式
  • 解析力学
• 数学関係
  • ド・ラームコホモロジー
  • テンソル
  • ベルソル（英語版）
• 人物
  • ジェームズ・クラーク・マクスウェル
  • ヘルマン・フォン・ヘルムホルツ
  • オリヴァー・ヘヴィサイド
  • ピーター・テイト
  • チャールズ・プロテウス・スタインメッツ