运动生物力学 × 功能表型知识库

运动生物力学 × 功能表型知识库

肌电信号与肌肉激活（Kinesiological Electromyography）

M24 · 肌电信号与肌肉激活（EMG and Muscle Activation） 阶段五 · 动力学与能量学 前面章节的净关节力矩（M21–M23）有一个致命盲区：看不到主动肌与拮抗肌的同时发力。本章补上这块拼图——EMG 是唯一能非侵入式观察”神经向肌肉下达了什么指令”的窗口。核心问题链：EMG 是什么（电生理）→ 怎么录干净（放大器四大指标）→ 怎么处理成可用曲线（整流/包络/积分）→ 能推断什么、不能推断什么（EMG-张力关系）。

0. 大白话定义

肌肉收缩时，肌纤维上会跑过一波波”电流涟漪”（运动单元动作电位，m.u.a.p.）。皮肤表面贴的电极就像放在湖面的浮标，测到的是附近所有涟漪叠加后的水面起伏——这就是 EMG。

几个关键直觉：

  • 电极只听得见附近的肌肉：小运动单元的”嗓门”只够传 0.5 cm，最大的也只传 1.5 cm。所以表面电极基本只反映它正下方那块浅表肌肉——隔一块肌肉”偷听”（串扰）虽存在，但比你想象的弱。
  • EMG 是”命令强度”不是”力量大小”：它反映神经激活了多少运动单元、放得多快，而肌腱上真实的力量要等 40–100 ms 才跟上（像按下开关后灯慢慢变亮）。
  • 离心收缩时”花小钱办大事”：同样的力，离心收缩的 EMG 明显更小——这就是为什么 EMG 不能直接当”力计”用。
  • 疲劳的指纹：EMG 振幅变大 + 频率变低（声音变”闷”），是肌肉疲劳在信号上的双重签名。

1. Definition 术语表（中英对照）

英文术语中文一句话定义
Electromyogram (EMG)肌电图伴随肌肉收缩的电信号，电极记录到的是所有 m.u.a.p. 的代数和
Motor Unit Action Potential (m.u.a.p.)运动单元动作电位一个运动单元被募集后沿其肌纤维传播的去极化波
Motor End Plate运动终板运动神经元与肌纤维间的突触接头（释放 ACh）
End Plate Potential (EPP)终板电位终板处的突触后膜去极化，达阈值才触发肌纤维动作电位
Twitch单收缩单个动作电位引起的收缩元冲量式张力响应（峰 50–110 ms，全程可达 300 ms）
Surface / Indwelling Electrode表面 / 肌内电极贴皮肤测浅表肌肉平均活动 / 针或细丝插入测深部或精细肌肉
Cross-talk串扰电极拾取到邻近肌肉 m.u.a.p. 的重叠区信号
Common-Mode Rejection Ratio (CMRR)共模抑制比差分放大器抑制工频干扰等共模信号的能力，优质 ≥ 80 dB
Full-Wave Rectification全波整流取 EMG 绝对值，使信号带随收缩强度波动的偏置
Linear Envelope线性包络整流后低通滤波得到的移动平均曲线，形状近似张力曲线
Median / Mean Power Frequency中位 / 平均功率频率功率谱的中位/加权平均频率，疲劳时下降
MVC最大随意收缩归一化 EMG 振幅的常用参照

2. 核心框架与公式

2.1 电生理链条：从神经冲动到 EMG

运动神经元（“最终共同通路”）→ 终板释放 ACh → 跨过 200–500 Å 突触间隙 → EPP 达阈值 → 肌纤维动作电位 → 沿横管系统内传 → 肌质网释放 Ca²⁺ → ATP 水解 → 横桥产生冲量力（twitch）。

  • 电极看到的：去极化波前 + 复极化波沿肌纤维传播形成的移动偶极子（dipole：源 +I 与汇 −I 相距 b = 0.5–2.0 mm）。
  • 单电极记录为双相波（波前接近时电位升至最大，恰过电极下方时过零，远离时变负）；双电极差分后变三相波。
  • 偶极子场强按 1/r 衰减（式 9.1–9.2）：拾取距离 = 小运动单元（50 纤维）约 0.5 cm，最大运动单元（2500 纤维）约 1.5 cm（Fuglevand et al., 1992）。
  • 电极面积越大，记录的 m.u.a.p. 时程越长、振幅被平均拉低 → 电极越小越好，且抗串扰更强。
  • 传导速度正常人约 4 m/s；疲劳与肌病时下降 → m.u.a.p. 时程变长 → 频谱向低频移动（1 mm 距离 MPF ≈ 160 Hz，20 mm 处仅 25 Hz）。
  • 神经肌肉传递障碍（如 ACh 耗竭）时，神经冲动与 m.u.a.p. 不再一一对应——终板阻滞可为完全性、高频性或间歇性。

