几何相同而动态相异：Uncontrolled Manifold (UCM) 时间平均分析对运动变异时间结构的盲区

《Experimental Brain Research》：Same geometry, different dynamics: what uncontrolled manifold (UCM) analysis misses about motor control

　　

Uncontrolled Manifold（UCM，无控制流形）分析利用关节运动与任务结果之间的几何关系——任务雅可比（task Jacobian）——将多关节运动变异分为影响任务目标的部分（任务相关）和不影响任务目标的部分（任务无关）。由此产生的分解被解释为神经控制策略的证据，这些策略稳定任务目标，同时允许执行中的灵活性。UCM通常以时间平均方差指数报告，而该时间平均步骤丢弃了运动变异的时间结构。使用三关节平面到达模型，研究人员生成具有三种时间结构（周期性振荡、白噪声和粉红噪声）且在空间协方差上匹配的变异，并发现时间平均几何指数——比值V_UCM/V_ORT和Z变换协同指数ΔV_z——对白噪声和粉红噪声保持统计学不可区分，而在同一试次上计算的时间指数（滞后1自相关、去趋势波动标度、谱斜率）显示出大的效应量。这种双重解离证实时间平均摘要对时间相关性仍然盲。投影操作在单试次和时间分辨成分中保留时间信息，但该信息对运动控制的含义仍然开放：从投影序列的时间属性到特定控制过程的映射本身尚未建立，因此将时间序列方法应用于投影开启了一个新的解释问题。因此，UCM识别协调是否由任务约束组织，同时使生成过程未定。确立这一边界是本工作的贡献：几何分区产生协调的空间解释，其时间起源超出该方法并需要专门的动力学理论。

研究背景与问题：运动变异曾被视为需最小化或平均掉的噪声，后被作为揭示控制原理的信号。Uncontrolled Manifold（UCM，无控制流形）分析成为从运动变异推断协调策略的主导框架。UCM利用任务雅可比（task Jacobian）将关节角、肌肉激活或力等基本变量变异分解为影响任务表现变量的成分和不影响的成分。V_UCM>V_ORT被解释为稳定任务结果而允许执行灵活性的控制组织，并应用于多指力、姿势、行走跑步、指向、到达和抓握。问题在于：UCM通常报告时间平均方差指数，该步骤丢弃运动变异的时间结构。两个探索相同姿态集合的执行者共享相同分区，即使一个慢漂移一个随机抖动。虽有时间分辨UCM，但未处理变异内部时间依赖。三十年来，没有研究系统检验时间平均指数能否区分对应机制不同控制过程的时间结构。理论测试和临床比较依赖几何证据，却未在生成过程层面检验。本研究通过模拟解决：需参数化控制时间动态，同时保持空间方差结构恒定，这在生物系统中不可行。

研究设计与方法（≤250字）：研究人员使用三关节平面臂模型（L₁=1.0、L₂=0.8、L₃=0.6，任意单位），标称关节角θ₀=[π/4,π/6,π/8]^T，目标端效器位置p₀=[1.948,1.330]^T。在10 s、N_t=100、N_trials=200下生成周期性（0.1 Hz）、白噪声和粉红噪声（Voss–McCartney）三种时间结构，并匹配任务空间与UCM空间变异。通过Jacobian的奇异值分解定义UCM零空间和ORT范围空间投影，计算V_UCM(t)、V_ORT(t)、r_UCM(t)。采用自相关、Welch功率谱密度（PSD）、去趋势波动分析（DFA）、端效器方差和交叉相关表征时间结构；逐试次计算几何指数（r_UCM、ΔV、ΔV_z）和时间指数（ρ(1)、α、β），用Kruskal–Wallis、Mann–Whitney U、Bonferroni校正和Cliff’s δ进行非参数统计。无生物样本队列，数据为模拟生成。

