try ai
科普
编辑
分享
反馈
  • 肌电信号处理:原理与应用

肌电信号处理:原理与应用

SciencePedia玻尔百科
核心要点
  • 通过整流和低通滤波处理原始肌电信号,是估算“神经驱动”的基本方法,而神经驱动是衡量向肌肉发出的神经指令的指标。
  • 恰当的滤波是一项关键的权衡,需要在去除运动伪影和工频干扰等噪声的同时,控制好时间延迟和信号失真等后果。
  • 高密度肌电图结合空间滤波等先进技术,可以增强聚焦,分离特定运动单元的活动,从而提高信号的特异性。
  • 肌电分析是一种功能多样的工具,具有深远的应用价值,从诊断运动障碍、引导外科医生,到个性化康复治疗和分析生物力学功能。

引言

我们肌肉的电活动,即肌电图 (EMG),就像体育场馆人群的呐喊——一个复杂、混沌但信息丰富的信号。虽然原始信号本身看起来充满噪声,但它包含了来自神经系统的精确指令,这些指令调控着我们的一举一动。研究人员和临床医生面临的核心挑战是如何解码这曲电气交响乐。我们如何将一个原始的、波动的波形,转化为一个清晰、可靠的衡量神经意图和肌肉功能的指标?

本文通过提供一份全面的肌电信号处理指南来应对这一挑战。首先,在“原理与机制”部分,我们将深入探讨采集、清理和解读肌电信号的基本步骤,探索从采样、滤波到时频分析等核心概念。随后,“应用与跨学科联系”一章将展示这些技术在现实世界中的应用,它们如何为临床神经学、外科学、康复学和生物力学等领域带来革命性的变化。读完本文,您将了解如何将身体的电气低语转化为关于运动、意图和健康的、有意义的故事。

原理与机制

人体的电气交响乐

想象一下,你正站在一个拥挤的体育场里。人群爆发出欢呼声。你无法分辨出任何一个人的声音,但那起伏的呐喊声却向你讲述着比赛的戏剧性时刻。我们从肌肉中测量到的电信号,即​​肌电图 (EMG)​​,与那呐喊声非常相似。它不是一个单一、干净的音符,而是由成千上万个微小的电脉冲叠加而成的混沌、复杂的信号。这些脉冲来自单个​​运动单元​​——执行我们大脑命令的神经和肌纤维的基本组合。

当你在屏幕上观察原始肌电信号时,它是一个在零点附近波动的、充满噪声且快速变化的波形。它的平均值是零,但其噼啪作响的纹理却蕴含着丰富的信息。这是一个干涉图样,是记录电极附近无数个运动单元放电时发出的微小电“呼喊”的总和。我们的首要任务是在不失真的情况下捕捉这曲交响乐。

我们用皮肤上的电极捕捉它,在离散的时间点上对电压进行快照——这个过程称为​​采样​​。但采样有一条基本规则,有点像制作电影的规则。如果你给一个旋转的车轮拍照的速度太慢,它可能会看起来像在倒转,甚至静止不动。这种高频伪装成低频的错觉,被称为​​混叠​​。为了避免这种情况,我们必须以至少两倍于信号中最高频率的速度进行采样——这个极限被称为​​奈奎斯特准则 (Nyquist criterion)​​。

这个概念不仅限于时间。借助现代​​高密度肌电图 (HD-EMG)​​,我们使用由数十或数百个微小电极组成的网格来创建肌肉活动的空间图谱。在这里,同样的原则也适用。如果我们的电极间距过大,就会产生空间混叠,从而误解了沿肌肉传播的复杂电波模式,就像我们误判了旋转车轮的方向一样。理解这个简单的原则是忠实记录人体电气乐谱的第一步。

从原始噪声到神经驱动

我们记录了信号之后,该如何理解它呢?原始的、均值为零的噼啪声并不是肌肉用力的直接度量。我们需要对其进行处理,以提取来自神经系统的潜在指令信息。这类似于将人群混乱的呐喊声转化为代表整体兴奋水平的平滑曲线。这个处理后的信号,作为大脑意图的代表,我们称之为​​神经驱动​​。

