M23 · 三维变量与事件检测（3D Variables and Event Detection） 阶段五 · 动力学与能量学 前面章节的逆向动力学是二维的；本章回答：相机拍到的只是一堆标记点的 x, y, z 坐标，如何一步步变成解剖坐标系下的关节角度、力矩和功率？ 核心链条是：全局坐标 → 标记坐标 → 解剖坐标 → Cardan 角 → 角速度/角加速度 → 3D 牛顿-欧拉方程 → 三维力矩与功率。事件检测（HC/TO）则是把所有数据”对齐到同一条时间轴”的前提。
0. 大白话定义
坐标变换像”翻译接力”。相机只懂实验室的语言（全局坐标系 GRS），肌肉骨骼只懂身体的语言（解剖坐标系）。中间请了一位随行翻译——绑在腿上的三个标记点组成的”标记坐标系”。先做一次校准，记住”标记系”和”解剖系”之间的固定关系（这本词典一次编好，全程不变）；之后每一帧，只要把相机语言翻成标记语言（每帧都要重翻），再查固定词典翻成解剖语言，就完成了。
Cardan 角像报一个姿态的”三步口令”：先绕 X 轴转 θ₁，再绕新的 Y 轴转 θ₂，最后绕新的 Z 轴转 θ₃。顺序不能乱——就像”先左转再前进”和”先前进再左转”结果完全不同（矩阵乘法不可交换）。
事件检测就是给运动录像”打拍子”：脚跟着地（HC）是”哒”，脚尖离地（TO）是”哒”。所有受试者的曲线都按这个拍子对齐成 0–100%，才能互相比较。健康人用测力台打拍子最准（力超过 20 N 就算着地）；但脑瘫患者可能根本没有”脚跟先着地”，就得换运动学方法打拍子。
1. Definition 术语表（中英对照）
| 英文术语 | 中文 | 一句话定义 |
|---|---|---|
| Global Reference System (GRS) | 全局参考系 | 固定在实验室空间的坐标系，与测力台轴向对齐 |
| Local Reference System (LRS) | 局部参考系 | 附于体节的坐标系，分标记坐标系与解剖坐标系两种 |
| Anatomical Axis System | 解剖坐标系 | 由骨性标志定义的体节主轴坐标系，原点通常在环节质心 |
| Cardan Sequence | Cardan 转序 | 按 x–y–z 顺序的三次连续旋转，生物力学常用 |
| Euler Angles / Sequences | 欧拉角 / 转序 | 12 种合法旋转序列的统称；z–x–z 为工程学常用 |
| Gimbal Lock | 万向锁 | 两个转轴平行导致一个自由度退化的数学奇点 |
| Tracking / Calibration Markers | 跟踪 / 校准标记点 | 全程贴附的标记 vs 仅静态校准时临时贴附的标记 |
| [G to M] / [M to A] / [G to A] | 全局→标记 / 标记→解剖 / 全局→解剖矩阵 | 三个 3×3 旋转矩阵；[G to A] = [M to A] × [G to M] |
| Heel Contact (HC) / Toe-Off (TO) | 脚跟着地 / 脚尖离地 | 划分支撑期与摆动期的步态事件 |
| Time Normalization | 时间归一化 | 将数据重采样至步态周期（或支撑/摆动期）的 0–100% |
| Induced Acceleration Analysis | 诱导加速度分析 | 研究一个关节力矩对全身所有相连环节加速度影响的方法 |
2. 核心框架与公式
2.1 三套坐标系与两级变换
- GRS：X 前后、Y 左右、Z 竖直，与测力台对齐；用 L 形标定架放在测力台上完成空间标定，相机一旦移动必须重新标定，即使不动也要定期标定防漂移。
- 标记坐标系 xm–ym–zm：每个环节至少 3 个不共线标记点构成跟踪平面；常取 mT1 为原点，mT1→mT3 为 zm，ym 垂直于跟踪平面，xm 按右手系补全。
- 解剖坐标系：原点设在环节质心，z 轴沿环节长轴（骨盆则由 PSIS/ASIS 等骨性标志定义）。
- 两级变换：[G to A] = [M to A] × [G to M]。其中 [M to A] 由静态校准得到、假设恒定；[G to M] 逐帧重算。解出 [G to A] 后按选定转序反解三个时变角 θ₁、θ₂、θ₃。
- 平移部分：环节质心全局位置 Rc = Rm + c（Rm = 标记原点全局坐标，c = 标记原点到质心的常矢量，需先转到全局系）。
2.2 Cardan x–y–z 旋转（式 7.1–7.5）
三次旋转各对应一个旋转矩阵 Φ₁（绕 x）、Φ₂（绕新 y）、Φ₃（绕新 z），合成：
[G to A] = Φ₃Φ₂Φ₁（顺序不可交换！）
展开式（7.5）的底行三个元素最简洁：s₂, −s₁c₂, c₁c₂——反解角度时从此行入手：
- s₂ = 矩阵 (3,1) 元素 → θ₂
- −s₁c₂ = (3,2)、c₁c₂ = (3,3) → θ₁
- c₂c₃ = (1,1)、−c₂s₃ = (2,1) → θ₃
- 每个角都有两个候选解，必须用其余元素回代验证唯一解。
算例（7.2.5.3，摆动期 frame 6）：解得 θ₁ = −48.7°、θ₂ = −17.21°、θ₃ = −69.62°。最大的 θ₃ 反映摆动期小腿相对实验室的大角度旋转。
理论上共 12 种合法转序（Euler 提出）；生物力学常用 Cardan（x–y–z），机械工程常用 Euler（z–x–z）。万向锁（gimbal lock）：当某次旋转使两个转轴平行时，3D 姿态无法唯一确定——上肢大活动度分析时尤其要按 ISB 建议选转序规避。
2.3 点积与叉积（式 7.6）
- 点积（标量积）：3D 功率 P = F·V = FxVx + FyVy + FzVz（第 6 章的推广）。
- 叉积（矢量积）：C = A × B，方向由右手定则、垂直于 AB 平面，|C| = |A||B|sinθ。建坐标系的基本工具：x = y × z，再用 x × z 修正 y，保证三轴严格正交（见 2.4 校准流程）。
- 所有矢量最后都要归一化为单位矢量再拼成旋转矩阵。
2.4 解剖校准流程（7.2.5.1 算例要点）
以小腿为例（Table 7.1 静态校准数据）：
- 关节中心 = 内外侧标记中点：踝 = (mT2 + mC1)/2，膝 = (mT3 + mC2)/2；
- 环节质心 = 0.567 × 膝 + 0.433 × 踝（体节参数）；
- z_an = 膝 − 踝；临时 y_an = mC1 − mT2（内外踝连线，仅近似垂直）；
- x_an = y_an × z_an；再修正 y_an = x_an × z_an，确保严格正交；
- 三个轴各自归一化 → 拼出 [LA to G]，转置得 [LG to A]（旋转矩阵的逆 = 转置）；
- 同样流程对三个跟踪标记建标记坐标系 → 得恒定的 [M to A]；校准标记（mC1、mC2）此后即可拆除。
临床适配：脑瘫、卒中患者无法摆标准解剖位 → 临床步态实验室改为舒适站立位校准 + anthropometric 测量（踝/膝宽度 + 通用 X 光体节数据推算关节中心偏移），标准流程见 Davis et al. (1991)、Õunpuu et al. (1996)。
2.5 标记点设计与标记组（Markerset）
- 尺寸由相机距离与分辨率决定：常规 10 mm；面部等密集场合 < 5 mm。
- 相机靠近红外反光的质心定位：两标记在视野中重叠会被误认为一个大反光点 → 两个都定位错误（图 7.3）。布点必须保证任意视角下可区分。
- 经典标记组：Helen Hayes（标记少、省时，但自由度受限、关节角误差较大）、Cleveland Clinic（下肢簇标记，利于横断面）、6DOF 标记组（大腿/小腿用簇板，融合前两者优点，Collins et al. 2009）。
- 实战坑：内侧标记互相碰撞脱落（解法：内侧标记仅校准时用，动态采集前拆除）；簇板布套沿肢体下滑 → 标记偏离校准参考位，产生显著误差；运动员出汗 + 高速碰撞需双面胶加固；反光涂层破损要及时更换。
2.6 步态事件检测（7.2.2）
- 测力台阈值法：垂直力升破/跌破 20 N 判定 HC/TO（20 N 是为排除噪声与摆动期脚擦台）。健康人好用；脑瘫等非典型着地模式常失效。
- 运动学法（Zeni et al., 2008）：
- 平地：HC = 足跟标记与骶骨标记前后向相对位移最大时刻；TO = 足尖标记与骶骨标记前后向相对位移最大时刻；
- 跑台： heel/toe 标记前后向速度由正转负 = HC，由负转正 = TO；
- 正常步速步行可靠，跑步分析可能不准。
- 采集纪律：整只脚必须完全落在单块测力台内、台面上不得有其他接触（拐杖、另一只脚、物品）；定期清零（tare）防电噪声/漂移制造假事件；运动学法对异常步态同样可能失效，需个案核验。
- 其他动作事件：跳起/落地——起跳与落地用垂直力阈值，起跳始动与落地完成用骨盆标记垂直位移偏离/恢复站立位；深蹲——以骨盆标记垂直位移的极大/极小值划分下蹲相与站起相（这正是我们 Squat MVP 可直接复用的事件定义）。
- 归一化：时间归一化到支撑/摆动期 0–100%；动力学量按研究惯例除以体重 / 身高 / 速度之一（各变量口径不同，参照领域标准）。
2.7 角速度与角加速度（式 7.7a–7.7b）
Cardan 角的时间导数不是环节角速度分量——必须经变换矩阵投影：
ω = [ c₂c₃, s₃, 0; −c₂s₃, c₃, 0; s₂, 0, 1 ] · [θ̇₁, θ̇₂, θ̇₃]ᵀ
角加速度 α 由 ω 数值微分（与 2D 相同的有限差分）。
2.8 3D 逆向动力学：牛顿-欧拉方程（式 7.8a–7.8c）
- 力在 GRS 中算（重力、测力台、质心加速度都天然在 GRS）：ΣF = ma，逐轴求近端反作用力；
- 力矩在解剖坐标系中算（转动惯量定义在解剖主轴上）：所有反作用力、远端力矩先用 [G to A] 变换；
- 欧拉方程（相对 2D 多出角速度耦合项 (I₁−I₂)ω₁ω₂）：
Ixαx + (Iz−Iy)ωyωz = ΣMx = Rzd·ld + Rzp·lp + Mxp − Mxd Iyαy + (Ix−Iz)ωxωz = ΣMy = Myp − Myd Izαz + (Iy−Ix)ωxωy = ΣMz = −Rxd·ld − Rxp·lp + Mzp − Mzd
- 算例（7.4.3，支撑期小腿 frame 6，60 Hz）：m = 3.22 kg，Ix = Iz = 0.0138、Iy = 0.0024 kg·m²。三步走：
- Step 1（GRS 反作用力）：RXP = 3.22 × 7.029 + (−119.04) = −96.41 N；RYP = 3.22 × 1.450 + (−12.03) + 31.60 = 24.24 N；RZP = 3.22 × (−0.348) + (−791.44) = −792.56 N；
- Step 2：全部力/力矩经 [G to A] 转到解剖系（Rzd = −781.34 N，Mxd = 58.29 N·m 等）；
- Step 3：角速度经式 7.7b 投影（ωz = −2.71 rad/s），角加速度由 frame 5/7 中心差分（αz = −17.08 rad/s²）；
- 解得膝关节三力矩：Mxp = +307.78 N·m（外展力矩，对抗上身体重作用于支撑腿内侧的重力负荷）[补充说明：此数值为教材原文结果，物理量级偏大，疑为教材排版错误，实际膝外展力矩通常为几十 N·m；复算时请以公式与原始数据自行验证]；Myp = −103.86 N·m（矢状面屈肌力矩，助力支撑末膝屈曲）；Mzp = −67.63 N·m（轴向力矩，反映髋内旋肌驱动骨盆前旋以增加步长）。
- 链条传递：算完环节的近端力矩后，需用 [A to G] 将其变回 GRS，作为上一环节（大腿）的远端输入，逐环节向上推进。
2.9 三维关节功率与诱导加速度
- 3D 功率 = 三个分量分别相乘再相加：P = Mxωx + Myωy + Mzωz（式 7.9），正 = 生成、负 = 吸收。
- 诱导加速度分析（Induced Acceleration Analysis）：任一关节力矩都会加速全身所有相连环节，无论多远。经典例子：肩屈时前臂因惯性相对肘伸；棒球投球的速度无法由单关节力矩解释，只能由环节间相互诱导的累积效应解释。[补充说明：这是 M25 肌骨建模中”肌肉功能 ≠ 其跨越关节的表观转动”的理论源头。]
2.10 正常步态 3D 力矩/功率曲线判读（Eng & Winter, 1995；图 7.5–7.6）
- 惯例：矢状面伸肌力矩为正、横断面外旋为正、额状面外翻为正；力矩按体重归一（N·m/kg），功率 W/kg。HC = 0%，TO ≈ 60%。
- 踝：步行中最大力矩在矢状面——HC 后小背屈力矩放平足，随后跖屈力矩增至约 50% 周期峰值 → A2-S 推离（push-off）正功爆发。
- 膝：8–25% 膝伸肌力矩控制承重屈膝（K1-S 吸收）；随后转屈肌力矩（腓肠肌贡献的副产品）；TO 前后小伸肌力矩限制屈膝（K3-S）；HC 前屈肌力矩减速摆动腿（K4-S）。
- 髋矢状面：前半支撑伸肌力矩（稳躯干 + 防膝塌陷 + “后方推送” H1-S 生成）；后半屈肌力矩（防躯干后仰 + 50–75% 周期”拉离（pull-off）“大腿入摆动 H3-S 生成）。
- 髋横断面：前半支撑外旋肌减速骨盆水平旋转（H1-T 吸收），后半内旋肌稳定前旋。
- 髋额状面：强外展力矩防止骨盆（连同上身，重力作用于支撑髋内侧约 10 cm）下垂——H1-F 吸收（骨盆下沉）→ H2-F/H3-F 生成（抬升骨盆与摆动腿，保障摆动期足趾廓清）。
- 膝额状面外展力矩不是肌肉产物：而是承重膝对重力负荷的被动响应——内侧髁加载、外侧髁卸载形成内源性外展力矩。韧带骨骼结构可以帮（或碍）肌肉的忙；假肢设计中的弹簧、阻尼、限位器同样会产生内力矩，康复工程必须计入。
3. Physical Meaning（物理意义）
- 旋转矩阵是正交矩阵：逆 = 转置，保长度、保角度——坐标变换本身不创造也不消灭任何物理量，只是换”语言”。
- 矩阵乘法不可交换 ↔ 三维旋转不可交换：这是”转序必须写进论文方法学”的物理根源，同一组姿态数据用不同转序会报告出不同的”关节角”。
- 欧拉方程中的 (I₁−I₂)ω₁ω₂ 项是陀螺耦合项：长杆状环节（Ix = Iz ≠ Iy）高速旋转时，一个平面内的角速度会在另一平面诱导力矩——2D 分析永远看不到这项。
- 万向锁是欧拉角参数化的固有奇点，不是测量误差。
4. Human Movement Meaning（人体运动意义）
- 二维分析只能看矢状面；但步行中最大的代偿与病理信号常在额状面与横断面：髋外展力矩（防骨盆下垂）、膝内/外翻负荷（内侧间室骨关节炎的核心力学指标）、髋旋转力矩（脑瘫内旋步态）。
- 正常步态的额状面力矩揭示一个反直觉事实：膝关节的外展力矩主要来自骨与韧带的被动承载，不是肌肉——说明人体把一部分稳定工作”外包”给了被动结构，这正是节能设计。
- 事件检测 + 时间归一化让”步态周期”成为可跨人、跨试验对齐的标准坐标，一切组间比较（病理 vs 健康）都建立其上。
5. How It Is Measured（如何测量）
- 空间标定（L 形架/3D 标定架对测力台）→ 建立 GRS；
- 贴标记：每环节 ≥ 3 个不共线跟踪标记（必要时加簇板），另贴内侧校准标记；
- 静态校准（解剖位或舒适站立位 ~1 s 平均）→ 算恒定 [M to A] 与质心偏移矢量 c；
- 拆除校准标记，动态采集；
- 逐帧：三标记 → [G to M] → [G to A] = [M to A][G to M] → 反解 θ₁θ₂θ₃（底行入手 + 回代验证）；
- 数值微分 + 式 7.7b → ω、α（解剖系）；
- 测力台 + ΣF = ma（GRS）→ 反作用力 → [G to A] 变换 → 欧拉方程解近端力矩 → [A to G] 回传 GRS 供上一环节使用；
- P = ΣMiωi 得三维功率；
- 事件检测（力阈值或运动学法）→ 时间归一化 → 按体重归一化力矩/功率 → 与常模比较。
6. Typical Applications（典型应用）
- 临床步态分析标准流程：脑瘫儿童手术决策（Davis 1991；Õunpuu 1996 是该领域操作蓝本）。
- 内侧间室膝骨关节炎：额状面膝内收力矩（KAM）是疾病进展的核心力学预测指标。[补充说明：本章给的”膝外展内力矩”即 KAM 的平衡方。]
- 跑步与跳跃的 3D 分析：着陆时横断面/额状面力矩与 ACL 损伤风险关联（M27/M28 论文的方法学基础）。
- 辅助器具研究：拐杖/助行器分担体重会改变 COP 与力台数据，研究器具效应需改用体重背心/悬吊替代支撑。
- 器械惯性：高尔夫挥杆必须连球杆一起建模——不规则器械需先实测惯性参数才能进模型。
7. Clinical Relevance（康复临床相关性）
- 校准体位妥协：患者摆不了标准解剖位是常态而非例外——“舒适站立位 + 维度测量推算关节中心”是临床可行的标准替代，但其关节中心定位误差大于研究级流程，解读病理数据时要知道精度底数。
- 事件检测失配：对足下垂（先足尖着地）、蹲伏步态患者用力阈值法可能完全检测不到 HC；对拖地步态运动学法也可能失效——事件错了，整段归一化曲线全错，这是临床报告里最常见的隐形错误源。
- 软组织伪影：簇板下滑、肥胖患者标记随软组织移动——标记系不再刚体，[M to A] 恒定假设被破坏，这是 3D 分析最大误差源之一。[补充说明：教材 7.2.1 的布套下滑警告即此。]
- 额状面信号：Trendelenburg 步态（髋外展肌无力）的力学本质是 H1-F 生成不足；内侧膝 OA 的力学本质是人体”被动外包”给内侧髁的负荷过大。
8. Performance Relevance（运动表现相关性）
- 鞭打类动作（投掷、挥拍、踢球）的速度来自环节链诱导加速度的累积效应，单关节力量训练无法替代链式协调训练。
- 高分辨率小标记 + 近距离相机使表情捕捉、手部精细动作分析成为可能；高速项目（短跑、跳跃）则要权衡标记尺寸与脱落风险。
- 6DOF 标记组可分辨横断面微小旋转，为短跑/变向中膝踝旋转负荷监测提供精度基础（Collins et al., 2009：重复性优于 Helen Hayes）。
9. Possible Interpretations（表现型 → 机制推理）
| 观察（表现型信号） | B 类：测量/技术因素 | C 类：功能性因素 | D 类：结构性/病理性因素 |
|---|---|---|---|
| 支撑期检测不到 HC | 力台未清零/漂移；脚未完全落台；20 N 阈值对低体重者过高 | 步速极慢；扶助器具分担负荷 | 足下垂/马蹄足（无脚跟着地模式） |
| 关节角曲线出现不合理跳变 | 标记重叠被合并质心；标记脱落重贴；万向锁附近转序退化 | 极端关节活动度（近活动边界） | 关节半脱位等（罕见，需影像确认） |
| 横断面旋转角噪声大 | 标记共线度差/标记间距太小；Helen Hayes 类稀疏标记组 | 步态变异性大 | 胫骨扭转、股骨前倾等骨性异常 |
| 额状面髋力矩减小 | 簇板下滑致 [M to A] 失效；步宽改变未记录 | 步速下降；躯干侧倾代偿（减小重力力臂） | 髋外展肌无力（Trendelenburg）；髋内翻 |
| 左右不对称的 3D 力矩 | 只标定了一侧/两侧标记组不对称 | 疼痛性避荷；单手杖使用 | 单侧关节病、下肢不等长、偏瘫 |
10. What It Cannot Tell Us（局限性）
- 皮肤不是骨头：所有 3D 分析建立在”标记相对骨骼刚性固定”的假设上，软组织伪影无法由模型本身消除。
- 欧拉角依赖转序：同一运动用不同转序得到不同角度曲线；跨研究比较必须核对转序与 ISB 建议的一致性。
- 净力矩依然看不到共收缩与被动结构贡献的分解：膝额状面力矩是骨/韧带/肌肉的混合账，逆向动力学无法拆分。
- 关节中心是推算的：尤其髋关节中心，误差直接传导为力矩误差（力臂错误）。
- 事件算法不是普适的：每种算法都有适用边界（步速、病种、场地），教材明确警告跑步与非典型步态。
- 2D 与 3D 结果不可直接混比：耦合项与平面外运动的差异意味着同一人 2D/3D 分析结论可能不同。
11. Related Metrics（相关概念双链）
- 上游基础：M21_力台与链段模型（2D 版方程）、M19b_平滑与微分问题（微分/滤波直接影响 θ̇、α 质量）、M19_运动学测量系统对比（标记/相机硬件）
- 横向衔接：M22_关节功率与能量流（3D 功率 = 三分量点积）
- 下游应用：M24_肌电信号与肌肉激活（净力矩的共收缩盲区）、M25_肌骨建模与平衡控制（诱导加速度 → 正向动力学）、M26_步态综合案例
- 论文延伸：M27_Bahr2005_损伤机制、M28_Krosshaug2007_落地机制（3D 额状面/横断面力矩是损伤机制研究的核心变量）
12. Evidence（证据与出处）
- 坐标系定义、Cardan 转序与反解算例：Winter 第 7 章 7.1–7.2 节（书页 145–158），Table 7.1–7.2。
- 标记组比较（Helen Hayes / Cleveland Clinic / 6DOF）：Collins et al. (2009), Gait Posture 30:173–180。
- 运动学事件检测算法：Zeni, Richards & Higginson (2008), Gait Posture 27:710–714。
- 临床实验室流程：Davis et al. (1991), Hum. Mov. Sci. 10:575–587；Õunpuu et al. (1996), Gait Posture 4:62–78。
- ISB 关节坐标系建议：Wu et al. (2002, Part 1), J. Biomech. 35:543–548；Wu et al. (2005, Part 2), J. Biomech. 38:981–992（Table 7.3–7.4）。
- 3D 牛顿-欧拉算例（frame 6）：第 7 章 7.4.1–7.4.3 节（书页 162–166），Table 7.5–7.6。
- 正常步态 3D 力矩/功率常模曲线：Eng & Winter (1995), J. Biomech. 28:753–758（图 7.5–7.6）。
13. My Understanding（我的理解，5 点）
- 3D 分析 = 一本固定词典 + 每帧一次翻译。[M to A] 恒定、[G to M] 时变，两者相乘得 [G to A]——这条主线的优雅之处在于把”难算的”（解剖系定义）压缩进一次校准，把”好算的”（标记坐标）留给每帧。
- 反解角度从矩阵底行入手，且必须回代验证：s₂、−s₁c₂、c₁c₂ 三个元素锁定 θ₁θ₂，再验 θ₃。每个角两个候选解——不验证就有 1/4 概率拿错解，曲线会莫名翻转。
- 事件检测是”地基中的地基”：事件错 → 归一化错 → 组间比较全错。而且健康人验证过的算法到病理人群可能系统性失效——我的系统在做表现型比较前，必须先有一层”事件有效性自检”。
- 额状面与横断面才是临床金矿：矢状面看功能（推离、承重），额状面看稳定（外展力矩、KAM），横断面看旋转控制。且膝外展力矩是被动结构的贡献——“内力矩 ≠ 肌肉力矩”，这对机制推理框架是重要警示。
- 角速度不是 Cardan 角导数：必须经式 7.7b 投影。直接对 θ₁θ₂θ₃ 求导再当 ω 用是初学者最易犯的隐形错误——三轴不共轴时误差巨大。任何信号→特征管线若扩展到 3D，这一步都必须单独实现。
模块自测（7 题）
Q1 为什么每个体节至少需要 3 个不共线的跟踪标记点？如果三个标记共线，会丢失什么信息？
参考答案：3 个不共线点才能唯一确定 3D 空间中的一个平面（坐标系需要两个面内轴 + 一个法向轴）。若三点共线，则只能确定一条直线（一个轴的方向），环节绕该直线（长轴）的旋转完全不可测——即丢失横断面（轴向旋转）自由度，6DOF 退化为不足。教材还给出实践推论：保证至少 3 个跟踪标记后，内侧校准标记可在动态采集前拆除，避免相互碰撞。
Q2 用 Table 7.1 的校准数据验证：zan = (3.855, −1.720, 31.730)，临时 yan = (−0.210, 7.670, 0.120)。计算 xan = yan × zan（保留一位小数），并说明为什么还要再算一次 yan = xan × zan。
参考答案：叉积展开：xx = (7.670)(31.730) − (0.120)(−1.720) ≈ 243.4 + 0.2 = 243.6；yx = (0.120)(3.855) − (−0.210)(31.730) ≈ 0.5 + 6.7 = 7.1；zx = (−0.210)(−1.720) − (7.670)(3.855) ≈ 0.4 − 29.6 = −29.2（教材值 246.7 / 7.1 / −29.6，微小差异来自教材用未舍入的关节中心坐标）。再算 yan = xan × zan 的原因：内外踝连线只是近似垂直于长轴，直接用它建系会得到非正交坐标系；用 xan × zan 重算 yan 可强制三轴严格正交（右手系），这是构造合法旋转矩阵的前提。
Q3 用算例 7.4.3 的数据求 frame 5 的全局近端反作用力（小腿 m = 3.22 kg，g = 9.814 m/s²；aX = 5.89、aY = 1.30、aZ = −1.66 m/s²；RXD = −102.71、RYD = −12.12、RZD = −790.27 N）。
参考答案：Step 1 三式： RXP = maX + RXD = 3.22 × 5.89 + (−102.71) = 18.97 − 102.71 = −83.74 N； RYP = maY + RYD + mg = 3.22 × 1.30 + (−12.12) + 3.22 × 9.814 = 4.19 − 12.12 + 31.60 = 23.67 N； RZP = maZ + RZD = 3.22 × (−1.66) + (−790.27) = −5.35 − 790.27 = −795.62 N。
Q4 一位脑瘫患儿步态分析中，测力台 20 N 阈值法始终检测不到可靠的 HC。请给出两种运动学替代方案（分别适用于平地与跑台），并指出各自的运动学判据与已知局限。
参考答案（Zeni et al., 2008）：① 平地：HC = 足跟标记与骶骨标记前后向相对位移最大的时刻；TO = 足尖标记与骶骨标记前后向相对位移最大时刻。② 跑台：HC = 足跟标记前后向速度由正转负的时刻；TO = 由负转正的时刻。局限：两种方法均在正常步速步行上验证，跑步时可能不准；且运动学法本身对非标准步态模式敏感——对严重异常步态需人工核验事件时刻。
Q5 什么是万向锁（gimbal lock）？为什么它在分析上肢大活动度运动（如投掷）时比分析步行时更容易成为问题？
参考答案：万向锁是欧拉角参数化的固有奇点：当中间那次旋转把第一个与第三个转轴转到互相平行时，两个旋转自由度退化重合，3D 姿态无法用三个角唯一表示，角度解出现奇异跳变。步行中下肢关节基本在中立位附近活动，中间角远离奇点；而上肢（肩、躯干）活动范围大且多平面耦合，极易进入中间角接近 ±90° 的区域——因此 ISB 对上肢各环节分别推荐了特定转序（Table 7.4，如肩 Y–X–Y），把奇点”推到”生理活动范围之外。
Q6 教材报告 frame 6 的膝关节近端力矩 Mxp = +307.78 N·m，并将其解读为外展力矩。请写出欧拉方程式 7.8a 中该力矩的三个主要贡献项（用算例数值），并评价该量级的合理性。
参考答案：Mxp = Ixαx + (Iz−Iy)ωyωz − Rzd·ld − Rzp·lp + Mxd，代入：惯性项 ≈ 0.0138 × 3.066 ≈ 0.04 N·m（可忽略）；−Rzd·ld = −(−781.336)(0.1386) ≈ +108.3 N·m；−Rzp·lp = −(−777.674)(0.1815) ≈ +141.2 N·m；+Mxd = +58.3 N·m，合计 ≈ 307.8 N·m。可见力矩几乎全部由反作用力 × 力臂贡献。物理解读为正外展力矩：对抗上身体重作用于支撑髋内侧约 10 cm 处产生的内收重力矩。但 308 N·m 的量级远超正常膝外展力矩常模（数十 N·m 量级），教材原文数据疑有排版/单位错误，复算时应以原始表 7.5–7.6 数据自行验证——这也提醒：逆向动力学结果必须过一遍”量级合理性”检查。
Q7 正常步态额状面膝关节存在内源性外展力矩，但它并非肌肉收缩的产物。请解释其真实来源、功能意义，以及这一事实对假肢/矫形器设计的启示。
参考答案：真实来源是被动承载：承重塑膝在重力负荷（作用于支撑膝内侧）下试图内翻，内侧髁软骨/骨加载增大、外侧髁卸载，内外侧负荷差形成被动外展力矩；韧带结构也参与约束。功能意义：人体把额状面稳定的一部分”外包”给被动结构，节省了肌肉代谢成本。启示：① 逆向动力学给出的净力矩是肌肉+骨+韧带的混合账，不能全算在肌肉头上；② 假肢膝关节内的弹簧、阻尼、机械限位同样会产生内力矩——设计者的”助力”若方向不当会变成肌肉必须对抗的”阻力”，康复工程设计必须显式建模这些被动元件力矩。
上一模块：M22_关节功率与能量流 ｜ 下一模块：M24_肌电信号与肌肉激活