Investigation on parametric analysis of dynamic EMG signals by amuscle-structured simulation model

In the analysis of electromyographic (EMG) signals during dynamic movement, we have proposed an estimation algorithm for the time-varying parameters of an autoregressive model. The parameters correspond to less biased time-varying reflection coefficients. We determined the less biased estimation using a locally quasi-stationary model and named these parameters "k parameters." We estimated k parameters up to the fifth order for the surface EMG signals of a masseter muscle during rapid open-close movement of the lower jaw, a ballistic contraction, and fatigue. According to the results, the time courses of the k parameters displayed remarkable properties. In order to study the behavior of k parameters physiologically, we produced a muscle-structured simulation model based on anatomical and physiological data. The simulation results suggested that the behavior of the third parameter is related to the number of active motor units (MU's) at the shallow layer of a muscle. The detailed recruitment mechanism in terms of the MU's types has not yet been solved. Although further study is required, the parametric analysis using k parameters offers a new perspective for evaluation of muscle dynamics during several movements.     

Modeling of surface myoelectric signals. II. Model-based signalinterpretation

Experimental electromyogram (EMG) data from the human biceps brachii were simulated using the model described in [10] of this work. A multichannel linear electrode array, spanning the length of the biceps, was used to detect monopolar and bipolar signals, from which double differential signals were computed, during either voluntary or electrically elicited isometric contractions. For relatively low-level voluntary contractions (10%-30% of maximum force) individual firings of three to four-different motor units were identified and their waveforms were closely approximated by the model. Motor unit parameters such as depth, size, fiber orientation and length, location of innervation and tendonous zones, propagation velocity, and source width were estimated using the model. Two applications of the model are described. The first analyzes the effects of electrode rotation with respect to the muscle fiber direction and shows the possibility of conduction velocity (CV) over- and under-estimation. The second focuses on the myoelectric manifestations of fatigue during a sustained electrically elicited contraction and the interrelationship between muscle fiber CV, spectral and amplitude variables, and the length of the depolarization zone. It is concluded that a) surface EMG detection using an electrode array, when combined with a model of signal propagation, provides a useful method for understanding the physiological and anatomical determinants of EMG waveform characteristics and b) the model provides a way for the interpretation of fatigue plots.     

Estimation of muscle fiber conduction velocity with a spectral multidip approach

We propose a novel method for estimation of muscle fiber conduction velocity from surface electromyographic (EMG) signals. The method is based on the regression analysis between spatial and temporal frequencies of multiple dips introduced in the EMG power spectrum through the application of a set of spatial filters. This approach leads to a closed analytical expression of conduction velocity as a function of the auto- and cross-spectra of monopolar signals detected along the direction of muscle fibers. The performance of the algorithm was compared with respect to that of the classic single dip approach on simulated and experimental EMG signals. The standard deviation of conduction velocity estimates from simulated single motor unit action potentials was reduced from 1.51 m/s [10 dB signal-to-noise ratio (SNR)] and 1.06 m/s (20 dB SNR) with the single dip approach to 0.51 m/s (10 dB) and 0.23 m/s (20 dB) with the proposed method using 65 dips. When 200 active motor units were simulated in an interference EMG signal, standard deviation of conduction velocity decreased from 0.95 m/s (10 dB SNR) and 0.60 m/s (20 dB SNR) with a single dip to 0.21 m/s (10 dB) and 0.11 m/s (20 dB) with 65 dips. In experimental signals detected from the abductor pollicis brevis muscle, standard deviation of estimation decreased from (mean +/- SD over 5 subjects) 1.25 +/- 0.62 m/s with one dip to 0.10 +/- 0.03 m/s with 100 dips. The proposed method does not imply limitation in resolution of the estimated conduction velocity and does not require any iterative procedure for the estimate since it is based on a closed analytical formulation.     

Assessment of average muscle fiber conduction velocity from surface EMG signals during fatiguing dynamic contractions

