模块定位：步态、骑行、划船、举重等重复性运动，每次循环的曲线都长得差不多却又不完全一样。本节回答三个工程问题：怎样把许多次循环合成一条”代表曲线”（集合平均）？每人每步时长不同，怎么对齐（时间归一化到 100%）？合成之后如何用一个数字说明”这个人动作稳不稳”（波形变异系数 CV）？——这是从”单次数据”走向”个体画像”和”人群参照”的关键一步，直接决定本系统功能表现型的可比性。
0. 大白话定义
你走了十步，每步的膝关节曲线都略有不同。哪条才算”你的”步态？答案是：都算上——把这十条曲线按步态周期对齐，逐点求平均，得到一条集合平均曲线；再把每条曲线偏离平均的程度平均一下，得到一个变异系数。前者是你的”典型动作”，后者是你的”动作稳定性”。
但十步的时长各不相同，直接对秒数逐点平均会对错位置。所以先把每条曲线的时间轴从”秒”拉伸/压缩成统一的 0–100% 周期（线性插值重采样），再逐点平均。代价是：真实时间信息被丢掉了——凡是要算”加载速率""总负荷时长”的指标，绝不能在这条归一化曲线上算。
临床上最有力量的用法：把 20 个同龄健康人的平均曲线 ±1 标准差画成一条”正常走廊”，把患者的曲线叠上去——曲线一出走廊，异常即刻可见。
1. Definition 定义（术语表）
| 术语 | 英文 | 定义 |
|---|---|---|
| 集合平均 | Ensemble averaging | 将重复性运动的多次循环（个体内）或多个受试者（个体间）的同变量曲线，按周期对齐后逐点平均，得到代表性波形 |
| 个体内/个体间平均 | Intra- / inter-subject averaging | 前者对同一受试者多次循环平均；后者对一群受试者（每人先平均）再平均 |
| 时间归一化 | Time normalization (to 100%) | 将各循环的时间轴从实际秒数线性插值重采样为统一的 0–100% 周期，消除时长差异以便平均 |
| 线性插值 | Linear interpolation | 重采样时按相邻两原始采样点的距离比例加权求值的方法 |
| 变异系数 | Coefficient of variation (CV) | 标准差与均值之比 σ/X̄；单点测量用样本标准差，波形用式 2.20（见下） |
| 波形变异系数 | Waveform CV (式 2.20) | 全周期 N 个点的平均标准差 ÷ 平均信号绝对值，度量平均波形周围的平均离散度 |
| 正常参照带 | Normative band | 健康人群集合平均 ± 1 标准差构成的带状区域，患者曲线叠加其上做临床对照 |
| 协方差分析 | Covariance analysis | 由各关节集合平均的方差求关节间协方差，分析多关节协同变异（支撑力矩相关，见第 11 章） |
| 倒摆模型 | Inverted pendulum model | 安静站立时将全身简化为单一质心绕踝关节摆动的模型（Gage et al., 2004），用于 COM-COP 比较 |
2. 核心公式与流程
2.1 集合平均的价值与依据
- 为什么平均：平均波形更可靠；平均周围的变异给出变量的”随机性”画像。
- 一个关键现象（步态）：个体内下肢关节角度变异很小，而同一关节的力矩曲线变异相当大——正是这一现象催生了协方差分析与支撑力矩研究（详见第 11 章）。这与 M16 的变异性层级（EMG > 力矩 > 运动学）互为印证。
- 临床用法：集合平均波形 + 正常人群参照带是临床评估的基础工具，患者曲线叠加于同龄健康人群平均之上，可直观定位具体运动异常。
2.2 时间归一化到 100%（2.3.2，含插值公式）
设某变量在一个周期内有 n 个采样点，要归一化为 N 个点（通常 N = 100）。归一化曲线上每个间隔对应原曲线的 n/N 个采样。
算例（书中原例）：n = 107，N = 100。
- 第 1 个归一化点对应原曲线第 1.07 个采样 → y₁ = 0.93x₁ + 0.07x₂
- 第 2 点对应第 2.14 个采样 → y₂ = 0.86x₂ + 0.14x₃
- 第 3 点对应第 3.21 个采样 → y₃ = 0.79x₃ + 0.21x₄
- ……
- 第 99 点对应第 105.93 个采样 → y₉₉ = 0.07x₁₀₅ + 0.93x₁₀₆
- 第 100 点 → y₁₀₀ = x₁₀₇
即：归一化点 k 对应原始位置 k·(n/N)，取整得左邻点，小数部分为右邻点权重。
⚠️ 临床警告（本节最重要的一句话）：归一化把时间（秒）从生理数据中抹掉了，任何以时间为因子的指标都不得在时间归一化曲线上计算——下肢承重关节的加载速率与总负荷时长被认为影响生物组织完整性，这类计算必须用真实时间轴的原始曲线。
2.3 波形变异系数（2.3.3，式 2.20）
(1/N)·Σᵢ σᵢ
CV = ──────────────
(1/N)·Σᵢ |Xᵢ|
- 分子：全周期 N 个点上各点标准差的平均（即平均离散度）；
- 分母：全周期信号绝对值的平均（即用整条曲线的”典型幅度”做归一）；
- 有研究者建议分母改用信号的 rms 值，但对 CV 分数影响很小。
- 解读：CV 越小，重复性/稳定性越高；可与单点测量的 CV = σ/X̄ 类比记忆。
3. Physical Meaning 物理意义
- 平均 = 提高信噪比的时域手段：若各次循环的随机扰动互不相关，平均 m 条曲线可使随机分量按 √m 缩小——集合平均与滤波（M17）是同一目标的两种实现：一个利用”重复”，一个利用”频率分离”。
- 归一化 = 坐标系变换：把”秒”换成”周期百分比”是一次非均匀的时间轴缩放（每人缩放率不同）。它保留了形状与时相比例（事件发生在周期的哪个阶段），牺牲了绝对时间与速率量纲。
- CV 是无量纲稳定性指标：用信号自身幅度做分母，使不同幅度、不同单位、不同关节之间的变异性可横向比较。
4. Human Movement Meaning 人体运动意义
- “僵直腿步态”的识别范式（Figure 2.20）：一位髋骨关节炎患者的矢状面髋关节角度（虚线）叠加于 20 名同龄健康人平均（实线）±1 SD（灰带）之上：支撑初期髋屈曲不足、支撑末期髋伸展缺失，整体关节活动幅度变小——“stiff-leg gait”。这一范例确立了临床读图顺序：看形态 → 看出带位置（周期的哪一段）→ 看方向（多还是少）→ 命名异常模式。
- 角度稳、力矩不稳：步态中运动学高重复性而力矩高变异，说明 CNS 允许多种力矩组合实现”同一个动作”（呼应 M16 的冗余与四级汇聚）——只测角度会高估一个人动作的”稳定性”，只测力矩会低估它。
- 归一化后的比较才是公平的：不同身高、步频者的支撑期长短不同，不归一化就把”第 0.3 秒”对在一起比较是错位比较；归一化把比较单位统一为”动作周期的阶段”。
5. How It Is Measured 如何实现
- 分割周期：按事件（足跟着地→下一次足跟着地，或支撑期起止）把连续数据切成单个循环；事件检测质量决定对齐质量。
- 逐条归一化：每条循环线性插值到 100 点（或 101 点含两端）。
- 逐点统计：对每个百分比位置求均值与标准差 → 平均曲线 + ±SD 带。
- 算 CV（式 2.20）：个体内给”动作稳定性”，群体给”人群离散度”。
- 临床出图：健康人群平均 ± 1 SD 带 + 患者曲线叠加（Figure 2.20 范式）。
- 红线：加载速率、冲量、负荷时长等指标回到未归一化的原始时间轴上计算。
6. Typical Applications 典型应用
| 场景 | 用法 |
|---|---|
| 临床步态报告 | 患者曲线 vs 同龄正常带，逐关节逐平面出图（髋 OA 僵直腿步态范例） |
| 康复进展监测 | 同一患者不同周次的集合平均曲线与 CV 追踪 |
| 跑步/骑行/划船/举重技术分析 | 周期性动作的代表波形与稳定性 |
| 人群参照数据库 | 按年龄/性别/速度分层的 normative band（本系统 Metric Dictionary 的参照形态） |
| 多关节协同研究 | 集合平均 → 方差 → 协方差 → 支撑力矩变异结构（第 11 章） |
7. Clinical Relevance 临床相关性
- 参照带诊断法的力量与边界：叠加图能定位”哪个关节、周期哪一段、偏多还是偏少”，但它仍是 M16_测量描述分析评估分层 中的描述层工具——看出异常 ≠ 知道原因（需动力学/EMG 分析层跟进）。
- 匹配对照人群是硬要求：范例使用 20 名年龄匹配成人——年龄、性别、步速不匹配，参照带本身就是错的。
- 归一化的临床陷阱：对组织负荷类问题（应力性损伤、骨关节炎负荷管理，衔接 M15_过用损伤与应力），加载速率与负荷时长恰恰是核心指标——在归一化曲线上算速率会系统性失真，属于方法学错误。
- CV 作为康复终点指标：动作稳定性（CV 下降）可作为”功能改善”的量化证据，但仍需 M16 的警告——CV 变化只记录变化，不证明治疗因果。
8. Performance Relevance 表现相关性
- 精英运动员的 CV 通常更低：高重复性是熟练技术的标志，CV 可作为技术稳定性与疲劳/状态波动的监控指标。
- 个体间平均用于”理想模型”要谨慎：把群体平均当技术模板可能抹掉个体化最优解——运动员的评估目标是自身纵向比较（M16），集合平均的正确用法是”自己 vs 自己”。
- 非周期动作不适用：投掷、跳远起跳等单次爆发动作无循环可平均，需改用多次试次的单次曲线比较或功能数据方法 [补充说明]。
9. Possible Interpretations 可能的解释层级
| 观察 | Level B（描述层） | Level C（机制层解释） | Level D（常见误读） |
|---|---|---|---|
| 患者髋角曲线在支撑末期低于正常带 | 髋伸展不足（形态学事实） | 可能为疼痛回避策略、关节囊挛缩或伸髋肌力不足——需力矩/EMG/查体区分 | ”角度异常=该关节有病变” |
| 个体内力矩 CV 高而角度 CV 低 | 力矩曲线重复性差、角度重复性好 | 正常冗余：CNS 以不同力矩组合实现同一运动学输出 | ”力矩不稳=控制能力差”（健康人也如此） |
| 归一化曲线上患者加载峰提前 | 峰值出现在周期更早阶段 | 需回到真实时间轴确认速率是否真改变 | 在归一化曲线上直接算”加载速率”并下结论 |
| 康复后 CV 从 12% 降到 8% | 动作重复性提高 | 可能反映神经肌肉控制改善，也可能是速度变慢等混杂 | ”CV 下降=治疗有效”（监测≠因果） |
10. What It Cannot Tell Us 这套方法不能告诉我们什么
- 平均抹掉单次特征：偶然出现但有诊断价值的单次异常（如某一步的失足）在平均中消失；平均曲线应与原始曲线叠加检查。
- 归一化丢弃绝对时间：速率、冲量、负荷时长一律失效（本节红线）。
- 参照带不解释机制：出带是”什么/哪里”，不是”为什么”。
- CV 不区分变异来源：生理性冗余、测量噪声、病理不稳都计入同一个 CV；需结合 M17 的噪声评估手段分解。
- 非周期动作不适用：无重复周期则无集合平均可言。
11. Related Metrics 相关概念与双链
- M17_相关分析与频谱：相关/频谱要求有意义的数据段，集合平均解决”多段重复数据如何合成代表”；循环相关与周期信号的归一化互为补充。
- M16_测量描述分析评估分层：参照带诊断法是”Describe + Data on normals”层级的标准实现；CV 追踪属”监测”，不得直接宣称因果。
- M15_过用损伤与应力：加载速率/总负荷时长是过用损伤核心变量，归一化红线直接保护这类指标的效度。
- M19_运动学测量系统对比（下一模块）：采集系统的噪声水平决定单条曲线质量，进而决定集合平均的可靠性。
- 支撑力矩（M16）：协方差分析（本节点到，详见 Winter 第 11 章，M25 前后）是其变异结构分析工具。
- 项目产出关联：Functional Profile Framework 的”个体典型曲线”= 集合平均，“动作稳定性”维度 = 波形 CV；正常参照带是 Metric Dictionary 各指标的标配参照形态。
12. Evidence 证据
- Winter 5th ed., Chapter 2, 2.3（书页 29–32）：集合平均的动机与”角度变异小、力矩变异大”现象；Figure 2.20 髋骨关节炎 vs 20 名年龄匹配对照的支撑期髋角叠加图（“stiff-leg gait”）；2.3.2 归一化插值算例（n=107 → N=100）及”加载速率/总负荷不得在归一化曲线上计算”的警告；2.3.3 波形 CV 公式（式 2.20）与 rms 分母讨论。
- 关联引用：Gage et al. (2004) 倒摆模型；Carpenter et al. (2001) COP 记录时长；Nelson-Wong et al. (2008) 站立任务臀中肌模式（第 11 章展开）。
- [补充说明] 现代步态报告软件（如 Vicon Nexus、Visual3D）的”normal gait band”图即 Figure 2.20 范式的工业实现；功能数据分析（FDA/fPCA）是集合平均的现代推广，超出本节范围。
13. My Understanding 我的理解
- 本节把”一个人”变成统计学对象：集合平均给出个体典型、CV 给出个体稳定性、参照带给出人群坐标——三者合起来就是功能表现型的最小统计单元。本系统 Functional Profile 的每个指标都应存”均值曲线 + SD 带 + CV”三件套。
- 归一化红线要写进数据管线：加载速率类指标在归一化曲线上计算是方法学错误而非误差——Squat MVP 的管线设计应为”原始时间轴算速率类指标，归一化轴只用于形态比较与分型”，两条轨道分开。
- “角度稳、力矩不稳”是机制推理的重要先验：运动学稳定不代表动力学稳定，功能分型若只用角度数据会系统性高估稳定性——MVP 若只有 IMU/视觉数据，报告中必须声明这一盲区。
- 参照带是描述层的上限：出带定位异常模式已是纯运动学能做的极限；再往前一步（为什么）就必须进入力矩、EMG、能量层——这正是 M21 之后动力学深化的存在理由。
- CV 的双重身份要小心使用：它既是”技术熟练度”指标（运动员低 CV 好），也可能是”病理代偿僵化”或”测量噪声大”的表象——CV 的解读必须绑定人群与任务背景，Metric Dictionary 中应为其标注双向解释。
模块自测（附参考答案）
Q1. 什么是集合平均？它对个体内与个体间数据分别如何实施？主要收益是什么？
参考答案：集合平均是将重复性运动的多次循环（个体内：同一受试者多条周期曲线）或多个体（个体间：每人先平均再对群体平均）按周期对齐后逐点平均。收益：平均波形更可靠（随机扰动按 √m 缩小），且平均周围的变异给出该变量的随机性画像，是临床参照带与协方差分析的基础。
Q2.（计算）某受试者一个步态周期采得 51 点膝角数据，要归一化到 100 点。求归一化曲线上第 10 个点 y₁₀ 的插值表达式。
参考答案：间隔 = n/N = 51/100 = 0.51；第 10 点对应原始位置 10×0.51 = 5.1，即第 5 与第 6 采样点之间，小数 0.1 为右邻权重：y₁₀ = (1−0.1)·x₅ + 0.1·x₆ = 0.9x₅ + 0.1x₆。
Q3.（计算）某患者步态周期 100 个归一化点上，平均标准差为 2.5°，平均信号绝对值为 25°。求波形变异系数 CV 并解释。
参考答案：CV = 2.5 / 25 = 0.10 = 10%。含义：该变量在平均波形周围的平均离散度约为其典型幅度的 10%，可与本人既往或人群的 CV 横向比较动作稳定性。
Q4. 为什么不能在时间归一化到 100% 的曲线上计算加载速率或总负荷？
参考答案：归一化把每个人的真实时间（秒）非均匀地拉伸/压缩为百分比，时间量纲被抹除且各人缩放率不同；而加载速率（力/秒）与总负荷时长（秒）以真实时间为因子，在归一化曲线上计算会得到系统性失真的结果。下肢承重组织完整性与加载速率和时长直接相关，这类指标必须回到原始时间轴计算。
Q5. 步态研究中发现”个体内关节角度变异很小，而同一关节力矩曲线变异很大”。这一现象说明什么？它与神经肌肉系统的哪一性质相关？
参考答案：说明 CNS 能以多种不同的力矩组合实现几乎相同的运动学输出——运动学层是最终”被约束”的目标，力矩层保留了大量冗余自由度。这与 M16 的四级汇聚/变异性层级（EMG > 力矩 > 运动学）及 Bernstein 冗余原理一致；正是这一现象催生了关节间协方差分析与支撑力矩研究。
Q6. 临床报告中将患者髋角曲线叠加于”20 名年龄匹配健康人平均 ± 1 SD”的参照带之上，发现支撑初期屈曲不足、末期伸展不足。按规范读图顺序写出结论，并指出该方法的边界。
参考答案：读图：形态上支撑初期髋屈曲减少、支撑末期髋伸展缺失，整体活动幅度下降，符合”僵直腿步态（stiff-leg gait）“模式。边界：① 这是描述层结论，只能说明”哪里、多少”，不能说明原因（疼痛回避？挛缩？肌力不足？需动力学/EMG/查体跟进）；② 参照带效度依赖人群匹配（年龄、性别、步速）；③ 偶然单次异常可能已被平均抹掉，需回看原始曲线。
上一模块：M17_相关分析与频谱｜下一模块：M19_运动学测量系统对比