精读笔记

Problem Setting

《Discovering the Kalman-Bucy-Koopman Filter》（arXiv preprint / 2026-06-30）处理的是连续时间非线性系统的状态估计，目标不是近似 Kalman filter 的一步更新，而是给 Mortensen maximum-likelihood / MAP 滤波找到一个具有 Kalman-Bucy 风格的有限维 Riccati 实现。真正困难在于：非线性滤波的精确 density evolution 是 SPDE，Mortensen/HJ 路线虽更接近 deterministic optimal control，但值函数一般不可闭合；而 EKF/UKF/PF 的近似要么局部、要么采样昂贵、要么缺少清晰代数结构。本文要解决的关键矛盾是，在保留非线性状态结构的同时，让估计增益仍来自一个结构化的 Riccati-type evolution。

Motivation

已有 Koopman observer 或 EDMD-based estimator 通常先学一个 lifted linear model，再在 lifted space 里做 observer/filter；这类方法的弱点是最优估计结构并不自然出现，lift 的选择也更多是 regression convenience。作者的出发点是另一个缺口：Mortensen HJ value function 和 Koopman eigenfunction spectral theory 之间存在结构连接，如果 value function 在 Koopman coordinates 中是二次的，那么非线性估计可以像线性 Kalman-Bucy 那样被 Riccati 化。换言之，缺的不是又一个非线性 observer，而是一个能把 HJ optimality、Koopman linearity 和 Riccati algebra 同时闭合的坐标系统。

Core Idea

论文最核心的思想是：不要在线性化原系统，也不要直接传播后验密度，而是把 MAP 滤波的负 log posterior / value function 写在 Koopman principal eigenfunctions 上。由于 principal eigenfunctions 满足 ∂Φ/∂x f(x)=ΛΦ(x)，非线性 drift 在 Φ 坐标中被谱线性化；如果测量 h(x) 也能由这些 eigenfunctions 或其乘积展开，那么 HJ 方程中的项可以按 Φ 的二次、一次、常数项匹配，得到 P、s、r 的 Riccati-type ODE。 本质区别在于建模对象变了：prior Koopman filtering 多在 observable dynamics 上做 linear prediction；KBK 在 value function 上做 Koopman spectral parameterization。它引入的 inductive bias 是“估计误差能量在流的谱坐标中近似二次”，而不是“状态局部近似线性”。这使得非线性被编码进 Φ 的几何形状，Riccati 结构则在谱坐标中保留。

Method

关键机制可以压缩为以下几步。第一，用 Mortensen minimum-energy cost 替代 density filtering，直接优化最可能轨迹的终点；这使问题进入 HJ PDE，而不是 Zakai/KS SPDE。第二，在 Koopman eigenfunction 坐标中重写初始 cost 和过程噪声 metric，引入由 Jacobian 诱导的 R1(x)；这一步保证坐标变换下的能量度量一致，但也把正则性和可逆性假设推到 Φ 上。第三，假设 h(x)=CΦ(x) 时，二次 value ansatz 精确满足 HJ，得到 KBK 的 P/s/r 方程和状态估计 ODE；这是理论闭合的核心。第四，当 h 不在 principal span 中，用 projection residual 给 value/state error bound，并用 higher-order eigenfunction products 扩充表示；这更像受 Koopman algebra 约束的输出 lift。第五，用 path-integral 和 characteristics-inspired basis 计算 eigenfunctions；这不是滤波本身的核心创新，但决定方法能否落地。

Key Insight / Why It Works

最重要的 insight 是：Kalman-Bucy 的核心不是线性状态方程本身，而是 value function/posterior log-density 在某个坐标中保持二次闭合。Koopman principal eigenfunctions 提供了一个候选坐标，使确定性流线性化；一旦 measurement 也能在该坐标中线性表示，HJ 方程就退化为 Riccati algebra。这解释了为什么方法在结构匹配的系统上有效：它不是近似非线性为线性，而是找到了让非线性系统“看起来线性”的 observables。 真正的贡献是把 Mortensen HJ filtering 和 Koopman spectral coordinates 对齐，而不是 path-integral eigenfunction computation 本身。higher-order products、basis construction、data-driven regression 更像辅助工程：它们扩大可表示性，但也可能只是把误差压到更大的 dictionary 和更好的数据覆盖里。尤其 data-driven 例子中，KBK 相对 mismatched EKF 的优势归因并不干净：EKF 被故意给错参数，而 KBK 使用来自真实系统的数据学习谱表示，增益可能主要来自 data/model access advantage，而非滤波方程本身更优。文中未充分说明 eigenfunction approximation error 如何进入最终滤波稳定性。

