2025-12-01T01:43:19.722729

Approximate Tracking Controllability of Systems with Quadratic Nonlinearities

Rissel, Tucsnak
Given a finite-dimensional time continuous control system and $\varepsilon>0$, we address the question of the existence of controls that maintain the corresponding state trajectories in the $\varepsilon$-neighborhood of any prescribed path in the state space. We investigate this property, called approximate tracking controllability, for linear and quadratic time invariant systems. Concerning linear systems, our answers are negative: by developing a systematic approach, we demonstrate that approximate tracking controllability of the full state is impossible even in a certain weak sense, except for the trivial situation where the control space is isomorphic to the state space. Motivated by these negative findings for linear systems, we focus on nonlinear dynamics. In particular, we prove weak approximate tracking controllability on any time horizon for a general class of systems with arbitrary linear part and quadratic nonlinear terms. The considered weak notion of approximate tracking controllability involves the relaxation metric. We underline the relevance of this weak setting by developing applications to coupled systems (including motion planning problems) and by remarking obstructions that would arise for natural stronger norms. The exposed framework yields global results even if the uncontrolled dynamics might exhibit singularities in finite time.
academic

二次非線形性を持つシステムの近似追跡可制御性

基本情報

  • 論文ID: 2511.12634
  • タイトル: Approximate Tracking Controllability of Systems with Quadratic Nonlinearities
  • 著者: Manuel Rissel (ShanghaiTech University)、Marius Tucsnak (Université de Bordeaux)
  • 分類: math.OC (最適化と制御)
  • 発表日: 2025年11月16日 (arXiv プレプリント)
  • 論文リンク: https://arxiv.org/abs/2511.12634

概要

本論文は、有限次元時間連続制御システムの近似追跡可制御性問題を研究する。与えられたε>0に対して、状態軌跡が状態空間内の任意の予定経路のε近傍内に留まるような制御が存在するかどうかを調査する。本論文は線形および二次時不変システムについて体系的に研究を行った。線形システムについては否定的な結果が得られた：制御空間と状態空間が同型であるという自明な場合を除き、弱い意味においても全状態の近似追跡可制御性は不可能である。線形システムのこれらの否定的な結果に基づいて、著者は非線形動力学へと転じ、任意の線形部分と二次非線形項を持つ一般的なシステムクラスに対して、任意の時間範囲における弱近似追跡可制御性を証明した。考慮される弱近似追跡可制御性の概念は、緩和度量(relaxation metric)を含む。結合システム（運動計画問題を含む）への応用を開発することで、著者はこの弱い設定の関連性を強調し、より強いノルムの下で現れる障害を指摘した。このフレームワークは、制御なし動力学が有限時間内に特異点を持つ可能性がある場合でも、大域的な結果を生成する。

研究背景と動機

1. 研究問題

従来の可制御性理論は主に点対点制御問題に焦点を当てており、入力信号が状態軌跡を任意の初期状態から予定された最終状態へ駆動できるかどうかを研究している。これとは対照的に、**追跡可制御性(tracking controllability)**は、時間区間0,τ全体にわたってシステムの全状態または出力関数を制御し、予定された軌跡を追跡することに関心を持つ。具体的には、本論文は以下の制御システムを研究する：

x˙(t)+Ax(t)+f(x(t))=Bu(t),t∈[0,τ]\dot{x}(t) + Ax(t) + f(x(t)) = Bu(t), \quad t \in [0,\tau]

ここで、xは状態、uは制御、Aは線形作用素、Bは制御作用素、fは（可能性のある）非線形項である。

2. 問題の重要性

追跡可制御性問題は複数の応用分野で重要な意義を持つ：

  • 複雑ネットワーク制御：ネットワーク状態を期望軌跡の近くに維持する必要性
  • ロボット工学：経路追跡と運動計画の基礎理論
  • 機械学習：ニューラルODEなどの連続深層学習モデルの可制御性分析

3. 既存方法の限界

  • 出力追跡の研究は比較的充実している（García-Planas と Domínguez-García、Zamorano と Zuazua など）が、全状態追跡可制御性は体系的な研究が不足している
  • 非線形システムの追跡可制御性に関する既存研究は主にドリフトのないアフィン制御システムに限定されている
  • 有限次元システム、特に二次非線形を伴うシステムについては、体系的な理論フレームワークが不足している