研究结果：1）关节空间与任务空间运动学在不同条件间不同但共享标称配置：周期性产生0.1 Hz平滑振荡；白噪声不规则无相关；粉红噪声持续平滑长程相关。端效器方差为周期性2.01×10^-5、白噪声4.03×10^-5、粉红噪声4.67×10^-5。2）白噪声与粉红噪声的时间平均方差比相等：r?_UCM分别为周期性2.86、白噪声8.50、粉红噪声8.72，白噪声与粉红噪声差0.22，即2.6%；V_UCM和V_ORT差小于2%。原因是跨试次在固定t汇聚独立试次，粉红结构是试次内属性，被割裂。3）方差成分自相关在不同条件间不同：周期性V_ORT的ρ(1)=0.98，V_UCM的ρ(1)=0.14；白噪声ρ_UCM(1)=-0.09，ρ_ORT(1)=0.10；粉红噪声ρ_UCM(1)=0.25，ρ_ORT(1)=0.17。该差异未反映在时间平均比值。4）跨试次谱斜率白噪声与粉红噪声差异小：周期性在0.1 Hz有峰；白噪声β_UCM=-3.43，粉红噪声β_UCM=-3.32，Δβ=0.11。但该结果是跨试次聚合，单试次PSD斜率恢复大分离，Cliff’s δ=+0.81，p<10^-44。5）单试次投影保留白噪声与粉红噪声时间区分：UCM成分白噪声α=0.49、粉红噪声α=0.79；ORT成分白噪声α=0.57、粉红噪声α=0.75；ACF白噪声ρ(1)=-0.05、粉红噪声ρ(1)=0.26。投影线性保留试次内变异，折叠发生在时间平均。6）DFA标度指数在投影成分中分离白噪声与粉红噪声：白噪声α_UCM=0.58、α_ORT=0.59；粉红噪声α_UCM=0.87、α_ORT=0.79；周期性α_UCM=0.81、α_ORT=1.81。UCM白噪声与粉红噪声Δα=0.29，约50%增加，对比时间平均比值2.6%差异。7）瞬时比值序列无白噪声与粉红噪声差异：周期性均值2.86振荡；白噪声均值8.58；粉红噪声均值8.83；ACF和PSD相似；V_UCM与V_ORT交叉相关近零（周期性-0.02、白噪声-0.05、粉红噪声-0.19）。8）逐试次指数显示几何对比可忽略、时间对比大：白噪声与粉红噪声r_UCM为8.69±1.66与8.67±0.83，δ=-0.04，p=0.48；ΔV_z为1.23±0.09与1.24±0.05，δ=-0.04；V_UCM的δ=-0.07，V_ORT的δ=+0.02。时间指数：ρ_UCM(1)为0.003±0.088与0.451±0.112，δ=-1.00，p=1.0×10^-66；α_UCM为0.60±0.09与0.88±0.12，δ=-0.92，p=7.0×10^-57；β的δ=+0.81，p=6.8×10^-45。ORT类似。该模式构成效应量空间中的双重解离。

讨论总结：标准UCM摘要不能区分仅时间不同的运动。时间平均比值和协同指数在根本不同时间结构上取等效值。逐试次检验给出白噪声与粉红噪声的关键双重解离：几何指数效应可忽略，而自相关和DFA效应大。ΔV_z是比值的单调变换，限制扩展到所有将分区沿时间平均形成的标量。该不变性是数学性的：时间平均将分区折叠为单一数字，移除时间信息。两个系统可有相同时间平均比值，却反映主动反馈、被动相关噪声积累等不同过程。该局限早已被认识，但此前未量化时间平均比值对不同时间结构的不变性。本研究表明，白噪声、粉红噪声和周期动力学比值仅差2.6%，而标度指数差约50%。这与Campolo et al. 2013的Riemannian坐标敏感性解决方案正交：后者去除坐标几何伪迹，本工作表明几何正确的投影仍丢弃时间组织。模拟直接规定时间结构以隔离聚合步骤，说明时间平均指数丢弃时间动态，无论其如何产生。主动控制与被动力学均可产生时间结构，时间与几何方法互补。高V_UCM/V_ORT与 Equilibrium Point Hypothesis（EPH，平衡点假说）的弹簧样松弛、最优反馈控制、最小干预策略或被动约束均一致，不能区分。多分形标度和临床比较也受相同限制。简化平面臂模型限制外推，但时间平均折叠分区、投影传递时间排序的结构性质不依赖复杂度。逐试次时间指数应优先于跨试次系列。

研究结论翻译：Uncontrolled Manifold（UCM）分析通过将基本变量偏差投影到Jacobian定义的子空间，将运动变异分为任务相关和任务无关成分。研究人员的核心结果涉及通过时间平均摘要报告该分区——比值r?_UCM和协同指数ΔV_z。这些摘要对时间结构保持不敏感：白噪声和粉红噪声在自相关和标度上相差数量级，但产生统计学不可区分的值。不敏感性源于时间平均步骤；投影将时间信息传递到V_UCM(t)、V_ORT(t)和单试次投影，其中信息对控制的含义仍为开放问题。时间平均比值区分确定性机制与随机机制，同时对共享瞬时协方差的随机过程赋予匹配值。因此，高比值确定变异性沿任务相关方向组织——这是真实且有用的发现——但使生成过程——主动控制、被动力学或生物力学滤波——未定。时间平均丢弃的时间信息存活于单试次和时间分辨投影中，这种存活标志着该方法未解决的边界。自相关、去趋势波动、谱和多分形分析可应用于投影成分，产生指数和相关结构；这些量对控制的含义是另一个问题，因为投影序列继承Jacobian和任务变量的几何，从它们的时间行为读取控制需要理论和实证论证。UCM描述变异性相对于任务需求的位置。随时间生成该变异性的过程超出其范围，定位它们需要其解释已被独立建立的方法。这是运动控制的几何、时间平均解释的持续局限：它精确映射变异性的空间组织，同时将这种组织的动力学起源留给尚待建立的理论。