转换过程包括两个关键步骤:

  1. ​​整流​​:由于原始肌电信号在零点上下振荡,其简单的平均值毫无意义。为了捕捉其强度,我们必须首先去除负值。标准方法是​​全波整流​​,即取信号的绝对值。这将所有负值部分翻转为正值,确保每一次电波动,无论其极性如何,都对我们测量的总活动量做出贡献。

  2. ​​平滑(低通滤波)​​:整流后的信号仍然是一个非常尖锐、高频的混乱信号。为了看到总体趋势——肌肉用力的缓慢增减——我们必须对其进行平滑处理。我们使用​​低通滤波器​​来完成这项工作。顾名思义,它允许低频内容(总体幅值的缓慢变化)通过,同时阻断高频内容(运动单元动作电位的单个尖峰)。结果是一条平滑的包络线,优美地代表了神经指令随时间变化的强度。这种神经驱动信号的典型带宽非常低,通常低于 10-15 赫兹,这反映了我们神经系统调节自主肌肉力量的速度。

这个“先整流,后平滑”的两步过程是肌电分析的基石,它将一个看似随机的信号转化为一个有意义的神经意图度量。

滤波的艺术与科学

当然,在现实世界中,情况要复杂一些。我们美丽的肌电信号常常被不必要的噪声污染,而滤波行为本身也可能引入自身的问题。选择正确的工具需要理解其中涉及的微妙权衡。

清理信号

在考虑神经驱动之前,我们必须先清理我们的记录。有两个主要元凶会破坏我们的数据:

  • ​​运动伪影​​:如果电极或其下的皮肤移动,会产生大的、缓慢的电压波动,通常低于 20 赫兹。
  • ​​工频干扰​​:我们建筑物中的电线会产生一个恒定的电磁场,在我们的灵敏记录中感应出嗡嗡声,通常在 50 赫兹或 60 赫兹。

我们用滤波器来对付这些噪声。​​带通滤波器​​是我们的第一道防线,旨在只保留我们关心的频率。对于表面肌电图,信号的大部分能量位于大约 20 赫兹到 450 赫兹之间,所以我们设置滤波器让这个频带通过,从而拒绝低频的运动伪影和一些高频噪声。

为了消除持续的工频干扰,我们可能会使用一种非常特定的​​陷波滤波器​​,它被设计成只消除 50 或 60 赫兹周围的一个窄频带。但在这里我们遇到了一个关于时间和频率之间深层联系的绝佳例证。一个在频域中非常尖锐的滤波器(一个窄陷波)必须有一个在时域中非常长且振荡的冲激响应。这意味着我们信号中的任何急剧瞬变都会导致滤波器像被敲响的钟一样“振铃”,从而增加可能扭曲我们想要研究的事件的伪影。一个巧妙的替代方案是​​自适应干扰消除​​,我们使用一个单独的天线直接记录工频噪声,然后简单地从我们的肌电信号中减去一个缩放版本的噪声,从而使底层的神经信号保持不变。

时间的暴政:因果性与延迟

滤波时,时间就是一切。一个关键的区别在于​​因果​​滤波器和​​非因果​​滤波器。顾名思义,因果滤波器只使用当前和过去的信息来产生其输出。它无法预见未来。任何实时运行的系统——比如必须立即响应用户意图的机器人外骨骼——必须使用因果滤波器。

但因果性是有代价的:​​延迟​​。每个因果滤波器,无论多么简单,都会引入一个时间滞后,称为群延迟。对于一个标准的线性相位滤波器,这个延迟是可预测的——对于一个长度为 NNN 的 FIR 滤波器,延迟恰好是 (N−1)/2(N-1)/2(N−1)/2 个样本。在一个以 1000 赫兹采样的外骨骼控制器中,一个看似无害的长度为 64 的滤波器会引入 31.5 毫秒的延迟。再加上后续 20 毫秒平滑窗口的延迟(另外 9.5 毫秒),总延迟超过 40 毫秒。这听起来可能不多,但在一个闭环人机系统中,这可能是无缝辅助与笨拙、不稳定交互之间的区别。

