The Scientific World Journal

Volume 2012 (2012), Article ID 478952, 5 pages

doi:10.1100/2012/478952

Research Article

Analysis of EMG Signals in Aggressive and Normal Activities by Using Higher-Order Spectra

Necmettin Sezgin

Department of Electrical and Electronics Engineering, Faculty of Architecture and Engineering, Batman University, 72060 Batman, Turkey

Received 2 September 2012; Accepted 12 October 2012

Academic Editors: T. Arendt and J. Benito-Leon

Copyright © 2012 Necmettin Sezgin. 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.

Abstract.

The analysis and classification of electromyography (EMG) signals are very important in order to detect some symptoms of diseases, prosthetic arm/leg control, and so on. In this study, an EMG signal was analyzed using bispectrum, which belongs to a family of higher-order spectra. An EMG signal is the electrical potential difference of muscle cells. The EMG signals used in the present study are aggressive or normal actions. The EMG dataset was obtained from the machine learning repository. First, the aggressive and normal EMG activities were analyzed using bispectrum and the quadratic phase coupling of each EMG episode was determined. Next, the features of the analyzed EMG signals were fed into learning machines to separate the aggressive and normal actions. The best classification result was 99.75%, which is sufficient to significantly classify the aggressive and normal actions.

1. Introduction

Electromyography (EMG) is the electrical activity of muscle cells and has been used for the classification of actions [1–7], disease detection [8], prosthetic hand control [9], and emotion detection [10]. In this study, 8 channels recorded the EMG signals of 10 aggressive and 10 normal actions of 3 males and 1 female, which were then analyzed for classifying normal and aggressive actions. With this purpose in mind, a composed model of higher-order spectra (HOS) and the learning machine algorithm was proposed. In signal processing, 2nd-order statistics methods such as the power spectrum have gained significant importance. However, many signals have nonlinearity and non-Gaussian behavior, and such signals cannot be examined properly by 2nd-order statistical methods. Thus, higher-order statistical methods have been proved. The HOS was first applied to real signal processing problems in the 1970s and since then it has continued to be applied in many different areas, such as economics, speech signal processing, noisy and artifact removal, biomedical signal processing, and optics. Since EMG signals are nonstationary and non-Gaussian signals, they should be examined by HOS methods. Bispectrum, which is the Fourier transform of the 3rd-order cumulant, can be applied to nonlinear and non-Gaussian signals to extract nonlinear information.

Bispectrum analysis reveals phase information called quadratic phase coupling (QPC). In the present study, the EMG signals analyzed using bispectrum and the QPCs were determined for all of the datasets, and then these QPCs were fed into the extreme learning machine (ELM) algorithm. The ELM is capable of training and testing data fast and with a high accuracy. The main advantage of ELM over the traditional learning methods is that it is very fast due to its algorithm. In the ELM algorithm, the weights between the input layer and the hidden layer and the hidden layer’s biases are selected randomly, while the weights between the hidden layer and the output layer are determined analytically. Therefore, considerable time saving is attained in the training stage. Moreover, the performance of the classification method (ELM) was compared with some other machine learning methods, such as support vector machine (SVM), logistic regression (LR), linear discriminant analysis (LDA), and artificial neural network (ANN). The proposed method is satisfactory due to the compared classification methods.

2. Material and Methods

2.1. Database

In this study, the dataset of the “EMG physical action data set” from the machine learning repository (UCI) [11] was used. 3 male and 1 female subjects took part in the experiment (aged 25 to 30 years), who have experienced aggression in scenarios such as physical fighting. Each subject had to perform 10 normal and 10 aggressive activities. The normal activities were bowing, clapping, handshaking, hugging, jumping, running, seating, standing, walking, and waving, while the aggressive activities were elbowing, front kicking, hammering, headering, kneeing, pulling, punching, pushing, side kicking, and slapping. There were 8 electrodes used, which corresponds to 8 input time series, one for each muscle channel (ch1–8): right bicep (ch1), right tricep (ch2), left bicep (ch3), left tricep (ch4), right thigh (ch5), right hamstring (ch6), left thigh (ch7), and left hamstring (ch8). Each time series contained about 10,000 samples, which were 10 s in length.

2.2. Bispectrum Analysis

Bispectrum analysis reveals the phase relation between components of a signal [12–14]. Unlike the power spectrum, the bispectrum is capable of extracting extra information from biological signals such as an EMG signal, which is non-Gaussian and nonlinear. The bispectrum is defined as the Fourier transform of the 3rd-order cumulant.

The 3rd-order cumulant of a discrete signal , which is stationary and has a 0 mean, is defined as [14]

where denotes the expected process.

The th degree moment of is defined as

Thus, (1) can be rewritten as:

Alternatively, the 3rd-order cumulant can be written as

where is the 3rd-order moment function of and is the 3rd-order moment function of a Gaussian random process with the same 1st- and 2nd-order characteristics of

An important result of (4) is that if is a Gaussian process, then its 3rd-order cumulant is 0 [14, 15]:

The 3rd-order cumulant and 3rd-order moment of a process, which has a 0 mean, are equal to each other. Thus, (4) becomes

