[논문 리뷰] Solutions of Higher Order Homogeneous Linear Matrix Differential Equations: Singular Case

Grigoris I. Kalogeropoulos, Charalambos P. Kontzalis|arXiv (Cornell University)|2015. 01. 22.

Matrix Theory and Algorithms참고 문헌 20인용 수 7

한 줄 요약

이 논문은 주어진 계수 행렬이 특이하거나 펜슬의 차원이 비정방형인 특이 케이스에서 고차수 동차 선형 행렬 미분방정식을 해결하기 위해 행렬 펜슬 방법을 개발한다. 해의 존재성과 유일성을 위한 필요 및 충분 조건을 제시하고, 최소 지수를 통해 해 공간을 규명하며, 기본 해 행렬의 열 생성공간 외부에 있는 초기 조건을 가질 경우 해가 유일하지 않을 수 있음을 수치 예제를 통해 설명한다. 명시적인 해 형태를 제시한다.

ABSTRACT

The main objective of this talk is to develop a matrix pencil approach for the study of an initial value problem of a class of singular linear matrix differential equations whose coefficients are constant matrices. By using matrix pencil theory we study the cases of non square matrices and of square matrices with an identically zero matrix pencil. Furthermore we will give necessary and sufficient conditions for existence and uniqueness of solutions and we will see when the uniqueness of solutions is not valid. Moreover we provide a numerical example.

연구 동기 및 목표

• 특이 행렬 펜슬을 가진 고차수 동차 선형 행렬 미분방정식의 초기값 문제를 해결하기 위해.
• 계수 행렬이 비정방형인 경우 또는 행렬식이 항상 0인 경우 발생하는 펜슬의 특이성에 대해 분석하기 위해.
• 특이 케이스에서 해의 존재성과 유일성을 위한 필요 및 충분 조건을 수립하기 위해.
• 해가 유일하지 않을 경우 해 공간의 차원을 최소 지수(c.m.i. 및 r.m.i.)를 사용하여 특성화하기 위해.
• 해를 계산하고 초기 조건의 일致성을 평가하기 위한 수치적 프레임워크를 제공하기 위해.

제안 방법

• 상태 변수를 사용하여 n차 미분방정식을 1차 시스템으로 변환하여 $ FY'(t) = GY(t) $ 형태로 표현한다.
• 행렬 펜슬 이론을 적용하여 $ sF - G $ 펜슬을 분석하며, 특히 $ \det(sF - G) \equiv 0 $ 또는 $ F $ 가 비정방형인 특이 케이스에 중점을 둔다.
• 구조적 불변량 기반으로 유한 및 무한 기본 인수, 최소 지수(c.m.i. 및 r.m.i.)를 고려하여 다섯 개의 하위시스템으로 시스템을 분해한다.
• 비특이 행렬 $ P $ 및 $ Q $ 를 통한 엄격 동치 변환을 적용하여 $ F $ 와 $ G $ 를 크로네커 표준형 $ F_K $ 와 $ G_K $ 으로 간소화한다.
• 해를 $ Y(t) = Q_p e^{J_p(t-t_0)} Z_p(t_0) $ 형태로 구성하며, 여기서 $ Q_p $ 는 해 공간을 생성하고 $ Z_p(t_0) $ 는 일致한 초기 조건으로 결정된다.
• 해의 존재성을 확인하기 위해 $ Y(t_0) \in \text{colspan}(Q_p) $ 를 검증하며, 이 조건을 만족하지 못할 경우 해는 비유일해진다.

실험 결과

연구 질문

• RQ1 특이 행렬 펜슬을 가진 고차수 동차 선형 행렬 미분방정식이 언제 해를 갖는가?
• RQ2 해는 언제 유일하며, 특이 케이스에서 해 공간의 차원은 무엇에 의해 결정되는가?
• RQ3 유한 및 무한 기본 인수, 그리고 최소 지수는 해의 구조에 어떤 영향을 미치는가?
• RQ4 열 최소 지수와 행 최소 지수는 해 공간을 특성화하는 데 어떤 역할을 하는가?
• RQ5 주어진 초기 조건이 시스템과 일치하는지 알고리즘적으로 어떻게 판단할 수 있는가?

주요 결과

• 초기 조건 벡터 $ Y(t_0) $ 가 유한 기본 인수와 관련된 기본 행렬 $ Q_p $ 의 열공간에 포함될 경우에만 해가 존재한다.
• 만약 $ Y(t_0) \notin \text{colspan}(Q_p) $ 라면 시스템은 무한히 많은 해를 가지며, 해 공간의 차원은 자유 매개변수의 수와 같으며, 예시 2에서는 2이다.
• 초기 조건이 기본 해 행렬 $ Q_p $ 와 일치할 때에만 해가 유일하며, 이는 유한 기본 인수와 최소 지수에 의해 결정된다.
• 예시 1에서 펜슬은 유한 기본 인수 $ s-2 $, $ s-1 $ 과 최소 지수 $ \epsilon_1 = 0 $, $ \epsilon_2 = 2 $ 를 가지며, 이로 인해 무한히 많은 해가 존재한다.
• 예시 2에서는 비정방형 $ F $ 와 $ G $ 를 가지며, 해 공간의 차원이 2인 무한 해를 가지며, 상수 $ c_1 $, $ c_2 $ 를 매개변수로 사용한다.
• 해 형태 $ Y(t) = Q_p e^{J_p(t-t_0)} Z_p(t_0) $ 는 $ Y(t_0) \in \text{colspan}(Q_p) $ 일 때에만 유효하며, 이 조건을 만족하지 못할 경우 해는 비유일해지고 $ Y(t) = Q_p e^{J_p t} C $ 형태로 임의의 $ C $ 를 포함하여 표현된다.