模块定位：这是整个测量阶段”最技术性”也最容易翻车的一节。动捕给出的原始坐标含高频噪声，而我们要的速度、加速度是对坐标的一次、二次微分——微分会把噪声按频率倍数放大（速度 ×n、加速度 ×n²）。本模块回答：噪声从哪来、为什么必须先平滑、四种平滑方法孰优孰劣、截止频率怎么选（残差分析）、以及平滑之后如何正确计算角度、速度、加速度。这些规则直接成为 Squat MVP 数据处理管线的技术规范。
0. 大白话定义
动捕相机拍下的坐标曲线，肉眼看很平滑，其实上面叠着一层细砂般的噪声。位移曲线上这层砂几乎看不见，但如果你拿这条曲线去算速度，砂粒就被放大 20 倍；算加速度，再放大 20 倍——位移上 0.5 毫米的噪声，到加速度上能比真实信号还大 4 倍。所以原始坐标必须先”打磨”（低通滤波）再微分。
打磨也有讲究：单向滤波会让曲线整体”迟到”（相位滞后），所以正向滤一遍、反向再滤一遍，时间失真归零。打磨力度的选择不能拍脑袋：用残差分析——不断收紧滤波，看”滤掉的成分”何时从纯噪声变成伤及信号，拐点就是最佳截止频率。
1. Definition 定义（术语表）
| 术语 | 英文 | 定义 |
|---|---|---|
| 原始数据 | Raw data | 动捕转换出的未处理坐标序列，含叠加性随机噪声（电子噪声/空间精度/数字化误差） |
| 微分噪声放大 | Differentiation noise amplification | 第 n 次谐波经一次微分幅值 ×n，二次微分 ×n²——高频噪声在速度/加速度中占比暴增 |
| 曲线拟合平滑 | Curve fitting | 假定轨迹服从多项式/谐波级数/样条，用最优拟合替代含噪数据 |
| 样条拟合 | Spline fitting | 分段多项式拟合（以拐点分段），段间特殊衔接 |
| 相位滞后 | Phase lag | 低通滤波使输出相对输入的时间延迟；二阶滤波器在 fc 处滞后 90° |
| 双向滤波/零滞后滤波 | Dual-pass / zero-lag filtering | 正向滤波后沿时间反向再滤一次，相位超前抵消滞后，净相移为零，等效四阶滤波器（Winter et al., 1974） |
| 截止频率修正因子 | Correction factor C | 多遍滤波使有效截止频率下移，设计单遍 fc 时须除以 C 补偿：Butterworth C=(2^(1/n)−1)^0.25；临界阻尼 C=(2^(1/2n)−1)^0.5 |
| 残差分析 | Residual analysis | 计算各 fc 下滤波与原始信号的均方根残差，作残差-fc 曲线，选”信号失真=噪声通过”的 fc（Wells & Winter, 1980） |
| 最优截止频率 | Optimal cutoff frequency | Yu et al. (1999) 公式 fc,2 = 0.06fs − 0.000022fs² + 5.95ε，同时考虑采样率与残差，加速度用值高于位移 |
| 有限差分 | Finite difference | 用相邻采样点差商近似导数：中心差分 (x_{i+1}−x_{i−1})/2Δt |
| 中心差分 | Central difference | 用 i−1 与 i+1 点求第 i 点斜率，使导数与位移在时间上对齐 |
| 环节角 | Segment angle | θij = arctan[(yj−yi)/(xj−xi)]，环节在绝对参照系中的角度 |
| 关节角 | Joint angle | 相邻两环节角之差（按各关节约定），如膝 θk = θ21 − θ43 |
| 半功率点 | Half-power point | 滤波器响应降至 0.707 的频率，即截止频率定义 |
2. 核心公式与推导
2.1 信号与噪声的频谱分布（3.5.2）
- 步行 120 步/分 → 步频 2 Hz、步态周期频率 1 Hz，谐波位于 1, 2, 3… Hz；跑动频率含量整体上移；远端环节（足跟/足尖）频率高于近端（膝/躯干）。
- 20 名受试者足尖轨迹谐波分析（Winter et al., 1974）：99.7% 信号功率 < 7 次谐波（< 6 Hz）；7–20 次谐波仍有小幅值成分，具噪声特征（与 M17 图 2.17 同源）。
- 噪声 = 非运动过程本身产生的成分：光电设备电子噪声、扫描/数字化空间精度、标记重建误差；总和为高频随机叠加。
2.2 微分噪声放大的数学（3.5.3，式 3.3–3.5）
设位移信号为谐波和：x = Σ Xn·sin(nω₀t + θn)，则
速度 Vx = Σ nω₀·Xn·cos(nω₀t + θn) → 第 n 次谐波幅值 × n
加速度 Ax = −Σ (nω₀)²·Xn·sin(nω₀t + θn) → 第 n 次谐波幅值 × n²
书中算例（必须记住的数量级）：步态轨迹中基波幅值 x₁ = 5 cm，第 20 次谐波噪声 x₂₀ = 0.5 mm（位移图上几乎不可见）——
- 速度中：20 次谐波放大 20 倍 → 达基波的 1/5；
- 加速度中：再放大 20 倍 → 达基波的 4 倍。
图 3.18 的实证：未滤波坐标二次差分得到的足尖加速度曲线呈随机噪声状，滤波后同一计算给出光滑、可重复的加速度。
2.3 四种平滑方法对比（3.5.4–3.5.5）
| 方法 | 原理 | 致命弱点 |
|---|---|---|
| 多项式拟合（式 3.6） | 假定轨迹为 n 阶多项式，最小均方误差拟合 | 高阶才勉强贴合，Pezzack 实验中 9 阶多项式加速度与实测严重不符；计算昂贵 |
| 谐波重构（式 3.7） | 取最低 N 次谐波重建 | 各谐波幅值/相位的平稳性假设不成立（M17 已证：足部摆动期高频、支撑期低频） |
| 样条拟合 | 以拐点分段的多项式 | 拐点须从含噪数据本身确定——被要消除的噪声污染 |
| 数字滤波（推荐） | 利用信号/噪声频谱差异衰减噪声 | 引入相位滞后——由双向滤波消除 |
Pezzack et al. (1977) 经典验证实验（同一摆臂运动，四种方法 vs 加速度计实测）：
- (a) 测角计 vs 动捕角度：吻合（测角计稍噪）；
- (b) 数字滤波 + 有限差分 vs 加速度计：极好吻合，且比直接实测更干净；
- (c) 9 阶多项式拟合加速度 vs 实测：严重偏差，足以否定多项式法；
- (d) 原始坐标直接差分 vs 实测：噪声大到”毫无意义（too noisy to mean anything）“。
2.4 双向零滞后滤波（3.5.4.2）
- 二阶低通在 fc 处相位滞后 90°，通带内高次谐波产生相位失真——可能比过渡带的幅值失真更严重。
- 解决：正向滤波后，沿时间反向再滤一次（Winter et al., 1974）→ 引入等量相位超前，净相移 = 0，等效四阶零滞后滤波器，且截止特性加倍陡峭。
- 代价：有效截止频率下移至约 0.80fc。设计时用修正因子补偿（式 3.8：ωc = tan(πfc/fs)/C）：
- Butterworth：C = (2^(1/n) − 1)^0.25，双遍 n=2 → C = 0.802
- 临界阻尼：C = (2^(1/2n) − 1)^0.5，双遍 → C = 0.435
- 二者时域特性权衡：Butterworth 阶跃响应有轻微过冲但上升快；临界阻尼无过冲但上升慢。人体运动数据罕有冲激型输入 → 优选 Butterworth。
- 系数设计（同 M17 式 2.19）：a₀ = K₂/(1+K₁+K₂)，a₁ = 2a₀，a₂ = a₀，b₁ = −(2a₀+K₃)，b₂ = 1−2a₀−K₃；系数只取决于 fs/fc 比值与 C。
2.5 截止频率的选择——残差分析（3.5.4.3，Wells & Winter, 1980）
对每个候选 fc 计算残差（式 3.9）：
R(fc) = √[ (1/N)·Σᵢ (Xᵢ − X̂ᵢ)² ]
Xᵢ 为原始数据，X̂ᵢ 为该 fc 下四阶零滞后滤波结果。作 R–fc 曲线（Figure 3.19）：
- 若数据为纯噪声，残差是一条从 0 Hz 处截距（= 噪声 rms 值）到奈奎斯特频率处零点的直线；
- 含真实信号时，fc 越低残差越向上偏离直线——偏离量 = 信号失真；
- 决策：从截距 a 水平投射交残差曲线于 b，对应频率 fc′，此时信号失真（bc 段）= 通过的噪声，二者平衡。
步态实测（Figure 3.20）：4 个标记点残差回归线 Residual = 1.8 − 0.05n（r = 0.79）→ 噪声 rms = 1.8 mm（相机 5 m 外、视场 2×3 m，噪声 < 视场 1/1000）；足跟/跖骨点信号功率达约 6 Hz，肋/髋点仅约 3 Hz → 不同标记点可用不同截止频率。
2.6 高阶导数的最优截止（3.5.4.4）
- Giakas & Baltzopoulos (1997)：最优 fc 依赖噪声水平与求导阶数（但其基准信号为谐波重构，受 M17 平稳性批评）。
- Yu et al. (1999) 公式（式 3.10）：
fc,2 = 0.06·fs − 0.000022·fs² + 5.95·ε
fs 为采样率，ε 为滤波/原始数据的相对平均残差。要点：算加速度的最优 fc 高于位移残差分析给出值；且 fs 越高 Δt 越小，有限差分算出的噪声越大（见式 3.17/3.18c 的分母），fc 需相应调整。
2.7 角度、速度、加速度的正确计算（3.6）
环节角（式 3.13/3.14）：θij = arctan[(yj − yi)/(xj − xi)]——两标记点不必在环节端点，只要与长骨轴共线；角度逆时针为正、水平为 0°。
关节角（3.6.2）：
- 膝：θk = θ21 − θ43（θ21 大腿角、θ43 小腿角）；θ21 > θ43 为屈曲；
- 踝：θa = θ43 − θ65 + 90°；正 = 跖屈，负 = 背屈；
- 关节角约定在研究者间差异极大，报告必须声明所用约定（绝对环节角约定统一后角速度/角加速度符号才能与后续动力学一致）。
速度——中心差分（式 3.15/3.16）：
Vxᵢ = (xᵢ₊₁ − xᵢ₋₁) / 2Δt ωᵢ = (θᵢ₊₁ − θᵢ₋₁) / 2Δt
为什么不用 (xᵢ₊₁−xᵢ)/Δt：那种前向差分的速度对应两点中点时刻，与位移数据错半拍，后续联合分析会错位。中心差分假设 i−1 到 i+1 连线斜率 = i 点切线斜率，时间对齐。
加速度（式 3.17 / 3.18a–c）：
- 对速度再中心差分（式 3.17）需 i±2 共 5 个点；
- 更优的三点式（式 3.18c，由两个”半点速度”相减推得）：
Axᵢ = (xᵢ₊₁ − 2xᵢ + xᵢ₋₁) / Δt²
角加速度以角数据代换即可。前提：所有微分都作用在已平滑数据上。
3. Physical Meaning 物理意义
- 微分是高通操作：乘以 nω₀ 等价于按频率加权——微分天然偏爱高频，而噪声恰在高频，这是矛盾的数学本质。
- 滤波的相位是时间轴的形变：幅值失真改变”多少”，相位失真改变”何时”——对时序敏感的变量（峰值时刻、事件同步）相位失真更致命，故零滞后不可省。
- 残差直线 = 噪声的指纹：白噪声在所有频率等功率，滤波移除的噪声量随带宽线性变化，故残差-fc 为直线；信号是低频集中的，收紧滤波时信号损失”加速”——直线与曲线的分岔就是信噪边界。
- Δt 的双刃剑：采样率翻倍使 Δt 减半，有限差分把固定大小的坐标噪声放大一倍——高采样率必须配相应的滤波策略，不是免费的午餐。
4. Human Movement Meaning 人体运动意义
- 步行数据 6 Hz 截止的经验法则来自真实谐波结构：信号功率 99.7% 在 7 次谐波以下；但远端快速环节（足尖摆动期）与近端慢速环节（髋/肋）频率上限不同（6 Hz vs 3 Hz）——统一 fc 是妥协，分标记点 fc 更优。
- 跑动、冲刺、投掷频率整体上移：步行的 6 Hz 法则不可外推，每个项目/动作必须重做谐波或残差分析。
- 病理步态：更慢、更不规则的步态基频更低、平稳性更差——谐波重构类方法风险更大，数字滤波 + 残差分析的稳健性优势更突出。
- 平滑后曲线”穿过噪声曲线中部”（图 3.18）：这是判断滤波合理性的直观标准——残差应像围绕信号的随机散布，而非系统性偏离。
5. How It Is Measured 标准处理管线
- 检查原始数据：绘制位移曲线，识别缺失点（先做 gap filling）与明显伪迹。
- 估计噪声水平：残差分析回归线截距即噪声 rms（或从静态标定段估计）。
- 选 fc：谐波分析（初步）→ 残差分析（精细，信号失真=噪声通过点）→ 若目标变量是加速度，用 Yu 公式上调。
- 设计滤波器：二阶 Butterworth，按 fs/fc 与双遍修正 C=0.802 算系数（式 3.8 及 M17 式 2.19）。
- 双向滤波：正向一遍 + 反向一遍 → 四阶零滞后。
- 微分：中心差分算速度（式 3.15/3.16），三点式算加速度（式 3.18c）——全部在平滑数据上。
- 验证：滤波曲线应穿过原始曲线中部；有条件时用加速度计等独立传感器交叉验证（Pezzack 范式）。
- 记录元数据：fc、滤波器类型/阶数/遍数、差分格式——写入指标字典元数据。
6. Typical Applications 典型应用
| 场景 | 处理要点 |
|---|---|
| 步态实验室标准流程 | 6 Hz（步行）四阶零滞后 Butterworth + 中心差分 |
| 跑动/冲刺 | 重做残差分析，fc 上移（常 8–15 Hz [补充说明]） |
| 加速度敏感分析（落地冲击、功率） | Yu 公式定 fc；优先三点式 3.18c |
| 多标记点不同频率特性 | 足端 6 Hz、躯干 3 Hz 分点处理 |
| 与 IMU/加速度计融合验证 | Pezzack 范式：独立传感器验证滤波+差分链 |
7. Clinical Relevance 临床相关性
- 临床报告的每个派生变量都依赖本节管线：步态报告的关节角速度/加速度若省略平滑说明，数值不可复现、不可跨实验室比较。
- 加速度类指标对处理最敏感：落地冲击、负荷速率（衔接 M15_过用损伤与应力 的过用损伤负荷指标）在 fc 选择不当时可差数倍——临床结论可能被处理方法而非病情驱动。
- 病理数据更脆弱：不规则步态使平稳性假设进一步破坏，谐波法与多项式法风险更高；残差分析应作为每个新人群的必经步骤。
- “滤波曲线穿过噪声中部”可培训为临床质控直觉：快速目检即可发现滤波过度（曲线跟不上真实快速事件）或不足（加速度曲线毛刺）。
8. Performance Relevance 表现相关性
- 快速动作是微分噪声的重灾区：投掷/短跑中远端环节速度高、事件窗口短，过度平滑会削掉峰值加速度——而峰值恰是表现与损伤评估的核心量。
- 峰值保真 vs 噪声压制的权衡在表现分析中最尖锐：宁可保留少量噪声换取峰值保真，并在报告中声明噪声水平，也不要把真实爆发力”滤没”。
- 技术对比（训练前后/运动员间）必须使用同一处理参数：fc 或差分格式不同，峰值差异可能纯属处理伪差。
9. Possible Interpretations 可能的解释层级
| 观察 | Level B（描述层） | Level C（机制层解释） | Level D（常见误读） |
|---|---|---|---|
| 加速度曲线高频毛刺 | 存在高频成分 | 大概率为差分放大的噪声（对照频谱/残差） | “真实的高频冲击” |
| 滤波后峰值下降 15% | 峰值被平滑 | fc 偏低致信号失真，或原峰值本含噪声——需残差分析仲裁 | ”滤波前后差异=误差”或”滤波后才是真相” |
| 峰值时刻前移 | 事件提前 | 单向滤波相位滞后未校正也会改变峰值时刻 | 在相位失真数据上解读时序 |
| 两实验室报告同一运动员峰值加速度不同 | 数值不一致 | 优先怀疑处理参数（fc/差分格式/采样率）差异 | 直接归因于状态变化 |
10. What It Cannot Tell Us 这套方法不能告诉我们什么
- 滤波不能恢复被采丢的信息：欠采样混叠（M17）在先，任何平滑都无力回天。
- 残差分析给的是平衡点不是真值：“信号失真=噪声通过”是决策准则，不证明该 fc 最优到每个变量。
- 平滑数据上的导数仍是估计：中心差分本身有截断误差，平滑+差分链条的累积不确定度需以独立传感器标定。
- 不能跨动作套参数：6 Hz 是步行的答案，不是人体运动的答案。
- 处理链不能修复系统性误差：软组织伪迹、标记错贴不在本节方法的能力范围内（见 M19_运动学测量系统对比）。
11. Related Metrics 相关概念与双链
- M19_运动学测量系统对比：噪声的物理来源（各系统精度）决定本节残差水平。
- M17_相关分析与频谱：采样定理、Butterworth 系数设计、白噪声谱——本节全部建立其上。
- M18_集合平均与归一化：平滑与集合平均是两条降噪路线（频域 vs 重复平均），管线中先平滑后平均。
- M15_过用损伤与应力：负荷速率类指标的效度直接受 fc 选择影响。
- M20_人体测量学与环节参数（下一模块）：平滑坐标 + 环节参数 → 质心轨迹，是逆动力学输入链的下一环。
- Squat MVP 管线规范：四阶零滞后 Butterworth（双遍 C=0.802）+ 中心差分速度 + 三点式加速度 + 残差分析定 fc + 元数据全记录。
12. Evidence 证据
- Winter 5th ed., Chapter 3, 3.5–3.6（书页 52–64）：式 3.3–3.5 微分放大推导与”20 次谐波加速度达基波 4 倍”算例；图 3.15 足尖谐波（99.7% < 7 次）；图 3.16 微分幅值-频率关系；四种平滑方法对比与 Pezzack et al. (1977) 验证实验（图 3.21，9 阶多项式式 3.11–3.12 失败案例）；双向滤波原理与 C 因子（Butterworth 0.802 / 临界阻尼 0.435，式 3.8）；图 3.17 二阶 vs 四阶零滞后响应；图 3.18 足尖加速度滤波前后；残差分析式 3.9、图 3.19 决策法、图 3.20 步态实测（噪声 rms 1.8 mm，足跟/跖 6 Hz vs 肋/髋 3 Hz）；Yu et al. (1999) 最优 fc 式 3.10；环节角式 3.13–3.14、膝/踝角约定、中心差分式 3.15–3.17、三点加速度式 3.18a–c；章末 Problems（足尖离地间隙答案 1.52 cm 等）。
- [补充说明] 现代实现：SciPy
signal.filtfilt即双向零滞后滤波；savgol_filter（Savitzky-Golay）为常用替代，但本章的残差分析决策逻辑仍然适用。
13. My Understanding 我的理解
- “先滤波后微分”不是建议而是数学必然：噪声 ×n² 的放大律意味着任何未经平滑论证的加速度数字都不可引用——这条应成为本系统数据管线的强制门禁。
- 残差分析是”可辩护的 fc”：拍脑袋定 6 Hz 是经验，残差曲线给出每个数据集自己的信噪分岔点——把经验变成了可审计的决策。Metric Dictionary 应要求每个派生指标附 fc 及其确定依据。
- 相位滞后比幅值失真更阴险：幅值失真肉眼可见，相位失真悄悄移动峰值时刻——时序类机制推理（如 M17 的 70 ms 提前量）若建立在单向滤波数据上，结论可能整个错向。零滞后是时序研究的生命线。
- Pezzack 实验是方法学自信的来源：滤波+差分链可以比直接加速度计更干净——信号处理不是”不得已的妥协”，而是能超越硬件的估计艺术；但前提是像他们一样用独立传感器验证。
- 对 MVP 的直接约束：Squat 动作慢（基频约 0.3–0.5 Hz [补充说明]），fc 可低至 2–4 Hz，但必须用自己的残差分析确认；且下蹲离心/向心转换点的角加速度峰值对 fc 敏感，报告需标注。
模块自测（附参考答案）
Q1. 为什么位移曲线上几乎不可见的高频噪声，会在二次微分后淹没真实加速度信号？用书中 x₁ = 5 cm、x₂₀ = 0.5 mm 的算例说明。
参考答案：第 n 次谐波一次微分幅值 ×n、二次微分 ×n²。位移上 0.5 mm 的 20 次谐波噪声仅为基波（50 mm）的 1/100，不可见；速度中放大 20 倍 → 10 mm，为基波速度贡献的 1/5；加速度中再放大 20 倍 → 200 mm 当量，为基波加速度贡献的 4 倍——噪声彻底淹没信号（图 3.18 的随机状曲线）。
Q2. 四种平滑方法中为什么数字滤波成为标准？其余三种的根本缺陷各是什么？
参考答案：数字滤波利用信号（低频）与噪声（高频）的频谱差异衰减噪声，不预设信号形状，稳健性最好。其余：① 多项式——Pezzack 实验中 9 阶仍严重失真，高阶计算昂贵；② 谐波重构——要求各谐波幅值/相位全程平稳，而足部等信号摆动期高频、支撑期低频，假设不成立；③ 样条——分段拐点需从含噪数据本身确定，被待消除的噪声污染。
Q3.（计算）用二阶 Butterworth 双向滤波处理 100 Hz 采样数据，目标有效截止频率 8 Hz。设计单遍滤波器时 fc 应设为多少？（C = 0.802）
参考答案：双遍后有效截止 ≈ 0.802 × 单遍设计值，故单遍 fc = 8 / 0.802 ≈ 9.98 ≈ 10 Hz。即在式 3.8 中以 fc = 10 Hz 计算 ωc 与滤波器系数，双向滤波后得到有效 8 Hz 截止。
Q4.（计算）平滑后坐标数据 x₉ = 1.204 m、x₁₀ = 1.218 m、x₁₁ = 1.233 m，采样 50 Hz。用中心差分求第 10 点速度，用三点式求第 10 点加速度。
参考答案：Δt = 1/50 = 0.02 s。速度 V₁₀ = (x₁₁ − x₉)/2Δt = (1.233 − 1.204)/0.04 = 0.725 m/s。加速度 A₁₀ = (x₁₁ − 2x₁₀ + x₉)/Δt² = (1.233 − 2×1.218 + 1.204)/0.0004 = 0.001/0.0004 = 2.5 m/s²。
Q5. 为什么计算速度时推荐中心差分 (xᵢ₊₁−xᵢ₋₁)/2Δt 而非前向差分 (xᵢ₊₁−xᵢ)/Δt？
参考答案：前向差分结果对应两采样点的中点时刻，与位移数据错开半个采样间隔，后续把速度与位移/力数据联合分析时产生时序错位；中心差分用 i−1 与 i+1 对称点，其连线斜率近似 i 点切线斜率，导数与原始数据时间上对齐。
Q6. 描述残差分析选择截止频率的完整逻辑。若残差回归线截距为 1.8 mm，这一数字的含义是什么？
参考答案：对一系列候选 fc 分别计算滤波与原始信号的 rms 残差并作 R–fc 曲线。纯噪声的残差为从 0 Hz 截距到奈奎斯特频率零点的直线；实际曲线随 fc 降低而上翘，上翘量即信号失真。从截距水平投射与曲线相交，交点处信号失真 = 通过的噪声，对应 fc 即平衡点。截距 1.8 mm = 0 Hz”滤波”（即取均值）时的残差 = 坐标噪声的 rms 值，本例中相对 2×3 m 视场不足 1/1000，且各标记点回归线相同说明噪声主要来自数字化过程。
Q7. 为什么 Yu et al. (1999) 给出的加速度最优截止频率高于位移残差分析的结果？采样率提高对差分噪声有何影响？
参考答案：加速度把噪声按 n² 放大，噪声在加速度中占比远高于位移，为保住信号的高频真实成分同时控制噪声，最优平衡点向更高 fc 移动。采样率提高使 Δt 减小，而差分公式分母为 Δt（速度）或 Δt²（加速度），同样大小的坐标噪声被除得更”狠”——差分噪声随之增大，故最优 fc 也是 fs 的函数（式 3.10 含 fs 项）。
Q8. 一名研究员比较患者术前/术后的膝关节峰值角加速度，术前数据用 fc = 6 Hz 处理、术后用 fc = 4 Hz，报告”峰值下降 30%，手术效果显著”。指出两处方法学错误并给出正确做法。
参考答案：① 处理参数不一致：fc 不同会直接造成峰值差异（4 Hz 平滑更强、峰值更低），“下降 30%“可能全是处理伪差——两组数据必须用完全相同的 fc、滤波器、差分格式重处理；② 把组间差值直接归因于手术属于监测≠因果（M16）：还需对照组或自然变异数据排除其他解释。正确做法：统一管线重处理 → 报告 fc 及残差分析依据 → 结合力矩/功能指标多层面验证 → 谨慎表述因果强度。
上一模块：M19_运动学测量系统对比｜下一模块：M20_人体测量学与环节参数