Forest background
バイオインフォの森へようこそ
ベクトルの導入と演算

ベクトルの導入と演算

目次

1. ベクトルの定義

(1) ベクトルの定義と表記

定義1.1: ベクトル

a1,a2,,ana_1, a_2, \cdots, a_n1,2,3+i1, 2, 3+iなどの数とすると、ベクトルとは以下のように数を縦に複数を並べたもので、記号としてはa\mathbf{a}のようにボールド体を用いて表します。

a=(a1a2an)\begin{align} \mathbf{a} = \begin{pmatrix} a_1 \\ a_2 \\ \vdots \\ a_n \end{pmatrix} \end{align}

なお、以下のように数を横に並べてベクトルとする場合もありますが、特に縦横の指定がなくベクトルが定義された場合は縦に数字を並べたベクトルと解釈することに注意して下さい。

a=(a1,a2,,an)\begin{align} \mathbf{a'} = (a_1, a_2, \cdots, a_n) \end{align}

a\mathbf{a}a\mathbf{a'}は数学的に異なり、両者を区別するためにa\mathbf{a}縦ベクトルa\mathbf{a'}横ベクトルと呼びます。後述する行列の考え方を用いると、a\mathbf{a}nn11列の行列で、a\mathbf{a'}11nn列の行列としてはっきり区別されます。

※ より正確に述べますと、単純に数を並べただけではベクトルとは呼びません。後述する和(差)とスカラー倍が決められている場合にベクトルと呼びます。

(2) ベクトルの次元と成分

定義1.2: 次元と成分

並べられた数がnn個のベクトルは nn次元のベクトルと呼びます。また、並べられたそれぞれの数はベクトルの成分と言います。

ベクトルの次元と成分を簡潔に表記する方法があります。まず、実数全体の集合をR\mathbb{R}、複素数全体の集合をC\mathbb{C}と表記します。また、成分が全て実数のnn次元ベクトル全体の集合をRn\mathbb{R}^n、成分に複素数を含むnn次元ベクトル全体の集合をCn\mathbb{C}^nと表記します。このとき、ベクトルa\mathbf{a}が実数成分のnn次元ベクトルであること、すなわちnn次元の実ベクトルであることをaRn\mathbf{a} \in \mathbb{R}^n、成分に複素数を含むnn次元のベクトル、すなわちnn次元の複素ベクトルであることをaCn\mathbf{a} \in \mathbb{C}^nと表記します。

(3) 様々なベクトル(1):ゼロベクトル

定義1.3: ゼロベクトル

すべての成分が0であるベクトルをゼロベクトルと呼び、0\mathbf{0}(00)\begin{pmatrix} 0 \\ \vdots \\ 0 \end{pmatrix} と表記します。

(4) 様々なベクトル(2):単位ベクトル

定義1.4: 単位ベクトル

単位ベクトルとは、大きさ(ノルム)が1のベクトルです。

※ ノルムの詳細については1章4節を参照して下さい。

(5) 様々なベクトル(3):幾何ベクトル

定義1.5: 幾何ベクトル

2次元または3次元のベクトルは原点からその成分が表す座標への矢印とみなすことが可能であり、その意味でベクトルは大きさと向きを持つ量とも言えます。

このようにベクトルの幾何学的側面を強調して記述する場合、そのベクトルを幾何ベクトルと呼び、OA\vec{OA}のように表します。

(6) 様々なベクトル(4):位置ベクトル

定義1.6: 位置ベクトル

原点を始点として平面や空間上の点 P(x,y,z)P(x, y, z) などへのベクトルを特に位置ベクトルと呼び、以下のように表現されます。

r=(xyz)\begin{align} \mathbf{r} = \begin{pmatrix} x \\ y \\ z \end{pmatrix} \end{align}

(7) 様々なベクトル(5):方向ベクトル

定義1.7: 方向ベクトル

直線や線分と平行なベクトルを方向ベクトルといいます。

2点 A(x1,y1,z1)A(x_1, y_1, z_1), B(x2,y2,z2)B(x_2, y_2, z_2) を通る直線の方向ベクトルは以下です。

AB=(x2x1y2y1z2z1)\begin{align} \vec{AB} = \begin{pmatrix} x_2 - x_1 \\ y_2 - y_1 \\ z_2 - z_1 \end{pmatrix} \end{align}

また、以下の直線の方程式においては方向ベクトルが(abc)\begin{pmatrix} a \\ b \\ c \end{pmatrix}で与えられます。

xx0a=yy0b=zz0c\begin{align} \frac{x-x_0}{a}=\frac{y-y_0}{b}=\frac{z-z_0}{c} \end{align}

※ 空間の直線の方程式の詳細については[2章][2]を参照して下さい。

(8) 様々なベクトル(6):法線ベクトル

定義1.8: 法線ベクトル

ある直線や平面に対して垂直なベクトルを法線ベクトルといいます。

例えば、直線 ax+by+c=0ax + by + c = 0 に対する法線ベクトルは

n=(ab)\begin{align} \mathbf{n} = \begin{pmatrix} a \\ b \end{pmatrix} \end{align}

です。

また、平面 ax+by+cz+d=0ax + by + cz + d = 0 に対する法線ベクトルは

n=(abc)\begin{align} \mathbf{n} = \begin{pmatrix} a \\ b \\ c \end{pmatrix} \end{align}

です。

※ 平面の方程式の詳細については[2章][2]を参照して下さい。

2. ベクトルの和と差

(1) ベクトル和の定義

定義1.9: ベクトルの和

aCn\mathbf{a} \in \mathbb{C}^nbCm\mathbf{b} \in \mathbb{C}^mの和は次元が等しい場合、つまりn=mn=m場合に定義されa+b\mathbf{a}+\mathbf{b}と表記し、その成分は対応する成分の和とします。すなわち、a=(a1an)\mathbf{a}=\begin{pmatrix}a_1\\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b}=\begin{pmatrix}b_1 \\ \vdots \\ b_n \end{pmatrix}とすると、a\mathbf{a}b\mathbf{b}の和は以下です。

a+b=(a1+b1an+bn)\begin{align} \mathbf{a} + \mathbf{b} = \begin{pmatrix} a_1 + b_1 \\ \vdots \\ a_n + b_n \end{pmatrix} \end{align}

(2) ベクトルの差の定義

定義1.10: ベクトルの差

aCn\mathbf{a} \in \mathbb{C}^nbCm\mathbf{b} \in \mathbb{C}^mの差は次元が等しい場合、つまりn=mn=m場合に定義されab\mathbf{a}-\mathbf{b}と表記し、その成分は対応する成分の差とします。すなわち、a=(a1an)\mathbf{a}=\begin{pmatrix}a_1\\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b}=\begin{pmatrix}b_1 \\ \vdots \\ b_n \end{pmatrix}とすると、a\mathbf{a}b\mathbf{b}の差は以下です。

ab=(a1b1anbn)\begin{align} \mathbf{a} - \mathbf{b} = \begin{pmatrix} a_1 - b_1 \\ \vdots \\ a_n - b_n \end{pmatrix} \end{align}

(3) 和(差)の性質

交換法則

a,bCn\mathbf{a}, \mathbf{b} \in \mathbb{C}^nについて、以下が成立します。

a+b=b+a\begin{align} \mathbf{a} + \mathbf{b} = \mathbf{b} + \mathbf{a} \end{align}
証明をみる

a=(a1an)\mathbf{a}=\begin{pmatrix}a_1\\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b}=\begin{pmatrix}b_1 \\ \vdots \\ b_n \end{pmatrix}とすると、

a+b=(a1+b1an+bn)=(b1+a1bn+an)=b+a\begin{align} \mathbf{a} + \mathbf{b} &= \begin{pmatrix}a_1+b_1\\ \vdots \\ a_n+b_n\end{pmatrix} \\ &= \begin{pmatrix}b_1+a_1\\ \vdots \\ b_n+a_n\end{pmatrix} \\ &= \mathbf{b} + \mathbf{a} \end{align}

より、成立します。

結合法則

a,b,cCn\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{C}^nについて、以下が成立します。

(a+b)+c=a+(b+c)\begin{align} (\mathbf{a} + \mathbf{b}) + \mathbf{c} = \mathbf{a} + (\mathbf{b} + \mathbf{c}) \end{align}
証明をみる

a=(a1an)\mathbf{a}=\begin{pmatrix}a_1\\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b}=\begin{pmatrix}b_1 \\ \vdots \\ b_n \end{pmatrix}, c=(c1cn)\mathbf{c}=\begin{pmatrix}c_1 \\ \vdots \\ c_n \end{pmatrix}とすると