然而,对于离线分析,我们可以“欺骗”时间。当我们将整个记录存储在计算机上时,我们可以使用​​非因果​​或​​零相位滤波器​​。这包括先在正向对数据进行一次滤波,然后在反向再进行一次。两次滤波的相位失真会完美抵消,从而产生零附加时间延迟。这为我们提供了关于信号时序的更准确的图像,但这是一种在实时世界中无法实现的奢侈,因为它需要预知未来 [@problem-id:4193649]。

更深层次的意义

有了一个干净、时序准确的神经驱动信号,我们就可以开始提出更深层次的问题。

什么是“幅值”?

我们将神经驱动定义为肌电信号的“幅值”,但我们应该如何测量它呢?两种最常用的方法是​​平均绝对值 (MAV)​​——这正是我们的“整流-平滑”过程所近似的——和​​均方根 (RMS)​​,后者与信号的功率有关。很长一段时间里,人们认为这两种方法在很大程度上是可以互换的。我们现在知道事实并非如此。

MAV 和 RMS 之间的关系取决于肌电信号幅值分布的统计形状。运动单元的募集方式或放电同步方式的改变可以改变这种形状。例如,一个信号可能从具有钟形的高斯分布变为更“尖峰”的拉普拉斯分布,即使其总功率(以及因此的 RMS 值)保持不变。当这种情况发生时,MAV 会发生变化!一个在一种条件下使用 MAV 校准力估计模型的研究者可能会发现,在另一种条件下,该模型会系统性地低估真实力,这仅仅是因为肌电信号的统计特性发生了变化。这揭示了我们的处理选择不仅仅是技术性的;它们与潜在的生理学紧密交织在一起。

聚焦镜头:空间滤波

回到高密度肌电图 (HD-EMG),电极网格不仅仅是为了制作漂亮的图片;它是一个强大的天线阵列,允许我们在空间上对信号进行滤波。通过简单地改变我们组合相邻电极信号的方式,我们可以极大地改变我们所“看到”的:

  • ​​单极记录​​(每个电极相对于一个远端参考点)为我们提供了原始的电位图谱。它的视野开阔,但容易受到邻近肌肉“串扰”的模糊影响。
  • ​​双极记录​​(两个相邻电极之间的差值)起到空间一阶导数的作用。它能抑制均匀的、远处的信号,并锐化图像,突显电极正下方的活动。
  • ​​拉普拉斯记录​​(从中心电极减去邻近电极的平均值)近似于空间二阶导数。这起到强大的空间高通滤波器的作用,极大地锐化焦点,以至于我们常常可以分离出单个运动单元的活动,有效地使系统对人群的喧嚣“充耳不闻”,而能听到单个的声音。

时频困境

有时我们不仅需要知道肌肉何时活跃,还需要知道该活动的频率内容是什么。经典的工具是​​短时傅里叶变换 (STFT)​​,它通过一个滑动窗口来分析信号。在这里,我们面临著名的时频不确定性原理:短窗口提供出色的时间分辨率但频率分辨率差,而长窗口则相反。这就像拥有一台固定镜头的相机——你要么获得广角视野,要么获得特写视野,但不能两者兼得。

​​连续小波变换 (CWT)​​ 提供了一个优雅的解决方案。它不使用固定窗口,而是使用一系列可以拉伸或压缩的“小波”。在分析高频事件(如肌肉爆发的突然开始)时,它使用短而压缩的小波,提供极佳的时间分辨率。在分析低频现象时,它使用长而伸展的小波,提供出色的频率分辨率。CWT 就像一台带有智能变焦镜头的相机,它会根据所观察的对象自动调整焦距,使其成为一个完美适应肌电信号丰富多变动态的工具。

从信号到科学:一堂谦逊的课

当这些原理帮助我们解决科学难题时,它们的真正威力就显现出来了。考虑一个生物力学中的经典悖论:​​负机电延迟 (EMD)​​的测量。EMD 是指肌电信号与肌肉力量产生之间的生理滞后——这个延迟必须是正的。然而,科学家有时测量到力量似乎在电信号之前就出现了!

