EMG Processing Based Measures of Fatigue Assessment during Manual Lifting

Abstract

Manual lifting is one of the common practices used in the industries to transport or move objects to a desired place. Nowadays, even though mechanized equipment is widely available, manual lifting is still considered as an essential way to perform material handling task. Improper lifting strategies may contribute to musculoskeletal disorders (MSDs), where overexertion contributes as the highest factor. To overcome this problem, electromyography (EMG) signal is used to monitor the workers' muscle condition and to find maximum lifting load, lifting height and number of repetitions that the workers are able to handle before experiencing fatigue to avoid overexertion. Past researchers have introduced several EMG processing techniques and different EMG features that represent fatigue indices in time, frequency, and time-frequency domain. The impact of EMG processing based measures in fatigue assessment during manual lifting are reviewed in this paper. It is believed that this paper will greatly benefit researchers who need a bird's eye view of the biosignal processing which are currently available, thus determining the best possible techniques for lifting applications.

E. F. Shair 1; 2 and S. A. Ahmad 1 and M. H. Marhaban 1 and S. B. Mohd Tamrin 3 and A. R. Abdullah 2

Academic Editor:Stefan Rampp

1, Department of Electrical and Electronic Engineering, Faculty of Engineering, Universiti Putra Malaysia (UPM), Selangor, Malaysia

2, Faculty of Electrical Engineering, Universiti Teknikal Malaysia Melaka (UTeM), Melaka, Malaysia

3, Department of Environmental and Occupational Health, Faculty of Medicine and Health Sciences, Universiti Putra Malaysia (UPM), Selangor, Malaysia

Received 28 November 2016; Accepted 31 January 2017; 19 February 2017

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Musculoskeletal disorders (MSDs) caused by manual lifting tasks have been perceived for quite some time as one of the primary work-related injuries which influences the personal satisfaction of industrial workers around the world [1, 2]. Common reasons of MSDs in manual lifting are due to improper lifting and muscle fatigue. National Institute of Occupational Safety and Health (NIOSH) has published a technical report entitled Work Practices Guide for Manual Lifting as a reference to perform correct lifting [3]. However, in terms of muscle fatigue monitoring, researchers are still in search of finding the best processing technique to assess the human muscle condition.

2. Manual Lifting

Despite the wide utilization of robots in modern industries, there are still numerous assignments in the industry that are performed physically by humans, such as carrying, lifting, pushing, or pulling. These physical movements of the body often cause MSDs that can be subdivided into more particular and perceived regions, such as back disorders (BDs), upper limb disorders (ULDs), and lower limb disorders (LLDs) [4]. ULDs cover distinctive work-related musculoskeletal dissensions in the arm, wrist, hand, elbow, neck, and shoulder, while LLDs are associated with the lower areas of the body such as feet, legs, and hips. The anatomical zones of the body mostly influenced by these MSDs issues due to manual handling are in the spine and lumbar (back area), and the arm, and wrist (upper limb area). These are further supported by various researches done in previous years on muscle areas investigated during manual lifting and are summarized in Table 1.

Table 1: Common muscles investigated during manual lifting.

Reference	Muscle areas
Kamarudin et al. [7]	Biceps brachii and triceps brachii
Voge and Dingwell [8]	Biceps, middeltoid, midtrapezius, and postdeltoid
Zawawi et al. [9]	Biceps brachii
Al-Ashaik et al. [10]	Biceps brachii, erector spinae, anterior deltoid, and trapezius
Granata and Marras [11]	Erector spinae, external abdominal obliques, internal abdominal obliques, rectus abdominis, and latissimus dorsi
Roy et al. [12]	Longissimus thoracis, iliocostalis lumborum, and multifidus
Seroussi and Pope [13]	Erector spinae and external oblique
Graham et al. [14]	Lumbar erector spinae, thoracic erector spinae, and rectus abdominis
Gagnon et al. [15]	Thoracic erector spinae, lumbar erector spinae, internal oblique, external oblique, rectus abdominis, and latissimus dorsi
Dolan and Adams [16]	Erector spinae
Bonato et al. [17]	Iliocostalis lumborum, longissimus thoracis, and multifidus
Potvin [18]	Thoracic erector spinae and lumbar erector spinae
Kingma and Van Dieën [19]	External oblique, internal oblique, rectus abdominis, iliocostalis lumborum, longissimus thoracis, and pars thoracis
Shin and Kim [20]	Erector spinae, rectus abdominis, latissimus dorsi, internal oblique, and external oblique
Cholewicki et al. [21]	Lumbar erector spinae, thoracic erector spinae, latissimus dorsi, internal oblique, external oblique, and rectus abdominis

Manual errands that are performed improperly by workers in the industries may most of the time bring pain in the lumbar spinae [5]. This situation happens when one individual performs forward flexion; the weight of abdominal area and the load affects the lumbar spinae, which is adjusted by trunk muscles (predominately by erector spinae muscles) [6]. This movement causes the erector spinae muscles to make incredible expansive strengths for adjusting the body. Even though there are numerous assistive gadgets such as lift trucks, trolleys, and derricks that can help reduce loads on the lumbar spinae, however in many cases, these instruments are costly and inconvenient.

Lifting strategies are an important measure during lifting and because physically correct lifting is significant, their impact on musculoskeletal well-being is broadly talked about [22]. The impact of lifting techniques including position width (characterized as the separation between the feet in the average sidelong heading with the sagittal symmetry of position) has been considered, recommending the utilization of a wide position to lessening the burden on the spine [23]. Lifting procedures (stoop or squat and semisquat) can also have a significant effect on spine stacking and steadiness during lifting. Past studies have researched lifting methods in biomechanical terms to give some mediation methodologies and recognize the correct lifting procedures [24, 25]. Squat lifting seems to have less lumbar shear stretch and place less weight on noncontractile connective tissues. Stoop lifting has all the earmarks of being more regular and less exhausting. Semisquat lifting might be a decent trade-off lift, despite the fact that early confirmation recommends perhaps higher lumbar minutes [22].

3. Muscle Fatigue Assessment

Muscle fatigue is characterized as the long lasting deterioration of the performance of the human operator to create force [24, 25]. A person's muscle capabilities may differ from one another, but the characteristics should have some resemblance. Past researchers have outlined several features of muscle fatigue such as an increment in amplitude and the transition from high frequency range to low frequency [26].

To assess this muscle fatigue characteristics, various techniques are available, depending on the usage. At this moment, even though the measurement of muscle fatigue by use of invasive means such as blood test (blood lactate level and blood oxygen level) and muscle biopsies (pH of muscle) offers higher accuracy compared to noninvasive techniques, it is unsuitable for some applications such as sports, ergonomics, and occupational therapy [27]. The use of noninvasive approach to acquire muscle signals such as electromyography (EMG), sonomyography (SMG), mechanomyography (MMG), and near infrared spectroscopy (NIRS) for research purposes is considered to be more convenient and easy to use [28]. Every technique has its own advantages and disadvantages based on its applications. For the purpose of assessing muscle physiology during manual lifting, EMG is the most commonly used tool and has been well reported by many researchers [29].

4. Electromyography Signal

EMG is the measure of electrical potential present on the skin in consequence to a muscle contraction that represents the neuromuscular activities [26]. It can be measured by two methods: (1) applying electrodes to the skin surface (noninvasive) or (2) intramuscular (invasive) within the muscle. Even though intramuscular EMG (imEMG) provides additional benefits to overcome the drawbacks of surface EMG (sEMG) such as maintaining robust electrode contact with the skin and the capacity to record from profound muscles with little EMG crosstalk [27], sEMG and imEMG have been proven to have equal classification performance for data of comparative nature (unmodulated). However the performance of imEMG diminished, compared to surface, when tested on modulated data [28, 29].

