Next Article in Journal
Location-Awareness for Failure Management in Cellular Networks: An Integrated Approach
Next Article in Special Issue
Measuring Gait Velocity and Stride Length with an Ultrawide Bandwidth Local Positioning System and an Inertial Measurement Unit
Previous Article in Journal
Validation of Using Smartphone Built-In Accelerometers to Estimate the Active Energy Expenditures of Full-Time Manual Wheelchair Users with Spinal Cord Injury
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 21
Issue 4
10.3390/s21041499
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Data-Driven Approach to Predict Fatigue in Exercise Based on Motion Data from Wearable Sensors or Force Plate

by
Yanran Jiang
1,*,
Vincent Hernandez
2,
Gentiane Venture
2,
Dana Kulić
1 and
Bernard K. Chen
1
1
Mechanical and Aerospace Department, Monash University, Melbourne, VIC 3800, Australia
2
Department of Mechanical Systems Engineering, Tokyo University of Agriculture and Technology, Tokyo 184-0012, Japan
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(4), 1499; https://doi.org/10.3390/s21041499
Submission received: 19 January 2021 / Revised: 6 February 2021 / Accepted: 9 February 2021 / Published: 22 February 2021
(This article belongs to the Special Issue Sensor-Based Measurement of Human Motor Performance)
Browse Figures
Versions Notes

Abstract

:
Fatigue increases the risk of injury during sports training and rehabilitation. Early detection of fatigue during exercises would help adapt the training in order to prevent over-training and injury. This study lays the foundation for a data-driven model to automatically predict the onset of fatigue and quantify consequent fatigue changes using a force plate (FP) or inertial measurement units (IMUs). The force plate and body-worn IMUs were used to capture movements associated with exercises (squats, high knee jacks, and corkscrew toe-touch) to estimate participant-specific fatigue levels in a continuous fashion using random forest (RF) regression and convolutional neural network (CNN) based regression models. Analysis of unseen data showed high correlation (up to 89%, 93%, and 94% for the squat, jack, and corkscrew exercises, respectively) between the predicted fatigue levels and self-reported fatigue levels. Predictions using force plate data achieved similar performance as those with IMU data; the best results in both cases were achieved with a convolutional neural network. The displacement of the center of pressure (COP) was found to be correlated with fatigue compared to other commonly used features of the force plate. Bland–Altman analysis also confirmed that the predicted fatigue levels were close to the true values. These results contribute to the field of human motion recognition by proposing a deep neural network model that can detect fairly small changes of motion data in a continuous process and quantify the movement. Based on the successful findings with three different exercises, the general nature of the methodology is potentially applicable to a variety of other forms of exercises, thereby contributing to the future adaptation of exercise programs and prevention of over-training and injury as a result of excessive fatigue.

1. Introduction

Physical activity can improve health and well-being, reduce the risk of many diseases, and improve the quality of life [1]. However, a large number of people suffer injuries during exercise [2]. A major contributing factor of exercise injuries is fatigue [2]. Fatigue caused by repeated movement accumulates over time and may exceed the muscle tissues’ tolerance, contributing to musculoskeletal disorders (MSDs) [2,3,4]. Thus, monitoring and predicting fatigue are important to reduce the risk of injuries. In the context of sports training, fatigue estimation can be used by coaches and physical therapists to avoid high levels of fatigue, which may adversely impact training and hinder performance in competition. In the context of rehabilitation, many patients are instructed to perform rehabilitation exercises at home by themselves. Without the therapist’s instructions and feedback on their movements, there are greater risks of secondary injury. Therefore, actively monitoring the onset of fatigue could provide important feedback in sports training, competition, and rehabilitation [4].
Human activity recognition (HAR) is a broad research field that involves the identification of various human activities or gestures and more detailed knowledge about human activities (e.g., the quality of motions, emotions, and gender) based on sensor data [5,6,7]. Recently, many novel data processing and analysis methods have been applied to HAR due to the introduction of wearable, low cost, low power sensors and live streaming of data [5]. Meanwhile, advances in computer vision, machine learning, and artificial intelligence have enabled HAR to be widely used in athletic competition, healthcare, and elderly care applications [5,8,9,10]. Even though much research has been conducted on human movement analysis and action recognition [5,6,7], the study of automatic fatigue prediction or estimation is somewhat limited.
Typically, biomechanical variables are measured using motion capture for the kinematics of body segments, electromyography for muscle activity, and plantar pressure measurements for detecting stepping [11]. Common techniques for detecting fatigue are to measure muscle activity, e.g., surface electromyogram (sEMG) [12], and the kinematics of joint angles, e.g., optical motion capture [1]. However, sEMG has limitations, e.g., the sensors may lose contact over time, particularly during vigorous exercise, and can only measure the activity of the particular muscle to which the sensor is attached. For optical motion capture, limited capture area and occlusion are key issues because reflective markers can be hidden from the camera and additional performers will increase occlusion in team sports. Even though there are some technologies that use structural lighting [13] and multi-camera motion capture systems [14,15] that can exclude blind spots, they are difficult to apply to outdoor activities. The multi-camera system usually takes longer to set up due to the large amount of equipment and is less flexible. On the contrary, wearable IMUs are small, lightweight, and robust to occlusions and interference. They do not restrict body movements and allow a participant to perform various tasks in arbitrary environments. On the other hand, force plates (FPs) are easy to use and do not require any equipment on the body, therefore saving much time in the setup of experiments. They also help to record the motion abnormalities of lower body segments during specific activities [16]. Foot plantar sensors are also seldom discussed in the application of fatigue estimation. Therefore, force plates and IMUs are selected to measure exercise-induced fatigue, and their performances are compared.
A number of studies [17,18,19,20] have been conducted to detect binary fatigue (fatigue vs. non-fatigue) by using various invasive or noninvasive devices including IMUs and sEMG. However, these models can only represent a simplistic process of fatigue development, and an essential early intervention is impossible before the athlete is deemed to be fatigued. Limited research [1,20,21] has been conducted to monitor gradual and continuous changes in fatigue levels. The change from low levels of fatigue to high levels of fatigue is a continuous process, which may take place gradually, and the variations of human motion due to fatigue are fairly small during exercise, making it difficult to detect continuous changes in fatigue [17]. Moreover, the associated literature [17,18,19,20] only focuses on one task or exercise type (running only or jumping only) and lacks the ability to generalize the fatigue models to different exercises.
The onset of muscle fatigue is complex and may depend on personal fitness level, health conditions, types of exercise, and gender [8]. The observed changes in movement or muscle electrical activity may not be consistent among all individuals. In other words, the expression of fatigue is more likely to be person-dependent.
In this work, a data-driven approach is investigated to estimate the onset of fatigue using force plate or IMU measurements in three different exercises (squat, high knee jack, and corkscrew toe-touch). The proposed framework is the first of its kind to model the continuous increase in fatigue based on each single repetition (rep) of exercise and the first to be evaluated on multiple exercises, showing its potential generalizability. Unlike most previous studies, which focused on extreme fatigue and non-fatigue detection, here, the approach is to monitor the fatigue accumulated status with continuous feedback. In the field of human motion recognition, the deep learning based framework makes a great contribution to understanding the quality of exercise for rehabilitation in a continuous process, automatically, without specific domain knowledge.