4. 研究の動機

  • 線形システムの追跡可制御性は非常に制限されている（本論文で証明）ことが、非線形システムの研究を促した
  • 二次非線形は多くの数学モデルに自然に現れる（例：ローレンツシステム、流体力学方程式の有限次元近似）
  • 二次展開は単なる線形化よりも非線形システムのより正確な近似を提供できる

核心的貢献

本論文の主な貢献は以下の通りである：

  1. 線形システムの否定的結果：線形システムの弱近似追跡可制御性は、制御作用素Bが全射である（すなわち、制御空間と状態空間が同型である）という自明な場合にのみ成立することを体系的に証明した（命題1.5）
  2. 二次システムの肯定的結果：特定の飽和仮定（仮定3.3）を満たす二次非線形システムに対して、任意の時間範囲における弱近似追跡可制御性を証明した（定理3.5）
  3. 緩和ノルムの重要な役割：緩和ノルム|||·|||τを測定ツールとして導入し、この弱位相の下で肯定的な結果を得ることが可能になり、同時に反例を通じてより強いL²ノルムの下では同様の結果が得られないことを示した
  4. 結合システムへの応用：主要な結果を結合ODE システムの制御問題に適用し、動的運動計画問題を含む（定理7.1）
  5. 大域的結果：制御なし動力学が有限時間内に爆発する可能性がある場合でも、このフレームワークは大域的可制御性を保証する

方法の詳細説明

タスク定義

定義1.1 (近似追跡可制御性)：τ>0とL²(0,τ;X)上のノルムNτが与えられたとき、システム(1.1)が0,τ上でNτに関して近似追跡可制御であるとは、すべてのψ∈W^{1,2}((0,τ);X)とε>0に対して、初期値x(0)=ψ(0)を持つ制御関数u∈L²(0,τ;U)と解xが存在して以下を満たすことである：

Nτ(x−ψ)≤εN_\tau(x-\psi) \leq \varepsilon

緩和ノルム（重要なツール）：v∈L¹(0,τ;E)に対して、以下のように定義される：

∣∣∣v∣∣∣τ:=sup⁡t∈[0,τ]∥∫0tv(s)ds∥E|||v|||_\tau := \sup_{t\in[0,\tau]} \left\|\int_0^t v(s)ds\right\|_E

このノルムはL²ノルムより弱いが、この弱い位相こそが二次システムの追跡可制御性を可能にする。

線形システムの分析（第2節）

中心的な考え方：双対性分析を通じて、追跡可制御性問題を作用素の零空間の特性付けに変換する。

主要なステップ：

  1. 作用素Fτ：L²(0,τ;U)→L²(0,τ;Y)を導入する： (Fτu)(t)=C∫0texp⁡((t−σ)A)Bu(σ)dσ(F_\tau u)(t) = C\int_0^t \exp((t-\sigma)A)Bu(\sigma)d\sigma
  2. その双対作用素Ψτを計算する： (Ψτg)(σ)=B∗[∫τσexp⁡((t−σ)A∗)C∗g(t)dt](\Psi_\tau g)(\sigma) = B^*\left[\int_\tau^\sigma \exp((t-\sigma)A^*)C^*g(t)dt\right]
  3. 命題2.5の証明：
    • Bが全射でない場合、η∈X{0}が存在してB*η=0
    • 特殊な関数g(t)=-α(t)A*η-α̇(t)ηを構成する。ここでα(t)=(τ-t)η
    • Ψτg=0だがg≠0であることを証明し、したがってFτの像集合は稠密でない
  4. 弱近似追跡可制御性については、拡張システムを考慮して同様の結論を得る（命題1.5）

二次システムの主要な結果（第3～6節）

飽和仮定（仮定3.3）：

  • 再帰的数列を定義する：E₀=Range(B)、Eᵢ=ℱ(Eᵢ₋₁)
  • ここでℱ(E)は最大部分空間であり、任意のγ∈ℱ(E)がγ=ξ₀-∑f(ξᵢ)、ξᵢ∈Eと表現できる
  • n_Xが存在してE_=Xであると仮定する

簡略仮定（仮定1.7）：各γ∈Xに対して、u,ξ∈Uが存在してγ=Bu-f(Bξ)