这是否违反了因果律?不。它教给我们的是关于系统复杂性和我们测量中假设的知识。我们经常使用逆动力学计算的“力”是关节处的净力矩,而不是来自单个肌肉的力。这个净力矩可以由以下因素产生:

  • 来自其他关节运动的​​交互力矩​​,这些力矩可以在局部肌肉激活之前就开始移动肢体。
  • 从预拉伸的肌腱中​​释放储存的弹性势能​​,这些肌腱像橡皮筋一样工作。
  • ​​关节角度的改变​​改善了肌肉的杠杆作用(力臂),即使在恒定力量下也能增加力矩。

此外,我们应用于肌电信号的处理本身也引入了延迟,将检测到的肌电信号起始点推向更晚的时间。这个悖论并没有揭示物理定律的违背,而是让我们对系统有了更深的理解。它提醒我们,我们处理后的“神经驱动”信号 u(t)u(t)u(t),与“激活度” a(t)a(t)a(t) 的内部生物物理状态不同,而后者本身又与关节产生的最终机械力矩截然不同。肌电信号处理是将身体的电气低语转化为关于运动、意图和生命美妙力学的真实而有意义的故事的、至关重要且错综复杂的艺术。

应用与跨学科联系

现在我们已经探索了肌电图背后的原理——我们如何捕捉和理解我们肌肉的电气低语——我们可以开始一段旅程,看看这些知识将我们带向何方。这是一段多么精彩的旅程!因为肌电图不仅仅是实验室里的一个奇观;它是一个强大的透镜,彻底改变了临床医学、外科学、康复学乃至人类语言研究等多个领域。它是我们探查神经系统的听诊器,让我们能够倾听大脑与身体之间的对话,揭示出卓越协调、细微功能障碍和强大恢复力的故事。

医生探查神经和肌肉的听诊器

想象一下,只通过远距离观察音乐家来诊断一首复杂的管弦乐曲。你可能会看到一个小提琴弓在错误的时间移动,或者一个喇叭手看起来很吃力,但错误的真正性质——一个错误的音符、糟糕的时机、缺乏和谐——仍然难以捉摸。这通常是神经病学面临的挑战。一个病人可能表现出震颤、痉挛或无力,虽然临床观察非常宝贵,但它可能是主观的。肌电图提供了一种客观的语言,一本人体交响乐的乐谱,让我们能够看到底层的神经指令。

思考一下运动障碍这个令人困惑的领域。像震颤、肌阵挛和肌张力障碍等病症,表面上可能看起来相似,但它们源于神经系统中非常不同的功能障碍。通过肌电图,我们可以将电极贴在主动肌-拮抗肌对上,比如屈伸手腕的肌肉,然后倾听。我们听到的声音在每种病症中都惊人地不同。

例如,典型的震颤表现为两条相对肌肉之间优美而有节奏的交替吟唱。一条肌肉收缩,然后另一条,来回往复,其规律性如此精确,以至于频率分析显示出一个尖锐的主峰——在帕金森震颤中通常在 555 赫兹左右。这两个信号是精确的反相关;当一个声音大时,另一个就安静。用信号处理的语言来说,我们说它们具有很高的相干性,但相位差约为 180∘180^{\circ}180∘,就像两个在跷跷板上完美不同步的孩子。

另一方面,肌阵挛不是一种节奏,而是一系列突然、尖锐、不自主的呼喊。肌电图的特征是在主动肌和拮抗肌中同时出现短暂、同步的活动爆发,仅持续几十毫秒。这是一个突然的、全员出动的命令,既短暂又出人意料。

肌张力障碍则呈现出另一种模式:一种持久的、不和谐的争论。在这里,主动肌和拮抗肌长时间同时收缩,相互对抗,处于持续的协同收缩状态。肌电图显示两个通道同时活跃,活动常常“溢出”到本应休息的邻近肌肉中。这是一个系统失去了指令特异性能力的标志,将一个简单的意图性动作变成了一场肌肉的斗争。通过将这些运动转化为定量特征——爆发持续时间、频率、节律性和同步性——肌电图为临床医生提供了一种强大的鉴别诊断工具,一种解读底层乐谱并精确理解音乐是如何出错的方法。