The greater selectivity of imEMG with respect to sEMG is due to wires exposed only at the tip but then again may be a hindrance, since the signal may provide local, rather than global information where imEMG recordings depend on the recruitment of motor units; it may be the case that insufficient information was captured at low amplitude/frequency [27].

5. Electromyography Preprocessing

When dealing with EMG signal, noise removal is a very important factor to be taken into consideration. This is due to the sensitive nature of EMG signal itself that different types of noises or artifacts unavoidably contaminate it. These noises are caused by various sources originates from either the skin-electrode interface, hardware that intensifies the signal, or other external sources. Therefore, before proceeding with the analysis and classification of EMG signal, a preprocessing is usually introduced to remove these noises. Some of the noises that may affect the signals are inherent noise in electronics equipment, ambient noise, motion artifact, ECG artifacts, inherent instability of signal, and crosstalk. However, this review only covers three types of noises: motion artifact, ECG noise, and device noise.

5.1. Motion Artifact

Cutting edge technology is considerably invulnerable to some of the EMG noises, but not to the motion artifact as it has frequency spectra that are affecting the low frequency part of EMG signal. This motion artifact originates at the electrode-skin interface, where there exist muscle movements underneath the skin, and force impulse that goes through the muscle and the skin underlying the electrodes causing movements [30]. These movements will result in a time-varying voltage across the two electrodes.

To remove the low frequency noise caused by the motion artifact, a study has been conducted to ascertain a reasonable value for the high-pass corner frequency filter. The result of the study recommends the use of a Butterworth filter with a 20 Hz corner frequency and a 12 dB/oct slope [31].

5.2. ECG Noise

When EMG signal is measure close to the heart, ECG bursts may pollute the EMG recording, which can affect the analysis and result in misinterpretations [32]. This is a natural relic that is often unavoidable. It can be diminished by proper skin preparation and adjusted position of the ground electrode. To date, other than the conventional ECG removal procedures such as the use of high past filtering, several ECG noise removal techniques have been introduced that can clean these ECG bursts but still maintaining the regular EMG characteristics.

In 2009, Lu et al. came out with a study to evaluate the performance of recursive-least-square adaptive filter in the removal of ECG interference from EMG signal [33]. The result of the study shows that due to fast convergence of the recursive-least-square algorithm, the filter is reliable to effectively remove ECG noise.

Later in 2012, independent component analysis (ICA) was proposed to remove the ECG contamination. The filtering effects are measured in terms of root mean square relative errors and linear envelopes correlation of uncontaminated and contaminated EMG. This technique excellently produces good results when EMG and ECG are statistically independent [34].

5.3. Device Noise

Another type of noise that usually affects EMG signal is the noise generated inside the EMG device itself that is commonly referred to as inherent noise in electronics equipment [35]. This noise cannot be completely eliminated as all electronic equipment generates noise [36]. However, the severity of the noise can be reduced using high-quality electronic components.

6. Electromyography Processing

Biosignal processing is a critical part in biomedical engineering in order to classify the frequency content of a signal. These signals include EMG, electrocardiography (ECG), and electroencephalography (EEG). Most of these biosignals are nonstationary signals due to its time-varying characteristics [37-39]. The problem of this nonstationary signal is the process of assuming it to be stationary over short-time intervals, where stationary analysis techniques are used. However, this assumption is not always suitable, and further methods for nonstationary processes are needed.

Past researchers have developed several techniques to process the nonstationary signal specifically the EMG signal and these include parametric and nonparametric approaches [40, 41]. The differences between these two approaches are mainly due to the parameters involved. Since most of the signals encounters are unknown signals without known parameters, the development of nonparametric approach, which involves a variety of transforms, is actively researched compared to the parametric approach that requires the modeling of the nonstationary process. Figure 1 provides an overview of the EMG signal processing methods that are reviewed in the manuscript.

Figure 1: Overall view of the electromyography signal processing methods that are reviewed in this manuscript.

[figure omitted; refer to PDF]

6.1. Parametric Approach

In the field of EMG processing, one type of parametric approach that has been explored is based on the time-varying linear predictive models. This includes time-varying autoregressive (TVAR) model where its parameters vary with time. The time-varying spectrum in TVAR is estimated from the time-varying model parameters and the nonstationary signal's instantaneous frequency (IF) can be extracted [42]. IF represents the spectral peak location of the EMG signal as it varies with time [43]. Even though the poles and zeros of the estimated model can be estimated directly from the EMG signal, it is not guaranteed that the time-varying poles remain inside the unit circle on the z-plane, thus making TVAR model temporarily unstable [44]. The constraints can be easily imposed by factorizing the denominator of TVAR in direct form and formulating it in cascaded form [45]. Computation of the mean frequency (MNF) is taken as the time average of the IF estimation based on the cascaded model in each interval [46].

Research conducted by Al Zaman et al. has compared TVAR with short-time Fourier transform (STFT) and concludes that even though both TVAR and STFT have produced similar decreasing muscle fatigue patterns from the MNF, TVAR has achieved a better performance with the slope and interception of the linear regression line is closer [47].

Another parametric approach that has been used is the time-varying autoregressive moving average (TVARMA). Studies have been made to compare the performance of TVAR and TVARMA in analyzing EMG signals, and it is concluded that AR models are more favorable than ARMA because of the following reasons: (a) estimation of AR models is moderately basic since they are represented by an arrangement of direct mathematical statements; (b) AR models are predominant in displaying the low recurrence sEMG segment; and (c) fluctuation of the AR model estimation on shorter fragments is lower than for ARMA models [40].

6.2. Nonparametric Approach

Another EMG processing method is the nonparametric approaches, which include STFT, spectrogram, and wavelet transform. These approaches appear to be more favorable to researchers since the fatigues indices can be obtained directly from the EMG signal without knowing its parameters. The analysis of the EMG signal based on the nonparametric approaches can be further divided into three types of distribution techniques: time domain, frequency domain, and time-frequency domain [48-50].

6.2.1. Time Distribution (TD)

Extraction of time domain fatigue indices is very simple and easy and involves low computational complexity compared to the other two techniques since it does not involve any transformation [51]. It can be directly obtained from the time representation of the raw EMG signal just by doing some simple mathematical statistics. The idea of fatigue is approximately connected with the amount and rate of change of some variables that reflect muscle alterations amid sustained contractions [52].

Listed in Table 2 are several fatigue indices in the time domain that have been explored by previous researchers. Even though the analysis based on time domain is already established and are widely used since the last decade, however it is considered less accurate since their calculations are based on EMG signal amplitude values, where there are much interference acquired through the recording, especially for features that are extracted from energy property [53]. Another major disadvantage of this time domain features that affects the accuracy comes from the nonstationary property of the EMG signal where the statistical properties changes over time, but time domain features assume the data as stationary signal [54].

Table 2: Fatigue indices in time domain.

Reference	Fatigue indices
Malinzak et al. [55]	Integrated EMG (IEMG)
Arabadzhiev et al. [56]	Root mean square (RMS)
Suetta et al. [57]	Mean absolute value (MAV)
Oskoei and Hu [58]	Modified mean absolute value type 1 (MAV1)
Villarejo et al. [59]	Modified mean absolute value type 2 (MAV2)
Phinyomark et al. [60]	Simple integral square (SSI)
Zardoshti-Kermani et al. [61]	Variance of EMG (VAR)
Tkach et al. [62]	v-order (V)
Zardoshti-Kermani et al. [61]	Log detector (LOG)
Kiguchi et al. [63]	Waveform length (WL)
Al Omari et al. [64]	Average amplitude change (AAC)
Siddiqi et al. [65]	Difference absolute standard deviation value (DASDV)
Kilbom et al. [66]	Zero crossing (ZC)
AlOmari and Liu [67]	Myopulse percentage rate (MYOP)
Phinyomark et al. [68]	Willison amplitude (WAMP)
Rogers et al. [69]	Slope sign change (SSC)
Buchenrieder [70]	Mean absolute value slope (MAVSLP)
Venugopal et al. [71]	Multiple hamming windows (MHW)
Du and Vuskovic [51]	Multiple trapezoidal windows (MTW)
Phinyomark et al. [72]	Histogram of EMG (HIST)
Al-Quraishi et al. [73]	Autoregressive coefficient (AR)
Chang et al. [74]	Cepstral coefficient (CC)

