# Mechanomyographic Analysis for Muscle Activity Assessment during a Load-Lifting Task

Matthieu Correa, Maxime Projetti, Isabelle A. Siegler, Nicolas Vignais

► To cite this version:

Matthieu Correa, Maxime Projetti, Isabelle A. Siegler, Nicolas Vignais. Mechanomyographic Analysis for Muscle Activity Assessment during a Load-Lifting Task. Sensors, 2023, 23 (18), pp.7969. 10.3390/s23187969. [hal-04448990](https://hal.science/hal-04448990v1)

HAL Id: hal-04448990

https://hal.science/hal-04448990v1

Submitted on 9 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Mechanomyographic Analysis for Muscle Activity Assessment during a Load-Lifting Task

Matthieu Correa<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,* ID, Maxime Projetti<fnref n="3" />, Isabelle A. Siegler<fnref n="1" />,<fnref n="2" /> ID and Nicolas Vignais<fnref n="1" />,<fnref n="2" /> ID

**Abstract:** The purpose of this study was to compare electromyographic (EMG) with mechanomyographic (MMG) recordings during isometric conditions, and during a simulated load-lifting task. Twenty-two males (age: 25.5 ± 5.3 years) first performed maximal voluntary contractions (MVC) and submaximal isometric contractions of upper limb muscles at 25%, 50% and 75% MVC. Participants then executed repetitions of a functional activity simulating a load-lifting task above shoulder level, at 25%, 50% and 75% of their maximum activity (based on MVC). The low-frequency part of the accelerometer signal (<5 Hz) was used to segment the six phases of the motion. EMG and MMG were both recorded during the entire experimental procedure. Root mean square (RMS) and mean power frequency (MPF) were selected as signal extraction features. During isometric contractions, EMG and MMG exhibited similar repeatability scores. They also shared similar RMS vs. force relationship, with RMS increasing to 75% MVC and plateauing to 100%. MPF decreased with increasing force to 75% MVC. In dynamic condition, RMS<sub>MMG</sub> exhibited higher sensitivity to changes in load than RMS<sub>EMG</sub>. These results confirm the feasibility of MMG measurements to be used during functional activities outside the laboratory. It opens new perspectives for future applications in sports science, ergonomics and human-machine interface conception.

**Keywords:** mechanomyography; electromyography; isometric contractions; dynamic contractions; load-lifting

## 1. Introduction

Today, measuring muscle activity in ecological situations is a major issue. It could allow a better understanding of the neuromuscular system and its adaptations to the environment compared to isolated contractions in the laboratory. It may also have significant applications. In sports science, assessing muscle load may enhance performance optimization strategies, injury prevention monitoring and rehabilitation. In physical ergonomics, muscle evaluation can help reduce work-related musculoskeletal disorders through workplace assessment and reconception [1]. In human-machine interface conception, recording muscle activity in real time can be useful to control assisting device such as active exoskeletons, cobots and myoelectronic prostheses.

To effectively measure muscle activity in functional and ecological situations, the method used should be easy to set up and to use, and robust in all environments and during prolonged duration acquisitions. Today, the reference method for muscle activity measurement is electromyography (EMG), which measures the electrical activity of muscles. However, it requires skin preparation and a precise sensor placement for reproducible results [2,3], making it hard to use and to set up outside of a laboratory. Moreover, EMG is

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Correa, M.; Projetti, M.; Siegler, I.A.; Vignais, N. Mechanomyographic Analysis for Muscle Activity Assessment during a Load-Lifting Task. *Sensors* **2023**, *23*, 7969. https://doi.org/10.3390/s23187969

Academic Editor: Ki H. Chon

Received: 31 July 2023

Revised: 15 September 2023

Accepted: 15 September 2023

Published: 19 September 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Laboratoire CIAMS (Complexité, Innovation, Activités Motrices et Sportives), Université Paris-Saclay, CEDEX, 91405 Orsay, France; isabelle.siegler@universite-paris-saclay.fr (I.A.S.); nicolas.vignais@universite-paris-saclay.fr (N.V.)</footnote>

<footnote n="2">² Laboratoire CIAMS (Complexité, Innovation, Activités Motrices et Sportives), Université d'Orléans, 45067 Orléans, France</footnote>

<footnote n="3">³ Moten Technologies, 92800 Puteaux, France</footnote>

<footnote>* Correspondence: matthieu.correa@universite-paris-saclay.fr</footnote>
sensitive to interference and changes in impedance, making it non-robust to the environ-
ment and during extended time acquisitions, such as when sweating occurs, for example [4].
Finally, EMG is sensitive to motion artifacts because of the electrodes' shift relative to mus-
cle fibers in dynamic muscle actions [5,6]. Altogether, these limitations make EMG hardly
suitable for recording muscle activity in functional and ecological conditions [3].

Other techniques have emerged for recording muscle activity such as mechanomyo-
graphy (MMG), the mechanical counterpart of EMG, which is the measurement of the
low-frequency lateral oscillations of active muscle fibers [7]. MMG reflects three physiolog-
ical phenomena: (i) the gross lateral movement of the muscle at the start of the contraction,
(ii) subsequent vibrations at the muscle's resonance frequency and (iii) dimensional changes
of active muscle fibers [8]. MMG can be obtained by multiple means: piezoelectric crystals,
condenser microphones, displacement lasers and accelerometers [9]. Yet, accelerometers
are often preferred because of their flat frequency range, ease of fixation on the skin, inex-
pensiveness, lightweightness and their measurement in physical units (m/s²). However,
accelerometers are highly sensitive to motion artifacts [10]. Thus, without proper filtering
or experimental design, the MMG signal could be polluted with motion artifacts, biasing
its interpretation. Most studies used a band-pass filter between 5 and 100 Hz to eliminate
both motion artifacts and high-frequency noise [8]. Recently, advances in MMG signal
processing methods have led to more efficient motion artifact filtering in dynamic condi-
tions [11,12]. Being sensitive to motion could also be viewed as a strong argument in favor
of accelerometer-based MMG, as the low-frequency part of the signal (motion accelerations)
can be used as a segmenting tool [13], making the most out of the overall signal.

MMG appears to have a number of advantages over EMG such as no skin preparation,
ease of fixation and less sensitivity to sensor placement [14], making it more convenient.
Moreover, it is more robust than EMG thanks to its non-sensitivity to changes in impedance
and interference, thus guaranteeing signal quality over long periods of time. They can
also be used in a "2 in 1" fashion, for both measuring muscle activity and segmenting
postures and motions. Altogether, these technical aspects make MMG a convenient and
suitable alternative to EMG for measuring and segmenting functional activities in controlled
environments and ecological situations.

Despite a great number of studies regarding isometric and isokinetic contractions, very
little MMG literature characterizes MMG under poly-articular functional activities such
as squatting [3,13]. To the best of our knowledge, no MMG studies were found regarding
load-lifting tasks, despite being omnipresent in sports, manual work and daily activities.
First of all, such research could assess the validity of MMG measurements under these
conditions. Moreover, it could provide significant insights about muscle synergies and
neuromuscular adaptations in functional situations.

Thus, the objectives of this study are twofold: (i) comparing a new accelerometer-based
MMG system to an EMG reference system during isometric contractions and (ii) analyzing
a load-lifting task with both MMG and EMG to explore MMG applications during complex
dynamic conditions.

**2. Materials and Methods**

2.1. *Participants*

Only male subjects were selected in this study for gender considerations [15]. Par-
ticipants should not have suffered from any musculoskeletal disorders in the 6 months
prior to the study. Given a significance criterion of α = 0.05 and an a priori large effect size,
the optimal number of participants for the current study was 18, with an actual statistical
power of 0.96.

Twenty-two active males (age: 25.5 ± 5.3 years; height: 1.79 ± 0.07 m; weight: 75.7 ± 11.9 kg) participated in this study. All subjects turned out to be right-handed. Participants provided written informed consent before taking part in this experiment. This study was approved by the Academic Research Ethics Committee (Université Paris-Saclay, 2021-287).
## 2.2. Sensors

EMG signal was recorded using sensors (Miniwave, Cometa, Montegranaro, Italy) sampled at 2000 Hz and positioned following the SENIAM recommendations [2]. Prior to EMG sensor placement, the skin was shaved and sanitized with 70% alcohol-soaked wipes. MMG signal was recorded using accelerometer sensors (Moten Technologies, Puteaux, France) with a 1000 Hz sampling rate. The MMG sensor specifications are displayed in Table 1. Each sensor was placed on the muscle's belly following recommendations from the literature [16]. A force sensor was used to measure external force at wrist level with a sampling rate of 2000 Hz (Sensor 2715-ISO, Sensy, Jumet, Belgium). Sensors were placed only on the dominant arm of each subject (see Figure 1 for EMG and MMG sensors' location).

