精读笔记
Problem Setting
【Nonsmooth High-Order Averaging Theory with Application to Extremum Seeking Optimization and Control】（arXiv preprint / 2026-06-02）
这篇论文实际解决的是：对于 \dot{x}=\epsilon f(x,t,\epsilon) 这类周期非线性系统，当 f 对状态 x 非光滑时，如何仍然得到类似 smooth averaging 的高阶平均系统，尤其是二阶平均公式，并把它用于控制仿射 ESC 的稳定性分析。
关键矛盾在于：高阶 averaging 不是简单积分平均。二阶及以上项需要对向量场做状态方向展开，需要 Jacobian、链式法则、Lie bracket 型组合和 near-identity coordinate transform。smooth 情况下这些是标准操作；nonsmooth 情况下 classical derivative 不存在，Clarke derivative 虽然集合值稳定，但不满足足够 sharp 的 calculus rules，实际推导和计算都会卡住。
以前方法的问题不在于完全不能处理非光滑，而是处理得太局部或太结构化：一阶 averaging 可以绕开很多导数问题；hybrid/discontinuous averaging 通常关注特定类系统；ESC 中的 non-C2 Lie bracket 多局限于原点非光滑或特定控制律。本文试图把问题提升到 general canonical form，并且给出用户能真正代入计算的二阶公式。
Motivation
作者的核心动机可以概括为：非光滑 averaging 缺的不是 stability language，而是可计算的高阶微分代数。
Clarke/Filippov 工具可以给 existence、inclusion 和 Lyapunov stability，但它们不天然给出 smooth averaging 那种干净的二阶项公式。尤其在 ESC 中，Lie bracket approximation 本质上就是 higher-order averaging；如果 objective 或 vector field 非光滑，传统 Lie bracket 推导失去合法性，而只靠 smoothing approximation 又会改变系统局部结构。
lexicographic differentiation 被引入的原因很直接：它在 L-smooth 函数类上既覆盖 abs/min/max/convex/PC1 及其复合，又有 sharp sum/product/chain rules。换句话说，它不是为了更一般的 generalized derivative，而是为了让 averaging 推导中的“符号微分流程”重新可用。这个选择很务实，也构成本文相对已有 nonsmooth averaging 的主要技术切入点。
Core Idea
论文真正的核心思想是：保留 smooth averaging 的几何/代数结构，只替换其中不可用的 Jacobian。具体说，用 lexicographic derivative J_L f 作为 nonsmooth Jacobian-like object，并据此定义 L-Lie bracket [f,g]_L=J_L g f-J_L f g。这样二阶平均项仍然可以写成 f2、[w,f1]_L、[w,g1]_L 的周期积分，而不是退化成难以操作的 set-valued expression。
这改变了非光滑系统的建模方式：不是把非光滑性看成必须 convexify 的不确定性，而是把它看成 piecewise/lexicographic 可解析结构。inductive bias 是“非光滑函数在合适方向序下仍有可计算的局部线性化”。因此方法比传统 Clarke-based route 更适合做高阶公式推导，也比针对某个 ESC 控制律的 ad hoc Lie bracket 更 general。
本质区别在于：prior 多数是在非光滑系统理论里寻找 averaging stability；本文是在 averaging calculus 本身里替换微分算子。这一点是实质性的。
Method
1）L-Lie bracket：解决高阶平均项中 Lie bracket 无法定义/计算的问题。它把 classical Jacobian 替换成 lexicographic derivative，并在 smooth 情况自动退化为标准 Lie bracket。核心变化是二阶 averaging 的 algebraic form 被保留下来。
2）nonsmooth near-identity transformation：解决 smooth averaging 中坐标变换在 Lipschitz 情况下是否合法的问题。作者用 Clarke nonsmooth implicit function theorem 证明 x=y+\epsilon w(y,t,\epsilon) 在小 \epsilon 下局部可逆。这个步骤是理论闭环的必要条件，否则平均系统只是形式展开。
3）L-smooth Taylor-like expansion：解决 f(y+\epsilon w,t) 在非光滑点无法 Taylor 展开的问题。这里用 lexicographic directional derivative 给出一阶 remainder 控制，使得 O(\epsilon^2) 项能被组织出来。它是二阶平均公式成立的关键技术桥。
4）二阶平均公式：在 A4=PWC∞ 条件下定义 g2(y)=1/T ∫[f2+[w,f1]_L+[w,g1]_L]dt。 这一步把非光滑性集中到 L-derivative/L-Lie bracket 中，而平均结构仍与 smooth theory 对齐。
5）Filippov 解与误差界：由于 g2 可能是 PWD∞，二阶平均系统不一定 classical，需要 Filippov solution。O(\epsilon) 误差界较直接；O(\epsilon^2) 需要 g1=0 和 crossing hypothesis，说明更强误差不是免费得到的。
6）ESC 应用：把控制仿射 ESC 经高频 rescaling 写成 averaging canonical form，通过二阶平均得到 \dot z=-A\nabla_L J(z)。这相当于把 smooth ESC 的 gradient flow 替换成 generalized-gradient flow。
Key Insight / Why It Works
最重要的 insight 是：高阶 averaging 的难点是 calculus compatibility，而不是 generalized derivative 的存在性。lexicographic differentiation 的价值在于它提供了 sharp chain rule，这正好匹配 near-identity transform 和 Lie bracket expansion 的需求。Clarke derivative 给的是集合，理论上稳健但计算上粗；L-derivative 给的是方向序相关的 single Jacobian-like object，牺牲某种 set-valued 完整性，换来可计算的高阶公式。
方法有效的真正原因有三层。
第一，L-smooth/PWC∞ 函数在几乎处处仍是 smooth 的，非光滑集合测度为零；因此大部分积分平均项实际上沿 smooth pieces 计算。L-derivative 主要负责跨 piece 或位于 kink 处时给出一致的导数选择。这意味着理论并不是在完全任意 Lipschitz 世界里工作，而是在 piecewise-smooth structure 上工作。
第二，periodic averaging 把高频时间依赖压缩成低频 drift；二阶项捕捉的是 fast oscillation 与 state-dependence 的耦合。ESC 中这个耦合正好产生 gradient/generalized-gradient flow。也就是说，ESC 的“优化能力”仍来自经典 dither-induced Lie bracket 机制，只是 gradient 被替换为 L-gradient。
第三，误差分析依赖小参数 \epsilon 和紧集有界性。这里的 closeness 很大程度来自 averaging scaling，而不是非光滑工具本身；非光滑工具保证的是展开合法和公式可计算。换言之，性能增益不是来自新的控制律，而来自把已有高频 averaging 机制扩展到 nonsmooth calculus。
最可能的核心贡献是二阶平均项的 closed formula 与 L-Lie bracket，而不是 ESC 仿真。ESC 部分更像是展示该 calculus 的自然应用。Law 1/2/3 推出相同 generalized-gradient flow 的结果有用，但本质上延续 smooth ESC 设计谱系；新增点是允许 J 或 b_i 非光滑且仍能形式化稳定性。
需要注意：L-derivative 在非光滑点依赖方向矩阵，文中默认 identity/natural choice。这个选择在计算上自然，但在动力系统解释上并非完全无代价。不同 lexicographic order 可能导致 kink 处不同 vector field selection。论文没有充分讨论这种 selection dependence 对 averaged dynamics 和 control behavior 的影响。
Relation To Prior Work
最接近的路线有三条：classical high-order averaging / near-identity transform；Lie bracket averaging for control-affine systems and ESC；nonsmooth/hybrid averaging with Clarke-Filippov tools。
相对 classical averaging，本文不是重写 averaging 理论，而是把其微分算子替换成 lexicographic derivative，并证明在 L-smooth/PWC∞ 条件下相同结构仍成立。smooth 情况自动恢复，这是一个好的 consistency check。
相对 Lie bracket ESC，本文的实质新增是说明 first-order Lie bracket approximation 等价的二阶 averaging 可以在 nonsmooth objective/vector field 下闭式计算。Scheinker-Krstic 的 non-C2 ESC 更像特定非光滑控制律的扩展；Suttner 的随机方向 nonsmooth optimization 是另一种绕开 gradient regularity 的方式。本文则给 deterministic generalized-gradient drift。
相对 Clarke/Filippov nonsmooth analysis，本文没有试图用 Clarke generalized Jacobian 直接构造高阶公式，因为那会丢失 sharp calculus。lexicographic differentiation 是真正的新工具输入。可以说本文属于“把 computational generalized differentiation 引入 averaging/control theory”的谱系，而不是单纯的 nonsmooth stability paper。
看似新的部分中，near-identity transform、averaging error、practical stability transfer 都是 classical theory 的非光滑化移植；实质创新集中在 L-Lie bracket、二阶公式，以及证明这些对象在 PWC∞/PWD∞ 下可控。
Dataset / Evaluation
这是一篇理论控制/数学论文，evaluation 主要是数值仿真而非 dataset benchmark。仿真覆盖两类：一个二维 nonsmooth NLTP 例子，目标是验证即使系统不稳定，二阶平均轨迹仍能跟随原高频系统；若干一维 ESC 例子，目标函数包含 abs/min/多段结构，用来展示 generalized-gradient averaged system 能捕捉原 ESC 行为。
这些实验支持的 claim 是有限但对口的：它们说明公式可操作，L-gradient 能处理 kink，且高频增大时 averaged approximation 变好。它们没有验证高维 scalability、复杂 switching geometry、实际测量噪声、真实 plant、约束优化或多极小非凸场景。
因此 evaluation 更像 proof-of-concept，而不是经验性能证明。对这类论文这并不是致命问题，因为主要贡献在理论；但如果把 ESC 部分理解为实际 model-free optimizer 的泛化能力，证据是不够的。文中未充分说明这些 nonsmooth ESC laws 在真实系统中面对噪声、delay、actuator saturation 时是否仍保留 generalized-gradient drift。
Limitation
主要限制不是“只做到二阶”这么表面，而是理论成立依赖的结构相当明确。
1）函数类限制：虽然 L-smooth 覆盖很广，但二阶部分实际加强到 PWC∞，并要求 piecewise 结构良好。任意 locally Lipschitz 函数并不在可操作范围内。所谓 nonsmooth generality 主要是 piecewise/lexicographic 可微结构上的 generality。
2）方向选择问题：L-derivative 依赖 directions matrix。论文后续基本使用 identity/natural derivative，但非光滑点处这种选择会决定 bracket 和 averaged vector field。若系统轨迹长期贴着 switching surface，selection 的物理含义可能很关键；文中未充分说明。
3）误差界依赖小参数与有界性：O(\epsilon) 或 O(\epsilon^2) closeness 都是在 t=O(1/\epsilon)、紧集内、有解存在的条件下。实际控制中如果 averaged system 只给 practical stability，原系统的 transient 可能仍受高频幅值、噪声和 actuator 限制影响。
4）二阶 O(\epsilon^2) 需要 g1=0 和 crossing hypothesis。crossing hypothesis 排除了很多 sliding / grazing / tangential switching 行为，而这些恰恰是非光滑系统中最麻烦的动力学。换句话说，最强结论避开了最难的非光滑现象。
5）ESC 部分维度和复杂度有限。仿真是一维或二维，objective 虽然形式复杂但结构可控。没有证据表明该方法在高维 nonsmooth nonconvex objective 上比已有随机/gradient-free ESC 更稳健。
6）高阶扩展未真正完成。标题是 high-order，但闭式公式到二阶；更高阶理论是否可计算、表达是否爆炸，文中没有充分说明。可能二阶是 sweet spot，更高阶会受到 lexicographic derivative 嵌套和 piecewise region combinatorics 限制。
Takeaway
- 1）最值得迁移的 insight：在非光滑控制理论里，选择 generalized derivative 时不应只看数学包容性，还要看它是否支持目标推导所需的 calculus rules。
- lexicographic differentiation 在需要 chain rule 的高阶展开中非常合适。
- 2）这篇真正推动的是 nonsmooth averaging 的“可计算化”：从 set-valued/stability-oriented analysis 走向可以写 closed-form high-order averaged vector field 的 calculus。
- 3）对 ESC 来说，本文说明很多 smooth Lie-bracket optimizer 可以系统地替换成 generalized-gradient optimizer，而不必为每个 nonsmooth objective 重新做 ad hoc 分析。
一句话总结
这篇论文把 lexicographic differentiation 引入高阶 averaging，用可计算的 L-Lie bracket 恢复了非光滑系统的二阶平均公式，并将 smooth ESC 的 Lie-bracket/gradient-flow 机制系统扩展为 nonsmooth generalized-gradient-flow 机制。