Nowadays, researches involving time distribution only focus on finding new fatigue indices with higher robustness, while most of other researchers nowadays opt for time-frequency analysis, which is considered to be more accurate.

6.2.2. Frequency Distribution (FD)

In many studies of fatigue muscle, frequency domain features are usually used to extract information from a signal. To date, there are more than 20 fatigue indices in the frequency domain, where two of the widely used parameters are the mean frequency (MNF) and median frequency (MDF) [75]. Other fatigue indices are listed in Table 3. Several techniques have been used to extract information from EMG signals such as fast Fourier transform (FFT), power spectral density (PSD), and parametric methods (i.e., AR model).

Table 3: Fatigue indices in frequency domain.

Reference	Fatigue indices
Merletti and Lo Conte [76]	Mean frequency (MNF)
De Luca [77]	Median frequency (MDF)
Khanam and Ahmad [78]	Peak frequency (PKF)
Khanam and Ahmad [78]	Mean power (MNP)
Khanam and Ahmad [78]	Total power (TTP)
Altaf et al. [79]	The 1st, 2nd, and 3rd spectral moments (SM1, SM2, and SM3)
Phinyomark et al. [80]	Power spectrum ratio (PSR)
González-Izal et al. [81]	Instantaneous frequency variance (Fvar )
Georgakis et al. [75]	Averaged instantaneous frequency (AIF)
Dimitrov et al. [82]	Dimitrov's spectral fatigue index (FInsm5 )

The EMG signal should undergo Fourier transform in order to represent the signal in frequency domain. The definition of Fourier transform in the form of power spectrum is shown in [figure omitted; refer to PDF] where x(t) is the signal in time domain.

However, the fluctuations of the EMG frequency component due to the adjustments in the muscle force, length and contraction speed throughout time, have caused challenges in the usage of FFT and other traditional processing methods. The fundamental confinement of a FFT is that it cannot give simultaneous time and frequency localization [83]. Therefore, investigation of the EMG signal in dynamic contraction (i.e., repetitive lifting) utilizing these methods may not be successful since it requires the signal to be stationary.

6.2.3. Time-Frequency Distribution (TFD)

Time-frequency representation (TFR) of a signal maps a one-dimensional signal of time into a two-dimensional of time and frequency. In analyzing, modifying, and synthesizing nonstationary signals, TFR are widely used since both representation of time and frequency are taken into account, thus leading to higher accuracy [84].

In general, TFDs consist of linear TFDs and bilinear (quadratic) TFDs. The most basic form of TFD technique is the short-time Fourier transform (STFT). STFT is one of the linear TFDs that have been used to analyze EMG signals along with spectrogram, wavelet transform, S-transform, and so forth. Examples of the bilinear TFDs include the Wigner-Ville distribution (WVD) and Choi-William distribution (CWD).

Linear TFD

Short-Time Fourier Transform (STFT) . To overcome the disadvantages of time distribution and frequency distribution and satisfy the stationary condition, it is common to separate long haul signals into blocks of narrow fragments [85]. This ought to be sufficiently narrow to be viewed as stationary and take the FT of every segment. Every FT gives the spectral information of a different time-slice of the signal, giving simultaneous time and frequency estimation.

This was proposed by Gabor who had developed the STFT, the extended version of FT [86]. [figure omitted; refer to PDF] where x(τ) is the EMG signal, w(τ-t) is the observation window, and t is the variable that slides the window over the signal, x(τ).

The choice of the window function in STFT is critical in order to get accurate results. The shape of the window can be either rectangular, Gaussian, or elliptic, depending on the shape of the signal. For sampled STFT using a Gaussian window, it is often called Gabor transform. Even though the window should be narrow enough to ensure that the portion of the signal that falls within the window is stationary, it must not be too narrow since it will lead to bad localization in the frequency domain. For infinitely long window, w(t)=1 will eventually cause STFT to turn into FT, providing excellent frequency localization, but no time localization. In contrast to infinitely short window, w(t)=δ will result in the time signal (with a phase factor), which will provide excellent time localization but no frequency localization.

Research by MacIsaac et al. uses MNF as fatigue index and has drawn three important conclusions. First, the nonstationarities in EMG do not seem to affect the MNF values obtained from Fourier transform. Secondly, the STFT method to measure MNF has successfully midpoints out the impacts of nonstationarities in dynamic compressions. Third, the STFT is capable of identifying a descending pattern with fatigue in dynamic contraction, as long as the range of motion remains consistent across sequential time interval [87].

Since the STFT is basic in idea and execution and is equipped for distinguishing a pattern in muscle fatigue, it is an undeniable choice for applications that require simple analysis with acceptable results.

Spectrogram . The squared magnitude of the STFT is called spectrogram. It can be expressed as [figure omitted; refer to PDF] where x(τ) is the EMG signal, w(τ-t) is the observation window, and t is the variable that slides the window over the signal, x(τ).

Spectrogram can be used to obtain the power distribution and energy distribution of the signal along the frequency direction at a given time. In EMG processing, the instantaneous energy is useful to separate muscle activation from the baseline, which is called the segmentation process. Segmentation is important in reducing the high computational complexity of TFDs. The muscle activation segmentation process by using spectrogram has been proposed by Shair et al. and resulted in a mean absolute percentage error (MAPE) of just 1.404% [88].

For both STFT and spectrogram, there is a compromise between the time-based and frequency-based perspective of a signal. Both time and frequency are represented in limited precision, where precision is controlled by the span of the window, and the size of the window chosen will be the same for all frequencies.

In 2010, Jankovic and Popovic have presented in their paper the use of spectrogram in the analysis of muscle fatigue. They highlighted that even though the traditional MDF technique provides good sensitivity for fatigue indication, however, there is a significant slope error particularly for low activity of coinitiated muscles. They proposed a hybrid of alternative methods (spectrogram and scalogram) that can successfully provide complete information but at the same time produce less error prediction for low-level coinitiated muscles [89]. As a recommendation for improvement, they also proposed further investigation on the robustness of this method in dynamic exercise contractions.

Later in 2016, Zawawi et al. came out with a detailed analysis of EMG signals by using spectrogram for manual lifting application. By taking the average instantaneous RMS voltage (_Vrms (t)) as fatigue index, she concludes that, as the lifting height is increased, the average _Vrms (t) is also increased. However, as the average _Vrms (t) decreased, the number of lifting repetitions is increased [90].

Wavelet Transform (WT) . A wavelet is a waveform of effectively limited duration that has an average value of zero. Some of its properties are short-time localized waves with zero integral value, the possibility of time shifting, and flexibility. Wavelet analysis produces a time-scale view of the signal. The result of a continuous wavelet transform (CWT) is wavelet coefficients. Multiplying each coefficient with the appropriately scaled and shifted wavelet yields the constituent wavelet of the original signal. Equation (4) shows the formula for CWT. [figure omitted; refer to PDF] where t is the translation, a is the scale parameter, and Ψ is the mother wavelet.