The correlation is a relation between 2 points, whereas the 3rd-order cumulant is a relation between combinations of 3 points in a time series. The 3rd-order cumulant has symmetry properties as

The Fourier transform of the 3rd-order cumulant is bispectrum and defined as

where is the 2-dimensional window function that decreases the variance of the bispectrum. In this study, a Hanning window, which is 0.05 s in duration, was used. Equation (9) can also be defined in the Fourier transform of as

where * denotes a complex conjugate.

is a symmetric function, such that a triangular region , could completely describe the whole bispectrum. The other regions in the bispectrum are the symmetry of the defined triangular region. A peak observed in the triangular region indicates that the energy component at frequency is produced, likely due to the quadratic nonlinearity dependence, called QPC [16]. On the contrary, a flat bispectrum at the 2 frequency components and suggests no such activities. Consequently, phase coupled components contribute extensively to the 3rd-order cumulant sequence of a process. This unique capability of bispectral analysis becomes a useful tool to detect and quantify the possible existence of QPCs in the EMG signals of aggressive activities. To quantify the QPC, one can take advantage of the quantification of non-Gaussianity, which has a direct relation to phase coupling, of a random process as the sum of the magnitudes of the estimated bispectrum given by [17]:

The bispectrum quantity of all of the episodes in the database was determined through (11) and fed as input into the ELM classifier in order to separate aggressive activities from normal activities.

2.3. Extreme Learning Machine Algorithm

In the ELM, the network has 3 layers: input, output, and 1 hidden layer. The weights between the input and hidden layers and the hidden layer biases are selected randomly and , the hidden layer output matrix, is determined analytically. In the ELM, the training of the network is to minimize the sum square error for the as

where is the minimum square form of and is the weights between the hidden and output layers. is the inverse of the generalization Moore-Penrose of . The estimations are the minimum of the solution of the sum square of .

The ELM does not only find the minimum error, but can also achieve the best performance with respect to conventional gradient based methods. This performance arises by the singularity of the matrix . is a singular solution. The ELM algorithm can be summarized in 3 steps as follows [18].

(i)The weights , which are between the input layer and the hidden layer, and the hidden layer biases , are selected randomly.(ii)The output of the hidden layer, , is determined. (iii)The weights , which are between the hidden layer and the output layer, are calculated as , where is the target vector.

3. Results

In the present study, the ELM was used in order to classify the EMG signals as either belonging to an aggressive action or a normal action. In the 1st stage of the QPC, each 10 s EMG episode was determined by bispectral analysis. After bispectral analysis of the EMG signal, in the 2nd stage, the extracted features, which are the QPC quantity, were fed into the input of the ELM classifier. For the ELM algorithm, the training-testing rate was randomly chosen as 50%-50% from the extracted features of the EMG.

An example of normal EMG activity (waving) and aggressive activity (frontkicking) is shown in Figures 1 and 2, respectively. In Figures 1 and 2, normal and aggressive EMG actions (Figures 1(a) and 2(a)) and their corresponding power spectrums (Figures 1(b) and 2(b)), bispectrums (Figures 1(c) and 2(c)), and bispectrums in 2 dimensions (Figures 1(d) and 2(d)) are shown, respectively.

(a)

(b)

(c)

(d)

Figure 1: (a) The EMG activity of normal action, (b) its power spectrum, (c) its bispectrum, and (d) its bispectrum in 2 dimensions.

(a)

(b)

(c)

(d)

Figure 2: (a) The EMG activity of aggressive action, (b) its power spectrum, (c) its bispectrum, and (d) its bispectrum in 2 dimensions.

As shown in Figure 1, the bispectrum (Figure 1(d)) is about 20 times higher than the power spectrum (Figure 1(b)), and in Figure 2, the bispectrum (Figure 2(d)) is about 100,000 times higher than the power spectrum (Figure 2(b)). This means that the nonlinearity and non-Gaussian signals are increased rapidly by aggressive actions. Accordingly, the bispectrum of aggressive activity (Figures 2(c) and 2(d)) is much higher than normal activity (Figures 1(c) and 1(d)). Furthermore, the nonlinearity and non-Gaussian concentration of aggressive activity (Figure 2(c)) are centered at lower frequencies than normal activities (Figure 1(c)).

The ELM algorithm built in this study has 1 node in the input layer, 40 nodes in the hidden layer, and 2 nodes in the output layer. The activation function was sigmoid since it allowed the best performance. The number of nodes in the hidden layer and the type of activation function were chosen by trial and error. The classification accuracy for the ELM was 99.75% and the duration of the training-testing phases was 0.07 s and 0.005 s, respectively. The same data were also classified by ANN, SVM, LR, and LDA classifiers. Their results and a comparison with the ELM is shown in Table 1. It is clear that the ELM is faster and has higher accuracy performance than conventional learning machines.

Table 1: Performances of the ANN, SVM, LR, LDA, and ELM learning machines.

Model
Training process time (s)
Testing process time (s)
Accuracy (%)

ANN
32.25
1.18
98.20

SVM
1.80
0.20
96.15

LR
0.10
0.05
97.50

LDA
0.09
0.04
97.25