2. Related Work

In this section, learning approaches for fatigue estimation are discussed.
Various research has been conducted to evaluate fatigue through data from different sensors and different methodologies, and the choice of sensors depends on the purpose or objectives of each study. Most studies (e.g., [12,22,23,24,25,26,27]) were conducted using surface electromyogram sensors (sEMG) to determine the physiological status of a muscle due to the activity. For example, Chattopadhyay et al. proposed an exhaustive set of features from the sEMG signals and analyzed the variability between subjects and between trials [25]. Dong et al. proposed a method to evaluate the overall fatigue of human body movement based on combined sEMG and accelerometer signals and introduced a “forgetting factor” and fatigue level fusion coefficient to combine different localized muscle fatigue estimates with the overall fatigue level [4]. However, sEMG may lose contact over time, particularly for dynamic exercise, and can only measure the particular muscle to which a sensor is attached. It cannot be used in the real world.
Inertial measurement units (IMUs) have also been used to analyze fatigue based on kinematic movement, especially in gait [18], running [19,20,28], and sprinting [29]. Those results indicated the capability of IMUs to provide reliable and accurate measurements of temporal parameters during exercise. The most relevant studies to our research are [19,20]. Buckley et al. predicted subject-dependent and subject-independent fatigue levels (non-fatigued and fatigue status) through data from a single IMU [19]. In the experiment, running four-hundred meters at a natural pace was considered as a non-fatigue status, and a subsequent beep test (also known as the “multi-stage fitness test”) was used to induce fatigue. The results showed that a single IMU on the right shank had better performance than on the lumbar spine when assessing subject-independent fatigue estimation, and the subject-dependent classifier had higher accuracy than the subject-independent classifier. Stohrmann et al. [20] monitored human fatigue by extracting kinematic parameters from wearable sensor data and investigated the kinematic changes evoked by fatigue during running. Twenty-one runners of different skill levels performed experiments on a treadmill and conventional outdoor track. Their findings showed that kinematic changes were related to fatigue for all runners, and fatigue was dependent on participants’ running technique [20]. To date, there are only a few studies focusing on continuous movement changes induced by fatigue, e.g., [1,20,21]. Ramos et al. presented a machine learning system to evaluate fatigue using electromyographic (EMG) and heart rate variability (HRV) measurements [21]. The approach showed a potential to implement a combination of a dimensionless (0–1) global fatigue descriptor to reflect the onset of fatigue.
Additionally, in the past research, a force plate was generally used to measure the postural stability performance [30,31], and some studies assessed counter-movement jump performance through the investigation of the ground reaction force-time profile [32,33,34]. The study [34] showed that fatigue due to the exercise of the calf-muscles of one leg could influence the body balance in the short term, and this can be measured by a force plate and an accelerator, which indicates the potential to distinguish human motion status through these two tools. In view of automatic recognition technology, different classifiers have been applied, such as the fuzzy logic (FL) classifier, the random forest classifier [19], the support vector machine (SVM) classifier [12,18], and the hidden Markov model (HMM) classifier [35]. However, these currently applied algorithms can only classify fatigue and non-fatigued status, and it it challenging to monitor a gradual and continuous changing of fatigue level. The methods also vary according to various exercises and muscle types. Different exercises may cause fatigue in different human body parts, and the measures of fatigue highly rely on the types of sensor and sensor positions. Thus, this paper aims to develop related algorithms to detect continuous fatigue changes and generalize fatigue models for various exercises.

3. Materials and Methods

Subject-dependent recognition of fatigue for three different exercises, squat, high knee jack, and corkscrew toe-touch, is investigated. Figure 1 presents an overview of the four steps in the proposed methodology for predicting fatigue levels.

3.1. Data Collection

The fatigue dataset was collected using motion analysis motion capture (Motion Analysis, Santa Rosa, CA, USA) (sample rate: 100 Hz), Xsens IMUs (Xsens Technologies B.V., El Segundo, CA, USA) (sample rate: 240 Hz), and AMTI force plates (AMTI, Watertown, MA, USA) (sample rate: 100 Hz) simultaneously while performing full body exercises with one or two stationary feet. Figure 2 shows the distribution of the fatigue levels reported by participants in the experiment. A set of 32 reflective markers and 17 IMU sensors was attached on the human body (Figure 3). The dataset was collected at the University of Waterloo with ethics approval for the experiments. Fourteen healthy adults participated in the study (Table 1). The body weights were extracted from the force plate measurements. Before conducting the experiment, each participant was given a detailed explanation of the study and provided informed consent.
The participants were asked to perform three different exercises: squats, high-knee jacks, and corkscrew toe-touch, as described in Table 2. During each exercise, participants were asked to perform repeated sets of exercises to exhaustion. Each set consisted of 5 repetitions of an exercise. After each exercise set, participants were instructed to assess their fatigue level on a scale of 1 to 10 with “10” as “you are so tired you can’t stand up anymore”. The total number of sets of each exercise for each participant is summarized in Table 3, ranging from 3 to 52. This variation may depend on fitness level, age, gender, emotion, and pre-exercise fatigue level. Meanwhile, after each set, participants could take a short break to recover and then continue to the next exercise set.

