659.47 쿼터니언의 수학적 정의

1. 개요

쿼터니언(quaternion)은 4차원 결합 나눗셈 대수(associative division algebra)의 원소로서, 실수의 확장 체계 중 하나이다. 로봇공학에서 쿼터니언은 3차원 회전을 특이점 없이 표현하는 수학적 도구로 활용된다. 본 절에서는 쿼터니언의 엄밀한 수학적 정의, 대수적 구조, 그리고 다양한 표현 형식을 체계적으로 기술한다.

2. 해밀턴 쿼터니언의 정의

2.1 기본 정의

쿼터니언은 실수체 $\mathbb{R}$ 위에 세 개의 허수 단위 $i, j, k$ 를 도입하여 구성한 4차원 수 체계이다. 쿼터니언 전체의 집합을 $\mathbb{H}$ 로 표기하며, 이는 해밀턴(Hamilton)의 이니셜에서 유래하였다. 임의의 쿼터니언 $\mathbf{q} \in \mathbb{H}$ 는 다음과 같이 정의된다:

$\mathbf{q} = q_0 + q_1 i + q_2 j + q_3 k, \quad q_0, q_1, q_2, q_3 \in \mathbb{R}$

여기서 $q_0$ 는 스칼라부(scalar part), $(q_1, q_2, q_3)$ 는 **벡터부(vector part)**라 한다. 허수 단위 $i, j, k$ 는 다음의 **해밀턴 기본 관계(Hamilton’s fundamental relations)**를 만족한다:

$i^2 = j^2 = k^2 = ijk = -1$

이 기본 관계로부터 허수 단위 간의 모든 곱셈 규칙이 유도된다.

2.2 허수 단위의 곱셈 규칙

해밀턴의 기본 관계 $i^2 = j^2 = k^2 = ijk = -1$ 로부터 다음의 곱셈표를 도출할 수 있다:

$\times$	$1$	$i$	$j$	$k$
$1$	$1$	$i$	$j$	$k$
$i$	$i$	$-1$	$k$	$-j$
$j$	$j$	$-k$	$-1$	$i$
$k$	$k$	$j$	$-i$	$-1$

이 곱셈표를 통해 확인할 수 있는 핵심적인 성질은 다음과 같다:

순방향 순환 곱셈:

$ij = k, \quad jk = i, \quad ki = j$

역방향 순환 곱셈:

$ji = -k, \quad kj = -i, \quad ik = -j$

이 규칙은 벡터 외적(cross product)의 순환 관계와 유사한 구조를 가진다. $i, j, k$ 를 각각 $\hat{\mathbf{x}}, \hat{\mathbf{y}}, \hat{\mathbf{z}}$ 단위 벡터에 대응시키면, 순방향 곱은 오른손 법칙에 따른 외적과 일치하고, 역방향 곱은 부호가 반전된다.

2.3 비가환성 (Non-commutativity)

쿼터니언 곱셈의 가장 중요한 특성은 비가환성이다:

$\mathbf{q}_1 \mathbf{q}_2 \neq \mathbf{q}_2 \mathbf{q}_1 \quad \text{(일반적으로)}$

이 성질은 $ij = k \neq ji = -k$ 에서 직접 확인된다. 비가환성은 쿼터니언이 3차원 회전을 표현하는 데 적합한 구조를 갖는 근본적인 이유이다. 3차원 회전 자체가 비가환적인 연산이기 때문이다. 예를 들어, $x$ 축 회전 후 $y$ 축 회전을 수행한 결과는 $y$ 축 회전 후 $x$ 축 회전을 수행한 결과와 일반적으로 다르다.

3. 쿼터니언의 대수적 구조

3.1 $\mathbb{H}$ 의 벡터 공간 구조

쿼터니언의 집합 $\mathbb{H}$ 는 기저 $\{1, i, j, k\}$ 를 가지는 $\mathbb{R}$ 위의 4차원 벡터 공간이다. 따라서 쿼터니언의 덧셈과 실수 스칼라 곱은 성분별(component-wise)로 정의된다:

