非直交基底の成分分解

(90d) 更新


非直交基底の成分分解

$n$ 次元数ベクトル空間 $K^n$ に、 必ずしも直交せず、正規化もされていない基底 $\{\bm b_i\}$ を取る。

基底であるから任意のベクトル $\bm x$ を成分分解して、

$$\bm x=\sum_{i=1}^n x_i\bm b_i\hspace{3cm}(*)$$

と表せるのであるが、この係数 $\{x_i\}$ をどうやって求めたらよいかというのがここでの問題。

もしも $\{\bm b_i\}$ が正規直交基底であったなら、$x_i=(\bm b_i, \bm x)$ として簡単に求まるのであるが、これは、

$$ \begin{aligned} (\bm b_i,\bm x)&=\sum_{j=1}^n x_j(\bm b_i,\bm b_j)\\ &=\sum_{j=1}^n x_j\delta_{ij}=x_i\\ \end{aligned} $$

すなわち、

$$(\bm b_i,\bm b_j)=\delta_{ij}$$

となることをうまく使っているのだった(他の項がゼロになり、求めたい項の係数が1になる)。

$\{\bm b_i\}$ が必ずしも正規直交でない場合にはこれが成り立たないため困るののだけれど、もし代わりに

$$(\bm \beta_i,\bm b_j)=\delta_{ij}$$

を満たすような $n$ 本のベクトル $\bm\beta_i$ を発見できれば、やはり $(\bm \beta_i, \bm x)=x_i$ とするだけで成分を求められる! というのが以下の道筋となる。

目次

逆基底

基底 $\{\bm b_i\}$ の「逆基底(仮)」をここでは $\{\bm \beta_i\}$ と書くことにする。任意の $i,j$ に対して

$$(\bm \beta_j,\bm b_i)=\delta_{ij}$$

を満たすことが逆基底(仮)の定義である。下で解説する通り、数学的には $\bm \beta_j$ のエルミート共役である $\bm \beta_j^\dagger$ が本来の逆基底であるが、ここでは線形代数 II からの導入のしやすさを優先して「逆基底(仮)」を定義しておく。

つまり、$\bm \beta_j$ は

  • $(\bm \beta_j,\bm b_j)=1$
  • $\bm b_j$ 以外と直交

つまり、$\bm b_j$ 以外で張られる空間の直交補空間の元で、その $\bm b_j$ 方向の成分が $1/\|\bm b_j\|$ であるようなベクトルである。そして、

inverse_basis.png
  • 一般に $\bm b_j$ と平行ではない
    • $\bm b_j$ が $\bm b_j$ 以外の基底と直交しないのだから当然
  • $\{\bm b_i\}$ が正規直交なら $\bm \beta_j=\bm b_j$ と取れる。
    • 一般には $\bm b_j$ に $|k|=1$ を満たす係数 $k$ を掛けて $\bm \beta_j=k\bm b_j$ と取ってよい。

たとえば $n=3$ のとき、$\bm\beta_1$ は $\bm b_2,\bm b_3$ が張る空間 $V_{2,3}$ の直交補空間 $V_{2,3}^\perp$ のベクトルで、その $\bm b_1$ 方向成分の長さが $1/\|\bm b_1\|$ であるようなものだ。

成分を求めるには

上述のとおり、逆基底(仮)が求まりさえすれば、

$$(\bm \beta_j,\bm x)=x_j$$

のように簡単に成分を求められる。

逆基底の求め方(逆行列を使う)

$$(\bm \beta_i,\bm b_j)={\bm \beta_i}^\dagger\bm b_i=\delta_{ij}$$

であるから、${\bm \beta_i}^\dagger$ を並べた行列 $\mathcal B^\dagger$ と、 $\bm b_i$ を並べた行列 $B$ とを掛け合わせれば、