3.2. Feature Extraction

The resultant accelerations and angular velocities were extracted from five IMU sensors on the lower body and selected as IMU features (sensors: right upper leg, right lower leg, left upper leg, left lower leg, and pelvis) for squat and high knee jack exercises. Since the corkscrew exercise also includes waist and hand motion, five sensors on the upper body were selected to measure their acceleration data (sensors: left upper arm, right upper arm, right forearm, left forearm, torso). In this study, the coordinates of the center of pressure (COP) were computed for a complete repetition according to Equations (1) and (2) and selected as input features instead of ground reaction force and moments (GRF&Ms). As the participant sways during exercise, the COP position varies over time and is one of the most commonly used features to quantify a person’s postural sway in both clinical and research contexts [36]. Comparing with the GRF&Ms, the COP features are more robust to inter-trial variability since they can remove the effect of the variability of different standing positions on the force plate. The complete feature set of the force plate used for training was the coordinates of COP and vertical ground reaction forces. Previous studies [36,37] have confirmed the reliability and informativeness of this choice of the feature set for postural stability.
C O P X = − M Y + ( F X × d Z ) F Z
C O P Y = − M X + ( F Y × d Z ) F Z
F X , F Y , and F Z respectively represent the components of the ground reaction force in which the X axis is in the mediolateral direction, the Y axis is in the anterior-posterior direction, and the Z axis in the superior-inferior direction. M X and M Y are components of the moments of force, and d Z is the depth of the force plate.
Other related COP features were also computed accordingly, including standard variance, variance, skew, kurtosis, area of the 95% ellipse, and first derivative and second derivative of COP positions. A explanation of the proposed features is provided in Table 4.

3.3. Data Preprocessing

Before using data in the conventional machine learning or deep learning model, it was necessary to preprocess the experimental data.
Segmentation: The continuous sequences of motion data obtained from the markers were divided manually into individual repetitions of squats by finding the segment points (i.e., the time step indicating the start and end of each repetition signal). Each repetition was labeled by the corresponding subjective fatigue level.
Synchronization: Since different sampling rates were used for the IMU data (240 Hz) and for the force data and marker data (100 Hz), a dynamic warping approach [38] was applied to synchronize the two different frequency time series data. This approach temporally aligns the two signals by identifying the alignment when the RMSE between two signals is minimized.
Temporal normalization: The lengths of the repetition segments were unevenly distributed due to the varying activity time of each participant. A fixed window size equal to 200 was considered, and each repetition was splined with a cubic spline interpolation.
Upsampling: Imbalanced datasets can skew the classifier towards the class that has the most samples [39]. Because only the beginning and end samples for each participant were in low and high fatigue while many samples were in between them, the highly unbalanced dataset was upsampled. The samples of each fatigue level were upsampled by adding additional samples to the minority class through random duplication, ensuring the final number of the minority class samples was equal to the number in the majority class.
Filtering: Butterworth low-pass filters with a cutoff frequency of 20 Hz were applied to all signals for noise reduction.
Standardization: The training, validation, and test data were standardized for each subject by computing the global mean μ and standard deviation σ from the training and validation sets. Then, each feature x m , n was standardized as follows:
x m , n = x m , n − μ n o r m σ n o r m
Training, validation, and test sets: The dataset of each participant was separated into a training set, a validation set, and a test set, which composed approximately 70%, 15%, and 15% of the motion data, respectively [40]. Before splitting the training/validation/test sets, the data sorted by time order were randomly shuffled, ensuring that the training/validation/test sets were representative of the overall distribution of the data. Then, nested K-fold cross-validation with a user-dependent approach was performed on the dataset, which consisted of a rotating K = 5 folds in the inner loop and K = 6 folds in the outer loop. The inner loops were used to stop the training when the accuracy of the training set increased, but the accuracy of the validation set decreased in case of overfitting of the model on the training set [41]. Moreover, they were also used to tune the hyperparameters of the classifier (e.g., number of layers, cells unit size, learning rate, etc.). Five trained models were created after K = 5 fold cross-validation. In each outer loop, the held-back test set was fed to each model (K = 5), and the majority vote decision was applied to select the final predicted class. The final performance of the classifier was averaged over 6 outer loops.

3.4. Data Analysis

Two regression models were trained to predict fatigue: random forest (RF) regression and convolutional neural network (CNN) regression. Random forest is more robust to hyperparameters compared to other traditional machine learning algorithms and is also computationally efficient [42]. In this study, the random forest regression hyperparameters were the number of estimators and the depth.
Compared with the commonly used recurrent neural network (RNN) for time series data, the convolutional layers of the CNN can better capture local and temporal patterns because the main assumption of the CNN model is that the same local patterns are relevant everywhere [43]. CNNs also tend to be more computationally efficient because there are fewer sequential calculations. Both random forest and CNN models are trained based on repetition based predictions.
The baseline CNN architecture has three 2D convolutional layers followed by a fully-connected layer (Figure 4). Rather than using hand-crafted features, the CNN directly learns features from the time series data without any prior knowledge of the features. Three convolutional layers were used to extract relevant features in the input signal. Each layer received the raw signals (or the one from the previous layer) and performed a convolution through a 2-dimensional kernel size of [10, 3] without pooling to extract relevant patterns. The filters for the three convolutional layers were 32, 64, and 128. The stride was [1, 1]. The rectified linear unit (ReLU) was employed as the activation function in the convolutional layers. To prevent overfitting, a dropout wrapper was added to each convolutional layer as a form of regularization to randomly select neurons (units) that were ignored at each epoch with a probability value of 0.8 [10]. Then, the output of the Conv layers was fed into a fully-connected layer composed of 128 units. Sigmoid was used as the activation function in the dense layer [10]. The Adam optimizer [44] was employed to minimize the cost function. The output layer had only one node. The CNN-based regression was trained with a batch size of 50 and a learning rate of η = 0.0001. The above hyperparameters were tuned by grid search. The model was trained by minimizing the loss function, where we used the root mean squared error (RMSE) between the predicted fatigue value and subjective fatigue level. The loss function (RMSE) is defined as:
R M S E = ∑ i = 1 N ( P r e d i c t e d i − A c t u a l i ) 2 N
Pearson’s correlation coefficient ρ between the estimated level of fatigue and self-reported fatigue level was used as the evaluation criterion (as shown in Equation (5)). The coefficient ρ can be categorized as weak ( ρ ≤ 0.35 ), moderate ( 0.35 < ρ ≤ 0.67 ), strong ( 0.67 < ρ ≤ 0.9 ), or excellent ( ρ > 0.9 ) [45,46]. For each K-fold cross-validation, the training phase was stopped when the accuracy of the validation set started to decrease, and the corresponding CNN-based regression was saved.
ρ = ∑ i = 1 n ( x i − x ¯ ) ( y i − y ¯ ) ∑ i = 1 n ( x i − x ¯ ) 2 ∑ i = 1 n ( y i − y ¯ ) 2
where x i is the predicted fatigue level, y i is the actual fatigue level, and x ¯ and y ¯ are the related mean values.

