9.65 수정된 로드리게스 파라미터(MRP)의 정의

1. MRP의 개념

수정된 로드리게스 파라미터(Modified Rodrigues Parameters, MRP)는 3차원 회전을 3개의 매개변수로 표현하는 방법 중 하나이다. 회전 축 $\hat{\mathbf{u}}$ 와 회전 각 $\phi$ 로부터 다음과 같이 정의된다.

$\boldsymbol{\sigma} = \tan\left(\frac{\phi}{4}\right)\hat{\mathbf{u}}$

이는 축-각도 표현의 변형이며, 회전 각을 $\tan(\phi/4)$ 로 압축한 형태이다. MRP는 짐벌 락이 없는 3매개변수 표현으로, 자세 추정과 자세 제어에 사용된다.

2. MRP의 역사

MRP는 항공 우주 공학과 우주선 자세 제어 분야에서 발전한 표현법이다. 슈스터(Malcolm D. Shuster)의 1993년 자세 표현 종합 논문 “A Survey of Attitude Representations“에서 표준적 형식으로 정리되었다. 이름의 “수정된“은 원래의 로드리게스 파라미터(Classical Rodrigues Parameters, CRP)와 구분하기 위해 사용된다.

3. CRP와 MRP의 비교

3.1 고전적 로드리게스 파라미터(CRP)

CRP는 깁스 벡터(Gibbs vector)라고도 하며 다음과 같이 정의된다.

$\mathbf{g} = \tan\left(\frac{\phi}{2}\right)\hat{\mathbf{u}}$

CRP는 $\phi = \pi$ 에서 발산한다(특이점). 따라서 180도 회전 근처에서 사용할 수 없다.

3.2 수정된 로드리게스 파라미터(MRP)

MRP는 다음과 같이 정의된다.

$\boldsymbol{\sigma} = \tan\left(\frac{\phi}{4}\right)\hat{\mathbf{u}}$

각도가 $\phi/4$ 로 더 압축되므로, MRP는 $\phi = 2\pi$ 에서야 발산한다. 따라서 정상적인 회전 범위 $[0, 2\pi)$ 에서 발산하지 않으며, $\phi = \pi$ 에서도 잘 정의된다.

표현법	정의	특이점
CRP	$\tan(\phi/2)\hat{\mathbf{u}}$	$\phi = \pi$
MRP	$\tan(\phi/4)\hat{\mathbf{u}}$	$\phi = 2\pi$

4. MRP의 매개변수 영역

MRP의 자유도는 3이며, 이는 회전 공간 $SO(3)$ 의 자유도와 일치한다. MRP의 크기는

$\lVert \boldsymbol{\sigma} \rVert = \tan(\phi/4)$

이다. $\phi \in [0, 2\pi)$ 에 대해 $\lVert \boldsymbol{\sigma} \rVert$ 는 $0$ 부터 $\infty$ 까지 변한다. 단, $\phi \in [0, \pi]$ 로 제한하면 $\lVert \boldsymbol{\sigma} \rVert$ 의 범위는 $[0, 1]$ 이다.

5. MRP에서 회전 행렬로의 변환

MRP $\boldsymbol{\sigma}$ 로부터 회전 행렬은 다음과 같이 계산된다.

$\mathbf{R}(\boldsymbol{\sigma}) = \mathbf{I} + \frac{8[\boldsymbol{\sigma}]_\times^2 - 4(1 - \lVert \boldsymbol{\sigma} \rVert^2)[\boldsymbol{\sigma}]_\times}{(1 + \lVert \boldsymbol{\sigma} \rVert^2)^2}$

여기서 $[\boldsymbol{\sigma}]_\times$ 는 MRP의 반대칭 행렬이다. 이는 닫힌 형태의 행렬 공식이며, 삼각 함수 계산이 필요 없다. 이것이 MRP의 주요 장점 중 하나이다.

6. 회전 행렬에서 MRP로의 변환

회전 행렬 $\mathbf{R}$ 로부터 MRP는 다음과 같이 추출된다. 먼저 회전 각 $\phi$ 와 회전 축 $\hat{\mathbf{u}}$ 를 추출한 후