덧셈:

$\mathbf{p} + \mathbf{q} = (p_0 + q_0) + (p_1 + q_1)i + (p_2 + q_2)j + (p_3 + q_3)k$

실수 스칼라 곱:

$\alpha \mathbf{q} = \alpha q_0 + \alpha q_1 i + \alpha q_2 j + \alpha q_3 k, \quad \alpha \in \mathbb{R}$

이 연산들은 교환법칙, 결합법칙, 분배법칙을 만족하며, $\mathbb{H}$ 는 $\mathbb{R}^4$ 와 동형(isomorphic)인 벡터 공간을 이룬다.

3.2 $\mathbb{H}$ 의 나눗셈 대수 구조

쿼터니언의 곱셈은 두 쿼터니언 $\mathbf{p} = p_0 + p_1 i + p_2 j + p_3 k$ 와 $\mathbf{q} = q_0 + q_1 i + q_2 j + q_3 k$ 에 대해, 분배법칙과 허수 단위의 곱셈 규칙을 적용하여 다음과 같이 전개된다:

$$
\mathbf{p}\mathbf{q} = (p_0 q_0 - p_1 q_1 - p_2 q_2 - p_3 q_3) \

(p_0 q_1 + p_1 q_0 + p_2 q_3 - p_3 q_2) i \
(p_0 q_2 - p_1 q_3 + p_2 q_0 + p_3 q_1) j \
(p_0 q_3 + p_1 q_2 - p_2 q_1 + p_3 q_0) k
$$

이 곱셈은 다음의 대수적 성질을 만족한다:

성질	수식	만족 여부
결합법칙	$(\mathbf{p}\mathbf{q})\mathbf{r} = \mathbf{p}(\mathbf{q}\mathbf{r})$	만족
분배법칙	$\mathbf{p}(\mathbf{q} + \mathbf{r}) = \mathbf{p}\mathbf{q} + \mathbf{p}\mathbf{r}$	만족
교환법칙	$\mathbf{p}\mathbf{q} = \mathbf{q}\mathbf{p}$	불만족
항등원 존재	$\mathbf{q} \cdot 1 = 1 \cdot \mathbf{q} = \mathbf{q}$	만족
역원 존재	$\mathbf{q}\mathbf{q}^{-1} = \mathbf{q}^{-1}\mathbf{q} = 1$ ( $\mathbf{q} \neq 0$ )	만족

$\mathbb{H}$ 는 교환법칙을 제외한 모든 체(field)의 공리를 만족하므로, 비가환 나눗셈 대수(non-commutative division algebra) 또는 **비가환체(skew field)**로 분류된다. 프로베니우스 정리(Frobenius theorem)에 의하면, $\mathbb{R}$ 위의 유한 차원 결합 나눗셈 대수는 $\mathbb{R}$ , $\mathbb{C}$ , $\mathbb{H}$ 세 가지뿐이다.

4. 쿼터니언의 다양한 표현 형식

4.1 순서쌍 표현 (Ordered Pair Representation)

쿼터니언을 스칼라부와 벡터부의 순서쌍으로 표현하는 방식이다:

$\mathbf{q} = (s, \mathbf{v}), \quad s \in \mathbb{R}, \quad \mathbf{v} \in \mathbb{R}^3$

여기서 $s = q_0$ 이고 $\mathbf{v} = (q_1, q_2, q_3)^T$ 이다. 이 표현에서 쿼터니언 곱셈은 다음과 같이 간결하게 표기된다:

$\mathbf{p}\mathbf{q} = (s_p, \mathbf{v}_p)(s_q, \mathbf{v}_q) = (s_p s_q - \mathbf{v}_p \cdot \mathbf{v}_q, \; s_p \mathbf{v}_q + s_q \mathbf{v}_p + \mathbf{v}_p \times \mathbf{v}_q)$