2.2 放大器四大指标（9.2）

  1. 增益与动态范围：表面 EMG 最大约 5 mV p-p（MVC），肌内可达 10 mV；放大器噪声应 < 50 μV（优选 20 μV）；增益可选 100–10,000；报告振幅时必须折算回电极处的 mV，与增益无关。削顶（overdriving）是最常见失真。
  2. 输入阻抗：电极-皮肤阻抗经良好备皮可降至 ≤ 1000 Ω，放大器输入阻抗应 ≥ 1 MΩ（肌内电极阻抗可达 50 kΩ → 需 ≥ 5 MΩ）。算例：Rs = 10 kΩ、Ri = 80 kΩ 时 2 mV 信号只剩 1.6 mV；备皮 + 1 MΩ 后保住 1.998 mV。
  3. 频率响应：EMG 谱 5–2000 Hz，能量集中在 20–200 Hz。推荐带宽：表面 10–1000 Hz，肌内 20–2000 Hz；模式识别需 5/10 kHz 上限。运动伪影在 0–10 Hz，工频干扰 60/50 Hz（陷波会连同该频段真实 EMG 一起去掉）。截止频率处增益降至 0.707（半功率点，−3 dB；dB = 20 log₁₀ 增益，式 9.3）。
  4. 共模抑制：人体是天线，工频 hum 可达 100 mV（淹没最大 EMG）。差分放大器做减法消除共模信号（式 9.4）；残余 hum = A × Vhum / CMRR（式 9.5）。算例 9.1：EMG 2 mV、hum 500 mV、CMRR 10,000:1、增益 2000 → 输出 EMG = 4 V，输出 hum = 100 mV。优质放大器 CMRR ≥ 80 dB（= 10,000:1，式 9.6）。

2.3 串扰的检验与量化（9.2.5）

  • 手法抵抗测试：实验前让受试者分别孤立收缩目标肌与邻肌（如外翻查 PER→TA/LG 串扰），可借助视觉反馈；详见 Winter et al. (1994)。
  • 互相关量化（式 9.7）：等长收缩约 10 s，Rxy(τ) 在 ±30 ms 内扫描；峰值在零相移处。Rxy² = 共享信号比例。实测（电极对平行肌纤维、等距运动点）：2 cm 间距 Rxy = 0.6 → 36% 串扰；2.5 cm → 23%；5 cm → 6%；7.5 cm → 2%。
  • 远处运动单元贡献主要串扰且其 m.u.a.p. 低频 → 微分处理可抑制串扰（双峰微分技术 DeLuca 1997，但临床少用）。
  • 预防胜于处理：电极放肌腹正中（横断+矢状面），检测面沿肌纤维方向排列（[补充说明：教材原文”检测面垂直于肌纤维方向”指两检测点连线沿纤维方向，即双极电极沿肌纤维走向放置]），偏离肌腹中心会使运动单元池随收缩改变、靠近肌腱则代表性差；超声可辅助定位。SENIAM 给出 27 块肌肉的标准位点。

2.4 五种处理方式（9.3，图 9.17）

  1. 全波整流：|EMG|，半定量看时相活动；单位仍为 mV。
  2. 线性包络：整流 + 低通 = 移动平均，形状近似张力曲线。不要叫 IEMG（与数学积分混淆）。
  • 生理依据：twitch 是二阶临界/微过阻尼系统，截止频率 2.3–7.8 Hz；fc = 1/(2πT)（式 9.8，T = twitch 达峰时间）。比目鱼肌 T ≈ 106 ms → fc = 1.5 Hz。
  1. 全收缩积分：曲线下面积，单位 mV·s；除以时间得真实平均 EMG（mV）。
  2. 定时重置积分（40–200 ms）：一串峰 = 移动平均；重置太快引入噪声，太慢跟不上波动。
  3. 定压重置积分：重置脉冲频率 ∝ 收缩强度（类似神经放电率），神经生理学家偏爱。

