物理のかぎしっぽ 記事ソース/面積素と微分形式

記事ソース/面積素と微分形式

これはrst2hooktailの記事ソース保存・変換用です(詳細).

コンバート

最近コンバートされた結果: HTMLPDFTeX

公開・更新メニュー ▼▲

記事ソースの内容

==========================================================
面積素と微分形式
==========================================================
曲面 $S$ が $\bm{r}=\bm{r}(u,v)$ で与えられているとき、その面積素は次式で与えられるのでした。面積素というのは、面積分するときに出て来る $dS$ ってやつです。ただし、 $S$ の法線ベクトルは $\bm{e_{u}} \times \bm{e_{v}}$ と同じ向きを正とします。( 面積素に関して補足_ も参照して下さい。)


<tex>
d\bm{S} = \left( \frac{\partial \bm{r}}{\partial u} \times \frac{\partial \bm{r}}{\partial v} \right) dudv 	\tag{1}
</tex>


.. image:: Joh-DiffSquare01.gif 


滑らかな曲面上の領域が、上図のように二つの独立なパラメーターで表わせるという事は直観的に了承して、先に進むことにしましょう。



.. [*] いま考えている点の近傍で曲面は滑らかだとしています。曲面上の領域が二つのパラメーターで表わせるという事をもう少し正確に言えば、この曲面上のどの点の近傍も、二次元ユークリッド空間と同じだと見ることが出来るということです。(このような曲面を、二次元多様体と呼びます。)ここでは、あまり深入りしませんが、詳しくは、 多様体_ カテゴリーを参考にして下さい。


.. image:: Joh-DiffSquare02.gif 



曲面 $S$ 上の一点 $P_{0}$ のまわりに微小面積を取ることを考えます。点 $P_{0}$ の位置ベクトルを $\bm{r}$ とし、上図のように $\bm{r}$ が $u$ とともに変化する方向( $\frac{\partial \bm{r}}{\partial u}$ の方向)に $P_{1}$ 、 $\bm{r}$ が $v$ とともに変化する方向( $\frac{\partial \bm{r}}{\partial v}$ の方向)に $P_{2}$ をとります。曲面は本当は曲がっているわけですが、いま微小量を考えているので、 $P_{0}P_{1}$ や  $P_{0}P_{2}$ はほぼ直線とみなし、 $dS$ として、 $\vec{P_{0}P_{1}}$ と $\vec{P_{0}P_{2}}$ が張る平行四辺形を考えます。


<tex>
d\bm{S} &= \vec{P_{0}P_{1}} \times \vec{P_{0}P_{2}} \\
&= (\vec{P_{1}}- \vec{P_{0}}) \times (\vec{P_{2}} - \vec{P_{0}}) \\ 
& = \left[ \left( \bm{r} + \frac{\partial \bm{r}}{\partial u}du \right) - \bm{r} \right] \times \left[ \left( \bm{r} + \frac{\partial \bm{r}}{\partial v}dv \right) - \bm{r} \right] \\ 
& = \frac{\partial \bm{r}}{\partial u}du \times \frac{\partial \bm{r}}{\partial v}dv \\ 
& =\left( \frac{\partial \bm{r}}{\partial u} \times \frac{\partial \bm{r}}{\partial v} \right) dudv \tag{2}
</tex>


偏導関数をいちいち分数の形に書くのは面倒なので、以後、微分変数を添字にして、例えば $\frac{\partial \bm{r}}{\partial u}$ を $\bm{r_{u}}$ のように略記することにします。 



.. [*] 細かいことを言えば、曲面上の微小面積は曲率を持っていて、平行四辺形による表現はあくまで近似のはずです。 $\Delta \bm{S} \sim (\bm{r_{u}} \times \bm{r_{v}}) dudv$ 。そして、 $\Delta \rightarrow 0$ の極限において初めて、 $\sim$ が $=$ になります。(つまり、曲率を持った面積素と平行四辺形の差が $0$ に収束します。)このような二段階で、 $\varepsilon -\delta$ 論法を使って示すのが正確ですが、式 $(2)$ は少しいい加減に最初から『面積素 $=$ 平行四辺形』として計算してしまいました。厳密に考えたい人は、詳しい教科書をあたって下さい。



蛇足ですが、曲面の第一基本量と呼ばれる量、 $E=\bm{r_{u}} \cdot \bm{r_{u}},  F=\bm{r_{u}} \cdot \bm{r_{F}},  G=\bm{r_{v}} \cdot \bm{r_{v}}$ を使えば、式 $(2)$ を次のように簡単に書くこともできます。



<tex>
dS = \sqrt{EG-F^{2}}dudv	\tag{3}
</tex>