(a+b)+c=((a1+b1)+c1(an+bn)+cn)=(a1+(b1+c1)an+(bn+cn))=a+(b+c)\begin{align} (\mathbf{a} + \mathbf{b}) + \mathbf{c} &= \begin{pmatrix}(a_1+b_1)+c_1\\ \vdots \\ (a_n+b_n)+c_n\end{pmatrix} \\ &= \begin{pmatrix}a_1+(b_1+c_1)\\ \vdots \\ a_n+(b_n+c_n)\end{pmatrix} \\ &= \mathbf{a} + (\mathbf{b} + \mathbf{c}) \end{align}

より、成立します。

3. ベクトルのスカラー倍

(1) スカラーの定義

定義1.11: スカラー

スカラーとは、1,3+i1, 3+iのようなベクトルではない1次元の普通の数を指す言葉で、ベクトルでないことを強調する際に用います。α\alphaをスカラーとした場合、α\alphaが実数であればαR\alpha \in \mathbb{R}α\alphaが複素数であればαC\alpha \in \mathbb{C}と表記します。

(2) スカラー倍の定義

定義1.12: ベクトルのスカラー倍

スカラー λC\lambda \in \mathbb{C} によるベクトル a=(a1an)Cn\mathbf{a}=\begin{pmatrix}a_1 \\ \vdots \\ a_n \end{pmatrix} \in \mathbb{C}^n のスカラー倍は以下のように定義されます。

λa=(λa1λan)\begin{align} \lambda \mathbf{a} = \begin{pmatrix} \lambda a_1 \\ \vdots \\ \lambda a_n \end{pmatrix} \end{align}

(3) スカラー倍の性質

結合法則(互換性)

α,βC\alpha, \beta \in \mathbb{C}aCn\mathbf{a} \in \mathbb{C}^n について、以下が成立します。

α(βa)=(αβ)a\begin{align} \alpha (\beta \mathbf{a}) = (\alpha \beta) \mathbf{a} \end{align}
証明をみる

a=(a1an)\mathbf{a} = \begin{pmatrix}a_1 \\ \vdots \\ a_n\end{pmatrix}とすると、

α(βa)=α(βa1βan)=(αβa1αβan)=((αβ)a1(αβ)an)=(αβ)a\begin{align} \alpha (\beta \mathbf{a}) &= \alpha \begin{pmatrix} \beta a_1 \\ \vdots \\ \beta a_n \end{pmatrix}\\ &= \begin{pmatrix} \alpha\beta a_1 \\ \vdots \\ \alpha\beta a_n \end{pmatrix}\\ &= \begin{pmatrix} (\alpha\beta) a_1 \\ \vdots \\ (\alpha\beta) a_n \end{pmatrix}\\ &= (\alpha \beta) \mathbf{a} \end{align}

より成立します。

分配法則

αC\alpha \in \mathbb{C}a,bCn\mathbf{a}, \mathbf{b} \in \mathbb{C}^nについて、以下が成立します。

α(a+b)=αa+αb\begin{align} \alpha(\mathbf{a} + \mathbf{b}) = \alpha \mathbf{a} + \alpha \mathbf{b} \end{align}
証明をみる

a=(a1an)\mathbf{a} = \begin{pmatrix}a_1 \\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b} = \begin{pmatrix}b_1 \\ \vdots \\ b_n\end{pmatrix}とすると、

α(a+b)=α(a1+b1an+bn)=(α(a1+b1)α(an+bn))=(αa1+αb1αan+αbn)=((αa1)+(αb1)(αan)+(αbn))=αa+αb\begin{align} \alpha (\mathbf{a} + \mathbf{b}) &= \alpha \begin{pmatrix} a_1+b_1 \\ \vdots \\ a_n+b_n \end{pmatrix}\\ &= \begin{pmatrix} \alpha(a_1+b_1) \\ \vdots \\ \alpha(a_n+b_n) \end{pmatrix}\\ &= \begin{pmatrix} \alpha a_1 + \alpha b_1 \\ \vdots \\ \alpha a_n + \alpha b_n \end{pmatrix}\\ &= \begin{pmatrix} (\alpha a_1) + (\alpha b_1) \\ \vdots \\ (\alpha a_n) + (\alpha b_n) \end{pmatrix}\\ &= \alpha \mathbf{a} + \alpha \mathbf{b} \end{align}

より成立します。

4. ベクトルの内積とノルム

本節では実ベクトルの内積を扱います。複素ベクトルに拡張した内積については[13章][13]を参照して下さい。また、本節以降では縦ベクトルa=(a1an)Rn\mathbf{a}=\begin{pmatrix}a_1 \\ \vdots \\ a_n\end{pmatrix} \in \mathbb{R}^nの成分を横に並べた横ベクトルa=(a1,,an)\mathbf{a'} = (a_1, \cdots, a_n)を、簡単のためa=aT\mathbf{a'}=\mathbf{a}^{\mathsf T}と表記することにします。aT\mathbf{a}^{\mathsf T}a\mathbf{a}の転置という概念であり、詳細については[3章][3]を参照して下さい。

(1) 内積の定義

定義1.13: 内積

aRn\mathbf{a} \in \mathbb{R}^nbRm\mathbf{b} \in \mathbb{R}^mの内積は同じ次元、すなわちn=mn=mの場合に定義され、ab\mathbf{a}\cdot\mathbf{b}と表記し、以下で定義されます。

ab=aTb\begin{align} \mathbf{a}\cdot\mathbf{b} = \mathbf{a}^{\mathsf T}\mathbf{b} \end{align}

aTb\mathbf{a}^{\mathsf T}\mathbf{b}は、a=(a1an)\mathbf{a}=\begin{pmatrix}a_1 \\ \vdots \\ a_n\end{pmatrix}, b=(b1bn)\mathbf{b}=\begin{pmatrix}b_1 \\ \vdots \\ b_n\end{pmatrix}とすると、以下で定義されます。

aTb=i=1naibi\begin{align} \mathbf{a}^{\mathsf T}\mathbf{b} = \sum_{i=1}^n a_ib_i \end{align}

内積ab=aTb\mathbf{a}\cdot\mathbf{b} = \mathbf{a}^{\mathsf T}\mathbf{b}は2つのベクトルからスカラーを与える演算であることに注意して下さい。また、aa\mathbf{a}\cdot\mathbf{a}a2\mathbf{a}^2と表記することがあることにも注意して下さい。

(2) 内積の性質

交換法則

a,bRn\mathbf{a}, \mathbf{b} \in \mathbb{R}^nについて、以下が成立します。

ab=ba\begin{align} \mathbf{a} \cdot \mathbf{b} = \mathbf{b} \cdot \mathbf{a} \end{align}
証明をみる

a,b\mathbf{a}, \mathbf{b}の第ii成分をそれぞれai,bia_i, b_iとすると、

ab=aTb=iaibi=ibiai=bTa=ba\begin{align} \mathbf{a} \cdot \mathbf{b} &= \mathbf{a}^{\mathsf T} \mathbf{b} \\ &= \sum_i a_i b_i \\ &= \sum_i b_i a_i \\ &= \mathbf{b}^{\mathsf T} \mathbf{a} \\ &= \mathbf{b} \cdot \mathbf{a} \end{align}
線形性

a,b,cRn\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{R}^nについて、以下が成立します。

a(b+c)=ab+ac\begin{align} \mathbf{a} \cdot (\mathbf{b} + \mathbf{c}) = \mathbf{a} \cdot \mathbf{b} + \mathbf{a} \cdot \mathbf{c} \end{align}
証明をみる

a,b,c\mathbf{a}, \mathbf{b}, \mathbf{c}の第ii成分をそれぞれai,bi,cia_i, b_i, c_iとすると、

a(b+c)=aT(b+c)=iai(bi+ci)=iaibi+i=1naici=aTb+aTc=ab+ac\begin{align} \mathbf{a} \cdot (\mathbf{b} + \mathbf{c}) &= \mathbf{a}^{\mathsf T} (\mathbf{b} + \mathbf{c}) \\ &= \sum_i a_i (b_i + c_i) \\ &= \sum_i a_i b_i + \sum_{i=1}^n a_i c_i \\ &= \mathbf{a}^{\mathsf T} \mathbf{b} + \mathbf{a}^{\mathsf T} \mathbf{c} \\ &= \mathbf{a} \cdot \mathbf{b} + \mathbf{a} \cdot \mathbf{c} \end{align}
正値性

