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

Approximate Tracking Controllability of Systems with Quadratic Nonlinearities

基本信息

  • 论文ID: 2511.12634
  • 标题: Approximate Tracking Controllability of Systems with Quadratic Nonlinearities
  • 作者: Manuel Rissel (ShanghaiTech University), Marius Tucsnak (Université de Bordeaux)
  • 分类: math.OC (Optimization and Control)
  • 发表时间: 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. 研究动机

  • 线性系统的跟踪可控性非常受限（本文证明了这一点），激发了对非线性系统的研究
  • 二次非线性在许多数学模型中自然出现（如Lorenz系统、流体力学方程的有限维近似）
  • 二阶展开可以提供比单纯线性化更精确的非线性系统近似

核心贡献

本文的主要贡献包括：

  1. 线性系统的否定结果：系统地证明了线性系统的弱近似跟踪可控性仅在控制算子B为满射（即控制空间与状态空间同构）这一平凡情形下才成立（Proposition 1.5）
  2. 二次系统的肯定结果：对于满足特定饱和假设（Assumption 3.3）的二次非线性系统，证明了在任意时间范围内的弱近似跟踪可控性（Theorem 3.5）
  3. 松弛范数的关键作用：引入松弛范数|||·|||τ作为度量工具，使得在该弱拓扑下可以获得正面结果，同时通过反例说明在更强的L²范数下无法获得类似结果
  4. 耦合系统的应用：将主要结果应用于耦合ODE系统的控制问题，包括动态运动规划问题（Theorem 7.1）
  5. 全局性结果：即使无控制动力学可能在有限时间内爆破，该框架仍能保证全局可控性

方法详解

任务定义

定义1.1 (近似跟踪可控性)：给定τ>0和L²(0,τ;X)上的范数Nτ，系统(1.1)在0,τ上关于Nτ是近似跟踪可控的，如果对于每个ψ∈W^{1,2}((0,τ);X)和ε>0，存在控制函数u∈L²(0,τ;U)使得初值x(0)=ψ(0)的解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²范数弱，但正是这种弱拓扑使得二次系统的跟踪可控性成为可能。

线性系统的分析（Section 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. Proposition 2.5的证明：
    • 若B非满射，存在η∈X{0}使得B*η=0
    • 构造特殊函数g(t)=-α(t)A*η-α̇(t)η，其中α(t)=(τ-t)η
    • 证明Ψτg=0但g≠0，故Fτ的像集不稠密
  4. 对于弱近似跟踪可控性，考虑扩展系统并得到类似结论（Proposition 1.5）

二次系统的主要结果（Sections 3-6）

饱和假设（Assumption 3.3）：

  • 定义递归序列：E₀=Range(B)，Eᵢ=ℱ(Eᵢ₋₁)
  • 其中ℱ(E)是最大子空间，使得任意γ∈ℱ(E)可表示为γ=ξ₀-∑f(ξᵢ)，ξᵢ∈E
  • 假设存在n_X使得E_=X

简化假设（Assumption 1.7）：对每个γ∈X，存在u,ξ∈U使得γ=Bu-f(Bξ)

证明策略（基于放大系统方法）

Step 1: 放大系统（Section 4） 引入带两个输入的系统： 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时退化为原系统。

关键引理（Proposition 4.2）：建立解对数据的连续依赖性，为后续逼近提供稳定性保证。

Step 2: 轨迹逼近（Section 5，核心技术）

Proposition 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. Lemma 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并利用周期性消去实现的

Step 3: 迭代构造（Section 6）

  1. 对于目标轨迹ψ，定义γ=ψ̇+Aψ+f(ψ)
  2. 应用Proposition 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
    • 每次迭代降低一层，最终达到可控子空间
    • 避免了对全局可控性的直接假设
  4. 稳定性分析：
    • Proposition 4.2保证了扰动下的连续依赖性
    • 允许用光滑控制逼近非光滑控制
    • 为迭代过程提供收敛性保证

实验设置

理论性质的验证实例

本文作为纯理论工作，没有数值实验，但提供了多个说明性例子：

Example 3.1 (Lorenz系统)： 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

验证了Assumption 1.7成立，因此该系统是弱近似跟踪可控的。

Example 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

  • 满足Proposition 1.8的假设，故弱近似跟踪可控
  • 但不是（强）近似跟踪可控

证明思路：选择目标轨迹ψ=ψ₁,0,0^T，其中ψ₁在1/3,2/3上恒等于1000。假设存在L²逼近序列，则x₁在该区间上≥999，从而ẋ₂≥998000，这与L²范数有界矛盾。

Example 3.4 (Assumption 3.3的非平凡性)：

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

可验证Assumption 3.3成立但Assumption 1.7不成立，说明前者确实更一般。

实验结果

主要理论结果