Yochum et al. have proposed a new fatigue index based on the continuous wavelet transform (CWT), named _ICWT , and have compared it to other fatigue indices from literature in terms of their sensitivity to noise. The results show that this new fatigue index quantifies the EMG signal elongation during a contraction and thus makes it a suitable fatigue index. _ICWT is less subjected to noise and less truncation dependent compared to other fatigue indices based on the frequency like MNF and MDF [91].

A comparison has been made by Dantas et al. between STFT and CWT in evaluating muscle fatigue in isometric and dynamic contractions. The after effects of this study exhibit that CWT and STFT analysis give comparable fatigue estimates (slope of MDF) in isometric and dynamic contractions. However, the aftereffects of CWT for both contractions indicate less variability (higher accuracy) in EMG signal analysis contrasted with STFT [92].

S-Transform . Another TFD technique exists is the S-transform. This technique has been widely used in diverse areas of telecommunication, power quality, geophysics, and biomedicine, due to its obvious advantages in processing nonstationary signals [93]. In the area of EEG and ECG, despite the use of S-transform in the time-frequency analysis, this technique is also a good alternative for denoising and removal of artifacts that exist in ECG and EEG signals [94]. Even though it has been explored in the areas of EEG and ECG, there are no available researches conducted in utilizing this technique for EMG signal.

S-transform was introduced in 1996 by Stockwell, Mansinha, and Lowe. It is a hybrid of two advanced signal processing techniques which are STFT and WT [95]. Due to this, it inherits the good qualities from both techniques. It provides good resolution in both time and frequency and allows users to assess any frequency component in the time-frequency domain without the need of using any digital filter. Even though S-transform uses a variable window length to maintain a good time-frequency resolutions for all frequencies, it still retains the phase information using a Fourier kernel.

Expression for S-transform is shown in [figure omitted; refer to PDF] where τ and f denote the time of the spectral localization and Fourier frequency, respectively.

The development of the S-transform is led by the urge to overcome the low resolution of STFT and the absence of the phase information in CWT. Since the features had existed in the S-transform, it is sometimes viewed as a variable sliding window STFT or a phase-corrected CWT [96].

Although S-transform has better time-frequency resolution than STFT, the resolution is still far from perfect and needs improvement. To date, there are several improved S-transform introduced by researchers such as the modified S-transform [96] and discrete orthonormal S-transform [94].

Bilinear TFD

Wigner-Ville Distribution (WVD) . Wigner first introduced WVD in the area of quantum mechanics in the year 1932, and later in 1948, Ville developed and applied the same transformation to signal processing and spectral analysis [97].

Compared to STFT and WT, WVD does not contain a windowing function and thus frees WVD from the smearing effect due to the windowing function. As a result, it provides best possible temporal and frequency resolution in the time-frequency plane. It possesses several interesting properties, such as energy conservation, real-valued, marginal properties, translation covariance, dilation covariance, instantaneous frequency, and group delay.

However, because of its quadratic nature, when dealing with signals with several frequency components, WVD suffers from the so-called cross components (inference terms), which represent significant defects in this method.

These interference terms are troublesome since they may overlap with autoterms (signal terms) and in this manner, it is hard to visually interpret the WVD image. It creates the impression that these terms must be available or the good properties of WVD (localization, group delay, instantaneous frequency, marginal properties, etc.) cannot be fulfilled. There is also a trade-off between the amount of interferences and the number of good properties. These are known issues with the WVD spectrum, and there are several ways to compensate these cross terms that has been discussed by previous researchers [98-100]. Equation (6) shows the formula for WVD. [figure omitted; refer to PDF] where x(t+τ/2)^{x[low *]} (τ-τ/2) is the instantaneous autocorrelation function and [low *] shows conjugate operation.

In the occasion that time smoothing window g(t) and a frequency smoothing window h(t) are connected to WVD, WVD will then transform into the smoothed-pseudo-Wigner-Ville distribution (SPWVD), as composed in [figure omitted; refer to PDF] where W is the WVD.

A research by Subasi and Kiymik in 2010 has compared the performance between STFT, SPWVD, and CWT in terms of accuracy, specificity, and sensitivity. The results suggested that all three methods provide almost similar performance and are found to be satisfactory, despite the fact that there are little differences between the results [99].

Choi-William Distribution (CWD) . CWD is a member of Cohen's class distribution function and was proposed in the year 1989 by [101]. It is able to avoid one of the main problems of WVD, which is the presence of interference in regions where one would expect zero power values. Adoption of the exponential kernel in the distribution helps to suppress the effect of cross term [102]. CWD is given by [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and the kernel function is given by [figure omitted; refer to PDF]

By studying the influence of muscle contraction towards the frequency content of EMG signal using WVD and CWD, Alemu et al. observed that cross terms existing in WVD are greatly reduced by CWD, thus resulting in easier interpretation with no loss of definition due to cross terms [103]. The reduction of the cross terms is due to the smoothing parameters introduced in CWD. However, this still depends on the smoothing parameter value selected. As observed from the result itself, when the smoothing parameter decreased, the reduction of the cross terms is better. This, however, affects the time and frequency resolutions by reducing it and leads to more signal loss. As the smoothing parameter approaches ∞, it resembles WVD more [103].

B-Distribution . Karthick et al. introduced B-distribution in 2015 to attenuate the cross terms exist in WVD and CWD. This is realized by introducing a smoothing function referred to as quadratic time-frequency distribution (QTFD) [104]. Expression of the QTFD is given as follows [105]: [figure omitted; refer to PDF] where G[n,m] is the smoothing kernel and operator [low *] represents convolution.

Kernel of the B-distribution is given by [figure omitted; refer to PDF] where β is a smoothing parameter.

Three fatigue indices, namely, instantaneous median frequency (IMDF), instantaneous mean frequency (IMNF), and instantaneous spectral entropy (ISPEn), were used alongside the B-distribution and it is found that all three features are distinct in both fatigue and nonfatigue conditions. The study concludes that B-distribution may be useful in EMG analysis under various clinical and normal conditions.

Modified B-Distribution . A year later, the same author came out with a study using modified B-distribution to monitor the progression of muscle fatigue. Since the performance of time-frequency distribution is assessed based on its ability to reduce cross term and provide closely spaced frequency components representation, this modified version successfully removes cross term interference, while maintaining high time-frequency resolution [106]. It is demonstrated that modified B-distribution based time-frequency distribution performs better in suppressing cross terms with good time and frequency resolution for multicomponent signals compared with other above-mentioned techniques [107]. Recently, this technique has been widely used to analyze nonstationarities related to EEG signals, heart rate signals, and accelerometer data based on fetal movements [108, 109].

7. Discussion and Conclusion

The use of EMG processing to assess muscle fatigue during manual lifting was reviewed and several fatigue indices were introduced. Research in the area of EMG signal currently evolves around the sensitivity, variability, and repeatability of the fatigue indices, as well as the best processing techniques with high efficiency and less computational complexity. Despite the use of traditional methods such as time distribution and frequency distribution, time-frequency distribution is found to be more superior in terms of monitoring since it enables the user to observe the progress of the signal in time and frequency. This is very important as at the point when muscle fatigue sets in, frequency compression can happen. Unlike the conventional frequency distribution technique, time-frequency analysis demonstrates the time when any changes in frequencies occur.

In recent years, researchers had identified that the bilinear TFD could perform better than the linear TFD such as STFT, spectrogram, and WT due to the fact that it does not suffer from the smearing effects cause by windowing function. Nonetheless, on account of its quadratic nature, bilinear TFD suffers the cross term effects. With the existence of high resolution Cohen class TFD such as the modified B-distribution introduced in 2016, the cross term effects can be removed and at the same time maintain the high time-frequency resolution.