aRn\mathbf{a} \in \mathbb{R}^nについて、以下が成立します。

aa0\begin{align} \mathbf{a} \cdot \mathbf{a} \geq 0 \end{align}
証明をみる

a\mathbf{a}の第ii成分をaia_iとすると、

aa=aTa=iaiai=iai20\begin{align} \mathbf{a} \cdot \mathbf{a} &= \mathbf{a}^{\mathsf T} \mathbf{a} \\ &= \sum_i a_i a_i\\ &= \sum_i a_i^2 \\ &\geq 0 \end{align}

(3) ノルムの定義

定義1.14: ノルム

a=(a1an)Rn\mathbf{a} = \begin{pmatrix}a_1 \\ \vdots \\ a_n \end{pmatrix} \in \mathbb{R}^nノルムa\|\mathbf{a}\|と表記し、自身との内積の平方根で以下のように定義されます。なお、ノルムは長さ大きさとも呼ばれます。

a=aa=i=1nai2\begin{align} \|\mathbf{a}\| &= \sqrt{\mathbf{a} \cdot \mathbf{a}} \\ &= \sqrt{\sum_{i=1}^n a_i^2} \end{align}

(4) ノルムの性質

以下の性質があります。

ノルムの性質

λR,aRn\lambda \in \mathbb{R}, \mathbf{a} \in \mathbb{R}^nについて、以下が成立します。

  • a0\|\mathbf{a}\| \geq 0
  • a=0a=0\|\mathbf{a}\| = 0 \Leftrightarrow \mathbf{a} = \mathbf{0}
  • λa=λa\|\lambda \mathbf{a}\| = |\lambda| \|\mathbf{a}\|

証明は割愛します。

(5) 内積とノルムに関する定理(1):単位ベクトル

定理1.15: 内積とノルムの関係

ゼロベクトルでないaRn\mathbf{a} \in \mathbb{R}^nについて、a/a\mathbf{a}/\|\mathbf{a}\|は単位ベクトルとなります。

証明をみる

a/a\mathbf{a}/\|\mathbf{a}\|のノルムが11であることを示せば十分です。a=n\|\mathbf{a}\|=nとおくと、

aa=an=an=nn=1\begin{align} \Big\|\frac{\mathbf{a}}{\|\mathbf{a}\|}\Big\| &= \Big\|\frac{\mathbf{a}}{n}\Big\| \\ &= \frac{\|\mathbf{a}\|}{n} \\ &= \frac{n}{n} \\ &= 1 \end{align}

よって、題意成立となります。

(6) 内積とノルムに関する定理(2):内積とノルムの関係

定理1.16: 内積とノルムの関係

a\mathbf{a}b\mathbf{b}が2次元または3次元の実ベクトルであるとき、その2つのベクトルの成す角をθ\thetaとすると、内積は以下にように表現できます。

ab=abcosθ\begin{align} \mathbf{a} \cdot \mathbf{b} = \|\mathbf{a}\|\|\mathbf{b}\|\cos\theta \end{align}
証明をみる

2次元平面または3次元空間上において、原点 OO、点 AA(位置ベクトル a\mathbf{a})、点 BB(位置ベクトル b\mathbf{b})をとります。三角形 OABOAB について、辺の長さはOA=aOA=\|\mathbf{a}\|OB=bOB=\|\mathbf{b}\|AB=abAB=\|\mathbf{a}-\mathbf{b}\|となります。また、AOB=θ\angle AOB=\thetaa\mathbf{a}b\mathbf{b} のなす角です。余弦定理より

ab2=a2+b22abcosθ\begin{align} \|\mathbf{a}-\mathbf{b}\|^2=\|\mathbf{a}\|^2+\|\mathbf{b}\|^2-2\|\mathbf{a}\|\|\mathbf{b}\|\cos\theta \end{align}

となります。一方、内積の性質から

ab2=(ab)(ab)=aa+bb2ab=a2+b22ab\begin{align} \|\mathbf{a}-\mathbf{b}\|^2 &= (\mathbf{a}-\mathbf{b})\cdot(\mathbf{a}-\mathbf{b}) \\ &= \mathbf{a}\cdot\mathbf{a}+\mathbf{b}\cdot\mathbf{b}-2\mathbf{a}\cdot\mathbf{b} \\ &= \|\mathbf{a}\|^2+\|\mathbf{b}\|^2-2\mathbf{a}\cdot\mathbf{b} \end{align}

となります。従って、

a2+b22ab=a2+b22abcosθ\begin{align} \|\mathbf{a}\|^2+\|\mathbf{b}\|^2-2\mathbf{a}\cdot\mathbf{b}=\|\mathbf{a}\|^2+\|\mathbf{b}\|^2-2\|\mathbf{a}\|\|\mathbf{b}\|\cos\theta \end{align}

両辺から a2+b2\|\mathbf{a}\|^2+\|\mathbf{b}\|^2 を消去して 2-2 で割ると

ab=abcosθ\begin{align} \mathbf{a}\cdot\mathbf{b}=\|\mathbf{a}\|\|\mathbf{b}\|\cos\theta \end{align}

となります。以上より、題意成立となります。

(7) 内積とノルムに関する定理(3):2つのベクトルが成す角度

定理1.17: 2つのベクトルが成す角度

a\mathbf{a}b\mathbf{b}が2次元または3次元の実ベクトルであるとき、その2つのベクトルの成す角をθ\thetaを以下のように求めることができます。

cosθ=abab\begin{align} \cos\theta = \frac{\mathbf{a} \cdot \mathbf{b}}{\|\mathbf{a}\| \|\mathbf{b}\|} \end{align}
証明をみる

内積とノルムの関係式を式変形することで得られます。

(8) 内積とノルムに関する定理(4):直交の判定

定理1.18: 直交の判定

a\mathbf{a}b\mathbf{b}が2次元または3次元の実ベクトルであるとき、「a\mathbf{a}b\mathbf{b}が直交」と「ab\mathbf{a}\cdot\mathbf{b}」は同値です。

証明をみる

内積とノルムの関係式でθ=π/2\theta=\pi/2を代入することで証明できます。

(9) 内積とノルムに関する定理(5):コーシー・シュワルツの不等式

定理1.19: 三角不等式

a,bRn\mathbf{a}, \mathbf{b} \in \mathbb{R}^nについて、以下が成立します。

abab\begin{align} \mathbf{a} \cdot \mathbf{b} \leq \|\mathbf{a}\|\|\mathbf{b}\| \end{align}
証明をみる

tRt\in\mathbb{R}について

atb2=b2t22abt+a20\begin{align} \|\mathbf{a}-t\mathbf{b}\|^2 &= \|\mathbf{b}\|^2t^2 -2\mathbf{a}\cdot\mathbf{b}t + \|\mathbf{a}\|^2 \\ & \geq 0 \end{align}

となります。これは、ttの2次式とみなした際の判別式が非正であることと同値なので、

(ab)2a2b20    (ab)2a2b2    (ab)ab\begin{align} &(\mathbf{a}\cdot\mathbf{b})^2 - \|\mathbf{a}\|^2\|\mathbf{b}\|^2 \leq 0 \\ &\iff (\mathbf{a}\cdot\mathbf{b})^2 \leq \|\mathbf{a}\|^2\|\mathbf{b}\|^2 \\ &\iff (\mathbf{a}\cdot\mathbf{b}) \leq \|\mathbf{a}\|\|\mathbf{b}\| \end{align}

となります。以上より、題意成立となります。

※ 以下の関係式においてcosθ\cos \theta1-1から11であることを用いる証明では、2次元または3次元ベクトルに対する証明にしかならないことに注意して下さい。

ab=abcosθ\begin{align} \mathbf{a} \cdot \mathbf{b} = \|\mathbf{a}\|\|\mathbf{b}\|\cos\theta \end{align}

(10) 内積とノルムに関する定理(6):三角不等式

定理1.20: 2つのベクトルが成す角度

a,bRn\mathbf{a}, \mathbf{b} \in \mathbb{R}^nについて、以下が成立します。

a+ba+b\begin{align} \|\mathbf{a} + \mathbf{b}\| \leq \|\mathbf{a}\| + \|\mathbf{b}\| \end{align}
証明をみる
(a+b)2=a2+b2+2ab\begin{align} (\|\mathbf{a}\| + \|\mathbf{b}\|)^2 = \|\mathbf{a}\|^2 + \|\mathbf{b}\|^2 + 2\|\mathbf{a}\|\|\mathbf{b}\| \end{align}