4. Results

The experiments were performed with open-source deep learning libraries, TensorFlow-based Keras [47] and Scikit-learn [48]. The computer used for the experiments was equipped with 3.6 GHz i7-8700 processors and an NVIDIA GTX1060 GPU.

4.1. Comparison of Different COP Features

Since random forest is very sensitive to the selection of features, different feature groups were trained in the random forest regression model with squat motion data. The results are summarized in Figure 5. The average Pearson coefficients were determined for all subjects, and the first COP feature set (i.e., the coordinates of COP and vertical ground reaction force) was found to perform best compared with other COP features. A one-way analysis of variance with repeated measures was performed to determine if the first COP feature set and the other COP feature sets differed significantly in their performance for predicting fatigue levels for each participant. The statistics software SPSS for Windows, Version 25 (IBM Corp., Armonk, NY, USA), was used in the analysis. It showed a significant main effect ( F ( 3 , 42 ) = 7.928 , p < 0.05 ) when the significance level was set at α = 0.05 , which showed that the first COP feature set performed significantly better than other COP feature sets at this 0.05 level.

4.2. Performance of Fatigue Estimation with Only IMU Data vs. Only FP Data

The selected COP and IMU feature sets were applied to train the random forest and CNN models separately for three different exercises. Then, the selected and trained models after cross-validation were fed with the corresponding test set. In Figure 6, the test results of each participant with different models are compared for the squat, high knee jack, and corkscrew toe-touch exercises, respectively. Examples of the regression relationships between the predicted and true fatigue levels are shown in Figure 7 (Participants 2, 5, and 8 based on force plate data for the squat exercise).
Paired t-tests with different models showed that there was no significant difference between the performances of the FP and IMU for the test sets at the 0.05 level in each random forest model and CNN model using the motion data from the squat, high knee jack, or corkscrew toe-touch exercises. For the squat and high knee jack motion exercises, the CNN model performed significantly better than the RF model, with both FP data and IMU data, while for the corkscrew toe-touch exercise, the performance of both methods did not differ significantly.
A Bland–Altman (B&A) plot was generated [49] to describe the difference between true fatigue values and predicted fatigue values. It was found that the majority of test points were within the 1.96 standard deviation (SD) agreement, and the mean difference between predicted values and true values was close to zero. In other words, most of the predicted fatigue values were close to the true values for each participant. For participants who achieved a strong Pearson coefficient, the agreement limits in the B&A analysis ranged from approximately −0.5 to 0.5; whereas for participants with a moderate or a weak Pearson coefficient, the corresponding limits of agreements were from around −1 to 1 and −2 to 2, respectively. Figure 8 shows example B&A plots of a participant with FP data with a moderate Pearson coefficient in squat motion and a strong Pearson coefficient in high knee jack and corkscrew toe-touch motions.

5. Discussion

The aim of the present work was to evaluate the methodologies for monitoring continuous fatigue levels automatically during exercise. The performances of the random forest and CNN classifiers using force plate and IMU data were evaluated in this study. By comparing the correlation coefficient in the RF model, the displacements of COP were found to be most highly associated with fatigue during exercise than other commonly used features of COP including standard variance, variance, skew, kurtosis, area of ellipse, and the first derivative and second derivative of COP positions.
By applying the selected feature set of the force plate and IMU in a participant-specific prediction model for continuous fatigue, the results showed that motion data from the FP had similar performance for fatigue prediction as the IMU in terms of the average Pearson correlation coefficient. Up to an 89%, 93%, and 94% Pearson coefficient for the squat, jack, and corkscrew exercises, respectively, could be achieved for a subject-dependent regression model for continuous fatigue prediction with motion data from the force plate. The proposed model translated the participants’ perceived exertion into numerical scores ranging from zero to 10, contributing to quantifying the movement performance via supervised learning. In addition, for most participants, the CNN performed better than the RF regression model. As expected, since the CNN is capable of extracting features automatically and no specific domain knowledge is required, the results were more accurate than those of random forest. Participants 2, 5, and 13 were found to consistently achieve strong correlation between the estimated fatigue level and their self-reported fatigue level in the squat exercise. The same conclusion can be drawn for Participants 5, 11, 13, and 14 in the high knee jack exercise and for Participants 2, 5, 11, and 13 in the corkscrew exercise. Considering the data and demographics (Table 3), the prediction performance was highly related to the number of sets that participants performed. The larger the number of repeated exercises per subject, the higher the performance of the prediction was, while the participants who performed a few reps of exercises had weaker correlations in the predicted model. This can be explained as large training data are expected to improve the model performance of the prediction. A further analysis with the same amount data for each participant will be conducted to examine the predictions in future work. Furthermore, when the experiments last for a long time, participants may lose interest and feel bored rather than really feel tired. In order to complete the experiments quicker, they may exaggerate fatigue levels, leading to unreliable results and low accuracy [50]. Therefore, it may be necessary to consider this respondent fatigue during complex and long surveys or experiments.
Unlike most previous work, which has been limited to studies on fatigue during one type of exercise [1,17,18], the present study was conducted on three different exercises. The test Pearson correlation coefficients in the RF and CNN models with FP data or IMU data demonstrated a good generalization ability of the proposed methods for different exercises. Since fatigue takes place gradually in a continuous fashion, the present work addressed the limitations associated with assuming only a two-class fatigue level prediction [17,18,19,20,29], thereby capturing the buildup of fatigue during the progress of exercise, which is important for preventing over-training at the early stage of fatigue, contributing to reducing the risk of injuries. Our model not only can predict the onset of fatigue, but also provide insight into the gradual change of exercise performance due to fatigue.
The most related to our study is the study [1], which conducted experiments to quantify the person-dependent continuous increase of fatigue over time, based on a squat dataset recorded with optical motion capture. It predicted fatigue for sets of squats assuming the foregoing squats were performed with the same fatigue level as the last squat within a set, which fit well with the subjective fatigue ratings since the ratings were given after sets of squats. In contrast, here, predicted fatigue was for every single repetition of squats without a prior assumption. In terms of accuracy, the predicted fatigue level in our study was also close to the real level based on the B&A analysis. The deviation of the regression line of the difference was small, which shows that the model can achieve good prediction whatever the level. In addition, optical motion capture was used in [1], which would pose challenges in real-world applications, since the capture space is limited and it is difficult to apply to outdoor activities. The proposed approach requires manual segmentation of the motion repetitions. For online deployment, automated segmentation needs to be implemented [51,52]. Alternatively, the approach could be modified to consider fixed-size windows, so that segmentation is not required [53]. Based on the successful findings of this study, there is potential to implement a real-time application for monitoring continuous fatigue-induced changes of motion based on data from wearable sensors or wireless insoles. In the future, the authors will explore the influence of sensor positions on movement detection across multiple exercises and a participant-independent approach for fatigue estimation by collecting a large-scale dataset through simulation techniques.