外科医生在微观世界中的向导

肌电图的力量从诊断诊所直接延伸到手术室。想象一位外科医生正在进行精细的喉返神经支配重建手术,这是一个为因神经损伤而失声的患者恢复声音的程序。目标是将一根健康的神经重新连接到控制声带的微小肌肉上。但在颈部错综复杂的解剖结构中,哪根神经纤维该去哪里?一毫米的错误可能意味着将一条本应用于闭合声带的神经连接到一条应该打开声带的肌肉上。

这就是术中神经监护发挥作用的地方。外科医生使用手持式刺激探头向神经纤维传递一个微小、精确的电流脉冲——就像敲门一样。与此同时,通常放置在患者呼吸管上紧邻声带的肌电电极正在倾听回应。如果外科医生刺激了正确的神经分支,相应的肌肉就会收缩,肌电系统将记录到一个清晰的反应,即一个复合肌肉动作电位。来自甲杓肌(一种声带内收肌)的反应发出的“音符”与来自环杓后肌(唯一的声带外展肌)的反应不同。

整个过程是我们讨论过的原理的一个完美应用。刺激脉冲宽度被选择为接近神经的时值(chronaxie)(通常为 100100100 到 200200200 微秒)以获得最大效率。电流保持在低水平(通常低于一毫安)以避免损伤脆弱的神经组织,这是一项由理解电极-组织界面电荷密度限制所支配的安全措施。通过这种方式“绘制”神经分支图,肌电图为外科医生提供了一份实时的功能路线图,以解剖学本身无法提供的确定性引导他们的手。

损伤后重训管弦乐队

一旦神经重新连接,或者在中风或其他损伤之后,工作远未结束。大脑必须重新学习如何控制新获得神经支配的肌肉。这个神经可塑性的过程,即大脑自我重塑的过程,可能缓慢而充满困难。肌电生物反馈提供了一个直接而直观的导师来指导这个过程。

考虑一个面部移植患者康复的非凡案例。第一个目标是简单地恢复力量。通过在目标肌肉上放置肌电电极,比如提升嘴角形成微笑的颧大肌,我们可以给患者实时反馈——一个声音变大或屏幕上的条形图上升——与他们肌肉的活动成正比。这使得康复团队能够设定客观、定量的目标。患者是否激活了肌肉,超过了背景噪声水平?激活是否稳定且可重复?是否存在来自邻近肌肉的过多污染,即*串扰*?基于这些肌电衍生的指标,治疗师可以安全地决定何时让患者从简单的等长收缩进展到更复杂的动态练习,在追求康复的动力与保护脆弱、愈合组织的需要之间取得平衡。

康复期间一个更微妙的挑战是联带运动,即不自主的关联运动的产生。例如,一个试图微笑的患者可能会发现他们的眼睛同时不受控制地眯起 [@problem-id:5170664]。这是因为再生的神经纤维“接错了线”。在这里,肌电生物反馈成为一种教导选择性的工具。我们可以在目标“微笑”肌肉上放置一个电极,在联带运动的“眼睛”肌肉上放置另一个电极。目标是只有当微笑肌肉活跃且眼睛肌肉安静时,才提供积极的奖励音调。通过定义一个“选择性指数”——期望活动与总(期望加不期望)活动的比率——我们可以创建一个复杂的学习环境,鼓励大脑分离正确的通路并修剪不正确的通路。患者实际上是在学习微调交响乐,鼓励小提琴演奏它们的部分,而不会让圆号意外地加入进来。

生物力学家的运动蓝图

在临床之外,肌电图是理解健康人体运动力学不可或缺的工具,从精英运动员到有背痛的办公室职员。它让生物力学家和物理治疗师能够透视皮肤之下,看到作用中的力量。

