·
肌电信号（EMG）传感器采集的信号包含大量噪声，有效的滤波方法是提取有用信息的关键。以下是肌电信号处理中常用的滤波方法，通常分为几个阶段：
一、 预处理（去除主要干扰）
- 工频陷波滤波器：
- 目的： 消除最主要的干扰 - 来自电网的 50Hz 或 60Hz 工频干扰及其谐波。
- 原理： 在工频频率点（如 50Hz, 60Hz, 100Hz, 150Hz 等）设置一个带阻滤波器，将这部分频率的能量滤除。
- 注意点： 设计陷波滤波器时需注意带宽不宜过宽，以免过多滤除有用信号（特别是基波50/60Hz附近的有用成分）。自适应陷波滤波器可以跟踪漂移的工频干扰。
- 高通滤波：
- 目的： 去除低频噪声（如运动伪迹 - 由导线或传感器运动引起）、电极-皮肤界面产生的缓慢变化的直流偏移（DC Offset）以及部分体动噪声。这些噪声通常在 1-20Hz 以下。
- 参数： 截止频率通常设置在 5Hz 到 20Hz 之间 (常见如 10Hz, 20Hz)。更高的截止频率（如 20Hz）能更有效地去除运动伪迹，但可能会损失部分有用的低频 EMG 能量。
- 类型： 常用一阶或二阶巴特沃斯高通滤波器。
- 直流阻断：
- 目的： 专门去除信号中的恒定直流分量（DC Offset）。
- 实现： 通常是一个截止频率非常低（< 0.5Hz）的高通滤波器，或者在软件中简单减去信号的平均值（静态基线）。
二、 主滤波阶段（去噪与特征提取）
- 带通滤波：
- 目的： 这是核心滤波环节，旨在保留 EMG 信号的主要能量区间（通常 10Hz-500Hz，主能量 50-150Hz），同时抑制带外的噪声（高频噪声如电子设备热噪声、环境射频干扰；预处理阶段未能完全滤除的低频噪声）。
- 参数：
- 低截止频率：通常在 10-30Hz 之间 (例如 20Hz)，与高通预处理相结合或替代（如果高通设置足够）。
- 高截止频率：通常在 250-500Hz 之间 (例如 400Hz, 500Hz)。限制高频噪声。
- 类型：
- 巴特沃斯滤波器： 通带内最平坦，常用，阶数选择（2-4阶常见）影响陡度。
- 切比雪夫滤波器： 通带内允许纹波，但过渡带更陡峭（在相同阶数下）。
- 椭圆滤波器： 通带和阻带都允许纹波，过渡带最陡峭（在相同阶数下），但设计复杂。
- FIR vs IIR：
- FIR滤波器： 本质稳定、线性相位（无相位畸变），但阶数通常较高，计算量大，实时应用有延迟。常用设计方法：窗函数法（Hamming, Blackman等）、等纹波法（Remez交换算法）。
- IIR滤波器： 阶数低，计算效率高，有相位非线性（在时序分析如神经传导速度测量中可能影响结果）。常用设计方法：基于模拟原型（Butterworth, Chebyshev, Bessel, Elliptic）的双线性变换。
- 选择： 如果后续分析对信号的相位信息敏感（如各通道间的相对时序、波形分析），优先选择 FIR。如果对实时性和计算效率要求高且相位误差可容忍，可选 IIR。
- 自适应滤波：
- 目的： 动态去除与参考信号相关的噪声，特别是对漂移的工频干扰特别有效。
- 原理： 使用一个参考输入（通常是纯净的工频正弦波或其测量值）和一个自适应算法（如 LMS, NLMS, RLS）来估计和减去噪声。
- 优势： 能跟踪噪声频率和幅度的变化。
- 小波去噪：
- 目的： 非平稳信号的有效去噪方法。适合去除突发性、幅度变化的噪声和抑制各种白噪声。
- 原理： 将信号分解到不同的“尺度”（频率带）和位置，对小波系数应用阈值（软阈值或硬阈值），然后重构信号。
- 步骤：
- 选择小波基函数（如 Daubechies, Symlets）。
- 选择分解层数。
- 为各层系数选择合适的阈值。
- 阈值处理系数。
- 小波重构。
- 优势： 能在不同频带上有效区分信号和噪声细节。
三、 特殊噪声处理
- 运动伪迹/体动噪声：
- 挑战： 低频分量强，频谱可能与 EMG 低频部分重叠。
- 应对方法：
- 使用相对高通（20Hz 或更高）。
- 加速度计辅助： 用加速度计测量运动信号作为参考输入，通过自适应滤波从 EMG 中减去运动伪迹。效果显著。
- 盲源分离： 如独立成分分析，假设伪迹源与 EMG 源统计独立，可以分离并去除。
- 基线漂移：
- 主要依靠有效的高通滤波或直流阻断。
四、 后处理（平滑与包络提取）
- 移动平均滤波器：
- 目的： 平滑信号，进一步减少高频波动/噪声，计算 EMG 信号包络（代表肌肉激活程度/幅值变化）。常用作全波整流后的低通平滑。
- 原理： 计算信号在一个滑动窗口内的平均值。
- 参数： 窗口长度是关键，通常在 100ms-300ms 之间，取决于信号的速率和应用。本质是一个低通滤波器。
- 包络提取专用低通滤波：
- 目的： 在信号经过整流（绝对值或平方）后，使用低通滤波器提取其平滑的包络线。
- 参数： 截止频率通常很低，用于反映肌肉收缩力随时间变化的速度，常用 2Hz - 5Hz。
- 希尔伯特变换：
- 目的： 直接计算信号的解析信号，其模值即为信号的瞬时包络。
总结与选择建议：
- 标准流程：
- (可选) 直流阻断 -> 高通滤波 (10-20Hz) -> 工频陷波 (50/60Hz + 谐波) -> 带通滤波 (主阶段: 如 20-500Hz, IIR/FIR) -> (可选: 特殊噪声处理) -> 整流 (绝对值或平方) -> 包络低通平滑 (2-5Hz)
- 关键考虑因素：
- 应用目标： 是表面 EMG 力估计？肌电控制（手势识别）？临床神经肌肉诊断？研究肌肉激活时序/模式？目标决定关注的特征（波形细节 vs. 幅度包络）以及对实时性的要求。
- 信号质量与噪声源： 环境电磁干扰强弱？受试者运动量大不大？传感器本身质量？
- 处理平台： 在线/实时系统（嵌入式设备）还是离线分析（PC）？计算能力直接影响算法选择（FIR vs IIR, 是否用计算量大的小波/自适应/ICA）。
- 相位敏感性： 若关注波形特征或通道间精确时序，必须选择 FIR 或有线性相位补偿的 IIR。
- 对信号畸变的容忍度： 如果目标是粗略估计肌肉激活水平，IIR 的非线性相位是可以接受的。
常用组合示例：
- 嵌入式/实时系统 (肌电控制)： 高通 IIR (20Hz) + 自适应工频陷波 + 带通 IIR (20-450Hz) + 全波整流 + 移动平均/LPF (100ms, ~3Hz)。
- 离线分析/研究 (表面肌电，肌肉活动模式)： 带通 FIR (Butterworth, 20-500Hz, 低阶FIR窗函数法) + 工频陷波 + (可选，如需要用小波去除运动伪迹) + 整流 + LPF (150ms, ~3Hz)。
- 临床EMG/NCS (神经传导速度，波形分析)： 强调线性相位，使用高通 FIR (10-20Hz) + 陷波 FIR + 带通 FIR (2-10kHz, 取决于信号类型，需专门设置) + 严格避免引入相位失真