ヤコビアンと面積素
=====================================================================
曲面 $S$ を表わすベクトル方程式 $\bm{r}(u,v)$ を用いて、式 $(2)$ もしくは $(3)$ の形で面積素 $dS$ を表現することが出来ました。曲面上のパラメーター二つで表わすのは簡潔で良いですが、応用上の便宜も考えて、次のようにデカルト座標系 $xyz$ で $\bm{r}(u,v)$ を表現することも考えてみましょう。

<tex>
\bm{r}(u,v) = x(u,v)\bm{e_{x}}+ y(u,v)\bm{e_{y}}+z(u,v)\bm{e_{z}}	\tag{4}
</tex>


これを使うと、式 $(2)$ は次のようになります。


<tex>
d\bm{S} & = ( \bm{r_{u}} \times \bm{r_{v}})dudv \\ 
&= \left( 
    \begin{array}{c}
\frac{\partial x}{\partial u} \\
\frac{\partial y}{\partial u} \\
\frac{\partial z}{\partial u} \\
    \end{array}
  \right) 
\times 
\left( 
    \begin{array}{c}
\frac{\partial x}{\partial v} \\
\frac{\partial y}{\partial v} \\
\frac{\partial z}{\partial v} \\
    \end{array}
  \right) dudv \\
&= 
  \left|
    \begin{array}{cc}
\frac{\partial y}{\partial u}  & \frac{\partial y}{\partial v} \\
\frac{\partial z}{\partial u}  & \frac{\partial z}{\partial v} \\
    \end{array}
  \right|
dudv  \bm{e_{x}} +
  \left|
    \begin{array}{cc}
\frac{\partial z}{\partial u}  & \frac{\partial z}{\partial v} \\
\frac{\partial x}{\partial u}  & \frac{\partial x}{\partial v} \\
    \end{array}
  \right|
dudv  \bm{e_{y}} + 
  \left|
    \begin{array}{cc}
\frac{\partial x}{\partial u}  & \frac{\partial x}{\partial v} \\
\frac{\partial y}{\partial u}  & \frac{\partial y}{\partial v} \\
    \end{array}
  \right|
dudv \bm{e_{z}}	\tag{5}
</tex>


なんだか、行列式がたくさん出てきてエライことになってますが、これらは $x,y,z$ のうちの二変数から $u,v$ へ座標変換する際の *ヤコビアン* になっていることが分かるかと思います。(ヤコビアンを忘れてしまっている人は、やばいですよ。すぐに復習して下さい!)




<tex>
\frac{\partial (x,y)}{\partial (u,v)} = 
  \left|
    \begin{array}{cc}
\frac{\partial x}{\partial u}  & \frac{\partial x}{\partial v} \\
\frac{\partial y}{\partial u}  & \frac{\partial y}{\partial v} \\
    \end{array}
  \right|
=\frac{\partial x}{\partial u}\frac{\partial y}{\partial v} - 
\frac{\partial x}{\partial v}\frac{\partial y}{\partial u}	\tag{6-1}
</tex>



<tex>
\frac{\partial (y,z)}{\partial (u,v)} = 
  \left|
    \begin{array}{cc}
\frac{\partial y}{\partial u}  & \frac{\partial y}{\partial v} \\
\frac{\partial z}{\partial u}  & \frac{\partial z}{\partial v} \\
    \end{array}
  \right|
=
\frac{\partial y}{\partial u} \frac{\partial z}{\partial v} - 
\frac{\partial y}{\partial v}\frac{\partial z}{\partial u} 	\tag{6-2}
</tex>



<tex>
\frac{\partial (z,x)}{\partial (u,v)} = 
  \left|
    \begin{array}{cc}
\frac{\partial z}{\partial u}  & \frac{\partial z}{\partial v} \\
\frac{\partial x}{\partial u}  & \frac{\partial x}{\partial v} \\
    \end{array}
  \right|
=
\frac{\partial z}{\partial u} \frac{\partial x}{\partial v} -
\frac{\partial z}{\partial v}\frac{\partial x}{\partial u}	\tag{6-3}
</tex>


この表記を使って $d\bm{S}$ をもう一度書き直せば、次のようになります。


<tex>
d\bm{S} = \frac{\partial (y,z)}{\partial (u,v)} dudv \bm{e_{x}} + 
\frac{\partial (z,x)}{\partial (u,v)} dudv \bm{e_{y}} + 
\frac{\partial (x,y)}{\partial (u,v)} dudv \bm{e_{z}} 	\tag{7}
</tex>


ずいぶん、すっきりしました。


