EMG Model Calibration and Trajectory Tracking by Remote Rehabilitation Exoskeleton

Preprint

Article

EMG Model Calibration and Trajectory Tracking by Remote Rehabilitation Exoskeleton

This version is not peer-reviewed.

Ubaldo Castro Jimenez,

Edgar Martínez-García^*

Ubaldo Castro Jimenez,

Edgar Martínez-García^*

This version is not peer-reviewed.

Altmetrics

Downloads

151

Views

Comments

Submitted:

18 January 2024

Posted:

19 January 2024

You are already at the latest version

Abstract

This study presents a novel kinematic tracking model, designed for a networked exoskeleton system that is asynchronously taught by a remote therapist. On the server side, the therapist’s rehabilitation exercises are quantitatively assessed using a monocular passive vision system. The resultant analytical metrics are then transmitted asynchronously over the network to patients equipped with exoskeletons. On the client side, the exoskeleton utilizes these analytical metrics as reference paths for exercises, complemented by electromyography (EMG) feedback. This work introduces a calibration approach aimed at estimating angular positions by utilizing EMG observations. The calibration model establishes real-time correlations between polynomial reference positions. This calibration mechanism is integrated into simulations of both upper and lower limb exoskeletons. We further explore redundant kinematics, incorporating an EMG observer for linear, time-variant rehabilitation tracking control. Our methodology is validated using vision-based metric data and experimental EMG measurements for various exercises, including shoulder flexion, elbow flexion, and rowing-like movements. This work also includes computer simulations for tracking the control of rehabilitation exercises, demonstrating the adaptability of the system in reliably, robustly, and effectively following desired trajectories.

Keywords:

Subject:

Computer Science and Mathematics - Robotics

Supplementary Material

supplementary.mp4 (148.72MB )

1. Introduction

From time to time, mankind has had to confront global infectious diseases with deadly impact, as has recently been the case with the Coronavirus, COVID-19. Family members have been compelled to take on in-home responsibilities due to quarantine, including patients following medical treatments. Pandemic-related robotic systematized solutions have been described for multiple scenarios [1]. Similarly, telesystems, communication networks, intelligent machines, and information technologies provide support for easy post-pandemic scenarios [2,3] by keeping us away from contagious scenarios [4,5]. The current technological market has provided a spectrum of modern communication technology such as distributed networked equipment, Internet-based telesystems, telepresence devices, either synchronous or asynchronous modalities [6]. Telesystems applications have evolved into Internet shopping, remote academic/job meetings, virtual scholarly classes, and healthcare support, such as tele-rehabilitation therapies.

An exoskeleton is a wearable robot fitted on a human’s body, designed to enhance biomechanical functions, with the human serving as the controller [12,13]. To be effective, an exoskeleton must provide a sense of lightness and agility, akin to the human bone structure. The advantages of exoskeleton technology extend to enhancing biomechanical capabilities, addressing diverse needs such as the substitution of body parts’ mobility, providing locomotion assistance [7,8], and aiding in rehabilitation [10,11,15,16]. Rehabilitation is particularly crucial in cases of post-corrective surgeries, where constant therapy is required to expedite recovery. Demand for exoskeletons spans various sectors, including athletics, sports, manual labor, and the military, where common surgeries involve conditions like rotator cuff injuries, acromioclavicular articulation injuries, elbow dislocations, hip dislocations, reconstruction of injured anterior cruciate ligaments, knee arthroscopy, and total knee replacement. Prolonged post-surgery rehabilitation plays a critical role in recovery, as prolonged limb immobilization can lead to a significant reduction in muscle strength and mobility. Additionally, various studies have explored control algorithms for lower-limb wearable robot systems [17], as well as for upper limbs [24], to provide power assistance [23]. Wearable robotic hands have been developed for purposes such as haptic feedback, assistive functions, and rehabilitation [34]. These devices also find applications in tasks like grasping [32], the development of multimodal rehabilitation exoskeletons [20,30], and walking assistance [27].

The work by [28] introduced an adaptive neural control system for rehabilitating lower limbs with unknown model dynamics, focusing on adapting to the complex dynamics of human-exoskeleton interactions. In this study, we propose a feedback kinematic tracking model for a networked exoskeleton, which is taught asynchronously by a remote therapist (Figure 1abc). The therapist’s rehabilitation exercises are visually assessed and transformed into analytical polynomial Cartesian functions through a least-square fitting approach. A noteworthy innovation presented in this work is the transmission of a concise stream of polynomial coefficients, as opposed to the entire set of Cartesian data. On the patient’s side, these Cartesian polynomials are employed to recreate the therapist’s experimental data with higher resolution. Additionally, by introducing a constant time period, denoted as T, we achieve the reconstruction of Cartesian components, x and y, as periodic functions of time using Fourier series. The integration of Cartesian and electromyography (EMG) data as feedback, similar to the approach in Kang et al. [9], represents another significant contribution of this work. Notably, our model is designed to infer the posture of Cartesian limbs from EMG signals, allowing for recalibration during new exercises. This feature offers a generalized approach, enabling the control law to adapt to the structure of different exoskeletons simply by readjusting their Jacobian matrices. Furthermore, a linear trajectory smoother model is introduced in cases where patients have biomechanical movement constraints. These smoothing functions facilitate a gradual convergence towards the desired reference trajectory shape. To evaluate the efficacy of the proposed exoskeleton, we conducted extensive simulations involving hundreds of tracking control rehabilitation exercises.

This study introduces a distinctive approach to remote therapy within the realm of tele-rehabilitation. Specifically, it proposes a feedback kinematic tracking model for a networked exoskeleton that undergoes asynchronous visual instruction from a remote therapist. The emphasis lies in transmitting concise streams of polynomial coefficients instead of the entire set of Cartesian data, subsequently transformed into analytical polynomial Cartesian functions. Additionally, the study underscores the integration of Cartesian and electromyography (EMG) data as feedback, facilitating the inference of limb posture from EMG signals and recalibration during new exercises. In the existing literature various approaches within closely related fields of study, including works by [36,37,38,39,40,41], partially intersect with the themes presented in this work. However, the authors remain uncertain as to whether fully analogous studies have been previously published.

The proposed solution specifically focuses on tele-rehabilitation, where a therapist’s exercises are visually assessed and transmitted to the patient’s exoskeleton. The key innovation of the proposed solution lies in the transmission of concise streams of polynomial coefficients, rather than transmitting the entire set of Cartesian data. This approach provides flexibility in tracking and control, making the rehabilitation process highly accurate and efficient. Overall, the proposed solution offers a generalized approach to tele-rehabilitation, making it more accessible and accommodating to a wide range of patient needs. It combines clinical precision with technological sophistication, providing an effective and efficient means of rehabilitation.

The proposed feedback kinematic tracking model for a networked exoskeleton taught asynchronously by a remote therapist is required to ensure accurate, personalized, and efficient rehabilitation for patients, regardless of their location or access to specialized facilities. In addition, it provides i) accuracy and precision: The feedback kinematic tracking model follows the desired trajectory and provides the necessary level of accuracy which is crucial for effective rehabilitation. ii) Adaptability and customization: The feedback model allows for adjustments and adaptations based on the patient’s specific biomechanical needs, and the use of electromyography (EMG) data as feedback allows for the inference of Cartesian posture from EMG signals, enabling recalibration during new exercises. iii) Remote guidance and supervision: Asynchronous teaching by a remote therapist regardless of their location, enables benefits for patients who may not have access to specialized rehabilitation facilities or therapists in their vicinity. iv) Flexibility and convenience: Asynchronous teaching allows patients to perform and follow the prescribed exercises and receive feedback from the therapist at a time that suits them best. v) Efficiency and cost-effectiveness: By utilizing a networked exoskeleton and asynchronous teaching, the rehabilitation process can save time and resources by avoiding frequent visits to rehabilitation centers, while still receiving expert guidance and support.

This paper is structured as follows: Section 2 delves into the therapist-side tele-rehabilitation, elucidating the teaching process facilitated by a vision system. In Section 3, we unveil the proposed EMG calibration and observation measurement models. Following this, Section 4 delves into the exoskeleton’s motion constraints and its kinematic control law. A comprehensive summary of the research’s results is provided in Section 5. Finally, Section 6 offers a set of conclusions drawn from this study.

2. Human therapist remote visual teaching

The proposed networked distributed system is structured as an asynchronous tele-system architecture, with the therapist’s side and the patients constituting the client side. The server plays a pivotal role by providing geometric data acquired visually from the therapist’s exercise sessions in front of a camera, transmitted over the network. An artificial vision system captures the data related to the training trajectories taught by the therapist (Figure 1d, left side).

Notably, the interaction session between the therapist and the patient occurs asynchronously, eliminating the need for simultaneous network connections. This approach stands in contrast to traditional tele-operation methods [19]. The exoskeleton has the flexibility to connect at any given date and time, allowing the download of therapy exercises in the form of analytic trajectories, represented by coefficients (Figure 1d, right side). The server’s computer is equipped with a monocular vision system that accurately measures the angles and positions of colored landmarks strategically placed at the therapist’s limb joints (Figure 3ab-left). Remarkably, the therapist does not require any specialized equipment during the exercise sessions.

In light of the evident benefits of rehabilitation exoskeletons ([22]), we selected commonly practiced exercises from a recognized exercise program ([18]) to serve as the basis for our study. These exercises, whether post-injury or post-surgery, were adopted as trajectories (Figure 1abc). For example, exercises such as the rotator cuff, shoulder and elbow flexion exercises were employed to teach various motions (Figure 1ac). Similarly, we incorporated exercises like knee extensions (Figure 1b) and the stand-up paddle exercise with an elastic band (Figure 1c) to broaden the scope of our research.

