― 講義ノート「(16) 行列式が 0 であることと 0 でないことは何を意味するのか?」の導出
行列式は,覚えるべき計算手順として教わってしまうことが多い量です. でも本当は,たった一言で言えます ―― 行列式は,ベクトルたちが張る図形の「大きさ」です. 2 次元なら面積,3 次元なら体積.そう思って眺めると, det = 0 も,線形独立も,rank も,連立方程式の解の形も, ぜんぶ同じ 1 枚の絵の別の言い方になります.
講義ノートの例をそのまま追います. , を列に並べて行列を作ると
一方,この 2 本が張る平行四辺形の面積を幾何で求めます. 全体を囲む長方形(横 4,縦 5)から,まわりの 2 つの長方形(合計 4)と 4 つの三角形(合計 10)を引くと
(1) = (2).行列式は,平行四辺形の面積そのものでした.
, として,同じことをやります.平行四辺形の面積は「底辺 × 高さ」= です. ここで と を使うと,
展開すると と が消え,
図 1.長方形から余りを引くと,行列式が残る. 講義ノートの数え方をそのまま絵にしたものです. 外側の長方形から,まわりの三角形と長方形を引くと平行四辺形の面積になり, それが と一致します.
実は,計算しなくても行列式の形は決まってしまいます. 「面積」に当然あってほしい性質を書き出すだけでよいのです.
| 性質 | 式 | 面積として当たり前なこと |
|---|---|---|
| 線形性 | ( についても同じ) |
底辺を 2 倍に伸ばせば,面積も 2 倍.底辺を 2 本に分ければ,面積も 2 つの和 |
| 交代性 | 同じ向きの 2 本では,平行四辺形がつぶれる | |
| 規格化 | 単位正方形の面積は 1 |
この 3 つだけで十分です. , を入れて線形性でばらすと 4 項出ますが,交代性のおかげで と が消え, さらに交代性から出る (下の囲みを見てください)でまとめると,
同じ 3 つの性質を 3 本のベクトルに拡張すると,出てくるのは
この形は,体積の定義そのままです. が底面の平行四辺形の面積, は底面に垂直なので, をその向きに射影した長さが高さ. 掛ければ底面積 × 高さ = 体積です. ただし (f × g)·h は符号つきの量で,体積そのものはその絶対値 |(f × g)·h| です. 符号は右手系か左手系かを表します.
講義ノートの例で確かめます. , , :
ここまでをまとめると,次の言い換えがすべて同じ意味になります (幾何の欄は n = 3,つまり 3 次元の場合です.2 次元なら「平面」を「直線」,「直線」を「原点」と読み替えてください).
| 幾何の言葉 | 代数の言葉 | 連立方程式の言葉 |
|---|---|---|
| 面積(体積)が 0 でない | det ≠ 0,線形独立,rank = n | 解は原点だけ(自明な解のみ) |
| つぶれて平面になった | det = 0,rank = n − 1 | 解が直線になる(自由度 1) |
| つぶれて直線になった | det = 0,rank = n − 2 | 解が平面になる(自由度 2) |
同次連立方程式 Mx = 0 を考えます.講義ノートの例②:
3 本目の式は,1 本目・2 本目から作れてしまいます(消去すると 0 = 0 になる). つまり本当に効いている式は 2 本だけ.未知数は 3 つあるので,
と置けば , . 解は x = k(−1, 1, 1) ―― 原点を通る 1 本の直線まるごとです. 例③(係数行列を M′ と書きます)では 3 本の式がすべて同じ式の定数倍なので rank = 1,自由度 2 で, 解は 2x1 − x2 + x3 = 0 という平面になります.
行列を「変換」と見ると,|det| は面積が何倍になるかの倍率です. だから 2 つ続けて変換すれば,面積の倍率は掛け算になり,向きが裏返るかどうかも掛け合わされます ―― それが です. 行列の積の定義から力ずくで証明することもできますが,面積の倍率だと思えば当たり前です.
図 2.行列は「面積を何倍にするか」. 黄色い単位正方形(面積 1)が B で青い平行四辺形に, さらに A で緑の平行四辺形に写ります. 面積は 倍, 倍と 2 段で変わり, 合わせて 倍になります.
ここまで来ると,永年方程式の意味がそのまま読めます.分子軌道の問題では
という同次連立方程式が出てきます. c = 0 では波動関数が消えてしまうので,自明でない解が必要です. そのためには「つぶれていなければならない」―― つまり
この条件から ε(エネルギー準位)が決まります. 「行列式を 0 と置く」は計算のおまじないではなく, 「解が原点だけにならないように,空間をわざとつぶす」という幾何の要請だったのです. 実際に動かして確かめたい方は Lagrange の未定乗数法のページへどうぞ.
シミュレーターに載っているコードそのものを走らせた結果です (32 項目,すべて通過).
| 確かめたこと | 結果 |
|---|---|
| (1)(2):講義ノートの |A| = 6 と,幾何で数えた面積 6 が一致するか | 両方 6.000000000000 |
| (3):|det| = 底辺 × 高さ(ランダムなベクトル 500 組) | 最大差 1.8×10−15 |
| (5):|det| = 底面積 × 高さ(3 次元,ランダム 500 組) | 最大差 1.4×10−14 |
| (6):講義ノートの |C| = 0,rank = 2,h = 0.6f + 0.2g | det = 0,法線 (5, 5, −5) に 3 本とも垂直 |
| (4):列の入れ替えで符号反転,同じ列で 0,単位行列で 1,回転で 1 | すべて厳密に成立 |
| §6:det(AB) = det A · det B | 5 = 5 × 1(厳密) |
| (7):例② rank 2 → 解 x = k(−1, 1, 1),M·x = 0 | 残差 0(厳密) |
| 例③ rank 1 → 自由度 2,解は 2x₁ − x₂ + x₃ = 0 の平面 | 基底 2 本とも M′·v = 0 |
| det ≠ 0 の例:rank 3,解空間の次元 0 | x = 0 のみ |
| 次元定理 rank + 自由度 = 3(ランダム 300 通り) | 不一致 0 件 |
| det = 0 ⟺ rank < 3(det = 0 の例を混ぜた 400 通り,うち 220 件が det = 0) | 食い違い 0 件 |
東京理科大学 望月研究室 / 講義ノート「(16) 行列式が 0 であることと 0 でないことは何を意味するのか?」より. §2 の「3 つの性質から一意に決まる」話と §6 の det(AB) は,ノートの外を補ったものです. 数値検証は倍精度の直接計算と,部分ピボット選択つき Gauss の消去法による. 動かして確かめたい方は シミュレーターへ. 数式はブラウザ標準の MathML で組んでいます(外部ライブラリ・通信は一切ありません).