精读笔记

Problem Setting

这篇论文实际处理的是一个很具体但工程上重要的问题：给定 IREM 模型

y_k = \sum_{\ell=-\infty}^{k} g_{k-\ell} u_\ell + \phi(x_k), x_k=x_{k-1}+f(u_k),

在当前时刻只观测到输出子集 \tilde y，而 MPC 需要预测未来全输出 y，包括未测内部变量。困难点有两个：一是卷积模型天然依赖无限过去，过去未知输入/扰动会以“状态”的形式残留；二是 equilibrium state x_0 进入非线性 \phi，而它也未必直接测得。

已有 state-space observer/Kalman/MHE 可以处理类似问题，但前提通常是有一个合适的低阶状态空间模型；而 IREM 的吸引力恰恰在于它可以从少量阶跃/脉冲仿真得到 g_k、\phi、f，不想再做 realization 或 reduction。因此关键矛盾是：模型表示是输入输出/冲激响应形式，但估计和预测需要某种状态可观测性。本文试图直接在冲激响应矩阵上给出这种可观测性和误差放大界。

Motivation

作者不是要提出一个新的 MHE algorithm，而是补 IREM 这类模型结构缺失的理论环节：如果用过去窗口拟合 measured outputs，然后拿估计出的 equivalent input 和 equilibrium state 去预测未来全输出，拟合得好是否意味着预测一定好？以及这个“好”由什么结构量决定？

传统路线的问题在这里很明确： 1. Kalman/filtering 路线需要 state-space realization，尤其电池高阶物理模型往往还要模型降阶； 2. generic nonlinear MHE 理论给的是抽象 observability/detectability 条件，不容易直接从实验得到的冲激响应判断； 3. 神经网络/Koopman 等数据路线虽然更 general，但对本文目标场景过重：作者假设系统类别较窄、模型结构较强、可通过少量仿真实验得到核心模型量。

核心观察是：对于 mildly nonlinear systems，很多复杂动态可以拆成“围绕当前平衡点的线性动态记忆”加“缓慢/积分驱动的 equilibrium shift”。只要线性卷积部分和 equilibrium 非线性在过去窗口中可分离，那么过去拟合误差就可以控制未来预测误差。

Core Idea

论文真正的核心是把状态估计问题改写成一个子空间几何问题。未知无限过去不显式建模，而是用一个有限长度 equivalent input U^e 表示它对当前及未来的影响；这相当于 Nerode equivalence 在冲激响应表示下的有限参数化。线性情形中，若所有在 measured past 中不可见的 equivalent input 方向也不会影响 future full outputs，即 Null(\tilde T_1^{p,r}) ⊆ Null(T_{p+1}^{p_f,r})，那么任何过去拟合误差都会通过一个可计算增益 σ_max(T_future \tilde T_past^+) 放大到未来误差。

IREM 的新增难点是 equilibrium state x_0 也未知，并且通过 \phi 非线性进入过去和未来输出。作者的做法不是构造 nonlinear observer，而是把 \Phi_p(\hat x_0)-\Phi_p(x_0) 写成 \Psi_p(\hat x_0,x_0,U_p)(\hat x_0-x_0)，然后要求这部分变化不能落入 equivalent input 的过去卷积子空间。也就是说，equilibrium state 的变化必须在 measured past 中留下无法被未知历史输入伪装的 signature。这个 signature 的强度就是 α。

本质区别在于：prior observer 理论通常从状态空间 dynamics 出发；这里从冲激响应 Toeplitz 子空间出发，直接刻画“历史记忆”和“平衡态漂移”的可分离性。这是一个更贴合 IREM 数据获取方式的 inductive bias。

Method

关键机制可以压缩为以下几项。

1. Finite equivalent input：解决无限历史不可用的问题。对于 n 维 LTI 卷积系统，长度 r≥n 的 U^e 可以代表 k≤-p 的无限过去对未来的全部影响。这一步把状态估计变成代数变量估计，而不是显式恢复 state-space state。

2. Past-fit to future-error gain：解决“拟合过去是否能保证预测未来”的问题。线性情形下，过去残差为 \tilde T_1^{p,r}(\hat U^e-U^e)，未来误差为 T_{p+1}^{p_f,r}(\hat U^e-U^e)。如果过去矩阵的 null space 被未来矩阵的 null space 包含，则未来误差可由过去残差通过 σ_max(T_future \tilde T_past^+) 控制。这个条件本质上是 partial-output observability 的冲激响应版本。

