線積分

概要

高校で学習する積分では、積分範囲が数直線の上に乗っていた。
つまり、1変数の関数f(x)があって、そのf(x)のグラフの曲線とx軸とに挟まれた領域の面積を求めるというのが高校で学習する積分である。

2変数関数h(x,y)を考える。
この関数は土地の起伏を表しているようなイメージである。地図上の位置を(x,y)で指定すると、h(x,y)がその地点の標高を返すと考えれば理解しやすい。
この関数h(x,y)を積分する。積分範囲は、このxy平面の上を走る自由な曲線コースとする。
自由な曲線コースの上を進みながら積分する、これが線積分である。
もちろん、それを計算するためには、その曲線を具体的に指定しなければならない。

まずは、イメージを説明する。
積分コースに沿ってうねるように立てられた衝立(ついたて)のようなものを考える。
この衝立の高さはその地点でのh(x,y)を意味している。
線積分で計算したいのはこの衝立の面積である。

スカラー場の線積分の計算方法

計算の基礎

xy平面上を走る曲線レールを細かく分割すると、その1つ1つはほぼ直線だとみなせる。
その微小な長さと、その近くでのh(x,y)の値とを掛け合わせたものを考えれば、それは衝立を微小な短冊状に切ったものの面積を表すことになる。
それらを全て加算すれば、望むものが得られる。

この計算をするには、コースの形がtの関数になったx(t)とy(t)で表されていると都合が良い。
コースのスタート地点が(x(a),y(a))であり、ゴール地点が(x(b),y(b))であり、a≦t≦bであるようなtによって(x(t),y(t))で示されるようにする。
(変数tは時刻であるかのようなイメージである)

時刻tからt+Δtまでの微小時間Δt内にコース上の点が動く距離を考える。
最初(x(t),y(t))にあった点が(x(t+Δt),y(t+Δt))にまで移動することになる。微小な時間だから、ほぼ直線的に移動したものと考える。
すると移動距離Δlは三平方の定理により、次式のように表せる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \Delta l=\sqrt{(x(t+ \Delta t)-x(t))^2+(y(t+ \Delta t)-y(t))^2}}

また、ほぼ直線的に変化しているので、x(t+Δt)−x(t)は構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{dx}{dt}\Delta t} 、y(t+Δt)−y(t)は${\frac {dy}{dt}}\Delta t$と近似できる。
$\Delta l={\sqrt {({\frac {dx}{dt}}\Delta t)^{2}+({\frac {dy}{dt}}\Delta t)^{2}}}={\sqrt {({\frac {dx}{dt}})^{2}+({\frac {dy}{dt}})^{2}}}\Delta t$

ここで、Δtが無限に小さいと考えれば、次式のようになる。
$dl={\sqrt {({\frac {dx}{dt}})^{2}+({\frac {dy}{dt}})^{2}}}dt$

このdlが微小時間dtに点がコース上を進む微小距離である。
x(t)とy(t)が具体的に分かっていれば、tのみの関数として表される。(1変数の関数f(x)を積分するときのdxに相当する部分である)
この地点でのh(x,y)の値は、tの時とt+dtの時とでほとんど同じ値なので、h(x(t),y(t))を使用する。

次のようにすれば、線積分の計算ができる。
$\int _{a}^{b}h(x(t),y(t)){\sqrt {({\frac {dx}{dt}})^{2}+({\frac {dy}{dt}})^{2}}}\,dt$

例題. 1

経路Cが以下のC1​、C2のそれぞれの場合、以下の線積分の値を求めよ。
$\int _{C}(2x+y)dx+(x-y)dy$

1. C1​ : 点A(0, 0)から点B(1, 1)へ直線 $y=x$ に沿う経路
2. C2​ : 点A(0, 0)から点B(1, 1)へ放物線 $y=x^{2}$ に沿う経路

$\int _{C}(2x+y)dx+(x-y)dy=\int _{C}(2x+y)dx+\int _{C}(x-y)dy$である。

1.では、$y=x$のため、第2項は$(x-y)=(x-x)=0$となる。
${\begin{aligned}\int _{C}(2x+y)dx+(x-y)dy&=\int _{C}(2x+y)dx\\&=\int _{0}^{1}3xdx\\&=3{\Big [}{\frac {x^{2}}{2}}{\Big ]}_{0}^{1}\\&={\frac {3}{2}}\end{aligned}}$

2.では、${\frac {dy}{dx}}=2x\quad {\mbox{よ り }}\quad dy=2xdx$である。
${\begin{aligned}\int _{C}(2x+y)dx+(x-y)dy&=\int _{0}^{1}(2x+x^{2})dx+\int _{0}^{1}(x-x^{2})2xdx\\&={\Big [}x^{2}+{\frac {x^{3}}{3}}{\Big ]}_{0}^{1}+2{\Big [}{\frac {x^{3}}{3}}-{\frac {x^{4}}{4}}{\Big ]}_{0}^{1}\\&={\Big (}1+{\frac {1}{3}}{\Big )}+2{\Big (}{\frac {1}{3}}-{\frac {1}{4}}{\Big )}\\&={\frac {4}{3}}+{\frac {1}{6}}\\&={\frac {3}{2}}\end{aligned}}$