ヤコビアンと微分形式
=====================================================================
前セクションでは、 $d\bm{S}$ を $xyz$ 成分に分解してみました。さて、式 $(7)$ から、第一成分の $\frac{\partial (y,z)}{\partial (u,v)} dudv$ だけを例に取って、さらに座標変換について考えてみましょう。このヤコビアンの符号には正負があり、 $dydz$ と $dudv$ を次式のように関係づけます。


<tex>
dydz= \pm \frac{\partial (y,z)}{\partial (u,v)} dudv	\tag{8}
</tex>

特に、正の符号を取る変換を *向きを保つ座標変換* と呼びます。座標系の向きを直観的に理解しやすい例は三次元直交座標で、右手系を右手系に、もしくは左手系を左手系に保つ座標変換が、正のヤコビアンによって表わされることになります。(後で、この話に戻って来ますから、ヤコビアンの符号と、座標変換の向きという話を覚えておいて下さい。)


それにしても、符号が $\pm$ となっているのは、いかにも面倒くさいですねぇ。



ヤコビアンの性質
`````````````````````````````````````````````````````````````````````
ここで、一般にヤコビアンの持つ性質を考えて見たいと思います。 行列式の性質_ により、ヤコビアンは次の性質を持つことが分かると思います。

<tex>
\frac{\partial (x,y)}{\partial (u,v)} =- \frac{\partial (y,x)}{\partial (u,v)} 
, \ \ \ 
\frac{\partial (y,z)}{\partial (u,v)} =- \frac{\partial (z,y)}{\partial (u,v)} 
, \ \ \ 
\frac{\partial (z,x)}{\partial (u,v)} =- \frac{\partial (x,z)}{\partial (u,v)} 
</tex>


.. important::  

	ヤコビアンの符号は、行か列を入れ替えると反転します。



また、次の性質もあることが言えます。


<tex>
\frac{\partial (x,x)}{\partial (u,v)} = 0 
, \ \ \ 
\frac{\partial (y,y)}{\partial (u,v)} = 0
, \ \ \ 
\frac{\partial (z,z)}{\partial (u,v)} = 0
</tex>



.. important:: 

	同じ行(もしくは同じ列)があると、ヤコビアンは $0$ になります。



座標の向きとの関係
```````````````````````````````````````````````````````````````
このセクションの最初に少し考えた、座標変換で向きを保つ、という話に戻ります。二つの座標軸を入れ替える座標変換のヤコビアンは負になりますので、向きを保ちません。このことは、例えば次のような座標変換のヤコビアンが $-1$ になることを見れば、明らかでしょう。

<tex>
\left( 
    \begin{array}{c}
\bm{e_{2}} \\
\bm{e_{1}} \\ 
    \end{array}
  \right) 
= 
\left( 
    \begin{array}{cc}
0 & 1 \\
1 & 0 \\
    \end{array}
  \right) 
\left( 
    \begin{array}{c}
\bm{e_{1}} \\
\bm{e_{2}} \\ 
    \end{array}
  \right) 	\tag{9-1}
</tex>

<tex>
\left( 
    \begin{array}{c}
\bm{e_{1}} \\
\bm{e_{3}} \\ 
\bm{e_{2}} \\
    \end{array}
  \right) 
= 
\left( 
    \begin{array}{ccc}
1 & 0  & 0 \\
0 & 0  & 1 \\
0 & 1  & 0 \\ 
    \end{array}
  \right) 
\left( 
    \begin{array}{c}
\bm{e_{1}} \\
\bm{e_{2}} \\ 
\bm{e_{3}} \\
    \end{array}
  \right) 	\tag{9-2}
</tex>

このように、『座標軸を入れ替える』という操作は、実はヤコビ行列の行(もしくは列)を入れ替える操作に等しくなります。式 $(9-1)(9-2)$ で考えたヤコビ行列の行列式が、単位行列の行を一回入れ替えることで得られることを確認して下さい。座標軸の入れ替えも、二回行えばヤコビアンの符号は元に戻ります。そこで、『座標変換が向きを保つか』という問題は、考えているヤコビ行列が、向きを保つ座標変換(例えば単位行列で表わされる恒等変換)の表現行列の行(もしくは列)を、何回入れ替えたものなのか、より具体的には、 *偶数回入れ替えたか奇数回入れ替えたか* という問題に帰着します。


なぜ、『向きを保つ』という問題にこだわっているのか、少し面倒臭く感じている人がいるかも知れません。私達が $dudv$ とか $dydz$ という面積素を考えているのは、それを使って面積分をしたい、という目論見が背後にあるからです。そして、面積分においては *曲面の向き* がとても重要でした。(例えば ガウスの発散定理_ を復習してみて下さい。)私達は、曲面上の領域を表わすのに、自由自在に座標変換したい訳ですが、このとき、面積素の面積だけでなく、曲面の向き(つまり曲面の法線ベクトルの方向)も自然に保つような座標変換の体系を考えたいのです。これが、向きにこだわっている動機です。(だって、積分のたびに向きをごちゃごちゃ考えるのは面倒ですよね。)


