精读笔记
Problem Setting
论文标题：Complete Trajectory Tracking for Polynomial Differential Variational Inequalities: A Moment-SOS Based Structure-Aware Approach（arXiv preprint / 2026-06-09）。
这篇论文实际解决的是：在非凸、多解、多项式 DVI 中，如何完整跟踪由不同 VI equilibrium selection 诱导出的轨迹分支。它不是普通的 DVI existence paper，也不是单纯把 Moment-SOS 套到每个时间步上；真正目标是把非凸 VI 的解集 multiplicity 作为动力系统状态的一部分来处理。
传统 DVI 通常把困难规避掉：K 凸、F 单调/强单调，SOL(K,F(x,·)) 基本单值，于是 ODE 右端可当作普通向量场处理。这里困难恰恰来自反面：K 可以是 annulus 这类非凸半代数集，F 是非单调多项式，静态 VI 对同一个 x 可能有多个 isolated equilibria，甚至在低维参数集合上出现 continuum solutions。于是“给定 x(t)，u(t) 是什么”不再是一个数值子问题，而是一个选择动力学问题。
关键矛盾是：如果只用 local solver，每个时间步得到哪个 equilibrium 取决于初始化和 solver 偏置，轨迹不可复现，也不代表完整系统行为；但如果枚举所有 equilibrium，又必须知道它们如何随 x 连续延拓、何时分裂/消失。论文的切入点就是把这个矛盾转化为半代数几何问题。
Motivation
作者的核心观察是：非凸 PVI 虽然失去凸 VI 的单值性，但并不是任意病态；由于所有数据是多项式，解图 Graph(S) 是半代数集，因此参数空间具有有限 stratification。换句话说，多解结构可以被几何化、有限化、可枚举化。
已有路线缺两件东西。第一，经典 DVI 依赖凸性/单调性来获得 well-posedness，一旦多解就没有机制解释 branch selection。第二，Moment-SOS / occupation-measure 在混合系统里更多用于全局验证、reachable set、optimal switching 等 ensemble-level 问题，不直接给出某条 Carathéodory trajectory 的 regular selection。
这篇论文试图补的缺口是 trajectory-wise completeness：不仅证明某个解存在，而是能在每个状态枚举所有 admissible controls，并沿着半代数分支连续跟踪；switching 不是外加逻辑，而是轨迹穿越参数分层边界的结果。
Core Idea
核心思想可以概括为：用 polynomial / semi-algebraic structure 取代 convexity 作为 DVI 的组织原则。对静态 parametric PVI，证明 S(x) 的 graph 是半代数集；借助 Hardt triviality，参数域可分解成有限个半代数 strata，每个 stratum 内 solution fiber 的拓扑类型不变。于是 PDVI 可以被视为 ODE 在一个有限分支 equilibrium manifold 上演化：stratum 内做连续选择，stratum 边界处发生离散切换。
这个视角改变了建模方式：prior 往往把 hybrid switching 当成额外 automaton / complementarity logic；这里 switching 来自同一个几何对象 Graph(S) 的 stratification。新的 inductive bias 是“解映射的几何连续性”：算法不再逐点独立求 VI，而是利用相邻时间点 equilibrium set 的连续延拓关系来选择 branch。
本质区别在于完整枚举而非局部求解。SOS-PDVIRK4 的重要性不在 RK4，而在 RK stage 的 VI solver 返回全体 isolated solutions，使得时间推进面对的是一个显式 finite choice set；这使多解动力学可以被系统性探索，而不是被 solver heuristic 隐式决定。
Method
1. 半代数解图与参数分层：解决“多解是否有有限结构”的问题。通过 Tarski-Seidenberg 将 VI 的 forall-z 条件转为半代数 definability，再用 Hardt triviality 得到有限 stratification。核心变化是把 equilibrium multiplicity 从不可控的 set-valued pathology 变成有限可分类的几何对象。
2. regular solution / finite switching 理论：解决“多值右端是否仍能产生可跟踪轨迹”的问题。论文用局部半代数 selection 在 full-dimensional strata 内定义连续控制，跨越 bifurcation set 时允许切换。核心变化是把 Carathéodory solution 扩展为 piecewise-continuous control trajectory，并用 bifurcation set 约束 discontinuity 的位置。
3. LME + Moment-SOS 枚举静态 PVI：解决“每个时间点如何不漏 equilibrium”的问题。LME 在 LICQ/nonsingular 条件下消去乘子，将 KKT/VI 候选点写成 polynomial constraints；Moment-SOS 提供全局 POP 求解与 isolated solutions recovery。核心变化是从 local VI solving 变成 certifiable global enumeration。
4. RK4 + branch selection：解决“枚举出多解后如何产生单条轨迹”的问题。每个 RK stage 都求当前状态的全部 VI 解，再用最近邻规则选择与上一控制最连续的分支。核心变化是把时间积分与 equilibrium consistency 绑定，避免 RK stage 使用不属于当前 S(x) 的控制。
5. boundary-aware hybrid solver：解决“结构边界附近 LME 不稳定/LICQ 可能失败”的问题。远离边界用 LME，靠近 Σ 切到 primal-dual KKT 并减小步长。这个机制更偏 engineering，但对避免 bifurcation 处漏分支或跳错分支是必要的。
Key Insight / Why It Works
最核心贡献是把 PVI solution map 的 semi-algebraic finiteness 作为 PDVI 可跟踪性的根基。方法有效不是因为 RK4 高阶，也不是因为 SOS 本身神奇，而是因为 polynomial structure 给了两种能力：一是解集结构有限，二是每个静态子问题可以用全局代数优化枚举。没有这两点，所谓 complete trajectory tracking 就没有可验证含义。
真正的机制是 latent structure exploitation：系统表面上是 nonsmooth / hybrid / set-valued，但底层 graph 是半代数集，有限 stratification 提供了隐式 mode graph。算法等价于在 test time 重建局部 mode graph，并沿连续分支做 tracking。这更像结构化 test-time compute，而不是传统数值积分方法的改进。
最可能的核心贡献：parametric PVI 分层 + PDVI finite-switching interpretation + 全局 equilibrium enumeration 的组合。单独看每一块都不是全新：Hardt triviality、semi-algebraic selections、Nie 的 LME-SOS、RK4 都已有；但把它们组织成 trajectory-wise 多解 DVI 框架，是这篇的实质新增。
可能只是辅助或 engineering 的部分：RK4 本身不是关键，换成其他高阶 integrator 也应成立；最近邻 selection 是合理 continuation heuristic，但不是由系统物理或最优性原则唯一推出；boundary threshold、step refinement、objective perturbation exclusion 都主要是数值稳定技巧。论文把这些包装成 structure-aware algorithm，有一定合理性，但理论核心不在这些细节。
需要警惕的地方：finite switching 的证明直觉依赖“轨迹与低维半代数集有限相交”，但 ODE 解一般不自动是半代数曲线；如果轨迹沿 bifurcation set 滑行、切触、或产生无限接近事件，文中论证不够扎实。收敛定理也在很强的 B1-B4 假设下退化为“如果 SOS 精确枚举且 branch selection 不错，那么 RK4 有四阶收敛”，这基本是 classical RK4 结论的条件化转述。
Relation To Prior Work
最接近的谱系有三条：经典 DVI / complementarity dynamics，PVI 的半代数解映射与 LME-Moment-SOS 求解，以及 occupation-measure/SOS 的 polynomial hybrid systems 分析。
相对经典 DVI，这篇的本质差异是放弃用凸性和单调性换单值 well-posedness，而是接受多解，并用半代数结构管理多解。它不是 classical DVI 的小扩展，而是把 well-posedness 从 uniqueness 转向 finite branch regularity。
相对 Nie 等 PVI 求解工作，这篇的新增不是静态 PVI solver，而是把静态全局枚举嵌入动态 tracking，并强调分支连续性和 switching geometry。也就是说，PVI solver 是底座，不是论文最原创的数学工具。
相对 Henrion-style occupation measures，这篇更 trajectory-wise：不是优化/验证一族轨迹的 measure relaxation，而是显式构造每条 admissible regular trajectory。这一点是实质区别，也决定了它更适合研究路径依赖、多稳定分支、branch selection，而不是大规模 reachability。
看似新的“structure-aware SOS-PDVIRK4”其实是已有技术重组：Hardt stratification 给结构，Moment-SOS 给全局求解，RK4 给时间推进，nearest-neighbor 给 continuation。实质创新在于把这些工具放进同一个 PDVI 完整轨迹跟踪问题中，并明确提出 polynomial structure 可以替代 convexity 成为非凸 DVI 的组织原则。
Dataset / Evaluation
实验是一个二维 cyber-physical control toy model：annular nonconvex actuation set，特定 polynomial F(x,u)，Duffing-like polynomial plant dynamics。它很好地服务于论文核心论点，因为该例子的解结构可解析、可视化，确实能展示同一 x0 下三条 equilibrium branches、跨越 G(K) 边界时 cardinality 跳变，以及控制轨迹的有限 discontinuities。
但 evaluation 覆盖范围很窄。没有高维系统，没有复杂半代数约束族，没有真实控制平台或真机，没有与 PATH、projection methods、普通 KKT local solver、occupation-measure 方法的定量对比。所谓 convergence study 主要说明在这个例子中 RK4 + 精确/近似 VI solve 表现正常，并不能证明 Moment-SOS 枚举在一般 PDVI 中可扩展。
实验确实支持“结构边界决定 switching”和“完整枚举可揭示多轨迹 dynamics”这两个 claim；但不支持更强的 claim，例如 scalable framework、long-term behavior exploration 的普适性、或复杂工程系统可部署性。长期收敛到共同 equilibrium 的结果也更像该构造系统的性质，不应过度外推。
Limitation
1. 依赖强结构假设：compact basic semi-algebraic K、polynomial f/F、LICQ/nonsingularity、Archimedean condition、generic isolated solutions。离开这些条件，LME 和 finite SOS recovery 都可能失效。
2. scalability 是硬上限：Moment-SOS SDP 规模随变量数、约束数、degree、relaxation order 快速爆炸。当前实验 m=n=2，不能说明高维控制或复杂约束下可用。未来用 TSSOS 是合理方向，但这篇没有实际解决。
3. continuum solution region 被基本绕开：论文承认初始化在 Σ 或 S(x) 无限时超出 finite-solution paradigm。可是非凸 VI 中 continuum equilibria 正是重要病态之一；当前框架更像“almost everywhere isolated branch tracker”。
4. finite-switching 论证存在潜在漏洞：文中把 trajectory image 当作 semi-algebraic curve 使用，但多项式 ODE 的解通常不是半代数函数。若轨迹沿 bifurcation set 滑动或出现无限次横穿，finite intersection 结论不自动成立。这里文中未充分说明。
5. branch selection 缺乏语义：最近邻控制选择保证数值连续，但不一定对应物理稳定性、能量最小、控制策略或 equilibrium selection principle。它把“选择哪条未来”从 solver 随机性转移到 heuristic continuity bias。
6. 收敛定理条件化较强：B3 的 uniform SOS exactness、B4 的 consistent branch selection 都是实际最难保证的部分。一旦 relaxation order 不够或解间距很小，误选分支造成的误差不是普通 RK4 truncation error 可以覆盖的。
7. boundary-aware 部分工程性较强：如何准确获得 Σ、如何设阈值、如何判断 cardinality change，在一般问题中并不清楚。增益来源不清：可能主要来自低维例子中边界结构已知，而不是算法自动理解了复杂 stratification。
Takeaway
- 1. 值得记住的核心 insight：非凸 DVI 不一定要回到凸性/单调性；如果数据是 polynomial，semi-algebraic geometry 可以提供另一种 well-posedness 语言——finite branch structure rather than uniqueness。
- 2. 对相关方向可迁移的点：在多解 equilibrium-coupled dynamics 中，先枚举完整 instantaneous solution set，再做 continuity-based branch tracking，比每步调用 local solver 更接近真实系统结构。
- 3. 未来真正值得做的是 scalable stratification-aware solvers：利用 sparsity、active-set continuation、certified branch tracking，而不是每个 RK stage 从头跑高阶 SOS。
- 4. 这篇推动的是问题表述和结构化计算范式，而不是一个已经工程可扩展的 solver。
一句话总结
这篇论文把非凸 PDVI 的多解动力学从局部求解器的偶然输出，重新表述为半代数 equilibrium manifold 上的有限分支跟踪问题，并用 Moment-SOS 枚举加连续选择给出一个理论上完整但目前主要适用于低维结构化系统的 trajectory-wise 框架。