$$ \mathcal B^\dagger B= \begin{pmatrix} \hspace{7mm}\bm \beta_1^\dagger\hspace{7mm}\\ \bm \beta_2^\dagger\\ \bm \beta_3^\dagger\\ \vdots\\ \end{pmatrix} \begin{pmatrix} \rule{0mm}{7mm}\\ \bm b_1&\bm b_2&\bm b_3&\cdots\\ \rule{0mm}{7mm}\\ \end{pmatrix}= \begin{pmatrix} (\bm \beta_1,\bm b_1)&(\bm \beta_1,\bm b_2)&\cdots&&\\ (\bm \beta_2,\bm b_1)&(\bm \beta_2,\bm b_2)\\ \vdots&&\ddots\\ \\ \end{pmatrix}= \begin{pmatrix} \delta_{11}&\delta_{12}&\cdots&&\\ \delta_{21}&\delta_{22}&\\ \vdots&&\ddots&\\ \\ \end{pmatrix}=E $$

のように単位行列となる。すなわち、

$$ \mathcal B^\dagger=B^{-1}\hspace{1cm}\text{つまり}\hspace{1cm} \begin{pmatrix} \hspace{7mm}{\bm \beta_1}^\dagger\hspace{7mm}\\ {\bm \beta_2}^\dagger\\ {\bm \beta_3}^\dagger\\ \vdots\\ \end{pmatrix}= \begin{pmatrix} \rule{0mm}{7mm}\\ \bm b_1&\bm b_2&\bm b_3&\cdots\\ \rule{0mm}{7mm}\\ \end{pmatrix}^{-1} $$

の関係がある。

つまり、基底ベクトルを列ベクトルに持つ行列の逆行列を求め、 その列ベクトルのエルミート共役を取れば逆基底ベクトルを求められる。
(基底ベクトルは一次独立であるから $B$ には必ず逆行列が存在し、逆基底ベクトルも必ず存在する)

というかまあカバーストーリーをまるっとすっ飛ばせば、

$$\bm x=\sum_{i=1}^n x_i\bm b_i\hspace{3cm}(*)$$

という式は、

$$ \bm x=\underbrace{\Bigg(\bm b_1\ \bm b_2\ \dots\ \bm b_n\Bigg)}_B\begin{pmatrix} x_1\\x_2\\\vdots\\x_n \end{pmatrix} $$

を表すから、

$$ \begin{pmatrix} x_1\\x_2\\\vdots\\x_n \end{pmatrix}= B^{-1}\bm x=\underbrace{\begin{pmatrix} \hspace{5mm}\bm \beta_1^\dagger\hspace{5mm}\\ \bm \beta_1^\dagger\\ \vdots\\ \bm \beta_n^\dagger\\ \end{pmatrix}}_{B^{-1}}\bm x $$

と書け、成分である $x_i$ を $B^{-1}$ の行ベクトル $\bm \beta_i^\dagger$ を使って

$$ x_i=\bm \beta_i^\dagger\bm x $$

と書け、これを内積で書けば

$$ x_i=(\bm \beta_i,\bm x) $$

となる、ということだ。

次の4つのベクトルは $\mathbb{R}^4$ の基底となる。

$$ \begin{pmatrix} 1\\2\\3\\4 \end{pmatrix}, \begin{pmatrix} 2\\2\\-1\\-1 \end{pmatrix}, \begin{pmatrix} 1\\1\\1\\1 \end{pmatrix}, \begin{pmatrix} 0\\1\\1\\0 \end{pmatrix} $$

これらを列ベクトルとする行列を作り、その逆行列を求めよう。