上記ではスカラー関数として考えてきたが、見方を変えてみる。
変位ベクトル${\overrightarrow {r}}=\langle dx,dy\rangle$とベクトル場${\overrightarrow {F}}=\langle 2x+y,x-y\rangle$を考える時、上記の線積分はベクトルの内積を用いて以下のように記述できる。
$\int _{C}(2x+y)dx+(x-y)dy=\int _{C}{\overrightarrow {F}}\cdot d{\overrightarrow {r}}$
$\nabla \times {\overrightarrow {F}}={\overrightarrow {0}}$(渦なし)の時、線積分は経路によらず始点と終点で決まる。

上記の例では、${\frac {\partial }{\partial x}}(x-y)-{\frac {\partial }{\partial y}}(2x+y)=1-1=0$となり、渦なしの条件を満たしているため、
1.および2.の2通りの線積分において、同様の計算結果となる。

線積分の制約

上記で説明した線積分の公式の中には、x(t)とy(t)を微分したものが出てくる。
ということは、x(t)とy(t)とはどちらも滑らかでないと微分ができない。
しかし、幾つかの地点で微分ができない場合は、そこでコースを分割して別々に計算して、後で足し合わせればよい。

3次元の線積分

上記の線積分を3次元に拡張することは簡単である。

3変数関数U(x,y,z)を考え、3次元空間を自由に飛び回るコースで積分する。
要するに、3次元を走る微小な長さと、その近くでのU(x,y,z)の値とを掛けせたものを端から端までの和をとる。

考え方は上記と同様で、次のようにする。
$\int _{a}^{b}U(x(t),y(t),z(t)){\sqrt {({\frac {dx}{dt}})^{2}+({\frac {dy}{dt}})^{2}+{\frac {dz}{dt}}^{2}}}\,dt$

ベクトル場の線積分の計算

計算の基礎

${\overrightarrow {F}}=F_{1}(x,y,z){\overrightarrow {i}}+F_{2}(x,y,z){\overrightarrow {j}}+F_{3}(x,y,z){\overrightarrow {k}}$のベクトル空間内の2点AからBを結ぶ曲線経路Cがある。

この曲線上で、ベクトル${\overrightarrow {F}}$の接線方向成分の大きさを表すスカラー関数$F_{t}(s)$を考える。(sはAからの孤の長さである)
※接線はタンジェントラインと呼ぶため、接線方向成分の意味でtを添字にしている。

$F_{t}(s)$の線積分は、次式となる。
$\int _{C}F_{t}(s)ds$

$F_{t}(s)$は、曲線経路Cの単位接線ベクトル${\overrightarrow {t}}(s)$とベクトル${\overrightarrow {F}}$の内積として求まるため、以下のように記述できる。
$\int _{C}F_{t}(s)ds=\int _{C}{\overrightarrow {F}}\cdot {\overrightarrow {t}}ds$

この量をベクトル${\overrightarrow {F}}$の線積分(接線線積分)という。
曲線経路Cをsのパラメータ(媒介変数)として、${\overrightarrow {r}}(s)=x(s){\overrightarrow {i}}+y(s){\overrightarrow {j}}+z(s){\overrightarrow {k}}$と表すと、
単位接線ベクトル${\overrightarrow {t}}$は、${\overrightarrow {t}}={\frac {d{\overrightarrow {r}}}{ds}}$となる。

したがって、ベクトル${\overrightarrow {F}}$の線積分は、以下のようにも記述できる。
${\begin{aligned}\int _{C}{\overrightarrow {F}}\cdot {\overrightarrow {t}}ds&=\int _{C}{\overrightarrow {F}}\cdot {\frac {d{\overrightarrow {r}}}{ds}}ds\\&=\int _{C}{\overrightarrow {F}}\cdot d{\overrightarrow {r}}\end{aligned}}$

記述方法が複数あるが、重要なことは、曲線経路Cにおいてベクトル場の接線成分を積分するということである。

例題.1

ベクトル関数${\overrightarrow {F}}(x,y)=-3x^{2}{\overrightarrow {i}}+5xy{\overrightarrow {j}}$について、経路Cにおける以下の線積分を求めよ。

ただし、経路Cは放物線$y=2x^{2}$の(0, 0)から(1, 2)に沿う曲線とする。
$\int _{C}{\overrightarrow {F}}\cdot d{\overrightarrow {r}}$

まず、経路Cをパラメータtで表示する。つまり、xやyをtの関数として記述する。
線積分を行う経路となる経路Cは放物線$y=2x^{2}$の(0, 0)から(1, 2)である。
$t=x$とすると、$0\leq t\leq 1$で$y=2t^{2}$である。
「経路Cについて考える」ということは、「xとyが経路Cの方程式を満たす」ということである。