In this paper, we propose techniques of surface electromyographic (EMG) signal detection and processing for the assessment of muscle fiber conduction velocity (CV) during dynamic contractions involving fast movements. The main objectives of the study are: 1) to present multielectrode EMG detection systems specifically designed for dynamic conditions (in particular, for CV estimation); 2) to propose a novel multichannel CV estimation method for application to short EMG signal bursts; and 3) to validate on experimental signals different choices of the processing parameters. Linear adhesive arrays of electrodes are presented for multichannel surface EMG detection during movement. A new multichannel CV estimation algorithm is proposed. The algorithm provides maximum likelihood estimation of CV from a set of surface EMG signals with a window limiting the time interval in which the mean square error (mse) between aligned signals is minimized. The minimization of the windowed mse function is performed in the frequency domain, without limitation in time resolution and with an iterative computationally efficient procedure. The method proposed is applied to signals detected from the vastus laterialis and vastus medialis muscles during cycling at 60 cycles/min. Ten subjects were investigated during a 4-min cycling task. The method provided reliable assessment of muscle fatigue for these subjects during dynamic contractions.     

Autoregressive Modeling of Surface EMG and Its Spectrum with Application to Fatigue

The following is an investigation of the ability of the autoregressive (AR) model to describe the spectrum of the processes underlying the recorded surface EMG. Surface EMG (SEMG) spectrum is influenced by two major factors; one attributed to the motor units (MU) firing rate and the second, the higher frequency one, to the morphology of the action potentials (AP) traveling along the muscle fiber. In the present paper, SEMG measurements were carried out on the biceps brachii muscle with fixed surface electrodes arrangement and isotonic conditions. Sufficient averaging of 0.5 s segments enabled the identification of the low-frequency peak related to the firing rates of the MU's.     

Direct On-Line Estimation of Muscle Fiber Conduction Velocity by Surface Electromyography

A new approach to the measurement of muscle fiber conduction velocity by surface electromyography (EMG) is developed, based upon a triple surface-electrode configuration. Unlike previously reported methods, it does not require the Fourier transformation of raw EMG data, and can be implemented on-line on a simple digital computer. Although programmed here on a PDP-11 machine, the algorithms developed are equally suited to a 16-bit microprocessor, without any need for hardware multiplication or floating-point facilities. Accurate conduction velocity estimates may be obtained every few seconds.     

Fatigue compensation of the electromyographic signal for prostheticcontrol and force estimation

During a sustained muscle contraction, the amplitude of electromyographic (EMG) signals increases and the spectrum of the EMG signal shifts toward lower frequencies. These effects are due to muscular fatigue and can cause problems in the control of myoelectric prostheses and in the estimation of contraction level from the EMG signal. It has been well known that the fatigue effects can be explained by the conduction velocity changes during the fatigue process and by the idea that the conduction velocity is linearly proportional to the median frequency of EMG signals. Hence the fatigue process can be monitored by measuring the median frequency. A fatigue compensation preprocessor has been developed. It uses the widely accepted power spectrum density model of EMG signals that contains the conduction velocity as a measure of fatigue. It was verified that the preprocessor scales down the amplitude of the fatigued EMG signal and decompresses the spectrum. Hence, the preprocessor eliminates the increase in amplitude and the shift in frequency and enables consistent EMG signals to be used to control prostheses     

Physiologically based simulation of clinical EMG signals

An algorithm that generates electromyographic (EMG) signals consistent with those acquired in a clinical setting is described. Signals are generated using a model constructed to closely resemble the physiology and morphology of skeletal muscle, combined with line source models of commonly used needle electrodes positioned in a way consistent with clinical studies. The validity of the simulation routines is demonstrated by comparing values of statistics calculated from simulated signals with those from clinical EMG studies of normal subjects. The simulated EMG signals may be used to explore the relationships between muscle structure and activation and clinically acquired EMG signals. The effects of motor unit (MU) morphology, activation, and neuromuscular junction activity on acquired signals can be analyzed at the fiber, MU and muscle level. Relationships between quantitative features of EMG signals and muscle structure and activation are discussed.     