$$ \begin{aligned} &\bigg(\ B\ \ \ E\ \ \bigg)=\begin{pmatrix} 1&2&1&0&1&0&0&0\\ 2&2&1&1&0&1&0&0\\ 3&-1&1&1&0&0&1&0\\ 4&-1&1&0&0&0&0&1\\ \end{pmatrix}\sim \begin{pmatrix} 1&2&1&0&1&0&0&0\\ 0&-2&-1&1&-2&1&0&0\\ 0&-7&-2&1&-3&0&1&0\\ 0&-9&-3&0&-4&0&0&1\\ \end{pmatrix}\\&\sim \begin{pmatrix} 1&2&1&0&1&0&0&0\\ 0&-2&-1&1&-2&1&0&0\\ 0&-14&-4&2&-6&0&2&0\\ 0&-2&-1&-1&-1&0&-1&1\\ \end{pmatrix}\sim \begin{pmatrix} 1&0&0&1&-1&1&0&0\\ 0&-2&-1&1&-2&1&0&0\\ 0&0&3&-5&8&-7&2&0\\ 0&0&0&-2&1&-1&-1&1\\ \end{pmatrix}\\&\sim \begin{pmatrix} 2&0&0&2&-2&2&0&0\\ 0&-6&-3&3&-6&3&0&0\\ 0&0&3&-5&8&-7&2&0\\ 0&0&0&-2&1&-1&-1&1\\ \end{pmatrix}\sim \begin{pmatrix} 2&0&0&2&-2&2&0&0\\ 0&-6&0&-2&2&-4&2&0\\ 0&0&6&-10&16&-14&4&0\\ 0&0&0&-2&1&-1&-1&1\\ \end{pmatrix}\\&\sim \begin{pmatrix} 2&0&0&0&-1&1&-1&1\\ 0&-6&0&0&1&-3&3&-1\\ 0&0&6&0&11&-9&9&-5\\ 0&0&0&-2&1&-1&-1&1\\ \end{pmatrix}\sim\bigg(\ E\ \ \ \mathcal B^\dagger\ \bigg) \end{aligned} $$

ここから、

$$ \bm \beta_1=\frac{1}{2}\begin{pmatrix}-1\\1\\-1\\1\end{pmatrix}, \ \bm \beta_2=\frac{1}{6}\begin{pmatrix}-1\\3\\-3\\1\end{pmatrix}, \ \bm \beta_3=\frac{1}{6}\begin{pmatrix}11\\-9\\9\\-5\end{pmatrix}, \ \bm \beta_4=\frac{1}{2}\begin{pmatrix}-1\\1\\1\\-1\end{pmatrix} $$

を得る。

例えば、$\bm x=\begin{pmatrix}-1\\ -4\\ 7\\ 14\end{pmatrix}$ と置くと、

$$ \begin{aligned} \bm x=&\sum_{i=1}^4(\bm\beta_i,\bm x)\bm b_i\\ =&\left[ \frac{1}{2} {\small \begin{pmatrix}-1\\1\\-1\\1 \end{pmatrix}\cdot \begin{pmatrix}-1\\ -4\\ 7\\ 14\end{pmatrix} }\right] \begin{pmatrix}1\\2\\3\\4\end{pmatrix}+ \left[ \frac{1}{6} {\small \begin{pmatrix}-1\\3\\-3\\1 \end{pmatrix}\cdot \begin{pmatrix}-1\\ -4\\ 7\\ 14\end{pmatrix} }\right] \begin{pmatrix}2\\2\\-1\\-1\end{pmatrix}\\ &+ \left[ \frac{1}{6} {\small \begin{pmatrix} 11\\-9\\9\\-5 \end{pmatrix}\cdot \begin{pmatrix}-1\\ -4\\ 7\\ 14\end{pmatrix} }\right] \begin{pmatrix}1\\1\\1\\1\end{pmatrix}+ \left[ \frac{1}{2} {\small \begin{pmatrix}-1\\1\\1\\-1 \end{pmatrix}\cdot \begin{pmatrix}-1\\ -4\\ 7\\ 14\end{pmatrix} }\right] \begin{pmatrix}0\\1\\1\\0\end{pmatrix}\\=& \underbrace{ 2\begin{pmatrix}1\\2\\3\\4\end{pmatrix}- 3\begin{pmatrix}2\\2\\-1\\-1\end{pmatrix}+ 3\begin{pmatrix}1\\1\\1\\1\end{pmatrix}- 5\begin{pmatrix}0\\1\\1\\0\end{pmatrix}}_\text{成分分解}= \begin{pmatrix}-1\\-4\\7\\14\end{pmatrix}\\ \end{aligned} $$

が成り立ち、正しく成分分解されることを確認できる。

逆基底の求め方(外積を使う)

$\{\bm b_i\}$ を順に横に並べ、$j$ 番目だけ $\bm x$ に置き換えた行列式を作り、これを $K^n\to K$ の写像 $T_j$ とみなす。

$$T_j(\bm x)= \Bigg| \,\bm b_1\ \bm b_2\ \dots\ \underset{\overset{\wedge}{j}}{\bm x}\ \dots\ \bm b_n \Bigg| $$