2.5 EMG 与生物力学变量的关系（9.4）

  • 等长张力：线性（Lippold 1952，小腿）与非线性（Zuniga & Simons 1969；Vredenbregt & Rau 1973，肘屈肌）都有报道——本质是特定长度下的静态标定。AP 计数率与张力近似线性（Close et al., 1960）。
  • 动态张力滞后 EMG 40–100 ms（twitch 达峰时间所致；图 9.20）；二阶低通模型可拟合（Gottlieb & Agarwal, 1971）。
  • 缩短与拉长（Komi 1973, 1987）：向心时张力降而 EMG 基本不变；离心时张力升而 EMG 反而更低 → EMG 反映的是收缩元的激活状态，不是肌腱张力；且负功的代谢成本低于正功（与 M22 呼应）。
  • 疲劳（9.4.3）：等张力维持需追加募集 → EMG 振幅升；频谱左移三机制（传导速度降、快大运动单元退出、运动单元同步化出现 8–10 Hz 震颤）。代表性数据（Krogh-Lund & Jørgensen, 1991）：传导速度 −10%、中位频率 −45%、RMS 振幅 +250%。
  • 中位频率 fm（式 9.9）与平均频率 MPF（式 9.10）几乎等价（Kerr & Callaghan, 1999）。
  • 三大悬而未决（动态多肌共跨关节时）：① 长度变化如何改变 EMG-张力关系；② 协同肌如何分担负荷；③ 拮抗肌活动引入的未知力 → 结论：经过适当标定的线性包络 EMG 只能作为长度变化不快肌肉的”粗粒度张力预测器”。

3. Physical Meaning（物理意义）

  • 体积导体中的移动偶极子：场强 1/r 衰减决定了”拾取半径”这个 EMG 一切应用的空间分辨率上限。
  • 差分放大 = 用空间对称性对抗电磁环境：共模信号（天线效应）在两端等幅，减法即消；差模信号（两电极下不同的肌肉活动）被保留放大。
  • 整流 + 低通 = 把”神经脉冲序列”卷积成”收缩元张力响应”的物理模拟——滤波器截止频率应匹配目标肌肉的 twitch 时间，这是包络处理唯一的生理学依据。
  • 频谱左移 = 传导速度下降 + 时程变宽的傅里叶对偶表现。

4. Human Movement Meaning（人体运动意义）

  • EMG 揭示了净力矩的盲区：步行中胫前肌与比目鱼肌存在可量化的共收缩（M22 的 COCON 指数即建立在 EMG 上）。
  • 肌肉激活模式是时相性的：正常步态中每块肌肉有精确的开关机时刻，病理状态下时相紊乱（延长、提前、持续激活）往往比振幅异常更有诊断价值。
  • 离心/向心的 EMG 不对称性解释了为什么”省力”的动作模式（如下坡走依赖离心刹车）代谢成本低于直觉。

5. How It Is Measured（如何测量）

  1. 选电极：浅表大肌肉 → Ag/AgCl 表面电极（直径约 1 cm）；深部肌肉（如骨间肌、髂腰肌、vastus intermedius、adductor magnus）→ 细 wire 肌内电极；
  2. 备皮至阻抗 ≤ 1000 Ω，肌腹正中放置、沿肌纤维方向，SENIAM 标准位点；
  3. 放大器：增益 100–10,000、输入阻抗 ≥ 1 MΩ、带宽 10–1000 Hz、CMRR ≥ 80 dB；
  4. 手法抵抗测试 + 互相关（Rxy²）排查串扰；
  5. 采样率 ≥ 2× 上限频率（表面 ≥ 2000 Hz）；
  6. 处理：全波整流 → 线性包络（fc 按目标肌 twitch 时间选 2–8 Hz）或积分；
  7. 归一化：MVC 或动态参照（[补充说明：教材强调振幅以电极处 mV 报告；MVC 归一化在动态任务中的适用性有争议，需注明参照方式]）；
  8. 疲劳分析：中位/平均功率频率 + RMS 振幅。