<table><caption>**Table 1.** MMG sensor technical specifications.</caption><tr><td>Type</td><td>3-axis capacitive accelerometer with digital output</td></tr><tr><td>Mass</td><td>13 g</td></tr><tr><td>Sensor resonant frequency</td><td>2.4 kHz</td></tr><tr><td>Output full scale range (FSR)</td><td>±2 g</td></tr><tr><td>Non linearity</td><td>0.1% FSR</td></tr><tr><td>Scale factor</td><td>3.9 µg/LSB (less significant bit)</td></tr><tr><td>Dynamic range</td><td>98 dB</td></tr></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** EMG and MMG sensors' location for the long portion of the Biceps Brachii (A), the lateral Deltoid (B) and the long head of the Triceps Brachii (C).</figcaption></figure>

## 2.3. Isometric Testing Procedure

The isometric testing procedure comprised two parts: (i) the measurement of the maximum voluntary contraction (MVC) of the Biceps Brachii (BB), Triceps Brachii (TRI) and the lateral Deltoid (DEL) muscles, and (ii) the submaximal contractions at 25%, 50% and 75% of MVC for each muscle.

For the MVC measurements, participants had to perform three maximal contractions of 3 s onto the force sensor attached to an adjustable metallic structure. Trials were separated by two-minute rest periods. For the BB and the TRI muscle, the participant had the elbow flexed at 90° against the body and the force sensor was placed on the wrist. The force sensor was facing downward for the BB and upward for the TRI. For the DEL muscle, the arm was extended forward parallel to the ground with the force sensor on top of the wrist facing downward. For each muscle, the MVC value was defined as the highest value obtained across the three trials [17].
The participant then had to maintain submaximal contractions at 25, 50 and 75% of MVC for 15 s, thanks to a visual feedback of force level displayed on a screen. Each submaximal contraction was repeated three times, with a 1.5-min rest period.

## 2.4. *Dynamic Testing Procedure*

The dynamic testing consisted of 10 repetitions of a load lifted above shoulder level with both hands. The different weights of the load were derived from the MVC obtained for the BB muscle in the isometric condition. More precisely, the force obtained during MVC was converted into a mass to lift (25%, 50% and 75% MVC). Indeed, BB was chosen over TRI and DEL muscles as a reference because it was the strongest muscle of the three and thus allowed to have the largest load differences between intensities. Each repetition started with lifting the charge from the mid-tibia level to above shoulder level and getting back into the initial position (see Figure 2). The participant was asked to maintain a constant speed throughout every repetition. Repetition rate was set at six seconds using a metronome. Prior to the experiment, participants had a training session with an empty bar (2 kg) to practice the movement, technique and required rate.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Example of the load-lifting task performed during the dynamic testing procedure. One repetition consisted of going from posture (a–d) and going back to initial posture. Postures were separated by one second from each other using a metronome. Participants were asked to keep a constant speed throughout every repetition.</figcaption></figure>