3. Nonlinear equilibrium decomposition：解决 x_0 未知的问题。作者把 nonlinear equilibrium contribution 的差分写成 \Psi_p(x,x_0,U_p)(x-x_0)，再把 \Psi_p 投影到过去卷积子空间正交补。若最小奇异值 α>0，则 equilibrium error 无法被 equivalent input error 完全吸收，因此过去拟合残差能同时约束 U^e error 和 x_0 error。

4. α>0 的充分条件：解决可验证性问题。论文给出两类条件：sector-bounded nonlinearity 接近一个 full-column-rank 线性映射 B；或 L-smooth nonlinearity 且 Jacobian 序列 band-limited，并与线性冲激响应模态在频谱上分离。技术上通过冲激响应分解为 k^l λ^k 模态，并证明不同 λ 模态生成的矩阵列独立。

5. No integrator in impulse-response part：通过要求 \tilde T_2^{n,n}(\tilde T_1^{n,n})^+ 没有 λ=1 或远离 λ=1，确保积分行为由 equilibrium state x_k 承担，而不是隐藏在线性卷积 g_k 中。否则 equilibrium drift 和 linear memory 的模态会混叠，可观测性会退化。

Key Insight / Why It Works

最重要的 insight 是：IREM 的可观测性不是一般 nonlinear observability，而是两个子空间之间的横向性。

第一层横向性是 linear memory：unknown past equivalent input 的某些方向如果在过去 measured outputs 中不可见，那么它们也必须对未来 full outputs 无影响。否则过去窗口无论怎么拟合都无法排除未来分叉。这正是 Null(\tilde T_past) ⊆ Null(T_future)。这个条件很干净，因为它完全从 impulse response 构造，不需要 realization。

第二层横向性是 equilibrium state：\phi(x) 的变化在过去窗口中必须有一部分不能被 linear convolution history 解释。α = σ_min((I-P_{col \tilde T})\Psi_p) 的含义非常直接：把 nonlinear equilibrium variation 投影到“不是未知历史能产生的输出模式”的方向上，看最弱方向还剩多少。如果 α 小，x_0 error 可以被 U^e error 伪装，过去拟合再好也不可靠；如果 α>0，过去残差同时限制 x_0 和 U^e，从而未来误差可控。

这篇最核心的贡献不是 bound 公式本身，而是把 IREM 的 observability 具体化成 impulse-response Toeplitz 子空间 + equilibrium nonlinear signature 的 separation。这个判断可以迁移：任何把慢变量/平衡变量和快线性动态叠加的模型，都应该检查慢变量 signature 是否落入快动态 memory 子空间。

哪些部分可能只是辅助？sector-bounded 与 band-limited Jacobian 条件更像证明 α>0 的 sufficient conditions，不一定是实际使用中最自然的诊断。实际工程里可能直接数值计算 α 或局部 α 更有意义。Theorem 6 这套 k^l λ^k 独立性证明是技术支撑，但不是核心建模 insight。

这不是 scaling，不是 data coverage，不是 retrieval，也不是 test-time compute。它本质上是 better inductive bias + latent structure separation：把系统记忆拆成 finite equivalent input，把 operating-point drift 拆成 integrator equilibrium state，并用子空间几何判断两者是否可分离。

Relation To Prior Work

最接近的是 Vincent et al. 2023 的 impulse-response prediction using equivalent input。本文在线性部分基本继承那条路线：用 finite equivalent input 表示 Nerode equivalence class，并通过 Toeplitz 矩阵做预测。新增的是两点：一是给出 past-fit 到 future-error 的确定性增益界；二是把该框架扩展到 IREM，即增加 nonlinear equilibrium/integrator state 并给出可观测性条件。

和 Dynamic Matrix Control 的关系也值得区分。DMC 同样使用冲激/阶跃响应做预测，但通常用有限脉冲响应近似和常值扰动校正，且默认输出可测。本文不是简单 DMC，而是试图在无限历史精确压缩和部分输出可测下恢复类似状态估计的功能。

和 classical observer/Kalman/MHE 的本质差异在表示层：传统方法先有 state-space model，再谈 observability；本文直接在 impulse response 上谈 observability。Theorem 3 说明 null-space 条件与 state-space observability 有联系，但实际判定不需要显式 realization。