여기서 $\mathbf{v}_p \cdot \mathbf{v}_q$ 는 벡터 내적(dot product), $\mathbf{v}_p \times \mathbf{v}_q$ 는 벡터 외적(cross product)이다. 이 공식은 쿼터니언 곱셈과 벡터 연산 간의 관계를 명확히 보여준다. 외적 항 $\mathbf{v}_p \times \mathbf{v}_q$ 이 곱셈의 비가환성을 유발하는 원인이다.

4.2 행렬 표현 (Matrix Representation)

쿼터니언은 실수 행렬로도 표현할 수 있다. 쿼터니언 $\mathbf{q} = q_0 + q_1 i + q_2 j + q_3 k$ 에 대응하는 $4 \times 4$ 실수 행렬 표현은 다음과 같다:

$\mathbf{q} \mapsto \begin{pmatrix} q_0 & -q_1 & -q_2 & -q_3 \\ q_1 & q_0 & -q_3 & q_2 \\ q_2 & q_3 & q_0 & -q_1 \\ q_3 & -q_2 & q_1 & q_0 \end{pmatrix}$

이 행렬 표현에서 쿼터니언의 곱셈은 행렬의 곱셈에 대응한다. 이 표현은 쿼터니언 연산을 행렬 연산으로 환원하여 분석할 수 있게 해준다.

또한 $2 \times 2$ 복소 행렬에 의한 표현도 가능하다:

$\mathbf{q} \mapsto \begin{pmatrix} q_0 + q_1 i_{\mathbb{C}} & q_2 + q_3 i_{\mathbb{C}} \\ -q_2 + q_3 i_{\mathbb{C}} & q_0 - q_1 i_{\mathbb{C}} \end{pmatrix}$

여기서 $i_{\mathbb{C}}$ 는 복소수의 허수 단위이다. 이 표현에서 단위 쿼터니언은 $SU(2)$ 군(특수 유니타리 군)의 원소에 대응하며, $SU(2)$ 와 $SO(3)$ 사이의 2대1 준동형(double covering homomorphism)이 성립한다.

4.3 차원 벡터 표현

계산적 관점에서 쿼터니언은 $\mathbb{R}^4$ 의 벡터로 표현할 수 있다:

$\mathbf{q} = \begin{pmatrix} q_0 \\ q_1 \\ q_2 \\ q_3 \end{pmatrix} \in \mathbb{R}^4$

이 표현에서 쿼터니언 곱셈 $\mathbf{r} = \mathbf{p}\mathbf{q}$ 는 왼쪽 곱셈 행렬 $L(\mathbf{p})$ 를 이용하여 다음과 같이 행렬-벡터 곱으로 표현된다:

$\mathbf{r} = L(\mathbf{p}) \mathbf{q} = \begin{pmatrix} p_0 & -p_1 & -p_2 & -p_3 \\ p_1 & p_0 & -p_3 & p_2 \\ p_2 & p_3 & p_0 & -p_1 \\ p_3 & -p_2 & p_1 & p_0 \end{pmatrix} \begin{pmatrix} q_0 \\ q_1 \\ q_2 \\ q_3 \end{pmatrix}$

유사하게 오른쪽 곱셈 행렬 $R(\mathbf{q})$ 를 이용하면:

$\mathbf{r} = R(\mathbf{q}) \mathbf{p} = \begin{pmatrix} q_0 & -q_1 & -q_2 & -q_3 \\ q_1 & q_0 & q_3 & -q_2 \\ q_2 & -q_3 & q_0 & q_1 \\ q_3 & q_2 & -q_1 & q_0 \end{pmatrix} \begin{pmatrix} p_0 \\ p_1 \\ p_2 \\ p_3 \end{pmatrix}$

두 행렬 $L(\mathbf{p})$ 와 $R(\mathbf{q})$ 의 차이는 곱셈의 비가환성을 반영한다.

5. 쿼터니언의 기본 연산

5.1 켤레 (Conjugate)

쿼터니언 $\mathbf{q} = q_0 + q_1 i + q_2 j + q_3 k$ 의 켤레는 벡터부의 부호를 반전한 것이다:

$\mathbf{q}^{*} = q_0 - q_1 i - q_2 j - q_3 k = (q_0, -\mathbf{v})$

