技术摘要：关于 Riccati 方程解的数据稀疏性及其在反馈控制中的应用

1. 问题陈述

本文探讨了求解大规模连续时间代数 Riccati 方程（CARE）的计算挑战，该方程形式为：
$$R(X) := A^\top X + XA - XFX + Q = 0$$
其中 $X \in \mathbb{R}^{n \times n}$ 是未知的稳定解，$F$ 和 $Q$ 是对称半正定矩阵。虽然依赖于 Hamiltonian 矩阵分解的标准方法对于小规模 $n$ 是可行的，但对于大规模系统（$n \gg 1$），这些方法会变得极其昂贵。

本研究的重点在于系数矩阵 $A, F, Q$ 具有**拟可分（quasiseparable）结构（即其非对角块为低秩）或为带状（banded）**结构的场景。作者研究了稳定解 $X$ 是否继承了这些结构特性，从而实现高效的近似与计算。这对于源自偏微分方程（PDE）控制和基于智能体模型（agent-based models）的无限时界最优控制问题尤为重要，因为在这些问题中，标准的低秩近似方法（如 RADI 或有理 Krylov 子空间方法）往往会失效，因为其二次项系数 $F$ 是满秩而非低秩的。

2. 研究方法

本文采用了结合对解结构进行理论分析与开发两种不同数值算法的双重方法。

2.1 理论分析：奇异值的衰减

作者证明了如果 $A, F, Q$ 是拟可分的，则解 $X$ 的非对角块其奇异值会呈现快速衰减。

• 与 Sylvester 方程的联系： 分析将 $X$ 的一般非对角块与一个具有低秩右侧项的 Sylvester 方程的解联系起来。
• Zolotarev 数： 衰减率由 Zolotarev 数 $Z_h(E, F)$ 限制，这类数值在涉及闭环矩阵或谱的有理逼近问题中出现。
• 改进的界限： 作者推导出了新的衰减界限，改进了现有关于具有拟可分系数的 Sylvester 方程的结论（特别是 [41, Theorem 2.7]）。这些界限被扩展到了 CARE 场景，涵盖了 $F$ 为单位矩阵、$F$ 为满秩拟可分以及 $F$ 为一般对称正定矩阵的情况。
• 张量列（Tensor Train, TT）秩： 这些衰减结果被用于推导相关值函数 $V(y) = y^\top X y$ 的 TT 秩更紧凑的上界，改进了以往仅限于低秩 $F$ 的估计。

2.2 算法开发

基于系数的结构特性，本文提出了两种求解器：

算法 2：针对一般拟可分系数的分治法（Divide-and-Conquer）

• 方法： 该方法利用基于层次半分离（HSS）矩阵格式的分治策略。
• 机制： 将 CARE 分解为一个块对角部分（递归求解）和一个低秩修正部分。修正项 $\delta X$ 使用投影到张量化 Krylov 子空间上的扩展 Krylov 子空间法（EKSM）进行近似。
• 复杂度： 在理想假设下，复杂度为 $O(n \log n \cdot r^2)$，其中 $r$ 是 HSS 秩。

算法 3：针对带状系数的截断不精确 Newton-Kleinman 法（tink）

• 方法： 专门针对系数为带状的特定情况。它采用了一种不精确的 Newton-Kleinman (NK) 迭代。
• 机制：
  • NK 迭代中产生的 Lyapunov 方程使用 GMRES 方法进行不精确求解。
  • 在每次迭代中应用阈值机制来截断迭代项的带宽，防止带宽不受控制地增长。
  • 线搜索策略确保了 Riccati 残差的单调递减以及闭环矩阵的稳定性。
• 收敛性： 作者提供了收敛性保证，表明通过适当选择内迭代和截断参数，该序列会收敛到稳定解。

3. 主要结果

3.1 理论发现

• 数值拟可分性： 本文证明了即使在二次项系数 $F$ 为满秩的情况下，具有拟可分系数的 CARE 的稳定解也是数值拟可分的。
• 改进的衰减界限： 推导出的非对角奇异值界限比现有文献更紧凑，特别是在满秩 $F$ 的情况下。
• TT 秩估计： 分析为值函数的 TT 秩提供了更精确的估计，显示其对近似误差 $\epsilon$ 的依赖关系为 $O(\log(\epsilon^{-1}))$，这在渐近意义上比之前的 $O(\log(\epsilon^{-1})^{7/2})$ 界限更紧凑。

3.2 数值实验

• 合成测试：
  • 算法 2： 证明了其相对于矩阵规模 $n$ 呈线性对数缩放，相对于拟可分阶数呈二次方缩放。对于较大的 $n$，它比稠密 MATLAB 求解器 icare 快 65 倍至 700 倍。
  • 算法 3 (tink)： 表明对于带状系数，该方法非常高效。然而，对于病态的 $F$，NK 迭代中的中间带宽可能会增长，从而影响性能。“1-LS”变体（仅在第一次迭代进行线搜索）被证明效率最高。
• PDE 应用（Allen-Cahn 方程）：
  • 这些方法被应用于 1D 和 2D Allen-Cahn 方程的最优控制。
  • 在 1D 中，tink 在保持精度的同时，比使用精确解公式的稠密评估实现了 16 倍的加速。
  • 在 2D 中，结合嵌套析散（nested dissection）排序的分治求解器（算法 4）允许解决大规模问题，而这些问题使用稠密求解器需要数千秒。随着网格分辨率的增加，每个控制步的计算成本保持在可控范围内。
• 智能体模型应用（Cucker-Smale）：
  • 使用 d&c_care 算法来控制相互作用的粒子系统。
  • 在求解前对智能体位置进行排序显著降低了解的 HSS 秩，从而降低了 CPU 时间。该方法成功地将系统稳定到了共识状态（consensus state）。

4. 重要性与主张

本文声称提供了首个理论结果，证明了在满秩二次项系数 $F$ 的情况下，CARE 解具有数值拟可分性。这一理论洞察为在比以往理解更广泛的控制问题中使用层次矩阵格式提供了依据。

作者强调，他们提出的算法能够解决目前对于依赖低秩 $F$ 的 RADI 或有理 Krylov 方法而言难以处理的大规模最优控制问题（特别是涉及 PDE 和基于智能体的模型）。这项工作弥合了反馈控制合成中理论衰减性质与实际可扩展算法之间的鸿沟。

5. 局限性与未来方向

作者谦虚地指出了一些有待进一步研究的领域：

• 通过改进分裂策略（如备注 1 和 2 所建议），可以进一步优化衰减界限中的参数 $t$。
• 理论界限中对条件数 $\kappa(F)$ 的依赖似乎是证明过程中的产物，因为数值实验表明，即使对于较大的 $\kappa(F)$，衰减行为依然稳健。
• 在带状情形下，tink 算法在收敛前有时会出现带宽大幅增长的瞬态阶段；提出替代的阈值技术以缓解这一问题被视为一个潜在的改进方向。