A Microprocessor-Based Polarity Correlator for the On-Line Measurement of Muscle Fiber Conduction Velocity

The measurement of the average conduction velocity of action potentials along muscle fibers has important applications in muscle fatigue study and ergonomics. A cross-correlation-based method for the estimation of this velocity from surface electromyography (EMG) is presented. It is implemented on a standard low-cost 8-bit microprocessor. The use of cross correlation for such an on-line application has been made possible by the use of an efficient polarity correlation algorithm.     

Simulation of surface EMG signals for a multilayer volume conductor with triangular model of the muscle tissue

This study analytically describes surface electromyogram (sEMG) signals generated by a model of a triangular muscle, i.e., a muscle with fibers arranged in a fan shape. Examples of triangular muscles in the human body are the deltoid, the pectoralis major, the trapezius, the adductor pollicis. A model of triangular muscle is proposed. It is a sector of a cylindrical volume conductor (with the fibers directed along the radial coordinate) bounded at the muscle/fat interface. The muscle conductivity tensor reflects the fan anisotropy. Edge effects have been neglected. A solution of the nonspace invariant problem for a triangular muscle is provided in the Fourier domain. An approximate analytical solution for a two plane layer volume conductor model is obtained by introducing a homogeneous layer (modeling the fat) over the triangular muscle. The results are implemented in a complete sEMG generation model (including the finite length of the fibers), simulating single fiber action potentials. The model is not space invariant due to the changes of the volume conductor along the direction of action potential propagation. Thus the detected potentials at the skin surface change shape as they propagate. This determines problems in the extraction and interpretation of parameters. As a representative example of application of the simulation model, the influence of the inhomogeneity of the volume conductor in conduction velocity (CV) estimation is addressed (for two channels; maximum likelihood and reference point methods). Different fiber depths, electrode placements and small misalignments of the detection system with respect to the fiber have been simulated. The error in CV estimation is large when the depth of the fiber increases, when the detection system is not aligned with the fiber and close to the innervation point and to the tendons.     

Nerve-bundle conduction velocity distribution measurement and transfer function analysis

Human peripheral nerves are composed of thousands of individual nerve fibers whose signal conduction velocities vary from 0.2 to 100 m/s. Though good correlation exist between fiber velocity and the physiological function subserved by these fibers, considerable variation is found for specific end organs (20 to 40 m/s for a single muscle). Though clinical neurophysiologists have routinely measured the maximum conduction velocity of peripheral motor nerves for 30 years, it has not been possible to determine the velocity distributions. This report details a model and noninvasive methodology for motor axon conduction latency distribution mesurement. This distribution, proportional to velocity dispersion, possesses a low-pass filter characteristic which agrees with theoretical predictions and may be significant in some sensory and motor functions.     

The sensory region of muscle spindles in the posterior cricoarytenoid muscle of the marmoset

We examined muscle spindles in the posterior cricoarytenoid (PCA) muscle of the marmoset. The spindle capsule contained only one intrafusal muscle fiber which did not form a typical aggregation of equatorial nuclei such as bags and chains. The intrafusal fiber was innervated by sensory endings which branched and arranged irregularly. These sensory endings lay on the surface of the fiber in shallow grooves or deeply penetrated into it. These findings suggest that the muscle spindles in the marmoset PCA muscle are strikingly different in the structure of intrafusal fibers in the arrangement of sensory endings from those in common mammalian skeletal muscles.     

A novel approach for estimating muscle fiber conduction velocity by spatial and temporal filtering of surface EMG signals

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.     

Estimation of Myoelectrc Conduction Yelocity Distribution

An estimator of muscle fiber conduction velocity probability density function is derived. It is based on the measurement of the cross-and autopower spectra of surface myoelectric signals. Sample experimental results are compared to the results of other methods.     