Theorem 3.5 (主要定理)： 设Assumption 3.3成立，τ>0。对任意ψ∈W^{1,2}((0,τ);X)和ε>0，存在控制u∈C^∞(0,τ;U)使得初值x₀=ψ(0)的解x满足： ∣x(τ)−ψ(τ)∣+∣∣∣x−ψ∣∣∣τ<ε|x(\tau)-\psi(\tau)| + |||x-\psi|||_\tau < \varepsilon

推论：

  1. 全局近似可控性（特殊情况）
  2. 即使无控制系统可能有限时间爆破，控制系统仍全局存在

耦合系统的应用（Section 7）

Theorem 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̄)满足第二个方程，且f,B满足Assumption 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. 动态控制：第一个方程作为第二个方程的动态控制器
  2. 运动规划（Example 7.3）：通过控制x来实现z的轨迹跟踪

Example 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²范数：即使对二次系统也可能失败（Example 1.9）
    • 松弛范数：允许利用高频振荡实现控制
  3. 饱和假设的充分性：
    • Assumption 3.3提供了可验证的代数条件
    • 通过有限步迭代即可判断
  4. 全局性vs局部性：
    • 传统小时间局部可控性需要复杂的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. 线性系统的局限性：弱近似跟踪可控性⟺控制算子满射（Proposition 1.5）
  2. 二次系统的可能性：在饱和假设下，任意时间范围内弱近似跟踪可控（Theorem 3.5）
  3. 范数的重要性：
    • 松弛范数：可实现
    • L²范数：一般不可实现（Example 1.9）
    • L^p范数(p≥1)：同样存在障碍（Remark 1.10）
  4. 应用价值：可应用于动态控制和运动规划（Theorem 7.1）

局限性

  1. 二次性质的本质性：
    • 方法高度依赖二次结构
    • 非二次扰动仅在"小"的假设下可处理（Section 1末尾）
    • 扩展到高阶多项式非平凡
  2. 弱拓扑的必要性：
    • 在更强的范数（如L²）下结果不成立
    • 物理应用中可能需要更强的控制
  3. 饱和假设的验证：
    • Assumption 3.3需要计算递归序列
    • 对高维系统可能计算复杂
  4. 光滑性要求：
    • 目标轨迹需要ψ∈W^{1,2}
    • 对非光滑目标的处理未涉及
  5. 控制成本：
    • 未讨论控制的L²范数或能量
    • 高频振荡可能导致高能量消耗

未来方向

论文暗示的研究方向：

  1. 高阶多项式：扩展到三次或更高阶非线性
  2. 无穷维情况：应用到PDE的Galerkin逼近
  3. 控制约束：考虑有界控制或其他约束
  4. 最优性：研究最小能量或最优控制问题
  5. 鲁棒性：对模型不确定性的鲁棒跟踪

深度评价

优点

  1. 理论深度：
    • 完整的正负两方面结果
    • 严格的数学证明
    • 深刻揭示了线性与非线性的本质差异
  2. 方法创新：
    • 巧妙地将PDE控制技术（Agrachev-Sarychev方法）适配到ODE
    • 松弛范数的引入既有数学意义又有物理直观
    • 分层逼近策略优雅且可操作
  3. 结构清晰：
    • 从简单到复杂：线性→二次
    • 从特殊到一般：Assumption 1.7→Assumption 3.3
    • 理论→应用：主定理→耦合系统
  4. 例子丰富：
    • Lorenz系统（Example 3.1）
    • 反例（Example 1.9）
    • 非平凡假设（Example 3.4）
    • 运动规划（Example 7.3）
  5. 技术细节：
    • Lemma 5.2的证明精巧（利用相对紧性+逐点收敛）
    • Proposition 4.2保证了逼近过程的稳定性
    • 对偶性分析（Section 2）清晰透彻

不足

  1. 实用性考虑：
    • 缺乏数值实验验证理论预测
    • 未讨论控制的实际可实现性（能量、频率等）
    • 高频振荡控制在实际系统中的可行性存疑
  2. 假设的限制性：
    • Assumption 3.3虽然可验证，但对某些系统可能不满足
    • 二次性质的限制较强
    • 未讨论假设的必要性
  3. 结果的强度：
    • 仅得到弱拓扑下的结果
    • Example 1.9表明无法推广到L²范数
    • 对终端状态的控制仅是近似的
  4. 分析的完整性：
    • 未讨论n_X（迭代次数）与系统维数的关系
    • 控制的L²范数估计缺失
    • 对ε的依赖性未明确
  5. 应用的深度：
    • Section 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，得到了深刻的正负两方面结果。论文的主要价值在于理论贡献和方法创新，但实用性和数值验证方面有待加强。对于从事非线性控制理论研究的学者，这是一篇值得深入阅读的文献；对于应用研究者，需要结合具体问题评估其适用性。