となります。ここで、コーシー・シュワルツの不等式を用いると、

(a+b)2=a2+b2+2aba2+b2+2ab=(a+b)2\begin{align} (\|\mathbf{a}\| + \|\mathbf{b}\|)^2 & = \|\mathbf{a}\|^2 + \|\mathbf{b}\|^2 + 2\|\mathbf{a}\|\|\mathbf{b}\| \\ &\geq \|\mathbf{a}\|^2 + \|\mathbf{b}\|^2 + 2\mathbf{a}\cdot\mathbf{b} \\ &= (\mathbf{a} + \mathbf{b})^2 \end{align}

となります。従って、

(a+b)2(a+b)2    a+ba+b    a+ba+b\begin{align} &(\|\mathbf{a}\| + \|\mathbf{b}\|)^2 \geq (\mathbf{a} + \mathbf{b})^2 \\ &\iff |\|\mathbf{a}\| + \|\mathbf{b}\|| \geq \|\mathbf{a} + \mathbf{b}\| \\ &\iff \|\mathbf{a}\| + \|\mathbf{b}\| \geq \|\mathbf{a} + \mathbf{b}\| \end{align}

となり、題意成立となります。

(11) 内積の公式 (1)

定理1.21: 内積の公式(1)

a,bRn\mathbf{a}, \mathbf{b} \in \mathbb{R}^nについて、以下が成立します。

a(aTb)=(aaT)b\begin{align} \mathbf{a} (\mathbf{a}^{\mathsf T}\mathbf{b}) = (\mathbf{a}\mathbf{a}^{\mathsf T})\mathbf{b} \end{align}

aaT\mathbf{a}\mathbf{a}^{\mathsf T}は行列です。行列については[3章][3]を参照してください。

証明をみる

a\mathbf{a}の第ii成分をaia_iとすると、

(aaT)b=(a1aTanaT)b=(a1aTbanaTb)=a(aTb)\begin{align} (\mathbf{a}\mathbf{a}^{\mathsf T})\mathbf{b} &= \begin{pmatrix}a_1\mathbf{a}^{\mathsf T} \\ \vdots \\ a_n\mathbf{a}^{\mathsf T}\end{pmatrix}\mathbf{b} \\ &= \begin{pmatrix}a_1\mathbf{a}^{\mathsf T}\mathbf{b} \\ \vdots \\ a_n\mathbf{a}^{\mathsf T}\mathbf{b}\end{pmatrix} \\ &= \mathbf{a} (\mathbf{a}^{\mathsf T}\mathbf{b}) \end{align}

となり、題意成立となります。

5. 外積

本節では定理や証明に[2章][2]で扱う正射影ベクトル、[3章][3]で導入する行列、[5章][5]で定義する行列式・余因子展開・置換などが登場します。適宜読み飛ばすか該当の章を参照して読み進めて下さい。

(1) 外積の定義

内積は2つのベクトルからスカラーを得る演算ですが、外積は2つのベクトルからベクトルを演算です。

定義1.22: 外積

a=(a1,a2,a3)T,b=(b1,b2,b3)TR3\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T} \in \mathbb{R}^3の外積a×b\mathbf{a} \times \mathbf{b}は以下のように定義されます。

a×b=(a2b3a3b2a3b1a1b3a1b2a2b1)\begin{align} \mathbf{a} \times \mathbf{b} = \begin{pmatrix} a_2 b_3 - a_3 b_2 \\ a_3 b_1 - a_1 b_3 \\ a_1 b_2 - a_2 b_1 \end{pmatrix} \end{align}

(2) 外積の行列式による表現

行列式の余因子展開に倣ったベクトル値ラプラス展開(成分にベクトルが含まれる行列の余因子展開)を用いると、外積は以下のように表現できます。

定理1.23: 外積と行列式

a=(a1,a2,a3)T,b=(b1,b2,b3)TR3\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T} \in \mathbb{R}^3とし、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T}とすると、

a×b=ijka1a2a3b1b2b3\begin{align} \mathbf{a}\times\mathbf{b} =\begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \end{align}

と書けます。

証明をみる

ベクトルが1行目に含まれた行列の行列式は、ベクトル値ラプラス展開によって定義されます。具体的には、第1行目に関する余因子展開として以下のように定義されます。従って、

ijka1a2a3b1b2b3=(1)1+1a2a3b2b3i+(1)1+2a1a3b1b3j+(1)1+3a1a2b1b2k=(a2b3a3b2)i(a1b3a3b1)j+(a1b2a2b1)k=(a2b3a3b2a3b1a1b3a1b2a2b1)\begin{align} \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} &= (-1)^{1+1}\begin{vmatrix}a_2 & a_3 \\ b_2 & b_3\end{vmatrix}\mathbf{i} +(-1)^{1+2}\begin{vmatrix}a_1 & a_3 \\ b_1 & b_3\end{vmatrix}\mathbf{j} +(-1)^{1+3}\begin{vmatrix}a_1 & a_2 \\ b_1 & b_2\end{vmatrix}\mathbf{k} \\ &= (a_2b_3 - a_3b_2)\mathbf{i} - (a_1b_3 - a_3b_1)\mathbf{j} + (a_1b_2 - a_2b_1)\mathbf{k} \\ &= \begin{pmatrix}a_2b_3 - a_3b_2 \\ a_3b_1 - a_1b_3 \\ a_1b_2 - a_2b_1\end{pmatrix} \end{align}

となり、題意成立となります。

ベクトル値ラプラス展開においても通常の行列式と同様に以下のような交代性・線形性・スカラー倍が成立します。これは、ベクトル値ラプラス展開の各成分が 2×2 行列式(=通常のスカラー行列式)になっているためです。2×2 行列式は行(列)ごとに線形で、行を入れ替えると符号が反転します。この性質がそのまま ベクトル値ラプラス展開の各成分に伝播することで成立します。

ベクトル値ラプラス展開の性質

a=(a1,a2,a3)T,b=(b1,b2,b3)TR3\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T} \in \mathbb{R}^3、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T}とし、

D(a,b)=ijka1a2a3b1b2b3=a2a3b2b3ia1a3b1b3j+a1a2b1b2k\begin{align} \mathcal{D}(\mathbf{a}, \mathbf{b}) &=\begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \\ &= \begin{vmatrix}a_2 & a_3 \\ b_2 & b_3\end{vmatrix}\mathbf{i} -\begin{vmatrix}a_1 & a_3 \\ b_1 & b_3\end{vmatrix}\mathbf{j} +\begin{vmatrix}a_1 & a_2 \\ b_1 & b_2\end{vmatrix}\mathbf{k} \end{align}

と定義すると、以下の性質が成り立ちます。ただし、α,βR\alpha, \beta \in \mathbb{R}、で、a1,a2,b1,b2R3\mathbf{a}_1, \mathbf{a}_2, \mathbf{b}_1, \mathbf{b}_2 \in \mathbb{R}^3とします。

  • 交代性:D(a,b)=D(b,a)\mathcal{D}(\mathbf{a}, \mathbf{b}) = -\mathcal{D}(\mathbf{b}, \mathbf{a})
  • 線形性・スカラー倍①:D(αa1+βa2,b)=αD(a1,b)+βD(a2,b)\mathcal{D}(\alpha\mathbf{a}_1+\beta\mathbf{a}_2, \mathbf{b}) = \alpha\mathcal{D}(\mathbf{a}_1, \mathbf{b}) + \beta\mathcal{D}(\mathbf{a}_2, \mathbf{b})
  • 線形性・スカラー倍②:D(a,αb1+βb2)=αD(a,b1)+βD(a,b2)\mathcal{D}(\mathbf{a}, \alpha\mathbf{b}_1+\beta\mathbf{b}_2) = \alpha\mathcal{D}(\mathbf{a}, \mathbf{b}_1) + \beta\mathcal{D}(\mathbf{a}, \mathbf{b}_2)

証明は割愛します。

(3) 外積の行列による表現

外積は行列とベクトルの積としても表現できます。

定理1.24: 外積と行列

ベクトル a=(a1,a2,a3)TR3\mathbf{a}=(a_1,a_2,a_3)^{\mathsf T} \in \mathbb{R}^3 に対し、以下の反対称行列

[a]×=(0a3a2a30a1a2a10)\begin{align} [\mathbf{a}]_\times=\begin{pmatrix} 0&-a_3&a_2\\ a_3&0&-a_1\\ -a_2&a_1&0 \end{pmatrix} \end{align}

