計算の基礎
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は三平方の定理により、次式のように表せる。
また、ほぼ直線的に変化しているので、 と近似できる。
ここで、Δtが無限に小さいと考えれば、次式のようになる。
dlが微小時間dtに点がコース上を進む微小距離である。(dtとは、1変数の関数f(x)を積分するときのdxに相当する部分である)
x(t)とy(t)が具体的に分かっていれば、tのみの関数として表される。
この地点での の値は、tの時と の時とでほとんど同じ値なので、 を使用する。
次のようにすれば、線積分の計算ができる。
例題. 1
経路Cが以下のC1、C2のそれぞれの場合、以下の線積分の値を求めよ。
1. C1 : 点A(0, 0)から点B(1, 1)へ直線 に沿う経路
2. C2 : 点A(0, 0)から点B(1, 1)へ放物線 に沿う経路
である。
1.では、 となるが、 のため、第2項は となる。
2.では、 である。
上記ではスカラー関数として考えてきたが、見方を変えてみる。
変位ベクトル とベクトル場 を考える時、上記の線積分はベクトルの内積を用いて以下のように記述できる。
(渦なし)の時、線積分は経路によらず始点と終点で決まる。
上記の例では、
渦なしの条件を満たしているため、
1.および2.の2通りの線積分において、同様の計算結果となる。
計算の基礎
のベクトル空間内の2点AからBを結ぶ曲線経路Cがある。
この曲線上で、ベクトル の接線方向成分の大きさを表すスカラー関数 を考える。(sはAからの孤の長さである)
※接線はタンジェントラインと呼ぶため、接線方向成分の意味でtを添字にしている。
の線積分は、次式となる。
は、曲線経路Cの単位接線ベクトル とベクトル の内積として求まるため、以下のように記述できる。
この量をベクトル の線積分(接線線積分)という。
曲線経路Cをsのパラメータ(媒介変数)として、 と表すと、
単位接線ベクトル は、 となる。
ただし、 は、曲線経路Cを細かく分割した微小ベクトルであり、線素ベクトルと呼ばれる。
したがって、ベクトル の線積分は、以下のようにも記述できる。
記述方法が複数あるが、重要なことは、曲線経路Cにおいてベクトル場の接線成分を積分するということである。
例題.1
ベクトル関数 について、経路Cにおける以下の線積分を求めよ。
ただし、経路Cは放物線 の(0, 0)から(1, 2)に沿う曲線とする。
まず、経路Cをパラメータtで表示する。つまり、xやyをtの関数として記述する。
線積分を行う経路となる経路Cは放物線 の(0, 0)から(1, 2)である。
とすると、 で である。
「経路Cについて考える」ということは、「xとyが経路Cの方程式を満たす」ということである。
したがって、xとyはtの関数として以下のように記述できる。
よって、xとyをtでそれぞれ微分すると以下のようになる。
式を変形するときのポイントは、 の意味から、 と記述する点である。
この例題では、ベクトル関数はz成分は無いため、 としている。
上記のポイントを踏まえて、線積分を計算する。
ここで、ドット(⋅)は内積を表している。
ここで、 を成分表示すると、 である。
同様に、 の成分表示は、 である。
次に、xとyをパラメータt(ここでは、 を代入する)で記述する。
これが、経路Cを満たすように式変形するということである。
例題. 2
ベクトル場 (aは正の定数)において、
経路Cを )とする時、
線積分 を求めよ。
をtで微分すると、
となり、
である。
を成分表示すると、 となる。
また、 を成分表示すると、 である。
内積はx成分、y成分、z成分それぞれ乗算した後に加算して求めるため、以下のように計算できる。
例題. 3
aを定数として、 というベクトル場を考える。
この時、次の問を求めよ。
(1) 下図左のように、点A : (1, 0, 0)、点B : (1, 1, 0)、点C : (0, 1, 0)の3点をつないだ閉曲線cを考える。
この閉曲線Cに沿った の線積分を求めよ。
ただし、線積分の向きは、下図左の矢印の方向(反時計方向)を正にとる。
(2) 下図右のように、点A : (1, 0, 0)から点C : (0, 1, 0)まで、原点Oを中心とする半径1の円の円周に沿って曲線ℓを引く。
この曲線lに沿った の線積分を求めよ。
ただし、線積分の向きは、下図右の矢印の方向(反時計方向)を正にとる。
(1)の求め方
経路C上の位置ベクトル を1つの変数で表す。
経路Cは3本の線分から構成されているため、各線分の線積分を計算して加算する。
(i) 線分AB上での線積分
線分AB上の位置ベクトルを とする時、 とおくと、 となる。
より、線素ベクトル は、 となる。
さらに、ベクトル場 に位置ベクトル を適用( を代入)すると、 となる。
よって、経路ABの矢印の向きを考慮すると、積分範囲は となるため、線分AB上での線積分は次式となる。
(ii) 線分BC上での線積分
線分BC上の位置ベクトルを とする時、 とおくと、 となる。
より、線素ベクトル は、 となる。
さらに、ベクトル場 に位置ベクトル を適用( を代入)すると、 となる。
よって、経路BCの矢印の向きを考慮すると、積分範囲は となるため、線分BC上での線積分は次式となる。
(iii) 線分CA上での線積分
線分CA上の位置ベクトルを とする時、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x = t, y = -t + 1}
とおくと、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{r_3} = \langle t, -t + 1, 0 \rangle}
となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{\vec{dr_3}}{dt} = \langle 1, -1, 0 \rangle = \begin{pmatrix} 1 \\ -1 \\ 0 \end{pmatrix}}
より、線素ベクトル構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{dr_3}}
は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{dr_3} = \langle 1, -1, 0 \rangle dt = \begin{pmatrix} 1 \\ -1 \\ 0 \end{pmatrix} dt}
となる。
さらに、ベクトル場構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{F}}
に位置ベクトル構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{r_3}}
を適用(構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x = t, y = -t + 1}
を代入)すると、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{F} = \langle at -a, at, 0 \rangle = \begin{pmatrix} at - a \\ at \\ 0 \end{pmatrix}}
となる。
よって、経路CAの矢印の向きを考慮すると、積分範囲は構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t:0 \rightarrow 1}
となるため、線分CA上での線積分は次式となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{CA} \vec{F} \cdot \vec{dr_3} &= \int_{0}^{1} \langle at - a, at, 0 \rangle \cdot \langle 1, -1, 0 \rangle \\ &= \int_{0}^{1} (at - a - at) dt \\ &= \int_{0}^{1} -a dt \\ &= -a \int_{0}^{1} dt \\ &= -a \Big[ t \Big]_{0}^{1} \\ &= -a \end{align} }
したがって、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int_{C} \vec{F} \cdot \vec{dr} = \int_{AB} \vec{F} \cdot \vec{dr_1} + \int_{BC} \vec{F} \cdot \vec{dr_2} + \int_{CA} \vec{F} \cdot \vec{dr_3} = a + a - a = a}
となる。
(2)の求め方
経路lは原点Oを中心とする半径1の円である。
経路l上の位置ベクトルを構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{r}}
とする時、極座標表示を用いて構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x = cos \theta, y = \sin \theta}
とおくと、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{r} = \langle \cos \theta, \sin \theta, 0 \rangle}
となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{\vec{dr}}{dt} = \langle - \sin \theta, \cos \theta, 0 \rangle = \begin{pmatrix} - \sin \theta \\ \cos \theta \\ 0 \end{pmatrix}}
より、線素ベクトル構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{dr}}
は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{dr} = \langle - \sin \theta, \cos \theta, 0 \rangle dt = \begin{pmatrix} - \sin \theta \\ \cos \theta \\ 0 \end{pmatrix} dt}
となる。
さらに、ベクトル場構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{F}}
に位置ベクトル構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{r}}
を適用(構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x = \cos \theta, y = \sin \theta}
を代入)すると、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \vec{F} = \langle -a \sin \theta, a \cos \theta, 0 \rangle = \begin{pmatrix} -a \sin \theta \\ a \cos \theta \\ 0 \end{pmatrix}}
となる。
よって、経路lの矢印の向きを考慮すると、積分範囲は構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t:0 \rightarrow \frac{\pi}{2}}
となるため、経路l上での線積分は次式となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \int_{l} \vec{F} \cdot \vec{dr} &= \int_{0}^{\frac{\pi}{2}} \langle -a \sin \theta, a \cos \theta, 0 \rangle \cdot \langle - \sin \theta, \cos \theta, 0 \rangle \\ &= \int_{0}^{\frac{\pi}{2}} (a \sin^2 \theta + a \cos^2 \theta) dt \\ &= a \int_{0}^{\frac{\pi}{2}} (\sin^2 \theta + \cos^2 \theta) dt \\ &= a \int_{0}^{\frac{\pi}{2}} dt \qquad \because \sin^2 \theta + \cos^2 \theta = 1 \\ &= a \Big[ t \Big]_{0}^{\frac{\pi}{2}} \\ &= \frac{a \pi}{2} \end{align} }