Figure 1. Tele-training rehabilitation exercises. a)-b) Elbow, knee flexion-extension. c) Shoulder sagittal rotation. d) Telesystem architecture.

The proposed tele-system was crafted on standard-capability computers running the Linux Gentoo operating system, version 4.4. The entire codebase was developed using the GNU C/C++ compiler, making extensive use of multiple open-source libraries, including openCV, Armadillo C++, the Open Dynamics Engine (ODE), and Gazebo. Communication between components was facilitated using the TCP/IP protocol via sockets. The vision sensor employed was a cost-effective USB camera (see Table 1). For the acquisition of EMG biosignals, we utilized the gMOBIlab+ device with two EMG channels. To bring our vision to life, we implemented a virtual exoskeleton prototype with original design mechanisms. This virtual exoskeleton operates under the control of human EMG signals, functioning as a master-slave system.

This work places some emphasis on vision-based teaching. To streamline the measurement process, we strategically position colored landmarks at the joints of the therapist’s limbs, as depicted in Figure 3. These landmarks serve as the key reference points for measuring angles and positions during exercise sessions. Our approach to measuring the metric trajectories of exercises involves a two-step process. First, we segment the joints’ landmarks based on their RGB color. Subsequently, we employ a nominal resolution geometric model to track the centroid positions of these metric landmarks. In this approach, we establish the reference coordinate system at the location of the visual sensor. In line with the extrinsic parameters outlined in Figure 2, we define an RGB image as a cubic matrix

I \in R^{m \times n \times 3}

, comprising the union of the three primary color channels: red

I_{R}

, green

I_{G}

, and blue

I_{B}

. To extract the shoulder’s landmark, we apply a threshold

ζ_{S}

to the red channel

I_{R}

and incorporate a zeros matrix

0 \in R^{m \times n}

. This process yields a binary segmented image denoted as

I_{S}

and is achieved through the following Boolean equation:

I_{S} = (I_{R} > ζ_{S}) \cup (I_{G} \cap 0) \cup (I_{B} \cap 0) .

(1)

Similarly, the isolation of the elbow’s joint landmark is achieved through the utilization of a binary image denoted as

I_{E}

. The effectiveness of this process hinges on the presence of a critical threshold denoted as

ζ_{E}

, which assumes a pivotal role in the Boolean criterion:

I_{E} = (I_{R} \cap 0) \cup (I_{G} < ζ_{E}) \cup (I_{B} \cap 0) .

(2)

In a similar manner, the identification of the wrist landmark is carried out through the utilization of a binary image represented as

I_{W}

The determination of this landmark relies on the precise specification of a threshold value denoted as

ζ_{W}

which is an integral part of the following criterion:

I_{W} = ((I_{R} \cap 0) \cup (I_{G} \cap 0)) \cup (I_{B} > ζ_{W}) .

(3)

For each segmented region in the images

I_{R}

I_{G}

, and

I_{B}

, the position of the landmarks in spatial coordinates is determined by their centroids, denoted as column

p_{c}

and row

p_{r}

, as expressed by:

p_{c} = \frac{1}{2} \sum_{i} c_{i}, p_{r} = \frac{1}{2} \sum_{i} r_{i},

(4)

In this context, the spatial coordinates

c_{i}

and

r_{i}

, where both

c_{i}

and

r_{i}

belong to a segmented region

R

, represent the pixels within the region. We will now proceed to offer a comprehensive description of the geometric visual model, as depicted in Figure 2.

Figure 2. Monocular vision-based geometrical model on the therapist’s side.

Geometrically, the equality for the vertical field of view ((5)) is valid.

tan (\frac{φ_{r}}{2}) = (\frac{h / 2}{z}),

(5)

Here, h represents the height of the object of interest, and z denotes the metric depth of the human therapist from the camera. Thus, building upon the previous expression, we can determine

φ_{r}

, which represents the vertical camera angle quantized in terms of image rows.

φ_{r} = 2 {tan}^{- 1} (\frac{h}{2 z}) .

(6)

With respect to (w.r.t.) the horizontal field of view,

tan (\frac{φ_{c}}{2}) = (\frac{w / 2}{z}),

(7)

then, dropping off the angle along image’s columns

φ_{c}

φ_{c} = 2 {tan}^{- 1} (\frac{w}{2 z}) .

(8)

We estimate the vertical metric resolution factor, denoted as

f_{h}

, as well as the horizontal metric resolution factor, denoted as

f_{w}

f_{h} = 2 tan (\frac{φ_{r}}{2} \cdot z), f_{w} = 2 tan (\frac{φ_{c}}{2} \cdot z) .

(9)

Subsequently, we determine the vertical metric position

Δ_{r}

and the horizontal metric position

Δ_{c}

of the center pixel

(p_{r}, p_{c})

of a landmark.

Δ_{r} = \frac{f_{h}}{p_{r}}, Δ_{c} = \frac{f_{w}}{p_{c}} .

(10)

Next, in the therapist’s side, the superior limb kinematic motion model in 3D space is provided by the expression

(\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} ℓ_{1} sin (ϕ_{0}) cos (ϕ_{1}) + ℓ_{2} sin (ϕ_{0}) cos (ϕ_{1} + ϕ 2) \\ ℓ_{1} sin (ϕ_{1}) + ℓ_{2} sin (ϕ_{1} + ϕ 2) \\ ℓ_{1} cos (ϕ_{0}) cos (ϕ_{1}) + ℓ_{2} cos (ϕ_{0}) cos (ϕ_{1} + ϕ 2) \end{matrix}),

where depth is affected by the yaw angle

ϕ_{0}

and double pitch joits

ϕ_{1}, ϕ_{2}

, and

ℓ_{1}

and

ℓ_{2}

to denote the actual metric lengths of the upper arm and forearm, respectively. When representing a 2D exercise scenario, the yaw angle

ϕ_{0} = π / 2

. Therefore, the depth dimension becomes irrelevant since the arm’s length is not projected, and the trajectory is drawn on a plane in front of the camera. The links

ℓ_{1}, ℓ_{2}

visual measurements,

ℓ_{1} = \sqrt[2]{{(x_{S} - x_{E})}^{2} + {(y_{S} - y_{E})}^{2}}

(11)

and

ℓ_{2} = \sqrt[2]{{(x_{E} - x_{W})}^{2} + {(y_{E} - y_{W})}^{2}} .

(12)

Similarly, we use the variables

θ_{1}

and

θ_{2}

to denote the joint angles of the limb.

θ_{1} = {tan}^{- 1} (\frac{y_{E} - y_{S}}{x_{E} - x_{S}}), θ_{2} = {tan}^{- 1} (\frac{y_{W} - y_{E}}{x_{W} - x_{E}})

(13)

We illustrate a shoulder exercise in Figure 3, and in Figure 9, we display the resulting Cartesian trajectory, with red points indicating the Cartesian positions of the blue wrist landmark.

Figure 3. Visual instruction of various exercises by the therapist.

In summary, the rehabilitation trajectory, denoted as

S = Θ_{0}, Θ_{1}, Θ_{2}, \dots, Θ_{N}

, is comprised of a set of joint angles

Θ_{i} = {(θ_{1}, θ_{2})}^{⊤}

and limb parameters

ℓ_{1}

and

ℓ_{2}

In Figure 4, a series of repetitions stemming from exercises performed by the therapist is presented. These exercises encompass a single experimental repetition involving shoulder flexion (depicted in Figure 4a), along with three experimental repetitions focusing on elbow flexion (illustrated in Figure 4b). Additionally, three experimental repetitions of a rowing-like motion are depicted (Figure 4c), wherein both the shoulder and elbow joints operate in unison during the exercise. It is imperative to emphasize that the reference coordinates, denoted as

(0, 0)

, are firmly anchored at the centroid of the shoulder.

Additionally, in Figure 5, we present angular motion measurements for all joints utilizing visual odometry, with a focus on the

x - y

plane. These angular measurements correspond to the Cartesian motion data depicted in Figure 4. In Figure 5a (left), we observe an experiment exclusively involving shoulder flexion, where a linear rigid arm is used, and there are a total of five repetitions. It’s important to note that this exercise is the simplest, resulting in data points with a higher degree of overlap. On the right side of Figure 5a, while still involving shoulder flexion, we observe slight angular perturbations in the wrist, averaging approximately

9^{\circ}

. These perturbations were primarily caused by minor limb vibrations from the therapist and were further exacerbated by contrast and brightness noise in the measurements. Figure 5b displays angular measurements for elbow flexion exercises. To ensure clarity without excessive data overlap, we have chosen to present three distinct experiments, each with multiple repetitions. On the left side of Figure 5b, we observe angular perturbations in the wrist during elbow flexion exercises. Conversely, the right side of Figure 5b depicts three separate experiments involving elbow-flexion exercises, displaying multiple repetitions that result in periodic oscillations. This presentation highlights variations in speed during the experiments, reflecting biomechanical consistency and differences in execution.

Figure 5c, in its entirety, represents three distinct experiments of rowing-like motions, each of which consists of six repetitions. This comprehensive view enables us to examine the angular dynamics in these exercises in greater detail. On the left side of Figure 5c, we can observe the angular magnitudes of elbow flexions. Meanwhile, the right side of the same figure portrays shoulder flexions, which synchronize their motion with the elbow flexions. On the therapist’s side, the process involves employing a least square polynomial fitting approach to derive a Cartesian function