ELM
0.07
0.005
99.75

4. Discussion and Conclusions

In the present study, aggressive and normal EMG signals were analyzed using bispectrum and the EMG signals were classified in aggressive and normal activities using learning machine algorithms. A best performance was obtained using the ELM algorithm, which has an accuracy of 99.75%. The bispectrum of the aggressive and normal activities of the EMG were analyzed and the QPC quantities of each episode were determined and fed into the input of the learning machines. Aside from the classification of the aggressive and normal activities, a comparison of the power spectrum and bispectrum was performed in the EMG signals. Thus, the bispectrum of the EMG signal is a candidate to separate aggressive and normal activities. The results obtained from the ELM classifier are acceptably high enough to differentiate aggressive activities from normal activities. This simple and effective method may help experts in defining aggressive activities and this can give important clues about some abnormalities related to EMG signals.

Conflict of Interests

The author declares that he has no conflict of interests.

Acknowledgments

The author would like to thank and acknowledge the UC Irvine Machine Learning Repository. This study was supported by Batman University, under Project number BTUBAP-2011-MMF-5.

References

K. Nazarpour, A. R. Sharafat, and S. M. P. Firoozabadi, “Surface EMG signal classification using a selective mix of higher order statistics,” in Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEE-EMBS '05), pp. 4208–4211, September 2005.

A. Phinyomark, S. Hirunviriya, C. Limsakul, and P. Phukpattaranont, “Evaluation of EMG feature extraction for hand movement recognition based on euclidean distance and standard deviation,” in Proceedings of the 7th Annual International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology (ECTI-CON '10), pp. 856–860, May 2010.

Y. Al-Assaf, “Surface myoelectric signal analysis: dynamic approaches for change detection and classification,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 11, pp. 2248–2256, 2006.

Z. Gao, J. Lei, Q. Song, Y. Yu, and Y. Ge, “Research on the surface EMG signal for human body motion recognizing based on arm wrestling robot,” in Proceedings of the IEEE International Conference on Information Acquisition (ICIA '06), pp. 1269–1273, August 2006.

K. Englehart, B. Hudgins, P. Parker, and M. Stevenson, “Time-frequency representation for classification of the transient myoelectric signal,” in Proceedings of the 20th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 20, no. 5, pp. 2627–2630, November 1998.

X. Zhang, Y. Yang, X. Xu, and M. Zhang, “Wavelet-based neuro-fuzzy classification for EMG control,” in Proceedings of the IEEE International Conference on Communications, Circuits and Systems and West Sino Expositions, vol. 2, pp. 1087–1089, 2002.

Z. Xizhi, “Study of surface electromyography signal based on wavelet transform and radial basis function neural network,” in Proceedings of the International Seminar on Future BioMedical Information Engineering (FBIE '08), pp. 160–163, December 2008.

M. Bodruzzaman, M. Wilkes, R. Shiavi, and A. Kilroy, “Classification of electromyographic signals by autoregressive modeling,” in Proceedings of the IEEE Southeastcon on Technologies Today and Tomorrow, vol. 2, pp. 508–510, April 1990.

Z. Jingdong, X. Zongwu, J. Li, C. Hegao, L. Hong, and G. Hirzinger, “EMG control for a five-fingered prosthetic hand based on wavelet transform and autoregressive model,” in Proceedings of the IEEE International Conference on Mechatronics and Automation (ICMA '06), pp. 1097–1102, June 2006.

M. Murugappan, “Electromyogram signal based human emotion classification using KNN and LDA,” in Proceedings of the IEEE International Conference on System Engineering and Technology, pp. 106–110, 2011.

A. Frank and A. Asuncion, UCI Machine Learning Repository, University of California, School of Information and Computer Science, Irvine, Calif, USA, 2010.

J. C. Sigl and N. G. Chamoun, “An introduction to bispectral analysis for the electroencephalogram,” Journal of Clinical Monitoring, vol. 10, no. 6, pp. 392–404, 1994.

M. J. Hinich and C. S. Clay, “The application of the discrete Fourier transform in the estimation of power spectra, coherence and bispectra of geophysical data,” Reviews of Geophysics, vol. 6, no. 3, pp. 347–363, 1968.

C. L. Nikias and A. P. Petropulu, Higher-Order Spectral Analysis: A Nonlinear Signal Processing Framework, Prentice-Hall, Engle-wood Cliffs, NJ, USA, 1993.

C. L. Nikias and M. R. Raghuveer, “Bispectrum estimation: a digital signal processing framework,” Proceedings of the IEEE, vol. 75, no. 7, pp. 869–891, 1987.

M. R. Raghuveer and C. L. Nikias, “Bispectrum estimation: a parametric approach,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 33, no. 5, pp. 1213–1230, 1985.

T. Ning and J. D. Bronzino, “Bispectral analysis of the rat EEG during various vigilance states,” IEEE Transactions on Biomedical Engineering, vol. 36, no. 4, pp. 497–499, 1989.

G. B. Huang, Q. Y. Zhu, and C. K. Siew, “Extreme learning machine: theory and applications,” Neurocomputing, vol. 70, no. 1–3, pp. 489–501, 2006.