を定めると、任意の bR3\mathbf{b}\in\mathbb{R}^3 に対して

a×b=[a]×b\begin{align} \mathbf{a}\times\mathbf{b}=[\mathbf{a}]_\times\,\mathbf{b} \end{align}

が成り立ちます。

証明をみる

b=(b1,b2,b3)T\mathbf{b}=(b_1, b_2, b_3)^{\mathsf T}とすると、

[a]×b=(0a3a2a30a1a2a10)(b1b2b3)=(a2b3a3b2a3b1a1b3a1b2a2b1)\begin{align} [\mathbf{a}]_\times \mathbf{b} &=\begin{pmatrix} 0&-a_3&a_2\\ a_3&0&-a_1\\ -a_2&a_1&0 \end{pmatrix}\begin{pmatrix}b_1\\ b_2 \\ b_3\end{pmatrix} \\ &= \begin{pmatrix}a_2b_3 - a_3b_2 \\ a_3b_1 - a_1b_3 \\ a_1b_2 - a_2b_1\end{pmatrix} \end{align}

となり、題意成立となります。

(4) 外積の性質

外積は以下の性質を満たします。いずれも成分を用いて証明ができます。

反交換法則

a,bR3\mathbf{a}, \mathbf{b} \in \mathbb{R}^3について以下が成り立ちます。

a×b=(b×a)\begin{align} \mathbf{a} \times \mathbf{b} = -(\mathbf{b} \times \mathbf{a}) \end{align}
証明をみる

a=(a1,a2,a3)T,b=(b1,b2,b3)T\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T}、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T}とすると、外積は

a×b=ijka1a2a3b1b2b3\begin{align} \mathbf{a}\times\mathbf{b} =\begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \end{align}

と書けます。ここで第2行と第3行(a(\mathbf{a}b\mathbf{b})を交換すると行列式は符号が反転するので

ijkb1b2b3a1a2a3=(a×b)\begin{align} \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ b_1 & b_2 & b_3\\ a_1 & a_2 & a_3 \end{vmatrix} =-(\mathbf{a}\times\mathbf{b}) \end{align}

ですが、左辺はまさに b×a\mathbf{b}\times\mathbf{a} です。したがって

a×b=(b×a)\begin{align} \mathbf{a}\times\mathbf{b}=-(\mathbf{b}\times\mathbf{a}) \end{align}

が成立します。

分配法則

a,b,cR3\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{R}^3について以下が成り立ちます。

a×(b+c)=a×b+a×c\begin{align} \mathbf{a} \times (\mathbf{b} + \mathbf{c}) = \mathbf{a} \times \mathbf{b} + \mathbf{a} \times \mathbf{c} \end{align}
証明をみる

a=(a1,a2,a3)T,b=(b1,b2,b3)T,c=(c1,c2,c3)T\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T}, \mathbf{c}=(c_1, c_2, c_3)^{\mathsf T}、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T}とすると、外積は

a×b=ijka1a2a3b1b2b3\begin{align} \mathbf{a}\times\mathbf{b} =\begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \end{align}

と書けます。第3行に関して行列式は線形なので

a×(b+c)=ijka1a2a3b1+c1b2+c2b3+c3=ijka1a2a3b1b2b3+ijka1a2a3c1c2c3=a×b+a×c.\begin{align} \mathbf{a} \times (\mathbf{b} + \mathbf{c}) &= \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1+c_1 & b_2+c_2 & b_3+c_3 \end{vmatrix} \\ &= \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} + \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ c_1 & c_2 & c_3 \end{vmatrix} \\ &=\mathbf{a}\times\mathbf{b}+\mathbf{a}\times\mathbf{c}. \end{align}

したがって所望の分配法則が成り立ちます。

結合法則

λR\lambda \in \mathbb{R}, a,bR3\mathbf{a}, \mathbf{b} \in \mathbb{R}^3について以下が成り立ちます。

(λa)×b=λ(a×b)\begin{align} (\lambda \mathbf{a}) \times \mathbf{b} = \lambda (\mathbf{a} \times \mathbf{b}) \end{align}
証明をみる

a=(a1,a2,a3)T,b=(b1,b2,b3)T\mathbf{a}=(a_1, a_2, a_3)^{\mathsf T}, \mathbf{b}=(b_1, b_2, b_3)^{\mathsf T}、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T}とすると、外積は

a×b=ijka1a2a3b1b2b3\begin{align} \mathbf{a}\times\mathbf{b} =\begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \end{align}

と書けます。第2行(a\mathbf{a} の行)に関して行列式は線形なので

(λa)×b=ijkλa1λa2λa3b1b2b3=λijka1a2a3b1b2b3=λ(a×b).\begin{align} (\lambda\mathbf{a})\times\mathbf{b} &= \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ \lambda a_1 & \lambda a_2 & \lambda a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \\ &=\lambda \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k}\\ a_1 & a_2 & a_3\\ b_1 & b_2 & b_3 \end{vmatrix} \\ &=\lambda(\mathbf{a}\times\mathbf{b}). \end{align}

従って、題意成立となります。

(5) 外積の公式 (1):スカラー三重積

定理1.25: スカラー三重積(混合積)

3つのベクトル a,b,cR3\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{R}^3 に対して、以下をスカラー三重積(または混合積)といいます。

a(b×c)\begin{align} \mathbf{a} \cdot (\mathbf{b} \times \mathbf{c}) \end{align}

スカラー三重積はスカラーであり、幾何学的には a,b,c\mathbf{a}, \mathbf{b}, \mathbf{c} によって張られる平行六面体の符号付き体積を表します。

証明をみる (幾何学的証明)

外積の大きさと方向(定理1.xx)より

b×c=bcsinθ\begin{align} \|\mathbf{b}\times\mathbf{c}\|=\|\mathbf{b}\|\,\|\mathbf{c}\|\,\sin\theta \end{align}

となります。ここでθ\thetab,c\mathbf{b},\mathbf{c} のなす角であり、b×c\mathbf{b}\times\mathbf{c}b,c\mathbf{b},\mathbf{c} に直交し、右手系で向きが定まります。ゆえに

a(b×c)=ab×ccosφ\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) =\|\mathbf{a}\|\,\|\mathbf{b}\times\mathbf{c}\|\,\cos\varphi \end{align}

となります。ここでφ\varphia\mathbf{a}b×c\mathbf{b}\times\mathbf{c} のなす角です。今、b×c\|\mathbf{b}\times\mathbf{c}\|底面(b,c\mathbf{b},\mathbf{c})の面積、acosφ\|\mathbf{a}\|\,\cos\varphia\mathbf{a}b×c\mathbf{b}\times\mathbf{c}方向への正射影ベクトルなので、

acosφb×c=底面積×高さ=体積.\begin{align} \|\mathbf{a}\|\,\cos\varphi \cdot \|\mathbf{b}\times\mathbf{c}\|=\text{底面積}\times\text{高さ}=\text{体積}. \end{align}

となります。符号は、{a,b,c}\{\mathbf{a},\mathbf{b},\mathbf{c}\} が右手系(正の向き)なら正、左手系(負の向き)なら負として与えられます。したがって

a(b×c)\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) \end{align}

は平行六面体の符号付き体積を与えます。

証明をみる (行列式を用いた証明)

a=(a1,a2,a3)T, b=(b1,b2,b3)T, c=(c1,c2,c3)TR3\mathbf{a}=(a_1,a_2,a_3)^{\mathsf T},\ \mathbf{b}=(b_1,b_2,b_3)^{\mathsf T},\ \mathbf{c}=(c_1,c_2,c_3)^{\mathsf T} \in \mathbb{R}^3、また、i=(1,0,0)T,j=(0,1,0)T,k=(0,0,1)T\mathbf{i}=(1, 0, 0)^{\mathsf T}, \mathbf{j}=(0, 1, 0)^{\mathsf T}, \mathbf{k}=(0, 0, 1)^{\mathsf T} とすると

