Farina Dario, Merletti Roberto
Dipartimento di Elettronica, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy.
IEEE Trans Biomed Eng. 2003 Dec;50(12):1340-51. doi: 10.1109/TBME.2003.819847.
We describe a new method for the estimation of muscle fiber conduction velocity (CV) from surface electromyography (EMG) signals. The method is based on the detection of two surface EMG signals with different spatial filters and on the compensation of the spatial filtering operations by two temporal filters (with CV as unknown parameter) applied to the signals. The transfer functions of the two spatial filters may have different magnitudes and phases, thus the detected signals have not necessarily the same shape. The two signals are first spatially and then temporally filtered and are ideally equal when the CV value selected as a parameter in the temporal filters corresponds to the velocity of propagation of the detected action potentials. This approach is the generalization of the classic spectral matching technique. A theoretical derivation of the method is provided together with its fast implementation by an iterative method based on the Newton's method. Moreover, the lowest CV estimate among those obtained by a number of filter pairs is selected to reduce the CV bias due to nonpropagating signal components. Simulation results indicate that the method described is less sensitive than the classic spectral matching approach to the presence of nonpropagating signals and that the two methods have similar standard deviation of estimation in the presence of additive, white, Gaussian noise. Finally, experimental signals have been collected from the biceps brachii muscle of ten healthy male subjects with an adhesive linear array of eight electrodes. The CV estimates depended on the electrode location with positive bias for the estimates from electrodes close to the innervation or tendon regions, as expected. The proposed method led to significantly lower bias than the spectral matching method in the experimental conditions, confirming the simulation results.
我们描述了一种从表面肌电图(EMG)信号估计肌肉纤维传导速度(CV)的新方法。该方法基于用不同空间滤波器检测两个表面EMG信号,并通过应用于信号的两个时间滤波器(以CV作为未知参数)对空间滤波操作进行补偿。两个空间滤波器的传递函数可能具有不同的幅度和相位,因此检测到的信号形状不一定相同。这两个信号首先进行空间滤波,然后进行时间滤波,当在时间滤波器中作为参数选择的CV值对应于检测到的动作电位的传播速度时,它们在理想情况下是相等的。这种方法是经典频谱匹配技术的推广。本文给出了该方法的理论推导,以及基于牛顿法的迭代方法的快速实现。此外,选择多个滤波器对得到的CV估计中的最低值,以减少由于非传播信号成分导致的CV偏差。仿真结果表明,所描述的方法对非传播信号的存在比经典频谱匹配方法更不敏感,并且在存在加性白高斯噪声的情况下,这两种方法具有相似的估计标准差。最后,用一个由八个电极组成的粘性线性阵列从十名健康男性受试者的肱二头肌收集了实验信号。正如预期的那样,CV估计值取决于电极位置,靠近神经支配或肌腱区域的电极估计值存在正偏差。在实验条件下,所提出的方法导致的偏差明显低于频谱匹配方法,证实了仿真结果。