Estimation of the input-to-state stability gain functions from finite-dimensional approximations

这是一份关于论文《ESTIMATION OF THE INPUT-TO-STATE STABILITY GAIN FUNCTIONS FROM FINITE-DIMENSIONAL APPROXIMATIONS》（从有限维近似估计输入到状态稳定性增益函数）的详细技术总结。

1. 研究背景与问题陈述 (Problem Statement)

背景：
输入到状态稳定性（Input-to-State Stability, ISS）是控制理论中的核心概念，用于量化外部输入对系统状态的影响。ISS 已在有限维系统中得到广泛研究，但在无限维系统（如偏微分方程 PDE）中，特别是涉及无界控制算子（unbounded control operators）的情况（例如边界控制），ISS 增益函数的数值估计仍是一个开放问题。

核心问题：
目前缺乏一种严格的数学框架，能够将无限维系统（特别是具有无界控制算子的系统）的 ISS 性质与其有限维数值近似（如有限差分或有限元离散化）联系起来。
具体而言，作者旨在解决以下两个问题：

1. 问题 2.18： 对于抽象的无限维控制系统 $\dot{x} = Ax + Bu$（其中 $B$ 是已知但无界的算子），如何通过有限维近似序列 $(A_n, B_n)$ 计算 ISS 增益函数 $\beta$ 和 $\gamma$？
2. 问题 2.19： 对于边界控制系统（形式为 $\dot{x}=Ax, Dx=u$），其中控制算子 $B$ 是通过 Fattorini 技巧隐含定义的（即 $B$ 未显式给出），如何仅通过边界算子 $D$ 的离散化来估计 ISS 增益？

动机：
该研究受到物理信息神经网络（PINNs）后验误差界分析的启发。在 PINNs 中，边界条件的不匹配被视为输入误差，需要 ISS 增益 $\gamma$ 来界定整体预测误差。现有的解析解仅适用于简单几何形状，而数值方法可以扩展到复杂几何形状，但前提是建立严格的收敛性理论。

2. 方法论 (Methodology)

作者利用半群理论（Semigroup Theory）和Trotter-Kato 逼近理论构建了理论框架。

2.1 数学基础

• 解析半群与分数阶幂： 假设算子 $-A$ 是扇形算子（sectorial operator），生成解析半群。利用分数阶幂空间 $X_\alpha$ 和 $X_{-1+\alpha}$ 来处理无界算子 $B$。
• Fattorini 技巧： 将边界控制系统转化为具有无界控制算子的抽象 Cauchy 问题：$B = AD_0 - A^{-1}D_0$，其中 $D_0$ 是边界算子 $D$ 的右逆。
• 逼近序列： 引入 Banach 空间序列 $(X_n, P_n, E_n)$，其中 $P_n$ 是投影算子，$E_n$ 是延拓算子，满足 Trotter-Kato 收敛条件。

2.2 核心策略

1. 针对抽象系统 (Problem 2.18)：

   • 将 ISS 增益分解为两部分：初始状态衰减项 $\beta$ 和输入影响项 $\gamma$。
   • $\beta$ 由半群的指数衰减界 $(M, \omega)$ 决定，可通过有限维半群的收敛性获得。
   • $\gamma$ 的估计依赖于积分 $\int_0^t S(t-s)Bu(s)ds$。利用解析半群性质，将积分核分解为 $S(t-s)(-A)^{1-\alpha} \cdot (-A)^{-1+\alpha}B$。
   • 分别估计算子范数：
     • $(-A)^{1-\alpha}S(t)$ 的衰减界（利用谱估计和半群类型）。
     • $(-A)^{-1+\alpha}B$ 的范数，通过证明有限维近似中的 $(-A_n)^{-1+\alpha}B_n$ 收敛到无限维情形。

2. 针对边界控制系统 (Problem 2.19)：

   • 引入预边界闭包离散化 (Pre-boundary closure discretization) 系统 (C)。这是一个中间步骤，保留了边界算子 $D_n$ 的结构，尚未转化为分布控制形式。
   • 在离散层面应用 Fattorini 技巧，构造离散控制算子 $B_n = A_n D_{n,0} - A_n R_n D_{n,0}$。
   • 证明在特定假设下（如 $D_{n,0}$ 的收敛性、算子的一致有界性），离散控制算子 $E_n B_n$ 在强算子拓扑下收敛到无限维算子 $B$。

