Open access peer-reviewed chapter

Trajectory Tracking Using Adaptive Fractional PID Control of Biped Robots with Time-Delay Feedback

By Joel Perez Padron, Jose P. Perez, C.F. Mendez-Barrios and V. Ramírez-Rivera

Submitted: February 24th 2018Reviewed: December 19th 2019Published: March 25th 2020

DOI: 10.5772/intechopen.90901

Downloaded: 153

Abstract

This paper presents the application of fractional order time-delay adaptive neural networks to the trajectory tracking for chaos synchronization between Fractional Order delayed plant, reference and fractional order time-delay adaptive neural networks. For this purpose, we obtained two control laws and laws of adaptive weights online, obtained using the fractional order Lyapunov-Krasovskii stability analysis methodology. The main methodologies, on which the approach is based, are fractional order PID the fractional order Lyapunov-Krasovskii functions methodology, although the results we obtain are applied to a wide class of non-linear systems, we will apply it in this chapter to a bipedal robot. The structure of the biped robot is designed with two degrees of freedom per leg, corresponding to the knee and hip joints. Since torso and ankle are not considered, it is obtained a 4-DOF system, and each leg, we try to force this biped robot to track a reference signal given by undamped Duffing equation. To verify the analytical results, an example of dynamical network is simulated, and two theorems are proposed to ensure the tracking of the nonlinear system. The tracking error is globally asymptotically stabilized by two control laws derived based on a Lyapunov-Krasovskii functional.

Keywords

  • biped robot
  • fractional time-delay adaptive neural networks
  • fractional order PID control
  • fractional Lyapunov-Krasovskii functions
  • trajectory tracking

1. Introduction

Fractional calculus is a generalization of differential and integral calculus which involves generalized functions. The first to work this new branch of mathematics was Leibniz. Due to the growing interest in the applications of fractional calculation, in this work we obtain conditions that guarantee the tracking of trajectories of nonlinear systems generated by differential equations of fractional order which we will call plants (This term is widely used in engineering), which in our case will be a mechanical arm, a helicopter, a plane or limbs of a humanoid, all of fractional order.

The problem of tracking control of trajectories is very important, since the control function allows the non-linear system to carry out a previously assigned task, work or trajectory, for example, a mechanical arm and its objective is to cut a piece with a previously generated form, or the coupling of two aircraft in space.

In this chapter we use adaptive recurrent neural networks with time delay, since its use allows us to work with systems whose mathematical model is unknown and with the presence of uncertainties, this is a well-known problem of robust control.

We include mathematical models with time delay, since the processing and transmission of information is important in this type of systems, which depending on the delay, these systems can generate undesirable oscillatory or chaotic dynamics, and cause instability in the mathematical model that describes the trajectory tracking error.

The chapter is organized as follows: first, the general mathematical model of non-linear systems is proposed, as a second part, the Neural Network is proposed that will adapt to the non-linear system and the reference signal that both must follow, as a third part obtains the dynamics of the tracking error between the non-linear system and the reference, after obtaining conditions in the laws of adaptation of weights in the Neural Network and obtaining the control law that guarantees that the tracking error converges to zero, so that the non-linear system will follow the indicated reference signal, which is what was wanted to be demonstrated. Finally simulations are presented, which illustrate the theoretical results previously demonstrated. The proposed new control scheme is applied via simulations to control of a 4-DOF Biped Robot [1].

We use the scheme of Figure 1 to indicate the procedure used in the obtaining of the laws of adaptation of weights and the laws of control that guarantee that the tracking error between the non-linear system, the neural network and the reference signal converges to zero.

Figure 1.

Adaptive recurrent control diagram.

2. Time-delay adaptive neural network and the reference

2.1 List of variables

Wis the matrix weights.

e=xpxr, error between the plant and the reference.

Ŵis part of the approach, given by W.

Ωu1, Ωu2, up, unare the controls

PIλDα=Kpet+KiaDtλet+KdaDtαetcontrol law

τis time delay

traDtαWTW=eTWσxtτ, learning law from the neural network weights

tτtØσTsŴTŴØσsds, Lyapunov-Krasovskii Function

Dqtq¨t+Cqtq̇tq̇t+Gqt=t, dynamics of the bipedal robot

D(qt) is the inertia matrix

Cqtq̇tis the matrix of Coriolis and centripetal forces

Gqtrepresents a matrix of gravitational effects

Bdefines the input matrix

There are several ways to define the fractional calculation, in this research we will use the well-known derivative of Caputo, which has the following notation:

aDtαft=1Γαnatfnτtταn+1E1
Forn1<α<n.

The nonlinear system, Eq. (2), which is forced to follow a reference signal:

aDtαxp=fptxpt+xptτ,t0T,0<α1,E2
xpt=gt
xp,fpRn,uRm,gpRnxn.

The differential equation will be modeled by the neural network [2]:

aDtαxp=Ax+WΓzx(tτ+Ωu.

The tracking error between these two systems:

wper=xxpE3

We use the next hypotheses.

aDtαwper=kwperE4

In this research we will use k=1, so that, Eq. (5), aDtαwper=aDtαxaDtαxp, so:

aDtαxp=aDtαx+wper

The nonlinear system is [3]:

aDtαxp=aDtαx+wper=Ax+WΓzxtτ+wper+ΩuE5

where the Wis the matrix weights.

3. Tracking error problem

In this part, we will analyze the trajectory tracking problem generated by

aDtαxr=frxrur,wr,xrRnE6

are the state space vector, input vector and fr, is a nonlinear vectorial function.

To achieve our goal of trajectory tracking, we propose the error between the plant and the reference as: e=xpxr=xpxn+xnxr=xpx+xxr.

Let ep=xpx, and en=xxr, be the trajectory tracking error and e=ep+er

en=xxrE7

The time derivative of the error is:

aDtαen=aDtαxaDtαxr=Ax+WΓzxtτ+wper+ΩufrxrurE8

Eq. (8), can be rewritten as follows, adding and subtracting, the next terms ŴΓzxrxtτ, αrtŴ, Aeand wper=xxp, then,

aDtαe=Ax+WΓzxtτ+xxp+Ωufrxrur+ŴΓzxrtτ
ŴΓzxrtτ+ΩαrtŴΩαrtŴ+AeAe
aDtαe=Ae+WΓzxtτ+Ωufrxrur+ŴΓzxrtτE9
+ΩαrtŴΩαrtŴexrAe+x+AxŴΓzxrtτ

The unknown plant will follow the fractional order reference signal, if:

Axr+ŴΓzxrtτ+xrxp+ΩαrtŴ=frxrur,

where

ΩαrtŴ=frxrurAxrŴΓzxrtτxr+xpE10
aDtαe=Ae+WΓzxtτŴΓzxrtτAe+A+Ixxr
+ΩuαrtŴE11

Now, Ŵis part of the approach, given by W. Eq. (11) can be expressed as Eq. (12), adding and subtracting the term ŴΓzxtτand if Γzxtτ=Γzxtτzxrtτ

aDtαe=Ae+WŴΓzxtτ+ŴΓzxtτzxrtτ
+A+IxxrAe+ΩuαrtŴE12

If

W=WŴandu=uαrtŴE13

And by replacing Eq. (13) in Eq. (12), we have:

aDtαe=Ae+WΓzxtτ+ŴΓzxtτzxrtτ
+A+IxxrAe+Ωu
aDtαe=Ae+W˜Γz(x(tτ))+ŴΓ(z(x(tτ))z(xp(tτ))+z(xp(tτ))z(xr(tτ)))+(A+I)(xxp+xpxr)Ae+Ωu˜E14

And:

u=u1+u2E15

So, the result for Ωu1is

Ωu1=ŴΓzxtτzxptτA+IxxpE16

and Eq. (14), is simplified:

aDtαe=Ae+WΓzxtτ+ŴΓzxptτzxrtτ
+A+IxpxrAe+Ωu

Taking into account that e=xpxr, shortening notation a little bit by setting σ=Γz, and defining Øσtτ=σxptτσxrtτ, the equation for aDtαeis

aDtαe=A+Ie+Wσxtτ+ŴØσtτ+Ωu2E17

Now, the problem is to find the control law Ωu2, which it stabilizes to the system Eq. (20). The control law, we will obtain using the fractional order Lyapunov-Krasovskii methodology.

4. Study of trajectory tracking error

Our mathematical model of the dynamics in the tracking error is described in (17). In this equation we can see that an equilibrium state of this system is eŴ=0.

Without loss of generality we can assume that the matrix Ais given A=λI,λ>0, where Iis the identity matrix of order nxn.

For the study of the stability of the tracking error we propose the following PID control law [4], widely used in science and engineering.

We will determine conditions in the parameters that guarantee that the tracking error converges to zero, and we will also use the following control law [5].

Ωu2=Kpe+KiaDtαe+KvaDtαeγ12+12Ŵ2Lϕ2eE18

We also include the following control law, PIλDα[6]:

ut=Kpet+KiaDtλet+KdaDtαet

Substituting Eq. (18) in Eq. (17):

aDtαe=A+Ie+Wσxtτ+ŴØσtτ
+Kpe+KiaDtαe+KvaDtαeγ12+12Ŵ2Lϕ2e,

then

1KvaDtαe=A+Ie+Wσxtτ+ŴØσtτ+Kpe+KiaDtαeγ12+12Ŵ2Lϕ2e.

If a=1Kv, then

aDtαe=1aA+Ie+1aWσxtτ+1aŴØσtτ+1aKpe+1aKiaDtαe
γa12+12Ŵ2Lϕ2eE19
aDtαe=1aλ1+Kpe+1aWσxtτ+1aŴØσtτ+1aKiaDtαe
γa12+12Ŵ2Lϕ2eE20

And if w=1aKiaDtαe, then aDtαw=1aKiet, [7], then Eq. (20) we rewrite as:

aDtαen=1aλ1+Kpe+1aWσxtτ+1aŴØσtτ+w
γa12+12Ŵ2Lϕ2eE21

We will show, the new state enwTis asymptotically stable, and the equilibrium point is enewT=00T, when Wσxrtτ=0, as an external disturbance.

Let Vbe, the next candidate Lyapunov function as [8, 9]:

V=12enTwTenwT+12atrWTWE22
+1atτtØσTsŴTŴØσsds

The fractional order time derivative of (22) along the trajectories of Eq. (21) is:

aDtαV=eTaDtαe+wTaDtαw+1atraDtαWTW+1aØσTtŴTŴØσtØσTtτŴTŴØσtτE23
aDtαV=eT(1a(λ1+Kp)e+1aW˜σ(x(tτ))+1aŴØσ(tτ)+wγa(12+12Ŵ2Lϕ2)e)+1aW˜TKie+1atr{aDtαW˜TW˜}+1a[ØσT(t)ŴTŴØσ(t)ØσT(tτ)ŴTŴØσ(tτ)]E24

In this part, we select the next learning law from the neural network weights as in [10, 11]:

traDtαWTW=eTWσxtτE25

Then Eq. (24) is reduced to

aDtαV=1aλ1+KpeTe+eTaŴØσtτ+1+KiaeTwγa12+12Ŵ2Lϕ2eTe+1aØσTtŴTŴØσtØσTtτŴTŴØσtτE26

Next, let us consider the following inequality proved in [12]

XTY+YTXXTΛX+YTΛ1YE27

Which holds for all matrices X,YRnxkand ΛRnxnwith Λ=ΛT>0. Applying (27) with Λ=Ito the term eTaŴØσtτfrom Eq. (26), where

eTaŴØσtτ1aeTe+ØσTtτŴTŴØσtτ

we get

aDtαV1aλ1+KpeTe+1aeTe2+12Ŵ2Lϕ2eTe+1+KiaeTwγa12+12Ŵ2Lϕ2eTeE28

Here, we select 1+Kia=0and Kv=Ki+1, with Kv0then Ki1, with this selection of the parameters from Eq. (28) is reduced to:

aDtαV1aλ1+KpeTeγ1a12+12Ŵ2Lϕ2eTE29

From the previous inequality, we need to guarantee that Eq. (29) is less than zero, for which we select,

λ1+Kp>0,a>0, γ1>0, so that: aDtαV0,e,w,Ŵ0, e0, is wanted to be demonstrate.

The control law is given by Eq. (30)

un=Ω[ŴΓzxntτzxptτ
A+Ixxp+Kpe+KiaDtαe+KvaDtαeE30
Γ12+12Ŵ2Lϕ2en+frxrurAxr
ŴΓzxrtτxr+xp]

Theorem: The control law Eq. (30) and the neuronal adaptation law given by Eq. (25) guarantee that the fractional tracking error converges to zero, by which the tracking of trajectories of the non-linear system is guaranteed Eq. (5).

Corollary 2: If aDtαV1aλ1+KpenTenγ1a12+12Ŵ2Lϕ2enTen<0, e0,Ŵ, where Vis decreasing and bounded from below by V0, and:

V=12enTwTenwT+12atrWTW+1atτtØσTsŴTŴØσsds,

then we conclude that e,WL1; this means that the weights remain bounded.

5. Modeling of the time-delay adaptive neural network and the delayed plant

The nonlinear delayed unknown plant and the neural network are given as:

aDtαxp=fpxptτ+gpxptup,
aDtαxn=Ax+WΓzxtτ+wper+Ωu

where xp,fpRn,uRm,gpRnxn. And fp, is unknown and gp=I,A=λI, with ΓLypschitz function, Ware the fixed weigths but unknown from the neural networks, which minimize the modeling error.

Theorem: We will show that epand entend to zero and therefore etends to zero, that is, the neural network follows the plant.

For this proposal, we first define the modeling error between the neural network and the plant: ep=xpxn, whose derivative in the time is

aDtαep=aDtαxpaDtαxnE31

Adding and subtracting, to the right hand side from (34) the terms ŴΓzxrxtτ, αptŴ

aDtαep=AepWΓzxntτ+ŴØeptτ+upE32

6. Identification of the unknown non-linear system by the neural network

First, it is easy to see that eŴ=0is a state of equilibrium (equilibrium point). Previous, so we propose to demonstrate that this point of equilibrium is asymptotically stable; for this, be:

up=γ12+12Ŵ2Lϕ2epE33

We will show, the feedback system is asymptotically stable. Replacing (36) in (35)

aDtαep=AepWΓzxntτ+ŴØeptτγ12+12Ŵ2Lϕ2epE34

We will show, the new state epis asymptotically stable, and the equilibrium point is ep0, when Ŵσxntτ=0, as an external disturbance.

Let Vbe, the next candidate Lyapunov function as

V=12epTwTepwT+12atrWTWE35
+1atτtØσTsŴTŴØσsds

Then, (35) is reduced to

aDtαV1aλ1+KpepTepγ1a12+12Ŵ2Lϕ2epTep<0E36

The previous inequality guarantees that the identification of the non-linear system is satisfied, that is, the approach error converges to zero asymptotically

up=Ω[ŴΓzxrtτŴΓzxntτzxptτ
A+Ixxp+Kpe+KiaDtαe+KvaDtαeE37
Γ12+12Ŵ2Lϕ2enΓ12+12Ŵ2Lϕ2ep+frxrur
fpxp+AxpAxr+ŴΓzxpxr+xp]

7. Simulation

The mathematical model, which describes the movement dynamics of the bipedal robot, is obtained using the Euler-Lagrange equations [1, 13] ( Figure 2 ).

D(qtq¨t+Cqtq̇tq̇t+Gqt=t

Figure 2.

Dynamic model of biped robot.

where qt=q31tq32tq41tq42tT, is the generalized coordinates vector. As usual, D(qt) is the inertia matrix, bounded and positive definite, and Cqtq̇tis the matrix of Coriolis and centripetal forces. Gqtrepresents a matrix of gravitational effects and Bdefines the input matrix. The vector τt=τ31tτ32tτ41tτ42tT, defines the applied joint torques of the robot.

To illustrate the theoretical results obtained, we propose an example, which, as can be seen in the simulations, trajectory tracking is guaranteed.

The neural network is described by:

aDtαxp=Ax+WΓzxtτ+Ωu, with τ = 25 s, A=20I,IR4x4, and, Wis estimated using the learning law given in (28).

Γzxtτ=tanhx1tτtanhx2tτtanhxntτT, Ω=00100001Tand the uis obtained using (33).

and the reference signal that they have to follow, both the non-linear system and the neural network is given by the Duffing equation [14].

x¨x+x3=0.114cos1.1t:x0=1,ẋ0=0.114
xtdt=yt
ytdt=xtx3tαyt+δcosωt

Here, the conventional derivatives are replaced by the fractional derivatives as follows:

aDtαxt=yt
aDtαxt=xtx3tαyt+δcosωt

where α,ω,δ, are the parameters of the Duffing differential equation, which we will use as a reference trajectory, that the non-linear system and the neural network have to follow.

As can be seen in Figures 3 6 , the tracking of trajectories in the states of the system are performed with satisfaction, Figure 7 shows the phase plane of the Duffing equation, while Figure 8 shows the plane phase of the same fractional order differential equation.

Figure 3.

Time evolution for the angular position Leg 1 and Leg 2 (rad) of link 1.

Figure 4.

Time evolution for the angular position Leg 1 and Leg 2 (rad) of link 2.

Figure 5.

Time evolution for the angular position Leg 1 and Leg 2 (rad) of link 1.

Figure 6.

Time evolution for the angular position Leg 1 and Leg 2 (rad) of link 2.

Figure 7.

A phase space trajectory of Duffing equation.

Figure 8.

A phase space trajectory of Duffing equation.

Figures 9 12 show the torques applied to the ends of the bipedal robot.

Figure 9.

Torque (Nm) applied to Leg 1 and Leg 2 of link 1.

Figure 10.

Torque (Nm) applied to Leg 1 and Leg 2 of link 2.

Figure 11.

Torque (Nm) applied to Leg 1 and Leg 2 of link 1.

Figure 12.

Torque (Nm) applied to Leg and Leg 2 of link 2.

Parameter values of the fractional order, alpha (0.001) and beta (0.0001) are included.

α=1,β=1
α=0.001,β=0.001

8. Conclusions

In this chapter we study the mathematical model and control of non-linear systems, which are modeled by differential equations of fractional order, where it is observed that these systems have a better performance than the systems modeled by ordinary differential equations, those of fractional order they produce responses, solutions at simulation level, softer, by varying the order of the derivative.

The magnitude of the fractional order systems are smaller than the responses of the systems of ordinary differential equations, and with smaller control signals, which implies, less energy in the control process.

In this research work, conditions have been obtained in the parameters of the adaptive recurrent neural network, as well as laws of control and laws of neuronal adaptation, which, together, guarantee that the tracking error of trajectories between the non-linear system and the reference signal converges asymptotically to zero, so that trajectory tracking is develops with satisfaction.

Acknowledgments

Authors thank, Mexican National Science and Technology Council, (CONACyT), Mexico, to the Autonomous University of Nuevo Leon, Dynamical Systems Group of the Department of Physical and Mathematical Sciences, (FCFM-UANL), Centro de Investigación y Estudios de Posgrado, Facultad de Ingeniería, Dr. Manuel Nava #8, Zona Universitaria, C.P. 78290, San Luis Potosí, S.L.P., México The research of V. Ramírez-Rivera has been supported by CONACyT-México under the project Problemas Nacionales (2015-01-786).

Conflict of interest

The first author of the reference manuscript, in his name and that of all authors, declares that there is no potential conflict of interest related to the article.

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Joel Perez Padron, Jose P. Perez, C.F. Mendez-Barrios and V. Ramírez-Rivera (March 25th 2020). Trajectory Tracking Using Adaptive Fractional PID Control of Biped Robots with Time-Delay Feedback, Becoming Human with Humanoid - From Physical Interaction to Social Intelligence, Ahmad Hoirul Basori, Ali Leylavi Shoushtari and Andon Venelinov Topalov, IntechOpen, DOI: 10.5772/intechopen.90901. Available from:

chapter statistics

153total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Motion Generation during Vocalized Emotional Expressions and Evaluation in Android Robots

By Carlos T. Ishi

Related Book

First chapter

Design of an Assistive Gait Device for Strength Endurance and Rehabilitation

By K. H. Low, Xiaopeng Liu and Haoyong Yu

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us