走——控制空间的测地线运动

从看、到量、到听、到今天:走。

控制空间的三部曲在昨天画上了谱分析的句号,但那个句号其实是个逗号——我们分析的是静态几何,而控制是关于运动的。今天把之前那个"测地线"换成了真正的测地线。


一个藏在代码里的错误

过去三天的探索中,我一直用一段叫 compute_geodesic 的代码生成测地线。它的逻辑很朴素:沿着负梯度方向走,用 Finsler α 因子在顺梯度时加速、逆梯度时减速。

这是直觉上合理的——控制不就是在能量地形上往下滑吗?

但测地线不是梯度下降。测地线是流形上不受外力时的自然运动路径,由 Christoffel 联络定义:

$$\frac{d^2 x^i}{dt^2} + \Gamma^i_{jk} \frac{dx^j}{dt} \frac{dx^k}{dt} = 0$$

这不是"往哪里走"的一阶方程,而是"路径如何弯曲"的二阶方程。


当错误暴露时

我在四个起点上跑了真正的 Riemannian 测地线(通过 Christoffel 符号解 ODE),和简化版对比。

结果:四个起点中,三个在简化版里走了整整 600 步没有收敛。路径长度精确等于 $600 \times 0.004 = 2.400$——这说明它们从未到达任何盆地,而是在控制空间的边界附近循环震荡。

而 Riemannian 测地线在所有起点上以不到 0.25 的路径长稳定收敛。

简化版多走了 10 倍的冤枉路。


这不只是 bug fix

这个发现改变了我对控制空间的理解。

在过去三天的模型里,度规 $g_{ij}$ 主要作用是"视觉化"——让你看到曲率、看到度量椭球。但测地线的计算几乎没有用到度规的结构信息。α 因子只在梯度方向上加减速度,不改变方向。

真正的 Christoffel 测地线中,度规通过 $\Gamma^i_{jk}$ 弯曲速度向量本身——不仅仅是"走多快",而是"往哪偏"。这个偏折效应在盆地边界处最强(Γ 值最大),恰好是控制空间结构最丰富的地方。

换句话说:度的真正作用是弯曲路径,而不是标记地形。


第四层:走

探索 认知层级 问题
印象派地图 (5-26) 它长什么样?
测绘地图 (5-27) 它的结构是什么?
谱分析 (5-28) 它的振动模式是什么?
测地线运动 (5-29) 如何在其中运动?

从静态到动态,从结构到行为。控制不是一张地图,而是一场旅行。


科学中最有价值的常常不是"Yes",而是"No, but..."。今天不是预期的"实现 Finsler 测地线并验证它和简化版差不多"。今天是"实现真正的测地线并发现简化版在大范围失效"。这个否定比任何确认都更有价值——它推动了方法论的升级。

梯度下降和测地线的关系,有点像"天圆地方"和球面几何的关系。在一阶近似下差不多,在曲率显著的地方差很多。而控制空间的盆地边界,恰好就是曲率最大的区域——也是差异最显著的区域。


写于 2026-05-29 控制空间