証明戦略（拡大システム方法に基づく）

ステップ1：拡大システム（第4節） 2つの入力を持つシステムを導入する： w˙(t)+A(w(t)+ζ(t))+f(w(t)+ζ(t))=γ(t)\dot{w}(t) + A(w(t)+\zeta(t)) + f(w(t)+\zeta(t)) = \gamma(t)

ここでζは乗法制御、γは加法制御である。ζ=0のとき、元のシステムに退化する。

主要な補題（命題4.2）：データに対する解の連続依存性を確立し、その後の近似に安定性保証を提供する。

ステップ2：軌跡近似（第5節、中核技術）

命題5.1の証明の概要：

  1. γ∈E_ℓ（E_ℓ内の値を取る区分定数関数）に対して、γ=η-f(ξ)、ξ,η∈E_{ℓ-1}と仮定する
  2. 二次性質を利用する(1.4)： f(w)−γ=12(f(w+ξ)+f(w−ξ))−ηf(w)-\gamma = \frac{1}{2}(f(w+\xi)+f(w-\xi))-η
  3. 高周波振動制御を構成する： ζn(t)=Enξ,ここで(Enξ)(t)=χ(nt/τ)ξ\zeta_n(t) = \mathcal{E}_n\xi, \quad \text{ここで} \quad (\mathcal{E}_n\xi)(t) = \chi(nt/\tau)\xi
    χは1周期関数：χ(t)=ξ（[0,1/2)上）、χ(t)=-ξ（[1/2,1)上）
  4. 補題5.2（主要技術補題）：連続φ：X→Xに対して、 hn:=ϕ(w1+ζn)−12(ϕ(w1+ξ)+ϕ(w1−ξ))+Aζnh_n := \phi(w_1+\zeta_n) - \frac{1}{2}(\phi(w_1+\xi)+\phi(w_1-\xi)) + A\zeta_n
    が以下を満たすことを証明する：lim⁡n→∞∥Ihn∥C([0,τ];X)=0\lim_{n\to\infty}\|\mathcal{I}h_n\|_{C([0,\tau];X)}=0
    ここで(Iz)(t)=∫0tz(s)ds(\mathcal{I}z)(t)=\int_0^t z(s)ds
  5. 緩和ノルムの収束を証明する： ∣∣∣ζn∣∣∣τ→0asn→∞|||\zeta_n|||_\tau \to 0 \quad \text{as} \quad n\to\infty
    これは積分∫0tζn(s)ds\int_0^t\zeta_n(s)dsを分析し、周期性を利用して消去することで実現される

ステップ3：反復構成（第6節）

  1. 目標軌跡ψに対して、γ=ψ̇+Aψ+f(ψ)を定義する
  2. 命題5.1を適用してE_ℓからE_{ℓ-1}へ近似する
  3. n_X回反復し、最終的にE₁=Range(B)内の制御を得る
  4. B:U→Range(B)が全射であるため、そのような制御を実現するu∈Uを見つけることができる

技術的革新点

  1. 緩和ノルムの巧妙な使用：
    • L²ノルムより弱く、高周波振動を許容する
    • 物理的な「平均効果」の概念と一致する
    • 積分的意味で誤差を制御する
  2. 高周波振動制御の設計：
    • Agrachev-Sarychev方法を借用（元々流体PDE用）
    • 二次非線形の特殊構造を利用する
    • 高速振動を通じて「凸化原理」を実現する
  3. 階層的近似戦略：
    • 再帰的に定義された部分空間列E_iを通じて
    • 各反復で1段階低下し、最終的に可制御部分空間に到達
    • 大域可制御性の直接的な仮定を回避する
  4. 安定性分析：
    • 命題4.2は摂動下の連続依存性を保証する
    • 滑らかでない制御を滑らかな制御で近似することを許容する
    • 反復プロセスの収束性を保証する

実験設定

理論的性質の検証例

本論文は純粋な理論的研究であり、数値実験はないが、複数の説明的な例を提供している：

