四面体の体積と行列式(volume of a tetrahedron and determinant)とは、四面体の体積を頂点の座標から求める式である。四面体 $\mathrm{OABC}$ の体積は、$\vec a=\overrightarrow{\mathrm{OA}}$、$\vec b=\overrightarrow{\mathrm{OB}}$、$\vec c=\overrightarrow{\mathrm{OC}}$ の 3 次の行列式(スカラー三重積)の絶対値の $\frac16$ に等しい。4 点が同じ平面の上にあれば両方とも $0$ である。絶対値を外した値は向きで符号の変わる符号つき体積で、各ベクトルについて 1 次なので、辺の上の点で切った四面体の体積比が倍率の積になる。分母の $6$ は、立方体が合同な $3!$ 個の四面体に分かれることから来る。
前提知識: 空間ベクトルと1次独立, ベクトルの内積(高校数学), 空間図形の計量, 錐体の体積とHilbertの第3問題
四面体の体積は、錐体の体積とHilbertの第3問題 で示したとおり $\dfrac13\times$ 底面積 $\times$ 高さ である。頂点の座標が分かっていても、高さが見えていないと、この公式を使うには手間がかかる。
$\mathrm O(0,0,0)$、$\mathrm A(1,0,0)$、$\mathrm B(0,2,0)$、$\mathrm C(0,0,3)$ とする。$\mathrm{OA}$、$\mathrm{OB}$、$\mathrm{OC}$ は座標軸の上にあって互いに垂直である。底面を直角三角形 $\mathrm{OAB}$ とみると、底面積は $\dfrac12\cdot1\cdot2=1$、高さは $\mathrm{OC}=3$ なので、体積は
$$
V=\frac13\cdot1\cdot3=1
$$
である。
$\mathrm O(0,0,0)$、$\mathrm A(2,1,0)$、$\mathrm B(0,2,1)$、$\mathrm C(1,1,3)$ とし、$\vec a=\overrightarrow{\mathrm{OA}}=(2,1,0)$、$\vec b=\overrightarrow{\mathrm{OB}}=(0,2,1)$、$\vec c=\overrightarrow{\mathrm{OC}}=(1,1,3)$ とおく。高校の方法で、底面 $\mathrm{OAB}$ の面積と、$\mathrm C$ から平面 $\mathrm{OAB}$ までの距離を求める。
ex-vtd-slanted では、平方根が途中に出てきて最後に消え、答えは $\dfrac{11}6$ という有理数になった。頂点の座標が整数なら、体積はいつも $\dfrac16$ の整数倍になるのだろうか。答えは「なる」で、その整数は 3 つのベクトルの成分から掛け算と引き算だけで求まる。この記事で答える問いは次の 4 つである。
| 高校の計算 | この記事の言葉 | 大学の言葉 |
|---|---|---|
| $\dfrac13\times$ 底面積 $\times$ 高さ | 行列式の絶対値の $\dfrac16$(thm-vtd-volume) | 単体の体積 $\dfrac{1}{n!}\lvert\det\rvert$ |
| 平面の三角形の面積 $\dfrac12\lvert x_1y_2-x_2y_1\rvert$ | 3 次の行列式(def-vtd-det) | 行列式=符号つき体積 |
| 4 点が同じ平面の上にあるか | 行列式が $0$ か | 1 次従属 |
| 頂点を底面に平行に動かしても体積は同じ | 行列式は行の足し算で変わらない(lem-vtd-det-rules) | 交代多重線形形式 |
| 辺の上の点で切った四面体の体積比 | 行列式の各行の 1 次性(prop-vtd-ratio) | 多重線形性 |
四面体 OABC と、ベクトル a、b、c の張る平行六面体(点線)。四面体の体積は平行六面体の体積の 6 分の 1 になる
空間のベクトル $\vec a=(a_1,a_2,a_3)$、$\vec b=(b_1,b_2,b_3)$、$\vec c=(c_1,c_2,c_3)$ に対して、
$$
\det(\vec a,\vec b,\vec c)=a_1b_2c_3+a_2b_3c_1+a_3b_1c_2-a_3b_2c_1-a_2b_1c_3-a_1b_3c_2
$$
を $\vec a$、$\vec b$、$\vec c$ の 行列式 という。$\vec a$、$\vec b$、$\vec c$ を上から順に行として並べた表
$$
\begin{pmatrix}a_1&a_2&a_3\\ b_1&b_2&b_3\\ c_1&c_2&c_3\end{pmatrix}
$$
の 行列式 と同じものである。
6 つの項は、どれも「$\vec a$ の成分 1 つ、$\vec b$ の成分 1 つ、$\vec c$ の成分 1 つ」の積で、3 つの成分の番号($1,2,3$)はすべて異なる。番号の並び方は $3!=6$ 通りあり、それが 6 つの項に対応する。表で右下がりの 3 本の斜めの線の積に $+$、右上がりの 3 本の積に $-$ を付けて足す、と覚える方法を Sarrus の方法 という。
この行列式は、空間ベクトルと1次独立 で 1 次独立の判定に使ったものと同じ式であり、外積 $\vec a\times\vec b$ を使うと $\det(\vec a,\vec b,\vec c)=(\vec a\times\vec b)\cdot\vec c$ と書ける(同じ記事の補題「行列式と外積の性質」)。この形 $(\vec a\times\vec b)\cdot\vec c$ を スカラー三重積 ともいう。
行列式の計算を楽にし、体積との関係を示すための性質をまとめる。
空間のベクトル $\vec a$、$\vec b$、$\vec c$、$\vec a'$ と実数 $s$、$t$ について、次が成り立つ。
方針:どの規則も、def-vtd-det の 6 つの項がどんな形をしているかから出る。4 は 1 と 3 から出す。
段 1(1 を示す)。6 つの項はどれも、$\vec a$ の成分をちょうど 1 つ因数にもつ。$\vec a$ を $s\vec a+t\vec a'$ に替えると、たとえば第 1 項は $(sa_1+ta_1')b_2c_3=s\,a_1b_2c_3+t\,a_1'b_2c_3$ と 2 つに分かれる。6 つの項すべてで同じように分けて集めると、$s\det(\vec a,\vec b,\vec c)+t\det(\vec a',\vec b,\vec c)$ になる。$\vec b$、$\vec c$ の位置についても、各項が $\vec b$ の成分、$\vec c$ の成分をちょうど 1 つずつもつので同じである。
段 2(2 を示す)。$\vec a$ と $\vec b$ を入れ替えた $\det(\vec b,\vec a,\vec c)$ は、def-vtd-det の式で文字 $a$ と $b$ を入れ替えたものなので
$$
b_1a_2c_3+b_2a_3c_1+b_3a_1c_2-b_3a_2c_1-b_2a_1c_3-b_1a_3c_2
$$
である。掛け算の順番を並べ替えると、正の項は $a_2b_1c_3$、$a_3b_2c_1$、$a_1b_3c_2$ で、これらはもとの式の負の項である。負の項は $a_2b_3c_1$、$a_1b_2c_3$、$a_3b_1c_2$ で、もとの式の正の項である。よって $\det(\vec b,\vec a,\vec c)=-\det(\vec a,\vec b,\vec c)$ である。ほかの 2 つの入れ替え($\vec b$ と $\vec c$、$\vec a$ と $\vec c$)も、同じように正の項と負の項がそっくり入れ替わる。
段 3(3 を示す)。たとえば $\vec a=\vec b$ とする。2 から $\det(\vec a,\vec a,\vec c)=-\det(\vec a,\vec a,\vec c)$ である(入れ替えても同じ組だから)。$D=-D$ なら $2D=0$ なので $D=0$ である。ほかの 2 つが等しい場合も同じである。
段 4(4 を示す)。1 と 3 により
$$
\det(\vec a+t\vec b,\vec b,\vec c)=\det(\vec a,\vec b,\vec c)+t\det(\vec b,\vec b,\vec c)=\det(\vec a,\vec b,\vec c)+t\cdot0
$$
である。足す相手や足す場所がほかの組み合わせでも、同じ計算になる。
段 5(5 を示す)。$\vec a=(p,x,y)$、$\vec b=(0,q,z)$、$\vec c=(0,0,r)$ では $b_1=c_1=c_2=0$ である。def-vtd-det の 6 つの項のうち、$a_1b_2c_3=pqr$ 以外の 5 つの項は、$c_1$、$b_1$、$c_2$ のどれかを因数にもつ($a_2b_3c_1$ と $a_3b_2c_1$ は $c_1$、$a_3b_1c_2$ と $a_2b_1c_3$ は $b_1$、$a_1b_3c_2$ は $c_2$)ので $0$ である。よって $\det=pqr$ である。$p=q=r=1$、$x=y=z=0$ とすると $\det(\vec e_1,\vec e_2,\vec e_3)=1$ である。$\square$
ex-vtd-det-calc の $\vec a=(2,1,0)$、$\vec b=(0,2,1)$、$\vec c=(1,1,3)$ の行列式を、lem-vtd-det-rules の 4 で上三角の形にしてから求める。
四面体 $\mathrm{OABC}$ で $\vec a=\overrightarrow{\mathrm{OA}}$、$\vec b=\overrightarrow{\mathrm{OB}}$、$\vec c=\overrightarrow{\mathrm{OC}}$ とすると、その体積は
$$
V=\frac16\,\bigl\lvert\det(\vec a,\vec b,\vec c)\bigr\rvert
$$
である。4 点 $\mathrm O$、$\mathrm A$、$\mathrm B$、$\mathrm C$ が同じ平面の上にある(四面体がつぶれている)ときは、体積を $0$ とみなせば、この式はそのまま成り立つ。頂点が一般の 4 点 $\mathrm A$、$\mathrm B$、$\mathrm C$、$\mathrm D$ のときは
$$
V=\frac16\,\bigl\lvert\det(\overrightarrow{\mathrm{AB}},\overrightarrow{\mathrm{AC}},\overrightarrow{\mathrm{AD}})\bigr\rvert
$$
である。
最後の式は、$\mathrm A$ が原点に来るように全体を平行移動すれば、はじめの式になる。平行移動で体積は変わらず、$\overrightarrow{\mathrm{AB}}$ などのベクトルも変わらないからである。つぶれた場合の主張は、4 点が同じ平面の上にあることと $\det=0$ が同値であること(空間ベクトルと1次独立 の定理「1 次独立と行列式」)を言い換えたものである。証明を 2 通り与える。1 つ目は高校の方法(ex-vtd-slanted)を文字で書いたもの、2 つ目は行列式の計算規則(lem-vtd-det-rules)と体積の性質を突き合わせるものである。
方針:ex-vtd-slanted と同じ手順を文字で行う。底面 $\mathrm{OAB}$ の面積と、$\mathrm C$ から平面 $\mathrm{OAB}$ までの距離を、外積 $\vec n=\vec a\times\vec b$ で表す。
段 0(つぶれた場合)。$\mathrm O$、$\mathrm A$、$\mathrm B$ が一直線上にあれば、$\vec a$、$\vec b$ は 1 次従属で、空間ベクトルと1次独立 の定理「1 次独立と行列式」の証明の段 1 と同じく $\det(\vec a,\vec b,\vec c)=0$ である。このとき 4 点は同じ平面の上にあり、体積も $0$ である。以下、$\vec a$、$\vec b$ は平行でないとする。
段 1(法線ベクトル)。$\vec n=\vec a\times\vec b=(a_2b_3-a_3b_2,\ a_3b_1-a_1b_3,\ a_1b_2-a_2b_1)$ とおく。$\vec n\cdot\vec a=0$、$\vec n\cdot\vec b=0$ で、$\vec a$、$\vec b$ が平行でないので $\vec n\ne\vec 0$ である(空間ベクトルと1次独立 の補題「行列式と外積の性質」の 3 と 5)。よって $\vec n$ は平面 $\mathrm{OAB}$ の法線ベクトルである。
段 2(底面積)。成分を展開して比べると
$$
\lvert\vec n\rvert^2=\lvert\vec a\rvert^2\lvert\vec b\rvert^2-(\vec a\cdot\vec b)^2
$$
が成り立つ(Lagrange の恒等式。証明は 空間の平面と法線ベクトル の定理「外積の性質」の 3 にあり、de Guaの定理 では $n$ 次元の形で展開して示す)。三角形の面積の公式 の $S=\dfrac12\sqrt{\lvert\vec a\rvert^2\lvert\vec b\rvert^2-(\vec a\cdot\vec b)^2}$ と合わせて、三角形 $\mathrm{OAB}$ の面積は $S=\dfrac12\lvert\vec n\rvert$ である。
段 3(高さ)。平面 $\mathrm{OAB}$ は原点を通り法線ベクトルが $\vec n$ なので、その方程式は $\vec n\cdot\vec x=0$ である。点と平面の距離の公式(空間の平面と法線ベクトル の定理「点と平面の距離」)により、$\mathrm C$ からこの平面までの距離は
$$
h=\frac{\lvert\vec n\cdot\vec c\rvert}{\lvert\vec n\rvert}
$$
である。
段 4(掛ける)。錐体の体積の公式に段 2・段 3 を代入すると
$$
V=\frac13Sh=\frac13\cdot\frac12\lvert\vec n\rvert\cdot\frac{\lvert\vec n\cdot\vec c\rvert}{\lvert\vec n\rvert}=\frac16\lvert\vec n\cdot\vec c\rvert
$$
である。$\vec n\cdot\vec c=(\vec a\times\vec b)\cdot\vec c=\det(\vec a,\vec b,\vec c)$ なので(空間ベクトルと1次独立 の補題「行列式と外積の性質」の 1)、$V=\dfrac16\lvert\det(\vec a,\vec b,\vec c)\rvert$ である。$\square$
1 つ目の証明は、$\lvert\vec n\rvert$ が分母と分子で打ち消し合うことを示した。これが ex-vtd-slanted で $\sqrt{21}$ が消えた理由である。2 つ目の証明では、外積を使わず、体積そのものの性質から出発する。鍵になるのは次の補題である。
頂点 A を、面 OBC の中の直線 OB に平行に A' まで動かしても、面 OBC からの高さ h は変わらない。底面 OBC も同じなので、体積は変わらない
四面体 $\mathrm{OABC}$ の体積を $V(\vec a,\vec b,\vec c)$ と書く(つぶれていれば $0$)。
方針:1 は四面体そのものが変わらないことから、2 は底面を面 $\mathrm{OBC}$ にとり、高さが変わらないことを法線ベクトルとの内積で確かめる(図 2)。
段 1(1 を示す)。並べ替えても、4 つの頂点 $\mathrm O$、$\mathrm A$、$\mathrm B$、$\mathrm C$ の集まりは同じなので、四面体も体積も同じである。
段 2(2 で、つぶれている場合)。$\mathrm O$、$\mathrm B$、$\mathrm C$ が一直線上にあれば、$\mathrm A$ をどこに動かしても 4 点は同じ平面の上にあり、体積は動かす前も後も $0$ である。
段 3(2 で、つぶれていない面を底面にとる)。$\mathrm O$、$\mathrm B$、$\mathrm C$ が一直線上にないとし、底面を三角形 $\mathrm{OBC}$ にとる。$\vec a$ を $\vec a'=\vec a+t\vec b$ に替えることは、頂点 $\mathrm A$ を $\mathrm A'$($\overrightarrow{\mathrm{OA'}}=\vec a'$)に動かすことで、底面 $\mathrm{OBC}$ は変わらない。
段 4(高さが変わらない)。平面 $\mathrm{OBC}$ の法線ベクトルを $\vec m$ とすると、平面 $\mathrm{OBC}$ は原点を通るので、点 $\mathrm X$ からこの平面までの距離は $\dfrac{\lvert\vec m\cdot\overrightarrow{\mathrm{OX}}\rvert}{\lvert\vec m\rvert}$ である(空間の平面と法線ベクトル の定理「点と平面の距離」)。$\vec b$ は平面 $\mathrm{OBC}$ の中の向きなので $\vec m\cdot\vec b=0$ であり、
$$
\vec m\cdot\vec a'=\vec m\cdot\vec a+t\,(\vec m\cdot\vec b)=\vec m\cdot\vec a
$$
である。よって $\mathrm A'$ の高さは $\mathrm A$ の高さと等しい。底面積も高さも同じなので、$\dfrac13\times$ 底面積 $\times$ 高さ は変わらない。
段 5(ほかの組み合わせ)。1 つのベクトル $\vec u$ に別のベクトル $\vec w$ の実数倍を足すときは、残りの 2 つのベクトル($\vec w$ を含む)で作る面を底面にとれば、段 3・段 4 と同じ議論になる。$\square$
補題の 2 は、lem-vtd-det-rules の 4「1 つのベクトルに別のベクトルの実数倍を足しても行列式は変わらない」と同じ形をしている。体積と行列式は、同じ操作で変わらない。これを使って、どんな四面体も直角をもつ四面体に変形する。
方針:lem-vtd-slide と lem-vtd-det-rules により、「順番の並べ替え」と「1 つのベクトルに別のベクトルの実数倍を足す」の 2 つの操作では、体積 $V$ も行列式の絶対値 $\lvert\det\rvert$ も変わらない(並べ替えで行列式は符号が変わることがあるが、絶対値は変わらない)。この 2 つの操作をくり返して、3 つのベクトルを座標軸の上に並べ、ex-vtd-corner の形にする。
段 1(第 1 成分をそろえる)。$a_1$、$b_1$、$c_1$ がすべて $0$ なら、3 つのベクトルはどれも平面 $x=0$ の中にあり、4 点 $\mathrm O$、$\mathrm A$、$\mathrm B$、$\mathrm C$ はこの平面の上にあるので $V=0$ である。一方、def-vtd-det の 6 つの項はどれも $a_1$、$b_1$、$c_1$ のどれかを因数にもつので $\det=0$ で、主張は成り立つ。そうでなければ、並べ替えて $a_1\ne0$ とし、$\vec b$ を $\vec b-\dfrac{b_1}{a_1}\vec a$ に、$\vec c$ を $\vec c-\dfrac{c_1}{a_1}\vec a$ に替える。すると新しい $\vec b$、$\vec c$ の第 1 成分は $0$ になる。
段 2(第 2 成分をそろえる)。いま $\vec b=(0,b_2,b_3)$、$\vec c=(0,c_2,c_3)$ である。$b_2=c_2=0$ なら、$\vec b$ と $\vec c$ はどちらも $z$ 軸の上にあり、$\mathrm O$、$\mathrm B$、$\mathrm C$ が一直線上にあるので $V=0$ である。このとき lem-vtd-det-rules の 5 で $q=0$ の場合にあたり $\det=0$ で、主張は成り立つ。そうでなければ、$\vec b$ と $\vec c$ を並べ替えて $b_2\ne0$ とし、$\vec c$ を $\vec c-\dfrac{c_2}{b_2}\vec b$ に替える。すると $\vec c=(0,0,c_3)$ になる。
段 3(上三角の形)。いま $\vec a=(a_1,a_2,a_3)$、$\vec b=(0,b_2,b_3)$、$\vec c=(0,0,c_3)$ で、$a_1\ne0$、$b_2\ne0$ である。$c_3=0$ なら $\mathrm C=\mathrm O$ なので $V=0$ で、lem-vtd-det-rules の 5 により $\det=a_1b_2\cdot0=0$ で、主張は成り立つ。
段 4(座標軸の上に並べる)。$c_3\ne0$ とする。$\vec b$ を $\vec b-\dfrac{b_3}{c_3}\vec c=(0,b_2,0)$ に替え、次に $\vec a$ を $\vec a-\dfrac{a_3}{c_3}\vec c-\dfrac{a_2}{b_2}\vec b$ に替えると $(a_1,0,0)$ になる(いまの $\vec b$ は $(0,b_2,0)$ である)。3 つのベクトルは $(a_1,0,0)$、$(0,b_2,0)$、$(0,0,c_3)$ になった。
段 5(直角をもつ四面体)。段 4 の四面体は、3 本の辺が座標軸の上にあって互いに垂直で、長さは $\lvert a_1\rvert$、$\lvert b_2\rvert$、$\lvert c_3\rvert$ である。ex-vtd-corner と同じく、底面積 $\dfrac12\lvert a_1\rvert\lvert b_2\rvert$、高さ $\lvert c_3\rvert$ なので $V=\dfrac16\lvert a_1b_2c_3\rvert$ である(空間図形の計量 の命題「直角をもつ四面体の面積と高さ」の 1)。一方 lem-vtd-det-rules の 5 により $\det=a_1b_2c_3$ である。よって $V=\dfrac16\lvert\det\rvert$ が成り立つ。
段 6(もとに戻す)。段 1〜段 4 の操作で $V$ も $\lvert\det\rvert$ も変わらなかったので、もとの 3 つのベクトルでも $V=\dfrac16\lvert\det(\vec a,\vec b,\vec c)\rvert$ である。$\square$
ex-vtd-slanted の $\vec a=(2,1,0)$、$\vec b=(0,2,1)$、$\vec c=(1,1,3)$ で、prf-vtd-volume-reduce の手順をたどる。
$\mathrm A(1,0,1)$、$\mathrm B(3,1,1)$、$\mathrm C(1,2,2)$、$\mathrm D(2,1,4)$ を頂点とする四面体の体積を求める。$\overrightarrow{\mathrm{AB}}=(2,1,0)$、$\overrightarrow{\mathrm{AC}}=(0,2,1)$、$\overrightarrow{\mathrm{AD}}=(1,1,3)$ で、これは ex-vtd-slanted の 3 つのベクトルと同じである(この四面体は ex-vtd-slanted の四面体を $(1,0,1)$ だけ平行移動したもの)。ex-vtd-det-calc により $\det=11$ なので、thm-vtd-volume から $V=\dfrac{11}6$ である。
thm-vtd-volume では行列式に絶対値を付けた。絶対値を外した値にも意味がある。
空間のベクトル $\vec a$、$\vec b$、$\vec c$ について、
$$
[\vec a,\vec b,\vec c]=\frac16\det(\vec a,\vec b,\vec c)
$$
を、四面体 $\mathrm{OABC}$ の(この順番の)符号つき体積 という。$\det(\vec a,\vec b,\vec c)>0$ のとき、3 つのベクトルの組 $(\vec a,\vec b,\vec c)$ は 右手系 である(基本ベクトルの組 $(\vec e_1,\vec e_2,\vec e_3)$ と同じ向きである)といい、$\det(\vec a,\vec b,\vec c)<0$ のとき 左手系 であるという。
$x$ 軸、$y$ 軸、$z$ 軸の正の向きを、右手の親指、人差し指、中指に合わせられるように座標をとるのが普通である。そのとき右手系の組は、右手の 3 本の指に(指の開きを変えながら)合わせられる組である。この記事では、この直観の説明には立ち入らず、行列式の符号を右手系の定義とする。
符号つき体積の良いところは、lem-vtd-det-rules の 1 により、各ベクトルについて 1 次式になることである。そのおかげで、次の体積比が一度に出る。
四面体 $\mathrm{OABC}$ で、半直線 $\mathrm{OA}$、$\mathrm{OB}$、$\mathrm{OC}$ の上($\mathrm O$ 以外)に点 $\mathrm P$、$\mathrm Q$、$\mathrm R$ をとる。四面体 $\mathrm{OPQR}$、$\mathrm{OABC}$ の体積をそれぞれ $V'$、$V$ とすると
$$
\frac{V'}{V}=\frac{\mathrm{OP}}{\mathrm{OA}}\cdot\frac{\mathrm{OQ}}{\mathrm{OB}}\cdot\frac{\mathrm{OR}}{\mathrm{OC}}
$$
である。
方針:$\overrightarrow{\mathrm{OP}}$ などを $\vec a$ などの正の実数倍で書き、行列式の各行の 1 次性で倍率をくくり出す。
段 1(実数倍で書く)。$\mathrm P$ は半直線 $\mathrm{OA}$ の上にあるので、$s=\dfrac{\mathrm{OP}}{\mathrm{OA}}>0$ として $\overrightarrow{\mathrm{OP}}=s\vec a$ である。同じく $t=\dfrac{\mathrm{OQ}}{\mathrm{OB}}$、$u=\dfrac{\mathrm{OR}}{\mathrm{OC}}$ として $\overrightarrow{\mathrm{OQ}}=t\vec b$、$\overrightarrow{\mathrm{OR}}=u\vec c$ である。
段 2(くくり出す)。lem-vtd-det-rules の 1 を $\vec a$ の位置、$\vec b$ の位置、$\vec c$ の位置で 1 回ずつ使うと
$$
\det(s\vec a,t\vec b,u\vec c)=s\det(\vec a,t\vec b,u\vec c)=st\det(\vec a,\vec b,u\vec c)=stu\det(\vec a,\vec b,\vec c)
$$
である。
段 3(体積の比)。thm-vtd-volume により $V'=\dfrac16\lvert stu\det(\vec a,\vec b,\vec c)\rvert=stu\cdot\dfrac16\lvert\det(\vec a,\vec b,\vec c)\rvert=stu\,V$ である($s,t,u>0$)。四面体 $\mathrm{OABC}$ はつぶれていないので $V>0$ で、両辺を $V$ で割れば主張が出る。$\square$
これは、三角形の面積比「$\dfrac{\triangle\mathrm{OPQ}}{\triangle\mathrm{OAB}}=\dfrac{\mathrm{OP}}{\mathrm{OA}}\cdot\dfrac{\mathrm{OQ}}{\mathrm{OB}}$」(三角形の面積比)の空間版である。平面の場合も、2 次の行列式の各行の 1 次性から同じように出る。
thm-vtd-volume の分母の $6$ は、証明の中では「$\dfrac13\times\dfrac12$」として現れた。$\dfrac13$ は錐体の体積の公式から、$\dfrac12$ は三角形の面積から来ている。ここでは、$6$ を「3 つの数の並べ方の個数 $3!$」として見る。
1 辺 1 の立方体を、3 つの座標の大小の順(6 通り)で 6 つの四面体に分ける。どれも体積は 6 分の 1 である
1 辺 $1$ の立方体 $0\le x\le1$、$0\le y\le1$、$0\le z\le1$ は、3 つの座標の大小の順で分けた 6 つの部分
$$
1\ge x\ge y\ge z\ge0,\quad 1\ge x\ge z\ge y\ge0,\quad 1\ge y\ge x\ge z\ge0,
$$
$$
1\ge y\ge z\ge x\ge0,\quad 1\ge z\ge x\ge y\ge0,\quad 1\ge z\ge y\ge x\ge0
$$
に分かれる(図 3)。どの部分も四面体で、体積は $\dfrac16$ である。たとえば $1\ge x\ge y\ge z\ge0$ の部分は、4 点 $(0,0,0)$、$(1,0,0)$、$(1,1,0)$、$(1,1,1)$ を頂点とする四面体である。
方針:$1\ge x\ge y\ge z\ge0$ の部分が四面体であることを、点を 3 つのベクトルの 1 次結合で表して示す。体積は thm-vtd-volume で求める。ほかの 5 つは座標の名前を入れ替えれば同じである。
段 1(1 次結合で表す)。$\vec u=(1,0,0)$、$\vec v=(1,1,0)$、$\vec w=(1,1,1)$ とおく。点 $(x,y,z)$ について
$$
(x-y)\vec u+(y-z)\vec v+z\vec w=(x-y+y-z+z,\ y-z+z,\ z)=(x,y,z)
$$
である。係数 $s=x-y$、$t=y-z$、$u=z$ は、$x\ge y\ge z\ge0$ のときどれも $0$ 以上で、和は $s+t+u=x\le1$ である。逆に、$s,t,u\ge0$、$s+t+u\le1$ なら、$(x,y,z)=s\vec u+t\vec v+u\vec w=(s+t+u,\ t+u,\ u)$ は $1\ge x\ge y\ge z\ge0$ を満たす。
段 2(四面体である)。$s,t,u\ge0$、$s+t+u\le1$ を満たす点 $s\vec u+t\vec v+u\vec w$ の全体は、原点と $\vec u$、$\vec v$、$\vec w$ の終点を頂点とする四面体(の内部と表面)である(点の存在範囲と凸結合)。段 1 から、$1\ge x\ge y\ge z\ge0$ の部分はこの四面体である。
段 3(体積)。$\det(\vec u,\vec v,\vec w)$ の 6 つの項のうち $0$ でないのは $u_1v_2w_3=1\cdot1\cdot1=1$ だけである($u_2=u_3=v_3=0$ なので、ほかの項は $u_2$、$u_3$、$v_3$ のどれかを因数にもつ)。よって $\det=1$ で、thm-vtd-volume により体積は $\dfrac16$ である。
段 4(6 つで立方体)。立方体のどの点でも、3 つの座標 $x$、$y$、$z$ を大きい順に並べる並べ方は少なくとも 1 つあるので、どれかの部分に入る。2 つの部分に同時に入るのは、2 つの座標が等しい点(平面 $x=y$、$y=z$、$z=x$ の上の点)だけで、そこは体積をもたない。6 つの部分の体積の和 $6\cdot\dfrac16=1$ は立方体の体積と一致する。$\square$
段 4 を逆に読むと、thm-vtd-volume を使わずに $\dfrac16$ が出る。$x$ と $y$ を入れ替える操作は平面 $x=y$ に関する鏡映(折り返し)で、図形を合同な図形に移す。座標の名前の入れ替えをくり返せば、6 つの部分はどれもほかのどれにも移るので、6 つとも合同で体積が等しい。和が $1$ なので、1 つは $\dfrac16$ である。分母の $6$ は、3 つの座標の並べ方の個数 $3!=6$ である。
| 外した仮定 | 崩れる主張 | ボックス |
|---|---|---|
| 行列式に絶対値を付けること | 式の値が体積($0$ 以上)になる | ex-vtd-negative |
| 3 本のベクトルが 1 つの頂点から出る 3 辺であること | 行列式の絶対値の $\dfrac16$ が四面体の体積になる | ex-vtd-face-edges |
| 原点が四面体の頂点であること | 頂点の位置ベクトル 3 つの行列式の絶対値の $\dfrac16$ が体積になる | ex-vtd-origin |
ex-vtd-signed の (2) のとおり、ex-vtd-slanted の四面体で順番を $(\vec b,\vec a,\vec c)$ にすると、$\dfrac16\det(\vec b,\vec a,\vec c)=-\dfrac{11}6$ で負になる。体積は $\dfrac{11}6$ である。行列式の符号は四面体の形ではなく、3 つのベクトルを並べる順番(向き)で決まる。
ex-vtd-four-points の四面体 $\mathrm{ABCD}$(体積 $\dfrac{11}6$)で、$\overrightarrow{\mathrm{AB}}=(2,1,0)$、$\overrightarrow{\mathrm{AC}}=(0,2,1)$、$\overrightarrow{\mathrm{BC}}=(-2,1,1)$ の行列式を計算すると、6 つの項は $2\cdot2\cdot1=4$、$1\cdot1\cdot(-2)=-2$、$0$、$0$、$1\cdot0\cdot1=0$、$2\cdot1\cdot1=2$ なので、$\det=4-2+0-0-0-2=0$ である。3 辺はどれも面 $\mathrm{ABC}$ の中にあり、$\overrightarrow{\mathrm{BC}}=\overrightarrow{\mathrm{AC}}-\overrightarrow{\mathrm{AB}}$ なので 1 次従属だからである。
一方、$\overrightarrow{\mathrm{AB}}$、$\overrightarrow{\mathrm{BC}}$、$\overrightarrow{\mathrm{CD}}$ のように「つながった 3 辺」をとると、$\overrightarrow{\mathrm{BC}}=\overrightarrow{\mathrm{AC}}-\overrightarrow{\mathrm{AB}}$、$\overrightarrow{\mathrm{CD}}=\overrightarrow{\mathrm{AD}}-\overrightarrow{\mathrm{AC}}$ なので、lem-vtd-det-rules の 4 を 2 回使って $\det(\overrightarrow{\mathrm{AB}},\overrightarrow{\mathrm{BC}},\overrightarrow{\mathrm{CD}})=\det(\overrightarrow{\mathrm{AB}},\overrightarrow{\mathrm{AC}},\overrightarrow{\mathrm{AD}})=11$ となり、正しい体積が出る。どの 3 辺でもよいわけではなく、3 辺が同じ平面の中にあると $0$ になる。
ex-vtd-four-points の $\mathrm A(1,0,1)$、$\mathrm B(3,1,1)$、$\mathrm C(1,2,2)$、$\mathrm D(2,1,4)$ で、原点から $\mathrm B$、$\mathrm C$、$\mathrm D$ への位置ベクトル $(3,1,1)$、$(1,2,2)$、$(2,1,4)$ の行列式は、6 つの項が $3\cdot2\cdot4=24$、$1\cdot2\cdot2=4$、$1\cdot1\cdot1=1$、$1\cdot2\cdot2=4$、$1\cdot1\cdot4=4$、$3\cdot2\cdot1=6$ なので、$24+4+1-4-4-6=15$ である。$\dfrac{15}6$ は四面体 $\mathrm{OBCD}$ の体積で、四面体 $\mathrm{ABCD}$ の体積 $\dfrac{11}6$ ではない。thm-vtd-volume の 1 つ目の式は、原点 $\mathrm O$ が四面体の頂点であることを使っている。頂点が一般の 4 点のときは、1 つの頂点から出る辺のベクトル $\overrightarrow{\mathrm{AB}}$、$\overrightarrow{\mathrm{AC}}$、$\overrightarrow{\mathrm{AD}}$ を使う。
prf-vtd-volume-reduce で使ったのは、体積 $V$ と $\dfrac16\lvert\det\rvert$ に共通する次の 3 つの性質だけである。(a) 3 本のベクトルの順番を並べ替えても変わらない。(b) 1 つのベクトルに別のベクトルの実数倍を足しても変わらない。(c) 座標軸の上に並んだ 3 本のベクトルでの値(ex-vtd-corner)。体積 $V$ は向きを考えない量なので、各ベクトルについて 1 次ではない(たとえば $\vec a$ を $-\vec a$ に替えても $V$ は変わらない)。行列式の側では、(a)〜(c) は行列式を決める次の 3 つの性質から出る。(i) 各ベクトルについて 1 次(lem-vtd-det-rules の 1)。(ii) 2 つのベクトルが等しければ $0$(同じ 3)。(iii) 基本ベクトルでの値が $\det(\vec e_1,\vec e_2,\vec e_3)=1$(同じ 5)。大学の線形代数では、(i)・(ii) を満たす関数を 交代 3 重線形形式 という。$\mathbb{R}^3$ の上の交代 3 重線形形式は、どれも $\det$ の定数倍であり(Axl24 の 9.37 の次元の主張、p. 351)、その定数は基本ベクトルでの値で決まる。$n$ 次の行列式も、$n$ 本のベクトルについての交代 $n$ 重線形形式である(Axl24 の 9.45、p. 355)。「体積の性質を満たすものは行列式しかない」ので、行列式は体積の性質を式にしたものと見られる。平面の場合の同じ議論は 1次変換と行列式:面積の拡大率 の補題「双線形で交代的な関数は符号付き面積の定数倍」にある。
$n$ 次元の空間 $\mathbb{R}^n$ で、$n+1$ 個の点 $\mathrm P_0,\dots,\mathrm P_n$ の凸結合の全体を $n$ 単体 といい(単体複体 を参照)、三角形($n=2$)と四面体($n=3$)の一般化である。$n$ 単体の体積は
$$
V=\frac1{n!}\bigl\lvert\det(\overrightarrow{\mathrm P_0\mathrm P_1},\dots,\overrightarrow{\mathrm P_0\mathrm P_n})\bigr\rvert
$$
である。$\dfrac1{n!}$ の理由は prop-vtd-cube-six と同じで、1 辺 $1$ の $n$ 次元立方体が、$n$ 個の座標の大小の順で合同な $n!$ 個の単体に分かれることから、基本ベクトルの単体の体積が $\dfrac1{n!}$ になる。一般の単体は、基本ベクトルを $\overrightarrow{\mathrm P_0\mathrm P_i}$ に移す 1 次変換で基本の単体を写したものであり、1 次変換で体積は $\lvert\det\rvert$ 倍になる(Axl24 の 9.61、p. 363。この記事では証明しない)。$n=2$ が三角形の面積の $\dfrac12$、$n=3$ が thm-vtd-volume の $\dfrac16$ である。
平面の多角形の面積を頂点の座標から求める 靴ひも公式 には、空間版がある。多面体の表面を三角形に分け、各三角形の頂点を多面体の外から見て反時計回りに $\mathrm P$、$\mathrm Q$、$\mathrm R$ と並べると、多面体の体積は、すべての三角形についての $\dfrac16\det(\overrightarrow{\mathrm{OP}},\overrightarrow{\mathrm{OQ}},\overrightarrow{\mathrm{OR}})$ の和である($\mathrm O$ は原点。一般の多面体での証明はこの記事ではしない)。各項は、原点と 1 つの三角形でできる四面体の符号つき体積(def-vtd-signed)で、原点から見て裏側の面の項は負になって、多面体の外の部分が打ち消される。ex-vtd-four-points の四面体 $\mathrm{ABCD}$ で 4 つの面についてこの和を計算すると、$\dfrac{11}6$ になる。
この記事の式は、頂点の座標から体積を求める。6 本の辺の長さだけから体積を求める式もあり、Cayley–Menger 行列式と呼ばれる $5\times5$ の行列式で $288V^2$ が表される。三角形の面積を 3 辺で表す Heron の公式の空間版で、Heronの公式とCayley–Menger行列式 で扱う。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する