a(b×c)=aijkb1b2b3c1c2c3=a(a2a3b2b3ia1a3b1b3j+a1a2b1b2k)=a1a2a3b2b3a2a1a3b1b3+a3a1a2b1b2=a1a2a3b1b2b3c1c2c3\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) &= \mathbf{a}\cdot \begin{vmatrix} \mathbf{i} & \mathbf{j} &\mathbf{k}\\ b_1 & b_2 & b_3\\ c_1 & c_2 & c_3 \end{vmatrix} \\ &= \mathbf{a}\cdot \Big( \begin{vmatrix}a_2 & a_3 \\ b_2 & b_3\end{vmatrix}\mathbf{i} -\begin{vmatrix}a_1 & a_3 \\ b_1 & b_3\end{vmatrix}\mathbf{j} +\begin{vmatrix}a_1 & a_2 \\ b_1 & b_2\end{vmatrix}\mathbf{k}\Big) \\ &= a_1\begin{vmatrix}a_2 & a_3 \\ b_2 & b_3\end{vmatrix} -a_2\begin{vmatrix}a_1 & a_3 \\ b_1 & b_3\end{vmatrix} +a_3\begin{vmatrix}a_1 & a_2 \\ b_1 & b_2\end{vmatrix}\\ &= \begin{vmatrix} a_1 & a_2 & a_3\\ b_1 & b_2 & b_3\\ c_1 & c_2 & c_3 \end{vmatrix} \end{align}

行列式は、3本のベクトルから張られる平行六面体の体積の符号付き版(3線形性・交代性・単位立方体で値1)を特徴づける唯一の関数であるため、

a(b×c)\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) \end{align}

は平行六面体の符号付き体積に等しくなります。

(6) 外積の公式 (2):スカラー三重積の巡回対称性

定理1.26: スカラー三重積(混合積)の巡回対称性

3つのベクトル a,b,cR3\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{R}^3 に対して、以下の性質が成り立ちます。

a(b×c)=b(c×a)=c(a×b)\begin{align} \mathbf{a} \cdot (\mathbf{b} \times \mathbf{c}) = \mathbf{b} \cdot (\mathbf{c} \times \mathbf{a}) = \mathbf{c} \cdot (\mathbf{a} \times \mathbf{b}) \end{align}
証明をみる

定理1.25の証明より、

a(b×c)=a1a2a3b1b2b3c1c2c3\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) = \begin{vmatrix} a_1 & a_2 & a_3\\ b_1 & b_2 & b_3\\ c_1 & c_2 & c_3 \end{vmatrix} \end{align}

行列の行を (1,2,3)(2,3,1)(1,2,3)\to(2,3,1)(1,2,3)(3,1,2)(1,2,3)\to(3,1,2) と巡回する操作は 3-巡回置換であり、これは偶置換なので行列式の値は不変。したがって

aTbTcT=bTcTaT=cTaTbT\begin{align} \begin{vmatrix} \mathbf{a}^{\mathsf T}\\ \mathbf{b}^{\mathsf T}\\ \mathbf{c}^{\mathsf T} \end{vmatrix} = \begin{vmatrix} \mathbf{b}^{\mathsf T}\\ \mathbf{c}^{\mathsf T}\\ \mathbf{a}^{\mathsf T} \end{vmatrix} = \begin{vmatrix} \mathbf{c}^{\mathsf T}\\ \mathbf{a}^{\mathsf T}\\ \mathbf{b}^{\mathsf T} \end{vmatrix} \end{align}

すなわち

a(b×c)=b(c×a)=c(a×b)\begin{align} \mathbf{a}\cdot(\mathbf{b}\times\mathbf{c}) = \mathbf{b}\cdot(\mathbf{c}\times\mathbf{a}) = \mathbf{c}\cdot(\mathbf{a}\times\mathbf{b}) \end{align}

となります。

(7) 外積の公式 (3):ベクトル三重積

定理1.27: ベクトル三重積

3つのベクトル a,b,cR3\mathbf{a}, \mathbf{b}, \mathbf{c} \in \mathbb{R}^3 に対して、以下をベクトル三重積といいます。

a×(b×c)\begin{align} \mathbf{a} \times (\mathbf{b} \times \mathbf{c}) \end{align}

ベクトル三重積はベクトルであり、以下の性質が成り立ちます。

a×(b×c)=(ac)b(ab)c\begin{align} \mathbf{a} \times (\mathbf{b} \times \mathbf{c}) = (\mathbf{a} \cdot \mathbf{c})\mathbf{b} - (\mathbf{a} \cdot \mathbf{b})\mathbf{c} \end{align}
証明をみる

ベクトル a=(a1,a2,a3)T\mathbf{a}=(a_1,a_2,a_3)^{\mathsf T}b=(b1,b2,b3)T\mathbf{b}=(b_1,b_2,b_3)^{\mathsf T}c=(c1,c2,c3)TR3\mathbf{c}=(c_1,c_2,c_3)^{\mathsf T} \in \mathbb{R}^3 とします。まず b×c\mathbf{b}\times\mathbf{c} の成分は外積の定義より

b×c=(b2c3b3c2b3c1b1c3b1c2b2c1)\begin{align} \mathbf{b}\times\mathbf{c}= \begin{pmatrix} b_2c_3-b_3c_2\\ b_3c_1-b_1c_3\\ b_1c_2-b_2c_1 \end{pmatrix} \end{align}

です。

次に,a×(b×c)\mathbf{a}\times(\mathbf{b}\times\mathbf{c}) を計算します。

a×(b×c)=(a2(b1c2b2c1)a3(b3c1b1c3)a3(b2c3b3c2)a1(b1c3b3c1)a1(b3c1b1c3)a2(b2c1b1c2))\begin{align} \mathbf{a}\times(\mathbf{b}\times\mathbf{c}) &=\begin{pmatrix} a_2(b_1c_2-b_2c_1)-a_3(b_3c_1-b_1c_3)\\ a_3(b_2c_3-b_3c_2)-a_1(b_1c_3-b_3c_1)\\ a_1(b_3c_1-b_1c_3)-a_2(b_2c_1-b_1c_2) \end{pmatrix} \end{align}

各成分を整理する。

(1) 第1成分

a2(b1c2b2c1)a3(b3c1b1c3)=a2b1c2a2b2c1a3b3c1+a3b1c3\begin{align} a_2(b_1c_2-b_2c_1)-a_3(b_3c_1-b_1c_3) =a_2b_1c_2-a_2b_2c_1-a_3b_3c_1+a_3b_1c_3 \end{align}

これを (ac)b(ab)c(\mathbf{a}\cdot\mathbf{c})\mathbf{b}-(\mathbf{a}\cdot\mathbf{b})\mathbf{c} の第1成分と比較します。

  • (ac)=a1c1+a2c2+a3c3(\mathbf{a}\cdot\mathbf{c})=a_1c_1+a_2c_2+a_3c_3
  • (ab)=a1b1+a2b2+a3b3(\mathbf{a}\cdot\mathbf{b})=a_1b_1+a_2b_2+a_3b_3

したがって、第1成分は

(ac)b1(ab)c1=b1(a1c1+a2c2+a3c3)c1(a1b1+a2b2+a3b3)=a2b1c2+a3b1c3a2b2c1a3b3c1\begin{align} (\mathbf{a}\cdot\mathbf{c})b_1-(\mathbf{a}\cdot\mathbf{b})c_1 &= b_1(a_1c_1+a_2c_2+a_3c_3)-c_1(a_1b_1+a_2b_2+a_3b_3) \\ &= a_2b_1c_2+a_3b_1c_3 - a_2b_2c_1 - a_3b_3c_1 \end{align}

となり、これは上で計算した第1成分と一致します。

(2) 第2成分

同様に計算すると

a3b2c3a3b3c2a1b1c3+a1b3c1\begin{align} a_3b_2c_3-a_3b_3c_2-a_1b_1c_3+a_1b_3c_1 \end{align}

となります。

一方,式 (ac)b2(ab)c2(\mathbf{a}\cdot\mathbf{c})b_2-(\mathbf{a}\cdot\mathbf{b})c_2 の第2成分は

b2(a1c1+a2c2+a3c3)c2(a1b1+a2b2+a3b3)=a3b2c3+a1b2c1a1b1c2a3b3c2\begin{align} &b_2(a_1c_1+a_2c_2+a_3c_3)-c_2(a_1b_1+a_2b_2+a_3b_3) \\ &= a_3b_2c_3+a_1b_2c_1 - a_1b_1c_2 - a_3b_3c_2 \end{align}

これも一致します。

(3) 第3成分

同様に計算すると

a1b3c1a1b1c3a2b2c1+a2b1c2\begin{align} a_1b_3c_1-a_1b_1c_3-a_2b_2c_1+a_2b_1c_2 \end{align}

となります。

一方,式 (ac)b3(ab)c3(\mathbf{a}\cdot\mathbf{c})b_3-(\mathbf{a}\cdot\mathbf{b})c_3 の第3成分は