例3.1 (ローレンツシステム)： x˙=[σ−σ0−ρ1000β]x+[0x2x3−x1x2]+[100101]u\dot{x} = \begin{bmatrix}\sigma & -\sigma & 0\\-\rho & 1 & 0\\0 & 0 & \beta\end{bmatrix}x + \begin{bmatrix}0\\x_2x_3\\-x_1x_2\end{bmatrix} + \begin{bmatrix}1 & 0\\0 & 1\\0 & 1\end{bmatrix}u

仮定1.7が成立することを検証し、したがってこのシステムは弱近似追跡可制御であることを示した。

例1.9 (否定的結果)： 以下のシステムを考える： x˙1+x1x2=u1,x˙2+x32−x12=0,x˙3−x3x2=u2\dot{x}_1 + x_1x_2 = u_1, \quad \dot{x}_2 + x_3^2 - x_1^2 = 0, \quad \dot{x}_3 - x_3x_2 = u_2

  • 命題1.8の仮定を満たすため、弱近似追跡可制御
  • しかし、（強）近似追跡可制御ではない

証明の概要：目標軌跡ψ=ψ₁,0,0^Tを選択する。ここでψ₁は1/3,2/3上で恒等的に1000である。L²近似列が存在すると仮定すると、x₁はこの区間上で≥999であり、したがってẋ₂≥998000となり、これはL²ノルムの有界性と矛盾する。

例3.4 (仮定3.3の非自明性)：

  • X=ℝ⁶、U=ℝ³
  • f(x)=0,0,0,x₁x₂,x₁x₃,x₂x₃^T
  • B=I₃,O₃^T

仮定3.3は成立するが仮定1.7は成立しないことを検証でき、前者が確かにより一般的であることを示す。

実験結果

主要な理論的結果

定理3.5 (主定理)： 仮定3.3が成立し、τ>0とする。任意のψ∈W^{1,2}((0,τ);X)とε>0に対して、初期値x₀=ψ(0)を持つ解xが以下を満たすような制御u∈C^∞(0,τ;U)が存在する： ∣x(τ)−ψ(τ)∣+∣∣∣x−ψ∣∣∣τ<ε|x(\tau)-\psi(\tau)| + |||x-\psi|||_\tau < \varepsilon

系：

  1. 大域近似可制御性（特殊な場合）
  2. 制御なしシステムが有限時間爆発の可能性がある場合でも、制御システムは大域的に存在する

結合システムへの応用（第7節）

定理7.1：結合システムを考える x˙+Ax+f(x)=Bu,z˙+Γ~(z,x)+B(z)=0\dot{x} + Ax + f(x) = Bu, \quad \dot{z} + \tilde{\Gamma}(z,x) + \mathcal{B}(z) = 0

(z̄,x̄)が第2の方程式を満たし、f,Bが仮定3.3を満たす場合、任意のε>0に対して、以下を満たすuが存在する： sup⁡t∈[0,τ]∣z(t)−zˉ(t)∣+∣x(τ)−xˉ(τ)∣+∣∣∣x−xˉ∣∣∣τ<ε\sup_{t\in[0,\tau]}|z(t)-\bar{z}(t)| + |x(\tau)-\bar{x}(\tau)| + |||x-\bar{x}|||_\tau < \varepsilon

応用シーン：

  1. 動的制御：第1の方程式が第2の方程式の動的コントローラとして機能
  2. 運動計画（例7.3）：xの制御を通じてzの軌跡追跡を実現

例7.3の具体的な設定：

  • X=Z=ℝ^d、Γ̃(a,b)=a₁b₁,...,a_db_d^T
  • 目標：特定の正交象限内に留まる参照曲線z_refを追跡
  • 構成：x_ℓ(t)=(𝒷_ℓ(z̄(t))-ż̄_ℓ(t))/z̄_ℓ(t)
  • 結果：sup_|z(t)-z_ref(t)|<ε

理論的発見

  1. 線形と非線形の根本的な違い：
    • 線形システム：追跡可制御性は極めて制限されている
    • 二次システム：弱位相の下で豊かな可制御性を持つ
  2. ノルム選択の重要性：
    • L²ノルム：二次システムでも失敗する可能性がある（例1.9）
    • 緩和ノルム：高周波振動の利用を許容する
  3. 飽和仮定の十分性：
    • 仮定3.3は検証可能な代数的条件を提供する
    • 有限ステップの反復で判定可能
  4. 大域性対局所性：
    • 従来の小時間局所可制御性は複雑なLie括弧条件を必要とする
    • 本論文の方法は直接大域的結果を与える