Simulation of surface EMG signals for a multilayer volume conductor with a superficial bone or blood vessel

This study analytically describes surface electromyogram (EMG) signals generated by a planar multilayer volume conductor constituted by different subdomains modeling muscle, bone (or blood vessel), fat, and skin tissues. The bone is cylindrical in shape, with a semicircular section. The flat portion of the boundary of the bone subdomain is interfaced with the fat layer tissue, the remaining part of the boundary is in contact with the muscle layer. The volume conductor is a model of physiological tissues in which the bone is superficial, as in the case of the tibia bone, backbone, and bones of the forearm. The muscle fibers are considered parallel to the axes of the bone, so that the model is space invariant in the direction of propagation of the action potential. The proposed model, being analytical, allows faster simulations of surface EMG with respect to previously developed models including bone or blood vessels based on the finite-element method. Surface EMG signals are studied by simulating a library of single-fiber action potentials (SFAP) of fibers in different locations within the muscle domain, simulating the generation, propagation, and extinction of the action potential. The decay of the amplitude of the SFAPs in the direction transversal to the fibers is assessed. The decay in the direction of the bone has a lower rate with respect to the opposite direction. Similar results are obtained by simulating motor unit action potentials (MUAPs) constituted by 100 fibers with territory 5 mm2. M waves and interference EMG signals are also simulated based on the library of SFAPs. Again, the decay of the amplitude of the simulated interference EMG signals is lower approaching the bone with respect to going farther from it. The findings of this study indicate the effect of a superficial bone in enhancing the EMG signals in the transversal direction with respect to the fibers of the considered muscle. This increases the effect of crosstalk. The same mathematical method used to simulate a superficial bone can be applied to simulate other physiological tissues. For example, superficial blood vessels (e.g., basilic vein, brachial artery) can influence the recorded EMG signals. As the electrical conductivity of blood is high (it is of the same order as the longitudinal conductivity in the muscle), the effect on EMG signals is opposite compared to the effect of a superficial bone.     

Polarity Correlator for Conduction Velocity Measurement

Techniques for the estimation of skeletal muscle fiber conduction velocity are of considerable interest. These techniques use, in general, some form of cross correlation or zero-crossing analysis. Cross correlation is a straightforward method of conduction velocity estimation, however, it is difficult to realize low-cost real-time processors. Polarity correlation offers a method which preserves the advantages of cross correlation and satisfies the speed-cost constraint. This paper describes the algorithm for implementation of a polarity cofrelator instrument using a microprocessor. The instrument is tested with deterministic and stochastic signals and used to estimate the conduction velocity of biceps brachii fibers of ten normal subjects. The performance of the instrument is compared to the zero-crossing technique described by Lynn.     

Selectivity of spatial filters for surface EMG detection from the tibialis anterior muscle

Many spatial filters have been proposed for surface electromyographic (EMG) signal detection. Although theoretical and modeling predictions on spatial selectivity are available, there are no extensive experimental validations of these techniques based on single motor unit (MU) activity detection. The aim of this study was to compare spatial selectivity of one- and two-dimensional (1-D and 2-D) spatial filters for EMG signal detection. Intramuscular and surface EMG signals were recorded from the tibialis anterior muscle of ten subjects. The simultaneous use of intramuscular wire and surface recordings (with the spike triggered averaging technique) allowed investigation of the activity of single MUs at the skin surface. The surface EMG signals were recorded with a grid of point electrodes (3 x 3 electrodes) and a ring electrode system at 15 locations over the muscle, with the wires detecting signals from the same intramuscular location. For most subjects, it was possible to classify, from the intramuscular recordings, the activity of the same MUs for all the contractions. The surface EMG signals were averaged with the intramuscularly detected MU action potentials as triggers. In this way, eight spatial filters--longitudinal and transversal, single and double differential (LSD, TSD, LDD, TDD), Laplacian (NDD), inverse binomial filter of the second order (IB2), inverse rectangle filter (IR), and differential ring system (C1)--could be compared on the basis of their spatial selectivity. The distance from the source (transversal with respect to the muscle fiber orientation) after which the surface detected potential did not exceed +/- 5% of the maximal peak-to-peak amplitude (detection distance) was statistically smaller for the 2-D systems and TDD than for the other filters. The MU action potential duration was significantly shorter with LDD and with the 2-D systems than with the other filters. The 2-D filters investigated (including C1) showed very similar performance and were, thus, considered equivalent from the point of view of spatial selectivity.     