## 2.5. *Signal Processing*

EMG and force signals were collected through the QTM software version 2022 (Qualisys, Göteborg, Sweden) allowing their synchronization. MMG signal was obtained using another software (Moten Launcher version 2.0, Moten Technologies, Puteaux, France). EMG and MMG signals were band-pass filtered between 20–500 Hz [5] and 5–100 Hz [16], respectively.

For the isometric testing procedure, each trial was segmented to select the 33% middle part of the contraction time according to previous studies [7,17–19].

The amplitudes of EMG and MMG signals were calculated by taking the root mean square (RMS) value of the segmented portion. The RMS values from EMG and MMG for each muscle were normalized from the average RMS values obtained during the three MVC trials [20]. The studied variables are referred to as *RMS<sub>EMG</sub>* and *RMS<sub>MMG</sub>*. To estimate the signal frequency content, a Fast Fourier Transform was performed over the segmented data and the mean power frequency (MPF) was calculated. The studied variables are referred to as *MPF<sub>EMG</sub>* and *MPF<sub>MMG</sub>*. Despite the non-stationarity of EMG and MMG signals during isometric and dynamic contractions, Fourier-based methods were shown to be acceptable during moderate velocity actions for estimating their frequency content [21].
$RMS_{MMG}$ and $MPF_{MMG}$ were computed following Equations (1) and (2), respectively [22]:

$$RMS_{MMG} = \sqrt{RMS_{X}^{2} + RMS_{Y}^{2} + RMS_{Z}^{2}} \qquad (1)$$

$$MPF_{MMG} = \sqrt{MPF_{X}^{2} + MPF_{Y}^{2} + MPF_{Z}^{2}} \qquad (2)$$

where X, Y and Z represent the 3 axes of the accelerometer.

For the dynamic testing, the accelerometer was first low-pass filtered at 5 Hz to keep only the motion data (MOT), which is composed of low frequency accelerations [10]. We used the MOT signal of the accelerometer located on the BB muscle to compute the sensor's pitch angle $\theta$ following Equation (3):

$$\theta = \arctan_{2}(Z, \sqrt{X^{2} + Y^{2}}) \times -180/\pi \qquad (3)$$

where X, Y and Z represent the 3 axes of the accelerometer.

$\theta$ was used to segment phases of (1) pulling the load from mid-tibia level to hip level (from posture a to b), (2) pulling the load from hip level to shoulder level (from posture b to c), (3) pushing the load above shoulder level (from posture c to d), (4) lowering the load to shoulder level (from posture d to c), (5) lowering the load to hip level (from posture c to b) and (6) lowering the load to mid-tibia level (from posture b to a), see Figures 2 and 3. Like isometric testing, we used the RMS and MPF to calculate signal amplitude and frequency content over the segmented data. EMG and MMG data from dynamic contractions were not normalized since there were no dynamic MVC measurements.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Example of the low-passed pitch angle of the accelerometer sensor located over the BB muscle (graph (A)) with corresponding MMG and EMG 1-second rolling RMS envelopes for the DEL muscle (graph (B) and graph (C), respectively) in function of time for one repetition of the load lifting task at 75% MVC (see Figure 2).</figcaption></figure>

## 2.6. Statistical Analysis

For the isometric testing, repeatability tests of absolute $RMS_{EMG}$ and $RMS_{MMG}$ between intensities and muscles were conducted based on intraclass correlation coefficients ICC (3, k) over the three trials for each muscle, following methodological recommendations [23]. Repeated-measures analyses of variances (ANOVA) were carried out to study the combined effects of muscles (BB/TRI/DEL) and intensities (25/50/75/100% MVC) on $RMS_{EMG}$, $RMS_{MMG}$, $MPF_{EMG}$, and $MPF_{MMG}$.
For the dynamic testing, an ICC (3, k) was conducted in order to quantify the repeatability of the MOT data over intensity levels for each muscle. Because EMG and MMG signals were not normalized in dynamic condition, two-way repeated measure ANOVAS (intensity × phase) were conducted for RMS<sub>EMG</sub>, RMS<sub>MMG</sub>, MPF<sub>EMG</sub> and MPF<sub>MMG</sub> separately for each muscle. We excluded phases 1 and 6 from the analysis because the recruitment of the three upper limb muscles during these phases was negligible. Post hoc analyses were conducted using the Bonferroni correction. When an intensity × phase interaction was found, we focused on the paired comparisons only between the concentric and eccentric phases of the motion where each muscle was recruited the most. For the BB muscle, phases 2 and 5 were identified as the concentric and eccentric phases, respectively (see Figure 2). For TRI and DEL muscles, phases 3 and 4 were identified as concentric and eccentric phases, respectively (see Figure 2). Average MMG and EMG RMS amplitudes in relation to % MVC were examined using a linear regression model and Pearson's correlation test for each muscle and subject. Independent t-tests were performed to compare EMG and MMG corresponding determination coefficients R<sup>2</sup>. All statistical analysis was conducted on JASP version 0.17. An alpha of 0.05 was considered significant for all comparisons and correlations. All data are presented as mean ± standard deviation.

# 3. Results

## 3.1. Isometric Testing

### 3.1.1. Time Domain Analysis

The two-way repeated-measures ANOVA revealed a significant effect of intensity level on the normalized force data ($F(3,63) = 1855, p < 0.001, \eta = 0.98$). Post hoc analysis showed that all force levels were significantly different from one another. All EMG and MMG isometric testing results are displayed in Table 2.