Other than the processing techniques listed in this paper, there are still several techniques such as Hilbert-Huang transform, which have been used to analyze the EMG signals. There are also already well established techniques in other research areas and are proven to be effective and have yet to be implemented in analyzing EMG signals such as the S-transform. This technique may be suitable for nonstationary signals like EMG itself, due to its good nature, including linearity, lossless reversibility, multiple resolution, good time-frequency resolution, and simple algorithm. In conclusion, it can be seen that there are many gaps to be filled in this area in either finding new and better fatigue indices or constantly developing new and improved processing techniques.

Acknowledgments

The authors acknowledge and thank the Control and Signal Processing Group of Universiti Putra Malaysia and Rehabilitation and Assistive Technology Research Group of Universiti Teknikal Malaysia Melaka for their continuous support in the research. This research is fully funded by the Ministry of Higher Education Malaysia (MOHE).

References

[1] S. Gallagher, J. R. Heberger, "Examining the interaction of force and repetition on musculoskeletal disorder risk: a systematic literature review,", Human Factors , vol. 55, no. 1, pp. 108-124, 2013.

[2] E. M. Eatough, J. D. Way, C.-H. Chang, "Understanding the link between psychosocial work stressors and work-related musculoskeletal complaints,", Applied Ergonomics , vol. 43, no. 3, pp. 554-563, 2012.

[3] T. R. Waters, V. Putz-Anderson, A. Garg, L. J. Fine, "Revised NIOSH equation for the design and evaluation of manual lifting tasks,", Ergonomics , vol. 36, no. 7, pp. 749-776, 1993.

[4] Y. Roquelaure, C. Ha, A. Leclerc, A. Touranchet, M. Sauteron, M. Melchior, E. Imbernon, M. Goldberg, "Epidemiologic surveillance of upper-extremity musculoskeletal disorders in the working population,", Arthritis Care and Research , vol. 55, no. 5, pp. 765-778, 2006.

[5] H. Arif, M. S. E. Kosnan, K. Jusoff, W. Sapto, T. J. S. Anand, Taufik, Sivarao, M. H. Abu Bakar, "Finite element analysis of conceptual lumbar spine for different lifting position,", World Applied Sciences Journal , vol. 21, pp. 68-75, 2013.

[6] C. J. Colloca, R. N. Hinrichs, "The biomechanical and clinical significance of the lumbar erector spinae flexion-relaxation phenomenon: a review of literature,", Journal of Manipulative and Physiological Therapeutics , vol. 28, no. 8, pp. 623-631, 2005.