6. Conclusions

  • In this paper, a data-driven approach was presented for predicting fatigue based on data from IMUs or a force plate using random forest and the CNN.
  • A high correlation (Pearson coefficient of up to 89%, 93%, and 94% for the squat, jack, and corkscrew exercises, respectively) can be achieved for a subject-dependent regression model for continuous fatigue prediction.
  • The proposed model can detect fairly small variations of human motion due to fatigue, which also indicates the ability to monitor the gradual decline in the quality of the movement during exercise execution. The results contribute to the area of human motion recognition for more detailed knowledge of human motion.
  • Fatigue prediction using force plate motion data had a similar performance as the data from IMU. The displacement of COP was found to be correlated to fatigue compared with other commonly used features of the force plate. Bland–Altman analysis also confirmed that the predicted fatigue levels were close to the true values.
  • These results show that it is feasible to predict the onset of fatigue and quantify continuous fatigue levels from IMU data or force plate data through machine learning and deep learning.
  • Based on the successful findings with three different exercises, the general nature of the methodology is potentially applicable to a variety of other forms of exercise and developing real-time biofeedback applications based on kinematic or kinetic data, thereby contributing to future adaptation of exercise programs and prevention of over-training and injury as a result of excessive fatigue.

Author Contributions

Y.J. designed and carried out the study, performed the statistical analysis, and wrote the manuscript; V.H. participated in the design and data analysis and reviewed the manuscript; G.V. and D.K. conceptualized and designed the study, participated in the data collection and analysis, and reviewed the manuscript; B.K.C. participated in the data analysis and reviewed the manuscript. 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 according to the guidelines of the Declaration of Helsinki and was approved by the Ethics Committee of the University of Waterloo.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors thanks Jonathan Feng-Shun Lin for sharing the IIT and UW dataset.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Karg, M.; Venture, G.; Hoey, J.; Kulic, D. Human movement analysis as a measure for fatigue: A hidden Markov-based approach. IEEE Trans. Neural Syst. Rehabil. Eng. 2014, 22, 470–481. [Google Scholar] [CrossRef]
  2. Wan, J.j.; Qin, Z.; Wang, P.y.; Sun, Y.; Liu, X. Muscle fatigue: General understanding and treatment. Exp. Mol. Med. 2017, 49, e384. [Google Scholar] [CrossRef]
  3. Bogdanis, G.C. Effects of physical activity and inactivity on muscle fatigue. Front. Physiol. 2012, 3, 142. [Google Scholar] [CrossRef] [Green Version]
  4. Dong, H.; Ugalde, I.; Figueroa, N.; El Saddik, A. Towards whole body fatigue assessment of human movement: A fatigue-tracking system based on combined semg and accelerometer signals. Sensors 2014, 14, 2052–2070. [Google Scholar] [CrossRef] [Green Version]
  5. Jobanputra, C.; Bavishi, J.; Doshi, N. Human Activity Recognition: A Survey. Procedia Comput. Sci. 2019, 155, 698–703. [Google Scholar] [CrossRef]
  6. Bevilacqua, A.; MacDonald, K.; Rangarej, A.; Widjaya, V.; Caulfield, B.; Kechadi, T. Human activity recognition with convolutional neural networks. In Joint European Conference on Machine Learning and Knowledge Discovery in Databases; Springer: Berlin/Heidelberg, Germany, 2018; pp. 541–552. [Google Scholar]
  7. Wang, J.; Chen, Y.; Hao, S.; Peng, X.; Hu, L. Deep learning for sensor-based activity recognition: A survey. Pattern Recognit. Lett. 2019, 119, 3–11. [Google Scholar] [CrossRef] [Green Version]
  8. Hassan, M.M.; Huda, S.; Uddin, M.Z.; Almogren, A.; Alrubaian, M. Human activity recognition from body sensor data using deep learning. J. Med. Syst. 2018, 42, 99. [Google Scholar] [CrossRef] [PubMed]
  9. Sousa Lima, W.; Souto, E.; El-Khatib, K.; Jalali, R.; Gama, J. Human Activity Recognition Using Inertial Sensors in a Smartphone: An Overview. Sensors 2019, 19, 3213. [Google Scholar] [CrossRef] [Green Version]
  10. Vincent, H.; Dana, K.; Gentiane, V. Adversarial autoencoder for visualization and classification of human activity: Application to a low-cost commercial force plate. J. Biomech. 2020, 103, 109684. [Google Scholar]
  11. Buldt, A.K.; Forghany, S.L.; Landorf, K.B.; Murley, G.S.; Levinger, P.; Menz, H.B. Centre of pressure characteristics in normal, planus and cavus feet. J. Foot Ankle Res. 2018, 11, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cifrek, M.; Medved, V.; Tonković, S.; Ostojić, S. Surface EMG based muscle fatigue evaluation in biomechanics. Clin. Biomech. 2009, 24, 327–340. [Google Scholar] [CrossRef] [PubMed]
  13. Bell, T.; Li, B.; Zhang, S. Structured light techniques and applications. Wiley Encycl. Electr. Electron. Eng. 1999, 1–24. [Google Scholar] [CrossRef]
  14. Wang, L.; Hu, W.; Tan, T. Recent developments in human motion analysis. Pattern Recognit. 2003, 36, 585–601. [Google Scholar] [CrossRef]
  15. Fiore, L.; Fehr, D.; Bodor, R.; Drenner, A.; Somasundaram, G.; Papanikolopoulos, N. Multi-camera human activity monitoring. J. Intell. Robot. Syst. 2008, 52, 5–43. [Google Scholar] [CrossRef]
  16. Whittle, M.W. Chapter 2—Normal gait. In Gait Analysis, 4th ed.; Whittle, M.W., Ed.; Butterworth-Heinemann: Edinburgh, UK, 2007; pp. 47–100. [Google Scholar]
  17. Karvekar, S.; Abdollahi, M.; Rashedi, E. A Data-Driven Model to Identify Fatigue Level Based on the Motion Data from a Smartphone. In Proceedings of the 2019 IEEE Western New York Image and Signal Processing Workshop (WNYISPW), New York, NY, USA, 4 October 2019; pp. 1–5. [Google Scholar]
  18. Zhang, J.; Lockhart, T.E.; Soangra, R. Classifying lower extremity muscle fatigue during walking using machine learning and inertial sensors. Ann. Biomed. Eng. 2014, 42, 600–612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Buckley, C.; O’Reilly, M.A.; Whelan, D.; Farrell, A.V.; Clark, L.; Longo, V.; Gilchrist, M.; Caulfield, B. Binary classification of running fatigue using a single inertial measurement unit. In Proceedings of the 2017 IEEE 14th International Conference on Wearable and Implantable Body Sensor Networks (BSN), Eindhoven, The Netherlands, 9–12 May 2017; pp. 197–201. [Google Scholar]
  20. Strohrmann, C.; Harms, H.; Kappeler-Setz, C.; Troster, G. Monitoring kinematic changes with fatigue in running using body-worn sensors. IEEE Trans. Inf. Technol. Biomed. 2012, 16, 983–990. [Google Scholar] [CrossRef]
  21. Ramos, G.; Vaz, J.; Mendonça, G.; Pezarat-Correia, P.; Rodrigues, J.; Alfaras, M.; Gamboa, H. Fatigue Evaluation through Machine Learning and a Global Fatigue Descriptor. J. Healthc. Eng. 2020, 2020, 6484129. [Google Scholar] [CrossRef]
  22. Shair, E.; Ahmad, S.; Marhaban, M.; Mohd Tamrin, S.; Abdullah, A. EMG processing based measures of fatigue assessment during manual lifting. BioMed Res. Int. 2017, 2017, 3937254. [Google Scholar] [CrossRef]
  23. Wang, L.; Wang, Y.; Ma, A.; Ma, G.; Ye, Y.; Li, R.; Lu, T. A comparative study of EMG indices in muscle fatigue evaluation based on grey relational analysis during all-out cycling exercise. BioMed Res. Int. 2018, 2018, 9341215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Shaw, D.; Huang, C.R. Assessing muscle fatigue by measuring the EMG of biceps brachii. In Proceedings of the 2010 3rd International Conference on Biomedical Engineering and Informatics, Yantai, China, 16–18 October 2010; pp. 773–777. [Google Scholar]
  25. Chattopadhyay, R.; Pradhan, G.; Panchanathan, S. Subject independent computational framework for myoelectric signals. In Proceedings of the 2011 IEEE International Instrumentation and Measurement Technology Conference, Hangzhou, China, 10–12 May 2011; pp. 910–913. [Google Scholar]
  26. Al-Mulla, M.; Sepulveda, F.; Colley, M.; Kattan, A. Classification of localized muscle fatigue with genetic programming on sEMG during isometric contraction. In Proceedings of the 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Minneapolis, MN, USA, 3–6 September 2009; pp. 2633–2638. [Google Scholar]
  27. Pilarski, P.M.; Qi, L.; Ferguson-Pell, M.; Grange, S. Determining the time until muscle fatigue using temporally extended prediction learning. In Proceedings of the 18th International Functional Electrical Stimulation Society Conference (IFESS), Donostia-San Sebastian, Spain, 6–8 June 2013; pp. 7–8. [Google Scholar]
  28. Op De Beéck, T.; Meert, W.; Schütte, K.; Vanwanseele, B.; Davis, J. Fatigue prediction in outdoor runners via machine learning and sensor fusion. In Proceedings of the 24th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining, London, UK, 19–23 August 2018; pp. 606–615. [Google Scholar]
  29. Schmidt, M.; Rheinländer, C.C.; Wille, S.; Wehn, N.; Jaitner, T. IMU-based determination of fatigue during long sprint. In Proceedings of the 2016 ACM International Joint Conference on Pervasive and Ubiquitous Computing: Adjunct, Heidelberg, Germany, 12–16 September 2016; pp. 899–903. [Google Scholar]
  30. Liu, T.; Inoue, Y.; Shibata, K. A wearable ground reaction force sensor system and its application to the measurement of extrinsic gait variability. Sensors 2010, 10, 10240–10255. [Google Scholar] [CrossRef] [Green Version]
  31. Lee, C.H.; Sun, T.L. Evaluation of postural stability based on a force plate and inertial sensor during static balance measurements. J. Physiol. Anthropol. 2018, 37, 27. [Google Scholar] [CrossRef] [Green Version]
  32. Wu, P.P.Y.; Sterkenburg, N.; Everett, K.; Chapman, D.W.; White, N.; Mengersen, K. Predicting fatigue using countermovement jump force-time signatures: PCA can distinguish neuromuscular versus metabolic fatigue. PLoS ONE 2019, 14, e0219295. [Google Scholar] [CrossRef] [Green Version]
  33. Raymond, F.; Lussier, B.; Dugas, F.; Charbonneau, M.; Croteau, F.; Kennedy, C.; Berryman, N. Using Portable Force Plates to Assess Vertical Jump Performance: A Metrological Appraisal. Sports 2018, 6, 149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Adlerton, A.K.; Moritz, U.; Moe-Nilssen, R. Forceplate and accelerometer measures for evaluating the effect of muscle fatigue on postural control during one-legged stance. Physiother. Res. Int. 2003, 8, 187–199. [Google Scholar] [CrossRef] [PubMed]
  35. Baghdadi, A.; Megahed, F.M.; Esfahani, E.T.; Cavuoto, L.A. A machine learning approach to detect changes in gait parameters following a fatiguing occupational task. Ergonomics 2018, 61, 1116–1129. [Google Scholar] [CrossRef]
  36. Dos Santos, D.A.; Fukuchi, C.A.; Fukuchi, R.K.; Duarte, M. A data set with kinematic and ground reaction forces of human balance. PeerJ 2017, 5, e3626. [Google Scholar] [CrossRef] [Green Version]
  37. Nakazato, K.; Scheiber, P.; Müller, E. A comparison of ground reaction forces determined by portable force-plate and pressure-insole systems in alpine skiing. J. Sport. Sci. Med. 2011, 10, 754. [Google Scholar]
  38. Kassidas, A.; MacGregor, J.F.; Taylor, P.A. Synchronization of batch trajectories using dynamic time warping. AIChE J. 1998, 44, 864–875. [Google Scholar] [CrossRef]
  39. He, H.; Garcia, E.A. Learning from Imbalanced Data. IEEE Trans. Knowl. Data Eng. 2009, 21, 1263–1284. [Google Scholar]
  40. Butt, A.H.; Rovini, E.; Esposito, D.; Rossi, G.; Maremmani, C.; Cavallo, F. Biomechanical parameter assessment for classification of Parkinson’s disease on clinical scale. Int. J. Distrib. Sens. Netw. 2017, 13. [Google Scholar] [CrossRef] [Green Version]
  41. Hernandez, V.; Rezzoug, N.; Gorce, P.; Venture, G. Force feasible set prediction with artificial neural network and musculoskeletal model. Comput. Methods Biomech. Biomed. Eng. 2018, 21, 740–749. [Google Scholar] [CrossRef]
  42. Huang, B.F.; Boutros, P.C. The parameter sensitivity of random forests. BMC Bioinform. 2016, 17, 331. [Google Scholar] [CrossRef] [Green Version]
  43. Sadouk, L. CNN Approaches for Time Series Classification. In Time Series Analysis-Data, Methods, and Applications; IntechOpen: London, UK, 2018. [Google Scholar]
  44. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. arXiv 2015, arXiv:1412.6980. [Google Scholar]
  45. Fleron, M.K.; Ubbesen, N.C.H.; Battistella, F.; Dejtiar, D.L.; Oliveira, A.S. Accuracy between optical and inertial motion capture systems for assessing trunk speed during preferred gait and transition periods. Sport. Biomech. 2019, 18, 366–377. [Google Scholar] [CrossRef]
  46. Evans, J.D. Straightforward Statistics for the Behavioral Sciences.; Thomson Brooks/Cole Publishing Co.: Pacific Grove, CA, USA, 1996. [Google Scholar]
  47. Abadi, M.; Agarwal, A.; Barham, P.; Brevdo, E.; Chen, Z.; Citro, C.; Corrado, G.S.; Davis, A.; Dean, J.; Devin, M.; et al. TensorFlow: Large-Scale Machine Learning on Heterogeneous Systems. 2015. Available online: https://arxiv.org/pdf/1603.04467.pdf (accessed on 19 February 2020).
  48. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  49. Bland, J.M.; Altman, D. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986, 327, 307–310. [Google Scholar] [CrossRef]
  50. O’Reilly-Shah, V.N. Factors influencing healthcare provider respondent fatigue answering a globally administered in-app survey. PeerJ 2017, 5, e3785. [Google Scholar] [CrossRef] [Green Version]
  51. Lin, J.F.S.; Bonnet, V.; Panchea, A.M.; Ramdani, N.; Venture, G.; Kulić, D. Human motion segmentation using cost weights recovered from inverse optimal control. In Proceedings of the 2016 IEEE-RAS 16th International Conference on Humanoid Robots (Humanoids), Cancun, Mexico, 15–17 November 2016; pp. 1107–1113. [Google Scholar]
  52. Kim, M.; Park, S. Golf swing segmentation from a single IMU using machine learning. Sensors 2020, 20, 4466. [Google Scholar] [CrossRef]
  53. Um, T.T.; Babakeshizadeh, V.; Kulić, D. Exercise motion classification from large-scale wearable sensor data using convolutional neural networks. In Proceedings of the 2017 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Vancouver, BC, Canada, 24–28 September 2017; pp. 2385–2390. [Google Scholar]
