1. 肌电信号基础：从肌肉到电信号
大家好，我是你们这门课的主讲。在正式开始讲非负矩阵分解之前，我觉得有必要先把肌电信号的基础捋一遍。你想想看，如果连信号是怎么来的都不清楚，后面做协同分析就容易心里没底。我个人习惯是，不管多复杂的算法，先回到物理源头看一看。
1.1 骨骼肌的生理结构
骨骼肌，说白了就是我们身上那些能听大脑指挥的肌肉。它由很多肌纤维捆在一起组成。每条肌纤维就是一个肌细胞，细长细长的，直径也就几十微米。
这里有个关键概念——运动单位。一个运动神经元，加上它支配的所有肌纤维，就构成一个运动单位。我刚开始接触这个领域时，总觉得这概念太简单。后来做实验才发现，运动单位的募集顺序和放电模式，直接决定了肌电信号的形态。嗯，这里要注意，不同运动单位的肌纤维数量差别很大。比如眼外肌一个运动单位可能只有十几条纤维，而大腿的股四头肌能有上千条。
核心要点：运动单位是肌肉收缩的基本功能单元。肌电信号，本质上就是这些运动单位动作电位的叠加。
1.2 动作电位的产生与传导
动作电位是怎么来的？简单说就是细胞膜内外离子分布变化导致的电位波动。静息状态下，膜内电位比膜外低约 -70mV。当收到刺激，钠离子通道打开，大量钠离子涌入，膜电位迅速上升，甚至能冲到 +30mV 左右。这就是去极化。
之后钾离子外流，膜电位回落，进入复极化。整个过程也就几毫秒。但就是这几毫秒，决定了我们能不能拿起一杯水。
动作电位沿着肌纤维传导，速度大概在 3-5 m/s。我在项目中遇到过一个问题：采集到的信号波形总是不对，后来发现是电极贴的位置离神经肌肉接头太远了，信号传导有延迟。所以，电极位置真的很重要。
小技巧：动作电位的传导速度受温度影响。冬天做实验时，肌肉温度低，传导速度会变慢，信号频谱会向低频偏移。记得做实验前让受试者充分热身。
1.3 表面肌电信号的采集原理
表面肌电，就是贴在皮肤表面测到的电信号。电极一般用Ag/AgCl材料，导电性好，极化电位低。我建议新手先别急着用高密度阵列电极，从双极导联开始，踩坑少一些。
采集系统的基本流程是这样的：
- 电极拾取信号——皮肤表面的电位变化
- 前置放大——信号很弱，只有毫伏级，需要放大1000倍左右
- 滤波——滤掉工频干扰和运动伪迹
- 模数转换——变成数字信号，采样率一般设1000-2000 Hz
我曾经踩过一个坑：采样率设得太低，只有500 Hz。结果分析出来的协同结构全是错的。为什么？因为肌电信号的能量主要集中在20-500 Hz，采样率低于1000 Hz就会发生频谱混叠。说白了，就是高频成分伪装成了低频，数据全废了。
避坑指南：我曾经因为贪图方便，用普通心电图电极来采集肌电。结果阻抗不匹配，信号质量一塌糊涂。肌电采集一定要用专用的肌电电极，或者至少保证电极阻抗在5 kΩ以下。
1.4 肌电信号的时频域特征
拿到一段肌电信号，我们怎么描述它？时域和频域是两个基本视角。
时域特征
时域就是看信号随时间怎么变。常用的指标有：
- 均方根值（RMS）——反映信号的能量大小，跟肌肉收缩力正相关
- 积分肌电值（iEMG）——一段时间内信号绝对值的积分
- 过零点数（ZC）——信号穿过零点的次数，跟肌肉疲劳有关
我个人的经验是，RMS最常用，但别只看RMS。有一次我分析抓握动作，RMS显示力量在增加，但实际受试者已经累了。后来看了频谱才发现，中位频率下降了，这才是疲劳的信号。
频域特征
频域分析，就是把信号从时间轴转到频率轴。肌电信号的频谱范围一般在20-500 Hz，能量主要集中在50-150 Hz。
常用指标：
| 指标 | 含义 | 应用场景 |
|---|---|---|
| 平均功率频率（MPF） | 频谱的加权平均 | 肌肉疲劳检测 |
| 中位频率（MF） | 将频谱分成两半的频率点 | 疲劳分析，比MPF更稳定 |
| 功率谱密度（PSD） | 各频率成分的能量分布 | 协同模式识别 |
你想想看，为什么疲劳时中位频率会下降？因为疲劳后动作电位传导速度变慢，信号中的高频成分衰减了。这个现象我在做康复评估项目时反复验证过，很可靠。
1.5 本章知识体系
下面这张图，是我自己梳理的本章知识结构。你可以把它当作一个地图，后面讲到非负矩阵分解时，随时回来对照。
好了，这一章的内容就到这里。肌电信号的基础，说白了就是搞清楚信号从哪里来、怎么采、怎么看。后面讲非负矩阵分解时，你会发现这些基础概念会反复出现。尤其是运动单位的募集和频谱特征，直接决定了协同结构的物理意义。
一句话总结：肌电信号是运动单位动作电位的时空叠加。理解了这个本质，后面的协同分析就有了根基。