b3(a1c1+a2c2+a3c3)c3(a1b1+a2b2+a3b3)=a1b3c1+a2b3c2a1b1c3a2b2c3,\begin{align} &b_3(a_1c_1+a_2c_2+a_3c_3)-c_3(a_1b_1+a_2b_2+a_3b_3) \\ &= a_1b_3c_1+a_2b_3c_2 - a_1b_1c_3 - a_2b_2c_3, \end{align}

これも一致します。

(1)(2)(3)より3つの成分が全て一致するので

a×(b×c)=(ac)b(ab)c\begin{align} \mathbf{a}\times(\mathbf{b}\times\mathbf{c})=(\mathbf{a}\cdot\mathbf{c})\mathbf{b}-(\mathbf{a}\cdot\mathbf{b})\mathbf{c} \end{align}

が成り立ちます。

(8) 外積の公式 (4):外積の恒等式

定理1.28: 外積の恒等式

u,v,x,yR3\mathbf{u}, \mathbf{v}, \mathbf{x}, \mathbf{y} \in \mathbb{R}^3について、以下の高等式が成り立ちます。

(u×v)(x×y)=(ux)(vy)(uy)(vx)\begin{align} (\mathbf u\times \mathbf v)\cdot(\mathbf x\times \mathbf y) = (\mathbf u\cdot \mathbf x)(\mathbf v\cdot \mathbf y)-(\mathbf u\cdot \mathbf y)(\mathbf v\cdot \mathbf x) \end{align}
証明をみる

ベクトル三重積(BAC–CAB)の公式

a×(b×c)=b(ac)c(ab)\begin{align} \mathbf a\times(\mathbf b\times \mathbf c) = \mathbf b(\mathbf a\cdot \mathbf c)-\mathbf c(\mathbf a\cdot \mathbf b) \end{align}

を用います。まず以下のように変形します。

(u×v)(x×y)=(x×y)(u×v)=u(v×(x×y)).\begin{align} (\mathbf u\times \mathbf v)\cdot(\mathbf x\times \mathbf y) &= (\mathbf x\times \mathbf y)\cdot(\mathbf u\times \mathbf v) \\ &= \mathbf u\cdot\bigl(\mathbf v\times(\mathbf x\times \mathbf y)\bigr). \end{align}

第2式から第3式への変形には、

c(a×b)=a(b×c)\begin{align} \mathbf{c} \cdot (\mathbf{a} \times \mathbf{b}) = \mathbf{a} \cdot (\mathbf{b} \times \mathbf{c}) \end{align}

のスカラー3重積の公式を、c=x×y\mathbf{c}=\mathbf{x}\times\mathbf{y}, a=u\mathbf{a}=\mathbf{u}, b=v\mathbf{b}=\mathbf{v}とおいて利用しています。ここで BAC–CAB を a=v, b=x, c=y\mathbf a=\mathbf v,\ \mathbf b=\mathbf x,\ \mathbf c=\mathbf y に適用すると

v×(x×y)=x(vy)y(vx)\begin{align} \mathbf v\times(\mathbf x\times \mathbf y)=\mathbf x(\mathbf v\cdot \mathbf y)-\mathbf y(\mathbf v\cdot \mathbf x) \end{align}

したがって

(u×v)(x×y)=u(x(vy)y(vx))=(ux)(vy)(uy)(vx)\begin{align} (\mathbf u\times \mathbf v)\cdot(\mathbf x\times \mathbf y) &=\mathbf u\cdot\bigl(\mathbf x(\mathbf v\cdot \mathbf y)-\mathbf y(\mathbf v\cdot \mathbf x)\bigr)\\ &=(\mathbf u\cdot \mathbf x)(\mathbf v\cdot \mathbf y)-(\mathbf u\cdot \mathbf y)(\mathbf v\cdot \mathbf x) \end{align}

求める等式が得られます。

(9) 外積の公式 (5):行列の恒等式

定理1.29: 行列の恒等式

aR3\mathbf{a} \in \mathbf{R}^3、また、IIR3×3R^{3\times 3}の単位行列としたとき、以下の行列[a]×[\mathbf{a}]_\timesについて、

[a]×=(0a3a2a30a1a2a10)\begin{align} [\mathbf{a}]_\times=\begin{pmatrix} 0&-a_3&a_2\\ a_3&0&-a_1\\ -a_2&a_1&0 \end{pmatrix} \end{align}

以下の恒等式が成り立ちます。

[a]×T[a]×=(aa)IaaT\begin{align} [\mathbf{a}]_\times^{\mathsf T}[\mathbf{a}]_\times=(\mathbf{a}\cdot\mathbf{a})I-\mathbf{a}\mathbf{a}^{\mathsf T} \end{align}
証明をみる

a,xR3\mathbf{a}, \mathbf{x} \in \mathbb{R}^3についてベクトル三重積の公式 a×(a×x)=a(ax)x(aa)\mathbf{a}\times(\mathbf{a}\times\mathbf{x})=\mathbf{a}(\mathbf{a}\cdot\mathbf{x})-\mathbf{x}(\mathbf{a}\cdot\mathbf{a}) を行列形で書くと

[a]×2x=a(aTx)(aa)x\begin{align} [\mathbf{a}]_\times^2 \mathbf{x}=\mathbf{a}(\mathbf{a}^{\mathsf T}\mathbf{x})-(\mathbf{a}\cdot\mathbf{a})\mathbf{x} \end{align}

となります。定理1.21の内積の公式(1)から、

a(aTx)=(aaT)x\begin{align} \mathbf{a}(\mathbf{a}^{\mathsf T}\mathbf{x})=(\mathbf{a}\mathbf{a}^{\mathsf T})\mathbf{x} \end{align}

となり、

[a]×2x=(aaT)x(aa)x\begin{align} [\mathbf{a}]_\times^2 \mathbf{x}=(\mathbf{a}\mathbf{a}^{\mathsf T})\mathbf{x}-(\mathbf{a}\cdot\mathbf{a})\mathbf{x} \end{align}

がわかります。これが任意のx\mathbf{x}について成立するため、

[a]×2=aaT(aa)I\begin{align} [\mathbf{a}]_\times^2=\mathbf{a}\mathbf{a}^{\mathsf T}-(\mathbf{a}\cdot\mathbf{a})I \end{align}

となります。ここで、半対称行列の性質である[a]×T=[a]×[\mathbf{a}]_\times^{\mathsf T}=-[\mathbf{a}]_\times を用いると

[a]×T[a]×=(aa)IaaT\begin{align} [\mathbf{a}]_\times^{\mathsf T}[\mathbf{a}]_\times=(\mathbf{a}\cdot\mathbf{a})I-\mathbf{a}\mathbf{a}^{\mathsf T} \end{align}

が得られます。

(10) 外積の公式 (6):ラグランジュの恒等式

定理1.30: ラグランジュの恒等式

a,bR3\mathbf{a}, \mathbf{b} \in \mathbb{R}^3について、以下の恒等式が成立します。

a×b2=a2b2(ab)2\begin{align} \|\mathbf{a}\times\mathbf{b}\|^2=\|\mathbf{a}\|^2\|\mathbf{b}\|^2-(\mathbf{a}\cdot\mathbf{b})^2 \end{align}
証明をみる

定理1.29の行列の恒等式より、

a×b2=(a×b)(a×b)=([a]×b)T([a]×b)=bT([a]×T[a]×)b=bT((aa)IaaT)b=(aa)bTIbbTaaTb=(aa)bTb(aTb)T(aTb)=(aa)(bb)(ab)2=a2b2(ab)2\begin{align} \|\mathbf{a}\times\mathbf{b}\|^2 &=(\mathbf{a}\times\mathbf{b})\cdot(\mathbf{a}\times\mathbf{b}) \\ &=\big([\mathbf{a}]_\times\mathbf{b}\big)^{\mathsf T}\big([\mathbf{a}]_\times\mathbf{b}\big) \\ &=\mathbf{b}^{\mathsf T}\big([\mathbf{a}]_\times^{\mathsf T}[\mathbf{a}]_\times\big)\mathbf{b} \\ &=\mathbf{b}^{\mathsf T}\big((\mathbf{a}\cdot\mathbf{a})I-\mathbf{a}\mathbf{a}^{\mathsf T}\big)\mathbf{b} \\ &=(\mathbf{a}\cdot\mathbf{a})\mathbf{b}^{\mathsf T}I\mathbf{b}-\mathbf{b}^{\mathsf T}\mathbf{a}\mathbf{a}^{\mathsf T}\mathbf{b} \\ &=(\mathbf{a}\cdot\mathbf{a})\mathbf{b}^{\mathsf T}\mathbf{b}-(\mathbf{a}^{\mathsf T}\mathbf{b})^{\mathsf T}(\mathbf{a}^{\mathsf T}\mathbf{b}) \\ &=(\mathbf{a}\cdot\mathbf{a})(\mathbf{b}\cdot\mathbf{b})-(\mathbf{a}\cdot\mathbf{b})^2 \\ &= \|\mathbf{a}\|^2\|\mathbf{b}\|^2-(\mathbf{a}\cdot\mathbf{b})^2 \end{align}