したがって、xとyはtの関数として以下のように記述できる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} x(t) &= t \quad (0 \le t \le 1) \\ y(t) &= 2t^2 \end{align} }

よって、xとyをtでそれぞれ微分すると以下のようになる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \frac{dx}{dt} &= 1 \quad \mbox{よ り } \quad dx = dt \\ \frac{dy}{dt} &= 4t \quad \mbox{よ り } \quad dy = 4t dt \end{align} }

式を変形するときのポイントは、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \overrightarrow{r}} の意味から、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \overrightarrow{r} = dx \overrightarrow{i} + dy \overrightarrow{j} + dz \overrightarrow{k}} と記述する点である。
この例題では、ベクトル関数はz成分は無いため、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \overrightarrow{r} = dx \overrightarrow{i} + dy \overrightarrow{j}} としている。

上記のポイントを踏まえて、線積分を計算する。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_C \overrightarrow{F} \cdot d \overrightarrow{r} &= \int_C (-3x^2 \overrightarrow{i} + 5xy \overrightarrow{j}) \cdot (dx \overrightarrow{i} + dy \overrightarrow{j}) \\ &= \int_C (-3x^2 dx + 5xy dy) \end{align} }
ドット(⋅)は内積を表している。
ここで、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \overrightarrow{F}(x, y) = -3x^2 \overrightarrow{i} + 5xy \overrightarrow{j}} を成分表示すると、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{pmatrix} -3x^2 \\ 5xy \end{pmatrix}} である。
同様に、 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle dx \overrightarrow{i} + dy \overrightarrow{j}} の成分表示は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{pmatrix} dx \\ dy \end{pmatrix}} である。

次に、xとyをパラメータt(ここでは、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x = t, \ y = 4t dt} を代入する)で記述する。これが、経路Cを満たすように式変形するということである。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_C (-3x^2 dx + 5xy dy) &= \int_{0}^{1} -3t^2 dt + 5t(2t^2) 4t dt \\ &= \int_{0}^{1} -3t^2 dt + 40 t^4 dt \\ &= \int_{0}^{1} (-3t^2 + 40t^4) dt \\ &= \Big [ -t^3 + 8t^5 \Big]_{0}^{1} \\ &= -1 + 8 \\ &= 7 \end{align} }

例題. 2

ベクトル場構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \overrightarrow{F} = a(\sin t \overrightarrow{i} + \cos t \overrightarrow{j})}
(aは正の定数)において、
経路Cを構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C: \overrightarrow{r} = \cos t \overrightarrow{i} + \sin t \overrightarrow{j} + bt \overrightarrow{k} \quad \Big(0 \le t \le \frac{\pi}{4} \Big)}
)とする時、
線積分構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_C \overrightarrow{F} \cdot d \overrightarrow{r}}
を求めよ。

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \overrightarrow{r} = \cos t \overrightarrow{i} + \sin t \overrightarrow{j} + bt \overrightarrow{k}} をtで微分すると、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d \overrightarrow{r}}{dt} = - \sin t \overrightarrow{i} + \cos t \overrightarrow{j} + b \overrightarrow{k}} となり、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \overrightarrow{r} = (-\sin t \overrightarrow{i} + \cos t \overrightarrow{j} + b \overrightarrow{k}) dt} である。

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d \overrightarrow{r}} を成分表示すると、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d\overrightarrow{r} = \begin{pmatrix} -\sin t \\ \cos t \\ b \end{pmatrix} dt} となる。
また、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \overrightarrow{F}} を成分表示すると、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \overrightarrow{F} = \begin{pmatrix} a \cos t \\ a \cos t \\ 0 \end{pmatrix}} である。

内積はx成分、y成分、z成分それぞれ乗算した後に加算して求めるため、以下のように計算できる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \overrightarrow{F} \cdot d \overrightarrow{r} &= \begin{pmatrix} a \sin t \\ a \cos t \\ 0 \end{pmatrix} \cdot \begin{pmatrix} - \sin t \\ \cos t \\ b \end{pmatrix} dt \\ &= a (- \sin^2 t + \cos^2 t) dt \\ &= a \cos 2t \qquad \because \cos{2 \theta} = \cos^2{\theta} - \sin^2{\theta} \end{align} }

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_C \overrightarrow{F} \cdot d \overrightarrow{r} &= \int_{0}^{\frac{\pi}{4}} a \cos 2t dt \\ &= \Big[ \frac{a \sin2t}{2}\Big]_0^{\frac{\pi}{4}} \\ &= \frac{a}{2} \Big\{ \sin(2 \times \frac{\pi}{4}) - \sin 0 \Big\} \\ &= \frac{a}{2} (1 - 0) \\ &= \frac{a}{2} \end{align} }

ベクトル量の線積分

物理学では、ベクトルの概念と組み合わさった形の別の線積分が出てくるので、これだけで線積分が完全に理解できない。
それについては、別のページで記載する。