<table><caption>**Table 2.** Isometric ICC scores and mean RMS values for each sensor, muscle and intensity level. ICC are displayed as point estimate—lower bound 95% confidence interval.</caption><thead><tr><th rowspan="2">Sensor</th><th rowspan="2">Muscle</th><th rowspan="2">Feature</th><th colspan="4">Intensity (%MVC)</th></tr><tr><th>25%</th><th>50%</th><th>75%</th><th>100%</th></tr></thead><tbody><tr><td rowspan="6">EMG</td><td rowspan="2">BB</td><td>ICC</td><td>0.96–0.93</td><td>0.97–0.94</td><td>0.99–0.98</td><td>-</td></tr><tr><td>RMS</td><td>22.1 ± 10.9 #†•</td><td>54.4 ± 23.7 *†•</td><td>80.0 ± 25.5 *#</td><td>90.2 ± 5.3 *#</td></tr><tr><td rowspan="2">TRI</td><td>ICC</td><td>0.98–0.97</td><td>0.99–0.98</td><td>0.996–0.992</td><td>-</td></tr><tr><td>RMS</td><td>15.0 ± 5.6 #†•</td><td>46.8 ± 17.6 *†•</td><td>81.7 ± 25.3 *#</td><td>89.5 ± 7.5 *#</td></tr><tr><td rowspan="2">DEL</td><td>ICC</td><td>0.99–0.98</td><td>0.99–0.98</td><td>0.99–0.985</td><td>-</td></tr><tr><td>RMS</td><td>33.8 ± 10.4 #†•</td><td>64.3 ± 18.5 *†•</td><td>92.2 ± 21.9 *#</td><td>93.6 ± 4.3 *#</td></tr><tr><td rowspan="6">MMG</td><td rowspan="2">BB</td><td>ICC</td><td>0.97–0.94</td><td>0.96–0.93</td><td>0.97–0.95</td><td>-</td></tr><tr><td>RMS</td><td>19.4 ± 8.6 #†•</td><td>60.8 ± 25.1 *†•</td><td>81.9 ± 24.7 *#</td><td>88.8 ± 6.6 *#</td></tr><tr><td rowspan="2">TRI</td><td>ICC</td><td>0.95–0.90</td><td>0.97–0.94</td><td>0.95–0.91</td><td>-</td></tr><tr><td>RMS</td><td>11.6 ± 5.3 #†•</td><td>45.2 ± 16.7 *†•</td><td>92.5 ± 20.9 *#</td><td>85.8 ± 7.2 *#</td></tr><tr><td rowspan="2">DEL</td><td>ICC</td><td>0.98–0.97</td><td>0.98–0.96</td><td>0.97–0.94</td><td>-</td></tr><tr><td>RMS</td><td>28.3 ± 12.7 #†•</td><td>75.1 ± 20.3 *†•</td><td>100.0 ± 29.4 *#</td><td>88.3 ± 7.6 *#</td></tr></tbody></table>

For both EMG and MMG, all ICC point estimates, lower and upper 95% confidence intervals were greater than 0.9. There was no significant effect of intensity level on ICC ratings. The ANOVA on RMS<sub>EMG</sub> exhibited a significant effect for muscle ($F(2,42) = 7.65, p < 0.001, \eta^2 = 0.03$), intensity ($F(3,63) = 270, p < 0.001, \eta^2 = 0.75$) and a significant interaction ($F(6,126) = 2.29, p = 0.039, \eta^2 = 0.01$). Similarly, the ANOVA on RMS<sub>MMG</sub> showed a significant effect for muscle ($F(2,42) = 6.98, p = 0.001, \eta^2 = 0.03$), intensity ($F(3,63) = 333, p < 0.001, \eta^2 = 0.74$) and a significant interaction ($F(6,126) = 7.6, p < 0.001, \eta^2 = 0.01$).

<footnote>*: significant difference with 25% MVC; #: significant difference with 50% MVC; †: significant difference with 75% MVC; •: significant difference with 100% MVC ($p < 0.05$).</footnote>
$p < 0.001$, $\eta^2 = 0.02$). Pairwise comparisons highlighted that, for all muscles, all $RMS_{EMG}$ and $RMS_{MMG}$ values were significantly different from one another, except between 75% and 100% MVC.

### 3.1.2. Frequency Domain Analysis

Two-way ANOVA on $MPF_{EMG}$ showed a significant effect of muscle ($F(2,42) = 4$, $p = 0.024$, $\eta^2 = 0.11$) and intensity ($F(3,63) = 9.7$, $p < 0.001$, $\eta^2 = 0.06$) on $MPF_{EMG}$. Post hoc analysis revealed that $MPF_{EMG}$ at 75% MVC was lower than MPF at other intensities.

The ANOVA on $MPF_{MMG}$ yielded significant effects of muscle ($F(2,42) = 4.2$, $p < 0.021$, $\eta^2 = 0.07$) and intensity ($F(3,63) = 11.3$, $p < 0.001$, $\eta^2 = 0.11$). Post hoc analysis revealed that $MPF_{MMG}$ at 25% and 50% MVC were larger than 75% and 100% MVC.

## 3.2. *Dynamic Testing*

### 3.2.1. Time Domain Analysis

The ICC results to assess the repeatability of MOT data between intensity levels for each muscle were: 0.96 for the BB, 0.99 for the DEL and 0.98 for the TRI.

All two-way repeated-measures ANOVA on $RMS_{MMG}$ and $RMS_{EMG}$ revealed a significant effect of intensity, phase and a significant intensity × phase interaction for each muscle (view Table 3).