y (x)

. This method is used to process the vision-based 3D measured data, resulting in the determination of polynomial coefficients of the n-th degree, as described in the following equation:

(\begin{matrix} a_{0} \\ a_{1} \\ ⋮ \\ a_{n} \end{matrix}) = {(\begin{matrix} n & \sum_{i} x_{i} & \dots & \sum_{i} x_{i}^{n} \\ \sum_{i} x_{i} & \sum_{i} x_{i}^{2} & \dots & \sum_{i} x_{i}^{n + 1} \\ ⋮ \\ \sum_{i} x_{i}^{j} & \sum_{i} x_{i}^{j + 1} & \dots & \sum_{i} x_{i}^{2 n} \end{matrix})}^{- 1} \cdot (\begin{matrix} \sum y_{i} \\ \sum x_{i} y_{i} \\ ⋮ \\ \sum_{i} x_{i}^{n} y_{i} \end{matrix}),

(14)

In this manner, we obtain a parametric function that accurately fits the exercise’s Cartesian data, following a polynomial form:

y (x) = a_{0} + a_{1} x + a_{2} x^{2} + a_{3} x^{3} + \dots + a_{n} x^{n} .

(15)

This approach significantly reduces the trade-off during data streaming to only a few coefficients and the time period T of one exercise oscillation. These coefficients require a minimal amount of bytes for transmission over the network, thereby mitigating significant time delays. Figure 6 illustrates experimental data from teaching trajectories that have been effectively fitted by various polynomials, exemplifying the practicality and efficiency of this approach.

By representing the exercise trajectory as a continuous polynomial function, denoted as

y (x)

, available on the exoskeleton’s side, we open the door to the replication of experimental data at higher metrical resolutions. With the time constant of a period T at our disposal, the exoskeleton’s side can autonomously convert the high-resolution Cartesian data into new components as functions of time. Using Fourier series, we generate new Cartesian components,

x (t)

and

y (t)

, as functions of time. Figure 7 effectively demonstrates this process, showcasing a Cartesian reference trajectory as

x (t)

versus

y (t)

, realistically reproducing three experiments involving elbow flexion. This approach, as presented in this work, offers sophisticated analytical functions that serve as effective substitutes for the therapist’s exercises, which typically involve large amounts of experimental data. In doing so, it efficiently addresses the challenges of long transmission time delays and the handling of extensive data tables, ultimately yielding parametric functions that streamline the entire process.

Moreover, in scenarios where biomechanical constraints affecting a patient’s movements are detected during exoskeleton deployment, we employ smoothing functions to facilitate gradual transitions that accommodate these constraints. One such smoothing function, based on a linear approach, is detailed below:

y_{r e f} = κ_{τ} + κ_{s} y (x),

(16)

In this equation, the original therapist’s trajectory is regarded as the ideal Cartesian reference motion, denoted as

y (x)

. Therefore, we introduce two key parameters:

κ_{τ}

, representing the translation factor of a family of functions, and

κ_{s}

, representing the scaling function’s coefficient. When

κ_{τ} = 0

and

κ_{s} = 1

, the exercise function aligns precisely with the therapist’s reference trajectory. However, by adjusting the values of these smoothing coefficients, we achieve gradual convergence toward the therapist’s ideal trajectory. Figure 8 purposely illustrates various smoothing teaching trajectories that progressively approach the original therapist’s reference trajectory, represented by the continuous line, underscoring the flexibility and efficacy of this approach.

Indeed, it’s important to note that both

κ_{τ}

and

κ_{s}

are coefficients specific to each exercise trajectory, differences of patient’s limbs size or patient’s motion conditions. The coefficients are manually adjusted individually for different exercises based on the patient’s comfort. This adaptability ensures that the smoothing coefficients are tailored to the characteristics of each exercise, optimizing their effectiveness.

3. Rehabilitation trajectory model

From the perspective of the patient-exoskeleton interface, the instructional trajectory

S

is inherently expressed in polar form. To facilitate processing, it is subsequently transformed into Cartesian points at the wrist location, represented as

x_{i}^{r e f} = ℓ_{1} cos (θ_{1}) + ℓ_{2} cos (θ_{1} + θ_{2})

and

y_{i}^{r e f} = ℓ_{1} sin (θ_{1}) + ℓ_{2} sin (θ_{1} + θ_{2})

, commonly referred to as the ideal or reference trajectory

p^{r e f}

. In lieu of managing an extensive dataset comprising hundreds or thousands of points, this study derives an adaptive analytical reference model expressed as functions of time, akin to the approach proposed by [25]. Given the periodic nature of a rehabilitation exercise assumed to be a Cartesian function, the suggested methodology involves obtaining xx and yy separately with respect to time. The utilization of a general Fourier series is deemed apt for representing a suitable function of time, providing a concise and analytically manageable characterization of the patient’s movement trajectory.

x (t) = α_{0}^{x} + α_{k}^{x} cos (k ω t) + β_{k}^{x} sin (k ω t) + \dots

(17a)

and

y (t) = α_{0} + α_{k}^{y} cos (k ω t) + β_{k}^{y} sin (k ω t) + \dots

(17b)

In the general formulation (18)-(20) the Fourier series estimate the coefficients for (17). The parameter n represents the number of data to be fitted through the Fourier series. However, the

k_{t h}

number of Fourier coefficients are predetermined with a sufficiently large value to appropriately accommodate the trajectory data. Equations 18, 19, and 20 comprise the Fourier series tailored to fit the trajectory data provided by the therapist, expressing it as a trigonometric function.

A C++ program was implemented to autonomously compute the discrete Fourier coefficients. Regarding the independent terms, the average of the discrete data is presented as follows:

α_{0}^{x} = \frac{1}{n} \sum_{k = 1}^{n} x_{k}, and α_{0}^{y} = \frac{1}{n} \sum_{k = 1}^{n} y_{k} .

(18)

Similarly, even though initially presented in continuous-time, the discrete Fourier series was formulated to derive the remaining coefficients for the Cartesian component x vs t by

α_{k}^{x} = \frac{2}{n} \int_{t} x (t) cos (k ω t) d t, β_{k}^{x} = \frac{2}{n} \int_{t} x (t) sin (k ω t) d t,

(19)

likewise, for the Cartesian component y vs t,

α_{k}^{y} = \frac{2}{n} \int_{t} y (t) cos (k ω t) d t, β_{k}^{y} = \frac{2}{n} \int_{t} y (t) sin (k ω t) d t .

(20)

Therefore, we have now acquired the Cartesian components of an

x y

trajectory as reference functions of time, denoted as

x (t)

and

y (t)

. The numerical coefficients related to time obtained to align with the experimental reference data of shoulder flexion-extension, illustrated in Figure 9, are:

\begin{matrix} x_{t}^{r e f} = 0.17 - 0.05 cos (ω t) - 0.06 cos (2 ω t) + 0.01 (sin (2 ω t) - sin (3 ω t)) \end{matrix}

(21)

and

\begin{matrix} y_{t}^{r e f} = - 0.30 - 0.22 cos (ω t) + 0.05 sin (ω t) - 0.02 sin (2 ω t) + 0.01 (cos (3 ω t) - sin (3 ω t)) . \end{matrix}

(22)

The experimental coefficients obtained are presented for illustrative purposes, specifically related to Figure 9. It’s essential to note that their numeric values may vary for different exercises or when administered by a human therapist, accounting for differences in limb size.

Figure 9. Periodic positions of shoulder flexion-extension fitted by Fourier series.

3.1. EMG sensing model calibration

This study suggests deriving a sensing model for estimating the Cartesian position of the limb through EMG measurements. Given the challenge of achieving repeatability in EMG data from the same exercises and subjects, an EMG calibration model is established at the onset of training. Real-time motion prediction becomes feasible using raw EMG data, typically encompassing combined muscle motions [33]. Table 2 provides a comprehensive overview, detailing the type of therapy and the associated muscles involved in obtaining EMG measurements.

Based on the electrodes illustrated in Figure 10c, the experimental raw EMG data is presented in Figure 14. This data pertains to the shoulder-elbow-wrist region and corresponds to multiple periodic flexion-extension exercises.

The series of plots illustrated in Figure 11 portrays three distinct experiments involving the same exercise, specifically elbow flexions, each comprising approximately five repetitions. In this configuration, Channel Ch5 was situated at the bicep brachii muscle, while concurrently, Channel Ch6 was positioned at the tricep brachial muscle. Notably, the magnitudes of amplitude (

μ V

) exhibit coherence and consistency across both channels.

The subsequent Figure 12 illustrates an experiment involving shoulder flexion. In this instance, Channel Ch5 was positioned at the Pectoralis Major muscle, and Channel Ch6 was situated at the Deltoid Anterior muscle. The plots represent approximately five oscillations (flexions) with amplitude magnitudes measured in

μ V

. It is noteworthy that, in this case, both frequency and amplitudes consistently differ from those observed in other exercises.

Figure 13 presents experimental data related to rowing-like movements, showcasing three distinct experiments involving the same exercise with approximately six rowing flexions each. Channel Ch5 was positioned at the Trapezius middle muscle, while Channel Ch6 was situated at the Deltoid anterior muscle. The amplitude magnitudes (

μ V

) and oscillation frequency exhibit distinct behavioral patterns, maintaining substantial consistency across both channels.

Figure 13. Three different experiments of rowing exercise using EMG channels Ch5 (Trapezius middle) and Ch6 (Deltoid anterior).

Figure 14. Native EMG signals

\tilde{E}

during simultaneous shoulder-elbow flexion extension and wrist.

Figure 14. Native EMG signals

\tilde{E}

during simultaneous shoulder-elbow flexion extension and wrist.

The activity level of muscles is inherently noisy due to the interaction of bi-joint muscles with movements in other joints. The EMG calibration measurement is acquired by initially filtering the noisy native EMG signal

\tilde{E}

(Figure 14) into a linear envelope

{\bar{E}}_{t}

within segments of

N_{k}

samples. Subsequently, a 2^nd order low-pass filter with a cutoff frequency

f_{c} = 20 H z

is applied,

{\bar{E}}_{t} = |\tilde{E} - \frac{1}{N_{k}} \sum_{1}^{N_{k}} {\tilde{E}}_{i}|, E_{t}^{c a l} = \frac{{\bar{E}}_{t}}{\sqrt[2]{1 + {(\frac{f_{t}}{f_{c}})}^{4}}} .

(23)

Consequently, an instantaneous EMG calibration measurement value, denoted as

E_{t}^{c a l}

, is deduced from the periodicity of the filtered EMG signal

{\bar{E}}_{t}

. To obtain an analytical representation of

E^{c a l}

over a period of a specified exercise, a fitness model employing a Fourier series in terms of time, denoted as

E_{t}^{c a l}

, is derived.

\begin{matrix} E_{t}^{c a l} = 56.60 - 29.65 cos (ω t) - 15.53 sin (ω t) \\ - 2.18 cos (2 ω t) + 7.73 sin (2 ω t) - 1.37 cos (3 ω t) \\ + 3.26 sin (3 ω t) + 1.76 cos (4 ω t) - 2.15 sin (4 ω t) \\ + 0.13 cos (5 ω t) + 2.58 sin (5 ω t) . \end{matrix}

(24)

The numerical coefficients in (24) are derived from the experimental shoulder-elbow flexion-extension exercise shown in Figure 14. Let’s introduce a polynomial function

\hat{E} (t)

as a time-dependent function capable of estimating EMG measurements over time. This function can then be compared against the calibrated model

E_{t}^{c a l}

, which serves as a reference function over time. To approximate

E_{t}^{c a l}

, we employ the Taylor series (Maclaurin), representing it as

\hat{E} (t)

\hat{E} (t) \approx \frac{E_{0}^{(0)}}{0!} t^{0} + \frac{E_{0}^{(1)}}{1!} t + \frac{E_{0}^{(2)}}{2} t^{2} + \dots \frac{E_{0}^{(m)}}{m!} t^{m} .

(25)

A polynomial approximation is advantageous for an observer due to its simplicity in terms of differentiation and integration with respect to time. In Figure 15a, a gradual polynomial fitting up to the 80^th degree aligns with our calibration model (24). Additionally, Figure 15b visually represents the correlation between raw EMG data, overlaid with calibrated signals (filtered out), and their corresponding Cartesian positions over time during periodic exercises.

3.2. Observability based on real-time EMG

Given that the polynomial observation model

\hat{E} (t)

is expressed as a function of time t, its inverse model

t ({\hat{E}}_{t})

provides a time value consistent with representing instantaneous Cartesian positions

x (t)

and

y (t)

. In the case of the experimental EMG data from Figure 16, an 80^th grade polynomial was obtained. As obtaining a non-linear analytical solution is not trivial, its inverse result for time is numerically calculated through the Newton-Raphson method for nonlinear equations. The instantaneous measurement

{\hat{E}}_{t}

substitutes the left side of the following expression,

\hat{E} (t) = a_{0} + a_{1} \hat{E} + a_{2} {\hat{E}}^{2} + a_{3} {\hat{E}}^{3} + \dots + a_{i} {\hat{E}}^{i}

(26)

where the time

t_{k}

is the unknown value. Therefore, the Newton-Raphson method necessitates the first-order derivative with respect to time

d \hat{E} / d t

\hat{\dot{E}} (t) = \frac{d \hat{E} (t)}{d t} = a_{1} + a_{2} {\hat{E}}_{t} + 2 a_{3} {\hat{E}}_{t}^{2} + \dots + a_{i} (i - 1) {\hat{E}}_{t}^{i - 1}

(27)

hence, a solution for

t_{k}

is provided in the following manner,

{\hat{t}}_{k} = {\hat{t}}_{k - 1} - \frac{\hat{E} (t)}{\hat{\dot{E}} (t)},

(28)

recursively, the best solution approximation is iteratively refined until the relative numeric error falls below

ϵ_{t}

|\frac{{\hat{t}}_{k} - {\hat{t}}_{k - 1}}{{\hat{t}}_{k}}| < ϵ_{t}

(29)

and now Cartesian models are satisfied as a functions of time

t_{k}

as well

\hat{x} (\hat{t}) = b_{0} + b_{1} \hat{t} + b_{2} {\hat{t}}^{2} + \dots + b_{i} {\hat{t}}^{i},

(30)

\hat{y} (\hat{t}) = c_{0} + c_{1} \hat{t} + c_{2} {\hat{t}}^{2} + \dots + c_{i} {\hat{t}}^{i} .

(31)

The solution for

t_{k}

establishes synchronicity between

{\hat{E}}_{t}

and the expected Cartesian positions (

x_{t}^{r e f}

y_{t}^{r e f}

) in (21) and (22), respectively. Figure 16a-b display the native data and filtered signals for two EMG channels. Similarly, Figure 16c presents only the filtered signals.

Figure 16. Shoulder flexion-extension EMG measurement.

4. Exoskeletons morphology and control law

In this study, a proposed exoskeleton model incorporates a well-designed mechanism to facilitate standard rehabilitation exercises (Figure 18). Alternative methodologies, such as adaptive geometrical design, have been documented in previous literature [35]. The design of speed reducers (e.g., the shoulder mechanism shown in Figure 17a) primarily incorporates pulse encoders for joint rotation sensing. These mechanisms are intended to enforce limb tracking along the reference trajectories.

The design assumes all encoders have a constant angular resolution R [pulses/rev]. The number of pulses

η_{t}

measures the motor’s instantaneous angle

θ_{0}

[rad], and its angular speed

w_{0}

[rad/s].

θ_{0} = \frac{2 π}{R} η_{t}, w_{0} = \frac{2 π}{R} (\frac{η_{t_{2}} - η_{t_{1}}}{t_{2} - t_{1}}) .

(32)

All gears’ rotary speeds are described by (33) and illustrated in Figure 17b. Gears with a radius

r_{i}

yield speeds

w_{i}

as functions of the motor’s velocity

w_{0}

(\begin{matrix} w_{1} \\ w_{23} \\ w_{5} \\ w_{7} \\ w_{8} \end{matrix}) = w_{0} (\begin{matrix} \frac{r_{0}}{r_{1}} \\ \frac{r_{0}}{r_{1}} \\ \frac{r_{0} r_{23}}{r_{1} r_{5}} \\ \frac{r_{0} r_{23}}{r_{1} r_{5}} \\ \frac{r_{0} r_{23} r_{7}}{r_{1} r_{5} r_{8}} \end{matrix}) .

(33)

By expressing (33) in vector notation, the time-invariant linear system (34) emerges,

w = w_{0} k .

(34)

This model encapsulates the functional representation of the joints, providing a framework to deduce Cartesian positions over time for feedback purposes. In Figure 18, the kinematic parameters of the upper and lower limbs are illustrated, featuring a structured organization of joint axes.

Figure 18. Exoskeleton mechanism structure.

Therefore, it follows a dynamic controller at the level of the exoskeleton’s joint gearing, rather than at the level of the limbs’ global kinematic-EMG calibration. Employing an Euler-Lagrange approach, the general Lagrangian model

L_{i}

for the

i^{t h}

gear, transmitting rotary kinetic and potential energy, is:

L_{i} = r_{i}^{2} m_{i} {\dot{ϕ}}_{i} - m_{i} g r_{i} sin (ϕ_{i}),

(35)

and developing the following differential equation

τ_{i} = \frac{d}{d t} (\frac{d L}{d {\dot{ϕ}}_{i}}) - \frac{d L}{d ϕ_{i}} .

(36)

Subsequently, by considering and substituting joint kinematics from (34), and simplifying in terms of controlled motor motion

ϕ_{0}

, the following dynamic vector

τ = {(τ_{1}, τ_{23}, τ_{5}, τ_{7}, τ_{8})}^{⊤}

is deduced,

τ = (\begin{matrix} r_{0} r_{1} m_{1} {\ddot{ϕ}}_{0} + m_{1} g r_{1} cos (\frac{r_{0}}{r_{1}} ϕ_{0}) \\ \frac{r_{23}^{2} r_{0}}{r_{1}} m_{23} {\ddot{ϕ}}_{0} + m_{23} g r_{23} cos (\frac{r_{0}}{r_{1}} ϕ_{0}) \\ \frac{r_{0} r_{23} r_{5}}{r_{1}} m_{5} {\ddot{ϕ}}_{0} + m_{5} g r_{5} cos (\frac{r_{0} r_{23}}{r_{1} r_{5}} ϕ_{0}) \\ \frac{r_{0} r_{23} r_{7}^{2}}{r_{1} r_{5}} m_{7} {\ddot{ϕ}}_{0} + m_{7} g r_{7} cos (\frac{r_{0} r_{23}}{r_{1} r_{5}} ϕ_{0}) \\ \frac{r_{0} r_{23} r_{7} r_{8}}{r_{1} r_{5}} m_{8} {\ddot{ϕ}}_{0} + m_{8} g r_{8} cos (\frac{r_{0} r_{23} r_{7}}{r_{1} r_{5}} ϕ_{0}) \end{matrix}) .

(37)

Expressing

τ

in terms of the explicit vectors for inertia

m

and gravity force

g

τ = (\begin{matrix} r_{0} r_{1} m_{1} {\ddot{ϕ}}_{0} \\ \frac{r_{23}^{2} r_{0}}{r_{1}} m_{23} \\ \frac{r_{0} r_{23} r_{5}}{r_{1}} m_{5} \\ \frac{r_{0} r_{23} r_{7}^{2}}{r_{1} r_{5}} m_{7} \\ \frac{r_{0} r_{23} r_{7} r_{8}}{r_{1} r_{5}} m_{8} \end{matrix}) {\ddot{ϕ}}_{0} + (\begin{matrix} m_{1} g r_{1} cos (\frac{r_{0}}{r_{1}} ϕ_{0}) \\ m_{23} g r_{23} cos (\frac{r_{0}}{r_{1}} ϕ_{0}) \\ m_{5} g r_{5} cos (\frac{r_{0} r_{23}}{r_{1} r_{5}} ϕ_{0}) \\ m_{7} g r_{7} cos (\frac{r_{0} r_{23}}{r_{1} r_{5}} ϕ_{0}) \\ m_{8} g r_{8} cos (\frac{r_{0} r_{23} r_{7}}{r_{1} r_{5}} ϕ_{0}) \end{matrix})

(38)

and expressing this relationship in vector notation, the following dynamic equation is formulated,

τ = m {\ddot{ϕ}}_{0} + g (ϕ_{0}) .

(39)

It follows that a controller is deduced from the dynamic equation (39). Thus, rewriting the torque vector with a differential equation notation,

τ = m \frac{d \dot{ϕ}}{d t} + g (ϕ_{0})

(40)

and by separating different differentials,

(τ - g (ϕ_{0})) d t = m d {\dot{ϕ}}_{0}

(41)

solving by definite integration,

(τ - g (ϕ_{0})) \int_{t_{2}}^{t_{1}} d t = m \int_{ϕ_{1}}^{ϕ_{2}} d {\dot{ϕ}}_{0},

(42)

then by algebraically arranging and changing subindex notation, a model-based controller is yielded

{\dot{ϕ}}_{0_{t + 1}} = {\dot{ϕ}}_{0_{t}} + m^{+} (τ^{r e f} - g (ϕ_{0})) (t_{2} - t_{1}),

(43)

where

m^{+}

is the inertia vector pseudoinverse2. Thus, the speeds of the exoskeleton’s joint mechanisms are recursively controlled according to a specified reference torque

τ^{r e f}

4.1. Limbs kinematic model

Based on the mechanical model presented in Figure 18, a set of Cartesian equations was deduced for kinematic analysis, forming the foundation for position control. The non-linear and non-square systems obtained in (44) and (45) evidently lack inverse analytical solutions. Therefore, an algebraic numeric recursive approach was developed. Expressions (44a)-(44c) describes the upper limb kinematics of the 5 rotary joints. Notice that

ϕ_{0}

and

ϕ_{3}

are rolling axis. Thus, for x,

x = l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{3} s_{123},

(44a)

for coordiante y

y = - s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{124}) + l_{3} s_{3} s_{124} s_{0}

(44b)

and for coordinate z,

z = - s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{124}) - l_{3} s_{3} s_{124} s_{0} .

(44c)

Expressions (45a)-(45c) are the kinematic equations of the lower limb with 4 rotary joints,

x = l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{123}

(45a)

and

y = - s_{0} (l_{1} s_{1} + l_{2} s_{12} - l_{3} s_{123})

(45b)

and

z = - s_{0} (l_{1} s_{1} + l_{2} s_{12} - l_{3} s_{123}) .

(45c)

The instantaneous Jacobian matrix represents the rate of change of coordinates

x y z

w.r.t. the joint-space vector

Φ = {(ϕ_{0}, ϕ_{1}, ϕ_{2}, \dots ϕ_{n})}^{⊤}

of length n. Such that,

J (Φ) = (\begin{matrix} \frac{\partial x}{\partial ϕ_{0}} & \frac{\partial x}{\partial ϕ_{1}} & \frac{\partial x}{\partial ϕ_{2}} & \dots & \frac{\partial x}{\partial ϕ_{n}} \\ \frac{\partial y}{\partial ϕ_{0}} & \frac{\partial y}{\partial ϕ_{1}} & \frac{\partial y}{\partial ϕ_{2}} & \dots & \frac{\partial y}{\partial ϕ_{n}} \\ \frac{\partial z}{\partial ϕ_{0}} & \frac{\partial z}{\partial ϕ_{1}} & \frac{\partial z}{\partial ϕ_{2}} & \dots & \frac{\partial z}{\partial ϕ_{n}} \end{matrix}) .

(46)

For the upper limb the Jacobian matrix elements are

\frac{\partial x}{\partial ϕ_{0}} = 0

(47)

\frac{\partial x}{\partial ϕ_{1}} = l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{3} s_{124}

(48)

\frac{\partial x}{\partial ϕ_{2}} = l_{2} s_{12} + l_{3} s_{3} s_{124}

(49)

\frac{\partial x}{\partial ϕ_{3}} = - l_{3} s_{3} s_{124}

(50)

\frac{\partial x}{\partial ϕ_{4}} = l_{3} s_{3} s_{124},

(51)

and for y,

\frac{\partial y}{\partial ϕ_{0}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{124}) + l_{3} s_{3} s_{124} s_{0}

(52)

\frac{\partial y}{\partial ϕ_{1}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{124}) + l_{3} s_{3} s_{124} s_{0}

(53)

\frac{\partial y}{\partial ϕ_{2}} = s_{0} (l_{2} s_{12} + l_{3} s_{124}) + l_{3} s_{3} s_{124} s_{0}

(54)

\frac{\partial y}{\partial ϕ_{3}} = l_{3} s_{3} s_{124} s_{0}

(55)

\frac{\partial y}{\partial ϕ_{4}} = l_{3} s_{124} s_{0} + l_{3} s_{3} s_{124} s_{0},

(56)

and for z,

\frac{\partial z}{\partial ϕ_{0}} = - s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{124}) + l_{3} s_{3} s_{124} s_{0}

(57)

\frac{\partial z}{\partial ϕ_{1}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{2} s_{12} + l_{3} s_{124}) - l_{3} s_{3} s_{124} s_{0}

(58)

\frac{\partial z}{\partial ϕ_{2}} = s_{0} (l_{2} s_{12} + l_{3} s_{124}) - l_{3} s_{3} s_{124} s_{0}

(59)

\frac{\partial z}{\partial ϕ_{3}} = - l_{3} s_{3} s_{124} s_{0}

(60)

\frac{\partial z}{\partial ϕ_{4}} = l_{3} s_{124} s_{0} - l_{3} s_{3} s_{124} s_{0}

(61)

Likewise, for the lower limb, the Jacobian matrix terms are

\frac{\partial x}{\partial ϕ_{0}} = 0; \frac{\partial x}{\partial ϕ_{1}} = l_{1} s_{1} + l_{2} s_{12} - l_{3} s_{123};

(62)

\frac{\partial x}{\partial ϕ_{2}} = l_{2} s_{12} - l_{3} s_{123}; \frac{\partial x}{\partial ϕ_{3}} = - l_{3} s_{123} .

(63)

and for y

\frac{\partial y}{\partial ϕ_{0}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} - l_{3} s_{123})

(64)

\frac{\partial y}{\partial ϕ_{1}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{123})

(65)

\frac{\partial y}{\partial ϕ_{2}} = s_{0} (l_{2} s_{12} + l_{3} s_{123}); \frac{\partial y}{\partial ϕ_{3}} = l_{3} s_{0} s_{123} .

(66)

and for z

\frac{\partial z}{\partial ϕ_{0}} = - s_{0} (l_{1} s_{1} + l_{2} s_{12} - l_{3} s_{123})

(67)

\frac{\partial z}{\partial ϕ_{1}} = s_{0} (l_{1} s_{1} + l_{2} s_{12} + l_{3} s_{123})

(68)

\frac{\partial z}{\partial ϕ_{2}} = s_{0} (l_{2} s_{12} + l_{3} s_{123}); \frac{\partial z}{\partial ϕ_{3}} = l_{3} s_{0} s_{123}

(69)

4.2. Recursive time-variant control law

Other exoskeleton control approaches have addressed nonlinearity using sliding surfaces [29] or adaptive methods, implementing impedance with feedback using biosignals [30], and neural approximation [31].

The present work takes an approach involving a model-based algebraic linear recursive system of equations. This work deduces a time-variant redundant kinematics controller since

J_{t}

is considered. The model incorporates the reference trajectory provided by the therapist, denoted as

p^{r e f}

. The real-time feedback utilizes the EMG observation

{\hat{p}}_{t} = {({\hat{x}}_{t}, {\hat{y}}_{t})}^{⊤}

. The innovation model, representing the local reference, is the inverse model

Φ_{t + 1}

, and the updating model involves the encoders’ observation

{\hat{Φ}}_{t}

Let us define the state linear equation (70) as the first-order derivative forward kinematics.

\dot{p} = J \cdot \dot{Φ}

(70)

and stating this expression in the form of differential equation, its derivative-order decreases

\frac{d p}{d t} = J \cdot \frac{Φ}{d t},

(71)

completing all differentials with their respective integrals at both sides of the equality

\int_{p_{1}}^{p_{2}} d p = J \cdot \int_{Φ_{1}}^{Φ_{2}} d Φ,

(72)

therefore, the expression (73) is the general recursive solution,

p_{t 2} - p_{t 1} = J \cdot (Φ_{t 2} - Φ_{t 1}) .

(73)

Therefore, from such an equality the expressions (74) lead to forward/backward general solutions.

The recursive forward:

p_{t + 1} = p_{t} + J_{t} \cdot (Φ_{t + 1} - Φ_{t})

(74a)

and the recursive backward:

Φ_{t + 1} = Φ_{t} + J_{t}^{+} \cdot (p_{f} - p_{t}) .

(74b)

Therefore, the inverse general model for the recursive control law is set by an ideal global trajectory

p_{t}^{r e f}

Φ_{t + 1} = Φ_{t} + J_{t}^{+} (p_{t}^{r e f} - {\hat{p}}_{t}) .

(75)

Some terms in the model are expanded. Since the Jacobian is non-squared, a pseudoinverse matrix solution is applied,

Φ_{t + 1} = Φ_{t} + {(J_{t} J_{t}^{⊤})}^{- 1} J_{t}^{⊤} [(\begin{matrix} x_{t}^{r e f} \\ y_{t}^{r e f} \end{matrix}) - (\begin{matrix} \hat{x} (\hat{t}) \\ \hat{y} (\hat{t}) \end{matrix})]

(76)

and the forward solution has the next local reference

Φ_{t + 1}

p_{t + 1} = p_{t} + J_{t} (Φ_{t + 1} - {\hat{Φ}}_{t}) .

(77)

Where the exoskeleton’s joints direct measurements are obtained from encoder devices,

p_{t + 1} = p_{t} + J_{t} (Φ_{t + 1} - \frac{2 π}{R} (\begin{matrix} η_{0} \\ η_{1} \\ ⋮ \\ η_{k} \end{matrix})) .

(78)

When the error condition

∥ p_{f} - p_{t} ∥ > ϵ

is met, the control loop achieves goal tracking. Figure 19a illustrates a complex trajectory and its tracking control performance concerning the joint angles (Figure 19b). Additionally, Figure 19c depicts the evolution of position metric errors throughout the entire reference trajectory.

In the context of numerical simulations for numerous controlled trajectories, the efficacy of the control approach was validated through a comparative analysis against a reference ideal function. The trajectory in Figure 20-above illustrates an

x y

path with an arbitrary reference defined as

y (x) = k_{5} x^{5} + k_{3} x^{3} - k_{0}

. Along the same row of the figure, the plots elucidate the error between the reference function and the controlled function. Ideally, a precise match is anticipated between the two functions. The observed errors in the x and y components, in general, are relatively small. As for the z component, its error is nearly zero since this trajectory is intentionally planar, confined solely to the

x y

plane. Consequently, any variation along the z axis should be virtually null. This illustration underscores the expectation that the controlled functions accurately align with a segment of a 5^th-grade polynomial function.

Moreover, Figure 20-below displays the resulting errors in joints’ angles induced by the trajectory in the same figure (expressed in physical units as radians).

A comprehensive depiction akin to Figure 20, showcasing the entire joint performance, is presented in Figure 21. A distinct range of coordinates from a different reference function is established as the reference trajectory, with the framework’s origin situated at the shoulder centroid. Similarly, Figure 21-above exhibits a Cartesian plot along with the errors along the three axes (reference versus controlled function). Figure 21-below illustrates the joints’ angle performance resulting from the controlled tracking trajectory.

5. Results and discussion

Summarizing the core focus of this work, we direct our attention to the following three figures: Figure 22, Figure 23, and Figure 24. These figures encapsulate the fundamental purpose and achievement of this study. They provide a visual representation of the work’s primary objective, which is the estimation of joint angles from EMG measurements. For the sake of clarity and simplicity, the plots in these figures display data for two key joint angles, highlighting the remarkable success of the approach. Figure 22 presents data from three distinct experiments involving elbow flexion. Within this figure, the shoulder angle is denoted as

ϕ_{1}

(depicted in red), and the elbow angle is represented as

ϕ_{2}

(in blue). Two specific channels, Ch5 and Ch6, correspond to the electromyography measurements of the biceps brachii and triceps brachii, respectively. The selection of these two channels reflects the authors’ recognition of their suitability for assessing electromyography activity during limb motion.

Following the calibration process, the results in Figure 22 demonstrate that the activity associated with elbow flexion (depicted in blue) reveals a strong correlation between the elevation angle and the measurements from Ch5, emphasizing the crucial role of the biceps brachii in this motion. Conversely, the decrease in the elbow’s angle is found to be correlated with the measurements from Ch6, highlighting the impact of the triceps brachii. Interestingly, the motion of the shoulder joint (

ϕ_{1}

) is not detected within this specific configuration, as indicated by the red-colored plot. This observation sheds light on the distinct contributions of the two channels to the control of different joint angles.

Figure 22. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in elbow-flexion performance.

Figure 22. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in elbow-flexion performance.

Figure 23 offers a close look at an experiment comprising five repetitions of shoulder flexion. Within this figure, the shoulder angle continues to be represented as

ϕ_{1}

, depicted in red, while the elbow angle is identified as

ϕ_{2}

and shown in blue. The electromyography measurements are captured by channels Ch5 and Ch6, corresponding to the activity of the pectoralis major and anterior pectoralis muscles, respectively. Following the calibration process, the results in Figure 23 reveal that during shoulder flexion activity (as depicted in the red-colored plot), the elevation angle is strongly influenced, with substantial contributions from both Ch6 and Ch5 measurements. This highlights the role of both the anterior pectoralis and the pectoralis major muscles in controlling the shoulder’s elevation angle.

However, as the exercise primarily involves shoulder motion and maintains the elbow at a fixed position, the motion of joint

ϕ_{2}

is not detected in this particular configuration, as indicated by the blue-colored plot. This distinction further underscores the remarkable accuracy of the approach in distinguishing the contributions of different muscle groups to specific joint movements. Overall, Figure 23 demonstrates the successful application of electromyography measurements in estimating joint angles, exemplifying the adaptability and precision of this approach for diverse exercises involving multiple joints.

Figure 23. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in shoulder-flexion performance.

Figure 23. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in shoulder-flexion performance.

Figure 24 offers insight into three distinct experiments involving rowing-like flexion, wherein both shoulder and elbow joints exhibit synchronized motion activity. As in previous figures, the shoulder angle is denoted as

ϕ_{1}

, represented by the red-colored plot, while the elbow angle is identified as

ϕ_{2}

and displayed in blue. The electromyography measurements for these experiments are derived from channels Ch5 and Ch6, corresponding to the activity of the middle trapezius and anterior deltoid muscles, respectively. After calibration, the results within Figure 24 reveal that elbow flexion activity (illustrated by the blue-colored plot) exhibits different oscillatory angular phases. The data suggests that the elevation of

ϕ_{1}

appears to be significantly impacted by Ch5 measurements, while

ϕ_{2}

is influenced by Ch6 data. The remarkable consistency observed across these three experiments serves to validate the proposed approach in this research, showcasing the method’s ability to accurately estimate joint angles during complex exercises that involve synchronized motion of multiple joints. This validation is a testament to the reliability and effectiveness of the approach, underscoring its potential for real-world applications in the field of tele-rehabilitation.

Figure 24. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in rowing-like flexions performance.

Figure 24. Calibrated EMG measurements, estimation joints’ angle

ϕ_{1}

and

ϕ_{2}

in rowing-like flexions performance.

Table 3. Physics-based computer simulation parameters^a.

Limbs:
m=1 kg	$I_{x} = 0.1666$ kg $m^{2}$	$I_{y} = 0.1666$ kg $m^{2}$	$I_{z} = 0.1666$ kg $m^{2}$
$ł_{1, 2} = 0.3$ m	limbs $r_{1, 2} = 0.15$ m	$τ_{r e f} = 19.62$ Nm	${\ddot{ϕ}}_{r e f} = 36.14$ $s^{- 2}$
Gears:
$r_{1} = 1.5 r_{0}$ m	$r_{2, 3} = 2.6 r_{0}$ m	$r_{4, 5} = 2 r_{0}$ m	$r_{6, 7} = 0.66 r_{0}$ m
$r_{8} = 1.16 r_{0}$ m	$n_{6, 7} = 2$ (threads)	$n_{8} = 40$ (teeth)	$m_{8} = 1.166 m_{0}$ kg
$m_{1} = 1.5 m_{0}$ kg	$m_{2, 3} = 3 m_{0}$ kg	$m_{4, 5} = 2 m_{0}$ kg	$m_{6, 7} = 1.5 m_{0}$ kg
$I_{1} = 3.37 r_{0}^{2} m_{0}$ kg $m^{2}$	$I_{2, 3} = 1.72 r_{0}^{2} m_{0}$ kg $m^{2}$	$I_{4, 5} = 8 r_{0}^{2} m_{0}$ kg $m^{2}$	$I_{6, 7} = 0.2875 r_{0}^{2} m_{0}$ kg $m^{2}$
$I_{8} = 1.58 r_{0}^{2} m_{0}$ kg $m^{2}$
^aGazebo enabled with Object Dynamic Object (ODE).

Although this study does not address more intricate rehabilitation types, the showcased exercises are representative of typical rehabilitation routines, serving as a testament to the feasibility of the proposed method. Currently, the demonstration of the model’s applicability is achieved using only two EMG channels. Given the limitation of relying on a singular visual sensor, the integration of various technological sensing devices becomes imperative. This broader array of sensing devices aims to encompass diverse perspectives surrounding the therapist’s tasks, thereby enhancing the overall accuracy of measurements. In the context of trajectory adjustment for the linear trajectory smoother model, a manual trial-and-error approach is employed to fine-tune the trajectory parameters, specifically

k_{t}

and

k_{s}

6. Conclusion

In this comprehensive research endeavor, a remarkable paradigm shift was introduced in the realm of tele-rehabilitation. The proposed approach has a feature: allowed for asynchronous interactions, removing the constraints of rigid scheduling, as exercises could be conducted at different dates and times, making rehabilitation more accessible and flexible for patients. The patient-exoskeleton interface was engineered to receive concise data structures, and the innovative use of coefficient sets emerged as a masterstroke. These sets proved invaluable in ensuring that data streaming remained swift and uninterrupted, unfazed by the passage of time. One of the standout features of this work was the treatment of all exercises as periodic functions. This approach enabled the seamless application of Fourier series fitting, providing an elegant and efficient means of data representation. The ability to translate complex exercises into Fourier series representations presented a new dimension of flexibility in tracking and control.

A cornerstone of this research was the method for reproducing rehabilitation exercises using Cartesian analytical polynomials. This method unlocked the door to producing high-resolution functions of time by synthesizing Cartesian functions of time with EMG data, with time itself serving as the universal synchronization variable. This synthesis represented the union of clinical precision and technological sophistication, making the rehabilitation process more accurate and efficient. A novel model for inferring Cartesian posture from EMG signals was an important issue of this work, showcasing remarkable reliability. The need for recalibration with each new exercise served as a testament to the adaptability and responsiveness of this approach, ensuring that the exoskeleton remained in sync with the patient’s unique characteristics and movements.

The proposed exoskeleton underwent rigorous simulation testing, encompassing hundreds of tracking control rehabilitation exercises. Its adaptability was underscored by a control law that could be seamlessly adjusted to suit various exoskeleton structures, simply by readjusting the Jacobian matrix. This adaptability introduced an exciting degree of versatility and expandability in the application of this technology, not limited by specific hardware constraints. In essence, this simulation prototype, anchored by a single reference trajectory, demonstrated remarkable potential in the calibration of numerous EMG experiments and the execution of multiple tracking control exercises. The system’s adaptability, efficiency, and the profound impact it could have on the field of tele-rehabilitation were evident throughout this research. The approach outlined in this study promises to redefine and enhance the way rehabilitation is conducted, making it not just more effective but also more accessible and accommodating to a wide range of patient needs. It represents a significant step forward in the integration of technology and healthcare, and its implications are bound to be felt for years to come.

As a future work, enhancing the visual teaching process can be accomplished by substituting the RGB camera with either an RGB-D or a ranging laser device. This alternative not only mitigates noise but is also impervious to constraints and illumination variations, offering the added advantages of infrared capabilities and 3D cloud data provision. In addition, to ensure the linear trajectory smoother model’s optimal performance, automatic optimization methods and tuning of rehabilitation trajectories, facilitated by parameters

k_{t}

and

k_{s}

, will enable a systematic convergence toward the desired reference trajectory shape. In our forthcoming work, there is an envisioned redesign of the exoskeleton mechanism. This redesign incorporates elastic soft materials, mimicking artificial muscles activated by electrical/optical signals, preserving EMG feedback. In addition, addressing the expansion of the EMG device array is pivotal, allowing for the comprehensive mapping of additional body muscles while utilizing more compact technological devices. We posit that exploring alternative controllers is imperative for achieving enhanced performance, tailored to different rehabilitation tasks, fostering a more nuanced understanding of optimal controller-task relationships.

Data Availability Statement

Experimental data obtained during the development of this research is accessible at, Martinez-Garcia, E. A., & Castro Jiménez, U. (2024). EMG model calibration and trajectory tracking by remote rehabilitation exoskeleton [Data set]. Zenodo. https://doi.org/10.5281/zenodo.10525565.

Acknowledgments

The authors are grateful for the technical guidance provided by the Robotics Laboratory team throughout the implementation of specific developments. Our sincere appreciation also goes to the involved professors for their invaluable suggestions, which played a crucial role in the remarkable success of this project.

References

Evgeni Magid, Aufar Zakiev, Tatyana Tsoy, Roman Lavrenov, Albert Rizvanov. 2021. Automating pandemic mitigation, Advanced Robotics, 35(9), 572-589. [CrossRef]
G. Yang, Z. Pang, M. J. Deen, M. Dong, Y. Zhang, N. Lovell, A.M. Rahmani. 2020. Homecare robotic systems for healthcare 4.0: visions and enabling technologies, IEEE J. of Biom & Health Informatics, 24(9):2535–2549. [CrossRef]
I. Haider, K.B. Khan, M.A. Haider, A. Saeed, K. Nisar. 2020. Automated Robotic System for Assistance of Isolated Patients of Coronavirus (COVID-19), IEEE 23rd Intl. Multitopic Conf., 5-7 Nov.
Z. Wu, L. Chitkushev, G. Zhang. 2020. A Review of Telemedicine in time of COVID-19, IEEE Intl. Conf. on Bioinformatics and Biomedicine, Seoul, Korea.
M. Touil, L. Bahatti, A.E. Magri. 2020. Telemedicine application to reduce the spread of Covid-19, In: IEEE 2nd Intl. Conf. on Electronics, Control, Optimization and Computer Science, Kenitra, Morocco, Dec.
J.R. Barr, D. D’Auria, F. Persia. 2020. Telemedicine, Homecare in the Era of COVID-19 & Beyond, In: Intl. Conf. on artificial intelligence for industries, CA, USA.
B. Samadi, S. Achiche, A. Parent, L. Ballaz, U. Chouinard, M. Raison. 2015. Custom sizing of lower limb exoskeleton actuators using gait dynamic modelling of children with cerebral palsy, Computer Methods in Biomech & Biomed Eng, 19(14):1519-1524. [CrossRef]
Z. Li et al. 2019. Hybrid brain/muscle signals powered wearable walking exoskeleton enhancing motor ability in climbing stairs activity, IEEE Trans. on Medical Robotics and Bionics, 1(4):218–227. [CrossRef]
I. Kang, P. Kunapuli, H. Hsu and A. J. Young. 2019. Electromyography (EMG) signal contributions in speed and slope estimation using robotic exoskeletons, In: IEEE 16th Intl. Conf. on Rehabilitation Robotics, Toronto, Canada, pp.548–553. [CrossRef]
D. Buongiorno, E. Sotgiu, D. Leonardis, S. Marcheschi, M. Solazzi and A. Frisoli. 2018. WRES: A novel 3 DoF wrist exoskeleton with tendon-driven differential transmission for neuro-rehabilitation and teleoperation, IEEE Rob and Autom Letters, 3(3):2152–2159. [CrossRef]
S. Pu, Y. Pei, J. Chang. 2020. Decoupling finger joint motion in an exoskeletal hand: a design for robot-assisted rehabilitation, IEEE Trans. on Ind Elec, 67(1)–686-697. [CrossRef]
Pons JL (Ed). 2008. Wearable robots: Biomechatronic Exoskeletons, 1st edn. Wiley & Sons.
C.J. Yang, J.F. Zhang, Y. Chen, Y.M. Dong, Y. Zhang. 2008. A Review of exoskeleton-type systems and their key technologies, Proc. of the Inst. of Mech. Engineers, Part C: J. of Mech Eng Sci, 222(8):1599–1612.
K. H. Low Xiaopeng Liu C. H. Goh Haoyong Yu. 2006. Locomotive control of a wearable lower exoskeleton for walking enhancement, J. of Vibration and Control, 12(12):1311–1336.
J. Charafeddine, D. Pradon, S. Alfayad, S. Chevallier, M. Khalil. 2019. Toward bio-kinematic for secure use of rehabilitation exoskeleton, Computer Methods in Biomechcanics and Biomedical Engineering, 22(sup1):s424-s426.
J.A. Díez, J.M. Catalán, L.D. Lledó, F. J. Badesa, N. Garcia-Aracil. 2016. Multimodal robotic system for upper-limb rehabilitation in physical environment, Adv in Mech Eng, 8(9).
R. Jiménez, O. Verlinden, Review of control algorithms for robotic ankle systems in lower-limb orthoses, prostheses and exoskeletons, J. of Med Eng & Physics, Elsevier.
G.S. Athwal. 2017. Shoulder Surgery Exercise Guide, American Academy of Orthopaedic Surgeons, USA, 2017. [Online] Available: https://orthoinfo.aaos.org/en/recovery/shoulder-surgery-exercise-guide/.
J. Lanini, T. Tsuji, P. Wolf, R. Riener and D. Novak. 2015. Teleoperation of two six-degree-of-freedom arm rehabilitation exoskeletons, In: IEEE Intl. Conf. on Rehabilitation Robotics, Singapore, 2015, pp. 514-519.
W. Sun et al., Research on Facial Emotion Recognition System Based on Exoskeleton Rehabilitation Robot, IEEE 11th Intl. Conf. on Software Engineering and Service Science, Beijing, China, 2020, pp. 481-484.
X. Li, Y. Pan, G. Chen, H. Yu. 2017. Multi-modal control scheme for rehabilitation robotic exoskeletons, The Intl. J. of Robot Res, 36(5-7):759–777.
C. Fisahn, M. Aach, O. Jansen, M. Moisi, A. Mayadev, K.T. Pagarigan, J.R. Dettori, T.A. Schildhauer. 2016. The Effectiveness and Safety of Exoskeletons as Assistive and Rehabilitation Devices in the Treatment of Neurologic Gait Disorders in Patients with Spinal Cord Injury: A Systematic Review, Global Spine Journal, 6(8):822–841.
Z. Li, B. Wang, F. Sun, C. Yang, Q. Xie and W. Zhang. 2014. sEMG-Based Joint Force Control for an Upper-Limb Power-Assist Exoskeleton Robot, IEEE J. of Biomedical and Health Informatics, 18(3):1043–1050.
G. Bauer and Y.J. Pan. 2020. Review of Control Methods for Upper Limb Telerehabilitation With Robotic Exoskeletons, IEEE Access, 2020(8):203382–203397.
M. V. Arteaga, J. C. Castiblanco, I. F. Mondragon, J. D. Colorado and C. Alvarado-Rojas. 2020. EMG-based adaptive trajectory generation for an exoskeleton model during hand rehabilitation exercises, In: 8th IEEE RAS/EMBS Intl. Conf. for Biomedical Robotics and Biomechatronics, New York, USA, 416–421.
H. Wang, X. Shi, H. Liu, L. Li, Z. Hou, Yu. 2011. Design, kinematics, simulation, and experiment for a lower-limb rehabilitation robot, Proc. of the IMEchE Part I: J. Sys and Cont Eng, 225(6):860–872. [CrossRef]
G. Aguirre-Ollinger, J.E. Colgate, M.A. Peshkin, A. Goswami. 2011. A One-Degree-of-Freedom Assistive Exoskeleton with Inertia Compensation: The Effects on the Agility of Leg Swing Motion, Proc. of the Inst.of Mech Eng, Part H: J. of Eng in Med, 225(3):228–245.
Wei He, Shuzhi Sam Ge, Yanan Li, Effie Chew, Yee Sien Ng. 2015. Neural Network Control of a Rehabilitation Robot by State and Output Feedback, J. of Intel and Robot Sys, vol. 80, p.15-31.
S. Zhu, X. Jin, B. Yao, Q. Chen, X. Pei, Z. Pan. 2016. Non-linear sliding mode control of the lower extremity exoskeleton based on human–robot cooperation, Intl J of Adv Rob Sys, 13(5).
Z. Li, Z. Huang, W. He, C. Y. Su. 2017. Adaptive Impedance Control for an upper limb robotic exoskeleton using biological signals, IEEE Trans. on Ind Elect, 64(2), pp. 1664-1674.
W. He, Z. Li, Y. Dong, T. Zhao. 2019. Design and Adaptive Control for an Upper Limb Robotic Exoskeleton in Presence of Input Saturation, IEEE Trans. on Neural Net and Learn Sys, 30(1), pp. 97-108. [CrossRef]
M. Gandolla, S. Ferrante, G. Ferrigno, D. Baldassini, F. Molteni, E. Guanziroli, M. Cotti Cottini, C. Seneci, A. Pedrocchi. 2016. Artificial neural network EMG classifier for functional hand grasp movements prediction, J. of Intl. Medical Research, 45(6):1831-1847.
T. Alenabi, M. Jackson, P. Tétreault, M. Begon, “Electromyographic activity in the shoulder musculature during resistance training exercises of the ipsilateral upper limb while wearing a shoulder orthosis”, J. of Shoulder and Elbow Surgery, vol.23, pp. 140-148, 2013.
M.Sarac, M. Solazzi, A. Frisoli, “Design requirements of generic hand exoskeletons and survey of hand exoskeletons for rehabilitation, assistive, or haptic use”, IEEE Trans. on Haptics, 12(4), pp. 400-413, 2019.
D. Wang, K. Lee, J. Guo, C. Yang, Adaptive Knee Joint Exoskeleton Based on Biological Geometries, IEEE/ASME Trans. on Mechatronics, 19(4), pp. 1268-1278, 2014.
Y. Liu, S. Guo, Z. Yang, H. Hirata, T. Tamiya, A Home-based Tele-rehabilitation System With Enhanced Therapist-patient Remote Interaction: A Feasibility Study, IEEE J. of Biomedical and Health Informatics, 26(8), p. 4176-4186, 2022.
Y. Bouteraa, I. B. Abdallah, Exoskeleton robots for upper-limb rehabilitation, IEEE 13th Intl. Multi-Conf. on Sys, Signals and Devices, 2016. [CrossRef]
M. Kim, H. Jeong, P. Kantharaju, D. Yoo, M. Jacobson, D. Shin, C. Han, J. L. Patton, Visual guidance can help with the use of a robotic exoskeleton during human walking, Nature Scientific Reports, 12:3881, 2022.
A. Akbari, F. Haghverd, S. Behbahani, Robotic Home-Based Rehabilitation Systems Design: From a Literature Review to a Conceptual Framework for Community-Based Remote Therapy During COVID-19 Pandemic, frontiers in Robotics and AI, 8:612331, 2021.
A. Cisnal, P. Gordaliza, J. Pérez Turiel, J. C. Fraile, Interaction with a Hand Rehabilitation Exoskeleton in EMG-Driven Bilateral Therapy: Influence of Visual Biofeedback on the Users’ Performance, Sensors, 23, 2048, MDPI, 2023.
S. Zhang, Q. Fu, S. Guo, Y. Fu, A Telepresence System for Therapist-in-the-Loop Training for Elbow Joint Rehabilitation, Applied science, 9, 1710, mdpi, 2019.