[7] N. H. Kamarudin, S. A. Ahmad, M. K. Hassan, R. M. Yusoff, S. Z. Dawal, "Muscle contraction analysis during lifting task," in Proceedings of the 3rd IEEE Conference on Biomedical Engineering and Sciences (IECBES '14), pp. 452-457, IEEE, Kuala Lumpur, Malaysia, December 2014.

[8] K. R. Voge, J. B. Dingwell, "Relative timing of changes in muscle fatigue and movement coordination during a repetitive one-hand lifting task," in Proceedings of the A New Beginning for Human Health: Proceddings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '03), pp. 1807-1810, IEEE, Cancun, Mexico, September 2003.

[9] T. N. S. T. Zawawi, A. R. Abdullah, E. F. Shair, I. Halim, S. M. Saleh, "EMG signal analysis of fatigue muscle activity in manual lifting,", Journal of Electrical Systems , vol. 11, no. 3, pp. 319-325, 2015.

[10] R. A. Al-Ashaik, M. Z. Ramadan, K. S. Al-Saleh, T. M. Khalaf, "Effect of safety shoes type, lifting frequency, and ambient temperature on subject's MAWL and physiological responses,", International Journal of Industrial Ergonomics , vol. 50, pp. 43-51, 2015.

[11] K. P. Granata, W. S. Marras, "An EMG-assisted model of trunk loading during free-dynamic lifting,", Journal of Biomechanics , vol. 28, no. 11, pp. 1309-1317, 1995.

[12] S. H. Roy, P. Bonato, M. Knaflitz, "EMG assessment of back muscle function during cyclical lifting,", Journal of Electromyography and Kinesiology , vol. 8, no. 4, pp. 233-245, 1998.

[13] R. E. Seroussi, M. H. Pope, "The relationship between trunk muscle electromyography and lifting moments in the sagittal and frontal planes,", Journal of Biomechanics , vol. 20, no. 2, pp. 135-146, 1987.

[14] R. B. Graham, M. J. Agnew, J. M. Stevenson, "Effectiveness of an on-body lifting aid at reducing low back physical demands during an automotive assembly task: assessment of EMG response and user acceptability,", Applied Ergonomics , vol. 40, no. 5, pp. 936-942, 2009.

[15] D. Gagnon, C. Larivière, P. Loisel, "Comparative ability of EMG, optimization, and hybrid modelling approaches to predict trunk muscle forces and lumbar spine loading during dynamic sagittal plane lifting,", Clinical Biomechanics , vol. 16, no. 5, pp. 359-372, 2001.

[16] P. Dolan, M. A. Adams, "The relationship between EMG activity and extensor moment generation in the erector spinae muscles during bending and lifting activities,", Journal of Biomechanics , vol. 26, no. 4-5, pp. 513-522, 1993.

[17] P. Bonato, P. Boissy, U. Della Croce, S. H. Roy, "Changes in the surface EMG signal and the biomechanics of motion during a repetitive lifting task,", IEEE Transactions on Neural Systems and Rehabilitation Engineering , vol. 10, no. 1, pp. 38-47, 2002.

[18] J. R. Potvin, "Mechanically corrected EMG for the continuous estimation of erector spinae muscle loading during repetitive lifting,", European Journal of Applied Physiology and Occupational Physiology , vol. 74, no. 1-2, pp. 119-132, 1996.

[19] I. Kingma, J. H. Van Dieën, "Lifting over an obstacle: effects of one-handed lifting and hand support on trunk kinematics and low back loading,", Journal of Biomechanics , vol. 37, no. 2, pp. 249-255, 2004.

[20] H.-J. Shin, J.-Y. Kim, "Measurement of trunk muscle fatigue during dynamic lifting and lowering as recovery time changes,", International Journal of Industrial Ergonomics , vol. 37, no. 6, pp. 545-551, 2007.

[21] J. Cholewicki, S. M. McGill, R. W. Norman, "Comparison of muscle forces and joint load from an optimization and EMG assisted lumbar spine model: towards development of a hybrid approach,", Journal of Biomechanics , vol. 28, no. 3, pp. 321-331, 1995.

[22] L. Straker, "Evidence to support using squat, semi-squat and stoop techniques to lift low-lying objects,", International Journal of Industrial Ergonomics , vol. 31, no. 3, pp. 149-160, 2003.

[23] C. K. Anderson, D. B. Chaffin, "A biomechanical evaluation of five lifting techniques,", Applied Ergonomics , vol. 17, no. 1, pp. 2-8, 1986.

[24] R. Burgess-Limerick, "Squat, stoop, or something in between?,", International Journal of Industrial Ergonomics , vol. 31, no. 3, pp. 143-148, 2003.

[25] S. M. Hsiang, G. E. Brogmus, T. K. Courtney, "Low back pain (LBP) and lifting technique--a review,", International Journal of Industrial Ergonomics , vol. 19, no. 1, pp. 59-74, 1997.

[26] A. Merlo, I. Campanini, "Technical aspects of surface electromyography for clinicians,", The Open Rehabilitation Journal , vol. 3, no. 1, pp. 98-109, 2010.

[27] E. N. Kamavuako, J. C. Rosenvang, R. Horup, W. Jensen, D. Farina, K. B. Englehart, "Surface versus untargeted intramuscular EMG based classification of simultaneous and dynamically changing movements,", IEEE Transactions on Neural Systems and Rehabilitation Engineering , vol. 21, no. 6, pp. 992-998, 2013.

[28] L. H. Smith, L. J. Hargrove, "Comparison of surface and intramuscular EMG pattern recognition for simultaneous wrist/hand motion classification," in Proceedings of the 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '13), pp. 4223-4226, Osaka, Japan, July 2013.

[29] L. J. Hargrove, K. Englehart, B. Hudgins, "A comparison of surface and intramuscular myoelectric signal classification,", IEEE Transactions on Biomedical Engineering , vol. 54, no. 5, pp. 847-853, 2007.

[30] A. Fratini, M. Cesarelli, P. Bifulco, M. Romano, "Relevance of motion artifact in electromyography recordings during vibration treatment,", Journal of Electromyography and Kinesiology , vol. 19, no. 4, pp. 710-718, 2009.

[31] C. J. De Luca, L. Donald Gilmore, M. Kuznetsov, S. H. Roy, "Filtering the surface EMG signal: movement artifact and baseline noise contamination,", Journal of Biomechanics , vol. 43, no. 8, pp. 1573-1579, 2010.

[32] H. L. Butler, R. Newell, C. L. Hubley-Kozey, J. W. Kozey, "The interpretation of abdominal wall muscle recruitment strategies change when the electrocardiogram (ECG) is removed from the electromyogram (EMG),", Journal of Electromyography and Kinesiology , vol. 19, no. 2, pp. e102-e113, 2009.

[33] G. Lu, J.-S. Brittain, P. Holland, J. Yianni, A. L. Green, J. F. Stein, T. Z. Aziz, S. Wang, "Removing ECG noise from surface EMG signals using adaptive filtering,", Neuroscience Letters , vol. 462, no. 1, pp. 14-19, 2009.

[34] N. W. Willigenburg, A. Daffertshofer, I. Kingma, J. H. van Dieën, "Removing ECG contamination from EMG recordings: a comparison of ICA-based and other filtering procedures,", Journal of Electromyography and Kinesiology , vol. 22, no. 3, pp. 485-493, 2012.

[35] M. I. Ibrahimy, M. R. Ahsan, O. O. Khalifa, "Design and optimization of Levenberg-Marquardt based neural network classifier for EMG signals to identify hand motions,", Measurement Science Review , vol. 13, no. 3, pp. 142-151, 2013.

[36] S. N. Sidek, A. J. Haja Mohideen, "Mapping of EMG signal to hand grip force at varying wrist angles," in Proceedings of the 2nd IEEE-EMBS Conference on Biomedical Engineering and Sciences (IECBES '12), pp. 648-653, IEEE, Langkawi, Malaysia, December 2012.

[37] Y. Zhan, D. Halliday, P. Jiang, X. Liu, J. Feng, "Detecting time-dependent coherence between non-stationary electrophysiological signals--a combined statistical and time-frequency approach,", Journal of Neuroscience Methods , vol. 156, no. 1-2, pp. 322-332, 2006.

[38] E. R. Ferrara, B. Widrow, "Fetal electrocardiogram enhancement by time-sequenced adaptive filtering,", IEEE Transactions on Biomedical Engineering , vol. 29, no. 6, pp. 458-460, 1982.

[39] R. H. Shumway, "Time-frequency clustering and discriminant analysis,", Statistics & Probability Letters , vol. 63, no. 3, pp. 307-314, 2003.

[40] D. Korosec, "Parametric estimation of the continuous non-stationary spectrum and its dynamics in surface EMG studies,", International Journal of Medical Informatics , vol. 58-59, pp. 59-69, 2000.

[41] G. Morren, S. Van Huffel, I. Helon, G. Naulaers, H. Daniëls, H. Devlieger, P. Casaer, "Effects of non-nutritive sucking on heart rate, respiration and oxygenation: a model-based signal processing approach,", Comparative Biochemistry and Physiology--A Molecular and Integrative Physiology , vol. 132, no. 1, pp. 97-106, 2002.

[42] R. Latif, E. Aassif, M. Laaboubi, A. Moudden, G. Maze, "Determination of the cut-off frequency of the anti-symmetric circumferential waves using the time-varying autoregressive TVAR," in Proceedings of the 3rd International Symposium on Communications, Control, and Signal Processing, (ISCCSP '08), pp. 357-361, IEEE, March 2008.

[43] R. Latif, E. Aassif, M. Laaboubi, G. Maze, "Determination of the group velocity using the time-varying autoregressive TVAR," in Proceedings of the 9th International Symposium on Signal Processing and Its Applications (ISSPA '07), pp. 1-4, IEEE, Sharjah, United Arab Emirates, February 2007.

[44] A. Al Zaman, X. Luo, M. Ferdjallah, A. Khamayseh, "A new TVAR modeling in cascaded form for nonstationary signals," in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '06), pp. 353-356, May 2006.

[45] A. Al Zaman, M. A. Ahad, M. Ferdjallah, J. J. Wertsch, "A new approach for muscle fatigue analysis in young adults at different MVC levels," in Proceedings of the IEEE International 48th Midwest Symposium on Circuits and Systems (MWSCAS '05), pp. 499-502, August 2005.

[46] Z. G. Zhang, H. T. Liu, S. C. Chan, K. D. K. Luk, Y. Hu, "Time-dependent power spectral density estimation of surface electromyography during isometric muscle contraction: methods and comparisons,", Journal of Electromyography and Kinesiology , vol. 20, no. 1, pp. 89-101, 2010.

[47] A. Al Zaman, M. Ferdjallah, A. Khamayseh, "Muscle fatigue analysis for healthy adults using TVAR model with instantaneous frequency estimation," in Proceedings of the 38th Southeastern Symposium on System Theory, pp. 44-47, March 2006.

[48] K. Englehart, B. Hudgins, P. A. Parker, M. Stevenson, "Classification of the myoelectric signal using time-frequency based representations,", Medical Engineering and Physics , vol. 21, no. 6-7, pp. 431-438, 1999.

[49] N. Selvaraj, J. Lee, K. H. Chon, "Time-varying methods for charac-terizing nonstationary dynamics of physiological systems,", Methods of Information in Medicine , vol. 49, no. 5, pp. 435-442, 2010.

[50] D. Tkach, H. Huang, T. A. Kuiken, "Study of stability of time-domain features for electromyographic pattern recognition,", Journal of NeuroEngineering and Rehabilitation , vol. 7, article 21, 2010.

[51] S. Du, M. Vuskovic, "Temporal vs. spectral approach to feature extraction from prehensile EMG signals," in Proceedings of the IEEE International Conference on Information Reuse and Integration (IRI '04), pp. 344-350, Las Vegas, Nev, USA, November 2004.

[52] R. M. Enoka, J. Duchateau, "Muscle fatigue: what, why and how it influences muscle function,", Journal of Physiology , vol. 586, no. 1, pp. 11-23, 2008.

[53] W. J. Williams, J. Jeong, "Kernel design for reduced interference distributions,", IEEE Transactions on Signal Processing , vol. 40, no. 2, pp. 402-412, 1992.

[54] M. B. I. Reaz, M. S. Hussain, F. Mohd-Yasin, "Techniques of EMG signal analysis: detection, processing, classification and applications,", Biological Procedures Online , vol. 8, no. 1, pp. 11-35, 2006.

[55] R. A. Malinzak, S. M. Colby, D. T. Kirkendall, B. Yu, W. E. Garrett, "A comparison of knee joint motion patterns between men and women in selected athletic tasks,", Clinical Biomechanics , vol. 16, no. 5, pp. 438-445, 2001.

[56] T. I. Arabadzhiev, V. G. Dimitrov, N. A. Dimitrova, G. V. Dimitrov, "Interpretation of EMG integral or RMS and estimates of 'neuromuscular efficiency' can be misleading in fatiguing contraction,", Journal of Electromyography and Kinesiology , vol. 20, no. 2, pp. 223-232, 2010.

[57] C. Suetta, P. Aagaard, A. Rosted, A. K. Jakobsen, B. Duus, M. Kjaer, S. P. Magnusson, "Training-induced changes in muscle CSA, muscle strength, EMG, and rate of force development in elderly subjects after long-term unilateral disuse,", Journal of Applied Physiology , vol. 97, no. 5, pp. 1954-1961, 2004.

[58] M. A. Oskoei, H. Hu, "Support vector machine-based classification scheme for myoelectric control applied to upper limb,", IEEE Transactions on Biomedical Engineering , vol. 55, no. 8, pp. 1956-1965, 2008.

[59] J. J. Villarejo, A. Frizera, T. F. Bastos, J. F. Sarmiento, "Pattern recognition of hand movements with low density sEMG for prosthesis control purposes," in Proceedings of the IEEE 13th International Conference on Rehabilitation Robotics (ICORR '13), Seattle, Wash, USA, June 2013.

[60] A. Phinyomark, C. Limsakul, P. Phukpattaranont, "A novel feature extraction for robust EMG pattern recognition,", Journal of Computing , vol. 1, no. 1, pp. 71-80, 2009.

[61] M. Zardoshti-Kermani, B. C. Wheeler, K. Badie, R. M. Hashemi, "EMG Feature Evaluation for Movement Control of Upper Extremity Prostheses,", IEEE Transactions on Rehabilitation Engineering , vol. 3, no. 4, pp. 324-333, 1995.

[62] D. Tkach, H. Huang, T. A. Kuiken, "Study of stability of time-domain features for electromyographic pattern recognition,", Journal of NeuroEngineering and Rehabilitation , vol. 7, no. 1, article no. 21, 2010.

[63] K. Kiguchi, S. Kariya, K. Watanabe, K. Izumi, T. Fukuda, "An exoskeletal robot for human elbow motion support--sensor fusion, adaptation, and control,", IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics , vol. 31, no. 3, pp. 353-361, 2001.

[64] F. Al Omari, J. Hui, C. Mei, G. Liu, "Pattern recognition of eight hand motions using feature extraction of forearm EMG signal,", Proceedings of the National Academy of Sciences India Section A--Physical Sciences , vol. 84, no. 3, pp. 473-480, 2014.

[65] A. R. Siddiqi, S. N. Sidek, A. Khorshidtalab, "Signal processing of EMG signal for continuous thumb-angle estimation," in Proceedings of the 41st Annual Conference of the IEEE Industrial Electronics Society (IECON '15), pp. 374-379, Yokohama, Japan, November 2015.

[66] Å. Kilbom, G. M. Hägg, C. Käll, "One-handed load carrying--cardiovascular, muscular and subjective indices of endurance and fatigue,", European Journal of Applied Physiology and Occupational Physiology , vol. 65, no. 1, pp. 52-58, 1992.

[67] F. AlOmari, G. Liu, "Analysis of extracted forearm sEMG signal using LDA, QDA, K-NN classification algorithms,", Open Automation and Control Systems Journal , vol. 6, no. 1, pp. 108-116, 2014.

[68] A. Phinyomark, G. Chujit, P. Phukpattaranont, C. Limsakul, H. Hu, "A preliminary study assessing time-domain EMG features of classifying exercises in preventing falls in the elderly," in Proceedings of the 9th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology (ECTI-CON '12), pp. 1-4, IEEE, Phetchaburi, Thailand, May 2012.

[69] D. R. Rogers, D. T. MacIsaac, "EMG-based muscle fatigue assessment during dynamic contractions using principal component analysis,", Journal of Electromyography and Kinesiology , vol. 21, no. 5, pp. 811-818, 2011.

[70] K. Buchenrieder, "Dimensionality reduction for the control of powered upper limb prostheses," in Proceedings of the 14th Annual IEEE International Conference and Workshops on the Engineering of Computer-Based Systems (ECBS '07), pp. 327-333, IEEE, Tucson, Ariz, USA, March 2007.

[71] G. Venugopal, M. Navaneethakrishna, S. Ramakrishnan, "Extraction and analysis of multiple time window features associated with muscle fatigue conditions using sEMG signals,", Expert Systems with Applications , vol. 41, no. 6, pp. 2652-2659, 2014.

[72] A. Phinyomark, C. Limsakul, P. Phukpattaranont, "EMG feature extraction for tolerance of white Gaussian noise," in Proceedings of the International Workshop and Symposium on Science and Technology (I-SEEC '08), pp. 178-183, December 2008.

[73] M. S. Al-Quraishi, A. J. Ishak, S. A. Ahmad, M. K. Hasan, "Impact of feature extraction techniques on classification accuracy for EMG based ankle joint movements," in Proceedings of the 10th Asian Control Conference (ASCC '15), pp. 1-5, Sabah, Malaysia, June 2015.

[74] G.-C. Chang, W.-J. Kang, L. Jer-Junn, C.-K. Cheng, J.-S. Lai, J.-J. J. Chen, T.-S. Kuo, "Real-time implementation of electromyogram pattern recognition as a control command of man-machine interface,", Medical Engineering and Physics , vol. 18, no. 7, pp. 529-537, 1996.

[75] A. Georgakis, L. K. Stergioulas, G. Giakas, "Fatigue analysis of the surface EMG signal in isometric constant force contractions using the averaged instantaneous frequency,", IEEE Transactions on Biomedical Engineering , vol. 50, no. 2, pp. 262-265, 2003.

[76] R. Merletti, L. R. Lo Conte, "Surface EMG signal processing during isometric contractions,", Journal of Electromyography and Kinesiology , vol. 7, no. 4, pp. 241-250, 1997.

[77] C. J. De Luca, "Use of the surface EMG signal for performance evaluation of back muscles,", Muscle and Nerve , vol. 16, no. 2, pp. 210-216, 1993.

[78] F. Khanam, M. Ahmad, "Frequency based EMG power spectrum analysis of Salat associated muscle contraction," in Proceedings of the 1st International Conference on Electrical and Electronic Engineering (ICEEE '15), Rajshahi, Bangladesh, November 2015.

[79] M. M. Altaf, B. M. Elbagoury, F. Alraddady, M. Roushdy, "Extended case-based behavior control for multi-humanoid robots,", International Journal of Humanoid Robotics , vol. 13, no. 2, 2016.

[80] A. Phinyomark, A. Nuidod, P. Phukpattaranont, C. Limsakul, "Feature extraction and reduction of wavelet transform coefficients for EMG pattern classification,", Elektronika ir Elektrotechnika , vol. 122, no. 6, pp. 27-32, 2012.

[81] M. González-Izal, A. Malanda, I. Navarro-Amézqueta, E. M. Gorostiaga, F. Mallor, J. Ibañez, M. Izquierdo, "EMG spectral indices and muscle power fatigue during dynamic contractions,", Journal of Electromyography and Kinesiology , vol. 20, no. 2, pp. 233-240, 2010.

[82] G. V. Dimitrov, T. I. Arabadzhiev, K. N. Mileva, J. L. Bowtell, N. Crichton, N. A. Dimitrova, "Muscle fatigue during dynamic contractions assessed by new spectral indices,", Medicine and Science in Sports and Exercise , vol. 38, no. 11, pp. 1971-1979, 2006.

[83] M. Vitor-Costa, H. Bortolotti, T. V. Camata, J. R. Cardoso, R. A. da Silva, T. Abrão, A. C. de Moraes, L. R. Altimari, "EMG spectral analysis of incremental exercise in cyclists and non-cyclists using fourier and wavelet transforms,", Revista Brasileira de Cineantropometria e Desempenho Humano , vol. 14, no. 6, pp. 660-670, 2012.

[84] M. Wacker, H. Witte, "Time-frequency techniques in biomedical signal analysis: a tutorial review of similarities and differences,", Methods of Information in Medicine , vol. 52, no. 4, pp. 279-296, 2013.

[85] S. Karlsson, J. Yu, M. Akay, "Time-frequency analysis of myoelectric signals during dynamic contractions: a comparative study,", IEEE Transactions on Biomedical Engineering , vol. 47, no. 2, pp. 228-238, 2000.

[86] A. L. Martinez-Herrera, L. M. Ledesma-Carrillo, M. Lopez-Ramirez, S. Salazar-Colores, E. Cabal-Yepez, A. Garcia-Perez, "Gabor and the Wigner-Ville transforms for broken rotor bars detection in induction motors," in Proceedings of the 24th International Conference on Electronics, Communications and Computers (CONIELECOMP '14), pp. 83-87, IEEE, Cholula. Mexico, February 2014.

[87] D. MacIsaac, P. A. Parker, R. N. Scott, "The short-time Fourier transform and muscle fatigue assessment in dynamic contractions,", Journal of Electromyography and Kinesiology , vol. 11, no. 6, pp. 439-449, 2001.

[88] E. F. Shair, A. R. Abdullah, T. N. S. Tengku Zawawi, S. A. Ahmad, S. Mohamad Saleh, "Auto-segmentation analysis of EMG signal for lifting muscle contraction activities,", Journal of Telecommunication, Electronic and Computer Engineering , vol. 8, no. 7, pp. 17-22, 2016.

[89] M. M. Jankovic, D. B. Popovic, "An EMG system for studying motor control strategies and fatigue," in Proceedings of the 10th Symposium on Neural Network Applications in Electrical Engineering (NEUREL-'10), pp. 15-18, Belgrade, Serbia, September 2010.

[90] T. N. S. T. Zawawi, A. R. Abdullah, I. Halim, E. F. Shair, S. M. Salleh, "Application of spectrogram in analysing electromyography (EMG) signals of manual lifting,", ARPN Journal of Engineering and Applied Sciences , vol. 11, no. 6, pp. 3603-3609, 2016.

[91] M. Yochum, T. Bakir, R. Lepers, S. Binczak, "Estimation of muscular fatigue under electromyostimulation using CWT,", IEEE Transactions on Biomedical Engineering , vol. 59, no. 12, pp. 3372-3378, 2012.

[92] J. L. Dantas, T. V. Camata, M. A. O. C. Brunetto, A. C. Moraes, T. Abrão, L. R. Altimari, "Fourier and Wavelet spectral analysis of EMG signals in isometric and dynamic maximal effort exercise," in Proceedings of the 32nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '10), pp. 5979-5982, IEEE, Buenos Aires, Argentina, September 2010.

[93] L. Ai, J. Nan, J. Shen, "Application of S-transform to driving fatigue in EEG analysis," in Proceedings of the IEEE International Conference on Intelligent Computing and Intelligent Systems (ICIS '10), pp. 814-817, IEEE, Guilin, China, October 2010.

[94] R. Upadhyay, P. K. Padhy, P. K. Kankar, "EEG artifact removal and noise suppression by discrete orthonormal S-transform denoising,", Computers and Electrical Engineering , vol. 53, pp. 125-142, 2016.

[95] C. Dora, P. K. Biswal, "Robust ECG artifact removal from EEG using continuous wavelet transformation and linear regression," in 2016 International Conference on Signal Processing and Communications (SPCOM), pp. 1-5, Bangalore, India, June 2016.

[96] S. Assous, B. Boashash, "Evaluation of the modified S -transform for time-frequency synchrony analysis and source localisation,", Eurasip Journal on Advances in Signal Processing , vol. 2012, article 49, 2012.

[97] A. L. Ricamato, R. G. Absher, M. T. Moffroid, J. P. Tranowski, "A time-frequency approach to evaluate electromyographic recordings," in Proceedings of the 5th Annual IEEE Symposium on Computer-Based Medical Systems, pp. 520-527, IEEE, Durham, NC, USA, 1992.

[98] M. Khezri, M. Jahed, "An inventive quadratic time-frequency scheme based on Wigner-Ville distribution for classification of sEMG signals," in Proceedings of the 6th International Special Topic Conference on Information Technology Applications in Biomedicine, pp. 261-264, IEEE, Tokyo, Japan, November 2007.

[99] A. Subasi, M. K. Kiymik, "Muscle fatigue detection in EMG using time-frequency methods, ICA and neural networks,", Journal of Medical Systems , vol. 34, no. 4, pp. 777-785, 2010.

[100] W. Martin, P. Flandrin, "Wigner-Ville Spectral Analysis of Nonstationary Processes,", IEEE Transactions on Acoustics, Speech, and Signal Processing , vol. 33, no. 6, pp. 1461-1470, 1985.

[101] H.-I. Choi, W. J. Williams, "Improved time-frequency representation of multicomponent signals using exponential kernels,", IEEE Transactions on Acoustics, Speech, and Signal Processing , vol. 37, no. 6, pp. 862-871, 1989.

[102] G. R. Pereira, L. F. de Oliveira, J. Nadal, "Reducing cross terms effects in the Choi-Williams transform of mioelectric signals,", Computer Methods and Programs in Biomedicine , vol. 111, no. 3, pp. 685-692, 2013.

[103] M. Alemu, D. K. Kumar, A. Bradley, "Time-frequency analysis of SEMG--with special consideration to the interelectrode spacing,", IEEE Transactions on Neural Systems and Rehabilitation Engineering , vol. 11, no. 4, pp. 341-345, 2003.

[104] P. A. Karthick, G. Venugopal, S. Ramakrishnan, "Analysis of surface EMG signals under fatigue and non-fatigue conditions using B-distribution based quadratic time frequency distribution,", Journal of Mechanics in Medicine and Biology , vol. 15, no. 2, 2015.

[105] V. Sucic, B. Barkat, B. Boashash, "Performance evaluation of the B-distribution," in Proceedings of the 5th International Symposium on Signal Processing and Its Applications (ISSPA '99), pp. 267-270, Queensland, Australia, August 1999.

[106] P. A. Karthick, S. Ramakrishnan, "Surface electromyography based muscle fatigue progression analysis using modified B distribution time-frequency features,", Biomedical Signal Processing and Control , vol. 26, pp. 42-51, 2016.

[107] B. Boashash, V. Sucic, "Resolution measure criteria for the objective assessment of the performance of quadratic time-frequency distributions,", IEEE Transactions on Signal Processing , vol. 51, no. 5, pp. 1253-1263, 2003.

[108] S. Dong, G. Azemi, B. Boashash, "Improved characterization of HRV signals based on instantaneous frequency features estimated from quadratic time-frequency distributions with data-adapted kernels,", Biomedical Signal Processing and Control , vol. 10, no. 1, pp. 153-165, 2014.

[109] B. Boashash, T. Ben-Jabeur, "Design of a high-resolution separable-kernel quadratic TFD for improving newborn health outcomes using fetal movement detection," in Proceedings of the 11th International Conference on Information Science, Signal Processing and their Applications (ISSPA '12), pp. 354-359, Québec, Canada, July 2012.