Open access peer-reviewed chapter

# Trajectory Tracking Error Using Fractional Order PID Control Law for Two‐Link Robot Manipulator via Fractional Adaptive Neural Networks

By Joel Perez P., Jose Paz Perez P. and Martha S. Lopez de la Fuente

Submitted: December 14th 2016Reviewed: June 5th 2017Published: December 6th 2017

DOI: 10.5772/intechopen.70020

Downloaded: 436

## Abstract

The problem of trajectory tracking of unknown nonlinear systems of fractional order is solved using fractional order dynamical neural networks. For this purpose, we obtained control laws and laws of adaptive weights online, obtained using the Lyapunov stability analysis methodology of fractional order. Numerical simulations illustrate the obtained theoretical results.

### Keywords

• fractional order PID control
• fractional adaptive neural networks
• fractional Lyapunov functions
• fractional nonlinear systems
• trajectory tracking

## 1. Introduction

The fractional calculus is a branch of mathematics that attracted attention since G.W. Leibnitz proposed it in the seventeenth century. However, the researchers were not attracted to this area because of the lack of applications and analytical results of the fractional calculus.

On the contrary, the fractional calculus currently attracts the attention of a large number of scientists for their applications in different fields of science, engineering, chemistry, and so on.

This chapter presents the design of a fractional order nonlinear identifier modeled by a dynamic neural network of fractional order.

Although PID controllers are introduced long time ago, they are widely used in industry because of their advantages such as low price, design simplicity, and suitable performance. While three parameters of design including proportional (Kp), integral (Ki), and derivative (Kd) are available in PID controllers, two more parameters exist in FOPID controllers for adjustment. These parameters are integral fractional order and derivative fractional order. In comparison with PID controllers, FOPID controllers have more flexible design that results in more precise adjustment of closed‐loop system. FOPID controllers are defined by FO differential equations. It is possible to tune frequency response of the control system by expanding integral and derivative terms of the PID controller to fractional order case. This characteristic feature results in a more robust design of control system, but it is not easily possible. According to nonlinearity, uncertainty, and confusion behaviors of robot arms, they are highly recommended for experimenting designs of control systems. Despite nonlinear behavior of robot arm, it is demonstrable that a linear proportional derivative controller can stabilize the system using Lyapunov. But, classic PD controller itself cannot control robot to reach suitable condition. Several papers and wide researches in optimizing performance of the robot manipulator show the importance of this issue.

There are several ways of defining the derivative and fractional integral, for example, the derivative of Grunwald‐Letnikov given by Eq. (1)

aDtαf(t)=limh01hαj=0[(tα)h](1)j(αj)f(tjh)E1

where [.] is a flooring operator, while the RL definition is given by:

aDtαf(t)=1Γ(nα)dndtnatf(τ)(tτ)αn+1dτE2

For (n1<α<n), Γ(x)is the well‐known Euler’s Gamma function.

Similarly, the notation used in ordinary differential equations, we will use the following notation, Eq. (3), when we are referring to the fractional order differential equations where αkR+.

which is:

g(t, x,aDtα1 x, aDtα2 x, )=0E3

The Caputo’s definition can be written as

aDtαf(t)=1Γ(αn)atf(n)(τ)(tτ)αn+1dτE4

For

(n1<α<n).E34

Trajectory tracking, synchronization, and control of linear and nonlinear systems are a very important problem in science and control engineering. In this chapter, we will extend these concepts to force the nonlinear system (plant) to follow any linear and nonlinear reference signals generated by fractional order differential equations.

The proposed adaptive control scheme is composed of a recurrent neural identifier and a controller (Figure 1). Figure 1.Recurrent neural network scheme.

We use the above scheme to model the unknown nonlinear system by means of a dynamic recurrent neural network of adaptable weights; the above is modeled by differential equations of fractional order. Also, the scheme allows us to determine the control actions, the error of approach of trajectories, as well as the laws of adaptation of adaptive weights and the interconnection of such systems.

## 2. Modelling of the plant

The nonlinear system (Eq. (5)) is forced to follow a reference signal:

aDtαxp=Fp(xp,u)fp(xp)+ gp(xp)upxp, fp.Rn, uRm, gpRnxn.E5

The differential equation will be modeled by:

aDtαxp=A(x)+W*Γz(x)+Ωu.E36

The tracking error between these two systems:

wper=xxpE6

We use the next hypotheses.

aDtαwper=kwperE7

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

aDtαxp= aDtαx+wperE37

The nonlinear system is :

aDtαxp=aDtαx+wper=A(x)+W*Γz(x)+wper+ΩuE8

Where the W*is the matrix weights.

## 3. Tracking error problem

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

aDtαxr=fr(xr,ur),  wr, xrRnE9

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

To achieve our goal of trajectory tracking, we propose

e=xpxrE10