この $T_j$ は線形なので、 その表現行列である $1\times n$ 行列($n$ 次元横ベクトル)を $\bm t_j^\dagger$ と書けば

$$T_j(\bm x)=\bm t_j^\dagger\bm x=(\bm t_j,\bm x)$$

のように標準内積で表せる。(行列式を展開すればこの形になるということ)

実はこの $\bm t_j$ から、

$$\bm \beta_j=\frac{1}{\overline{T_j}(\bm b_j)}\bm t_j$$

とすることで $\bm \beta_j$ を求められる。

  • $i\ne j$ のとき $(t_j,\bm b_i)=T_j(\bm b_i)=0$ ← 行列式の内に同じ行ベクトルが2回現れるから
  • $i=j$ のときは $(t_j,\bm b_i)=T_j(\bm b_i)=T_j(\bm b_j)$

に注意せよ($\beta_j$ は内積の左側に入れるので複素共役の $1/\overline T_j(\bm b_j)$ になっている)。

どこに外積が出てきた?

$n-1$ 個の $n$ 次元ベクトル $\bm a_1,\bm a_2,\dots,\bm a_{n-1}$ のベクトル積(外積) $\bm a_1\times\bm a_2\times\dots\times\bm a_{n-1}$ は、

$$\Bigg|\,\bm x\ \bm a_1\ \bm a_2\ \dots\ \bm a_{n-1}\Bigg|\,=\,(\,\bm a_1\times\bm a_2\times\dots\times\bm a_{n-1},\,\bm x\,)$$

を満たすベクトルとして定義される

この定義を用いると、

$$\bm t_j=(-1)^{j-1}\bm b_1\times\bm b_2\times\dots(j番目欠)\dots\times\bm b_n$$

$$ \bm \beta_j=\frac{(-1)^{j-1}\bm b_1\times\bm b_2\times\dots(j番目欠)\dots\times\bm b_n}{\overline{{}^{\ }\Bigg|\,\bm b_1\ \bm b_2\ \dots\ \bm b_n\Bigg|}^{\ }} $$

と表せる。

  • $n$ 次元のベクトル積は $n-1$ 個のベクトルの積のみが定義される。$n-1$ 個以外は定義されない。
  • 符号の $(-1)^{j-1}$ は、行列式内で $\bm x$ の列を一番左に持ってくるのに必要な互換の回数に対応する。

上の例の $\bm b_1,\bm b_2,\bm b_3,\bm b_4$ に対して、

$$T_1(\bm b_1)=T_2(\bm b_2)=T_3(\bm b_3)=T_4(\bm b_4)= \begin{vmatrix} 1&2&1&0\\ 2&2&1&1\\ 3&-1&1&1\\ 4&-1&1&0 \end{vmatrix}=6\hspace{1cm}\text{← 分母は共通になる} $$

$\displaystyle\frac{T_1(\bm x)}{T_1(\bm b_1)}=\frac{1}{6}\begin{vmatrix} x&2&1&0\\ y&2&1&1\\ z&-1&1&1\\ w&-1&1&0 \end{vmatrix}=\frac{1}{2}(-x+y-z+w)=\frac{1}{2}\begin{pmatrix} 1&-1&1&-1 \end{pmatrix} \begin{pmatrix}x\\y\\z\\w\end{pmatrix} $ より $\bm \beta_1=\frac{1}{2}\begin{pmatrix}-1\\1\\-1\\1 \end{pmatrix}$

$\displaystyle\frac{T_2(\bm x)}{T_2(\bm b_2)}=\frac{1}{6}\begin{vmatrix} 1&x&1&0\\ 2&y&1&1\\ 3&z&1&1\\ 4&w&1&0 \end{vmatrix}=\frac{1}{6}(-x+3y-3z+w)=\frac{1}{6}\begin{pmatrix}-1&3&-3&1 \end{pmatrix} \begin{pmatrix}x\\y\\z\\w\end{pmatrix} $ より $\bm \beta_2=\frac{1}{6}\begin{pmatrix}-1\\3\\-3\\1 \end{pmatrix}$