<table><caption>**Table 3.** MMG and EMG RMS amplitudes two-way repeated measures ANOVA results for each muscle from dynamic testing.</caption><thead><tr><th>Sensor</th><th>Muscle</th><th>Factor</th><th>df</th><th>F</th><th>p</th><th>η<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="12">MMG</td><td rowspan="3">BB</td><td>Intensity</td><td>(2, 36)</td><td>96</td><td>&lt;0.001</td><td>0.40</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>36</td><td>&lt;0.001</td><td>0.29</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>5.5</td><td>&lt;0.001</td><td>0.02</td></tr><tr><td rowspan="3">TRI</td><td>Intensity</td><td>(2, 36)</td><td>84</td><td>&lt;0.001</td><td>0.30</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>33</td><td>&lt;0.001</td><td>0.35</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>7.8</td><td>&lt;0.001</td><td>0.02</td></tr><tr><td rowspan="3">DEL</td><td>Intensity</td><td>(2, 36)</td><td>96</td><td>&lt;0.001</td><td>0.40</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>36</td><td>&lt;0.001</td><td>0.29</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>5.5</td><td>&lt;0.001</td><td>0.02</td></tr><tr><td rowspan="12">EMG</td><td rowspan="3">BB</td><td>Intensity</td><td>(2, 36)</td><td>37</td><td>&lt;0.001</td><td>0.16</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>21</td><td>&lt;0.001</td><td>0.32</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>5.5</td><td>&lt;0.001</td><td>0.04</td></tr><tr><td rowspan="3">TRI</td><td>Intensity</td><td>(2, 36)</td><td>5.2</td><td>0.01</td><td>0.05</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>3.8</td><td>0.015</td><td>0.13</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>4.58</td><td>&lt;0.001</td><td>0.01</td></tr><tr><td rowspan="3">DEL</td><td>Intensity</td><td>(2, 36)</td><td>32</td><td>&lt;0.001</td><td>0.12</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>58</td><td>&lt;0.001</td><td>0.59</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>13</td><td>&lt;0.001</td><td>0.02</td></tr></tbody></table>

### Intensity Levels

Post hoc analysis highlighted that for the BB muscle, $RMS_{MMG}$ values were significantly different between intensity levels except between 25% and 50% MVC in phase 5 (eccentric contraction). For the TRI muscle, $RMS_{MMG}$ values were different between intensities except between 50% and 75% MVC in phase 3 (concentric contraction). For the DEL muscle, $RMS_{MMG}$ values were different between all intensities for both concentric and eccentric phases of contraction.

For EMG, RMS values of the BB muscle were different between 25% and 75% MVC for phase 2 and 5 (concentric and eccentric contractions); however, there were no RMS differences between 50% and 75% MVC for phase 2 and between 25% and 50% MVC for phase 5. For TRI, RMS values were only different between 25% and 75% MVC for phase
3 (concentric contraction). For DEL, RMS values were different from each other except between 50% and 75% MVC for phase 3 (concentric contraction) (see Figure 4).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Absolute *RMS*<sub>EMG</sub> and *RMS*<sub>MMG</sub> values as a function of phase (between 2 and 5) for each intensity level and muscle (row graphics).</figcaption></figure>

### Concentric vs. Eccentric

For MMG, there was no difference in RMS between concentric and eccentric phases for the BB muscle for all intensity levels. For the TRI and DEL muscles, RMS during the concentric contraction (phase 3) was lower than the eccentric part (phase 4) for 50% and 75% MVC.

For EMG, there was no difference in RMS between concentric and eccentric phases for the BB and TRI muscles for all intensity levels. For the DEL muscle, RMS values during its concentric contraction was greater than the eccentric contraction for all intensity levels (*p* < 0.001).

### EMG and MMG RMS vs. % MVC Relationship

For BB and TRI muscles, MMG $R^2$ were significantly greater than EMG (*p* = 0.015 and *p* < 0.01, respectively). For the DEL muscle, there was no difference between EMG and MMG mean determination coefficients (see Figure 5).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Mean EMG and MMG RMS amplitudes vs. % MVC linear relationship for each muscle during the load lifting-task. Y-axis RMS amplitude unit is in mV for EMG and in mg for MMG.</figcaption></figure>

### 3.2.2. Frequency Domain Analysis

All two-way repeated-measures ANOVA on *MPF<sub>MMG</sub>* and *MPF<sub>EMG</sub>* are presented in Table 4).

<table><caption>**Table 4.** MMG and EMG MPF two-way repeated measures ANOVA results for each muscle from dynamic testing. Only significant effects are displayed.</caption><thead><tr><th>Sensor</th><th>Muscle</th><th>Factor</th><th>df</th><th>F</th><th>p</th><th>η<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="12">MMG</td><td rowspan="3">BB</td><td>Intensity</td><td>(2, 36)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>61</td><td>&lt;0.001</td><td>0.67</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">TRI</td><td>Intensity</td><td>(2, 36)</td><td>14</td><td>&lt;0.001</td><td>0.04</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>53</td><td>&lt;0.001</td><td>0.61</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>10</td><td>&lt;0.001</td><td>0.04</td></tr><tr><td rowspan="3">DEL</td><td>Intensity</td><td>(2, 36)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>4</td><td>&lt;0.001</td><td>0.05</td></tr><tr><td rowspan="12">EMG</td><td rowspan="3">BB</td><td>Intensity</td><td>(2, 36)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>9</td><td>&lt;0.001</td><td>0.24</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">TRI</td><td>Intensity</td><td>(2, 36)</td><td>56</td><td>&lt;0.001</td><td>0.20</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>17</td><td>&lt;0.001</td><td>0.28</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>7</td><td>&lt;0.001</td><td>0.04</td></tr><tr><td rowspan="3">DEL</td><td>Intensity</td><td>(2, 36)</td><td>14</td><td>&lt;0.001</td><td>0.05</td></tr><tr><td>Phase</td><td>(3, 54)</td><td>37</td><td>&lt;0.001</td><td>0.55</td></tr><tr><td>Intensity × Phase</td><td>(6, 108)</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

#### Intensity Levels

Post hoc analysis indicated that for MMG, TRI MPF decreased with increasing load. MPF at 75% MVC was significantly lower than 25 and 50% MVC. EMG MPF for the TRI and DEL highlighted similar results where MPF decreased significantly between all intensity levels.

### 4. Discussion

