[論文レビュー] A boundary value problem of a generalised linear discrete time system with no solutions and infinitely many solutions

Nicholas Apostolopoulos, Fernando Ortega+1 more|arXiv (Cornell University)|Oct 26, 2016

Differential Equations and Numerical Methods参考文献 25被引用数 3

ひとこと要約

本稿では、特異な先行係数行列を有する一般化線形離散時間システムにおける非一貫境界値問題（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.

研究の動機と目的

• 先行係数行列が特異で、解が存在しないか無限に多くの解が存在する可能性がある一般化線形離散時間システムにおける境界値問題（BVP）を扱う。
• 従来の解法が失敗する非一貫BVPにおいて、最適解を体系的に計算する手法を提供する。
• 不適切なBVPに対する正則化技術を組み込んだことで、一般化線形システムの既存理論を拡張する。
• さまざまな一貫性条件下での最適解を明示的かつ計算可能にテスト可能な式として導出する。
• ランク不足または一貫性の欠如がある系において、理論的枠組みの妥当性を数値例により検証する。

提案手法

• 行列ペンシル理論を用いて、システム $ sF - G $ の正則性を分析し、非特異変換 $ P $ および $ Q $ を用いてワイエルシュトラス標準形に分解する。
• ワイエルシュトラス形を用いて、システムの一般解を $ 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) $ の場合、モーア・ペントンス逆行列に類似した式 $ (K^*K)^{-1}K^*L $ を用い、残差 $ \|L - K\hat{C}\|_2^2 $ を最小化する。
• ランク不足の状況における逆行列の安定化のため、 $ \|E\|_2 = \theta \ll 1 $ を満たす小さな正則化パラメータ $ \theta $ を行列 $ E $ を通じて導入する。

実験結果

リサーチクエスチョン

• RQ1 特異な先行係数行列を有する一般化線形離散時間システムにおいて、どのような条件下で解が存在しないか、または無限に多くの解が存在するのか。
• RQ2 境界値問題が非一貫である場合、最適解を体系的に導出する方法は何か。
• 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 $ となり、残差の最小化が正確に達成された。
• 本手法は、解が存在しない場合と無限に多くの解が存在する場合の両方を効果的に処理でき、一元的な正則化フレームワークを提供する。
• 理論的結果は数値的に検証され、元のシステムが不適切な場合であっても、導出された式が安定的かつ最適な解をもたらすことが確認された。