The time derivative of the error is:

aDtαe=aDtαxpaDtαxr=A(x)+W*Γz(x)+wper+Ωufr(xr,ur)E11

The Eq. (11) can be rewritten as follows, adding and subtracting the next terms W^Γz(xr), αr(t,W^), Aeand wper=xxp,; then,

aDtαe=A(x)+W*Γz(x)+xxp+Ωufr(xr,ur)+W^Γz(xr)W^Γz(xr)+Ωαr(t,W^)Ωαr(t,W^)+AeAeaDtαe=Ae+W*Γz(x)+Ωufr(xr,ur)+W^Γz(xr)+Ωαr(t,W^)W^Γz(xr)Ωαr(t,W^)Axrxr+x+A(x)E12

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

Axr+W^Γz(xr)+xrxp+Ωαr(t,W^)=fr(xr,ur), where

Ωαr(t,W^)=fr(xr,ur)AxrW^Γz(xr)xr+xpE13
aDtαe=Ae+W*Γz(x)W^Γz(xr)Ae+(A+I)(xxr)+Ω(uαr(t,W^))E14

Now, W^is part of the approach, given by W*. The Eq. (14) can be expressed as Eq. (15), adding and subtracting the term W^Γz(x)and if Γz(x)=Γ(z(x)z(xr))

aDtαe=Ae+(W*W^)Γz(x)+W^Γ(z(x)z(xr))+(A+I)(xxr)Ae+Ω(uαr(t,W^))E15

If

W=W*W^andu=uαr(t,W^)E16

And by replacing Eq. (16) in Eq. (15), we have:

aDtαe=Ae+WΓz(x)+W^Γ(z(x)z(xr))+(A+I)(xxr)Ae+Ωu
aDtαe=Ae+WΓz(x)+W^Γ(z(x)z(xp)+z(xp)z(xr))+(A+I)(xxp+xpxr)Ae+ΩuE17

And:

u=u1+u2E18

So, the result for Ωu1is

Ωu1=W^Γ(z(x)z(xp))(A+I)(xxp),E19

and Eq. (17) is simplified:

aDtαe=Ae+WΓz(x)+W^Γ(z(xp)z(xr))+(A+I)(xpxr)Ae+ΩuE38

Taking into account that e=xpxr, the equation for aDtαeis

aDtαe=(A+I)e+W˜Γz(x)+W^Γ(z(e+xr)z(xr))+Ωu2=(A+I)e+Wσ(x)+W^(σ(e+xr)σ(xr))+Ωu2E39

If φ(e)=(σ(e+xr)σ(xr)), then

aDtαe=(A+I)e+Wσ(x)+W^φ(e)+Ωu2E20

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

## 4. Asymptotic stability of the approximation error

From Eq. (20), we consider the stability of the tracking error, for which we first observe that (e,W^)=0, is an equilibrium state of dynamical system from Eq. (20), and we consider a particular case when A=λI,λ>0

For such stability analysis of the trajectory tracking (Eq. (20)), we propose the following FOPID control law :

Ωu2=Kpe+KiaDtαe+KvaDtαeγ(12+12W^2Lφ2)e E21

Our objective is to find Kp, Ki, Kv,W^, Lφz2, and this guarantees that the tracking error given by Eq. (20) is asymptotically stable, for which we will later propose a Lyapunov function, with γ>0; this control law (Eq. (21)) is similar to .

A FOPID controller, also known as a PIλDαcontroller, takes on the form :

u(t)=Kpe(t)+KiaDtλe(t)+KdaDtαe(t)E40

where λand αare the fractional orders of the controller and e(t)is the system error, where λ=α. Note that the system error e(t)replaces the general function f(t).

We will show that the feedback system is asymptotically stable. Replacing Eq. (21) in Eq. (20), we have

aDtαe=(A+I)e+Wσ(x)+W^φ(e)+Kpe+KiaDtαe+KvaDtαeγ(12+12W^2Lφ2)e, then

(1Kv)aDtαe=(A+I)e+Wσ(x)+W^φ(e)+Kpe+KiaDtαeγ(12+12W^2Lφ2)e. If

a=(1Kv), then

aDtαe=1a(A+I)e+1aWσ(x)+1aW^φ(e)+1aKpe+1aKiaDtαeγa(12+12W^2Lφ2)eE22
aDtαe=1a(λ1+Kp)e+1aWσ(x)+1aW^φ(e)+1aKiaDtαeγa(12+12W^2Lφ2)eE23

And if w=1aKiaDtαe, then aDtαw=1aKie(t); then, we rewrite Eq. (23) as:

aDtαe=1a(λ1+Kp)e+1aWσ(x)+1aW^φ(e)+wγa(12+12W^2Lφ2)eE24

We will show that the new state (e,w)Tis asymptotically stable and the equilibrium point is (e,w)T=(0,0)T, when Wσ(xr)=0, as an external disturbance.

Let Vbe, the next candidate Lyapunov function as :

V=12(eT,wT)(e,w)T+12atr{WTW}E25

The fractional order time derivative of Eq. (25) along with the trajectories of Eq. (24) is

aDtαV=eTaDtαe+wTaDtαw+1atr{aDtα WTW}E26
aDtαV=eT(1a(λ1+Kp)e+1aWσ(x)+1aW^φ(e)+wγa(12+12W^2Lφ2)e)+1aWTKie+1atr{aDtα WTW}E27

In this part, we select the next learning law from the neural network weights as in  and :

tr{aDtα WTW}=eTWσ(x)E28

Then, Eq. (27) is reduced to

aDtαV=1a(λ1+Kp)eTe+eTaW^φ(e)+(1+Kia)eTwγa(12+12W^2Lφ2)eTeE29

Next, lets consider the following inequality proved in 

XTY+YTXXTΛX+YTΛ1YE30

which holds for all matrices X,YRnxkand ΛRnxnwith Λ=ΛT>0. Applying (30) with Λ=Ito the term eTaW^φ(e)from Eq. (29), where

eTW^φ(e)12e2+12Lφ2W^2e2=12(1+Lφ2W^2)e2E43

we get

aDtαV1a(λ1+Kp)eTe+1a(eTe2+12W^2Lφ2)eTe+(1+Kia)eTwγa(12+12W^2Lφ2)eTeE31

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

aDtαV1a(λ1+Kp)eTe(γ1)a(12+12W^2Lφ2)eTeE32

Of the previous inequality, Eq. (32), we need that the fractional order Lyapunov derivative, aDtαV0, to ensure that the trajectory tracking error is asymptotically stable, that is, limte(t)=0, which means that the nonlinear system follows the reference signal.

To achieve this purpose, we select:

λ1+Kp>0,a>0,(γ1)>0,aDtαV0, e, w, W^0,e0E9000

With the above Eq. (32), the control law that guarantees asymptotic stability of the tracking error is given by Eq. (33)