켤레는 다음의 성질을 만족한다:

$(\mathbf{q}^{*})^{*} = \mathbf{q}, \quad (\mathbf{p}\mathbf{q})^{*} = \mathbf{q}^{*}\mathbf{p}^{*}$

켤레의 곱셈 순서가 뒤집히는 성질 $(\mathbf{p}\mathbf{q})^{*} = \mathbf{q}^{*}\mathbf{p}^{*}$ 은 전치 행렬의 곱 $(AB)^T = B^T A^T$ 과 유사하다.

5.2 노름 (Norm)

쿼터니언의 노름은 다음과 같이 정의된다:

$\|\mathbf{q}\| = \sqrt{\mathbf{q}\mathbf{q}^{*}} = \sqrt{q_0^2 + q_1^2 + q_2^2 + q_3^2}$

노름은 다음의 중요한 성질을 갖는다:

$\|\mathbf{p}\mathbf{q}\| = \|\mathbf{p}\| \|\mathbf{q}\|$

이 곱셈적 노름(multiplicative norm) 성질은 단위 쿼터니언의 곱이 다시 단위 쿼터니언이 됨을 보장한다. 따라서 단위 쿼터니언의 집합은 곱셈에 대해 닫혀 있으며, 이는 회전의 합성이 다시 유효한 회전이 되는 물리적 사실과 일치한다.

5.3 역 쿼터니언 (Inverse Quaternion)

영이 아닌 쿼터니언의 역은 다음과 같이 정의된다:

$\mathbf{q}^{-1} = \frac{\mathbf{q}^{*}}{\|\mathbf{q}\|^2}$

이 정의는 $\mathbf{q}\mathbf{q}^{-1} = \mathbf{q}^{-1}\mathbf{q} = 1$ 을 만족한다. 단위 쿼터니언( $\|\mathbf{q}\| = 1$ )의 경우 역은 켤레와 동일하다:

$\|\mathbf{q}\| = 1 \implies \mathbf{q}^{-1} = \mathbf{q}^{*}$

이 성질은 역회전(inverse rotation)의 계산을 매우 효율적으로 만든다. 행렬의 역행렬 계산이 $O(n^3)$ 의 복잡도를 가지는 것과 비교하면, 쿼터니언의 역은 부호 반전만으로 $O(1)$ 에 완료된다.

6. 순수 쿼터니언과 3차원 벡터

6.1 순수 쿼터니언의 정의

스칼라부가 0인 쿼터니언을 **순수 쿼터니언(pure quaternion)**이라 한다:

$\mathbf{p} = 0 + p_1 i + p_2 j + p_3 k = (0, \mathbf{v})$

3차원 공간의 벡터 $\mathbf{v} = (v_x, v_y, v_z)^T$ 는 순수 쿼터니언 $(0, \mathbf{v})$ 와 자연스럽게 대응한다. 이 대응을 통해 3차원 벡터에 대한 쿼터니언 회전 연산을 정의할 수 있다.

6.2 순수 쿼터니언의 곱

두 순수 쿼터니언 $\mathbf{p} = (0, \mathbf{u})$ 와 $\mathbf{q} = (0, \mathbf{v})$ 의 곱은 다음과 같이 전개된다:

$\mathbf{p}\mathbf{q} = (-\mathbf{u} \cdot \mathbf{v}, \; \mathbf{u} \times \mathbf{v})$

이 결과의 스칼라부는 내적의 부호가 반전된 값이고, 벡터부는 외적이다. 따라서 두 3차원 벡터의 내적과 외적이 쿼터니언 곱셈이라는 하나의 연산으로 통합된다. 이 관계는 쿼터니언 대수와 벡터 대수 사이의 심층적 연결을 보여준다.

특히, 동일한 순수 쿼터니언의 제곱은:

$\mathbf{p}^2 = (-\|\mathbf{u}\|^2, \mathbf{0})$

이므로, 단위 순수 쿼터니언의 제곱은 $-1$ 이 된다. 이는 허수 단위 $i, j, k$ 가 모두 단위 순수 쿼터니언이라는 사실과 일치한다.