Influence of tissue inhomogeneities on noninvasive muscle fiberconduction velocity measurements-investigated by physical and numericalmodeling

The determination of conduction velocity in the muscle fibers of single motor units from noninvasive recordings of single motor unit action potentials can be improved by the method of spatially filtering multielectrode EMG. The use of this conduction velocity as a diagnostic tool requires a high reliability of the detected values. However, experiments did reveal that the measured conduction velocity values showed remarkably high fluctuations depending on the recording site along the muscle fibers which could not be attributed to the influence of the endplate and tendon region. The present work examines the hypothesis that the observed fluctuations in propagation velocity were caused by electrically inhomogeneous tissue, regions of different electrical conductivity which are located between the excited muscle fibers and the recording electrodes and which cause a deformation of the extracellular electric current field. The investigation was performed by means of a physical model as well as by finite element model calculations. In both models single, simple shaped (cylindrical) inhomogeneity regions with a conductivity of 0.1 to 10 times that of the surrounding medium and diameters ranging between 1.6 and 2.7 mm were placed between excitation sources and recording site. The results indicate that the observed conduction velocity fluctuations of up to some 10% can be well attributed to inhomogeneity effects of the tissue conductivity. Based on these results, one may look for signal processing methods to cut down such fluctuations in conduction velocity measurements.     

Noninvasive estimation of motor unit conduction velocity distribution using linear electrode arrays

Determining the conduction velocity of motor unit action potentials is one of the most important problems in surface electromyography. The estimate of one average conduction velocity value depends on a variety of uncontrollable factors. More meaningful information is obtained from the estimation of the distribution of the different delays in the myoelectric signals. A solution to the problem is the separation and characterization of the individual components propagating at different velocities. A technique, based on surface electrode array recording, is proposed to estimate motor unit conduction velocity distribution. The method consists in the identification of the single action potentials in the time scale domain (with the continuous wavelet transform) and in the estimation of their conduction velocities based on the beamforming algorithm. The performances of the technique have been evaluated using simulated and real myoelectric signals. The results demonstrate that the technique is accurate and reliable. The method may be useful for the diagnosis of neuromuscular disorders, for the monitoring of muscle fatigue and for noninvasive investigation of individual motor units.     

Simulation of surface EMG signals generated by muscle tissues with inhomogeneity due to fiber pinnation

Surface electromyographic (EMG) signal modeling has important applications in the interpretation of experimental EMG data. Most models of surface EMG generation considered volume conductors homogeneous in the direction of propagation of the action potentials. However, this may not be the case in practice due to local tissue inhomogeneities or to the fact that there may be groups of muscle fibers with different orientations. This study addresses the issue of analytically describing surface EMG signals generated by bi-pinnate muscles, i.e., muscles which have two groups of fibers with two orientations. The approach will also be adapted to the case of a muscle with fibers inclined in the depth direction. Such muscle anatomies are inhomogeneous in the direction of propagation of the action potentials with the consequence that the system can not be described as space invariant in the direction of source propagation. In these conditions, the potentials detected at the skin surface do not travel without shape changes. This determines numerical issues in the implementation of the model which are addressed in this work. The study provides the solution of the nonhomogenous, anisotropic problem, proposes an implementation of the results in complete surface EMG generation models (including finite-length fibers), and shows representative results of the application of the models proposed.