$\boldsymbol{\sigma} = \tan(\phi/4)\hat{\mathbf{u}}$

또는 쿼터니언 $\mathbf{q} = (q_w, q_x, q_y, q_z)$ 를 경유하여

$\boldsymbol{\sigma} = \frac{1}{1 + q_w}(q_x, q_y, q_z)^T$

이 식이 MRP와 단위 쿼터니언의 관계이다.

7. MRP와 쿼터니언의 관계

MRP는 단위 쿼터니언으로부터 직접 도출된다. 단위 쿼터니언 $\mathbf{q} = (q_w, \mathbf{q}_v) = (\cos(\phi/2), \sin(\phi/2)\hat{\mathbf{u}})$ 로부터

$\boldsymbol{\sigma} = \frac{\mathbf{q}_v}{1 + q_w} = \frac{\sin(\phi/2)\hat{\mathbf{u}}}{1 + \cos(\phi/2)} = \tan(\phi/4)\hat{\mathbf{u}}$

여기서 반각 공식 $\tan(\phi/4) = \sin(\phi/2)/(1 + \cos(\phi/2))$ 가 사용되었다.

이 관계는 MRP를 쿼터니언의 사영으로 해석할 수 있게 한다. 4차원 단위 쿼터니언을 3차원 공간으로 입체 사영(stereographic projection)한 결과가 MRP이다.

8. 그림자 MRP (Shadow MRP)

MRP는 매니폴드 $SO(3)$ 전체를 매개화할 수 있지만, $\phi = 2\pi$ 에서 발산한다. 이를 해결하기 위해 그림자 MRP가 도입된다. 동일한 회전을 두 가지 MRP로 표현할 수 있다.

8.1 원래 MRP

$\boldsymbol{\sigma} = \tan(\phi/4)\hat{\mathbf{u}}$

8.2 그림자 MRP

$\boldsymbol{\sigma}^S = -\frac{1}{\lVert \boldsymbol{\sigma} \rVert^2}\boldsymbol{\sigma} = -\tan((\pi - \phi/2)/2)\hat{\mathbf{u}}$

또는 동등하게

$\boldsymbol{\sigma}^S = -\cot(\phi/4)\hat{\mathbf{u}}$

원래 MRP와 그림자 MRP는 같은 회전을 나타내지만, 한쪽이 발산하면 다른 쪽이 잘 정의된다. 자세 추정에서 $\lVert \boldsymbol{\sigma} \rVert$ 가 1을 초과하면 그림자 MRP로 전환한다.

8.3 전환 규칙

$\lVert \boldsymbol{\sigma} \rVert \leq 1$ : 원래 MRP 사용 ( $\phi \leq \pi$ )
$\lVert \boldsymbol{\sigma} \rVert > 1$ : 그림자 MRP 사용 ( $\phi > \pi$ )

이 규칙으로 모든 회전을 항상 유한한 매개변수로 표현할 수 있다.

9. MRP의 운동학

MRP의 시간 미분은 각속도와 다음과 같이 연결된다.

$\dot{\boldsymbol{\sigma}} = \frac{1}{4}\left[(1 - \lVert \boldsymbol{\sigma} \rVert^2)\mathbf{I} + 2[\boldsymbol{\sigma}]_\times + 2\boldsymbol{\sigma}\boldsymbol{\sigma}^T\right]\boldsymbol{\omega}$

이 식은 MRP가 각속도로부터 어떻게 시간에 따라 변화하는지 기술하며, 자세 적분에 사용된다.

10. MRP의 합성

두 MRP의 합성은 단순한 벡터 합산이 아니라 다음의 식을 따른다.

$\boldsymbol{\sigma} = \frac{(1 - \lVert \boldsymbol{\sigma}_2 \rVert^2)\boldsymbol{\sigma}_1 + (1 - \lVert \boldsymbol{\sigma}_1 \rVert^2)\boldsymbol{\sigma}_2 - 2\boldsymbol{\sigma}_1 \times \boldsymbol{\sigma}_2}{1 + \lVert \boldsymbol{\sigma}_1 \rVert^2\lVert \boldsymbol{\sigma}_2 \rVert^2 - 2\boldsymbol{\sigma}_1^T\boldsymbol{\sigma}_2}$