7. 쿼터니언과 수 체계의 위치

쿼터니언은 수 체계의 확장 계층에서 다음과 같은 위치를 차지한다:

$\mathbb{R} \subset \mathbb{C} \subset \mathbb{H} \subset \mathbb{O}$

수 체계	차원	교환법칙	결합법칙	나눗셈	노름 곱셈성
실수 $\mathbb{R}$	1	만족	만족	가능	만족
복소수 $\mathbb{C}$	2	만족	만족	가능	만족
쿼터니언 $\mathbb{H}$	4	불만족	만족	가능	만족
팔원수 $\mathbb{O}$	8	불만족	불만족	가능	만족

케일리-딕슨 구성(Cayley-Dickson construction)에 의해, 각 단계에서 차원은 2배가 되고 하나의 대수적 성질이 상실된다. 쿼터니언은 교환법칙을 잃지만 결합법칙은 보존하며, 이 성질 덕분에 행렬 표현이 가능하고 회전군의 표현론에 활용될 수 있다.

8. 단위 쿼터니언과 회전

8.1 단위 쿼터니언의 정의

노름이 1인 쿼터니언을 **단위 쿼터니언(unit quaternion)**이라 한다:

$\|\mathbf{q}\| = \sqrt{q_0^2 + q_1^2 + q_2^2 + q_3^2} = 1$

단위 쿼터니언의 집합은 $\mathbb{R}^4$ 내의 3차원 초구면(3-sphere, $S^3$ )을 이루며, 쿼터니언 곱셈에 대해 군(group)을 형성한다. 이 군은 위상적으로 $SU(2)$ 와 동형이며, $SO(3)$ 의 이중 피복(double cover)이다.

8.2 지수 형식 (Exponential Form)

단위 쿼터니언은 다음의 지수 형식으로 표현할 수 있다:

$\mathbf{q} = e^{\frac{\theta}{2}\hat{\mathbf{n}}} = \cos\frac{\theta}{2} + \sin\frac{\theta}{2}(n_x i + n_y j + n_z k)$

여기서 $\hat{\mathbf{n}} = (n_x, n_y, n_z)$ 는 단위 벡터(회전축)이고, $\theta$ 는 회전 각도이다. 이 형식은 오일러 공식 $e^{i\alpha} = \cos\alpha + i\sin\alpha$ 의 자연스러운 일반화이다.

쿼터니언 지수 함수는 다음과 같이 정의된다. 순수 쿼터니언 $\mathbf{v} = (0, \mathbf{u})$ 에 대해:

$e^{\mathbf{v}} = \cos\|\mathbf{u}\| + \frac{\mathbf{u}}{\|\mathbf{u}\|}\sin\|\mathbf{u}\|$

역으로 쿼터니언 로그 함수는 단위 쿼터니언 $\mathbf{q} = (s, \mathbf{v})$ 에 대해:

$\ln\mathbf{q} = \frac{\mathbf{v}}{\|\mathbf{v}\|}\arccos(s)$

지수 형식은 쿼터니언의 거듭제곱, 보간, 미분 등의 연산에서 유용하게 활용된다. 특히, 쿼터니언의 $t$ 제곱은 다음과 같이 정의된다:

$\mathbf{q}^t = e^{t \ln\mathbf{q}} = \cos\left(\frac{t\theta}{2}\right) + \sin\left(\frac{t\theta}{2}\right)\hat{\mathbf{n}}$

이는 $t = 0$ 에서 항등 회전, $t = 1$ 에서 원래 회전을 나타내며, $0 \leq t \leq 1$ 범위에서 부드러운 회전 보간을 제공한다.

8.3 회전 연산의 정의

단위 쿼터니언 $\mathbf{q}$ 에 의한 벡터 $\mathbf{p} \in \mathbb{R}^3$ 의 회전은 다음과 같이 정의된다:

$L_{\mathbf{q}}(\mathbf{p}) = \mathbf{q}\hat{\mathbf{p}}\mathbf{q}^{*}$

