Next Article in Journal
Cooperative Detection of Multiple Targets by the Group of Mobile Agents
Previous Article in Journal
Neural Computing Enhanced Parameter Estimation for Multi-Input and Multi-Output Total Non-Linear Dynamic Models
Previous Article in Special Issue
Particle Swarm Optimisation: A Historical Review Up to the Current Developments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Emotional State Recognition from Peripheral Physiological Signals Using Fused Nonlinear Features and Team-Collaboration Identification Strategy

1
School of Mechanical Engineering, Changzhou University, Changzhou 213164, China
2
Remote Measurement and Control Key Lab of Jiangsu Province, School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(5), 511; https://doi.org/10.3390/e22050511
Submission received: 25 March 2020 / Revised: 25 April 2020 / Accepted: 27 April 2020 / Published: 30 April 2020
(This article belongs to the Special Issue Intelligent Tools and Applications in Engineering and Mathematics)

Abstract

:
Emotion recognition realizing human inner perception has a very important application prospect in human-computer interaction. In order to improve the accuracy of emotion recognition, a novel method combining fused nonlinear features and team-collaboration identification strategy was proposed for emotion recognition using physiological signals. Four nonlinear features, namely approximate entropy (ApEn), sample entropy (SaEn), fuzzy entropy (FuEn) and wavelet packet entropy (WpEn) are employed to reflect emotional states deeply with each type of physiological signal. Then the features of different physiological signals are fused to represent the emotional states from multiple perspectives. Each classifier has its own advantages and disadvantages. In order to make full use of the advantages of other classifiers and avoid the limitation of single classifier, the team-collaboration model is built and the team-collaboration decision-making mechanism is designed according to the proposed team-collaboration identification strategy which is based on the fusion of support vector machine (SVM), decision tree (DT) and extreme learning machine (ELM). Through analysis, SVM is selected as the main classifier with DT and ELM as auxiliary classifiers. According to the designed decision-making mechanism, the proposed team-collaboration identification strategy can effectively employ different classification methods to make decision based on the characteristics of the samples through SVM classification. For samples which are easy to be identified by SVM, SVM directly determines the identification results, whereas SVM-DT-ELM collaboratively determines the identification results, which can effectively utilize the characteristics of each classifier and improve the classification accuracy. The effectiveness and universality of the proposed method are verified by Augsburg database and database for emotion analysis using physiological (DEAP) signals. The experimental results uniformly indicated that the proposed method combining fused nonlinear features and team-collaboration identification strategy presents better performance than the existing methods.

1. Introduction

Emotions play an important role in human daily life. Currently, emotion recognition has received more and more attention in many fields, such as safe driving, distance education, health care and rehabilitation medical treatment [1]. Reliable and accurate emotion recognition system is one of the key problems of achieving natural human-machine interaction (HMI) [2].
In the last few decades, a variety of approaches for detecting human emotion have been performed by using speech, facial expression and behavior (gesture/posture) or physiological signals [3]. Leila et al. [4] proposed a global optimal feature fusion method for speech emotion recognition based on empirical mode decomposition and Teager-Kaiser Energy Operator (EMD-TKEO), according to the fact that the EMD combined with the TKEO gives an efficient time-frequency analysis of the non-stationary signals. In order to increase the accuracy rate of emotion recognition, unsupervised deep belief network (DBN) was proposed for depth level feature extraction from fused observations of Electro-Dermal Activity (EDA), Photoplethysmogram (PPG) and Zygomaticus Electromyography (zEMG) sensors signals and the fine Gaussian support vector machine was employed [5]. Hu et al. [6] proposed an integrated framework of two networks—a local network and a global network, which were based on local enhanced motion history image (LEMHI) and CNN-long short term memory (LSTM) cascaded networks respectively, for facial emotion recognition in video sequences. Picard et al. [7] proved that physiological signals can be used to classify emotions and successfully identified eight emotional states of calm, anger, disgust, sadness, peace, excitement, happiness and awe, using four physiological signals—electromyogram (EMG), pulse wave, skin conductivity (SC) and respiration (RSP). Many researchers have carried out researches on emotion recognition based on multi-modal physiological signals and EEG signals, many methods of emotion recognition were proposed, meanwhile, free emotion recognition databases with electroencephalogram (EEG) and peripheral neurophysiological signals were set up [8,9,10,11]. However, the informations of speech, facial expression and postures are more subjective [8]. Studies have shown that the signals controlled by the human nervous system [12,13,14,15,16,17,18], are closer to present people’s real emotions. Currently, using physiological signals for emotion recognition is one of the research hotspots.
In terms of using physiological signals to identify emotions, the core work of emotion identification mainly includes feature extraction and classification. In the aspect of feature extraction, numerous different features have been investigated for emotion recognition, from time domain, frequency domain, time-frequency domain and others. The extracted features based on the above mentioned analysis methods are considered inadequate when comes to use signals associated with the central nervous system (CNS) for emotion identification. On the one hand, a large number of features are needed to be extracted to represent the characteristics of each signal. On the other hand, the CNS related signals are non-stationary signals with characteristics varying over time and the traditional extraction methods can only achieve some superficial information [19,20,21]. However, the use of nonlinear analysis has the potential to provide discriminant features of CNS related signals that can enhance the classification accuracy of emotional states. In fact, nonlinear features have been successfully employed in a variety of ways. He et al. [22] employed wavelet packet entropy (WpEn) of speech signals as the feature to represent the emotion of angry, neutral and soft and achieved the accuracy of emotion recognition with 76%. Wanget al. [23] adopted three features, the approximate entropy (ApEn), hurst exponent and fractal dimension as extracted features from EEG and the performances of each feature with support vector machine (SVM) presented with 65.12% (ApEn), 71.38% (hurst exponent) and 70.63% (fractal dimension), respectively. Vayrynen et al [24] studied four type of emotional states (neutral, sad, angry and happy) with speech signals and employed two nonlinear features and k-nearest neighbour (KNN) classifier to realize emotion identification. In emotion identification, classifier also plays a very important role. The classifiers, SVM, DT, KNN and so on, are usually selected to realize emotion identification. Hosseini et al. [25] applied two entropy features and SVM to perform two types of emotion recognition with EEG. The sample entropy and SVM were adopted to realize two binary emotion recognition from EEG, with classification performances for positive and negative emotion by 80.43% and 79.11%, respectively [26]. Li et al. [27] adopted SVM to identify eight kinds of positive and negative emotions, combining the extracted features wavelet entropy (WaEn), ApEn and husrt exponent of EEG. Experimental results showed that the method achieved 85% recognition accuracy for two types of emotion recognition. Mohammadi et al [28] calculated the WaEn of each frequency band from EEG as the input feature of KNN classifier and achieved the classification accuracy of 86.75% for arousal level and 84.05% for valence level. According to the current research status, EEG signals are mainly selected for emotion recognition using signals related to the nervous system and the application of nonlinear features is relatively simple and not well integrated. In addition, when performing emotion recognition, a single classifier is usually used for recognition and determination.
As metioned above, physiological signals controlled by the human nervous system can also directly reflect emotional states. Because physiological signals are easy to obtain, this paper conducts research on emotion recognition from the perspective of using physiological signals. Previous researches of emotion recognition have largely been conducted on single-modality-based methods [12,17,24,25,26]. Compared with single-modality-based methods, fusing multiple sensors’ data is helpful to reflect emotional states from multiple perspectives, because the use of more than two modalities can exploit the complementary nature of different modalities so as to increase the credibility of measurements. In order to better represent the emotional state, this paper proposes the feature representation method of multi-nonlinear feature integration and multi-channel information feature fusion. Each classifier has its own performance characteristics. Samples that are misclassified by one classifier may easily be correctly classified by another classifier. In order to improve the accuracy of emotion identification, a team-collaboration dentification strategy is proposed, which can use different identification strategies to identify samples according to the characteristics of samples and the designed decision-making mechanism.

2. Materials and Methods

In this investigation, in order to effectively realize emotion identification, feature fusion of multi-channel signals is suggested to represent the emotional states and team-collaboration identification strategy is proposed to achieve emotion recognition. Two public databases, Augsburg dataset and DEAP dataset, are employed to verify the effectiveness and performances of the proposed methods.

2.1. Database

2.1.1. Augsburg Dataset and Data Pre-Processing