$\displaystyle\frac{T_3(\bm x)}{T_3(\bm b_3)}=\frac{1}{6}\begin{vmatrix} 1&2&x&0\\ 2&2&y&1\\ 3&-1&z&1\\ 4&-1&w&0 \end{vmatrix}=\frac{1}{6}(11x-9y+9z-5w)=\frac{1}{6}\begin{pmatrix} 11&-9&9&-5 \end{pmatrix} \begin{pmatrix}x\\y\\z\\w\end{pmatrix} $ より $\bm \beta_3=\frac{1}{6}\begin{pmatrix} 11\\-9\\9\\-5 \end{pmatrix}$

$\displaystyle\frac{T_4(\bm x)}{T_4(\bm b_4)}=\frac{1}{6}\begin{vmatrix} 1&2&1&x\\ 2&2&1&y\\ 3&-1&1&z\\ 4&-1&1&w \end{vmatrix}=\frac{1}{2}(-x+y+z-w)=\frac{1}{2}\begin{pmatrix}-1&1&1&-1 \end{pmatrix} \begin{pmatrix}x\\y\\z\\w\end{pmatrix} $ より $\bm \beta_4=\frac{1}{2}\begin{pmatrix}-1\\1\\1\\-1 \end{pmatrix}$

このように、上記と同じ結果が得られる。

ここでは省いたが実は行列式を求めるところは手計算でまじめにやるとかなり大変なので、 いっぺんにすべてを求めるのであれば逆行列を使う方がずっといい。

直交しない空間への成分分解

全体空間が $K^n=V\dot +W$ のように必ずしも直交しない2つの部分空間の直和で書けるとき、

$V$ の基底を $\bm b_1,\bm b_2,\dots,\bm b_m$

$W$ の基底を $\bm b_{m+1},\bm b_{m+2},\dots,\bm b_n$

とすれば、両者を合わせた $\bm b_1,\bm b_2,\dots,\bm b_n$ は $K^n$ の基底となる。

対応する逆基底を $\bm \beta_1,\bm \beta_2,\dots,\bm \beta_n$ とすると、 任意の $\bm x\in K^n$ は

$$ \begin{aligned} \bm x &=\sum_{i=1}^n (\bm \beta_i,\bm x)\bm b_i\\ &=\underbrace{\sum_{i=1}^m (\bm \beta_i,\bm x)\bm b_i}_{\in\,V}\ \,+ \underbrace{\sum_{i=m+1}^n (\bm \beta_i,\bm x)\bm b_i}_{\in\,W}\\ &=\bm x_V+\bm x_W \end{aligned} $$

のように $V$ 成分と $W$ 成分とに一意に分解できる。

さらに、

$$ \begin{aligned} \bm x_V &=\sum_{i=1}^m (\bm \beta_i,\bm x)\bm b_i\\ &=\sum_{i=1}^m \bm b_i(\bm \beta_i,\bm x)\\ &=\underbrace{\sum_{i=1}^m \bm b_i{\bm \beta_i}^\dagger}_{P_V}\ \bm x\\ &=P_V\bm x \end{aligned} $$

として $V$ への射影演算子が得られる。

$$P_V=\sum_{i=1}^m \bm b_i{\bm \beta_i}^\dagger$$

正規直交基底の時には $\bm \beta_i=\bm b_i$ と取れるため、 以前に見た通り

$$P_V=\sum_{i=1}^m \bm b_i\bm b_i^\dagger$$

に戻って、つじつまが合ってる。

この式を、「逆基底 $\bm\beta_i$ を掛けて $\bm b_i$ の係数を求めてから $\bm b_i$ を掛ければ成分そのものが得られるから、必要な成分を全部集めれば射影が求まる」と読めるようになっていればOKだ。

逆格子の基本ベクトル

固体物理において3次元空間の「逆格子の基本ベクトル」を求めるための公式として、

$${\bm \beta_1}'=2\pi\cdot\frac{\bm b_2\times\bm b_3}{\bm b_1\cdot(\bm b_2\times\bm b_3)}\hspace{1cm}=\ \ 2\pi\frac{(\bm b_1\,\text{を除いた外積})}{\ (\bm b_1\,\text{を除いた外積と}\,\bm b_1\,\text{の内積)}\ }$$

のようなものを学ぶが、これは係数 $2\pi$ を除けば上記の逆基底そのものである。