となり、求める恒等式が得られます。

(11) 外積の向き

定理:1.31: 外積の向き

a,bR3\mathbf{a}, \mathbf{b} \in \mathbb{R}^3について、a×b\mathbf{a} \times \mathbf{b}の方向はa\mathbf{a}b\mathbf{b}へ重ねるように右ねじを回した際の方向です (回す角度はπ\pi以下)。外積はa\mathbf{a}b\mathbf{b}が含まれる平面に垂直な方向になります。

証明をみる (成分計算による証明)

a=(a1,a2,a3)T\mathbf{a}=(a_1,a_2,a_3)^{\mathsf T}b=(b1,b2,b3)T\mathbf{b}=(b_1,b_2,b_3)^{\mathsf T} とします。外積の成分は

a×b=(a2b3a3b2, a3b1a1b3, a1b2a2b1)T\begin{align} \mathbf{a}\times\mathbf{b}=\bigl(a_2b_3-a_3b_2,\ a_3b_1-a_1b_3,\ a_1b_2-a_2b_1\bigr)^{\mathsf T} \end{align}

まず、(1)a(a×b)=0\mathbf{a}\cdot(\mathbf{a}\times\mathbf{b})=0 を示します。 左辺を成分で展開すると

a(a×b)=a1(a2b3a3b2)+a2(a3b1a1b3)+a3(a1b2a2b1)=a1a2b3a1a3b2+a2a3b1a1a2b3+a1a3b2a2a3b1=0\begin{align} \mathbf{a}\cdot(\mathbf{a}\times\mathbf{b}) &=a_1(a_2b_3-a_3b_2)+a_2(a_3b_1-a_1b_3)+a_3(a_1b_2-a_2b_1) \\ &=a_1a_2b_3-a_1a_3b_2+a_2a_3b_1-a_1a_2b_3+a_1a_3b_2-a_2a_3b_1 \\ &=0 \end{align}

続いて (2)b(a×b)=0\mathbf{b}\cdot(\mathbf{a}\times\mathbf{b})=0 を示します。

b(a×b)=b1(a2b3a3b2)+b2(a3b1a1b3)+b3(a1b2a2b1)=a2b1b3a3b1b2+a3b1b2a1b2b3+a1b2b3a2b1b3=0\begin{align} \mathbf{b}\cdot(\mathbf{a}\times\mathbf{b}) &=b_1(a_2b_3-a_3b_2)+b_2(a_3b_1-a_1b_3)+b_3(a_1b_2-a_2b_1) \\ &=a_2b_1b_3-a_3b_1b_2+a_3b_1b_2-a_1b_2b_3+a_1b_2b_3-a_2b_1b_3 \\ &=0 \end{align}

(1)(2)より

a(a×b)=0,b(a×b)=0\begin{align} \mathbf{a}\cdot(\mathbf{a}\times\mathbf{b})=0,\quad \mathbf{b}\cdot(\mathbf{a}\times\mathbf{b})=0 \end{align}

したがって a×b\mathbf{a}\times\mathbf{b}a\mathbf{a}b\mathbf{b} の両方に直交します。

証明をみる (行列式を用いた証明)

スカラー三重積の性質を用いると

a(a×b)=det[a,a,b]\begin{align} \mathbf{a}\cdot(\mathbf{a}\times\mathbf{b})=\det[\mathbf{a},\mathbf{a},\mathbf{b}] \end{align}

となります。行列式は列(または行)が同一なら 00 なので

det[a,a,b]=0\begin{align} \det[\mathbf{a},\mathbf{a},\mathbf{b}]=0 \end{align}

同様に

b(a×b)=det[b,a,b]=0\begin{align} \mathbf{b}\cdot(\mathbf{a}\times\mathbf{b})=\det[\mathbf{b},\mathbf{a},\mathbf{b}]=0 \end{align}

となります。よって直交性が示されました。

(12) 外積のノルム

定理1.32: 外積のノルム

a,bR3\mathbf{a}, \mathbf{b} \in \mathbb{R}^3、また、これらのベクトルの成す角をθ\thetaとすると、外積のノルムは以下です。

a×b=absinθ\begin{align} \|\mathbf{a} \times \mathbf{b}\| = \|\mathbf{a}\| \|\mathbf{b}\| \sin \theta \end{align}

外積の大きさはa\mathbf{a}b\mathbf{b}が作る平行四辺形の面積に相当します。

証明をみる

定理1.30のラグランジュの恒等式より、

a×b2=a2b2(ab)2\begin{align} \|\mathbf{a}\times\mathbf{b}\|^2=\|\mathbf{a}\|^2\|\mathbf{b}\|^2-(\mathbf{a}\cdot\mathbf{b})^2 \end{align}

です。定理1.16の内積とノルムの関係式ab=abcosθ\mathbf{a}\cdot\mathbf{b}=\|\mathbf{a}\|\,\|\mathbf{b}\|\cos\theta を代入すると

a×b2=a2b2(abcosθ)2=a2b2(1cos2θ)=a2b2sin2θ\begin{align} \|\mathbf{a}\times\mathbf{b}\|^2 &=\|\mathbf{a}\|^2\|\mathbf{b}\|^2-\big(\|\mathbf{a}\|\,\|\mathbf{b}\|\cos\theta\big)^2 \\ &=\|\mathbf{a}\|^2\|\mathbf{b}\|^2\big(1-\cos^2\theta\big) \\ &=\|\mathbf{a}\|^2\|\mathbf{b}\|^2\sin^2\theta \end{align}

両辺は非負なので平方根を取って

a×b=absinθ\begin{align} \|\mathbf{a}\times\mathbf{b}\|=\|\mathbf{a}\|\,\|\mathbf{b}\|\sin\theta \end{align}

が得られます。

6. 1次独立と1次従属

(1) 1次独立

定義1.33: 1次独立

ベクトル a1,a2,,akCn\mathbf{a}_1, \mathbf{a}_2, \dots, \mathbf{a}_k \in \mathbb{C}^n が 1次独立であるとは、以下の線形結合がゼロベクトルになるとき、すべての係数が 0 に限ることを言います。

c1a1+c2a2++ckak=0c1=c2==ck=0\begin{align} c_1 \mathbf{a}_1 + c_2 \mathbf{a}_2 + \cdots + c_k \mathbf{a}_k = \mathbf{0} \Rightarrow c_1 = c_2 = \cdots = c_k = 0 \end{align}

すなわち、ベクトル同士が互いに冗長でなく、それぞれが独立した方向を持つこと(=どんな2つのベクトルも平行ではないこと)を意味します。

なお、2つのベクトルが1次独立の場合、a/ ⁣/ ⁣ ⁣ ⁣ ⁣\b\mathbf{a} /\!/\!\!\!\!\backslash\, \mathbf{b}と表記することがあります。

(2) 1次従属

定義1.34: 1次従属

ベクトル a1,a2,,akCn\mathbf{a}_1, \mathbf{a}_2, \dots, \mathbf{a}_k \in \mathbb{C}^n が 1次従属であるとは、ゼロでない係数を用いた線形結合がゼロベクトルになる場合です。

c1a1+c2a2++ckak=0,ただし i such that ci0\begin{align} c_1 \mathbf{a}_1 + c_2 \mathbf{a}_2 + \cdots + c_k \mathbf{a}_k = \mathbf{0}, \quad \text{ただし } \exists i \text{ such that } c_i \neq 0 \end{align}

つまり、少なくとも1つのベクトルが他のベクトルの線形結合として表されることを意味します。

なお、2つのベクトルが1次独立の場合、a/ ⁣/b\mathbf{a} /\!/ \mathbf{b}と表記することがあります。

7. 関連書籍

線形代数 改訂版

クイズ

広告を取得中...
前の記事なし