関連研究

主要な研究方向

  1. 古典的可制御性理論：
    • Coronの専著：非線形システム可制御性の体系的論述
    • Jurdjevic-Kupkaの拡大方法：本論文の方法論的基礎
  2. 出力追跡可制御性：
    • García-Planas & Domínguez-García 12：線形システムの出力追跡
    • Zamorano & Zuazua 26：有限次元線形システムの追跡可制御性と双対特性付け
    • 本論文：出力追跡から全状態追跡への拡張
  3. アフィンシステムの追跡：
    • Liu & Sussmann 25：ドリフトなしアフィンシステムの近似追跡
    • Haynes & Hermes 13：Lie理論を通じた非線形可制御性
    • 本論文：ドリフト項（線形+二次）を許容
  4. 小時間局所可制御性：
    • Beauchard & Marbach 5：スカラー入力システムの二次障害
    • 本論文：大域全時間可制御性に焦点
  5. 無限次元システム：
    • Agrachev & Sarychev 2,3：Eulerおよび Navier-Stokes方程式の可制御性
    • Shirikyan 24、Nersesyan 20,22：流体方程式の近似可制御性
    • Koike et al. 15：非圧縮流体における緩和強化
    • 本論文：これらの技術を有限次元二次システムに適応

本論文の利点

  1. 体系性：有限次元システムの全状態追跡可制御性を初めて体系的に研究
  2. 完全性：否定的結果（線形）と肯定的結果（二次）の両方を提供
  3. 検証可能性：代数的に検証可能な十分条件を提供
  4. 応用性：結合システムと運動計画への応用を開発

結論と議論

主要な結論

  1. 線形システムの限界：弱近似追跡可制御性⟺制御作用素が全射（命題1.5）
  2. 二次システムの可能性：飽和仮定の下で、任意の時間範囲における弱近似追跡可制御性（定理3.5）
  3. ノルムの重要性：
    • 緩和ノルム：実現可能
    • L²ノルム：一般に実現不可能（例1.9）
    • L^p ノルム(p≥1)：同様の障害が存在（注釈1.10）
  4. 応用価値：動的制御と運動計画に適用可能（定理7.1）

限界

  1. 二次性質の本質性：
    • 方法は二次構造に高度に依存している
    • 非二次摂動は「小さい」仮定の下でのみ処理可能（第1節末尾）
    • 高次多項式への拡張は非自明
  2. 弱位相の必要性：
    • より強いノルム（例えばL²）の下では結果が成立しない
    • 物理的応用ではより強い制御が必要な場合がある
  3. 飽和仮定の検証：
    • 仮定3.3は再帰列の計算を必要とする
    • 高次元システムでは計算が複雑になる可能性がある
  4. 滑らかさの要件：
    • 目標軌跡はψ∈W^{1,2}を必要とする
    • 非滑らかな目標への対応は未検討
  5. 制御コスト：
    • 制御のL²ノルムまたはエネルギーについては未検討
    • 高周波振動は高エネルギー消費につながる可能性がある

今後の方向

論文が示唆する研究方向：

  1. 高次多項式：3次以上の非線形への拡張
  2. 無限次元の場合：PDE のGalerkin近似への応用
  3. 制御制約：有界制御またはその他の制約を考慮
  4. 最適性：最小エネルギーまたは最適制御問題の研究
  5. ロバスト性：モデル不確実性に対するロバスト追跡

深い評価

利点

  1. 理論的深さ：
    • 完全な正負両方の結果
    • 厳密な数学的証明
    • 線形と非線形の本質的な違いを深く明らかにする
  2. 方法的革新：
    • PDE制御技術（Agrachev-Sarychev方法）をODEに巧妙に適応
    • 緩和ノルムの導入は数学的意義と物理的直観の両方を持つ
    • 階層的近似戦略は優雅で操作可能
  3. 構造の明確性：
    • 単純から複雑へ：線形→二次
    • 特殊から一般へ：仮定1.7→仮定3.3
    • 理論→応用：主定理→結合システム
  4. 豊富な例：
    • ローレンツシステム（例3.1）
    • 反例（例1.9）
    • 非自明な仮定（例3.4）
    • 運動計画（例7.3）
  5. 技術的詳細：
    • 補題5.2の証明は精巧（相対コンパクト性+点別収束を利用）
    • 命題4.2は近似プロセスの安定性を保証
    • 双対性分析（第2節）は明確透徹