The public emotional dataset [29] from University of Augsburg records the different emotional states caused by corresponding musics (https://www.informatik.uni-augsburg.de/en/chairs/hcm/projects/tools/aubt/). To induce the subject to unaffectedly feel four different emotions, joy, anger, sadness and pleasure, some music songs that the subject himself carefully handpicked were employed. The chosen songs can bring back some spectial memories to the subject in respect of targeted emotion classes. When the participant listened to music songs, four kinds of physiological signals were recorded, including electromyogram (ECG), electrocardiogram (EMG), respiration change (RSP) and skin conductivity (SC). The sampling frequency of the ECG was 256 Hz and the other three signals were 32 Hz. Each recording of physiological signal was chosen to be 120s in length. The experiment was performed for 25 days and 25 samples of each emotional state were obtained, totally 100 samples (four types of motional states). The content summary of Augsburg Database is shown in Table 1. In this investigation, each 120s signal was segmented into ten 12s-samples without overlapping. Finally, a total of 1000 samples were obtained.

2.1.2. DEAP Dataset and Data Pre-Processing

Another publicly available database—DEAP [30] (http://www.eecs.qmul.ac.uk/mmv/datasets/deap/download.html)—was also adopted to investigate the universality of proposed methods. The DEAP dataset consists of 32 EEG channels and 8 peripheral physiological signals (electro-oculogram (EOG), galvanic skin response (GSR), blood volume pulse (BVP), RSP, EMG, skin temperature (SKT)) of 32 healthy subjects (labeled from s01 to s32) with half females and aged between 19 and 37. The EEG and physiological signals were recorded while each subject was watching 40 one-minute videos which were carefully selected to elicit different emotional states. In DEAP, each video clip was rated from 1 to 9 for arousal and valence by each subject after watching 40 trials of videos and the discrete rating value can be used as a classification label in emotion recognition. The total number of trials within the DEAP dataset was 1280 trials (40 trials per subject).
In this research, four dimensional emotional states (high valence-high arousal (HVHA), high valence-low arousal (HVLA), low valence-low arousal (LVLA), low valence-high arousal (LVHA)) were considered as target emotions based on the scales of valence and arousal. Four kinds of peripheral physiological signals GSR, BVP, RSP, EMG from DEAP dataset were employed for emotion identification. The content summary of preprocessed DEAP database is shown in Table 2. The preprocessed signals in DEAP dataset were down-sampled to 128Hz and the length of each trial was 63 seconds, including 3s baseline signal and 60s watching recodings. In this paper, after removing the 3s baseline signal, each 60s trial was segmented into ten 6s samples without overlapping. Finally, each subject presented with a total of 400 (40 trials × 10 samples) samples and the samples of five subjects (s01, s02, s03, s04, s05) were chosen to verify the effectiveness of the proposed methods. The sample distribution is shown in Table 3.

2.2. Emotion Labeling Schemes

In this research, two common types of emotion models were taken into consideration. One is to divide emotions into discrete categories. Ekman [31] regarded emotions as discrete and physiology related. He diveded emotions into six basic emotions with happy, sad, anger, fear, surprise and disgust and the other emotions were viewed as the production of reaction and combination of the six basic emotions. Izard [32] presented ten basic emotions—joy, sadness, angry, surprise, fear, interest, shyness, guilt, contempt and disgust. He suggested that each basic emotion corresponded to a simple brain circuit and there was no complex cognitive component involved. When using Augsburg Dataset, this discrete emotion model was employed to classify four different emotions, including joy, anger, sadness and pleasure, because the experimental emotional states were desigend with discrete categories. The other emotion model is based on the 2D arousal-valence emotion description model. In Russell’s theroy of emotion model [33], the emotional states are characterized by two dimension, valence and arousal. The emotion can be mapped to a plane with valence as the horizontal and arousal as the vertical axes. Arousal map emotions ranging from passtive to active while valence rangesfrom negative to positive. In this study, the DEAP dataset was divided into four-dimensional emotion groups according to Reference [30], including HVHA, HVLA, LVLA and LVHA based on the levels of valence and arousal, as shown in Figure 1. While valence > 5 and arousal > 5, it belongs to HVHA; while valence > 5 and arousal ≤ 5, it belongs toHVLA; while valence ≤ 5 and arousal > 5, it belongs to LVHA; while valence ≤ 5 and arousal ≤ 5, it belongs to LVLA.

2.3. Feature Extraction

Prior to using the classifier for sample identification, features which can represent the emotional states usually need be extracted. The traditional linear analysis methods can only achieve some superficial information, therefore numerous features need be extracted from physiological signals to display the characteristics of signals, which reduces the efficiency of emotion classification. In addition, compared with periodic signals, emotion-related physiological signals are nonlinear time-varying signals. Therefore, in this investigation, four nonlinear features are employed to represent the emotional states from physiological signals. The four selected nonlinear features, approximate entropy (ApEn), sample entropy (SaEn), fuzzy entropy (FuEn) and wavelet packet entropy (WpEn), can reveal the internal structure of the time series and provide a prior knowledge about the intrinsic similarity, deep-seated information and predictability of the signal.

2.3.1. Approximate Entropy

Approximate entropy (ApEn) is a nonlinear analysis method proposed by S.M. Pincus [34,35,36] to measure the complexity or irregularity of time series. Its physical meaning is to measure the probability of mode generation in a new state when the signal sequence changes in its dimensional space. The method is based on the theory of phase space reconstruction. The core idea of this theory is that embedding the signal into the phase space, when the embedding dimension of the phase space is increased from m to m + 1, predicting the probability of generating the new mode. The greater the probability of generating a new mode, the more complex the signal and the higher the corresponding value of ApEn. The computed steps of ApEn are as follows:
(1) For time series x i of finite length N ( 1 < i < N ) , x i is re-constructed into m-dimensional vector X m ( i ) :
{ X m ( i ) = { x ( i ) ,   x ( i + 1 ) ,   x ( i + 2 ) ,   ,   x ( i + m 1 ) } i = 1 ,   2 ,   , N m + 1 ,
where m is the embedded dimension.
(2) Calculate the distance d { X m ( i ) , X m ( j ) } between the components of X m ( i ) and X m ( j ) and define the maximum distance as D { X m ( i ) , X m ( j ) } :
D { X m ( i ) , X m ( j ) } = m a x { | x ( i + k ) x ( j + k ) | } .
(3) Calculate the probability of C i m ( r ) , which measures the regularity of x i , that is, the ratio of D { X m ( i ) , X m ( j ) } < r to the total number of N m + 1 :
C i m ( r ) = N m ( i ) N m + 1 ,
where r is the vector comparison distance and N m ( i ) represents the number of D { X m ( i ) , X m ( j ) } < r .
(4) Define Φ m ( r ) as:
Φ m ( r ) = i = 1 N m + 1 ln C i m ( r ) N m + 1 .
(5) When the dimension is increased to m + 1 , repeats the above steps to get Φ m + 1 ( r ) :
Φ m + 1 ( r ) = i = 1 N m ln C i m + 1 ( r ) N m .
(6) Finally, for fixed m, r and N, ApEn can be expressed as:
A p E n ( m , r , N ) = Φ m ( r ) Φ m + 1 ( r ) .
The value of ApEn is related to vector comparison distance r, the embedded dimension m and the data length N. According to the practice of Pincus et al. [34], it is suggested that a relatively stable estimate can be obtained by using shorter data; r ranges from 0.1std to 0.25std, where std is the standard deviation of the data; the value of m is from 2 to 5.

2.3.2. Sample Entropy

Sample entropy (SaEn), proposed by Richman et al [37], is used to measure self-similarity and complexity of time series signals, being an improved algorithm based on ApEn. In view of the shortcomings of ApEn, SaEn mainly made two improvements—(a) Eliminating the deviation of ApEn when matching with itself; (b) In order to avoid the condition of ln0 caused by eliminating self-matching, SaEn calculates the total number of matching boards before logarithmic operation and when calculates m + 1 dimensional statistics, the number of templates matching itself is used to calculate. The specific computing steps of SaEn are as follows:
(1) The first 3 steps are similar to the calculation of ApEn and define Q i m ( r ) as the probability that D { X m ( i ) , X m ( j ) } < r .
(2) After averaging all obtained Q i m ( r ) , gets the total number of template matches, Q m ( r ) :
Q m ( r ) = i = 1 N m ln Q i m ( r ) N m .
(3) When the dimension is increased to m + 1 , the total number of template matches is Q m + 1 ( r ) :
Q m + 1 ( r ) = i = 1 N m ln Q i m + 1 ( r ) N m .
(4) SaEn can be expressed as:
S a E n ( m , r , N ) = ln Q m + 1 ( r ) Q m ( r ) .

2.3.3. Fuzzy Entropy

Chen et al. [38,39] employed fuzzy theory to measure the complexity of signal sequence and proposed fuzzy entropy (FuEn) algorithm, which was successfully applied to the extraction and classification of EMG signals. Studies [40,41] show that FuEn algorithm has lower sensitivity and dependence on phase space dimension (m), similar tolerance limit (r) and the time series length (N) and contributes to the calculation efficiency. FuEn is computed as follows:
(1) The first 2 steps are similar to the calculation of ApEn. Define the maximum distance as D i j m :
D i j m = m a x { | x ( i + k ) x ( j + k ) | } .
(2) Imports fuzzy membership function A i j m :
A i j m = exp [ ln ( 2 ) ( D i j m r ) 2 ] .
(3) Define C i m ( r ) as:
C i m ( r ) = j = 1 , j i N m + 1 A i j m N m .
(4) Define Φ m ( r ) as:
Φ m ( r ) = i = 1 N m + 1 ln C i m ( r ) N m + 1 .
(5) When the dimension is increased to m + 1, repeats the above steps and finally FuEn is:
F u E n ( m , r , N ) = ln Φ m ( r ) ln Φ m + 1 ( r ) .

2.3.4. Wavelet Packet Entropy

Wavelet packet entropy (WpEn) is an algorithm that combines wavelet packet transform with information entropy, which takes the advantages of wavelet packet in accurately describing signals of different frequency bands and the information measurement of non-stationary signals based on information entropy. WpEn reflects the spectrum energy distribution of signals in different frequency band and can quantitatively describe the order or disorder degree of information distribution [42]. The specific algorithm steps of WpEn are as follows:
(1) The raw signal is decomposed into different signal components of different frequency bands by wavelet packet, the energy Ei,j for each frequency range in each time window can be computed as:
E i , j ( t ) = k = 1 L i , j ( x i , k ) 2 ,
where i denotes the number of layers of wavelet packet decomposition; j denotes the j-th frequency band; k is the summation index; Li, j denotes the coefficient energy.
(2) Total energy E t o t a l of the signal in each time window is calculated as:
E t o t a l = i = 1 N E i , j .
(3) According to Shannon’s information entropy theory and the definition of wavelet packet energy, WpEn is defined as follows:
W p E n = i p i ln ( p i ) ,
where P i = E i E t o t a l is computed as the ratio between the energy of each level.

2.3.5. Multimodal Feature Fusion

In this research, physiological signals are employed for emotion identification. The primary purpose of feature fusion is to improve the classification results by exploiting the complementary nature of different features [43,44,45]. Each of the four nonlinear features mentioned above can reflect the state characteristics of the signal to some extent from a certain point of view. Therefore, the early fusion of four nonlinear features combining as a single representation can reflect emotional states more effectively. Meanwhile, the feature fusion of multimodal physiological signals will be more comprehensive to represent emotional states from multiple perspectives. In this investigation, in order to represent emotional states effectively, we propose the feature representation method of multi-nonlinear feature integration and multi-channel information feature fusion. The four nonlinear features, are employed and extracted from every physiological signal, then the extracted features of various physiological signals are fused in order to represent the emotional states more effectively.

2.4. Team-Collaboration Identification Strategy Based on SVM-DT-ELM

Some traditional classifiers, such as support vector machine (SVM), extreme learning machine (ELM) and decision tree (DT), have been well applied in emotion recognition, while each classifier has his own shortcomings inevitably. For instance, SVM has better diagnostic performance under small sample conditions, poor performance for large sample conditions and emerge multiple categories with the same number of votes when voting. The initial input parameters of ELM are generated randomly, which requires a large number of training samples and cannot guarantee the optimal parameters. DT is inconsistent with the data of different samples and the information gain tend to those features with more values. In view of the requirements for accuracy and reliability of emotion recognition system and the uncertainty caused by a single classifier, a team-collaboration identification strategy based on the fusion of SVM, DT and ELM, is proposed, which exerts the function of collaborative diagnosis with multiple classifiers, thus eliminating the uncertainty brought by a single classifier and improving the recognition accuracy.

2.4.1. Support Vector Machine

Support vector machine (SVM) is a machine learning method based on the principle of structural risk minimization in statistical learning theory, which seeks the best performance between model complexity and learning ability to achieve the best generalization ability based on limited sample information [46,47,48]. The core idea is to realize nonlinear classification or regression fitting by mapping nonlinear classification or regression fitting problems into high-dimensional space by using kernel function to obtain the better classification or regression result. When making a decision in a classification problem, the voting method is usually adopted and the category with the most votes is the class to which the sample belongs.
For a data set { ( x 1 , y 1 ) , ( x 2 , y 2 ) ,   ,   ( x k , y k ) } of two classes with k as the number of samples, where x i R n represents the sample; y i { + 1 , 1 } is the class label;   i = 1 , 2 , , k is the training sample number. SVM seeks an optimal hyper-plane in the n dimensional data feature space by constructing the following function:
{ min 1 2 w T w + C i = 1 n ξ i y i ( w T ϕ ( x i ) + b ) 1 ξ i ξ i 0 , i = 1 , 2 , 3 , , n ,
where ϕ is a mapping function from low-dimensional space to high-dimensional space; ξ i is a slack variable to ensure the correctness of the classification in the case of inseparable samples; C is a penalty factor and a larger C indicates a greater penalty for misclassification; w and b are the weight vector and classification threshold of the decision function f ( x ) = ( w · x ) + b ; xi is the input vector and yi is the output vector.
Introducing the Lagrange function to get the dual optimization problem:
{ min 1 2 α T Q α e T α y T α = 0 0 α i C , i = 1 , 2 , 3 , , n ,
where e = ( 1 , 1 , , 1 ) T is the column vector; Q is the semi-positive definite matrix of n × n ; Q i j = y i y j K ( x i , x j ) , K ( x i , x j ) = ϕ ( x i ) T ϕ ( x j ) T is the kernel function; α i is the Lagrange multiplier; y is the sample label vector; α is the Lagrange multiplier vector.
Computing Equation (19), the optimal solution is:
w = i = 1 n y i α i ϕ ( x i ) .
The optimal hyper-plane decision function is:
f ( x ) = sgn ( i = 1 n y i α i K ( x i , x ) + b ) .
SVM can be extended to multi-classification problems by constructing multiple SVM two-class classifiers which include direct method, one-to-one and one-to-rest. Among them, the one-to-one method is used to classify the k classes of sample data by constructing k(k –1)/2 binary classifiers, which has a fast solving speed and is widely used in practice. The classification principle is the “voting mechanism,” that is, each classifier votes for its preference and the final result is based on the category with the most votes. This method can be expressed as:
{ min [ 1 2 ( x i j ) T w i j + C i = 1 n ξ t i j ] ( w i j ) T ϕ ( x i ) + b i j 1 ξ t i j ( y t = i ) ( w i j ) T ϕ ( x i ) + b i j 1 + ξ t i j ( y t = j ) ξ t i j 0 ,
where w i j and b i j are the weight vector and threshold obtained when designing the two-class classifier for the i-th sample and the j-th sample respectively; ξ t i j is the slack variable; x i j is the training sample vector; y t is the sample label; S is the sum of the i-th class samples and j-th class samples.

2.4.2. Decision Tree

Decision tree (DT) classifier is an instance-based inductive learning algorithm that uses inductive algorithm to generate readable decision trees and rules and then uses the decision tree to classify new data [49]. DT is an inverted tree structure similar to the flow chart, which mainly focuses on the two core problems of growth and pruning. The structure diagram of DT is shown in Figure 2. The knowledge acquired by DT is a formal representation of the tree, including the regression tree and the classification tree. The results of classification or prediction are reflected in the leaf nodes of DT. The average value of the output variables is the prediction result in the samples contained in the leaf nodes of the regression tree, while the mode of the output variable is the classification result in the samples contained in the leaf nodes of the classification tree.
Each none-leaf node in the figure represents the input attribute of the training data set, attribute value represents the value corresponding to the attribute and the leaf node means the value of the target category attribute. Yes and No represent positive and negative examples respectively.
DT classifier is computed as follows:
Input—training set D, feature set A and threshold ε ;
(1) If all samples in D belong to the same category of Ck, then T is a single node tree and Ck is used as the class of the node and returns T.
(2) If A is not an empty set, then T is a single node tree and the class Ck with the largest number in D is taken as the class of the node and returns T.
(3) Otherwise, calculate the information gain ratio (GA) of each feature in A according to Equation (23) and select feature Ag with the largest:
{ G R = G S I G = E n t r o p y ( S ) = v V ( A ) | S v | | S | E n t r o p y ( S v ) S I = i = 1 c | S i | | S | log 2 S i S E n t r o p y ( S ) = i = 1 c p i log 2 p i ,
where p i is the proportion of the sample of the i-th attribute value in the subset; V(A) is the range of the attribute A; SV is the subset of D whose value is V on the attribute A; E n t r o p y ( S ) is the entropy of D relative to C states.
(4) If Ag is less than ε , then T is a single node tree and the class Ck with the largest number in D is taken as the class of the node and T is returned;
(5) Otherwise, for each possible value ai of Ag, divide D into several non-empty subsets Di according to Ag = ai, mark the class with the largest number in Di as a mark to construct sub nodes and form a tree T by the node and return T.
(6) For node i, Di is used as training set and A-{Ag} as feature set. Step (1) to (5) are called recursively to get subtree Ti and return Ti.

2.4.3. Extreme Learning Machine

In order to improve the traditional learning algorithms (such as back propagation neural network), which easily fall into local minimum, slow model training speed and difficulty in adjusting learning rate, reference [50] proposed the Extreme Learning Machine (ELM), which consisted of only an input layer, a hidden layer and an output layer. The brief network structure of the algorithm is shown in Figure 3. According the inputs, ELM randomly generates the connection weights and the threshold of hidden layer neurons between input layer and hidden layer, which need not be adjusted during the training process. Users do not need to know the hidden layer because Gaussian kernel was applied. The optimal solution can be obtained by setting the number of hidden layer neurons, which is related to the number of features. Best values for positive regularization coefficient and Gaussian kernel parameter were found empirically after several experiments.
Suppose there are N training samples ( x i , y i ) R N × R m ( i = 1 ,   2 ,   ,   N ) , where x i R N is the input and y i R m the output. The mathematical model of a standard single hidden layer feed forward neural network with M hidden layer nodes is:
i = 1 M β i g ( ω j x i + b j ) = o i , i = 1 , 2 , 3 , , N ,
where ω j is the input weight vector connecting the input neuron and the j-th hidden layer neuron; β j is the output weight vector connecting the j-th hidden layer neuron and the output neuron; o i is the actual output vector; b j is the bias of the hidden layer neurons; g ( x ) is the activation function of the hidden layer neurons.
If the model can approximate the output y j of the training sample with zero error, which means i = 1 N o i y i = 0 , then β j , ω j and b j make the following formula hold:
i = 1 M β i g ( ω j x i + b j ) = y i , i = 1 , 2 , 3 , , N .
Simplifies (25):
H β = Y ,
where H is called the output matrix of the hidden layer of the neural network.
When the activation function of the neuron is arbitrarily differentiable, the training error of the single hidden layer feed forward neural network can approach infinitely small positive number ε . At this time, the input weight vector ω j and the hidden layer offset b j can remain unchanged during training process and can also be randomly assigned. Therefore, the training process is equivalent to finding the least squares solution β ^ of the linear system:
H β ^ T = min β H β T .
The solution is β ^ = H + Y and H + is the Moore-Penrose generalized inverse of the hidden layer output matrix H .

2.4.4. Team-Collaboration Identification Strategy

As mentioned above, each classifier is based on a different method principle, so each classifier has its own advantages and disadvantages. For a sample, it may be easily misclassified by one classifier but easily identified by other classifiers. In order to reduce the limitation of a single classifier and improve the accuracy of recognition, a team-collaboration identification strategy model, which combines SVM, DT and ELM, is proposed. For this proposed team-collaboration identification strategy, the SVM model is regarded as a major decision expert and the DT and ELM are emloyed to provide the decision suggestions for the samples which are easily misclassified by SVM. The core idea of SVM-DT-ELM is that selecting the possibly misclassified samples for the SVM and then employing DT and ELM to conduct referral for these samples and finally confirming the emotion class of the sample according to the designed decision-making mchanism. Main procedures of the suggested SVM-DT-ELM algorithm are as follows:
(1) Firstly, the training sets are used to train SVM, DT, ELM classification model respectively. In this research, during training SVM model, the radial basis function (RBF kernel) is selected and the grid search method is used to optimize the SVM model parameters to achieve better performances.
(2) Selecting samples that may be misclassified by SVM. According to the training model of SVM and the self-classification accuracy, selection conditions where the possibly misclassified samples belong to are determined. Analysis shows that samples distributed near the support vector or with the same number of votes during voting process are easily misclassified when SVM is used for classification. The SVM supports one-versus-one multiclassification. If k is the number of classes, k(k-1)/2 models will be generated, each model involves only two class. Focusing on above problems, set the following conditions:
a) When using SVM for classification, if the highest number of votes are equal during the voting process, the sample will be regarded as a possibly misclassified sample and the sample will be picked out for referral.
b) If the test sample satisfies the following condition after input to SVM, this sample is selected for re-diagnosis.
{ u min > v min h max u min > s max v min t 1 < h max v min u min s max < t 2 ,
where umin represents the smallest absolute value in the case of three votes; hmax is the largest absolute value in the case of three votes; vmin represents the smallest absolute value in the case of two votes; smax represents the largest absolute value in the case of two votes. The t 1 and t 2 are conditional parameters, which are determined by the performances of the trained SVM model, in this research, t 1 = 1.5 and t 2 = 3.0 .
(3) Decision principles. When applying SVM-DT-ELM team-collaboration strategy for emotion identification, this paper follows the following principles:
ⅰ) If the test sample is classified by SVM with full votes and the sample is outside the set conditions, then DT and ELM are not employed for further consultation. The output emotional categories are based on the results of SVM.
ⅱ) Samples except those satisfying condition ⅰ) should be classified by DT and ELM. If any one of the results between DT and ELM are the same to SVM, the output emotional classes are based on the principle of minority obeying majority.
ⅲ) If the results of DT, ELM are different from the category of highest ranked vote of SVM and any one of the referral results between DT and ELM is consistent with the SVM’s second highest ranked vote, the final diagnosis category is based on the referral result.
(4) According to the principle of step (3), the emotional categories of the test samples are confirmed.
The flow chart of emotion recognition algorithm based on SVM-DT-ELM is shown in Figure 4.

3. Results and Discussions

The Augsburg Dataset and DEAP dataset were employed in order to fully verify the effectiveness of the proposed method. The performances of the proposed method were analyzed and compared with the exisiting studies from multiple perspectives, which highlights the significant ability of the proposed methods to recognize emotions through peripheral physiological signals.

3.1. Experiment Environment

All implementations are performed using MATLAB (R2015b) running on Windows 10 Laptop PC with Intel(R)Core (TM) i7-8750H CPU @ 2.21GHz processor with 16 GB RAM. Table 4 shows the hardware and software for the experiments.

3.2. Procedure of Emotion Recognition

Figure 5 illustrates the architecture of emotion recognition from physiological signals. Firstly, the raw emotional physiological signals need to be preprocessed. Then, the nonlinear features are extracted from four types of signals. Next, the extracted multimodal features are fused, the training samples are used to train the classifiers and test samples are classified with the proposed team-collaboration identification strategy.

3.3. Model Performance Evaluation Method

In order to quantify the performance of the proposed approach, the Hold-Out method is adopted to train and test the constructed classifier model, where the dataset is divided into two mutually exclusive sets, one is the training set and the other is the testing set. The 60%–80% of the dataset are usually randomly selected for training and the remaining are used for testing. Generally, the experiments need be repeated several times with random division and the average value is as the final result. Besides, the standard evaluation metric, accuracy (Acc), was used as a measurement to evaluate the performance of the above classification models and the values under different conditions are reported as mean ± standard. The calculation formula of average recognition rate and standard deviation is shown in (29) and (30):
A c c = i N A c c i N
σ = i N ( A c c A c c i ) 2 N 1 ,
where Acci represents the recognition accuracy of the i-th experiment, that is, the number of correctly classified samples divided by the total number of samples; N represents the number of experiments.
According to the size of the Augsburg dataset, 80% of the samples are used for training, while the remaining samples are used for testing. This procedure is repeated 10 times to ensure the results more reasonable. The accuracy values of 10 times was averaged as the final classification performance. When the experiments are conducted on the DEAP dataset, for each subject, 75% of the samples are used for training, while the remaining 25% are used for testing. Then, for each subject, the average classification performance is computed over the ten train-test repetition.

3.4. Emotion Classification in Augsburg Dataset

In order to fully verify the effectiveness of the proposed method, the experimental results of emotion recognition are presented from three dimensions—(1) To show that the proposed feature level fusion is effective, the recognition results obtained by using features of a single signal are compared with that based on feature fusion; (2) To present the effectiveness of team-collaboration identification strategy, the recognition results obtained by using a single classifier to classify the samples with feature fusion are compared with the results based on the proposed team-collaboration identification strategy; (3) To further demonstrate the advantage of the proposed method, compare the results of this paper with that of related researches.

3.4.1. Feature Level Fusion

Prior to classification of samples, four nonlinear features, ApEn, SaEn, FuEn and WpEn, were extracted from each physiological signal, including ECG, EMG, RSP and SC. Then the extracted features were fused and fed into classifiers for classification.
In order to show that the feature level fusion mechanism in this paper is effective, the recognition results of single physiological signal were compared with the recognition results of multimodal signals (the classifier uses the SVM, DT and ELM respectively). Table 5 shows the comparison of the 10-times average recognition rate of four types of emotional states between single signal features and feature fusion with multitype signals (Randomly selecting 500 samples each time). As can be seen the recognition accuracy can be effectively improved by the fusion at the feature level with each classifier. The high accuracy identification with fused features combining different classifiers clearly indicate that the features obtained from different signals containing complementary or supplementary information.

3.4.2. Team-Collaboration Identification Strategy

In order to further improve the emotion recognition accuracy, the proposed SVM-DT-ELM team-collaboration identification strategy, is employed for classification. Firstly, 80 percent samples were randomly selected as training samples to establish the SVM, DT and ELM classification model respectively. According to the training model of SVM, selection conditions where the samples possibly misclassified were determined. The rest samples were used to test the performance of the proposed classification model and the experiments were repeated 10 times (Randomly selecting 500 samples each time). In addition, we also compared the classification accuracy of using SVM, DT and ELM with the recognition results of proposed method, as shown in Table 6.
Table 6 presents the accuracies of 10 experiments of using SVM, DT, ELM classifier and proposed team-collaboration identification strategy for four emotions when the features obtained from various physiological signals are taken in a fused manner. Compared with the accuracy of SVM classifier, the accuracy rate of the proposed SVM-DT-ELM reaches 98.6%, 3.1% higher than the result of SVM (95.5%). The proposed team-collaboration identification strategy presents better performances with 8.1% and 9.2% improved accuracy than DT and ELM, respectively. The experimental results demonstrate that the team-collaboration identification strategy is able to further improve the recognition accuracy and make the recognition results more reliable.

3.4.3. Comparison with Existing Methods

Table 7 provides an overview of the accuracy rate comparisons of existing studies and the proposed method with Augsburg Dataset. Some studies for emotion recognition reported their performances based on valence and arousal measures only, which are not included here for comparison. As we can see, the recognition accuracy of this paper is increased to 98.6% with 16 features based on the proposed nonlinear features fusion and team-collaboration identification strategy. From the perspective of feature dimension and recognition accuracy, the extracted features have better performance than other features [29,51,52,53,54] and effectively reduce the feature dimension and improve emotion recognition accuracy. In addition, the proposed team-collaboration identification strategy can integrate the advantages of other classifications to effectively improve the classification accuracy.

3.5. Emotion Classification in DEAP Dataset

In this section, we present the results obtained from DEAP dataset to demonstrate the effectiveness of proposed methods.

3.5.1. Feature Level Fusion

In order to perform a more reliable classification process, we constructed a training set and a test set for each subject (s01, s02, s03, s04, s05). The number of training set is 300 and the test set is 100 for each subject. The experiments were carried out ten times and the accuracy rates of four dimensional emotions identification using single signals and multimodal signals with SVM, DT and ELM are shown in Table 8, Table 9 and Table 10, respectively.
It can be seen in Table 8, Table 9 and Table 10 that the classification performances of RSP and EMG are good on the whole, presenting their certain advantages on detecting different emotions. These two physiological signals also perform well in Augsburg dataset, which is where the two databases are consistent. Compared with the other three signals, the classification accuracy of GSR is not good all in all. The results of BVP is more satisfying. In addition, it can be seen that the results of multi sensors fusion present much improvement comparing with the results of single signal for each subject. Performance analysis shows that multi-type information fusion can effectively improve the accuracy of emotion identification, which is consistent with the experimental results with the Augsburg dataset.

3.5.2. Team-Collaboration Identification Strategy

In this research, the team-collaboration identification strategy is proposed to avoid the limitations of a single classification method. The results of ten times experiments are shown in Table 11, comparing between SVM, DT, ELM and the proposed method for the five subjects individually. Due to individual differences, the accuracy of emotion recognition varies among different subjects. However, as can be seen, the results illustrate the improvement in the classification accuracy for each subject after using the proposed classification model. Taking the subject s01 as an example, when employing the proposed strategy, the average accuracy improved by 6% than the SVM classfication method, from 73.5% to 79.5%. And there were 19.2% and 18% improved accuracy than the rusults of DT and ELM, respectively. Overall, the average identification accuracies are 70.4% and 76.46% for SVM and proposed strategy, respectively, improving with 6.06%. Meanwhile, when the proposed method was employed, there were smaller mean square deviation, which present the better stability of the team-collaboration identification strategy than the other methds. The compared results for each subject with SVM, DT, ELM and proposed method are demonstrated in Figure 6.

3.5.3. Comparison with Existing Methods

In recent years, various studies of emotion recognition have been conducted on DEAP dataset. Generally, The research work mainly focuses on two dimensional classification (HA/LA, HV/LV) or four dimensional classification (HVHA, HVLA, LVLA, LVHA).
The previous studies [55,56,57,58,59] based on two dimensional emotion classes using DEAP dataset are shown in Table 12. As can be seen from Table 12, the methods adopted [47,48,49,50,51] are able to identify the two types of emotions to a certain extent. However, binary emotion model is difficult to fully express emotional states. Therefore, this paper takes four dimensional emotion classes (HVHA, HVLA, LVLA, LVHA) into consideration.
In this investigation, the emotional states were identified from subject-dependent perspective. The training and testing were performed on the same subject, which was same as the Ref [58,60,61]. An overview comparison of accuracy rate based on four dimensional emotion classes with different methods is shown in Table 13. In this paper, the identification accuracy of the four types of emotions was 76.46%. It can been seen that the proposed method with fused nonlinear features and team-collaboration identification strategy presents better identification performance than the existing methods in References [58,60,61]. Although only five experimental subjecs (s01~s05) were employed for analysis in this paper, the analysis results are not affected, because the subjects were selected sequentially, not deliberately. In addition, the studies of References [58,60,61] have been carried out based on EEG signals, while this paper is based on physiological signals. The experimental results show that as long as the method is proper, the use of physiological signals can also achieve good emotion recognition. In this investigation, the effectiveness of proposed methods of feature representation and classification were verified by comprehensive analysis with Augsburg dataset and DEAP dataset. The experimental analysis results of two datasets are consistant that the proposed methods present good performances of emotion identification, which indicates the universality of the proposed methods.

3.6. Discussions

In this investigation, the method of nonlinear features extraction and multi-signal feature fusion was proposed to effectively characterize emotional responses and team-collaboration identification strategy was suggested to improve the accuracy of emotion recognition. The effectiveness of the proposed feature extraction fusion method and classification method were verified by single factor analysis comparision on Augsburg and DEAP datasets respectively. The results before and after feature fusion were compared to verify the effectiveness of the proposed feature fusion method. The effectiveness of the proposed classification method was verified by comparing the results of the single classification method with that of the team-collaboration identification strategy. Meanwhile, the validity of the proposed methods were fully verified by comparing with the results of existing research methods based on the same database. Whether it was the comparison of the effect before and after the fusion or the comparison with the methods of other researchers, the methods proposed in this article have shown superior performances.
Various methods of feature extraction and classification [51,52,53,54,55,56,57,58,59,60,61] were proposed and employed by researchers, in terms of using physiological signals to identify emotions. No matter which method was adopted, the purpose was to improve the accuracy of emotion identification. Recently, deep learning methodologies have become popular to analyse physiological signals and employ to realize emotion classification [5]. Kwon et al [59] employed deep learning method to conduct two dimensional emotion identification on DEAP dataset. At present, there is no application of deep learning methods for the same research goal as this article on the Augsburg and DEAP datasets. Hence, the results of this paper were not compared with that of deep learning method. In this research, at present, we have studied the emotion recognition of a single person using physiological signals. In the future work, further investigation is needed on how to extract more discriminative features to make cross-subject emotion classification and how to construct and optimize classification models with higher accuracy for emotion recognition and how to effectively adopt deep learning method to identify emotions with peripheral neurophysiological signals and EEG signals.

4. Conclusions

In this paper, the methods of emotional state identification based on physiological signals were investigated. In order to represent emotional states effectively, the method of nonlinear features extraction and multi-signal feature fusion was suggested. Meanwhile, team-collaboration identification strategy was proposed for avoiding the limitations of a single classifier. The four nonlinear features, namely ApEn, SaEn, FuEn and WpEn were employed and extracted from each physiological signal, then the extracted features of different physiological signals were fused, for example fusion of ECG, EMG, RSP, SC with Augsburg dataset and EMG, RSP, BVP, GSR with DEAP dataset. Nonlinear features can represent emotion characteristics in depth and multi-signal feature fusion can represent emotion from multiple perspectives. Hence, the proposed feature represent method is more effective. In order to improve the accuracy of emotion identification, team-collaboration identification strategy was proposed. According to the performances and characteristics of SVM, DT and ELM classifiers, the team-collaboration decision-making mechanism was developed, which according to the difficulty of sample classification, single classification decision and team-collaboration identification can be adopted respectively. The proposed strategy can choose the appropriate decision method according to the characteristics of the samples and can effectively integrate the advantages of other classifiers to avoid the limitations of single classifier. The Augsburg dataset and DEAP dataset were employed to verify the validity and universality of the proposed method. The experimental results from the two public databases uniformly indicated that the proposed method combining fused nonlinear features and team-collaboration identification strategy has better performances than the existing methods.

Author Contributions

Formal analysis, L.P.; Software, L.P. and Z.Y.; Validation, L.P. and A.S.; Writing – original draft, L.P. and Z.Y.; and Writing – review and editing, L.P. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by National Natural Science Foundation of China: 61773078; Open Foundation of Remote Measurement and Control Key Lab of Jiangsu Province: YCCK201303; Industrial Technology Project Foundation of Chang Zhou Government: CE20175040.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kwakkel, G.; Kollen, B.J.; Krebs, H.I. Effects of Robot-Assisted Therapy on Upper Limb Recovery After Stroke: A Systematic Review. Neurorehabilit. Neural Repair 2008, 22, 111–121. [Google Scholar] [CrossRef] [PubMed]
  2. Long, X.; Fonseca, P.; Foussier, J.; Haakma, R.; Aarts, R. Promoting interactions between humans and robots using robotic emotional behavior. IEEE Trans. Cyber. 2016, 46, 2911–2923. [Google Scholar] [CrossRef]
  3. Calvo, R.A.; D’Mello, S.K. Affect Detection: An Interdisciplinary Review of Models, Methods, and Their Applications. IEEE Trans. Affect. Comput. 2010, 1, 18–37. [Google Scholar] [CrossRef]
  4. Kerkeni, L.; Serrestou, Y.; Raoof, K.; Mbarki, M.; Mahjoub, M.A.; Cleder, C. Automatic speech emotion recognition using an optimal combination of features based on EMD-TKEO. Speech Commun. 2019, 114, 22–35. [Google Scholar] [CrossRef]
  5. Hassan, M.M.; Alam, G.R.; Uddin, Z.; Huda, S.; Almogren, A.; Fortino, G. Human emotion recognition using deep belief network architecture. Inf. Fusion 2019, 51, 10–18. [Google Scholar] [CrossRef]
  6. Hu, M.; Wang, H.; Wang, X.; Yang, J.; Wang, R. Video facial emotion recognition based on local enhanced motion history image and CNN-CTSLSTM networks. J. Vis. Commun. Image Represent. 2019, 59, 176–185. [Google Scholar] [CrossRef]
  7. Picard, R.W.; Vyzas, E.; Healey, J. Toward machine emotional intelligence: Analysis of affective physiological state. IEEE Trans. Pattern Anal. Mach. Intell. 2001, 23, 1175–1191. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, J.; Yin, Z.; Chen, P.; Nichele, S. Emotion recognition using multi-modal data and machine learning techniques: A tutorial and review. Inf. Fusion 2020, 59, 103–126. [Google Scholar] [CrossRef]
  9. Zheng, W.-L.; Lu, B.-L. Investigating Critical Frequency Bands and Channels for EEG-Based Emotion Recognition with Deep Neural Networks. IEEE Trans. Auton. Ment. Dev. 2015, 7, 162–175. [Google Scholar] [CrossRef]
  10. Chen, P.; Zhang, J. Performance Comparison of Machine Learning Algorithms for EEG-Signal-Based Emotion Recognition. In Proceedings of the International Conference on Artificial Neural Networks (ICANN2017), Alghero, Italy, 11–14 September 2017; pp. 208–216. [Google Scholar]
  11. Yang, Z. Emotion Recognition Based on Nonlinear Features of Skin Conductance Response. J. Inf. Comput. Sci. 2013, 10, 3877–3887. [Google Scholar] [CrossRef]
  12. Song, T.; Zheng, W.; Lu, C.; Zong, Y.; Zhang, X.; Cui, Z. MPED: A Multi-Modal Physiological Emotion Database for Discrete Emotion Recognition. IEEE Access 2019, 7, 12177–12191. [Google Scholar] [CrossRef]
  13. Campbell, E.; Phinyomark, A.; Scheme, E. Feature Extraction and Selection for Pain Recognition Using Peripheral Physiological Signals. Front. Mol. Neurosci. 2019, 13, 437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cheng, X.F.; Wang, Y.; Dai, S.C.; Zhao, P.J.; Liu, Q.F. Heart sound signals can be used for emotion recognition. Sci. Rep. 2019, 9, 6486. [Google Scholar] [CrossRef]
  15. Dissanayake, T.; Rajapaksha, Y.; Ragel, R.; Nawinne, I. An Ensemble Learning Approach for Electrocardiogram Sensor Based Human Emotion Recognition. Sensors 2019, 19, 4495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Levenson, R.W. The Autonomic Nervous System and Emotion. Emot. Rev. 2014, 6, 100–112. [Google Scholar] [CrossRef]
  17. Ali, M.; Al Machot, F.; Mosa, A.H.; Jdeed, M.; Al Machot, E.; Kyamakya, K. A Globally Generalized Emotion Recognition System Involving Different Physiological Signals. Sensors 2018, 18, 1905. [Google Scholar] [CrossRef] [Green Version]
  18. Abadi, M.K.; Subramanian, R.; Kia, S.M.; Avesani, P.; Patras, I.; Sebe, N. DECAF: MEG-Based Multimodal Database for Decoding Affective Physiological Responses. IEEE Trans. Affect. Comput. 2015, 6, 209–222. [Google Scholar] [CrossRef]
  19. Choi, J.; Ahmed, B.; Gutierrez-Osuna, R. Development and Evaluation of an Ambulatory Stress Monitor Based on Wearable Sensors. IEEE Trans. Inf. Technol. Biomed. 2012, 16, 279–286. [Google Scholar] [CrossRef] [Green Version]
  20. Kappeler-Setz, C.; Arnrich, B.; Schumm, J.; La Marca, R.; Troster, G.; Ehlert, U. Discriminating Stress From Cognitive Load Using a Wearable EDA Device. IEEE Trans. Inf. Technol. Biomed. 2010, 14, 410–417. [Google Scholar] [CrossRef]
  21. Greco, A.; Valenza, G.; Lanata, A.; Scilingo, E.P.; Citi, L. cvxEDA: A Convex Optimization Approach to Electrodermal Activity Processing. IEEE Trans. Biomed. Eng. 2016, 63, 797–804. [Google Scholar] [CrossRef] [Green Version]
  22. He, L.; Lech, M.; Zhang, J.; Ren, X.; Deng, L. Study of wavelet packet energy entropy for emotion classification in speech and glottal signals. In Proceedings of the Fifth International Conference on Digital Image Processing, Beijing, China, 21–22 April 2013; Volume 8878. [Google Scholar] [CrossRef]
  23. Wang, X.-W.; Nie, D.; Lu, B.-L. Emotional state classification from EEG data using machine learning approach. Neurocomputing 2014, 129, 94–106. [Google Scholar] [CrossRef]
  24. Jie, X.; Cao, R.; Li, L. Emotion recognition based on the sample entropy of EEG. BioMed. Mater. Eng. 2014, 24, 1185–1192. [Google Scholar] [CrossRef] [PubMed]
  25. Hosseini, S.A.; Naghibi-Sistani, M.-B. Emotion recognition method using entropy analysis of EEG signals. Int. J. Image Graph. Signal Process. 2011, 3, 30–36. [Google Scholar] [CrossRef] [Green Version]
  26. Vayrynen, E.; Kortelainen, J.; Seppanen, T. Classifier-based learning of nonlinear feature manifold for visualization of emotional speech prosody. IEEE Trans. Affect. Comput. 2013, 4, 47–56. [Google Scholar] [CrossRef]
  27. Li, X.; Cai, E.; Tian, Y.; Sun, X.; Fan, M. An improved electroencephalogram feature extraction algorithm and its application in emotion recognition. J. Biomed. Eng. 2017, 34, 510–517. [Google Scholar] [CrossRef]
  28. Mohammadi, Z.; Frounchi, J.; Amiri, M. Wavelet-based emotion recognition system using EEG signal. Neural Comput. Appl. 2017, 28, 1985–1990. [Google Scholar] [CrossRef]
  29. Wagner, J.; Kim, J.; André, E. From Physiological Signals to Emotions: Implementing and Comparing Selected Methods for Feature Extraction and Classification. In Proceedings of the 2005 IEEE International Conference on Multimedia and Expo, Amsterdam, The Netherlands, 6 July 2005; pp. 940–943. [Google Scholar]
  30. Koelstra, S.; Muhl, C.; Soleymani, M.; Lee, J.-S.; Yazdani, A.; Ebrahimi, T.; Pun, T.; Nijholt, A.; Patras, I. DEAP: A Database for Emotion Analysis; Using Physiological Signals. IEEE Trans. Affect. Comput. 2012, 3, 18–31. [Google Scholar] [CrossRef] [Green Version]
  31. Ekman, P. An argument for basic emotions. Cogn. Emot. 1992, 6, 169–200. [Google Scholar] [CrossRef]
  32. Izard, C.E. Basic Emotions, Natural Kinds, Emotion Schemas, and a New Paradigm. Perspect. Psychol. Sci. 2007, 2, 260–280. [Google Scholar] [CrossRef]
  33. Russell, J.A. A circumplex model of affect. J. Pers. Soc. Psychol. 1980, 39, 1161–1178. [Google Scholar] [CrossRef]
  34. Pincus, S.M. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Pincus, S. Approximate entropy (ApEn) as a complexity measure. Chaos 1995, 5, 110–117. [Google Scholar] [CrossRef] [PubMed]
  36. Pincus, S.M. Approximate entropy as a measure of irregularity for psychiatric serial metrics. Bipolar Disord. 2006, 8, 430–440. [Google Scholar] [CrossRef] [PubMed]
  37. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [Green Version]
  38. Chen, W.; Wang, Z.; Xie, H.; Yu, W. Characterization of Surface EMG Signal Based on Fuzzy Entropy. IEEE Trans. Neural Syst. Rehabil. Eng. 2007, 15, 266–272. [Google Scholar] [CrossRef]
  39. Chen, W.; Zhuang, J.; Yu, W.; Wang, Z. Measuring complexity using FuzzyEn, ApEn, and SampEn. Med Eng. Phys. 2009, 31, 61–68. [Google Scholar] [CrossRef]
  40. Chen, X.-J.; Li, Z.; Bai, B.-M.; Pan, W.; Chen, Q.-H. A New Complexity Metric of Chaotic Pseudorandom Sequences Based on Fuzzy Entropy. J. Electron. Inf. Technol. 2011, 33, 1198–1203. [Google Scholar] [CrossRef]
  41. Sun, K.H.; He, S.B.; Yin, L.Z.; Duo, L.K. Application of FuzzyEn algorithm to the analysis of complexity of chaotic sequence. Acta Phys. Sin. 2012, 61, 130507. [Google Scholar] [CrossRef]
  42. Cheng, B.; Liu, G.Y. Emotion recognition based on wavelet packet entropy of surface EMG signal. Comp. Eng. Appl. 2008, 44, 214–216. [Google Scholar] [CrossRef]
  43. Snoek, C.; Worring, M.; Smeulders, A.W.M. Early versus late fusion in semantic video analysis. In Proceedings of the 13th Annual ACM International Conference on Multimedia, Hilton, Singapore, 6–11 November 2005. [Google Scholar]
  44. Turk, M. Multimodal Human–Computer Interaction. Real-Time Vision for Human–Computer Interaction; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  45. Guironnet, M.; Pellerin, D.; Rombaut, M. Video Classification based on low-level feature fusion model. In Proceedings of the European Signal Processing Conference, Antalya, Turkey, 4–8 September 2005. [Google Scholar]
  46. Vapnik, V.N. The Nature of Statistical Learning Theory; Springer: Berlin/Heidelberg, Germany, 1995; pp. 123–179. [Google Scholar]
  47. Vapnik, V. Statistical Learning Theory; Wiley: New York, NY, USA, 1998. [Google Scholar]
  48. Chang, C.-C.; Lin, C.-J. LIBSVM: A library for support vector machines. ACM Trans. Intell. Syst. Technol. 2011, 2. [Google Scholar] [CrossRef]
  49. Sorel, L.; Viaud, V.; Durand, P.; Walter, C. Modeling spatio-temporal crop allocation patterns by a stochastic decision tree method, considering agronomic driving factors. Agric. Syst. 2010, 103, 647–655. [Google Scholar] [CrossRef]
  50. Huang, G.-B.; Zhu, Q.-Y.; Siew, C.-K. Extreme learning machine: Theory and applications. Neurocomputing 2006, 70, 489–501. [Google Scholar] [CrossRef]
  51. Guendil, Z.; Lachiri, Z.; Maaoui, C.; Pruski, A. Emotion recognition from physiological signals using fusion of wavelet based features. In Proceedings of the 2015 7th International Conference on Modelling, Identification and Control (ICMIC), Sousse, Tunisia, 18–20 December 2015; pp. 1–6. [Google Scholar]
  52. Wong, W.M.; Tan, A.W.; Loo, C.; Liew, W.S. PSO optimization of synergetic neural classifier for multichannel emotion recognition. In Proceedings of the 2010 2nd World Congress on Nature and Biologically Inspired Computing (NaBIC), Fukuoka, Japan, 15–17 December 2010; pp. 316–321. [Google Scholar] [CrossRef]
  53. Zong, C.; Chetouani, M. Hilbert-Huang transform based physiological signals analysis for emotion recognition. In Proceedings of the 2009 IEEE International Symposium on Signal Processing and Information Technology (ISSPIT), Ajman, UAE, 14–17 December 2009; pp. 334–339. [Google Scholar] [CrossRef]
  54. Gong, P.; Ma, H.T.; Wang, Y. Emotion recognition based on the multiple physiological signals. In Proceedings of the 2016 IEEE International Conference on Real-time Computing and Robotics (RCAR), Angkor Wat, Cambodia, 6–10 June 2016; pp. 140–143. [Google Scholar]
  55. Chen, J.; Hu, B.; Moore, P.; Zhang, X.; Ma, X. Electroencephalogram-based emotion assessment system using ontology and data mining techniques. Appl. Soft Comput. 2015, 30, 663–674. [Google Scholar] [CrossRef]
  56. Zhuang, N.; Zeng, Y.; Tong, L.; Zhang, C.; Zhang, H.; Yan, B. Emotion Recognition from EEG Signals Using Multidimensional Information in EMD Domain. BioMed Res. Int. 2017, 2017, 8317357. [Google Scholar] [CrossRef]
  57. Yin, Z.; Zhao, M.; Wang, Y.; Yang, J.; Zhang, J. Recognition of emotions using multimodal physiological signals and an ensemble deep learning model. Comput. Methods Programs Biomed. 2017, 140, 93–110. [Google Scholar] [CrossRef]
  58. Alazrai, R.; Homoud, R.; Alwanni, H.; Daoud, M.I. EEG-Based Emotion Recognition Using Quadratic Time-Frequency Distribution. Sensors 2018, 18, 2739. [Google Scholar] [CrossRef] [Green Version]
  59. Kwon, Y.-H.; Shin, S.-B.; Kim, S.-D. Electroencephalography Based Fusion Two-Dimensional (2D)-Convolution Neural Networks (CNN) Model for Emotion Recognition System. Sensors 2018, 18, 1383. [Google Scholar] [CrossRef] [Green Version]
  60. Zubair, M.; Yoon, C. EEG Based Classification of Human Emotions Using Discrete Wavelet Transform. In Proceedings of the Conference on IT Convergence and Security 2017, Seoul, Korea, 25–28 September 2017; pp. 21–28. [Google Scholar]
  61. Zheng, W.L.; Zhu, J.Y.; Lu, B.L. Identifying stable patterns over time for emotion recognition from EEG. IEEE Trans. Affect. Comput. 2019, 10, 417–429. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Emotion model classification based on the scale of valence and arousal.
Figure 1. Emotion model classification based on the scale of valence and arousal.
Entropy 22 00511 g001
Figure 2. The structure diagram of decision tree (DT).
Figure 2. The structure diagram of decision tree (DT).
Entropy 22 00511 g002
Figure 3. The network structure of Extreme Learning Machine (ELM).
Figure 3. The network structure of Extreme Learning Machine (ELM).
Entropy 22 00511 g003
Figure 4. The flow chart of SVM-DT-ELM algorithm.
Figure 4. The flow chart of SVM-DT-ELM algorithm.
Entropy 22 00511 g004
Figure 5. The architecture of emotion recognition with fused features and team-collaboration identification strategy.
Figure 5. The architecture of emotion recognition with fused features and team-collaboration identification strategy.
Entropy 22 00511 g005
Figure 6. Average accuracy of each subject using SVM, DT, ELM and proposed method for four types of emotional identification.
Figure 6. Average accuracy of each subject using SVM, DT, ELM and proposed method for four types of emotional identification.
Entropy 22 00511 g006
Table 1. The content summary of Augsburg Database.
Table 1. The content summary of Augsburg Database.
Elicitation MaterialMusic
Emotional statesJoy, anger, sadness, pleasure
Number of subjects1
Collected signalsECG, EMG, RSP, SC
Length120 seconds
Sampling frequencyECG: 256 Hz; EMG, RSP and SC: 32 Hz;
Collected days25
Table 2. The preprocessed DEAP database content summary.
Table 2. The preprocessed DEAP database content summary.
Elicitation MaterialVideos
Emotion labelsArousal, valence
Number of subjects32
Collected signalsEEG, EOG, GSR, BVP, RSP, EMG, SKT
Length60 seconds
Sampling frequency128 Hz;
Rating valuesArousal: 1–9
Valence: 1–9
Table 3. The number of samples for four dimensional emotional states of each subject.
Table 3. The number of samples for four dimensional emotional states of each subject.
SubjectHVHAHVLALVLALVHATotal
s0113060100110400
s021606010080400
s031021011070400
s041204020040400
s0513011010060400
Total5504806103602000
Table 4. The hardware and framework specifications.
Table 4. The hardware and framework specifications.
CPUIntel Core i7-8750H
GPUNVIDIA GeForce GTX1050Ti 4GB
OSWindows 10
RAMDDR4 16GB
FrameworksMATLAB (R2015b)
Table 5. Comparison of the performances between single signal features and feature fusion of multitype signals.
Table 5. Comparison of the performances between single signal features and feature fusion of multitype signals.
Physiological SensorAcc*(%)
SVMDTELM
Single sensorECG65.7 ± 1.5562.1 ± 1.9258.4 ± 2.72
EMG72.1 ± 0.9760.1 ± 1.8462.3 ± 2.47
RSP66.4 ± 1.2266.2 ± 1.7559.7 ± 2.62
SC70.9 ± 1.4364.5 ± 1.5460.6 ± 2.33
Multi sensorsECG + EMG + RSP + SC95.5 ± 0.8590.5 ± 1.2789.4 ± 1.78
Table 6. An overview of the comparison of the classification accuracy using different classifiers.
Table 6. An overview of the comparison of the classification accuracy using different classifiers.
Number of ExperimentsClassification Methods (Acc*/%)
SVMDTELMSVM-DT-ELM
196918698
296918998
395928998
495899199
594908898
695909098
795939299
8968990100
996908899
1097909199
Acc* (%)95.5 ± 0.8590.5 ± 1.2789.4 ± 1.7898.6 ± 0.70
Table 7. Accuracy comparison of various studies.
Table 7. Accuracy comparison of various studies.
Classification MethodFeature DimensionAcc* (%)
LDF [29]3292.05
SVM [51]6495
PSO-SNC [52]3286
SVM [53]2876
C4.5 DT [54]15593
This paper1698.6
Table 8. The classification accuracy for each subject using single signal and multimodal signals with SVM.
Table 8. The classification accuracy for each subject using single signal and multimodal signals with SVM.
SubjectPhysiological Sensors
Single Sensor (Acc*/%)Multi Sensors (Acc*/%)
GSRRSPBVPEMGGSR + RSP + EMG + BVP
s0138.0 ± 1.4643.3 ± 2.7353.5 ± 3.0650.6 ± 1.5673.5 ± 2.07
s0234.0 ± 2.2553.1 ± 1.7343.2 ± 1.6752.1 ± 3.0365.1 ± 2.69
s0354.3 ± 1.5160.8 ± 2.7265.6 ± 1.8363.2 ± 2.3381.5 ± 1.35
s0448.6 ± 4.1352.3 ± 3.6756.5 ± 3.9956.1 ± 2.1762.7 ± 2.21
s0532.8 ± 2.5643.1 ± 3.8847.6 ± 3.6942.5 ± 3.0469.2 ± 2.70
Table 9. The classification accuracy for each subject using single signal and multimodal signals with DT.
Table 9. The classification accuracy for each subject using single signal and multimodal signals with DT.
SubjectPhysiological Sensors
Single Sensor (Acc*/%)Multi Sensors (Acc*/%)
GSRRSPBVPEMGGSR + RSP + EMG + BVP
s0132.8 ± 2.9539.8 ± 3.1141.6 ± 4.3446.4 ± 2.4160.3 ± 1.95
s0235.6 ± 1.8233.0 ± 2.2350.0 ± 2.3544.2 ± 3.7053.3 ± 2.36
s0340.8 ± 1.3049.4 ± 2.7055.2 ± 0.8457.4 ± 1.8260.3 ± 2.00
s0440.8 ± 2.1740.2 ± 3.0340.4 ± 4.1645.4 ± 1.9559.4 ± 2.01
s0528.6 ± 1.5241.4 ± 2.4037.2 ± 1.4841.6 ± 2.3052.7 ± 2.83
Table 10. The classification accuracy for each subject using single signal and multimodal signals with ELM.
Table 10. The classification accuracy for each subject using single signal and multimodal signals with ELM.
SubjectPhysiological Sensors
Single Sensor (Acc*/%)Multi Sensors (Acc*/%)
GSRRSPBVPEMGGSR + RSP + EMG + BVP
s0129.2 ± 1.6447.0 ± 2.5546.8 ± 1.9243.4 ± 3.6561.5 ± 2.37
s0230.4 ± 1.6734.4 ± 2.7049.6 ± 3.5042.6 ± 1.8255.4 ± 2.37
s0340.0 ± 3.0854.6 ± 3.3654.4 ± 2.9754.6 ± 2.3062.9 ± 1.29
s0439.6 ± 2.7047.8 ± 2.5642.6 ± 3.2145.8 ± 3.5650.1 ± 2.28
s0534.8 ± 3.3544.4 ± 2.9739.8 ± 3.8343.2 ± 3.1153.9 ± 2.73
Table 11. Comparison of the results using SVM, DT, ELM and proposed classification strategy for each subject.
Table 11. Comparison of the results using SVM, DT, ELM and proposed classification strategy for each subject.
SubjectMethodThe Identification Accuracy of Each Experiment (%)Average (%)
12345678910
s01DT6461586058595960626260.3 ± 1.95
ELM6260645862586460646361.5 ± 2.37
SVM7572727675777273717273.5 ± 2.07
Proposed8080808180827979767879.5 ± 1.65
s02DT5551535350565157525553.3 ± 2.36
ELM5854595253565357575555.4 ± 2.37
SVM6866716363646466636365.1 ± 2.69
Proposed7270767069706870687070.3 ± 2.31
s03DT6258596362615759606260.3 ± 2.00
ELM6563616362626363656262.9 ± 1.29
SVM8482808283818081808281.5 ± 1.35
Proposed8888878986868786878687 ± 1.05
s04DT5862605761585958635859.4 ± 2.01
ELM4851494753504752535150.1 ± 2.28
SVM6366656465606160626162.7 ± 2.21
Proposed7070727273687070686770 ± 1.94
s05DT5052525550515951525552.7 ± 2.83
ELM5758525056535451525653.9 ± 2.73
SVM6967686676697070696869.2 ± 2.70
Proposed7573787382757574757575.5 ± 2.68
Overall averageDT
ELM
SVM
57.2 ± 2.23
56.8 ± 2.21
70.4 ± 2.20
Proposed76.46 ± 1.93
Table 12. Accuracy comparison of various studies in two-dimensional classification.
Table 12. Accuracy comparison of various studies in two-dimensional classification.
MethodAcc* (%)
ArousalValence
Chen et al. [55]69.0967.89
Zhuang et al. [56]71.969.1
Yin et al. [57]77.176.1
Alazrai et al. [58]86.685.8
Kwon et al [59]76.5680.46
Table 13. Accuracy comparison of various studies.
Table 13. Accuracy comparison of various studies.
MenthodEmotionsAcc*(%)
M Zubair and C Yoon [60]HVHA, HVLA,
LVLA, LVHA
49.7
Alazrai et al [58]HVHA, HVLA,
LVLA, LVHA
75.1
Zheng et al [61]HVHA, HVLA,
LVLA, LVHA
69.67
This paperHVHA, HVLA,
LVLA, LVHA
76.46

Share and Cite

MDPI and ACS Style

Pan, L.; Yin, Z.; She, S.; Song, A. Emotional State Recognition from Peripheral Physiological Signals Using Fused Nonlinear Features and Team-Collaboration Identification Strategy. Entropy 2020, 22, 511. https://doi.org/10.3390/e22050511

AMA Style

Pan L, Yin Z, She S, Song A. Emotional State Recognition from Peripheral Physiological Signals Using Fused Nonlinear Features and Team-Collaboration Identification Strategy. Entropy. 2020; 22(5):511. https://doi.org/10.3390/e22050511

Chicago/Turabian Style

Pan, Lizheng, Zeming Yin, Shigang She, and Aiguo Song. 2020. "Emotional State Recognition from Peripheral Physiological Signals Using Fused Nonlinear Features and Team-Collaboration Identification Strategy" Entropy 22, no. 5: 511. https://doi.org/10.3390/e22050511

APA Style

Pan, L., Yin, Z., She, S., & Song, A. (2020). Emotional State Recognition from Peripheral Physiological Signals Using Fused Nonlinear Features and Team-Collaboration Identification Strategy. Entropy, 22(5), 511. https://doi.org/10.3390/e22050511

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