波数ベクトル $\bm k$ が基本逆格子ベクトルの整数倍の和であり、

$$\bm k=l{\bm \beta_1}'+m{\bm \beta_1}'+n{\bm \beta_1}'$$

実ベクトル $\bm x$ が基本格子ベクトルの整数倍の和であるとき、

$$\bm x=l'\bm b_1+m'\bm b_2+n'\bm b_3$$

両者の積は、

$$ \begin{aligned} \bm k\cdot\bm x &=(l{\bm \beta_1}'+m{\bm \beta_2}'+n{\bm \beta_3}')\cdot(l'\bm b_1+m'\bm b_2+n'\bm b_3)\\ &=2\pi(l{\bm \beta_1}+m{\bm \beta_2}+n{\bm \beta_3})\cdot(l'\bm b_1+m'\bm b_2+n'\bm b_3)\\ &=2\pi(ll'\bm \beta_1\cdot\bm b_1+mm'\bm \beta_2\cdot\bm b_2+nn'\bm \beta_3\cdot\bm b_3\underbrace{+\ 0}_\text{他は消える})\\[-10px] &=2\pi(ll'+mm'+nn')\\ \end{aligned}$$

となって、

$$e^{i\bm k\cdot\bm x}=1$$

を得る。

すなわち、波数空間に於いて逆格子ベクトルの格子点に来る波数の平面波は 実空間の格子点上で位相がぴったり揃うことになり、 これが固体物理のいろんなところで役に立つことになる。

本来の逆基底

上では「直交するベクトル」という感覚や逆格子ベクトルとの整合性を重視して、 縦ベクトルである逆基底(仮) $\bm \beta_i$ を考えたが、 数学的には逆基底はそのエルミート共役、すなわち横ベクトルとして定義される。

すなわち、$\bm b_i$ に対応する逆基底 $\bm b^i$ は、

$$\bm b^i=\bm \beta_i^\dagger$$

このとき、$(\bm \beta_i,\bm x)={\bm \beta_i}^\dagger\bm x=\bm b^i\bm x$ であるから、 逆基底の定義は

$$\bm b^i\,\bm b_j=\delta_{ij}$$

逆行列や外積を用いた定義は、

$$\begin{aligned} \bm b^i &=\underbrace{(0\ 0\ \dots\ 1\ \dots 0)}_{i\,\text{番目だけが}\,1}\Bigg(\bm b_1\ \bm b_2\ \dots\ \bm b_n\Bigg)^{-1}\\[6mm] &=\frac{(-1)^{i-1} }{\bigg|\bm b_1\ \bm b_2\ \dots\ \bm b_n\bigg|}\,\bm b_1\times\bm b_2\times\dots(i番目欠)\dots\times\bm b_n \end{aligned}$$

となる。このように定義された逆基底は「双対基底」とも呼ばれる。

双対空間について

ある $K$ 上のベクトル空間 $V$ に対して、 「$V$ の元を $K$ の元へ移す線形写像をすべて集めた集合 $U$」

$$U=\{\,f:V\to K\,|\,\forall x,y\in V,\forall a,b\in K,f(ax+by)=af(x)+bf(y)\,\}$$

を $V$ の双対空間と呼び $U=V^*$ と書く。

この集合内での和とスカラー倍を $f,g\in V^*, x\in V, a\in K$ に対して、

$$ f+g:x\mapsto f(x)+g(x) $$ $$ af:x\mapsto af(x) $$

のように定義すると $V^*$ はベクトル空間となる。

「線形写像が作るベクトル空間」と聞くと難しそうだが、$V$ を $n$ 次元空間とすれば $U$ の元 $f$ は $n$ 次元から $1$ 次元への線形写像であるから、適当な基底の下で $1$ 行 $n$ 列の「横ベクトル空間」と同一視できる。

上で見た「双対基底」はこの「双対空間」の基底になっており、元の空間の基底を $\bm b_i$ と書くとき、 対応する双対基底は上付き数字で $\bm b^i=\bm\beta_i^\dagger$ と書くのが一般的である。

質問・コメント





添付ファイル: fileinverse_basis.png 354件 [詳細]

Counter: 9784 (from 2010/06/03), today: 3, yesterday: 6