精读笔记

Problem Setting

论文标题：Continuous Time Identification of Linear Systems: Extended Version（arXiv preprint / 2026-06-10）。

这篇论文不是在解决一般 continuous-time identification，而是在一个更强约束下做 SISO 稳定 LTI 辨识：计算过程必须尽量连续时间、递归、低存储，不能依赖离散采样数据矩阵。已知阶数下，Kreisselmeier/Narendra-Annaswamy 类型 adaptive observer 已经基本解决参数估计；真正卡住的是阶数未知。

关键矛盾是：adaptive observer 需要一个阶数匹配的参数化才能让 prediction error 对参数线性、让 PE regressor 成立；但阶数本身又需要从输入输出数据中学习。传统 subspace/Ho-Kalman 通过 Hankel/rank 从采样数据里恢复阶数，这和本文的 neuromorphic 约束不兼容。本文要找的是一个连续时间对象，其几何维数直接编码 plant order。

Motivation

已有路线的缺口很明确：连续时间辨识文献通常仍以 sampled input-output data 为入口；而 adaptive observer 虽然连续时间，但默认 plant order 已知。作者的核心观察是：如果能构造一个阶数可能大于 plant 的输入输出等价模型，那么未知阶数可以不再作为 observer 结构的先验，而作为该过参数化模型内部不可控/可控分解的结果出现。

这也是为什么论文会回到 Kreisselmeier-Lozano 的 overmodeled adaptive control。那里的 overparameterization 主要用于处理 adaptive control 中的 overmodeling；本文把它重新解释为 system identification 中的阶数探测工具。缺的不是另一个参数估计算法，而是一个在阶数不匹配时仍有稳定几何含义的状态空间参数化。

Core Idea

核心思想可以压缩成一句：不要直接估计 plant order，而是让一个固定阶数 \bar n 的连续时间滤波器生成一个过参数化 realization，然后通过其可控子空间维数读出真实阶数。

当 \bar n\ge n 时，高阶模型可以通过 \alpha(s)\delta(s), \beta(s)\delta(s) 与 plant 输入输出等价，多余阶数变成 cancellation。关键不是这个事实本身，而是所有 admissible \theta 诱导同一个 controllable subspace，且 dim(C)=n+\bar n。于是 plant order 变成 rank(continuous-time covariance) - \bar n。这个转换很漂亮：它把辨识问题从“找到正确模型结构”变成“估计一个由 excitation 暴露出来的 invariant subspace”。

和 prior 的本质差异在于信息流重组：传统方法先采样、存数据、构造 Hankel 或回归矩阵，再做 rank；本文先用稳定连续滤波器把输入输出嵌入状态，再用时间演化和 PE 让子空间从协方差流中浮现。inductive bias 是强烈的 LTI/SISO/稳定/可激励结构，而不是数据驱动泛化。

Method

1. 过参数化 observer：对 y 和 u 分别通过同一个稳定可控滤波器 (F,G)，得到 \hat\xi_1,\hat\xi_2。这一步解决的是阶数未知时无法直接构造 plant observer 的问题。核心变化是，真实 plant 被表示为滤波状态上的静态输出映射，而不是显式 plant state。

2. overmodeled 几何定理：若 \bar n\ge n，则存在 \theta 使模型输入输出等价，并且可控子空间维数为 n+\bar n。它解决的是阶数如何从过高维模型中被区分出来的问题。真正重要的是 uncontrollable modes 正好对应多余 cancellation polynomial \delta(s)。

3. PE 子空间估计：通过 sufficiently rich / multisine 输入，让特殊解 \xi_\star(t) 的 PE 子空间等于可控子空间。连续时间滤波 \dot\Sigma=-\epsilon\Sigma+\epsilon\xi\xi^T 相当于递归估计 autocovariance 的像空间。它解决的是无法直接访问 C 的问题。

4. 欠建模判别：当 \bar n<n 时不存在 overparameterized equivalent model，作者改看 observer 状态本身；随机多正弦输入下，几乎处处 rank(R[\hat\xi])=2\bar n。于是欠建模时估计阶数等于 \bar n，不会出现 plateau。

5. 递增阶数算法：从 \bar n=1 开始重复实验，直到 \hat n(\bar n)=\hat n(\bar n-1) 出现固定点。这个过程本质是 continuous-time 版本的 rank-growth test，形式上类似 Ho-Kalman 递增 Hankel rank。

6. 参数识别：阶数确定后令 \bar n=n，使用标准 gradient law \dot{\hat\theta}=-\gamma e\hat\xi。这里贡献主要是把 PE 论证建立在 overparameterized special solution 上，而不是提出新型参数估计器。

Key Insight / Why It Works

最核心的技术贡献是 Theorem 6，而不是算法本身。它说明 plant order 被编码在一个输入输出等价但非最小 realization 的可控维数中：dim(C)=n+\bar n。这个结果之所以成立，是因为过建模引入的 \delta(s) 只产生 uncontrollable modes，而真实 plant poles 和滤波器 f(s) 的 modes 留在 controllable part。换言之，多余阶数不是被噪声般混在状态里，而是结构性地落入不可控子空间。

第二个关键 insight 是 PE 子空间 = controllable subspace。对于稳定 LTI，足够丰富输入会让可控方向在状态协方差中显现，而不可控方向不会由输入激发。这个机制把 order estimation 转成 rank estimation。它本质上不是 scaling，也不是数据覆盖意义上的泛化，而是强模型结构下的 latent subspace recovery。

欠建模部分的作用是让递增算法闭环：如果 \bar n<n，rank 满为 2\bar n，于是不会误停；一旦 \bar n>n，rank 缺口出现，plateau 暴露真实 n。这类似 curriculum / incremental model growth，但不是学习 curriculum 带来的泛化，而是 rank signature 的离散搜索。

辅助成分包括 SVD、阈值、dimension plot、多正弦输入选择。这些更多是 engineering / numerical rank recovery。尤其噪声下的 \hat\sigma_{tol}\approx0.01\lambda_{min}(\hat\Sigma_2) 明显是 heuristic，增益来源不清；理论 rank gap 在有限时间、有限 SNR 下是否可靠，文中未充分说明。

Relation To Prior Work

最接近的谱系有三条：Kreisselmeier adaptive observer / adaptive control；Ho-Kalman/subspace order estimation；continuous-time sampled-data identification。本文不是提出新的最小 realization 算法，而是把 Kreisselmeier-Lozano overparameterized model 的老思想重新用于阶数识别。

和 Ho-Kalman/subspace 方法的本质差异：那些方法的 rank 对象是离散采样构造的 Hankel/data matrix；本文的 rank 对象是连续时间滤波状态的 autocovariance / PE subspace。形式上都靠 rank plateau，信息来源也都是输入输出行为，但计算 substrate 不同。

和传统 adaptive observer 的差异：传统方法要求 order matched，本文允许 observer order 与 plant order mismatch，并给出 overmodeled/undermodeled 两种 rank 行为。这一点是实质创新。

和 continuous-time identification 文献的差异：许多方法输出 continuous-time model，但仍从 sampled data 开始；本文严格强调 processing in continuous time。不过 SVD 和实验递增仍是离散 orchestration，因此“完全 neuromorphic”这一 claim 需要谨慎看待。

Dataset / Evaluation

evaluation 主要是数值仿真：低阶三阶系统、高阶八阶系统、输出测量噪声下的 Monte Carlo singular value gap。任务覆盖范围窄，基本都在作者理论假设内：稳定 SISO LTI、可设计多正弦输入、仿真可控噪声。

这些实验足以说明 dimension plot 机制在干净设置下能工作，也能显示噪声会压缩 singular value gap。但它们没有真正验证 neuromorphic claim：没有真实神经数据、模拟硬件、真实连续传感器系统或闭环 motor task。也没有和 sampled-data continuous-time identification、RLS mismatch 方法、subspace rank estimator 做系统比较。

因此 evaluation 支持的是“理论机制可数值实现”，而不是“这是实际更鲁棒/更 scalable 的辨识框架”。benchmark 没有明显 leakage 问题，但场景选择高度贴合理论假设，外推性有限。

Limitation

1. 假设很强：SISO、LTI、稳定、严格真/相对阶隐含在给定传函形式中、coprime、输入可自由设计且 sufficiently rich。任何一个放松都会影响可控子空间维数公式或 PE 论证。

2. rank estimation 是真正脆弱点。理论上 rank 有清晰 gap；实际中 finite time、噪声、滤波器 f(s) 选择、输入频带、plant attenuation 都会改变 singular value separation。论文的阈值策略有经验性，文中未充分说明如何给出可靠有限样本误差界。

3. scalability 不是免费的。协方差矩阵维度是 2\bar n×2\bar n，SVD 成本随阶数增长；高阶滤波状态对输入频率和幅值敏感，实验中八阶系统已经需要放大输入以克服高阶衰减。这说明方法可能更适合 low-order model learning，而非高阶复杂系统。

4. “连续时间”并不等价于完全无离散计算。算法需要逐次实验、终止判据、SVD/rank threshold，这些是离散元控制。若作为 brain-level model，这部分生物可实现性仍偏概念化。

5. 参数识别部分不是主要创新。阶数确定后基本回到标准 adaptive observer + gradient law；若输入不 PE 或只能自然激励，参数收敛保证会弱很多。

6. MIMO/非线性扩展不是平滑推广。SISO polynomial cancellation 给了非常干净的 \delta(s) 结构；MIMO 下 Smith-McMillan、非最小 realization、方向性零点和输入激励子空间都会让 dim law 复杂得多。

Takeaway

• 1. 最值得记住的是：阶数识别可以通过“过参数化输入输出等价模型的可控子空间维数”来做，而不是必须依赖 Hankel 数据矩阵。
• 2. PE subspace 是一个可迁移 insight：如果某个 latent subspace 是系统结构真正承载信息的地方，足够丰富输入 + 连续时间协方差流可以把它显现出来。
• 3. 这篇真正推动的是 adaptive observer 与 subspace identification 的接口：用 observer 生成状态嵌入，用 subspace/rank 读出结构，再回到 adaptive law 估参数。
• 4. 未来最值得做的不是再调仿真，而是给噪声/有限时间 rank 判定、输入设计、MIMO 情形的维数公式，以及闭环条件下的可辨识性建立硬保证。

一句话总结

这篇论文把 Kreisselmeier 式连续时间 adaptive observer 改造成一种基于过参数化可控子空间 rank signature 的阶数识别框架，是 adaptive identification 向结构/子空间几何辨识融合的一步，而非单纯的新参数估计算法。