Relation To Prior Work

它最接近三条线：Mortensen minimum-energy filtering、Koopman eigenfunction-based nonlinear coordinates、以及 Koopman observer / EDMD estimator。与 EKF 的差别是它不在当前估计点线性化 f,h，而是预先构造谱坐标；与 UKF/PF 的差别是它不传播样本分布，而传播一个结构化 value function；与 Zakai/KS 的差别是它瞄准 MAP mode，不是完整 posterior evolution。 看似新的部分中，Koopman lift、eigenfunction products、least-squares output projection 都有已有谱系；实质新增在于把 HJ 方程的二次 ansatz 放在 principal eigenfunction 坐标里，并给出 exact-when-h-in-span 的 Riccati 闭合。它属于 operator-theoretic nonlinear estimation，而不是 generic learned filtering。严格说，这篇不是把 nonlinear filtering 完全解决了，而是指出了一类在 Koopman spectral coordinates 中可以恢复 Kalman-Bucy 结构的系统/近似。

Dataset / Evaluation

评估主要是仿真，覆盖解析系统、quadrotor 垂直动力学参数 mismatch、车辆摩擦估计三个低到中低维场景。它们能够展示方法在结构友好、数据可得、eigenfunction 可近似的情况下可工作，也展示了相对 EKF 的鲁棒性。但 benchmark 没有充分验证全局 nonlinear filtering claim：没有高维系统，没有真实硬件，没有与 UKF/PF/ensemble Kalman/neural filters 的系统比较，也没有对 eigenfunction 误差、domain violation、Hessian conditioning 的消融。 quadrotor 例子中 EKF 使用 mismatched 参数，而 KBK 使用真实系统数据学习，这个设置更像验证“data-informed spectral model 比 wrong analytic model 好”，不等价于证明 KBK filter structure 本身优越。车辆例子支持 data-driven feasibility，但更接近 controlled simulation 下的 representation fitting；真实 deployment 中摩擦突变、激励不足、未覆盖工况会不会破坏 Φ 表示，文中未充分说明。

Limitation

最核心限制是方法把非线性滤波难题转移成 Koopman eigenfunction 构造难题。principal eigenfunctions 不一定全局存在，也不一定在感兴趣状态域内可逆或数值稳定；一旦轨迹离开 eigenfunction 有效子域，KBK 的理论基础会塌。ΠΦ 需要 Hessian 信息且要可逆，这在高维、有噪数据或 learned eigenfunctions 下非常脆弱。 第二个限制是输出表示假设强。exact result 依赖 h∈span{Φ}；不满足时的误差界依赖 sup projection error、measurement residual bound 和 local strong convexity，这些在实际系统中很难验证。higher-order products 能降低 projection error，但可能带来维数增长和 conditioning 问题，本质上接近结构化 polynomial lifting。 第三，MAP 目标不是 MMSE；posterior 非高斯/多峰时，mode 可能不是最有用估计。文中虽承认 bias/CRLB 问题，但实验没有真正检验多峰 posterior 或强非高斯 regime。 第四，data-driven claim 的边界不清。若 KBK 需要大量覆盖真实 dynamics 的轨迹来学 Φ 和 C，那么泛化能力可能主要来自数据覆盖；所谓 robustness 可能是模型 mismatch 设置下的 access advantage。增益来源不清。

Takeaway

• 第一，最值得迁移的 insight 是：非线性估计可以从“传播分布/线性化状态”转为“寻找使 value function 二次闭合的坐标”。
• 这比常规 Koopman observer 更接近最优控制/滤波结构。
• 第二，Koopman eigenfunctions 在这里的价值不是 generic lifting，而是作为 HJ value-function coordinates；如果某个领域能找到类似的 intrinsic spectral coordinates，也可能恢复 Riccati-like inference。
• 第三，未来真正关键的问题不是再堆 higher-order basis，而是可验证地学习局部可逆、稳定、带 Hessian 控制的 eigenfunctions，并把 eigenfunction approximation error 接入滤波稳定性理论。

一句话总结

这篇论文把 Mortensen/HJ 非线性 MAP 滤波重新表述为 Koopman eigenfunction 坐标中的 Kalman-Bucy 型 Riccati 估计，是一次有理论野心的 operator-theoretic reparameterization，但实际能力上限主要受 eigenfunction 存在性、可学习性和输出可表示性制约。