和 generic nonlinear MHE stability literature 的区别是目标不同。MHE 文献常给扰动/噪声条件下的收敛稳定性；本文给的是一个 deterministic algebraic gain：过去拟合误差乘以模型相关增益控制未来预测误差。它更局部、更结构化，也更容易针对 IREM 检查，但泛化性明显更窄。

看似新但其实是重组的部分：pseudo-inverse bound、Hankel/Toeplitz rank、Ho-Kalman 相关模态分解都不是新思想。实质创新是把这些线性系统工具组织到 IREM 的“linear convolution + nonlinear equilibrium”结构上，并指出 λ=1 模态必须从 g_k 中剥离给 integrator equilibrium state。

Dataset / Evaluation

本文没有常规意义上的 dataset/evaluation，也没有实验 benchmark。它是理论预印本，evaluation 主要是定理条件是否可从模型参数计算、bound 是否成立。

这意味着它并没有实证验证几个工程 claim：例如在电池 fast charging 场景中，α 的实际大小是否足够好；sector-bounded 或 band-limited Jacobian 条件是否在真实模型上成立；预测误差界是否紧；在测量噪声、模型误差、gain scheduling、冲激响应截断下是否仍然有可用 performance。

从支持核心 claim 的角度看，理论部分足以支持“存在可计算的 deterministic bound under assumptions”，但不足以支持“这个 bound 在实际 MPC/MHE 中有效”或“IREM 比 state-space observer 更实用”。后者仍需要数值实验或真实电池案例。

Limitation

最根本限制是结构假设强。IREM 假设系统可以拆成固定 linear convolution g_k 加 equilibrium function \phi(x_k)，且 x_k 只有 integrator dynamics。对 mild nonlinearity 和单一 operating-point coordinate 很合理，但对强非线性、路径依赖、hysteresis、多时间尺度耦合、输入相关 dynamics 的系统，这个拆分会漏掉主要误差。

第二，理论把 g_k、f、\phi 当作已知且精确。现实里 g_k 来自有限仿真/实验，通常还会随 state-of-charge 或温度调度。文中提到实践可 schedule g_k，但分析只覆盖 unscheduled model。这个 gap 很大：gain scheduling 会让 Toeplitz 矩阵随估计状态变化，null-space 条件和 α 都变成轨迹相关甚至估计相关。

第三，α 是理论上最关键但也最脆弱的量。α>0 是 qualitative observability；实际性能取决于 α 是否远离 0。若 \phi 的某些方向很平坦，或其过去窗口 signature 与 convolution memory 高度对齐，bound 会爆炸。文中未充分说明如何选择输入窗口 U_p 来增强 α，也未给出可操作的 experiment design。

第四，sector-bounded 条件偏强，要求 \phi 在全域内接近一个 full-column-rank B；L-smooth 条件虽然更自然，但需要 Jacobian sequence band-limited 且与 impulse-response eigenvalues 分离。这更像证明技术条件，而非一般 nonlinear system 中容易验证的性质。

第五，误差界没有显式处理噪声、扰动、模型失配。虽然结论可以把测量噪声当作 past-fit residual 传递，但这不是 stochastic filtering theory；没有给出 estimator design、regularization、bias-variance、constraint handling，也没有讨论 U^e 和 x_0 在噪声下的数值病态。

第六，bound 可能保守。它用 worst-case operator norm、pseudo-inverse 和 infimum singular value，实际预测误差可能远小于上界。增益来源清楚但实用 tightness 文中未说明。

Takeaway

• 1. IREM 的关键不是“用冲激响应加非线性函数”这么简单，而是要保证 linear memory 和 equilibrium drift 在过去测量中可分离；否则 estimation 会把 operating-point error 和 history error 混在一起。
• 2. 对输入输出模型做预测时，最有用的 observability 语言可能不是 state-space rank condition，而是 Toeplitz 子空间的 null-space inclusion 和 cross-subspace gain。
• 3. λ=1 模态的归属很关键：如果 linear impulse response 里藏着 integrator，而 equilibrium state 也在积分，模型结构会不可辨识。
• 把积分动态显式放进 equilibrium variable 是 IREM 成立的核心 inductive bias。

一句话总结

这篇论文把 IREM 从一个工程上好用的冲激响应/平衡态模型，推进为一个有确定性预测误差界和冲激响应级可观测性判据的结构化估计框架。