[논문 리뷰] Polynomial Time Nondimensionalisation of Ordinary Differential Equations via their Lie Point Symmetries

Évelyne Hubert, Alexandre Sedoglavic|arXiv (Cornell University)|2006. 04. 13.

Nonlinear Waves and Solitons참고 문헌 11인용 수 5

한 줄 요약

이 논문은 Lie 점 대칭을 사용하여 상미분방정식(OED)의 비차원화를 다항시간 알고리즘으로 제시한다. 특히 비율과 이동 대칭을 통해 매개변수의 수를 줄인다. 이러한 대칭으로부터 유리형 불변량을 계산함으로써, 매개변수 수가 줄어든 차원 없는 좌표계로 ODE 시스템을 재작성한다. 이 방법의 계산 복잡도는 ℓ(매개변수의 수)에 대해 O(ℓ³)로 유 bounds된다.

ABSTRACT

Lie group theory states that knowledge of a $m$-parameters solvable group of symmetries of a system of ordinary differential equations allows to reduce by $m$ the number of equation. We apply this principle by finding dilatations and translations that are Lie point symmetries of considered ordinary differential system. By rewriting original problem in an invariant coordinates set for these symmetries, one can reduce the involved number of parameters. This process is classically call nondimensionalisation in dimensional analysis. We present an algorithm based on this standpoint and show that its arithmetic complexity is polynomial in input's size.

연구 동기 및 목표

• 물리적 및 생물학적 모델에서 유도된 ODE 시스템의 매개변수 수를 체계적으로 줄이는 알고리즘을 개발하기 위해.
• 비차원화 과정을 Lie 점 대칭을 기반으로 한 변환으로 체계화하기 위해.
• 유리형 불변량 계산과 시스템 재작성의 계산 복잡도가 입력 크기의 다항시간 내에 보장됨을 증명하기 위해.
• 수동적인 수작업 조작을 피하는 체계적이고 알고리즘 기반의 매개변수 축소 프레임워크를 제공하기 위해.
• 차원 분석을 Lie 군 이론과 불변량 이론을 활용한 계산 기반 대수적 프레임워크로 확장하기 위해.

제안 방법

• ODE 시스템을 변화시키지 않는 일차 파라미터 그룹인 비율과 이동 대칭을 식별하기 위해.
• 편미분방정식의 결정 시스템을 사용하여 이러한 대칭 그룹의 무한소 생성자(infinitesimal generators)를 구성하기 위해.
• 대칭 작용의 계수 행렬에 대해 선형대수학, 특히 가우스 소거법을 활용하여 대칭 그룹의 유리형 불변량을 계산하기 위해.
• 기존 ODE 시스템을 이러한 유리형 불변량을 기반으로 재작성함으로써, 매개변수 수가 줄어든 차원 없는 시스템을 얻기 위해.
• 대칭에 의해 생성된 리 대수의 가비타이를 활용하여 알고리즘에서 이동 대칭을 비율 대칭보다 우선순위를 높이기 위해.
• 대칭 지수에 대한 선형대수 기법을 활용하여 알고리즘의 산술 복잡도가 다항시간에 유지되도록 보장하기 위해, O(ℓ³)로 유 bounds된다.

실험 결과

연구 질문

• RQ1 ODE 시스템에 대한 비차원화 과정을 다항시간 복잡도를 보장하는 대칭 기반 알고리즘으로 체계화할 수 있는가?
• RQ2 Lie 점 대칭에서 유도된 유리형 불변량을 체계적으로 계산하여 ODE의 매개변수 수를 줄일 수 있는가?
• RQ3 이러한 불변량을 계산하고 시스템을 불변량 좌표계로 재작성하는 데 드는 계산 복잡도는 무엇인가?
• RQ4 다양한 유형의 대칭(이동 대칭 대비 비율 대칭)을 적용할 때 매개변수 축소를 극대화하기 위해 어떤 순서로 적용해야 하는가?
• RQ5 이 방법을 비율과 이동 대칭을 초월한 더 일반적인 대칭 유형, 예를 들어 비선형 대칭에도 확장할 수 있는가?

주요 결과

• 알고리즘은 가우스 소거법을 사용하여 Lie 점 대칭에서 유리형 불변량을 계산함으로써 O(ℓ³)로 유 bounds되는 다항시간 복잡도를 보장한다.
• 대칭 그룹을 활용하여 ODE 시스템의 매개변수 수를 성공적으로 줄였으며, 이는 베르훌스트의 로지스틱 모델과 포식자-피식자 시스템에서 검증되었다.
• 선형 포식이 있는 베르훌스트 모델의 경우, 유도된 불변량(t̃ = (a−c)t, x̃ = b/(a−c) x)을 통해 세 매개변수에서 하나의 매개변수로 축소되었으며, 차원 없는 형태로 dx̃/dt̃ = x̃(1−x̃)를 얻었다.
• 포식자-피식자 모델에서는 스케일링 변환을 통해 세 매개변수(r, k₁, k₂)를 제거하여 나머지 세 매개변수(s, e, h)만 남는 시스템을 얻었다.
• 이 방법은 대칭 기저의 선택에 대해 강건하며, 가우스 소거법 중 변수 순서를 적절히 조정함으로써 사용자가 어떤 매개변수를 제거할지 제어할 수 있다.
• 비표준 대칭, 예를 들어 피츠휴-나구모 모델에서의 한 파라미터 군 ∂/∂y + b∂/∂a + ∂/∂d에 대해서도 프레임워크가 확장되며, 불변량 ỹ = y−d 및 â = a+bd를 통해 매개변수 축소가 가능하다.