타원 궤도 목표‑공격자‑방어자 게임 분석

초록

본 논문은 비조종 목표가 타원 궤도를 따라 움직이는 상황에서, 공격자와 방어자 사이의 3인 게임을 선형‑이차(LQ) 형태로 모델링하고, 행렬 리카티 방정식에 대한 완전한 해석적 해를 도출한다. 이를 통해 Nash 균형 제어법을 분석적으로 얻고, 공격자의 승리 조건을 기하학적으로 정의한 뒤, 타원 궤도 이심률이 게임 결과에 미치는 영향을 조사한다. 해석법은 수치법 대비 99.9 % 이상의 CPU 시간 절감과 0.004 % 수준의 비용 함수 오차를 보인다.

상세 분석

이 논문은 기존의 두 플레이어(추격‑회피) 차동 게임을 세 플레이어(공격‑방어‑목표) 구조로 확장하고, 특히 목표가 비조종(non‑maneuvering)이며 타원 궤도를 따라 움직이는 상황을 다룬다. 이를 위해 저자들은 Tschaüner‑Hempel(TH) 방정식을 사용하여 선형 상대역학 모델을 구축하고, 진위이상(真異常) 변수 f (진위이상) 를 독립 변수로 삼아 상태 전이 행렬(STM) Ω₁₁(f,f₀) 를 명시적으로 구한다. 변환 행렬 ρ=1+e cos f 를 도입해 좌표를 정규화하고, 제어 입력 u 을 포함한 상태 방정식 ˜X′=A˜X+Bu 로 정리한다.

LQ 게임의 비용 함수는 공격자와 방어자의 상태 가중치 Sₐ, S_d 와 제어 가중치 Rₐ, R_d 를 포함하며, 목표는 고정점(˜X_t=0)으로 가정한다. 해밀턴ian을 구성하고, 최적 제어 조건 ∂H/∂uₐ=0, ∂H/∂u_d=0 로부터 Nash 균형 제어법
uₐ* = –Rₐ⁻¹ Bᵀ(λ–ν), u_d* = R_d⁻¹ Bᵀ ν
을 얻는다. 여기서 λ, ν는 각각 공격자와 방어자의 비용ate 변수이며, λ′=–Aᵀλ, ν′=–Aᵀν 로 전이된다.

핵심 난제는 λ와 ν 사이의 관계를 정의하는 행렬 리카티 방정식
P′ = W₂₂ P – P W₁₁ – P W₁₂ P
이다. 기존 연구는 이 방정식을 수치적으로 뒤로 적분했지만, 저자는 TH 방정식의 해석적 전이 행렬과 변수 변환 기법을 이용해 P(f) 를 명시적으로 구한다. 구체적으로, W₁₁, W₁₂, W₂₂ 를 A, B, Rₐ, R_d 로 구성하고, 전이 행렬 Ω₁₁, Ω₂₂ 를 이용해 P(f) = Ω₂₂(f,f₀) P(f₀) Ω₂₂⁻¹(f,f₀) 형태로 전개한다. 초기 전이 조건 P(f_f)=diag(Sₐ, –S_d) 를 적용하면, 전체 게임에 대한 해석적 Nash 전략이 완성된다.

이 해석적 전략을 바탕으로 공격자의 승리 조건을 도출한다. 공격자는 초기 상대 위치 ˜xₐ와 방어자의 초기 위치 ˜x_d 사이의 거리와 방향이 특정 집합 𝔇 내에 있을 때만 목표를 잡을 수 있다. 𝔇는 P(f₀) 와 Ω₁₁, Ω₂₂ 에 의해 정의되는 6차원 타원형(또는 타원체) 영역이며, 기하학적으로는 “방어자 초기 위치가 공격자-목표 선을 둘러싼 원뿔형(또는 타원형) 안에 있어야 한다”는 형태로 해석된다.

또한 저자는 궤도 이심률 e 가 증가할수록 Ω₁₁·Ω₂₂ 의 비선형성도 증가해 𝔇의 크기가 축소됨을 수치·해석적으로 보여준다. 즉, 고이심률 궤도에서는 방어자가 더 좁은 영역에 위치해야만 공격자가 승리할 가능성이 높아진다.

시뮬레이션에서는 동일한 초기 조건에 대해 해석적 해와 수치적 백워드 적분 해를 비교했으며, 비용 함수 차이는 0.004 % 이하, CPU 시간은 해석적 방법이 0.1 % 수준(99.9 % 절감)으로 크게 우수함을 입증한다. 이러한 결과는 실시간 우주 작전, 특히 빠른 결정을 요구하는 위협 회피 시나리오에 직접적인 적용 가능성을 시사한다.