The aim of this study was to assess the feasibility of MMG measurements in complex dynamic situations by comparing MMG to EMG, which is the reference system to measure
muscle activity, during isometric contractions and a load lifting task at different intensities. Twenty-two male participants took part in this study, consisting of submaximal isometric contractions at 25%, 50% and 75% MVC of the BB, TRI and DEL muscles. Participants then performed repetitions of a lifting task from the ground and raising it above shoulder level at three relative intensity levels. MMG and EMG signals were acquired for the Biceps, Triceps and Deltoid of the dominant arm.

## 4.1. *Isometric Testing*

### 4.1.1. Time Domain Analysis

Overall, MMG repeatability scores obtained in this study can be considered excellent [24] and were higher than results previously related to MMG studies [25,26]. As an example, a relatively recent study obtained an ICC score of 0.79 with a microphone MMG [27].

$RMS_{MMG}$ results obtained in this study agree with previous investigations and the most recent MMG review article [28]. In this review, it has been shown that $RMS_{MMG}$ increases with force up to 80% MVC and reaches a plateau or even decreases to 100% MVC in isometric condition. Authors have suggested that an increase in $RMS_{MMG}$ with force is due to the increasing number of active motor units to produce force. However, at higher intensities, an increase in muscle stiffness, which is a function of attached cross-bridges and fusion of motor-unit twitches, causes an MMG amplitude plateau or a drop [28]. Similar behavior was shared with EMG, as the amplitude increases to 75% MVC and plateaus at MVC in the current study. EMG results were also in agreement with previous studies and reviews [28,29]. These results suggest that at approximately 75% MVC, the major part of motor units (MU) is recruited [28].

The evolution of MMG and EMG amplitude as a function of force level has been showed to be muscle-specific [17] but accounted only for a very small portion of the total variance ($\eta^2$ = 0.02 and $\eta^2$ = 0.01 in the current results for MMG and EMG, respectively).

### 4.1.2. Frequency Domain Analysis

$MPF_{EMG}$ decreased from 25% to 75% MVC and increased to the initial value at 100% MVC. This result agrees with previous investigations in step isometric contraction for which authors found a similar $MPF_{EMG}$ decreasing trend from 20% to 80% MVC [30]. $MPF_{MMG}$ decreased from low intensities (25–50% MVC) to high intensities (75–100% MVC).

For both EMG and MMG, MPF correlates with the motor unit firing rate. As MPF decreased with increasing intensity, we can conclude that the increase in force was modulated by the increasing number of active motor units and by the decrease in the global motor units firing rate. Other authors suggested that muscle stiffness and/or intramuscular pressure, with increasing intensity, could impair the firing rate of motor units [31]. This impairment has been shown to stabilize or even decrease the $MPF_{MMG}$ with increasing force in step contractions.

## 4.2. *Dynamic Testing*

### 4.2.1. Time Domain Analysis

Thanks to the repeatability analysis of the MOT data and the ICC scores, it was verified a posteriori that there were no significant kinematic differences between the different-load trials that could lead to misinterpretations of the MMG and EMG signals.

In the current study, we observed a significant effect of intensity, where $RMS_{EMG}$ and $RMS_{MMG}$ increased linearly with load. This result agrees with those obtained in isometric testing and in previous investigations during mono-articular isokinetic condition [7,19,32].

It was found that the $EMG_{RMS}$ of the solicited muscle was greater in the concentric part of the contraction than the eccentric part, although an inverse behavior was found for MMG RMS. These results agree with previous findings in isokinetic condition for which authors concluded that muscular tremor was greater during the eccentric part of the contraction [33]. These results are also consistent with those obtained in an unloaded
squatting task for which the RMS<sub>MMG</sub> of the rectus femoris was greater entering (eccentric phase) compared to exiting the squat (concentric phase) [3].

Finally, the main result of this study concerns the fact that MMG was more sensitive to changes in load than EMG in dynamic condition. MMG amplitude exhibited greater systematic effect sizes of intensity in ANOVAs and significantly greater correlations with load for the BB and TRI muscles. A similar result has been found during incremental cycling, for which authors exhibited that “MMG amplitude more closely reflects changes in power output” [34].

### 4.2.2. Frequency Domain Analysis

In this study, $MPF_{EMG}$ significantly decreased with increasing load for the TRI and DEL muscles. The same MPF vs. load relationship was found for MMG but only for the TRI muscle. These results agree with those obtained in isometric condition where MPF at 75% MVC was lower than 25% MVC. Moreover, MPF vs. force responses are muscle-specific, as there was no significant increase nor decrease for the BB with both EMG and MMG. This result agrees with previous studies in isokinetic condition [7,19,32].

### 4.3. *Limits*

One of the limits of this study is the absence of EMG and MMG normalization in dynamic conditions. Because MMG amplitude is greater in dynamic condition compared to isometric condition [35], we could not have used RMS<sub>MMG</sub> from MVC values as a reference to normalize signals. Another limit is that in dynamic condition, MMG may be polluted with motion artifacts despite high-pass filtering above 5 Hz [10]. Even if an ulterior verification was done on kinematics between intensity levels, artifacts such as muscular tremor or higher-frequency motion harmonics might have altered the MMG signal. In future work, to address these limits, we could improve MMG and EMG signal processing to better filter motion artifacts in dynamic conditions using EMD (Empirical Mode Decomposition)-based algorithms, which showed more effective results than band-pass filtering recently [11,12,36]. Other limits of the current study are related to gender as only male subjects have been selected. Moreover, the age of the participants (25.5 ± 5.3) may also reduce the scope of this work. To extend outcomes of this study in the future, it would be interesting to include more participants of different ages, both men and women.

## 5. Conclusions

