精读笔记
Problem Setting
这篇论文解决的是稳定线性控制系统在静态非线性输出下的 GOE 判定问题：给定周期输入 v(t) 与同均值常值输入 \bar v，比较周期稳态输出均值和常值稳态输出。表面上这是周期控制是否有收益的问题，实质上是在问：当状态动力学本身完全线性、平均响应没有任何一阶优势时，非线性读出是否足以制造平均性能差异。
关键矛盾是线性系统的平均运算会把周期输入压回到同均值平衡点，因此任何线性输出都不可能产生 GOE；但输出非线性又使“先输出再平均”和“先平均再输出”不交换。以前在一般非线性/收缩系统里的 GOE 分析容易把增益来源混在动力学非线性、轨道变形和输出非线性里，本文刻意拿掉动力学非线性，把问题压缩到最小机制。
Motivation
已有路线不够的地方在于，GOE 在非线性系统中通常是高阶效应，很难判断增益到底来自 entrainment、本体动力学、输入波形，还是输出函数本身。作者的核心观察是：即使状态动力学完全线性，周期稳态轨道仍会在平均平衡点周围产生波动；只要输出函数有曲率，这个波动就会通过 Jensen gap 变成非零平均输出差。
因此本文的动机不是提出更强的 entrainment 理论，而是隔离 GOE 的一个基本来源：输出几何。缺口在于此前缺少一个既能给出全局符号判据、又能给出频域设计规则的最简模型。Wiener 系统正好是这个缺口的干净载体。
Core Idea
核心思想可以概括为：把 GOE 看成“周期轨道围绕其均值的输出几何偏差”，而不是看成“周期输入带来的动态性能提升”。Hurwitz 线性系统保证周期轨道 \gamma^v 的时间平均等于 e^{\bar v}，所以所有一阶平均效应都被精确消掉；剩下的就是 h 在该均值附近如何响应状态波动。凸性给正 Jensen gap，凹性给负 Jensen gap，二阶局部极限由 Hessian 与轨道协方差决定。
这和 prior 的本质区别是建模归因的改变：不是从非线性动力学中寻找 entrainment 后的性能改善，而是把动力学当作一个稳定线性滤波器，把周期输入转换成状态空间中的波动分布，再由输出非线性的曲率决定收益。这个 inductive bias 很强：系统部分负责生成 covariance，输出部分负责决定 covariance 的价值。
Method
第一，论文证明稳定线性系统的周期稳态平均恒等式：\bar \gamma^v=e^{\bar v}。这一步解决了公平比较中的基准问题，也解释了为什么线性输出 GOE 必为零；它是全文所有结论的结构支点。
第二，对光滑输出做小幅周期扰动展开。由于扰动轨道均值为零，一阶项消失，主导项是 (1/2)tr(\nabla^2 h(e^{\bar v})\Sigma^w)。这一步说明 GOE 的局部符号由输出曲率和线性系统诱导的状态波动协方差共同决定，而不是由输入均值决定。
第三，用 Jensen 和 Bregman divergence 给出全局几何刻画。凸/凹输出在 controllable subspace 上即可保证任意周期输入的 GOE 符号；Bregman 表达把 GOE 写成平均 D_h(\gamma^v(t),e^{\bar v})，说明它本质是非线性输出相对切平面的平均偏离。
第四，在二次输出 h(x)=x^TQx 下，GOE 变成 tr(Q\Sigma^v)，并可写成频域和式 1/2 \sum_k c_k^T Re(G^*(ik\omega)QG(ik\omega)) c_k。这里解决的是可计算性和设计问题：符号判据变成频率逐点半正定性，能量约束下的最优输入变成最大特征值对应的单频单方向激励。
Key Insight / Why It Works
最重要的 insight 是：GOE 在这个模型里不是 entrainment 的神秘收益，而是 Jensen gap。线性稳定动力学只负责保证存在唯一周期轨道，并把输入频率成分滤波到状态空间；真正决定正负的是输出函数在可达状态子空间上的曲率。换句话说，系统的“动态性”只决定波动往哪里投影、投影多大，输出的“几何性”决定这些波动是否有价值。
论文最核心的贡献是把 GOE 拆成两个可分离对象：状态波动协方差 \Sigma^v 与输出几何 h/Q。这比一般非线性 GOE 分析更干净，因为它消除了动力学非线性带来的归因混乱。Bregman divergence 不是装饰性解释，它准确说明了为什么凸性给全局符号结论；而二阶展开只是这个全局几何事实的局部版本。
频域部分更偏 engineering，但很有用。对二次输出，不同谐波独立贡献不是深层非线性发现，而是线性系统与正交 Fourier 基的直接结果；真正有迁移价值的是把“周期控制设计”转成“选择频率-方向以最大化输出加权传递增益”。这部分的增益不是 scaling，也不是数据覆盖，而是更好的结构归因和经典频域工具的重用。
需要注意：所谓最优周期激励只在非常理想的能量约束和二次输出下成立。它不是一般周期控制最优性的证据，而是该简化模型下的闭式解。
Relation To Prior Work
最接近的是 GOE/entrainment 在 contractive systems、ribosome flow model、weakly contractive bilinear systems 中的分析，以及经典 Wiener 系统和频域线性系统理论。与 contractive/nonlinear GOE 工作相比，本文不试图解决周期轨道存在和全局吸引的难题，而是把这部分交给 Hurwitz 线性系统，从而专注于平均性能差异的来源。
与线性系统理论相比，本文新增的信息不是 transfer function 本身，而是把 transfer function 与非线性输出曲率连接成 GOE 判据。与 Jensen 类结果相比，本文的实质创新在于把 Jensen gap 嵌入到稳定线性控制系统的周期稳态结构中，并指出只需在 controllable subspace 上要求凸/凹。
一些看似新的部分其实是已有思想重组：Bregman divergence、KL 例子、谐波正交分解、最大特征值选频都是成熟工具。实质新增在于它们被组织成一个 GOE 归因框架：平均恒等式 → 几何 gap → 协方差/频域设计。
Dataset / Evaluation
这不是数据驱动论文，evaluation 主要是两个应用例子而非 benchmark。RLC 电路覆盖的是凸二次能量输出，验证正定 Q 下周期输入必然提升平均储能，并展示频率响应公式能解释最优频率不等于共振频率。药效模型覆盖的是正系统加凹 Emax 输出，验证周期给药在同平均剂量下不能优于常值给药。
这些例子支持论文的核心理论 claim：输出凸性/凹性支配 GOE 符号，频域公式可用于设计。但它们不构成跨场景验证，也没有真实系统实验或带约束的部署测试。尤其药效例子忽略峰值毒性、给药离散性、药代非线性和临床目标函数，不能直接推出实际周期给药无价值；它只说明该线性 compartment + 凹 readout 模型内 GOE 为负。
Limitation
核心限制是成立条件非常强。A 必须 Hurwitz，状态方程必须线性，输出只依赖 x，且比较只固定输入均值。只要动力学中出现非线性项，周期轨道平均就不再等于平均输入对应平衡点，GOE 会多出平均漂移项；此时输出曲率不再是唯一来源，本文的干净分解会被破坏。
凸/凹判据是充分条件，不是必要条件；二次情形给了必要充分频域条件，但这只覆盖 h=x^TQx。对于一般非凸输出，GOE 的全局符号仍可能高度依赖输入幅度、轨道覆盖区域和系统滤波结构。
最优激励结论的上限也明显：能量约束下单频最优来自线性叠加和二次目标，一旦加入输入非负性、幅值上限、状态安全约束、稀疏脉冲给药、执行器带宽或输出饱和，最优波形可能完全改变。文中未充分说明这些现实约束下频域结论的鲁棒性。
此外，文中应用只验证了理论方向，没有展示真实工程/生物系统中的模型误差敏感性。增益来源在本文模型内很清楚，但在实际系统里可能被未建模非线性、约束和噪声主导。
Takeaway
- 第一，周期控制相对常值控制的收益可以被干净地归因为“状态波动 × 输出曲率”，这比泛泛讨论 entrainment 是否有益更可迁移。
- 第二，controllable subspace 是关键：输出函数在全空间是否凸不重要，真正相关的是周期轨道实际能访问的子空间上的几何。
- 第三，对二次性能指标，GOE 可以直接转成频域加权增益问题；这给设计周期激励提供了比时域仿真更透明的工具。
- 第四，未来真正值得做的是把这个分解推广到弱非线性动力学，弄清平均漂移项和输出曲率项如何竞争；这比继续在线性 Wiener 系统里扩展示例更有价值。
一句话总结
这篇论文把稳定 Wiener 系统中的 GOE 从一个周期控制现象还原为输出几何诱导的 Jensen/Bregman gap，并在二次输出下将其闭式化为经典频域加权能量设计问题。