Sensors 21 01499 g001 550
Figure 1. An overview of the proposed method. COP, center of pressure.
Figure 1. An overview of the proposed method. COP, center of pressure.
Sensors 21 01499 g001
Sensors 21 01499 g002 550
Figure 2. Distributions of the fatigue levels reported by participants in the experiment: (a) squat; (b) high knee jack; (c) corkscrew toe-touch.
Figure 2. Distributions of the fatigue levels reported by participants in the experiment: (a) squat; (b) high knee jack; (c) corkscrew toe-touch.
Sensors 21 01499 g002
Sensors 21 01499 g003 550
Figure 3. An illustration of the location of the markers. (a) Reflective markers for optical motion capture; (b) Xsens IMU sensors.
Figure 3. An illustration of the location of the markers. (a) Reflective markers for optical motion capture; (b) Xsens IMU sensors.
Sensors 21 01499 g003
Sensors 21 01499 g004 550
Figure 4. CNN structure (N represents the number of features. In this study, the value of N is 6 and 10 when using the force plate and IMU features as the input, respectively; 1, 32, 64, and 128 are the extracted channels at the corresponding convolution layer, and the stride size is [1, 1].).
Figure 4. CNN structure (N represents the number of features. In this study, the value of N is 6 and 10 when using the force plate and IMU features as the input, respectively; 1, 32, 64, and 128 are the extracted channels at the corresponding convolution layer, and the stride size is [1, 1].).
Sensors 21 01499 g004
Sensors 21 01499 g005 550
Figure 5. Correlation results with different COP features.
Figure 5. Correlation results with different COP features.
Sensors 21 01499 g005
Sensors 21 01499 g006 550
Figure 6. Results with IMU only and force plate (FP) only in the RF and CNN models: (a) squat; (b) high knee jack; (c) corkscrew toe-touch.
Figure 6. Results with IMU only and force plate (FP) only in the RF and CNN models: (a) squat; (b) high knee jack; (c) corkscrew toe-touch.
Sensors 21 01499 g006
Sensors 21 01499 g007 550
Figure 7. The predicted results of example subjects based on squat motion data from the FP representing the correlation between the predicted fatigue level and self-reported fatigue level (a) Subject 2 (strong correlation), (b) Subject 5 (moderate correlation), (c) Subject 8 (weak correlation).
Figure 7. The predicted results of example subjects based on squat motion data from the FP representing the correlation between the predicted fatigue level and self-reported fatigue level (a) Subject 2 (strong correlation), (b) Subject 5 (moderate correlation), (c) Subject 8 (weak correlation).
Sensors 21 01499 g007
Sensors 21 01499 g008 550
Figure 8. Example Bland–Altman (B&A) plot of a participant (subject = 5) with FP data for (a) squat, (b) high knee jack, and (c) corkscrew toe-touch. The blue line represents the mean differences. The red lines represent the 95% CI of the limits of agreement (the 1.96 SD). The pink line is the regression line of the differences. The orange dotted line is the line of equality (difference = 0).
Figure 8. Example Bland–Altman (B&A) plot of a participant (subject = 5) with FP data for (a) squat, (b) high knee jack, and (c) corkscrew toe-touch. The blue line represents the mean differences. The red lines represent the 95% CI of the limits of agreement (the 1.96 SD). The pink line is the regression line of the differences. The orange dotted line is the line of equality (difference = 0).
Sensors 21 01499 g008
Table
Table 1. Participants’ characteristics (SD: standard deviation).
Table 1. Participants’ characteristics (SD: standard deviation).
CharacteristicsRange
Number%MinMaxMean (SD)
Total Number14
Female2
Male12
Age (years) 203527.4 (4.15)
Anthropometric
Height (m) 1.621.871.76 (0.076)
Weight (kg) 579974.13 (12.13)
Exercise hours per week (h) 0124.53 (3.44)
0213.3
[1,5)653.8
[5,10)326.7
>10213.3
Table
Table 2. Exercise descriptions.
Table 2. Exercise descriptions.
NameInitial PoseDescription
SquatStandingBend knees to lower torso vertical position while keeping upper body upright and head forward
High knee jackStandingStart (arms above head; legs straight), ext1(right knee raise; arm comes down so hand claps under the leg), return to rest, ext2 (left knee raise; arm comes down so hand claps under leg)
Corkscrew toe-touchStandingStart (arm straight out to the side), ext1 (bend at waist, right hand touching left leg), return to rest, then ext2 (bend at waist, left hand touching right leg)
Table
Table 3. Number of recorded exercise sets.
Table 3. Number of recorded exercise sets.
Participant1234567891011121314
No. of sets of squat95231303126379201720153724
No. of sets of high knee jack6101514441331773221112215
No. of sets of corkscrew toe-touch420101631121864321593615
Table
Table 4. Explanation of generated features.
Table 4. Explanation of generated features.
NumberNameExplanation
1COPThe displacements of the center of pressure of two feet during each repetition
2Variance of COPDeviation of COP velocities in anteroposterior (AP) and mediolateral (ML) directions during each repetition
3STD of COPStandard deviation of COP velocities in AP and ML directions during each repetition
4Mean velocityMean velocity, which represents the total distance traveled by COP for each repetition in AP and ML directions
5First derivative of COPThe change rates of COP displacements (i.e., velocity) for each repetition in AP and ML directions
6Second derivative of COPThe change rates of COP velocities (i.e., acceleration) for each repetition in AP and ML directions
7Skew of COPThe value of skewness of COP displacements for each repetition in AP and ML directions
8Kurtosis of COPThe value of kurtosis of COP displacements for each repetition in AP and ML directions
9Ellipse of COPThe area of the ellipse covered by the trajectory of the COP for each repetition with a 95% confidence interval
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, Y.; Hernandez, V.; Venture, G.; Kulić, D.; K. Chen, B. A Data-Driven Approach to Predict Fatigue in Exercise Based on Motion Data from Wearable Sensors or Force Plate. Sensors 2021, 21, 1499. https://doi.org/10.3390/s21041499

AMA Style

Jiang Y, Hernandez V, Venture G, Kulić D, K. Chen B. A Data-Driven Approach to Predict Fatigue in Exercise Based on Motion Data from Wearable Sensors or Force Plate. Sensors. 2021; 21(4):1499. https://doi.org/10.3390/s21041499

Chicago/Turabian Style

Jiang, Yanran, Vincent Hernandez, Gentiane Venture, Dana Kulić, and Bernard K. Chen. 2021. "A Data-Driven Approach to Predict Fatigue in Exercise Based on Motion Data from Wearable Sensors or Force Plate" Sensors 21, no. 4: 1499. https://doi.org/10.3390/s21041499

APA Style

Jiang, Y., Hernandez, V., Venture, G., Kulić, D., & K. Chen, B. (2021). A Data-Driven Approach to Predict Fatigue in Exercise Based on Motion Data from Wearable Sensors or Force Plate. Sensors, 21(4), 1499. https://doi.org/10.3390/s21041499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop