IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS, VOL. 62, NO. 7, JULY 2015 4267
A State-Space EMG Model for the Estimation
of Continuous Joint Movements
Jianda Han, Member, IEEE, Qichuan Ding, Anbin Xiong, and Xingang Zhao, Member, IEEE
Abstract—A state-space electromyography (EMG) model
is developed for continuous motion estimation of human
limb in this paper. While the general Hill-based muscle
model (HMM) estimates only joint torque from EMG signals
in an “open-loop” form, we integrate the forward dynamics
of human joint movement into the HMM, and such an
extended HMM can be used to estimate the joint motion
states directly. EMG features are developed to construct
measurement equations for the extended HMM to form a
state-space model. With the state-space HMM, a normal
closed-loop prediction–correction approach such as the
Kalman-type algorithm can be used to estimate the continu-
ous joint movement from EMG signals, where the measure-
ment equation is used to reject model uncer tainties and
external disturbances. Moreover, we propose a new nor-
malization approach for EMG signals for the purpose of re-
jecting the dependence of the motion estimation on varying
external loads. Comprehensive experiments are conducted
on the human elbow joint, and the improvements of the
proposed methods are verified by the comparison of the
EMG-based estimation and the inertial measurement unit
measurements.
Index Terms—Closed-loop estimation, continuous joint
motion, electromyography (EMG), muscle model.
I. INTRODUCTION
U
SING ELECTROMYOGRAPHY (EMG) as control sig-
nals to realize a “friendly” human–robot interface (HRI)
has been there for assistive robots [1]–[3]. Classification, which
distinguishes different patterns of motion from EMG, is one
of the key techniques for such an HRI [4]. The classification
involves two steps: 1) extracting feature sets from EMG and
2) classifying different motions based on the selected feature
sets [5]. For step 1, the EMG feature sets include EMG am-
plitude, autoregressive coefficients, waveform length, cepstrum
coefficients, and the wavelet packet transform [6]. Certain
feature sets, such as the wavelet packet transform, have to
be used in conjunction with dimension reduction algorithms
such as principal components analysis [7] or linear discriminant
Manuscript received February 15, 2014; revised June 11, 2014 and
October 4, 2014; accepted November 7, 2014. Date of publication
January 1, 2015; date of current version May 15, 2015. This work was
supported in part by the National Natural Science Foundation of China
under Grant 61273355, Grant 61273356, and Grant 61035005.
J. Han and X. Zhao are with the State Key Laboratory of Robotics,
Shenyang Institute of Automation, Chinese Academy of Sciences,
Shenyang 110016, China (e-mail: jdhan@sia.cn).
Q. Ding and A. Xiong are with the State Key Laboratory of Robotics,
Shenyang Institute of Automation, Chinese Academy of Sciences,
Shenyang 110016, China, and also with the University of the Chinese
Academy of Sciences, Beijing 100049, China.
Color versions of one or more of the figures in this paper are available
online at http://ieeexplore.ieee.org.
Digital Object Identifier 10.1109/TIE.2014.2387337
analysis [8] to yield proper signal representation. For step 2,
commonly used classification algorithms are support vector
machine [9], artificial neural networks [10], hidden Markov
model [11], Gaussian mixture model [12], and neuro-fuzzy
algorithm [13].
The aforementioned algorithms are able to identify a limited
number of discrete motion patterns from EMG. However, re-
cently, instead of the discrete patterns, how to determine the
human’s intent of a continuous motion from EMG has become
an active issue. This is due to the requirements of several possi-
ble applications, such as powered prosthesis, exoskeletons, and
rehabilitative robots [14], [15]. In such a system, the human
motion intents have to be continuously recognized from EMG
and further reconstructed as control commands to a robotic
device, so that the robot could match the human’s intent and
then perform efficient assistance.
Involving a physiological muscle model into classification
algorithm is a way to achieve continuous EMG recognition,
and the Hill-based muscle model (HMM) is the most frequently
utilized one [16]–[19]. In [16], an EMG-based forward dynam-
ics model was proposed, which consists of muscle activation
dynamics, Hill-based muscle contraction dynamics, muscu-
loskeletal geometry, and joint forward dynamics. This model
is very complex and involves many unknown physiological
parameters, which limits its applications in real mechatronic
systems. A simplified model for controlling lower extremity
exoskeleton was developed in [18], where a step calibration
process was constructed to optimize the estimation of unknown
parameters.
However, there exist two problems while the normal Hill-
type models are used to estimate continuous motion states of
human limbs. First, a full Hill model always involves many
physiological parameters, which are very difficult to be ac-
curately identified. The complex model will also introduce
extra computational burden. Thus, a simplified model is often
used for real-time applications. However, a simplified model
inherently has modeling errors that will degrade the estimation
accuracy. Normally, the measurement equations in a state-
space model will reduce the effects of modeling errors and
uncertainties. However, the Hill-type models are in “open-loop”
structure; there are no other measurement equations that can be
used as a kind of feedback to improve the estimation accuracy.
Second, the Hill-type models are often used to directly
calculate the joint torque from EMG signals [18]–[21]. For the
applications where the continuous estimations of human joint
are required, for example, to control a robotic device assisting
the movement of human limb [22], we have to further calculate
the motion states, such as angular velocity and position, from
0278-0046 © 2015
IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission.
See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.