不足

  1. 実用性の考慮：
    • 理論予測を検証する数値実験の欠如
    • 制御の実際の実現可能性（エネルギー、周波数など）について未検討
    • 実際のシステムにおける高周波振動制御の実行可能性に疑問
  2. 仮定の限定性：
    • 仮定3.3は検証可能だが、特定のシステムでは満たされない可能性がある
    • 二次性質の制限は強い
    • 仮定の必要性については未検討
  3. 結果の強度：
    • 弱位相の下での結果のみ
    • 例1.9はL²ノルムへの一般化が不可能であることを示す
    • 終端状態の制御は近似的のみ
  4. 分析の完全性：
    • n_X（反復回数）とシステム次元の関係について未検討
    • 制御のL²ノルム推定が欠落
    • εへの依存性が明確でない
  5. 応用の深さ：
    • 第7節の応用は比較的単純
    • 実際の運動計画アルゴリズムとの比較がない
    • 計算複雑性分析の欠如

影響力

分野への貢献：

  1. 有限次元システムの全状態追跡可制御性を初めて体系的に研究
  2. 線形システムの不可能性の完全な理論を確立
  3. 二次システムに対して操作可能な十分条件を提供
  4. 有限次元ODE と無限次元PDE の制御理論を連結

実用的価値：

  1. 運動計画に理論的基礎を提供
  2. 動的コントローラ設計の新しい視点
  3. ロボット工学と自動制御への潜在的応用

再現性：

  • 理論的結果は完全に再現可能（純粋数学的証明）
  • 例の検証は直接的（代数計算）
  • 数値実装にはさらなる作業が必要（高周波制御の数値近似）

限界：

  • 弱位相の実際的な意義はさらなる説明が必要
  • 二次性質は直接応用の範囲を制限
  • 理論的結果を支持する数値実験が必要

適用可能なシーン

理論的に適切なシーン：

  1. 二次非線形を持つ機械システム：平方項を伴う減衰など
  2. 化学反応ネットワーク：二次反応動力学
  3. 生物システム：Lotka-Volterra型モデルの制御
  4. 流体システムの低次元モデル：Galerkin切断

慎重に扱うべきシーン：

  1. 強い制御要件：L²またはL^∞追跡が必要な場合は不適切
  2. 高速応答システム：高周波振動が実行不可能な場合
  3. エネルギー制限システム：制御コストが過度に高い可能性
  4. 高次非線形が支配的：二次近似が不十分

潜在的な拡張方向：

  1. 数値最適化と結合し、実用的なアルゴリズムを設計
  2. 有界制御制約を考慮
  3. ロバスト性と感度を研究
  4. 時変または確率的システムへの拡張

参考文献

論文が引用する主要な文献：

  1. 2,3 Agrachev & Sarychev: 2D Euler および Navier-Stokes方程式の可制御性、本論文の方法の主要な着想源
  2. 26 Zamorano & Zuazua: 有限次元線形システムの追跡可制御性、双対性分析フレームワークを提供
  3. 15 Koike et al.: 制御を通じた非圧縮流体の緩和強化、PDE における緩和ノルムの応用を示す
  4. 24 Shirikyan: 3D Navier-Stokes方程式の近似可制御性、凸化原理の体系的説明
  5. 14 Jurdjevic & Kupka: 群作用下の制御システム、拡大方法の理論的基礎
  6. 8 Coron: 非線形可制御性の古典的専著

総合評価：これは高品質な理論制御論文であり、有限次元システムの全状態追跡可制御性の理論フレームワークを体系的に確立している。緩和ノルムの導入と二次非線形の特殊構造の利用を通じて、無限次元PDE制御の技術を有限次元ODE に巧妙に適応させ、深刻な正負両方の結果を得ている。論文の主な価値は理論的貢献と方法的革新にあるが、実用性と数値検証の側面では強化の余地がある。非線形制御理論の研究に従事する学者にとって、これは深く読む価値のある文献である；応用研究者にとっては、具体的な問題に対する適用可能性を評価する必要がある。