6. Typical Applications（典型应用）

  • 步态肌肉时相分析：判断某肌”何时开、何时关”，与力矩/功率曲线联合解释（如 A2-S 推离爆发对应腓肠-比目鱼 EMG 高峰）。
  • 共收缩量化：COCON 指数（M22）的输入信号。
  • 疲劳监测：劳动工效、运动训练中频谱左移作为客观疲劳指标。
  • 生物反馈：让患者看着自己的 EMG 学习孤立收缩（手法抵抗测试中也用作辅助）。
  • 运动单元分解：四极针电极 + 模板匹配可追踪单个运动单元 0–100% MVC 的放电率（DeLuca 1993；Erim et al., 1996）。

7. Clinical Relevance（康复临床相关性）

  • 痉挛/挛缩鉴别：被动牵伸时的 EMG 反射性升高提示神经性成分，无 EMG 变化的僵硬提示结构性挛缩。[补充说明：这是临床常用推理，教材 9.0 提及反射活动对信号的影响。]
  • 病理步态时相紊乱：脑瘫患儿常见的”持续激活”（该关时不关）是步态能量低效的可测来源，与 M22 四大低效机制中的共收缩直接对应。
  • 终板阻滞类疾病（如重症肌无力）：神经-肌肉传递不一一对应，EMG 解读必须考虑传递障碍。
  • 疲劳评估：低频成分增加 + 振幅增加的双重签名可用于下背痛患者竖脊肌耐力评估等场景。
  • 串扰误判的临床风险：把邻肌活动误判为目标肌激活，会直接导致错误的”哪块肌该练/该松”决策——手法抵抗测试应成为临床 EMG 的标配前奏。

8. Performance Relevance（运动表现相关性）

  • 力量训练中 EMG 振幅常被用作”激活程度”的代理指标比较不同动作——但受长度-张力与离心/向心不对称影响，跨动作比较需谨慎。
  • 频谱分析监测训练疲劳与恢复状态。
  • 短跑等高速运动中 EMG 时相分析可揭示预激活（触地前激活）策略——预激活是落地稳定的关键机制（与 M28 落地损伤论文衔接）。

9. Possible Interpretations（表现型 → 机制推理）

观察（表现型信号）B 类：测量/技术因素C 类：功能性因素D 类：结构性/病理性因素
目标肌 EMG”异常升高”串扰（邻肌同相活动）；电极偏离肌腹中心；增益设置变化拮抗共收缩（稳定策略）；负荷增加；疲劳致追加募集痉挛；肌张力障碍；代偿性过度使用
EMG 振幅低但力矩正常电极-皮肤阻抗高；放大器饱和削顶离心收缩（EMG 天然低于同力向心）；肌腱储能贡献大肌肉激活不足（神经驱动下降）——需排除前两类才能下此结论
频谱左移电极间距过大；低温真实疲劳（传导速度降 + 同步化）；血流受限肌病（传导速度病理性下降）
基线 hum 明显CMRR 不足；接地不良；附近电源/荧光灯——（几乎总是技术问题，先查设备再谈机制）——
支撑期全程”持续激活”包络滤波 fc 过低把时相抹平疼痛性保护性共收缩；新任务学习期脑瘫/卒中后的时相控制障碍；锥体外系强直

10. What It Cannot Tell Us（局限性）

  1. EMG ≠ 力：长度变化、收缩速度、收缩类型（离心/向心）、疲劳都使 EMG-张力关系漂移；只有”长度变化不快的肌肉经适当标定后可作粗预测”（教材原话）。
  2. 看不到单块深部肌肉（表面电极）：拾取半径 0.5–1.5 cm，深部肌必须肌内电极（有创）。
  3. 无法分解协同肌负荷分担：多肌共跨一关节时，各自贡献多少力是未知数（M25 建模的冗余性问题）。
  4. 振幅受一篮子非神经因素调制：电极位置、面积、间距、皮下脂肪厚度、串扰——跨人/跨日比较必须归一化且承认残余不确定性。
  5. 时程失真：电极面积与距离使 m.u.a.p. 变宽，频率信息系统性偏低。
  6. 不显示被动张力：挛缩、肌腱弹性贡献在 EMG 上是”静音”的。