そこで、いっそ式 $(8)$ を次のように書いてしまえば便利でしょう。

<tex>
dy\land dz=  \frac{\partial (y,z)}{\partial (u,v)} dudv	\tag{10}
</tex>


いきなりウェッジ積が出てきましたが、『積の順序を(奇数回)変えれば符号が変わる』というウェッジ積の性質と、『二つの行や列を(奇数回)入れ替えれば符号が変わる』というヤコビアンの性質が、見事にマッチしている訳です。これで、最初に基準となる $uv$ 座標系の向きさえ決めておけば(つまり、曲面の表・裏を最初に決めるということです。どちらが表でどちらが裏か、という決定には任意性があります)、その後は座標変換に伴う向きの変化を気にしなくて済むことになりますo(^ー^)o。



ヤコビアンの二番目の性質、『同じ行や列があると零になる』という性質ですが、これもウェッジ積の『同じもののウェッジ積は零になる』という性質 $dx \land dx =0$ に、見事にマッチしています。どうやら、二次微分形式を使うと、 *向きに関するヤコビアンの性質を織り込んだ面積素* を考えることができそうです。



.. important:: 

	二次微分形式は、曲面の向き(表・裏)に関する情報も織り込んだ面積素だと考えても良いです。


これは、面積分で重宝しそうですね!微分形式を使えば、次のような表現が可能です。


<tex>
dx \land dy = \frac{\partial (x,y)}{\partial (u,v)}dudv
</tex>

<tex>
dy \land dz = \frac{\partial (y,z)}{\partial (u,v)}dudv
</tex>

<tex>
dz \land dx = \frac{\partial (z,x)}{\partial (u,v)}dudv
</tex>

これを使って、式 $(7)$ を次のように簡単な形にできます。もう一度、式 $(7)$ と比べてみて下さい。

<tex>
d\bm{S} = dy \land dz + dz \land dx + dx \land dy \tag{11}
</tex>



例えば、ガウスの発散定理に出てくる面積分は、式 $(11)$ を使って表現すると次のようになります。

<tex>
\int \limits_{S} \bm{A} \cdot d\bm{S}
= 
\int \limits _{S}
A_{x}dy \land dz + A_{y}dz \land dx  + A_{z}dx \land dy	\tag{12}
</tex>


ここで式 $(12)$ の右辺は二次微分形式になっていることが分かると思います。図形的に、式 $(12)$ は『単位時間あたりに面 $S$ を通過する $\bm{A}$ の *流量* 』という意味を持っていましたが、これはそのまま微分二形式の物理的意味として重要な解釈です。




.. [*] 代数的に、ベクトル空間→テンソル空間→外積代数と勉強して来た人は、外積空間なんていうものは抽象的思考の産物で、実際の役に立たないのではないかと感じていたかも知れません。(実際、グラスマン( $\text{Hermann}$  ${\rm G\ddot{u}nter}$  $\text{Grassmann (1809-1877)}$ )が外積代数を考え出したときは、これが何の役に立つのか、誰も分かりませんでした。)この記事では曲面の向きと面積分という観点だけから微分形式を見たわけですが、外積代数の基底を $dx,dy,dz$ にしてみただけで、『曲面の向き』『座標変換』『面積分の向き』といった(面倒な)話題が、綺麗に、微分形式の表現に収まってしまったのを見ました。外積代数も、やはり役に立たない抽象思考の産物では無かったわけです。なんと言っても、微分形式は美しいです。



.. _面積素に関して補足: http://www12.plala.or.jp/ksp/vectoranalysis/SurfaceIntegralApp/
.. _多様体: http://www12.plala.or.jp/ksp/manifold/index.html
.. _行列式の性質: http://www12.plala.or.jp/ksp/linearalgebra/DeterminantsProperties/
.. _行列の階数: http://www12.plala.or.jp/ksp/linearalgebra/DeterminantsProperties/MtrxRank/
.. _ガウスの発散定理: http://www12.plala.or.jp/ksp/vectoranalysis/GaussDivTheorem/


@@author:Joh@@
@@accept:2006-11-13@@
@@category: 微分形式@@
@@id: DiffFormsArea@@
トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Modified by 物理のかぎプロジェクト PukiWiki 1.4.6 Copyright © 2001-2005 PukiWiki Developers Team. License is GPL.
Based on "PukiWiki" 1.3 by yu-ji Powered by PHP 5.3.29 HTML convert time to 0.011 sec.