10.49 수정 로드리게스 매개변수(MRP)의 정의

1. MRP의 정의

수정 로드리게스 매개변수(Modified Rodrigues Parameters, MRP)는 고전 로드리게스 매개변수의 개선된 형태로, 회전 각의 1/4의 탄젠트를 사용한다.

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

여기서

• \hat{\mathbf{u}}: 회전 축 (단위 벡터)
• \phi: 회전 각

이는 고전 로드리게스 매개변수의 \tan(\phi/2)\hat{\mathbf{u}}와 비교하여 절반의 각도를 사용한다.

2. MRP의 도입 동기

2.1 고전 로드리게스 매개변수의 한계

고전 로드리게스 매개변수(Classical Rodrigues Parameters, CRP)는 \phi = \pi에서 발산하여 180도 회전을 표현할 수 없다.

2.2 MRP의 개선

MRP는 \tan(\phi/4)를 사용하여 특이점 위치를 \phi = 2\pi로 이동시킨다. \phi = 2\pi는 실제로 \phi = 0과 같은 회전이므로, MRP는 거의 완전한 회전 범위 [0, 2\pi)를 매개화할 수 있다.

3. 매개변수의 크기

MRP의 크기는 다음과 같이 회전 각에 의존한다.

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

3.1 구체적 값

• \phi = 0: \lVert\boldsymbol{\sigma}\rVert = 0
• \phi = 90°: \lVert\boldsymbol{\sigma}\rVert = \tan(22.5°) \approx 0.414
• \phi = 180°: \lVert\boldsymbol{\sigma}\rVert = \tan(45°) = 1
• \phi = 270°: \lVert\boldsymbol{\sigma}\rVert = \tan(67.5°) \approx 2.414
• \phi = 360°: \lVert\boldsymbol{\sigma}\rVert \to \infty

CRP의 \phi = 180°에서 \lVert\mathbf{g}\rVert \to \infty와 비교하면, MRP는 180도 회전에서 \lVert\boldsymbol{\sigma}\rVert = 1의 유한 값을 가진다.

4. MRP와 쿼터니언의 관계

쿼터니언 \mathbf{q} = (\cos(\phi/2), \sin(\phi/2)\hat{\mathbf{u}})로부터 MRP는

\boldsymbol{\sigma} = \frac{\mathbf{q}_v}{1 + q_w}

여기서 \mathbf{q}_v는 쿼터니언의 벡터 부분이고 q_w는 스칼라 부분이다.

4.1 유도

쿼터니언의 스칼라 부분이 \cos(\phi/2)이고 벡터 부분의 크기가 \sin(\phi/2)이다. 반각 공식

\tan(\phi/4) = \frac{\sin(\phi/2)}{1 + \cos(\phi/2)}

을 사용하면 위의 관계가 직접 유도된다.

4.2 역방향

MRP로부터 쿼터니언은

\mathbf{q} = \frac{1}{1 + \lVert\boldsymbol{\sigma}\rVert^2}(1 - \lVert\boldsymbol{\sigma}\rVert^2, 2\boldsymbol{\sigma})

이는 정규화를 통해 단위 쿼터니언을 보장한다.

5. 대체 MRP (Shadow MRP)

MRP의 핵심 기능은 대체 매개화(shadow MRP)이다. 매개변수 크기가 커지면 대체 MRP로 전환하여 크기를 제한한다.

5.1 대체 MRP의 정의

\boldsymbol{\sigma}^S = -\frac{\boldsymbol{\sigma}}{\lVert\boldsymbol{\sigma}\rVert^2}

5.2 대체의 특성

• \boldsymbol{\sigma}와 \boldsymbol{\sigma}^S가 같은 회전을 나타낸다.
• \lVert\boldsymbol{\sigma}\rVert\cdot\lVert\boldsymbol{\sigma}^S\rVert = 1
• \lVert\boldsymbol{\sigma}\rVert > 1이면 \lVert\boldsymbol{\sigma}^S\rVert < 1

5.3 전환 규칙

매개변수 크기가 1을 초과하면 대체 MRP로 전환한다.

if norm(sigma) > 1:
    sigma = -sigma / (norm(sigma) ** 2)

이렇게 하면 매개변수 크기가 항상 [0, 1] 범위 내에 있다.

6. MRP의 회전 범위

대체 MRP 메커니즘을 통해 MRP는 완전한 회전 범위를 매개화할 수 있다.

6.1 단위 공 내부

\lVert\boldsymbol{\sigma}\rVert < 1이면 \phi < 180°에 해당한다. 이 영역에서 매개변수가 자연스럽게 사용된다.

6.2 단위 공 바깥

\lVert\boldsymbol{\sigma}\rVert > 1이면 \phi > 180°에 해당한다. 이 경우 대체 MRP로 전환하여 다시 단위 공 내부의 값으로 변환한다.

6.3 실용적 효과

이 메커니즘으로 MRP는 모든 회전을 항상 크기 1 이하의 매개변수로 표현할 수 있다. 이는 수치적으로 매우 안정적이다.

7. MRP의 운동학

