精读笔记
Problem Setting
《From the Linear Quadratic Regulator (LQR) to the (Deterministic) Kalman Filter in Two Easy Steps》(arXiv preprint / 2026)处理的是确定性 Kalman filtering 的结构化推导问题:给定线性系统 ˙x=Ax+w, y=Cx+v,只观测 y[0,T],在所有与动力学和测量方程一致的轨迹中,选择过程扰动、测量误差和初始状态加权 L2 能量最小的那条轨迹。难点不是最小二乘本身,而是这个问题不像标准 LQR:目标函数因 y-Cx 含有线性项和常数项,初始状态 xi 也是优化变量,当前/终端状态也未知。换句话说,估计问题同时破坏了 LQR 的两个标准结构:pure quadratic objective 和 fixed initial condition。已有 deterministic derivation 通常借 tracking duality、time reversal、value-function optimization 混在一起解释,正确但不够模块化;这篇 note 的价值在于把这些障碍拆开,使 Kalman observer 的来源在 Riccati block structure 中变得透明。
Motivation
作者不试图证明 stochastic Kalman filter 的最优性,而是绕开概率叙事,把问题放回 Willems 式 deterministic least-squares filtering:观测 y 是事实,其余 w、v、xi 是解释,最优估计就是最小能量解释。已有 stochastic derivation 的问题是 covariance / white noise 叙事对很多控制背景读者并非必要,且会遮蔽一个更基本的事实:Kalman filter 本质上是一个 cost-to-arrive 递推。作者的核心观察是 affine-quadratic optimal control 中的线性项不是 nuisance;在齐次坐标后,它成为扩展 Riccati 矩阵的 off-diagonal block。这个 block 在 LQ tracking 中对应 anti-causal feedforward,在 estimation 中则对应 causal information vector。缺口正是:如何从标准 LQR 的 algebra 中无缝读出这个 observer dynamics。
Core Idea
论文真正核心的思想是把 deterministic filtering 看成 final-condition affine-quadratic LQR,而不是把 Kalman filter 当作概率估计器来推。首先,测量 residual cost (y-Cx)^*V(y-Cx) 造成的 affine-quadratic 结构通过增加一个恒等于 1 的状态 α 转成高一维纯二次 LQR;这一步不是技术花招,而是把外生信号 y 对优化的影响编码进扩展状态代价矩阵的 off-diagonal block。其次,用 final-condition LQR 的 cost-to-arrive 而不是 standard LQR 的 cost-to-go:因为估计的自然问题是“到达当前候选状态需要付出多少解释能量”,而不是“从已知初态出发未来要付出多少控制代价”。
和 prior 的本质差别在于组织信息流的方向:tracking 问题给定初态、参考信号 r 影响未来控制,因此线性项产生 backward-in-time 的 feedforward;filtering 问题当前估计应当随观测在线更新,因此 final-condition formulation 使同样的 off-diagonal block forward-in-time 演化,成为 observer 的内部状态。这个视角把 Kalman filter 的 causality 归因于 cost-to-arrive,而不是概率条件化本身。
Method
1. Final-condition LQR:作者先并列 standard LQR with initial condition 和 LQR with final condition。后者给出 forward DRE:S˙=-A^*S-SA-SBR^{-1}B^*S+Q,S(0)=X,并把 value function 解释为 cost-to-arrive。它解决的问题是估计中“初值未知但当前候选状态可被评价”的结构;核心变化是从 backward cost-to-go 切换到 forward accumulated evidence。
2. Homogeneous coordinates:对 affine-quadratic functional x^*Hx+h_1^*x+x^*h_1+h_0,引入 ξ=[x;α], α≡1,将其写成 ξ^*Ĥξ。它解决的问题是 LQR 只能直接处理纯二次目标,而测量 tracking / filtering cost 含线性项;核心变化是线性项变成扩展 Riccati 的 off-diagonal block p1 或 s1。
3. Riccati block partition:扩展 DRE 的 (1,1) block 给出常规 Riccati,(1,2) block 给出线性系统,(2,2) block 只影响 value scalar。关键不是 p0,而是 s1:在 final-condition case 中 s1 forward 演化,并通过 x_hat=-S^{-1}s1 与估计状态绑定。
4. 对未知终端状态优化:先假设 x(t)=xf 已知,得到 cost-to-arrive V_t(xf)=xf^*Sxf+2s1^*xf+s0;再对 xf 最小化,得到 x_hat(t)=-S^{-1}s1。把这个关系对时间求导并代入 S 与 s1 的方程,恢复 observer:x_hat˙=Ax_hat+S^{-1}C^*V(y-Cx_hat)。这一步是从 batch minimum-energy smoothing-like statement 到 online filtering dynamics 的关键桥。
Key Insight / Why It Works
最重要的 insight 是:Kalman filter 的动态 observer 结构并不是额外假设出来的,而是 affine-quadratic LQR 中线性项的 Riccati off-diagonal dynamics 在 final-condition / cost-to-arrive 方向上的表现。对于纯二次 LQR,optimal control 是 memoryless state feedback,因为 value function 只有 quadratic curvature;一旦目标含外生信号导致线性项,value function 必须携带一个一阶项。这个一阶项就是“信息向量”的前身。tracking 中它由未来 reference 驱动,所以 anti-causal;filtering 中它由过去 measurement 驱动,所以 causal。
方法成立的根本原因是线性-二次结构的闭合性:quadratic value function 在 Hamilton-Jacobi / dynamic programming 下保持 quadratic,affine-quadratic value 在齐次扩展后也保持 quadratic。因此所有未知轨迹优化都被压缩进有限维的 S 和 s1。进一步对当前状态做最小化,相当于把 cost-to-arrive 的 bowl 的中心作为估计值,而 bowl 的曲率 S 表示确定性意义下的 certainty / information。
我认为实质贡献在于概念重排而非新滤波算法:final-condition LQR + homogeneous coordinates + block partition 这三者都不是新技术,但组合后清晰解释了为什么 deterministic minimum-energy filtering 与 stochastic Kalman-Bucy information filter 完全同构。辅助部分是 LQ tracking 对比;它不是必要推导,但很好地说明同一 affine term 在不同边界条件下产生完全不同的时间方向。这里没有 scaling、data、retrieval、benchmark trick;全部增益来自更好的数学组织和 inductive bias 暴露。
Relation To Prior Work
这篇属于 deterministic least-squares / minimum-energy filtering 与 LQR duality 的谱系,最接近 Willems 的 deterministic Kalman filtering、Sontag 中通过 LQ tracking/time reversal 的推导,以及 MHE 里的 cost-to-arrive 解释。它不是在 stochastic Kalman theory 上新增定理,也不是提出替代 filter;真正不同点是把已有推导中纠缠在一起的步骤拆开:time reversal 只用于理解 final-condition LQR,homogeneous coordinates 只负责 affine-to-quadratic,terminal-state minimization 只负责从 constrained final state 到 filter estimate。
看似新的“two easy steps”本质上是已有思想的重新组合:齐次坐标是优化中的标准 trick,final-condition DRE 是 LQR time reversal,information filter 也是经典形式。实质创新是解释层面的:它把 observer dynamics 的来源定位到扩展 Riccati 的 off-diagonal block,并用 tracking vs filtering 的对照说明 causality 差异。这比传统“stochastic covariance inversion”叙事更适合连接 LQR、MHE 和 deterministic estimation。
Dataset / Evaluation
没有 dataset、实验或数值 evaluation;作为 tutorial note,这是合理的。论文的 claim 不是 empirical performance,而是解析等价与概念澄清,因此 evaluation 应该看推导是否闭合、符号是否一致、是否恢复标准 Kalman-Bucy equations。就这一点而言,文中通过 S^{-1}=Σe、X^{-1}=Σi、V^{-1}=Σv、W^{-1}=Σw 对齐 stochastic Riccati 与 gain,基本支持核心 claim。
但它没有验证任何工程层面的 claim:没有离散时间实现、没有数值稳定性比较、没有高维 sparse information filter 的计算优势展示,也没有 MHE 或 nonlinear extension 的案例。因此不能把这篇看成算法性能论文;它只证明“这个建模视角能导出已知 filter”。
Limitation
主要限制是结构前提很硬:线性动力学、二次能量、正定权重、L2 扰动、无约束连续时间系统。S^{-1} 的存在在推导中被直接使用,但文中未充分说明在半正定 X、不可观方向、奇异噪声权重或有限 horizon 退化情形下如何处理。对于实际滤波,很多关键问题被绕开了:模型误差、非高斯异常值、约束状态、离散采样、数值 conditioning、稳态极限、detectability/stabilizability 条件等。
另一个上限是它没有产生比 Kalman filter 更强的估计器。方法只是把问题“转移”为一个更清晰的 LQR 表述;如果原问题不满足 LQ 闭合性,齐次坐标不会神奇地保留 tractability。推广到 nonlinear minimum-energy filtering 时,value function 不再二次,observer 也不再由有限维 Riccati 精确刻画;此时这篇的 insight 可能仍有建模价值,但算法上需要 approximation。所谓 information / certainty 的解释也依赖 S>0 的 bowl curvature,在退化方向上会变得不完整。
Takeaway
- 1. 最值得记住的是:Kalman filter 可以被看成 final-condition LQR 的 cost-to-arrive center,而不是必须从 stochastic conditional expectation 出发。
- 2. affine-quadratic optimal control 的线性项携带外生信号的记忆;在 initial-condition tracking 中形成 anti-causal feedforward,在 final-condition filtering 中形成 causal observer。
- 这是一个可迁移的结构性 insight。
- 3. information matrix S 的 deterministic 解释不是 covariance inverse 的附庸;它是 cost-to-arrive 对当前状态的曲率,表示偏离最小能量解释的惩罚强度。
一句话总结
这篇论文在 LQR—deterministic filtering 谱系中贡献的是一种干净的结构化推导:用齐次坐标和 final-condition cost-to-arrive 把 Kalman information filter 解释为 affine-quadratic LQR 的 Riccati block dynamics,而不是提出新的滤波算法。