12. Evidence（证据与出处）

  • 拾取距离与电极几何影响：Fuglevand, Winter, Patla & Stashuk (1992), Biol. Cybern.（图 9.3、9.5；书页 192–195）。
  • 串扰量化（互相关法与距离梯度 36%→2%）：Winter, Fuglevand & Archer (1994), J. Electromyogr. Kinesiol. 4:15–26。
  • 电极放置与报告标准：SENIAM（Hermens et al., 2000, J. Electromyogr. Kinesiol. 10:361–374）；ISEK Ad Hoc 报告（Winter et al., 1980）。
  • 线性包络的 twitch 匹配依据：Crosby (1978)；Milner-Brown et al. (1973)（截止 2.3–7.8 Hz）。
  • 离心收缩 EMG 更低：Komi (1973)；Komi, Kaneko & Aura (1987)。
  • 疲劳频谱三机制与代表数据：Kadefors et al. (1973)；Krogh-Lund & Jørgensen (1991)（传导速度 −10%、中位频率 −45%、RMS +250%）。
  • 运动单元分解：DeLuca (1993)；Erim et al. (1996)。

13. My Understanding（我的理解，5 点）

  1. EMG 是”神经指令”的模糊传真：它忠实地反映激活状态，但从激活到肌腱力之间隔着 twitch 滤波、长度-张力、力-速度三重变换——把 EMG 包络直接当力曲线用，是康复文献中最常见的过度解读。
  2. “录得干净”是一门工程学：阻抗、CMRR、带宽、电极几何，每一项都有可量化的验收标准（≤1 kΩ、≥80 dB、10–1000 Hz、面积最小化）。我的系统在接收 EMG 数据时应先跑一遍”信号质量清单”。
  3. 串扰可测、可防、不可假设不存在：Rxy² 给出共享比例，7.5 cm 间距才降到 2%。小肌肉群（前臂、小腿）几乎必然有串扰，结论措辞必须留余地。
  4. 疲劳的双签名（振幅↑ + 频率↓）是可计算的：中位频率与 RMS 都是现成公式，可直接成为功能表现型中的”耐力维度”指标。
  5. EMG + 力矩 + 功率三位一体才是完整画像：力矩说”净效果”，功率说”能量流向”，EMG 说”谁在发力、何时发力、有没有内耗”——这正是机制推理框架中区分 C 类（共收缩策略）与 D 类（时相控制障碍）的关键数据组合。

模块自测（7 题）

Q1 为什么表面电极记录的是双相波，而双电极差分后变成三相波？拾取距离大约是多少？

参考答案：m.u.a.p. 是沿肌纤维传播的偶极子（去极化源 +I 在前、复极化汇 −I 在后）。波前接近单电极时 r₁ < r₂，电位为正并升至最大；偶极子恰在电极正下方时 r₁ = r₂，两贡献抵消过零；远离时 r₁ > r₂，电位转负——形成双相波。双电极各记录一个时间上略有错位的双相波，差分相减后正负峰错位叠加，形成三相波。拾取距离：小运动单元（50 纤维）约 0.5 cm，最大运动单元（2500 纤维）约 1.5 cm（Fuglevand et al., 1992）。

Q2 受试者皮肤 EMG 为 2 mV，工频 hum 达 500 mV。放大器增益 2000、CMRR = 10,000:1。求输出端的 EMG 与残余 hum 各为多少？该 CMRR 换算成分贝是否达到优质放大器标准？

参考答案：输出 EMG = A × (emg₁ − emg₂) = 2000 × 2 mV = 4 V；残余 hum = A × Vhum / CMRR = 2000 × 500 mV / 10,000 = 100 mV。CMRR(dB) = 20 log₁₀(10,000) = 80 dB——恰好达到优质生物放大器 ≥ 80 dB 的下限标准（算例 9.1）。

