2025-11-16T06:37:11.852603

On the equivalence of model-based and data-driven approaches to the design of unknown-input observers

DisarÃ², Valcher

In this paper we investigate a data-driven approach to the design of an unknown-input observer (UIO). Specifically, we provide necessary and sufficient conditions for the existence of an unknown-input observer for a discrete-time linear time-invariant (LTI) system, designed based only on some available data, obtained on a finite time window. We also prove that, under weak assumptions on the collected data, the solvability conditions derived by means of the data-driven approach are in fact equivalent to those obtained through the model-based one. In other words, the data-driven conditions do not impose further constraints with respect to the classic model-based ones, expressed in terms of the original system matrices.

academic

未知入力オブザーバ設計におけるモデルベース法とデータ駆動法の等価性について

基本情報

論文ID: 2311.00673
タイトル: On the equivalence of model-based and data-driven approaches to the design of unknown-input observers
著者: Giorgia Disarò, Maria Elena Valcher (パドヴァ大学、イタリア)
分類: math.DS (力学系)
発表時期: 2023年11月 (arXiv プレプリント)
論文リンク: https://arxiv.org/abs/2311.00673

要約

本論文は、未知入力オブザーバ(UIO)設計のデータ駆動法を研究している。具体的には、著者は離散時間線形時不変(LTI)システムに対して、有限時間窓データに基づくUIO設計の存在性に関する必要十分条件を提供している。同時に、収集データに対する弱い仮定の下で、データ駆動法から導出される可解条件が、実はモデルベース法から得られた条件と等価であることを証明している。言い換えれば、原システム行列で表現された古典的なモデルベース条件に対して、データ駆動条件は追加的な制約を課さない。

研究背景と動機

問題記述: 多くの制御工学応用において、システムの内部状態を把握することは、状態フィードバック安定化と故障検出などの基本的な問題を解決するために重要である。しかし、システム状態は通常直接取得できず、状態ベクトルを推定するための適切なオブザーバを設計する必要がある。
課題: 実際のシステム動態は、しばしば干渉、測定誤差、または状態進化の識別に利用できない他の未知信号の影響を受ける。従来のオブザーバ設計はシステムモデルが既知であると仮定しているが、場合によっては、利用可能なデータのみではシステムを一意に識別できない。
既存方法の限界:
- モデルベース法はシステム行列が既知であることを必要とする
- 既存のデータ駆動法(例:18)は十分条件のみを提供し、得られたシステムが漸近安定でない場合、UIoを設計できるかどうかが不明確である
研究動機: 利用可能な大量のデータを活用して、制御工学問題を解決するためのデータ駆動技術を開発すること。特に、正確なシステムモデルが得られない場合の状態推定問題に対応すること。

核心的貢献

必要十分条件の提供: データに基づくUIO存在性を検証するための検証可能な必要十分条件を提供
完全なパラメータ化: すべての候補UIoの完全なパラメータ化表現を提供
等価性の証明: 特定のデータ仮定の下で、データ駆動法の可解条件とモデルベース法の条件が完全に等価であることを証明
理論的ブレークスルー: データ駆動法とモデルベース法の間に全単射対応関係を確立

方法の詳細

タスク定義

離散時間LTIシステムΣを考える：

x(t+1) = Ax(t) + Bu(t) + Ed(t)  (1)
y(t) = Cx(t)                    (2)

ここで：

x(t) ∈ ℝⁿ: 状態ベクトル
u(t) ∈ ℝᵐ: 既知制御入力
y(t) ∈ ℝᵖ: 出力ベクトル
d(t) ∈ ℝʳ: 未知入力(例:干渉)

UIO定義: 以下の形式のLTIシステムΣ̂：

z(t+1) = A_UIO z(t) + B^u_UIO u(t) + B^y_UIO y(t)  (3)
x̂(t) = z(t) + D_UIO y(t)                           (4)

推定誤差e(t) = x(t) - x̂(t)が初期条件と未知入力動態に依存せず漸近的にゼロに収束するようなもの。

モデルベース法

必要十分条件(定理2)：UIoが存在するための必要十分条件は以下のいずれかが成立することである：

行列(A_UIO, B^u_UIO, B^y_UIO, D_UIO)が存在して以下を満たす：
- A_UIoはSchur安定
- D_UIO CE = E
- B^u_UIO = (I - D_UIO C)B
- A_UIO(I - D_UIO C) + B^y_UIO C = (I - D_UIO C)A
行列条件：
- rank(CE) = rank(E) = r
- rankzIₙ - A -E; C 0 = n + r, ∀z ∈ ℂ, |z| ≥ 1
三つ組(A,E,C)は強*検出可能である

データ駆動法

データ組織: 履歴データを行列形式に組織化：

Uₚ = uₐ(0) ... uₐ(T-2) ∈ ℝᵐˣ⁽ᵀ⁻¹⁾
Xₚ = xₐ(0) ... xₐ(T-2) ∈ ℝⁿˣ⁽ᵀ⁻¹⁾
Xf = xₐ(1) ... xₐ(T-1) ∈ ℝⁿˣ⁽ᵀ⁻¹⁾
Yₚ = yₐ(0) ... yₐ(T-2) ∈ ℝᵖˣ⁽ᵀ⁻¹⁾
Yf = yₐ(1) ... yₐ(T-1) ∈ ℝᵖˣ⁽ᵀ⁻¹⁾

核心仮定: 行列Uₚᵀ Dₚᵀ Xₚᵀᵀが完全行ランクm+r+nを有する。

主要結果(定理9)：以下の条件は等価である：

(3)-(4)の形式のUIoが存在してT_Σ ⊆ T_Σ̂
行列T₁,T₂,T₃,T₄が存在してXf = [T₁ T₂ T₃ T₄]Uₚ; Yₚ; Yf; XₚかつT₄がSchur安定
データ条件：ker(Xf) ⊇ ker(Uₚ; Yₚ; Yf; Xₚ)かつrankzXₚ-Xf; Uₚ; Yₚ = n+m+r, ∀z∈ℂ,|z|≥1
三つ組(A,E,C)は強*検出可能である

技術的革新点

等価性の確立: データ駆動条件とモデルベース条件の完全な等価性を初めて厳密に証明
全単射対応: UIoパラメータ行列とデータ駆動パラメータ間の全単射関係を確立(系10)
実用的アルゴリズム: 問題をより扱いやすい形式に変換する簡略化された計算方法を提供

実験設定

数値例

文献4の例を使用：

A = [-1 -1  0; -1  0  0; 0 -1 -1]
C = [1  0  0; 0  0  1]
E = [-1; 0; 0]

実験パラメータ

時間窓：T = 20
未知入力：(-2,2)区間内で均一にランダム生成
初期条件：ランダムに選択
干渉：(-10,10)区間内でランダム生成

モデル検証: rank(CE) = rank(E) = rと強*検出可能性条件を検証
データ駆動設計: データから行列Cを正常に再構成し、定理9の条件を検証
性能比較:
- 本論文の方法で得られたA_UIO行列は冪零行列
- 文献18の方法で得られた固有値は{0, 0, -0.2}
- 両方法とも機能するが、本論文の方法は収束速度の選択を許可