3. 主要贡献 (Key Contributions)

1. 建立了 ISS 增益的数值计算框架： 首次为具有无界控制算子的无限维系统提供了从有限维近似计算 ISS 增益的严格理论依据。
2. 提出了预边界闭包离散化方法： 针对边界控制问题，提出了一种新的离散化策略（系统 C），作为连接边界形式（System A）和分布控制形式（System D）的桥梁。这解决了直接离散化边界条件时难以处理无界算子 $B$ 的难题。
3. 去除了光滑性假设： 之前的某些方法（如 [10]）要求输入 $u$ 关于时间可微。本文的方法消除了这一限制，恢复了标准的 ISS 估计形式，适用于 $L^\infty$ 输入。
4. 理论结果的推广： 证明了在满足 Trotter-Kato 条件及额外的算子收敛假设下，有限维近似的 ISS 增益序列收敛到无限维系统的真实增益。

4. 主要结果 (Results)

4.1 理论定理

• 定理 3.7 (抽象系统)： 给出了 ISS 增益 $\hat{\beta}$ 和 $\hat{\gamma}$ 的显式公式。
  • $\hat{\beta}(s, t) = \hat{M} e^{-\hat{\omega}t} s$，其中 $\hat{M}, \hat{\omega}$ 来自有限维半群类型的极限。
  • $\hat{\gamma}(s) = \mu_e \kappa \lim_{n\to\infty} \|(-A_n)^{-1+\alpha} B_n\| s$。
  • 其中 $\kappa$ 是由半群衰减率和 resolvent 界导出的常数。
• 定理 4.1 (边界系统)： 证明了在假设 (B1)-(B7) 下（包括 $D_{n,0}$ 的收敛性和算子的一致有界性），离散控制算子 $E_n B_n$ 强收敛于 $B$。这确保了从边界离散化直接计算 ISS 增益的有效性。

4.2 数值算例：一维热传导方程

• 系统： 具有 Dirichlet 边界控制的一维热方程。
• 离散化： 使用有限差分法。
• 验证：
  • 验证了所有理论假设（如扇形性、Trotter-Kato 收敛、分数阶幂空间的收敛性）。
  • 数值计算了增长界 $\omega$、resolvent 界 $D$ 以及算子范数 $\|(-A_n)^\alpha B_n\|$。
  • 结果： 计算出的 ISS 增益 $\gamma(s) \approx 0.8989 s$。
  • 对比： 该结果与文献 [15] 中通过严格解析方法得到的一侧控制增益 $\gamma(s) = \sqrt{2/3}s \approx 0.8165 s$ 非常接近，验证了该方法的有效性和精度。

5. 意义与展望 (Significance & Future Work)

意义：

• 工程应用价值： 使得在复杂几何形状和实际工程问题中，利用数值模拟（如 FEM/FDM）计算 ISS 增益成为可能，从而为 PINNs 等数据驱动方法的误差分析提供了理论工具。
• 理论桥梁： 在数值分析（离散化）与抽象控制理论（ISS、半群）之间建立了严格的联系，特别是处理了无界算子这一难点。
• 方法论创新： “预边界闭包离散化”为处理边界控制问题提供了一种系统化的数值分析视角，不仅限于 ISS，也可用于其他控制理论问题。

局限与未来方向：

• 算子类型限制： 目前仅适用于生成解析半群的扇形算子（主要是抛物型方程），尚未扩展到双曲型问题（如波动方程）。
• 边界条件： 目前主要针对 Dirichlet 边界，Neumann 或混合边界条件的假设验证仍需进一步研究。
• 非线性与非自治系统： 未来计划将框架扩展至半线性系统和非自治系统。
• Riccati 方程结合： 探索将离散化框架与基于 Riccati 方程的 ISS 增益计算方法相结合。

总结：
该论文通过引入预边界闭包离散化和严格的半群逼近理论，成功解决了无限维边界控制系统 ISS 增益的数值估计难题，为复杂控制系统的稳定性分析和误差界计算提供了强有力的数学工具。