To the best of our knowledge, this study was the first to investigate MMG and EMG during a functional load lifting task at different intensities. This study widens the scope of MMG knowledge in the context of functional activities. During isometric conditions, EMG and MMG had excellent repeatability scores, as well as similar RMS and MPF vs. force relationship. During a dynamic situation corresponding to a load-lifting task, RMS from the MMG signal was more sensitive to changes in load than EMG signal. Conversely, MPF from the EMG signal was more sensitive to changes in load. Altogether, these results sustain the feasibility of using MMG measurement in complex functional activities, such as load lifting and carrying. Moreover, EMG and MMG signals may be used together to give complementary information concerning strategies of motor unit recruitment, in both time and frequency domains. This study highlighted the fact that MMG signal may be used for practical applications such as in sports science or physical ergonomics assessment. Further work is now necessary in terms of signal processing to improve the reliability of MMG signal during ecological conditions.
**Author Contributions:** Conceptualization, M.C., M.P., I.A.S. and N.V.; methodology M.C., M.P., I.A.S. and N.V.; software M.C.; validation M.P., I.A.S. and N.V.; formal analysis M.C., M.P., I.A.S. and N.V.; investigation M.C., M.P., I.A.S. and N.V.; resources M.P., N.V. and I.A.S.; data curation M.C.; writing—original draft preparation M.C.; writing—review and editing M.P., N.V. and I.A.S.; visualization M.C., M.P., I.A.S. and N.V.; supervision N.V., I.A.S. and M.P.; project administration M.C., M.P., I.A.S. and N.V.; funding acquisition M.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Université Paris-Saclay, (2021-287).

**Informed Consent Statement:** Written informed consent has been obtained from the participants to publish this paper.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author. The data are not publicly available due to commercial privacy policy.

**Acknowledgments:** We would like to thank Dorian Verdel and Jean Jeuvrey for their help in the experimental design.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Gazzoni, M.; Afsharipour, B.; Merletti, R. Surface EMG in Ergonomics and Occupational Medicine. In *Surface Electromyography: Physiology, Engineering, and Applications*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2016; pp. 361–391.

2. Hermens, H.J.; Freriks, B.; Disselhorst-Klug, C.; Rau, G. Development of recommendations for SEMG sensors and sensor placement procedures. *J. Electromyogr. Kinesiol. Off. J. Int. Soc. Electrophysiol. Kinesiol. 2000*, *10*, 361–374. [CrossRef]

3. Woodward, R.B.; Stokes, M.J.; Shefelbine, S.J.; Vaidyanathan, R. Segmenting Mechanomyography Measures of Muscle Activity Phases Using Inertial Data. *Sci. Rep.* **2019**, *9*, 5569.

4. Laferriere, P.; Chan, A.D.C.; Lemaire, E.D. Surface electromyographic signals using a dry electrode. In *Proceedings of the 2010 IEEE International Workshop on Medical Measurements and Applications*, Ottawa, ON, Canada, 30 April–1 May 2010; pp. 77–80. [CrossRef]

5. De Luca, C.J.; Donald Gilmore, L.; Kuznetsov, M.; Roy, S.H. Filtering the surface EMG signal: Movement artifact and baseline noise contamination. *J. Biomech.* **2010**, *43*, 1573–1579. [CrossRef]

6. Farina, D. Interpretation of the surface electromyogram in dynamic contractions. *Exerc. Sport Sci. Rev.* **2006**, *34*, 121–127. [CrossRef]

7. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic and electromyographic time and frequency domain responses during submaximal to maximal isokinetic muscle actions of the biceps brachii. *Eur. J. Appl. Physiol.* **2004**, *92*, 352–359. [CrossRef]

8. Ibitoye, M.O.; Hamzaid, N.A.; Zuniga, J.M.; Abdul Wahab, A.K. Mechanomyography and muscle function assessment: A review of current state and prospects. *Clin. Biomech.* **2014**, *29*, 691–704. [CrossRef]

9. Islam, M.A.; Sundaraj, K.; Ahmad, R.B.; Ahamed, N.U.; Ali, M.A. Mechanomyography Sensor Development, Related Signal Processing, and Applications: A Systematic Review. *IEEE Sens. J.* **2013**, *13*, 2499–2516.

10. Posatskiy, A.O.; Chau, T. The effects of motion artifact on mechanomyography: A comparative study of microphones and accelerometers. *J. Electromyogr. Kinesiol.* **2012**, *22*, 320–324. [CrossRef]

11. Wang, D.; Wu, H.; Xie, C.; Gao, L. Suppression of Motion Artifacts in Multichannel Mechanomyography Using Multivariate Empirical Mode Decomposition. *IEEE Sens. J.* **2019**, *19*, 5732–5739.

12. Li, Z.; Gao, L.; Lu, W.; Wang, D.; Cao, H.; Zhang, G. A Novel Noise Suppression and Artifact Removal Method of Mechanomyography Based on RLS, IGWO-VMD, and CEEMDAN. *J. Sens.* **2022**, *2022*, e4239211.

13. Yu, J.; Zhang, Y.; Xia, C. Study of gait pattern recognition based on fusion of mechanomyography and attitude angle signal. *J. Mech. Med. Biol.* **2020**, *20*, 1950085.

14. Silva, J.; Heim, W.; Chau, T. A Self-Contained, Mechanomyography-Driven Externally Powered Prosthesis. *Arch. Phys. Med. Rehabil.* **2005**, *86*, 2066–2070. [CrossRef]

15. Evetovich, T.K.; Housh, T.J.; Johnson, G.O.; Smith, D.B.; Ebersole, K.T.; Perry, S.R. Gender comparisons of the mechanomyographic responses to maximal concentric and eccentric isokinetic muscle actions. *Med. Sci. Sport. Exerc.* **1998**, *30*, 1697–1702. [CrossRef]

16. Ibitoye, M.O.; Hamzaid, N.A.; Zuniga, J.M.; Hasnan, N.; Wahab, A.K.A. Mechanomyographic parameter extraction methods: An appraisal for clinical applications. *Sensors* **2014**, *14*, 22940–22970. [CrossRef]