여기서 $\hat{\mathbf{p}} = (0, \mathbf{p})$ 는 $\mathbf{p}$ 에 대응하는 순수 쿼터니언이다. 이 연산이 실제로 3D 회전임을 증명하려면 다음을 보여야 한다:

순수 쿼터니언 보존: $\text{scalar}(\mathbf{q}\hat{\mathbf{p}}\mathbf{q}^{*}) = 0$ (결과가 순수 쿼터니언이므로 3D 벡터로 해석 가능)
노름 보존: $\|\mathbf{q}\hat{\mathbf{p}}\mathbf{q}^{*}\| = \|\hat{\mathbf{p}}\|$ (벡터의 길이가 변하지 않음)
선형성: $L_{\mathbf{q}}(\alpha\mathbf{p}_1 + \beta\mathbf{p}_2) = \alpha L_{\mathbf{q}}(\mathbf{p}_1) + \beta L_{\mathbf{q}}(\mathbf{p}_2)$
행렬식 보존: 대응하는 선형 사상의 행렬식이 $+1$ (반사가 아닌 순수 회전)

이 네 가지 성질이 모두 만족됨을 보이면, $L_{\mathbf{q}}$ 는 $SO(3)$ 의 원소, 즉 고유 회전(proper rotation)에 해당한다.

9. 로봇공학에서의 수학적 의의

쿼터니언의 수학적 정의가 로봇공학에서 갖는 실질적 의의는 다음과 같다:

매개변수 효율성: 4개의 파라미터와 1개의 구속 조건( $\|\mathbf{q}\| = 1$ )으로 3자유도 회전을 표현한다. 회전 행렬(9개 파라미터, 6개 구속 조건)보다 간결하다.
특이점 부재: $S^3$ 위의 점으로 표현되므로, 오일러 각과 같은 좌표 특이점이 존재하지 않는다.
곱셈적 노름: $\|\mathbf{p}\mathbf{q}\| = \|\mathbf{p}\|\|\mathbf{q}\|$ 성질에 의해, 단위 쿼터니언의 곱은 항상 단위 쿼터니언이다. 이는 회전 합성의 수학적 정합성을 보장한다.
효율적 역연산: 단위 쿼터니언의 역이 켤레와 동일하므로, 역회전 계산이 $O(1)$ 에 완료된다.
수치 안정성: 부동소수점 연산의 누적 오차로 인해 노름이 1에서 벗어나더라도, 단순한 정규화로 복원할 수 있다.

10. 요약

쿼터니언 $\mathbb{H}$ 는 해밀턴에 의해 정의된 4차원 비가환 나눗셈 대수로서, 기저 $\{1, i, j, k\}$ 와 기본 관계 $i^2 = j^2 = k^2 = ijk = -1$ 에 의해 완전히 결정된다. 쿼터니언은 순서쌍 표현, 행렬 표현, 4차원 벡터 표현 등 다양한 형식으로 나타낼 수 있으며, 단위 쿼터니언은 지수 형식을 통해 축-각 표현과 직접적으로 연결된다. 곱셈적 노름, 결합법칙, 비가환성 등의 대수적 성질은 3차원 회전의 물리적 특성과 정확히 부합하며, 이러한 수학적 기반이 쿼터니언을 로봇공학의 표준 회전 표현 방식으로 확립하는 근거가 된다.

11. 참고 문헌

Hamilton, W. R. (1844). “On Quaternions; or on a New System of Imaginaries in Algebra.” Philosophical Magazine, 25(3), 489-495.
Hamilton, W. R. (1853). Lectures on Quaternions. Hodges and Smith, Dublin.
Kuipers, J. B. (1999). Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality. Princeton University Press.
Altmann, S. L. (1986). Rotations, Quaternions, and Double Groups. Oxford University Press.
Shoemake, K. (1985). “Animating Rotation with Quaternion Curves.” ACM SIGGRAPH Computer Graphics, 19(3), 245-254.
Diebel, J. (2006). “Representing Attitude: Euler Angles, Unit Quaternions, and Rotation Vectors.” Stanford University Technical Report.