이 합성 공식은 회전 행렬의 곱에 대응한다.

11. MRP의 장점

11.1 매개변수의 최소성

MRP는 정확히 3개의 매개변수로 회전을 매개화한다. 쿼터니언(4)이나 회전 행렬(9)보다 매개변수가 적다.

11.2 자유 매개변수

MRP의 매개변수에는 정규화 제약이 없다. 임의의 3차원 벡터가 유효한 MRP이다(다만 $\lVert \boldsymbol{\sigma} \rVert > 1$ 이면 그림자 MRP 사용 권장).

11.3 닫힌 형태 변환

MRP에서 회전 행렬로의 변환이 닫힌 형태이며, 삼각 함수 계산이 필요 없다.

11.4 짐벌 락 부재

오일러 각과 달리 MRP는 짐벌 락이 없다.

12. MRP의 단점

12.1 합성의 복잡성

MRP의 합성은 단순한 행렬 곱이 아니라 위의 복잡한 식을 따른다. 회전 행렬이나 쿼터니언보다 합성이 어렵다.

12.2 그림자 MRP 전환의 처리

큰 회전을 다룰 때 그림자 MRP로의 전환이 필요하며, 전환 시점에서 매개변수가 비연속적으로 변할 수 있다.

12.3 보간의 복잡성

MRP를 직접 보간하면 일정한 각속도를 보장하지 않는다. SLERP와 같은 부드러운 보간을 위해서는 쿼터니언으로 변환 후 보간하는 것이 일반적이다.

13. MRP의 응용

13.1 우주선 자세 제어

NASA와 ESA의 우주선 자세 제어 시스템에서 MRP가 표준적으로 사용된다. 짐벌 락이 없고 매개변수가 적다는 장점이 우주선 자세 추정에 적합하다.

13.2 자세 추정 칼만 필터

확장 칼만 필터의 자세 상태로 MRP를 사용하는 변형(MRP-EKF)이 존재한다. 쿼터니언 기반 EKF와 달리 정규화 제약을 처리할 필요가 없다.

13.3 매니퓰레이터 제어

일부 매니퓰레이터 제어 알고리즘이 MRP를 자세 표현으로 사용한다. 자세 오차의 PID 제어에 자연스럽게 적용된다.

13.4 비선형 최적화

MRP의 자유 매개변수 성질로 인해, 비선형 최적화의 변수로 사용될 수 있다. 단, 큰 회전에서는 그림자 MRP 처리가 필요하다.

14. 다른 매개변수 표현과의 비교

표현법	매개변수 수	제약	특이점	합성
회전 행렬	9	6	없음	행렬 곱
쿼터니언	4	1	$\pm$ 이중성	쿼터니언 곱
오일러 각	3	0	짐벌 락	어색함
회전 벡터	3	0	경미	BCH 공식
CRP	3	0	$\phi = \pi$ 발산	복잡
MRP	3	0	$\phi = 2\pi$ (그림자로 해결)	복잡

15. 참고 문헌

Shuster, M. D. (1993). “A Survey of Attitude Representations.” Journal of the Astronautical Sciences, 41(4), 439–517.
Marandi, S. R., & Modi, V. J. (1987). “A Preferred Coordinate System and the Associated Orientation Representation in Attitude Dynamics.” Acta Astronautica, 15(11), 833–843.
Schaub, H., & Junkins, J. L. (2014). Analytical Mechanics of Space Systems (3rd ed.). AIAA.
Crassidis, J. L., & Junkins, J. L. (2011). Optimal Estimation of Dynamic Systems (2nd ed.). CRC Press.
Tsiotras, P. (1996). “Stabilization and Optimality Results for the Attitude Control Problem.” Journal of Guidance, Control, and Dynamics, 19(4), 772–779.

version: 1.0