17. Ebersole, K.T.; Housh, T.J.; Johnson, G.O.; Evetovich, T.K.; Smith, D.B.; Perry, S.R. MMG and EMG responses of the superficial quadriceps femoris muscles. *J. Electromyogr. Kinesiol.* **1999**, *9*, 219–227. [CrossRef]
18. Evetovich, T.K.; Boyd, J.C.; Drake, S.M.; Eschbach, L.C.; Magal, M.; Soukup, J.T.; Webster, M.J.; Whitehead, M.T.; Weir, J.P. Effect of moderate dehydration on torque, electromyography, and mechanomyography. *Muscle Nerve* **2002**, *26*, 225-231. [CrossRef]

19. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Comparison of a piezoelectric contact sensor and an accelerometer for examining mechanomyographic amplitude and mean power frequency versus torque relationships during isokinetic and isometric muscle actions of the biceps brachii. *J. Electromyogr. Kinesiol.* **2006**, *16*, 324–335. [CrossRef]

20. Coburn, J.W.; Housh, T.J.; Weir, J.P.; Malek, M.H.; Cramer, J.T.; Beck, T.W.; Johnson, G.O. Mechanomyographic responses of the vastus medialis to isometric and eccentric muscle actions. *Med. Sci. Sport. Exerc.* **2004**, *36*, 1916–1922. [CrossRef]

21. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Comparison of Fourier and wavelet transform procedures for examining the mechanomyographic and electromyographic frequency domain responses during fatiguing isokinetic muscle actions of the biceps brachii. *J. Electromyogr. Kinesiol.* **2005**, *15*, 190–199. [CrossRef]

22. Nogueira-Neto, G.; Scheeren, E.; Krueger, E.; Nohama, P.; Button, V.L.S.N. The Influence of Window Length Analysis on the Time and Frequency Domain of Mechanomyographic and Electromyographic Signals of Submaximal Fatiguing Contractions. *Open J. Biophys.* **2013**, *3*, 178.

23. Shrout, P.E.; Fleiss, J.L. Intraclass Correlations: Uses in Assessing Rater Reliability. *Psychol. Bull.* **1979**, *86*, 420. [CrossRef]

24. Koo, T.K.; Li, M.Y. A Guideline of Selecting and Reporting Intraclass Correlation Coefficients for Reliability Research. *J. Chiropr. Med.* **2016**, *15*, 155–163. [CrossRef] [PubMed]

25. Akataki, K.; Mita, K.; Itoh, Y. Repeatability study of mechanomyography in submaximal isometric contractions using coefficient of variation and intraclass correlation coefficient. *Electromyogr. Clin. Neurophysiol.* **1999**, *39*, 161–166. [PubMed]

26. Herda, T.J.; Ryan, E.D.; Beck, T.W.; Costa, P.B.; DeFreitas, J.M.; Stout, J.R.; Cramer, J.T. Reliability of mechanomyographic amplitude and mean power frequency during isometric step and ramp muscle actions. *J. Neurosci. Methods* **2008**, *171*, 104–109. [CrossRef] [PubMed]

27. Meagher, C.; Franco, E.; Turk, R.; Wilson, S.; Steadman, N.; McNicholas, L.; Vaidyanathan, R.; Burridge, J.; Stokes, M. New advances in mechanomyography sensor technology and signal processing: Validity and intrarater reliability of recordings from muscle. *J. Rehabil. Assist. Technol. Eng.* **2020**, *7*, 2055668320916116. [CrossRef] [PubMed]

28. Talib, I.; Sundaraj, K.; Lam, C.K.; Sundaraj, S. A systematic review of muscle activity assessment of the biceps brachii muscle using mechanomyography. *J. Musculoskelet. Neuronal Interact.* **2018**, *18*, 446–462.

29. Madeleine, P.; Bajaj, P.; Søgaard, K.; Arendt-Nielsen, L. Mechanomyography and electromyography force relationships during concentric, isometric and eccentric contractions. *J. Electromyogr. Kinesiol.* **2001**, *11*, 113–121. [CrossRef]

30. Bilodeau, M.; Arsenault, A.B.; Gravel, D.; Bourbonnais, D. EMG power spectra of elbow extensors during ramp and step isometric contractions. *Eur. J. Appl. Physiol. Occup. Physiol.* **1991**, *63*, 24–28. [CrossRef]

31. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Cramer, J.T.; Weir, J.P.; Coburn, J.W.; Malek, M.H. Does the frequency content of the surface mechanomyographic signal reflect motor unit firing rates? A brief review. *J. Electromyogr. Kinesiol. Off. J. Int. Soc. Electrophysiol. Kinesiol.* **2007**, *17*, 1–13. [CrossRef]

32. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic and electromyographic amplitude and frequency responses during fatiguing isokinetic muscle actions of the biceps brachii. *Electromyogr. Clin. Neurophysiol.* **2004**, *44*, 431–441.

33. Qi, L.; Wakeling, J.M.; Ferguson-Pell, M. Spectral properties of electromyographic and mechanomyographic signals during dynamic concentric and eccentric contractions of the human biceps brachii muscle. *J. Electromyogr. Kinesiol.* **2011**, *21*, 1056–1063. [CrossRef] [PubMed]

34. Perry, S.R.; Housh, T.J.; Weir, J.P.; Johnson, G.O.; Bull, A.J.; Ebersole, K.T. Mean power frequency and amplitude of the mechanomyographic and electromyographic signals during incremental cycle ergometry. *J. Electromyogr. Kinesiol.* **2001**, *11*, 299–305. [CrossRef] [PubMed]

35. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic amplitude and mean power frequency versus torque relationships during isokinetic and isometric muscle actions of the biceps brachii. *J. Electromyogr. Kinesiol.* **2004**, *14*, 555–564. [CrossRef] [PubMed]

36. Zhang, Y.; Xu, P.; Li, P.; Duan, K.; Wen, Y.; Yang, Q.; Zhang, T.; Yao, D. Noise-assisted multivariate empirical mode decomposition for multichannel EMG signals. *Biomed. Eng. Online* **2017**, *16*, 107. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.