Q3 电极-皮肤阻抗 Rs₁ = Rs₂ = 10 kΩ，放大器输入阻抗 Ri = 80 kΩ，EMG 信号 2 mV。输入端实际得到多少？若备皮使 Rs 降至 1 kΩ 且 Ri 升至 1 MΩ，又得到多少？这说明了什么？

参考答案：分压比 = Ri/(Ri + Rs₁ + Rs₂)。前者 = 80/(80 + 20) = 0.8 → 2 × 0.8 = 1.6 mV（两个电极各损失 0.2 mV）。后者 = 1000/(1000 + 2) ≈ 0.998 → 1.998 mV。说明：备皮（降电极阻抗）与高输入阻抗同等重要——目标指标为皮肤阻抗 ≤ 1000 Ω、输入阻抗 ≥ 1 MΩ（肌内电极需 ≥ 5 MΩ）。

Q4 比目鱼肌的 twitch 达峰时间约 106 ms。若要让线性包络在物理上模拟该肌的张力曲线，临界阻尼二阶低通滤波器的截止频率应设为多少？写出依据公式。

参考答案：依据 twitch 的二阶冲激响应匹配关系 fc = 1/(2πT)（式 9.8）：fc = 1/(2π × 0.106) ≈ 1.5 Hz。这就是”包络滤波不是随便选 fc”的原因——它应对应目标肌肉 twitch 的力学响应（文献报道范围 2.3–7.8 Hz，快肌取高端、慢肌取低端）。

Q5 两名相邻肌肉（电极间距 2 cm）的互相关峰值 Rxy = 0.6（零相移处）。两通道信号有多大比例是共同的？若把电极间距加大到 5 cm，串扰约降到多少？

参考答案：共享方差比例 = Rxy² = 0.6² = 0.36，即 36% 串扰。间距 5 cm 时 Rxy ≈ 0.24 → Rxy² ≈ 0.06，串扰约 6%（Winter et al., 1994 实测梯度：2 cm→36%、2.5 cm→23%、5 cm→6%、7.5 cm→2%）。零相移出现峰值是因为两电极对平行肌纤维且等距运动点，共同 m.u.a.p. 几乎同相到达。

Q6 受试者在等速肌力测试仪上分别以同等主观努力做最大向心与离心收缩。预期 EMG 振幅与肌腱张力在两种收缩中各如何表现？这一不对称对”用 EMG 估计肌肉力”意味着什么？

参考答案（Komi 1973, 1987）：向心收缩时张力下降但 EMG 振幅基本不变；离心收缩时张力升高但 EMG 振幅反而低于同等正功水平。说明 EMG 反映的是收缩元的激活状态而非肌腱张力——离心时单个运动单元可产生更大张力（力-速度关系），用较少激活即可达成大力。含义：任何 EMG→力的换算都必须按收缩类型分开标定；跨向心/离心比较 EMG 振幅没有意义；同时这也与负功代谢成本低于正功互为印证。

Q7 一名腰痛患者竖脊肌耐力测试中，记录到 EMG 中位频率下降 45%、RMS 振幅上升 250%、传导速度下降 10%，且信号中出现 8–10 Hz 节律。请逐一解释这四个现象的机制，并说明为何”振幅上升”不能解读为”肌肉变强了”。

参考答案（Krogh-Lund & Jørgensen, 1991 为代表性数据）：① 中位频率 ↓45%——疲劳使 m.u.a.p. 传导速度下降、时程变宽，频谱整体左移；② 振幅 ↑250%——为维持目标张力，神经系统追加募集新运动单元以补偿已募集单元放电率与峰张力下降（同时同步化也推高振幅）；③ 传导速度 ↓10%——肌内代谢环境改变（缺血/底物耗竭）直接减慢膜传导；④ 8–10 Hz 节律——疲劳下运动单元由独立随机放电转为同步成簇放电的神经性震颤，在 EMG 与力记录上均可见。“振幅上升”恰恰反映的是效率下降（需要更多电活动才能维持同样张力），是疲劳而非增强的证据。

上一模块：M23_3D变量与事件检测 ｜ 下一模块：M25_肌骨建模与平衡控制

反向链接

本页提及