MRP의 시간 미분과 각속도의 관계는 다음과 같다.

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

(본체 좌표계 각속도)

이는 CRP보다 복잡하지만 모든 회전 범위에서 잘 정의되어 있다.

8. MRP의 회전 행렬

MRP로부터 회전 행렬을 직접 계산할 수 있다.

\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}

이는 MRP의 대수적 성질에서 유도된 닫힌 형태 공식이다.

9. MRP의 합성

두 MRP의 합성은 다음과 같다.

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

이는 복잡한 공식이지만 닫힌 형태이다.

10. MRP의 장점

10.1 3 매개변수

오일러 각이나 CRP와 같은 최소 매개변수이다.

10.2 거의 완전한 회전 범위

대체 MRP 메커니즘으로 거의 모든 회전을 매개화한다.

10.3 자유 매개변수

단위 노름 같은 제약이 없다.

10.4 크기 제한

대체 MRP 전환으로 매개변수 크기를 [0, 1]로 제한할 수 있다.

10.5 수치적 안정성

CRP보다 훨씬 안정적이다. 큰 회전에서도 잘 작동한다.

11. MRP의 단점

11.1 \phi = 360° 특이점

매개변수 자체의 특이점이 \phi = 0과 같으므로 실용적 문제는 없지만 이론적으로 존재한다.

11.2 대체 전환의 복잡성

대체 MRP 메커니즘이 구현 복잡성을 증가시킨다. 전환 시점의 결정이 중요하다.

11.3 운동학 방정식의 복잡성

쿼터니언의 운동학 방정식보다 복잡하다.

11.4 보간의 어려움

SLERP와 같은 자연스러운 보간 공식이 쿼터니언에 비해 복잡하다.

12. MRP의 응용

12.1 우주선 자세 제어

MRP는 주로 우주선 자세 제어에서 사용된다. 자유 매개변수와 거의 완전한 회전 범위가 우주선 제어 시스템에 적합하다.

12.2 자세 추정

일부 자세 추정 알고리즘이 MRP를 사용한다. 특히 특이점 없는 자세 표현이 필요한 경우이다.

12.3 비선형 최적화

자세를 자유 매개변수로 최적화하는 경우에 MRP가 선택될 수 있다.

12.4 교육적 가치

MRP는 회전 표현의 진화와 한계를 이해하는 데 교육적으로 유용하다.

13. MRP vs. 쿼터니언

쿼터니언이 일반적으로 선호되지만, 우주 항공 특수 응용에서는 MRP가 사용된다.

14. MRP의 수학적 특징

14.1 스테레오그래픽 사영

MRP는 4차원 단위 구면(쿼터니언 공간)에서 3차원 공간으로의 스테레오그래픽 사영(stereographic projection)에 해당한다. 이는 MRP의 대체 매개화가 자연스럽게 나타나는 이유이다.

14.2 S^3에서 \mathbb{R}^3로의 사영

쿼터니언의 “북극” (scalar = 1)에서 방사형으로 3차원 공간으로 사영한다. 북극에 가까운 쿼터니언(작은 회전)은 원점 근처로, 적도에 있는 쿼터니언(180도 회전)은 단위 구면 위로 사영된다.

14.3 대체의 기하학적 의미

쿼터니언 \mathbf{q}와 -\mathbf{q}는 S^3의 대척점이다. 대체 MRP는 이 대척점에서의 사영에 해당한다. 북극과 남극에서 각각 사영하여 두 MRP를 얻는다.

15. MRP의 역사

MRP는 1996년 Tsiotras 등에 의해 도입되었다. 우주 항공 분야에서 CRP의 한계를 극복하기 위한 대안으로 제안되었다. 이후 우주선 자세 제어의 표준 표현 중 하나가 되었다.

16. 결론

수정 로드리게스 매개변수(MRP)는 고전 로드리게스 매개변수의 180도 회전 특이점 문제를 해결한 3 매개변수 회전 표현이다. 대체 MRP 메커니즘을 통해 거의 완전한 회전 범위를 매끄럽게 매개화하며, 자유 매개변수이므로 최적화와 제어에 적합하다. 우주 항공 분야에서 주로 사용되지만, 쿼터니언이 일반적으로 더 많이 채택된다.

17. 참고 문헌

• Tsiotras, P. (1994). “New Control Laws for the Attitude Stabilization of Rigid Bodies.” IFAC Symposium on Automatic Control in Aerospace, 316–321.
• Schaub, H., & Junkins, J. L. (1996). “Stereographic Orientation Parameters for Attitude Dynamics: A Generalization of the Rodrigues Parameters.” Journal of the Astronautical Sciences, 44(1), 1–19.
• Shuster, M. D. (1993). “A Survey of Attitude Representations.” Journal of the Astronautical Sciences, 41(4), 439–517.
• Crassidis, J. L., & Markley, F. L. (2003). “Unscented Filtering for Spacecraft Attitude Estimation.” Journal of Guidance, Control, and Dynamics, 26(4), 536–542.
• Markley, F. L., & Crassidis, J. L. (2014). Fundamentals of Spacecraft Attitude Determination and Control. Springer.

version: 1.0