许多常见的肌肉骨骼问题,如肩痛,并非由单一的灾难性故障引起,而是由于肌肉协同工作的微妙、长期失衡所致。例如,当你举起手臂时,肩胛骨的平稳上旋依赖于你侧面的前锯肌和你背部的下斜方肌之间精细调谐的“力偶”。如果其中一块肌肉无力或激活过迟,其他肌肉,如上斜方-肌,可能会试图代偿,从而变得过度劳累和占主导地位。这种情况,被称为肩胛骨运动障碍,通常肉眼无法察觉。但通过肌电图,情况就变得清晰了。我们可以看到活动不足的原动肌延迟启动和幅值减小,以及协同肌的过度代偿性活动。肌电图提供了揭示设计缺陷的蓝图,指导针对性的锻炼以恢复适当的平衡。

肌电图还可以告诉我们肌肉如何随时间响应挑战,为我们提供了一个观察疲劳的窗口。如果你让某人举起一个中等重量的物体或保持一个姿势一分钟,他们的肌肉会开始疲劳。虽然他们产生的力保持不变,但肌电信号会发生显著变化。随着大脑募集更多的运动单元来分担负荷,总幅值通常会增加。同时,信号的频率成分会发生变化。中值频率,一个衡量信号“音高”的指标,开始下降 [@problem-id:5104351]。这种频谱偏移是肌肉中代谢副产物积累的直接生理后果,这减慢了电脉冲沿肌纤维的传导速度。这是外周疲劳的标志性特征,清楚地表明管弦乐队的演奏者们正在疲劳,即使他们仍在努力维持那个音符。

解码言语的机制

也许肌电图最复杂的应用之一是在言语科学领域。言语产生是我们拥有的最复杂的运动技能之一,依赖于喉部、下颌和舌头中数十块肌肉的闪电般快速的协调。舌头尤其是一个生物工程的奇迹——一个没有骨骼的“肌肉静水压骨骼”,由交织的内在肌和外在肌组成,使其能够呈现出令人难以置信的各种形状。

大脑是如何编排这场肌肉芭蕾来产生 /t/、/s/ 和 /n/ 之间细微差别的呢?这些音都是在舌尖接触同一位置(牙齿后的齿龈脊)时发出的,但它们在声学上是截然不同的。使用细如发丝的线电极进行的肌内肌电图,使我们能够倾听负责这些动作的单个肌肉。对于 /t/,我们看到舌上纵肌短暂而尖锐的活动爆发以抬起舌尖,以及舌横肌的活动以使其变硬,形成一个牢固、不透气的封闭。对于 /s/,同样的肌肉也活跃,但是以一种持续、强直的方式来维持一个精确、狭窄的通道,让空气冲过,产生湍流。而对于 /n/,模式再次相似,但保持时间更长,对应于持续的鼻音。肌电图使我们能够解构这种复杂的行为,揭示构成口语基石的特定肌肉协同作用。

不速之客:当信号即是噪声

在我们的旅程中,我们一直将肌电图视为感兴趣的信号。但在最后一个转折中,当肌电图是噪声时,其原理同样重要。考虑记录脑电图 (EEG) 的挑战,即来自大脑皮层的微弱电信号。这些信号非常微小——几十微伏——而来自头皮和面部肌肉的肌电信号可能高达几百微伏或更多。当一个人咬紧牙关或皱起眉头时,产生的肌电活动可以完全淹没我们试图测量的精细脑信号 [@problem-id:4477166]。

我们如何区分它们?再一次,频谱特性来拯救我们。脑信号的功率谱通常随频率急剧下降。而肌肉信号,正如我们所见,其频谱要平坦得多,富含高频能量。通过计算谱指数——一个衡量功率下降速度的指标——我们可以为信号创建一个指纹。陡峭的斜率表明信号源于皮层;平坦的斜率则指向肌肉污染。正是这个原理使我们能够识别并(通过仔细滤波)去除不必要的肌电伪影,擦净窗户,以便我们能再次看到下方的脑活动。这是一个美丽的科学统一性的证明,即让我们在一个情境中理解信号的相同原理,也赋予我们在另一个情境中消除它的力量。

从手术室到运动实验室,从婴儿的第一声啼哭到句子的复杂细微之处,我们肌肉的电气生命讲述了一个丰富多彩的故事。肌电图赋予我们倾听的能力,提供了与人类运动、健康和表达机制的深刻联系。它远不止是一条弯曲的线;它是一种等待被阅读的语言。