[논문 리뷰] A boundary value problem of a generalised linear discrete time system with no solutions and infinitely many solutions

Nicholas Apostolopoulos, Fernando Ortega외 1명|arXiv (Cornell University)|2016. 10. 26.

Differential Equations and Numerical Methods참고 문헌 25인용 수 3

한 줄 요약

이 논문은 특이 주도 계수 행렬을 가진 일반화된 선형 이산 시간 시스템에서 비일관성 있는 경계값 문제(BVPs)에 대해 명시적이고 최적의 정규화된 해를 제안한다. 매트릭스 펜슬 이론과 티코노프 유형 정규화를 사용하여, BVP에 해가 없거나 무한히 많은 해가 존재할 경우의 최적 해에 대한 닫힌 형식의 공식을 유도하며, 질량 결함 또는 비일관성 있는 시스템에서도 이론적 접근이 검증된다.

ABSTRACT

In this article we study a class of generalised linear systems of difference equations with given boundary conditions and assume that the boundary value problem is non-consistent, i.e. it has infinite many or no solutions. We take into consideration the case that the coefficients are square constant matrices with the leading coefficient singular and provide optimal solutions. Numerical examples are given to justify our theory.

연구 동기 및 목표

• 주도 계수 행렬이 특이하고 시스템의 해가 존재하지 않거나 무한히 많은 해를 가질 수 있는 일반화된 선형 이산 시간 시스템에서 경계값 문제(BVPs)를 다루는 것.
• 기존의 해법이 실패하는 비일관성 있는 BVP에서 최적의 해를 체계적으로 계산하는 방법을 제공하는 것.
• 기존의 일반화된 선형 시스템 이론을 확장하여, 부정확한 BVP에 대한 정규화 기법을 통합하는 것.
• 다양한 일관성 조건 하에서 최적의 해에 대한 명시적이고 계산 가능하게 검증 가능한 공식을 도출하는 것.
• 질량 결함 또는 비일관성 있는 시스템에서도 이론적 프레임워크의 타당성을 수치 예제로 검증하는 것.

제안 방법

• 비특이 변환 $ P $ 와 $ Q $ 를 통해 Weierstrass 표준형으로 분해하는 매트릭스 펜슬 이론을 사용하여 시스템 $ sF - G $ 의 정규성 분석을 수행한다.
• Weierstrass 형식을 활용하여 시스템의 일반해를 $ Y_k = Q_p J_p^k C $ 로 표현하며, 여기서 $ C $ 는 경계 조건으로부터 결정되는 상수 벡터이다.
• 경계 조건으로부터 행렬 $ K $ 와 벡터 $ L $ 을 정의하여 해의 존재성과 일관성 여부를 판단하는 선형 시스템 $ K C = L $ 을 구성한다.
• 행렬 $ K $ 가 질량 결함이고 $ L \notin \text{col}(K) $ 인 경우, $ (K^*K + E^*E)^{-1}K^*L $ 을 통한 티코노프 정규화를 적용하여 유일하고 안정적인 해를 확보한다.
• 행렬 $ K $ 가 전순위이고 $ L \notin \text{col}(K) $ 인 경우, 잔차 $ \|L - K\hat{C}\|_2^2 $ 를 최소화하기 위해 모어-펜로즈 역행렬 유사 공식 $ (K^*K)^{-1}K^*L $ 을 사용한다.
• 질량 결함인 경우 역행렬의 안정성을 확보하기 위해 행렬 $ E $ 를 통해 작은 정규화 매개변수 $ \theta $ 를 도입하며, $ \|E\|_2 = \theta \ll 1 $ 이다.

실험 결과

연구 질문

• RQ1 주도 계수 행렬이 특이한 일반화된 선형 이산 시간 시스템에서 언제 해가 존재하지 않거나 무한히 많은 해가 존재하는가?
• RQ2 비일관성 있는 BVP에서 최적의 해를 어떻게 체계적으로 유도할 수 있는가?
• RQ3 질량 결함 또는 비일관성 있는 BVP에서 유일하고 안정적이며 계산 가능성이 있는 해를 보장하기 위한 정규화 전략은 무엇인가?
• RQ4 닫힌 형식의 공식으로서 명시적이고 쉽게 검증 가능한 해를 어떻게 표현할 수 있는가?
• RQ5 매트릭스 펜슬 이론은 이러한 일반화된 시스템에서 해의 구조를 어떻게 특성화하는가?

주요 결과

• 비일관성 있는 BVP에서 $ p < r_1 + r_2 $, $ \text{rank}(K) = p $, $ L \notin \text{col}(K) $ 라면 최적의 해는 $ \hat{Y}_k = Q_p J_p^k (K^*K)^{-1}K^*L $ 으로 주어진다.
• 행렬 $ K $ 가 질량 결함이고 $ L \notin \text{col}(K) $ 라면 최적의 해는 $ \hat{Y}_k = Q_p J_p^k (K^*K + E^*E)^{-1}K^*L $ 로 주어지며, $ \|E\|_2 = \theta \ll 1 $ 이여서 역행렬이 존재함을 보장한다.
• 예제 4.1에서 $ \theta = 0.00001 $ 인 경우 최적의 해는 $ \hat{Y}_k = \left[0, 0, 0, \frac{11}{14 \cdot 4^k}, \frac{11}{14 \cdot 4^k}\right]^T $ 로 수렴하며, 정규화에 의한 안정성이 입증된다.
• 예제 4.2에서 $ K $ 가 전순위일 경우 최적의 해는 $ \hat{Y}_k = \left[0.0349 \cdot 2^{-k}, 0, 0, 0.5006 \cdot 4^{-k}, 0.5006 \cdot 4^{-k}\right]^T $ 로 주어지며, 잔차 최소화가 정확하게 이루어짐을 보여준다.
• 이 방법은 해가 존재하지 않는 경우와 해가 무한히 많은 경우를 모두 다루며, 통합된 정규화 프레임워크를 제공한다.
• 이론적 결과는 수치적으로 검증되었으며, 유도된 공식이 원래 시스템이 부정확한 경우에도 안정적이고 최적의 해를 제공함을 확인한다.