==========================================
テンソル楕円体
==========================================
この記事は、 テンソルの主軸１_ 、 テンソルの主軸２_ の続きです。しばらく工学上の応用のため二階のテンソルばかりを扱い、行列の固有値や対角化など、行列の話題が出てきましたので、もう一度注意を喚起しておきます。


.. important:: 

	テンソルを行列のことだと思ってはいけません


テンソルとは、あくまでも テンソルの概念_ に紹介した定義を満たす量で、代数的には、多重線形性を満たすようにベクトル空間からテンソル積によって構成されるものでした。二階のテンソルを行列表記することが便利なため、どうしても応用的なテンソルの計算には行列がよく出て来ることになり、その結果『テンソルとは（ $3 \times 3$ の）行列のことだ』と思い込んでいる人さえいるようなのです。恐ろしいことです。


注意を喚起したところで、再び二階のテンソルの話に戻ります。これも工学上重要な応用的話題です^^


テンソル楕円体
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
対称テンソル $T_{ij}$ を使って表わされる次の二次曲面を考えます。

<tex>
T_{ij}x^{i}x^{j}=1		\tag{1}
</tex>


二階の対称テンソルですから、 $T_{ij}=T_{ji}$ に注意すると、式 $(1)$ は一般に次のように展開できることが分かるでしょう。 二次曲線の係数_ を参照して下さい

<tex>
(1) \ \rightarrow \  &  T_{11}x^{1}x^{1} + T_{22}x^{2}x^{2} + T_{33}x^{3}x^{3} \\  & \ \ \  2T_{12}x^{1}x^{2} + 2T_{23}x^{2}x^{3} + 2T_{31}x^{3}x^{1}=1 \tag{2}
</tex>


テンソル $T_{ij}$ が固有値 $\lambda_{1},\lambda_{2},\lambda_{3}$ を持ち、主方向に座標系 $(x_{1}^{*},x_{2}^{*},x_{3}^{*})$ を取るとき、式 $(2)$ を次式のように簡単化することが出来るでしょう。

<tex>
T_{11}(x_{1}^{*})^2 + T_{22}(x_{2}^{*})^2 + T_{33}(x_{3}^{*})^2 = \lambda_{1} (x_{1}^{*})^2 + \lambda_{2} (x_{2}^{*})^2 + \lambda_{3} (x_{3}^{*})^2 =1  \tag{3}
</tex>


特に、 $T_{ij}$ を対称テンソルだと仮定しているため、主軸は互いに直交しますので、 $(x_{1}^{*},x_{2}^{*},x_{3}^{*})$ 座標系として正規直交座標系を取ることが出来ます。これは大変に便利な話です。式 $(3)$ を次のように書き換えると、これが軸半径 $\frac{1}{\sqrt{\lambda_{1}}},\frac{1}{\sqrt{\lambda_{2}}},\frac{1}{\sqrt{\lambda_{3}}}$ の楕円体を表わす式だということが分かると思います。


<tex>
\left( \frac{x_{1}^{*}}{1/\sqrt{\lambda_{1}}} \right)^2 + \left( \frac{x_{2}^{*}}{1/\sqrt{\lambda_{2}}} \right)^2 + \left( \frac{x_{3}^{*}}{1/\sqrt{\lambda_{3}}} \right)^2 = 1	\tag{4}
</tex>


.. image:: Joh-TensorEllipsoid1.gif 




.. [*] テンソルの話題というより、単に行列の固有値と対角化の問題であるような気もします。冒頭の注意をもう一度思い出してください。応用上、応力テンソルや慣性モーメントを楕円体の形で表わすと便利な場合があります。例えば、自由に運動する剛体は、慣性楕円体が中心の高さを一定に保ちつつ、平面に滑らずに接して動くときの運動として幾何学的に捉えることが可能です。（楕円体の形は慣性モーメントによって、平面からの高さは剛体の持つ運動量によって決まります。楕円体を使って剛体の運動を考える方法を $polehode$ 法と言います。）また、その他の楕円体の応用として、結晶の屈折率に異方性を示す屈折率楕円体というものもあります。これは、様々な方向から入射した直線偏光の屈折率をテンソル楕円体で表したもので、光が入射する方向からの楕円体の断面を見ることによって、幾何的に屈折率を知ることができます。


補足：対称テンソルではない場合
---------------------------------------------------------------------------
最初にテンソル $T_{ij}$ は対称テンソルだと仮定しましたが、これを一般のテンソルだとするとどうなるのでしょうか。この場合、 $T_{ij}$ は対称部分 $S_{ij}$ と反対称部分 $A_{ij}$ に分離可能です。（ 対称テンソルと反対称テンソル_ 参照。）

<tex>
T_{ij} = S_{ij} +A_{ij}		\tag{5}
</tex>


式 $(1)$ に式 $(5)$ を代入してみます。


<tex>
T_{ij}x^{i}x^{j}=S_{ij}x^{i}x^{j}+A_{ij}x^{i}x^{j}=1
</tex>



ここで、 $A_{ij}=-A_{ji}$ かつ $x_{i}x_{j}=x_{j}x_{i}$ であることを考えると、 $A_{ij}x^{i}x^{j}=0$ となってしまい、対称部分しか残らないことになります。すなわち、一般の二階のテンソルでも、互いに直交する主軸を持つテンソル楕円体を描くことは出来ますが、その楕円体の方程式はテンソルの対称部分のみによって規定されることになり、対称部分が共通なテンソルのテンソル楕円体は全て同一になります。



.. _対称テンソルと反対称テンソル: http://www12.plala.or.jp/ksp/vectoranalysis/SymtensorAntisymtensor/
.. _テンソルの概念: http://www12.plala.or.jp/ksp/vectoranalysis/TensorConcept/
.. _テンソルの主軸１: http://www12.plala.or.jp/ksp/vectoranalysis/TensorPrincipalAxis1/
.. _テンソルの主軸２: http://www12.plala.or.jp/ksp/vectoranalysis/TensorPrincipalAxis2/
.. _二次曲線の係数: http://www12.plala.or.jp/ksp/vectoranalysis/QuadricCoeff/

@@author:Joh・黒子@@
@@accept: 2006-08-25@@
@@category: ベクトル解析@@
@@id: TensorEllipsoid@@