1	Experimental estimations $f_{h} = 1.10249$ m, $f_{w} = 1.48499$ m.
2	The Moore-Penrose right-pseudoinverse is applied $m^{+} = m^{⊤} {(m m^{⊤})}^{- 1}$

Figure 4. Experimental measurements in flexion-extension exercises of shoulder, elbow and wrist.

Figure 5. Angular measurements for flexion-extension exercises of shoulder, elbow and wrist.

Figure 6.

y (x)

models fitting real trajectories. Depiction of

x y

positions and joints’ angle.

Figure 6.

y (x)

models fitting real trajectories. Depiction of

x y

positions and joints’ angle.

Figure 7. Analytical

x (t)

y (t)

functions fitting elbow-flexion experimental data of Figure 4.

Figure 7. Analytical

x (t)

y (t)

functions fitting elbow-flexion experimental data of Figure 4.

Figure 8. Family of smoothed reference trajectories adjusted by

κ_{τ}

and

κ_{s}

Figure 8. Family of smoothed reference trajectories adjusted by

κ_{τ}

and

κ_{s}

Figure 10. The upper limbs’ EMG electrode placement is delineated in Channel 1 (left) and Channel 2 (right).

Figure 11. Three elbow-flexion EMG measurements. Channels are Ch5 (bicep braquial) and Ch6 (tricep braquial).

Figure 12. Shoulder flexion EMG measurements. Channels Ch5 are (Pectoralis Major) and Ch6 (Deltoid Anterior).

Figure 15.

E (t)

polynomial fitting with

i^{t h}

coefficients by McLaurin series.

Figure 15.

E (t)

polynomial fitting with

i^{t h}

coefficients by McLaurin series.

Figure 17. The joint’s mechanism.

Figure 19. Tracking control of a complex reference trajectory.

Figure 20. Trajectory tracking metric errors, reference

p^{r e f} - p_{t}

, tracking

y (x) = k_{5} x^{5} + k_{3} x^{3} - k_{0}

and joints’ error.

Figure 20. Trajectory tracking metric errors, reference

p^{r e f} - p_{t}

, tracking

y (x) = k_{5} x^{5} + k_{3} x^{3} - k_{0}

and joints’ error.

Figure 21. Metric errors of controlled planar trajectory versus a reference (

p^{r e f} - p_{t}

). Tracking

y (x) = κ_{1} cos (\frac{κ_{2} x}{x + κ_{3}})

Figure 21. Metric errors of controlled planar trajectory versus a reference (

p^{r e f} - p_{t}

). Tracking

y (x) = κ_{1} cos (\frac{κ_{2} x}{x + κ_{3}})

Table 1. Experimental visual sensor’s intrinsic/extrinsic parameters^a.

Resolution	$640 \times 480$	height	$1.1024982$ m
Focal lenght	1.5 m	width	$1.4849993$ m
Frames/s	30	z (distance)	2.0 m
$φ_{h o r i z o n t a l}$	$52 . 6708^{o}$	$φ_{v e r t i c a l}$	$40 . 3567^{o}$
H. resol $f_{h}$	1.10249 m	V. resol. $f_{v}$	1.48499 m

^aLogitech USB 2.0 camera, views of Figure 3.

Table 2. EMG electrodes placement for different therapy.

Exercise^a	Channel 1	Channel 2
Elbow	Biceps brachii	Triceps brachii
Shoulder	Pectoralis major	Anterior pectoralis
Shoulder and elbow	Middle trapezius	Anterior deltoid
Knee	Rectur femoris	Vastus medialis

^aFlexion-extension exercise. See Figure 10 for channels placement.

MOBIlab+ sampling 256Hz, high/low pass filters

100 / 0.5

Hz.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Alerts

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

EMG Model Calibration and Trajectory Tracking by Remote Rehabilitation Exoskeleton

Abstract

Keywords:

Subject:

1. Introduction

2. Human therapist remote visual teaching

3. Rehabilitation trajectory model

3.1. EMG sensing model calibration

3.2. Observability based on real-time EMG

4. Exoskeletons morphology and control law

4.1. Limbs kinematic model

4.2. Recursive time-variant control law

5. Results and discussion

6. Conclusion

Data Availability Statement

Acknowledgments

References

MDPI Initiatives

Important Links

Subscribe