u=Ω[W^Γ(z(x)z(xp))(A+I)(xxp)+Kpe+KiaDtαe+KvaDtαeγ(12+12W^2Lφ2)e+fr(xr,ur)AxrW^Γz(xr)xr+xpE33

Theorem: The control laws (Eq. (33)) and the adaptive weights (Eq. (28)) ensure that the trajectory tracking error between the fractional nonlinear system (Eq. (8)) and the fractional reference signal (Eq. (9)) satisfies limte(t)=0

Remark 2: From Eq. (32), we have

aDtαV1a(λ1+Kp)eTe(γ1)a(12+12W^2Lφ2)eTe<0,  e0,  W^, where Vis decreasing and bounded from below by V(0), and:

V=12(eT,wT)(e,w)T+12atr{WTW},; then we conclude that e, WL1; this means that the weights remain bounded.

## 5. Simulation

The manipulator used for simulation is a two revolute joined robot (planar elbow manipulator), as shown in Figure 2. Figure 2.Diagram of the prototype planar robot with two degrees of freedom.

The dynamics of the robot is established by [10, 11], Mij(q),  i, j=1,2of the inertia matrix M(q)as

M11(q)=m1lc12+m2(l12+lc22+2l1lc2cos(q2))+I1+I2;E45
M12(q)=m2(lc22+l1lc2cos(q2))+I2;E46
M21(q)=m2(lc22+l1lc2cos(q2))+I2;E47
M22(q)=m2lc22+I2.E48
C11(q,q˙)=m2l1lc2sin(q2)q˙2;E49
C21(q,q˙)=m2l1lc2sin(q2)(q˙1+q˙2);E50
C21(q,q˙)=m2l1lc2sin(q2)q˙1E51
C22(q,q˙)=0E52

And the torque vector g(q):

g1(q)=(m1lc1+m2l1)gsin(q1)+m2lc2gsin(q1+q2);E53
g2(q)=m2lc2gsin(q1+q2);E54

Thus, it is possible to write the equations of motion using the Lagrange equations for fractional manipulator system as :

(m1+m2)l12θ¨1+m2l1l2θ¨2cos(θ2θ1)m2l1l2θ˙22sin(θ2θ1)+(m1+m2)gl1sin(θ1)+(α1)(tτ)[(m1+m2)l12θ˙1+m2l1l2θ˙2cos(θ2θ1)]=Q1E55
m2l22θ¨2+m2l1l2θ¨1cos(θ2θ1)+m2l1l2θ˙12sin(θ2θ1)+m2gl2sin(θ2)+(α1)(tτ)[m2l22θ˙2+m2l1l2θ˙1cos(θ2θ1)]=Q2E56

The terms containing αindicate the additional terms resulting from the fractional order model and the right‐hand sides denote the generalized force terms resulting from the forcing functions, and there is a specific set of values for Q1and Q2for each case.

With the end of supporting the effectiveness of the proposed controller, we have used a Duffing equation.

The fractional order neural network is modelling by the differential equation:

aDtαxp=A(x)+W*Γz(x)+Ωu, with A=λI,IR4x4, and λ=20, W*is estimated using the learning law given in Eq. (28).

Γz(x)=(tanh(x1), tanh(x2),,tanh(xn))T, Ω=(00   1000  01)Tand the uis calculated using Eq. (33). The plant is stated in  and , and it is given by:

D(q)q¨+C(q,q˙)q˙+G(q)=τE57

We try to force this manipulator to track a reference signal  given by undamped Duffing equation:

x¨x+x3=0.114cos(1.1t):x(0)=1,x˙(0)=0.114E58

To get the fractional order Duffing’s system, this equation can be rewritten as a system of the first‐order autonomous differential equations in the form :

x(t)dt=y(t)E59
y(t)dt=x(t)x3(t)αy(t)+δcos(ωt)E60

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

aDtαx(t)=y(t)E61
aDtαx(t)=x(t)x3(t)αy(t)+δcos(ωt)E62

where αis the fractional orders and α, δ, ωare the system parameters.

Illustrated, the response in the time, angular position and torque applied to the fractional nonlinear system are shown in Figures 37. As can be observed, the trajectory tracking objective is obtained Figure 3.A phase space trajectory of Duffing equation. Figure 4.Time evolution for the angular position q 1 (rad) of link 1. Figure 5.Time evolution for the angular position q 2 (rad) of link 2. Figure 6.Torque (Nm) applied to link 1. Figure 7.Torque (Nm) applied to link 2.
α=1,  β=1E63

Its phase space trajectory is given in Figure 8, and the time evolution for the position angles and applied torque are shown in Figures 912. As can be seen in Figures 9 and 10, the trajectory tracking is successfully obtained where plant and reference signals are the same. Figure 8.A phase space trajectory of Duffing equation. Figure 9.Time evolution for the angular position q 1 (rad) of link 1. Figure 10.Time evolution for the angular position q 2 (rad) of link 2. Figure 11.Torque (Nm) applied to link 1. Figure 12.Torque (Nm) applied to link 2.
α=0.99,  β=0.99E64

Its phase space trajectory is given in Figure 13, and the time evolution for the position angles and applied torque are shown in Figures 1417. As can be seen in Figures 14 and 15, the trajectory tracking is successfully obtained where plant and reference signals are the same. Figure 13.A phase space trajectory of Duffing equation. Figure 14.Time evolution for the angular position q 1 (rad) of link 1. Figure 15.Time evolution for the angular position q 2 (rad) of link 2. Figure 16.Torque (Nm) applied to link 1. Figure 17.Torque (Nm) applied to link 2.
α=0.001,  β=0.001E65

As can be observed, in the graphs of the trajectory tracking, the experimental results obtained in this chapter show a good experimental performance. The laws of control are obtained online, as well as the laws of adaptive weights in the fractional order neural network.

The control laws obtained are robust to modeling errors and nonmodeled dynamics (unknown nonlinear systems).

## 6. Conclusions

We have discussed the application of the stability analysis by Lyapunov of fractional order to follow trajectories of nonlinear systems whose mathematical model is unknown. The convergence of the tracking error is established by means of a Lyapunov function, as well as a control law based on Lyapunov and laws of adaptive weights of fractional order dynamical neural networks.

The results show a satisfactory performance of the fractional order dynamical neural network with online learning.

## Acknowledgments

Authors thank Mexican National Science and Technology Council, (CONACyT), Mexico.

The first and second authors thank the Autonomous University of Nuevo Leon, Dynamical Systems Group of the Department of Physical and Mathematical Sciences (FCFM‐UANL), Mexico.

The third author thanks the University of Monterrey (UDEM), Mexico.

## More

© 2017 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

### Cite this chapter Copy to clipboard

Joel Perez P., Jose Paz Perez P. and Martha S. Lopez de la Fuente (December 6th 2017). Trajectory Tracking Error Using Fractional Order PID Control Law for Two‐Link Robot Manipulator via Fractional Adaptive Neural Networks, Robotics - Legal, Ethical and Socioeconomic Impacts, George Dekoulis, IntechOpen, DOI: 10.5772/intechopen.70020. Available from:

### chapter statistics

436total chapter downloads

### More statistics for editors and authors

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

### Related Content

Next chapter

#### Robots Liability: A Use Case and a Potential Solution

By Alejandro Zornoza, José C. Moreno, José L. Guzmán, Francisco Rodríguez and Julián Sánchez-Hermosilla

First chapter

#### Introductory Chapter: Drones

By George Dekoulis

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.

View all books