SlideShare a Scribd company logo
International Journal of Electrical and Computer Engineering (IJECE)
Vol. 9, No. 2, April 2019, pp. 1078~1089
ISSN: 2088-8708, DOI: 10.11591/ijece.v9i2.pp1078-1089  1078
Journal homepage: http://iaescore.com/journals/index.php/IJECE
Using real interpolation method for adaptive identification of
nonlinear inverted pendulum system
Phu Tran Tin1
, Tran Hoang Quang Minh2
, Tran Thanh Trang3
, Nguyen Quang Dung4
1
Faculty of Electronics Technology, Industrial University of Ho Chi Minh City, Ho Chi Minh City, Vietnam
2
Optoelectronics Research Group, Faculty of Electrical and Electronics Engineering,
Ton Duc Thang University, Ho Chi Minh City, Vietnam
3
Faculty of Electrical and Electronics Engineering, Ho Chi Minh City University of Food Industry,
140 Le Trong Tan, Ho Chi Minh City, Vietnam
4
Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam
Article Info ABSTRACT
Article history:
Received Jun 3, 2018
Revised Oct 10, 2018
Accepted Dec 8, 2018
In this paper, we investigate the inverted pendulum system by using real
interpolation method (RIM) algorithm. In the first stage, the mathematical
model of the inverted pendulum system and the RIM algorithm are
presented. After that, the identification of the inverted pendulum system by
using the RIM algorithm is proposed. Finally, the comparison of the linear
analytical model, RIM model, and nonlinear model is carried out. From the
results, it is found that the inverted pendulum system by using RIM
algorithm has simplicity, low computer source requirement, high accuracy
and adaptiveness in the advantages.
Keywords:
Identification
Inverted pendulum
Non-linear system
Real interpolation method
RIM Copyright © 2019 Institute of Advanced Engineering and Science. All rights
reserved.
Corresponding Author:
Tran Hoang Quang Minh,
Optoelectronics Research Group,
Faculty of Electrical and Electronics Engineering,
Ton Duc Thang University,
19 Nguyen Huu Tho, District 7, Ho Chi Minh City, Vietnam.
Email: tranhoangquangminh@tdtu.edu.vn
1. INTRODUCTION
Inverted pendulum system, which is an underactuated mechanical system and exhibits strong
nonlinearity, has been considered as an exciting and challenging problem in automation. Notably, their
features such as swing up, stabilization problem, higher order, nonlinear, robust coupling, and multivariate
have been attracted more and more attention from control researchers for the last decades. Moreover,
inverted pendulum system is widely used in the field of robotics and aerospace field, and so has essential
significance both in the field of the theoretical study research and practice [1]-[3]. There are a lot of
researchers about most popular types inverted pendulums such as the rotational single arm pendulum, the cart
inverted pendulum, and the double inverted investigated the inverted pendulum by switching two different
control laws [3], and in [4] LQR, double-PID and pole placement control techniques are used to control a cart
inverted pendulum system. Moreover, fractional order PID controllers are proposed in [5] and [6] for control
of the inverted pendulum system. In [7], a simulation study on PID and LQR control of a cart inverted
pendulum system is carried out with and without disturbance input. Furthermore, in [8], control of a cart
inverted pendulum on an inclined surface is performed using fuzzy and PID controllers. An ANFIS controller
is proposed in [9] as an intelligent control technique to control an inverted pendulum. In addition, in [10],
sliding mode controller is proposed to control a cart inverted pendulum system in both simulation and real-
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1079
time implementation. From this point of view, research in the field of inverted pendulum system is still a hot
direction in both theory and practice automation area.
In this paper, we investigate the inverted pendulum system by using real interpolation method
(RIM) algorithm. In the first stage, the mathematical model of the inverted pendulum system and the real
interpolation method (RIM) algorithm are presented. After that, the identification of the inverted pendulum
system by using the real interpolation method (RIM) algorithm is proposed. Investigated and derived.
Finally, the analytical algorithm derived by Simulink Mat Lab Software. From the results, it is found that the
inverted pendulum system by using real interpolation method (RIM) algorithm is simplicity, low computer
source requirement, high accuracy. The main contributions of this paper are presented as the following:
1) The mathematical model of the cart inverted pendulum system is proposed.
2) The identification of the inverted pendulum system by using the real interpolation method (RIM)
algorithm is presented and derived.
The rest of the paper is proposed as follows. The cart inverted pendulum system, the real
interpolation method (RIM) and the identification of the inverted pendulum system by using the real
interpolation method (RIM) algorithm is introduced in Section 2. After that, some results and discussions are
given in Section 3. Finally, conclusions are provided in Section 4.
2. RESEARCH METHOD
2.1. Mathematical model of the cart inverted pendulum
The cart inverted pendulum model is presented in Figure 1 [11].
Figure 1. The cart inverted pendulum model
The equations of pendulum motion (a nonlinear system) can be formulated as:
     2
cos sin 1M m ml x mlx b x F        
, (1)
     2
cos 2 sinI ml mlx b mgl        . (2)
Furthermore, we try to linearize the system differential equations. It is a linearization around the
equilibrium point, in our case about vertical position . The prerequisite is that the system remains in the
vicinity of this position. Let presents a deviation from the equilibrium point, we get . When
pendulthe um is in upright position     2 2
cos 1, sin , 0cos sin                  .
Using above relation and from the nonlinear Equations (1) and (2) we obtained linearized equations
of the proposed system as in (3) and (4).
 2
2I ml b mgl mlx     
, (3)
  1x M m b x ml U    . (4)
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089
1080
Note that U has been substituted for the input force F. To obtain the transfer functions of the
linearized system equations, we must first take the Laplace transform of the system equation assuming zero
initial conditions. The resulting Laplace transforms are shown below:
         2 2 2
2Φ Φ ΦI ml s s b s s mgl s mlX s s   
, ( (5)
         2 2
1 ΦM m X s s b X s s ml s s U s    . (6)
In this analysis, we set
2
( )q I M m Mml   . (7)
The transfer function of pendulum’s position can determine from systems of Equations (5), (6) and
(7) as the following:
 
 
 
     
2
2 2
2
1 2 14 3 2
X
2 1
I ml b mgl
s s
s q q q
U s M m b I ml b mgl M m b b mglb
s s s s
q q q

 

          
   (8)
Here transfer function of pendulum’s angle can be formulated as:
 
       2
13 2 1
Φ
2 2 1
ml
s
s q
U s M m b I ml b mgl M m b b mglb
s s s
q q q

        
(9)
We consider the equilibrium point , the pendulum become the normal form. The transfer
functions of the pendulum’s angle and position are presented as in the formula (10) and (11).
 
 
 
     
2
2 2
2
1 2 14 3 2
X
2 1
I ml b mgl
s s
s q q q
U s M m b I ml b mgl M m b b mglb
s s s s
q q q

 

          
   , (10)
 
       2
13 2 1
Φ
2 2 1
ml
s
s q
U s M m b I ml b mgl M m b b mglb
s s s
q q q


        
. (11)
For the transfer function above are valid these units, we assume that:
 
 
 
 
 
Φ
,k v
s X srad m
P P s
U s N U s N
   
    
   
, or
 
 
 
 
 
Φ 180
,k v
s X sdeg m
P P s
U s N U s N
   
     
   
.
According to the transfer functions (8) and (9), the inverted pendulum with the equilibrium point
is unstable system. Due to unstable behavior of the upright equilibrium position, identification of the
inverted pendulum represents a difficult task. However, the normal pendulum with the equilibrium point
is a stable system with respect to the transfer functions (10) and (11). Looking at the details, it is
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1081
clearly seem that the absolute values of the coefficients of the transfer functions are the same for the
pendulum’s upright and downright positions. The difference is only in the sign of the coefficients. Based on
this conclusion, the identification task of the nonlinear pendulum system can be applied on direct pendulum
characteristics. In this paper, the response characteristics of the pendulum with the downright equilibrium
point, will be consider as input data for the identification task in the next subsection [11]-[12].
2.2. Real interpolation method (RIM)
The real interpolation method is one of the methods, which works on mathematical descriptions of
the imaginary area. The method is based on real integral transform as follow [13]:
     
0
, , , 0t
F f t e dt C C
 


   
(12)
which assigns to the original the image as a function of the real variable . Formula of direct
transform (12) can be considered as a special case of the direct Laplace transform by replacing the complex
variable for real variable. Another step towards the development of the instrumentation
method is the transition from continuous functions to their discrete form, using the computing resources
and numerical methods. For these purposes, real interpolation method is represented by numerical
characteristics { } . They are obtained as a set of values of the function in the nodes ̅̅̅̅̅,
where is the number of elements numerical characteristics, called its dimension.
Selecting of interpolations is a primary step in the transition to a discrete form, which has a
significant impact on the numerical computing and accuracy of problem solutions. Distribution of nodes in
the simplest variant is uniform. Another important advantage of real interpolation method is cross-conversion
properties [11]. It dues to the fact that the behavior of the function for large values of the argument is
determined mainly by the behavior of the original for small values of the variable . In the opposite
case, the result is the same: the behavior of the function for small values of the argument is
determined mainly by the behavior of the original for large values of the variable [13].
When considering the original of dynamic characteristics of dynamic systems, formula (12)
leads to an operator model, which under certain conditions can be considered as special cases of the models
based on the Laplace transform. Thus, in (12) replacing of the function by - the impulse transient
function of the dynamic system, we obtain its transfer function. From here we can find the elements of a
discrete model of the system, and its transfer function by performing the discretization procedures for nodes
̅̅̅̅̅
   
0
, 1,it
iW k t e dt i
 


 
(13)
Function is a real transfer function of control automatic systems, having an impulse transient
response . Function could be received based on determination of transfer function such as a
relationship of the imaginary of output and input signals
 
 
 
Y
W
X




(14)
In which the imaginary of the output signal and the input signal is calculated from the original functions of
the input and output signals:
   
0
t
Y y t e dt




, (15)
   
0
t
X x t e dt




, (16)
where input signal, and output signal of the system.
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089
1082
The input-output relationship of the system has a form like below:
     Y W X  
. (17)
2.3. Identification of the inverted pendulum based on the RIM
The inverted pendulum is an unstable system, investigation of this system based on input-output
signal use two kinds of the input signal, namely impulse and steps responses. For unstable systems,
identification of these systems using step response is not entirely suitable, because of the unstable transient
response under the effect of the step response. Using impulse response can be used, however, in practice, the
idea impulse response has challenging to implement. In this paper, the conversion impulse response is used
for identification of nonlinear system. The inverted pendulum is used conversion impulse response having a
form of transfer function:
 
 1 exp *s
X s M
s
 

, (18)
where is magnitude, delay time.
Under the effect of a conversion impulse response (18) on the output of the system, response
function is represented as a function . From formula (14) and (18), in the real form, the transfer function
of the inverted pendulum has a form.
 
 
 
1 exp *
M
W Y

 
 

 
, (19)
where    
0
t
Y k t e dt



 if response characteristics represents as a function , or    
1
i
N
t
i
i
Y k t e t
 

 
, is a sample of ..
The linear transfer function of the angle of the inverted pendulum can be formulated as:
  1
3 2
3 2 1
Φ  LR
B s
s
s A s A s A


  
(20)
The linear transfer function of the position of the inverted pendulum can be calculated as:
 
2
2 1 0
4 3 2
3 2 1
LR
B s B s B
s
s s s s
X
A A A
 

  
(21)
The linear transfer function of angle and position of the inverted pendulum has a form of real interpolation
transformation can be formulated as (22) and (23):
  1Φ
3 2
3 2 1
ΦLR
B
A A A


  


  
, (22)
 
2
2 1 0
4 3 2
3 2 1
X X X
LRX
B B B
A A A
 

   
 

  
(23)
Interestingly, the coefficients of the transfer function of the angle and position of the pendulum have
the same number in the denominator. The identification process can be determined the coefficients of the
denominator and separately the coefficients of the numerator. Firstly, identification of the transfer function of
pendulum angle. Chosen structure and order of the linear identification like formula (20) and (21) or (22) and
(23) in the real form. The number of the unknown coefficients for the angle of the pendulum is four; it means
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1083
that the dimension of the elements numerical characteristics is four. In a simple case, the value of nodes
given be calculated by nominal distribution in the range of , where
 Φ 1
1 Φ
Φ
1 , 2,
1
i i i
 
  


   

(24)
Discretization of by real interpolation method in ̅̅̅̅̅̅
 
  1Φ
Φ3 2
3 2 1
1 , 1,
1 exp *
i i
i
i i i i
M B
Y i
A A A
 
 
    

 
    
, (25)
where ∫ , for the case, the transient characteristics is given a form of function with
time. For the case, the transient characteristics is given a table or massive of values with sample time .
   Φ
1
1 Δi i
N
t
i i
i
Y k t e t
 

 
(26)
where, –sampling time, N-integrator. Here, we set replacement as:
 
 
  Φ1 1 , 1,
1 exp *
i
i i
i
M
F Y i

  
 
 
 
(27)
From (14) and (16) we have a system of the equation like:
  1Φ
Φ3 2
3 2 1
1 , 1,i
i
i i i
B
F i
A A A

 
  

 
  
(28)
or
  3 2
3 2 1 1Φ Φ1 , 1,i i i i iF A A A B i          
(29)
     2 3
3 2 1 1Φ Φ1 1 1 , 1,i i i i i i iF A F A A B F i             (30)
Rewrite the system of equations of the unknown coefficients
Matrix M1:
   
   
   
   
2
1 1 1 1 1
2
2 2 2 2 2
2
3 3 3 3 3
2
4 4 4 4 4
1 1 1
1 1 1
1
1 1 1
1 1 1
F F
F F
M
F F
F F
    
    
    
    
 
 
 
 
 
   (31)
Matrix N1
 
 
 
 
3
1 1
3
2 2
3
3 3
3
4 4
1
1
1
1
1
F
F
N
F
F
 
 
 
 
 
 
 
 
 
   (32)
Matrix X1 of the unknown coefficients
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089
1084
3
2
1
1Φ
1
A
A
X
A
B
 
 
 
 
 
  (33)
If the matrixes M1, N1 are known, the matrix X1 can be determined as the following:
1
1 1 1X M N
 (34)
The identification task leads to solving a system of equation of unknown coefficients of tthe
ransfer function. Identification of transfer function of pendulum’s position. We already determined the
coefficient of denominator n above step, in terms of simplicity, requirement determines unknown coefficients
of the numerator of the transfer function of pendulum’s position. The number of characteristics is .
Transfer function of pendulum’s position has the form as (23):
 
2
2 1 0
4 3 2
3 2 1
X X X
LRX
B B B
A A A
 

   
 

  
(35)
Obviously, in order to reduce identification task, the coefficient has already
determined in the first section, the identification task remains only determine coefficients of the numerator of
the transfer function of pendulum’s position . The identification procedure of transfer function
of pendulum’s position illustrates following steps. If response characteristics are given a massive of time,
with the sample time , value of pendulum position at and N-number of samples, numerical
characteristics are calculated by following form:
   
1
2 Δ  i i
N
t
i X i
i
Y k t e t
 

 
(36)
3. RESULTS AND DISCUSSION
The inverted pendulum has parameter as in Table 1.
Table 1. Parameters of the Pendulum System
Parameter Value Unit
Mass of the cart (M) 0.5 Kg
Mass of the pendulum (m) 0.2 Kg
Length of the pendulum (l) 0.3 M
Coefficient of friction for the cart (b1) 0.2 Nm-1s-1
Coefficient of friction for the pendulum (b2) 0.002 Nrad-1s-1
Mass moment of inertia of the pendulum (I) 0.006 kg/m2
Gravitation force (g) 9.81 m/s2
Force applied to the cart (F) - N
Position of the cart (x) - m
Pendulum angle (φ) - deg.
In this example we will compare results of the nonlinear modeling, describing as the formulas (1)
and (2), with a linear model, solving by the analytical approximation as formula (10) and (11) and with a
RIM model, which is identified by RIM. The pendulum has a set of parameters as the Table 1. On input of
the models, the rectangular input acceleration signal has 1 N in time 1 s in the Figure 2(a).
In the Figure 3, obviously, the RIM model and analytical linear model highly fit with nonlinear
model in the beginning of the transient process. However, in the end of the transient process, the linear model
represents higher in fitness comparing with the RIM model. Moreover, Figure 4 shows that the RIM model
has more accuracy at the beginning of the span time [0-7 s]. However, the fitness of RIM model is reducing
with the increase of time in the range 7-70 s. From the results, the maximum error of linear is 0.24
approximation by RIM is 0.45. Generally, the linear model represents a higher accuracy model than RIM
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1085
model when the angle deviation from the downright equilibrium position is lower than 15o
. The position of
the pendulum versus time is presented in Figure 5. Interestingly, the RIM model of position of pendulum
illustrates higher fitness level of nonlinear model than the analytical linear model. According to the absolute
errors of the analytical linear model and the RIM model presented in the Figure 6. The error of the RIM
model is considerably lower than of the analytical linear model in the time range [0-10 s]. In the rest of the
considered settling time, the analytical model is around 1.5 times as high as the RIM model. The maximum
error of the analytical and RIM models account for 0.038 and 0.026, respectively
(a) (b)
Figure 2. The input acceleration signal
Figure 3. Response characteristics of nonlinear – green-dot line, analytic - blue-dot-line, and identification by
RIM –red line
The input acceleration force with 3 N, illustrated in Figure 2(b) is active on a cart during 1 second,
based on the samplings of the nonlinear model, we find the approximation model by RIM for the new set of
data. For identification of the deviation angle of pendulum, the nodes considered in the range
, sample time and settling time 70 s. Then making a comparison with the
original nonlinear model, linear model, and RIM model, we show the results as Figure 7 and Figure 8.
Overall, the Figure 7 show that the RIM model represents more accuracy than the analytical linear model,
especially in the range 0-10s of the settling time. Interestingly, with a high displacement of the angle from the
equilibrium point , the RIM model represents more accuracy then linear model. Maximum error of the
linear model and the RIM model is 9.17 and 8.62, respectively.
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089
1086
Figure 4. Approximation error: Linearization error-red line, approximation error by RIM - blue line
Figure 5. The position of the pendulum: Nonlinear model red-line, linearization model- blue line, approximation
model by RIM – green line
Figure 6. The error of models: linearization model – red line, approximation model by RIM – blue line
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1087
Figure 7. The angle of the pendulum: nonlinear model – blue dot line, linear model – violet dot line,
approximation by RIM – red line
Figure 8. The error of angle – linear model – red line, RIM model – blue line
Position of the pendulum under acting by 3 N during 1 second presents in Figure 9. Obviously, the
RIM model has a higher fitness than the analytical linear model in the beginning of the settling time [0-10 s].
In the end of the settling time the fitness of both models is approximately the same. Looking at more details
in the Figure 10, the fitness of the RIM model is higher than of the analytical model in time range [0-20 s].
Figure 9. The position of the pendulum: nonlinear model –red line, linear model – blue line, RIM model –
green line
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089
1088
Figure 10. The error of model- linearization model –red line; approximation model by RIM – blue line
However, in the further transient time, the RIM model is slightly less accuracy than the analytical
linear model. The maximum error of the deviation of the pendulum’s position of both models is the same at
about 0.83. According to the two careful samples, the models identified by RIM represent higher accuracy
than model of analytical linear approach in cases of that the deviation of pendulum accounts for higher 45o
.
RIM provides adaptively tool to obtain the linear transfer function of non-linear system based on transient
characteristics, which has significant advantage compare analytical methods. The method which is based on
expression from physical laws of dynamic pendulum. Another advantage of RIM is that the method is
conducted in real domain, which is more effective compared to other methods implement in imaginary and
frequency domains.
4. CONCLUSION
In this paper, the inverted pendulum system by using real interpolation method (RIM) algorithm is
investigated. After the identification of the inverted pendulum system by using the real interpolation method
(RIM) algorithm is analyzed, the analytical algorithm derived by Simulink Mat Lab Software. From the
results, it is found that the inverted pendulum system by using real interpolation method (RIM) algorithm is
simplicity, low computer source requirement, high accuracy. It can be found like the prospective solution for
the inverted pendulum research direction in the near future.
REFERENCES
[1] Glück, Tobias, Andreas Eder, and Andreas Kugi, “Swing-up Control of a Triple Pendulum on a Cart with
Experimental Validation,” Automatica 49, no. 3 (03 2013): 801-08. doi:10.1016/j.automatica.2012.12.006.
[2] Aracil, J., J.á. Acosta, and F. Gordillo, “A Nonlinear Hybrid Controller for Swinging-up and Stabilizing the Furuta
Pendulum,” Control Engineering Practice 21, no. 8 (08 2013): 989-93. doi:10.1016/j.conengprac.2013.04.001.
[3] Åström, K.j., and K. Furuta, “Swinging up a Pendulum by Energy Control,” Automatica 36, no. 2 (02 2000): 287-
95. doi:10.1016/s0005-1098(99)00140-5.
[4] Shehu, Muhammad, Mohd Ridzuan Ahmad, Auwal Shehu, and Ahmad Alhassan, “LQR, Double-PID and Pole
Placement Stabilization and Tracking Control of Single Link Inverted Pendulum,” 2015 IEEE International
Conference on Control System, Computing and Engineering (ICCSCE), 11 2015.
doi:10.1109/iccsce.2015.7482187.
[5] S., Gopikrishnan, Ameya Anil Kesarkar, and N. Selvaganesan, “Design of Fractional Controller for Cart-pendulum
SIMO System,” 2012 IEEE International Conference on Advanced Communication Control and Computing
Technologies (ICACCCT), 08 2012. doi:10.1109/icaccct.2012.6320764.
[6] Mishra, Sunil Kumar, and Dinesh Chandra, “Stabilization and Tracking Control of Inverted Pendulum Using
Fractional Order PID Controllers,” Journal of Engineering 2014 (2014): 1-9. doi:10.1155/2014/752918.
[7] Prasad, Lal Bahadur, Barjeev Tyagi, and Hari Om Gupta, “Optimal Control of Nonlinear Inverted Pendulum
Dynamical System with Disturbance Input Using PID Controller & LQR,” 2011 IEEE International
Conference on Control System, Computing and Engineering, 11 2011. doi:10.1109/iccsce.2011.6190585.
[8] Kharola, Ashwani, Pravin Patil, Suyashi Raiwani, and Deepak Rajput, “A Comparison Study for Control and
Stabilisation of Inverted Pendulum on Inclined Surface (IPIS) Using PID and Fuzzy Controllers,” Perspectives in
Science 8 (09 2016): 187-90. doi:10.1016/j.pisc.2016.03.016.
Int J Elec & Comp Eng ISSN: 2088-8708 
Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin)
1089
[9] Tatikonda, Ravi Chandra, Venkata Praveen Battula, and Vijay Kumar, “Control of Inverted Pendulum Using
Adaptive Neuro Fuzzy Inference Structure (ANFIS),” Proceedings of 2010 IEEE International Symposium on
Circuits and Systems, 05 2010. doi:10.1109/iscas.2010.5537234.
[10] Elsayed, Belal A., Mohsen A. Hassan, and Saad Mekhilef, “Fuzzy Swinging-up with Sliding Mode Control for
Third Order Cart-inverted Pendulum System,” International Journal of Control, Automation and Systems 13, no. 1
(12, 2014): 238-48. doi:10.1007/s12555-014-0033-4.
[11] Strakos, Premysl, and Jiri Tuma, “Mathematical Modelling and Controller Design of Inverted Pendulum,” 2017
18th International Carpathian Control Conference (ICCC), 05 2017. doi:10.1109/carpathiancc.2017.7970431.
[12] Peker, Fuat, and Ibrahim Kaya, “Identification and Real Time Control of an Inverted Pendulum Using PI-PD
Controller,” 2017 21st International Conference on System Theory, Control and Computing (ICSTCC), 10 2017.
doi:10.1109/icstcc.2017.8107130.
[13] Goncharov, Nguyen Quang Dung, Than Viet Dung, Identification of object with distributed parameters, Theory and
practice in the physical, mathematical and technical sciences-XXIV International Scientific and Practice
Conference, London, May, 2012. doi:10.18007/gisap:tsca.v0i8.1421.
[14] Tin, Phu Tran, Tran Hoang Quang Minh, Tan N. Nguyen, and Miroslav Voznak, “System Performance Analysis of
Half-Duplex Relay Network over Rician Fading Channel,” TELKOMNIKA (Telecommunication, Computing,
Electronics and Control). 16, no. 1 (02, 2018): 189. doi:10.12928/telkomnika.v16i1.7491.
[15] Rashid, Tarique, Sunil Kumar, Akshay Verma, Prateek Raj Gautam, and Arvind Kumar, “Pm-EEMRP: Postural
Movement Based Energy Efficient Multi-hop Routing Protocol for Intra Wireless Body Sensor Network (Intra-
WBSN),” TELKOMNIKA (Telecommunication, Computing, Electronics and Control) 16, no. 1 (02, 2018): 166.
doi:10.12928/telkomnika.v16i1.7318.
[16] A. F. Morabito, “Power Synthesis of Mask-Constrained Shaped Beams Through Maximally-Sparse Planar
Arrays,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control), vol. 14, n. 4, pp. 1217-1219,
2016.

More Related Content

PDF
The International Journal of Engineering and Science (The IJES)
PDF
PID control dynamics of a robotic arm manipulator with two degrees of freedom.
PDF
Robot2L_IEEE00493506
PDF
Performance measurement and dynamic analysis of two dof robotic arm manipulator
PDF
Dynamics and control of a robotic arm having four links
PDF
AN INVESTIGATION ON FUZZY LOGIC CONTROLLERS (TAKAGI-SUGENO & MAMDANI) IN INVE...
PDF
Kinematics Analysis of Parallel Mechanism Based on Force Feedback Device
PDF
Ips csmc 14.06.2016
The International Journal of Engineering and Science (The IJES)
PID control dynamics of a robotic arm manipulator with two degrees of freedom.
Robot2L_IEEE00493506
Performance measurement and dynamic analysis of two dof robotic arm manipulator
Dynamics and control of a robotic arm having four links
AN INVESTIGATION ON FUZZY LOGIC CONTROLLERS (TAKAGI-SUGENO & MAMDANI) IN INVE...
Kinematics Analysis of Parallel Mechanism Based on Force Feedback Device
Ips csmc 14.06.2016

What's hot (20)

PDF
Fir 05 dynamics
PDF
Ball and beam
PDF
Image Based Visual Servoing for Omnidirectional Wheeled Mobile Robots in Volt...
PDF
Mechatronics design of ball and beam system education and research
PDF
Instrumentation and Automation of Mechatronic
PDF
Modal mass
PDF
OPTIMAL TRAJECTORY OF ROBOT MANIPULATOR FOR ENERGY MINIMIZATION WITH QUARTIC ...
PDF
PID Controller Design for a Real Time Ball and Beam System – A Double Integra...
PDF
Kinematic Model vs Dynamic Model
PDF
SLIDING MODE CONTROLLER DESIGN FOR GLOBAL CHAOS SYNCHRONIZATION OF COULLET SY...
PDF
Evaluation of Vibrational Behavior for A System With TwoDegree-of-Freedom Und...
PDF
Fractional-order sliding mode controller for the two-link robot arm
PDF
Robotics
PDF
Integral Backstepping Approach for Mobile Robot Control
PDF
SLAM of Multi-Robot System Considering Its Network Topology
PDF
PDF
IRJET-A Review Paper on using Mineral Admixture Coated Pet Fibres to Make Con...
PDF
PDF
Iaetsd modelling and controller design of cart inverted pendulum system using...
Fir 05 dynamics
Ball and beam
Image Based Visual Servoing for Omnidirectional Wheeled Mobile Robots in Volt...
Mechatronics design of ball and beam system education and research
Instrumentation and Automation of Mechatronic
Modal mass
OPTIMAL TRAJECTORY OF ROBOT MANIPULATOR FOR ENERGY MINIMIZATION WITH QUARTIC ...
PID Controller Design for a Real Time Ball and Beam System – A Double Integra...
Kinematic Model vs Dynamic Model
SLIDING MODE CONTROLLER DESIGN FOR GLOBAL CHAOS SYNCHRONIZATION OF COULLET SY...
Evaluation of Vibrational Behavior for A System With TwoDegree-of-Freedom Und...
Fractional-order sliding mode controller for the two-link robot arm
Robotics
Integral Backstepping Approach for Mobile Robot Control
SLAM of Multi-Robot System Considering Its Network Topology
IRJET-A Review Paper on using Mineral Admixture Coated Pet Fibres to Make Con...
Iaetsd modelling and controller design of cart inverted pendulum system using...
Ad

Similar to Using real interpolation method for adaptive identification of nonlinear inverted pendulum system (20)

PDF
Robust control theory based performance investigation of an inverted pendulum...
PDF
Min Max Model Predictive Control for Polysolenoid Linear Motor
PDF
Hh3312761282
PDF
STABILIZATION AT UPRIGHT EQUILIBRIUM POSITION OF A DOUBLE INVERTED PENDULUM W...
PDF
Experimental verification of SMC with moving switching lines applied to hoisti...
DOC
Analysis of Automobile Suspension
PDF
Linear quadratic regulator and pole placement for stabilizing a cart inverted...
PDF
Modern Control System (BE)
PDF
Research Inventy : International Journal of Engineering and Science
PDF
Navigation of Mobile Inverted Pendulum via Wireless control using LQR Technique
PDF
SLIDING MODE CONTROLLER DESIGN FOR GLOBAL CHAOS SYNCHRONIZATION OF COULLET SY...
PDF
Optimal FOPI-FOPD controller design for rotary inverted pendulum system using...
PDF
SLIDING MODE CONTROLLER DESIGN FOR SYNCHRONIZATION OF SHIMIZU-MORIOKA CHAOTIC...
PDF
SLIDING MODE CONTROLLER DESIGN FOR SYNCHRONIZATION OF SHIMIZU-MORIOKA CHAOTIC...
PDF
Modeling and State Feedback Controller Design of Tubular Linear Permanent Mag...
PDF
Comparison of a triple inverted pendulum stabilization using optimal control ...
PDF
Hybrid Chaos Synchronization of Hyperchaotic Newton-Leipnik Systems by Slidin...
PDF
Adaptive Neuro-Fuzzy Control Approach for a Single Inverted Pendulum System
PDF
Multi parametric model predictive control based on laguerre model for permane...
Robust control theory based performance investigation of an inverted pendulum...
Min Max Model Predictive Control for Polysolenoid Linear Motor
Hh3312761282
STABILIZATION AT UPRIGHT EQUILIBRIUM POSITION OF A DOUBLE INVERTED PENDULUM W...
Experimental verification of SMC with moving switching lines applied to hoisti...
Analysis of Automobile Suspension
Linear quadratic regulator and pole placement for stabilizing a cart inverted...
Modern Control System (BE)
Research Inventy : International Journal of Engineering and Science
Navigation of Mobile Inverted Pendulum via Wireless control using LQR Technique
SLIDING MODE CONTROLLER DESIGN FOR GLOBAL CHAOS SYNCHRONIZATION OF COULLET SY...
Optimal FOPI-FOPD controller design for rotary inverted pendulum system using...
SLIDING MODE CONTROLLER DESIGN FOR SYNCHRONIZATION OF SHIMIZU-MORIOKA CHAOTIC...
SLIDING MODE CONTROLLER DESIGN FOR SYNCHRONIZATION OF SHIMIZU-MORIOKA CHAOTIC...
Modeling and State Feedback Controller Design of Tubular Linear Permanent Mag...
Comparison of a triple inverted pendulum stabilization using optimal control ...
Hybrid Chaos Synchronization of Hyperchaotic Newton-Leipnik Systems by Slidin...
Adaptive Neuro-Fuzzy Control Approach for a Single Inverted Pendulum System
Multi parametric model predictive control based on laguerre model for permane...
Ad

More from IJECEIAES (20)

PDF
Redefining brain tumor segmentation: a cutting-edge convolutional neural netw...
PDF
Embedded machine learning-based road conditions and driving behavior monitoring
PDF
Advanced control scheme of doubly fed induction generator for wind turbine us...
PDF
Neural network optimizer of proportional-integral-differential controller par...
PDF
An improved modulation technique suitable for a three level flying capacitor ...
PDF
A review on features and methods of potential fishing zone
PDF
Electrical signal interference minimization using appropriate core material f...
PDF
Electric vehicle and photovoltaic advanced roles in enhancing the financial p...
PDF
Bibliometric analysis highlighting the role of women in addressing climate ch...
PDF
Voltage and frequency control of microgrid in presence of micro-turbine inter...
PDF
Enhancing battery system identification: nonlinear autoregressive modeling fo...
PDF
Smart grid deployment: from a bibliometric analysis to a survey
PDF
Use of analytical hierarchy process for selecting and prioritizing islanding ...
PDF
Enhancing of single-stage grid-connected photovoltaic system using fuzzy logi...
PDF
Enhancing photovoltaic system maximum power point tracking with fuzzy logic-b...
PDF
Adaptive synchronous sliding control for a robot manipulator based on neural ...
PDF
Remote field-programmable gate array laboratory for signal acquisition and de...
PDF
Detecting and resolving feature envy through automated machine learning and m...
PDF
Smart monitoring technique for solar cell systems using internet of things ba...
PDF
An efficient security framework for intrusion detection and prevention in int...
Redefining brain tumor segmentation: a cutting-edge convolutional neural netw...
Embedded machine learning-based road conditions and driving behavior monitoring
Advanced control scheme of doubly fed induction generator for wind turbine us...
Neural network optimizer of proportional-integral-differential controller par...
An improved modulation technique suitable for a three level flying capacitor ...
A review on features and methods of potential fishing zone
Electrical signal interference minimization using appropriate core material f...
Electric vehicle and photovoltaic advanced roles in enhancing the financial p...
Bibliometric analysis highlighting the role of women in addressing climate ch...
Voltage and frequency control of microgrid in presence of micro-turbine inter...
Enhancing battery system identification: nonlinear autoregressive modeling fo...
Smart grid deployment: from a bibliometric analysis to a survey
Use of analytical hierarchy process for selecting and prioritizing islanding ...
Enhancing of single-stage grid-connected photovoltaic system using fuzzy logi...
Enhancing photovoltaic system maximum power point tracking with fuzzy logic-b...
Adaptive synchronous sliding control for a robot manipulator based on neural ...
Remote field-programmable gate array laboratory for signal acquisition and de...
Detecting and resolving feature envy through automated machine learning and m...
Smart monitoring technique for solar cell systems using internet of things ba...
An efficient security framework for intrusion detection and prevention in int...

Recently uploaded (20)

PPTX
Internet of Things (IOT) - A guide to understanding
PDF
Model Code of Practice - Construction Work - 21102022 .pdf
PPTX
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PPTX
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
PPTX
Geodesy 1.pptx...............................................
PDF
Embodied AI: Ushering in the Next Era of Intelligent Systems
PPTX
bas. eng. economics group 4 presentation 1.pptx
PPTX
Construction Project Organization Group 2.pptx
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PDF
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
PPTX
OOP with Java - Java Introduction (Basics)
PDF
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
PPTX
Welding lecture in detail for understanding
PPTX
additive manufacturing of ss316l using mig welding
PDF
composite construction of structures.pdf
PPTX
Foundation to blockchain - A guide to Blockchain Tech
PPTX
Lecture Notes Electrical Wiring System Components
PPTX
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
PDF
Operating System & Kernel Study Guide-1 - converted.pdf
Internet of Things (IOT) - A guide to understanding
Model Code of Practice - Construction Work - 21102022 .pdf
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
Geodesy 1.pptx...............................................
Embodied AI: Ushering in the Next Era of Intelligent Systems
bas. eng. economics group 4 presentation 1.pptx
Construction Project Organization Group 2.pptx
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
OOP with Java - Java Introduction (Basics)
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
Welding lecture in detail for understanding
additive manufacturing of ss316l using mig welding
composite construction of structures.pdf
Foundation to blockchain - A guide to Blockchain Tech
Lecture Notes Electrical Wiring System Components
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
Operating System & Kernel Study Guide-1 - converted.pdf

Using real interpolation method for adaptive identification of nonlinear inverted pendulum system

  • 1. International Journal of Electrical and Computer Engineering (IJECE) Vol. 9, No. 2, April 2019, pp. 1078~1089 ISSN: 2088-8708, DOI: 10.11591/ijece.v9i2.pp1078-1089  1078 Journal homepage: http://iaescore.com/journals/index.php/IJECE Using real interpolation method for adaptive identification of nonlinear inverted pendulum system Phu Tran Tin1 , Tran Hoang Quang Minh2 , Tran Thanh Trang3 , Nguyen Quang Dung4 1 Faculty of Electronics Technology, Industrial University of Ho Chi Minh City, Ho Chi Minh City, Vietnam 2 Optoelectronics Research Group, Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam 3 Faculty of Electrical and Electronics Engineering, Ho Chi Minh City University of Food Industry, 140 Le Trong Tan, Ho Chi Minh City, Vietnam 4 Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam Article Info ABSTRACT Article history: Received Jun 3, 2018 Revised Oct 10, 2018 Accepted Dec 8, 2018 In this paper, we investigate the inverted pendulum system by using real interpolation method (RIM) algorithm. In the first stage, the mathematical model of the inverted pendulum system and the RIM algorithm are presented. After that, the identification of the inverted pendulum system by using the RIM algorithm is proposed. Finally, the comparison of the linear analytical model, RIM model, and nonlinear model is carried out. From the results, it is found that the inverted pendulum system by using RIM algorithm has simplicity, low computer source requirement, high accuracy and adaptiveness in the advantages. Keywords: Identification Inverted pendulum Non-linear system Real interpolation method RIM Copyright © 2019 Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Tran Hoang Quang Minh, Optoelectronics Research Group, Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, 19 Nguyen Huu Tho, District 7, Ho Chi Minh City, Vietnam. Email: tranhoangquangminh@tdtu.edu.vn 1. INTRODUCTION Inverted pendulum system, which is an underactuated mechanical system and exhibits strong nonlinearity, has been considered as an exciting and challenging problem in automation. Notably, their features such as swing up, stabilization problem, higher order, nonlinear, robust coupling, and multivariate have been attracted more and more attention from control researchers for the last decades. Moreover, inverted pendulum system is widely used in the field of robotics and aerospace field, and so has essential significance both in the field of the theoretical study research and practice [1]-[3]. There are a lot of researchers about most popular types inverted pendulums such as the rotational single arm pendulum, the cart inverted pendulum, and the double inverted investigated the inverted pendulum by switching two different control laws [3], and in [4] LQR, double-PID and pole placement control techniques are used to control a cart inverted pendulum system. Moreover, fractional order PID controllers are proposed in [5] and [6] for control of the inverted pendulum system. In [7], a simulation study on PID and LQR control of a cart inverted pendulum system is carried out with and without disturbance input. Furthermore, in [8], control of a cart inverted pendulum on an inclined surface is performed using fuzzy and PID controllers. An ANFIS controller is proposed in [9] as an intelligent control technique to control an inverted pendulum. In addition, in [10], sliding mode controller is proposed to control a cart inverted pendulum system in both simulation and real-
  • 2. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1079 time implementation. From this point of view, research in the field of inverted pendulum system is still a hot direction in both theory and practice automation area. In this paper, we investigate the inverted pendulum system by using real interpolation method (RIM) algorithm. In the first stage, the mathematical model of the inverted pendulum system and the real interpolation method (RIM) algorithm are presented. After that, the identification of the inverted pendulum system by using the real interpolation method (RIM) algorithm is proposed. Investigated and derived. Finally, the analytical algorithm derived by Simulink Mat Lab Software. From the results, it is found that the inverted pendulum system by using real interpolation method (RIM) algorithm is simplicity, low computer source requirement, high accuracy. The main contributions of this paper are presented as the following: 1) The mathematical model of the cart inverted pendulum system is proposed. 2) The identification of the inverted pendulum system by using the real interpolation method (RIM) algorithm is presented and derived. The rest of the paper is proposed as follows. The cart inverted pendulum system, the real interpolation method (RIM) and the identification of the inverted pendulum system by using the real interpolation method (RIM) algorithm is introduced in Section 2. After that, some results and discussions are given in Section 3. Finally, conclusions are provided in Section 4. 2. RESEARCH METHOD 2.1. Mathematical model of the cart inverted pendulum The cart inverted pendulum model is presented in Figure 1 [11]. Figure 1. The cart inverted pendulum model The equations of pendulum motion (a nonlinear system) can be formulated as:      2 cos sin 1M m ml x mlx b x F         , (1)      2 cos 2 sinI ml mlx b mgl        . (2) Furthermore, we try to linearize the system differential equations. It is a linearization around the equilibrium point, in our case about vertical position . The prerequisite is that the system remains in the vicinity of this position. Let presents a deviation from the equilibrium point, we get . When pendulthe um is in upright position     2 2 cos 1, sin , 0cos sin                  . Using above relation and from the nonlinear Equations (1) and (2) we obtained linearized equations of the proposed system as in (3) and (4).  2 2I ml b mgl mlx      , (3)   1x M m b x ml U    . (4)
  • 3.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089 1080 Note that U has been substituted for the input force F. To obtain the transfer functions of the linearized system equations, we must first take the Laplace transform of the system equation assuming zero initial conditions. The resulting Laplace transforms are shown below:          2 2 2 2Φ Φ ΦI ml s s b s s mgl s mlX s s    , ( (5)          2 2 1 ΦM m X s s b X s s ml s s U s    . (6) In this analysis, we set 2 ( )q I M m Mml   . (7) The transfer function of pendulum’s position can determine from systems of Equations (5), (6) and (7) as the following:             2 2 2 2 1 2 14 3 2 X 2 1 I ml b mgl s s s q q q U s M m b I ml b mgl M m b b mglb s s s s q q q                   (8) Here transfer function of pendulum’s angle can be formulated as:          2 13 2 1 Φ 2 2 1 ml s s q U s M m b I ml b mgl M m b b mglb s s s q q q           (9) We consider the equilibrium point , the pendulum become the normal form. The transfer functions of the pendulum’s angle and position are presented as in the formula (10) and (11).             2 2 2 2 1 2 14 3 2 X 2 1 I ml b mgl s s s q q q U s M m b I ml b mgl M m b b mglb s s s s q q q                   , (10)          2 13 2 1 Φ 2 2 1 ml s s q U s M m b I ml b mgl M m b b mglb s s s q q q            . (11) For the transfer function above are valid these units, we assume that:           Φ ,k v s X srad m P P s U s N U s N              , or           Φ 180 ,k v s X sdeg m P P s U s N U s N               . According to the transfer functions (8) and (9), the inverted pendulum with the equilibrium point is unstable system. Due to unstable behavior of the upright equilibrium position, identification of the inverted pendulum represents a difficult task. However, the normal pendulum with the equilibrium point is a stable system with respect to the transfer functions (10) and (11). Looking at the details, it is
  • 4. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1081 clearly seem that the absolute values of the coefficients of the transfer functions are the same for the pendulum’s upright and downright positions. The difference is only in the sign of the coefficients. Based on this conclusion, the identification task of the nonlinear pendulum system can be applied on direct pendulum characteristics. In this paper, the response characteristics of the pendulum with the downright equilibrium point, will be consider as input data for the identification task in the next subsection [11]-[12]. 2.2. Real interpolation method (RIM) The real interpolation method is one of the methods, which works on mathematical descriptions of the imaginary area. The method is based on real integral transform as follow [13]:       0 , , , 0t F f t e dt C C         (12) which assigns to the original the image as a function of the real variable . Formula of direct transform (12) can be considered as a special case of the direct Laplace transform by replacing the complex variable for real variable. Another step towards the development of the instrumentation method is the transition from continuous functions to their discrete form, using the computing resources and numerical methods. For these purposes, real interpolation method is represented by numerical characteristics { } . They are obtained as a set of values of the function in the nodes ̅̅̅̅̅, where is the number of elements numerical characteristics, called its dimension. Selecting of interpolations is a primary step in the transition to a discrete form, which has a significant impact on the numerical computing and accuracy of problem solutions. Distribution of nodes in the simplest variant is uniform. Another important advantage of real interpolation method is cross-conversion properties [11]. It dues to the fact that the behavior of the function for large values of the argument is determined mainly by the behavior of the original for small values of the variable . In the opposite case, the result is the same: the behavior of the function for small values of the argument is determined mainly by the behavior of the original for large values of the variable [13]. When considering the original of dynamic characteristics of dynamic systems, formula (12) leads to an operator model, which under certain conditions can be considered as special cases of the models based on the Laplace transform. Thus, in (12) replacing of the function by - the impulse transient function of the dynamic system, we obtain its transfer function. From here we can find the elements of a discrete model of the system, and its transfer function by performing the discretization procedures for nodes ̅̅̅̅̅     0 , 1,it iW k t e dt i       (13) Function is a real transfer function of control automatic systems, having an impulse transient response . Function could be received based on determination of transfer function such as a relationship of the imaginary of output and input signals       Y W X     (14) In which the imaginary of the output signal and the input signal is calculated from the original functions of the input and output signals:     0 t Y y t e dt     , (15)     0 t X x t e dt     , (16) where input signal, and output signal of the system.
  • 5.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089 1082 The input-output relationship of the system has a form like below:      Y W X   . (17) 2.3. Identification of the inverted pendulum based on the RIM The inverted pendulum is an unstable system, investigation of this system based on input-output signal use two kinds of the input signal, namely impulse and steps responses. For unstable systems, identification of these systems using step response is not entirely suitable, because of the unstable transient response under the effect of the step response. Using impulse response can be used, however, in practice, the idea impulse response has challenging to implement. In this paper, the conversion impulse response is used for identification of nonlinear system. The inverted pendulum is used conversion impulse response having a form of transfer function:    1 exp *s X s M s    , (18) where is magnitude, delay time. Under the effect of a conversion impulse response (18) on the output of the system, response function is represented as a function . From formula (14) and (18), in the real form, the transfer function of the inverted pendulum has a form.       1 exp * M W Y         , (19) where     0 t Y k t e dt     if response characteristics represents as a function , or     1 i N t i i Y k t e t      , is a sample of .. The linear transfer function of the angle of the inverted pendulum can be formulated as:   1 3 2 3 2 1 Φ  LR B s s s A s A s A      (20) The linear transfer function of the position of the inverted pendulum can be calculated as:   2 2 1 0 4 3 2 3 2 1 LR B s B s B s s s s s X A A A       (21) The linear transfer function of angle and position of the inverted pendulum has a form of real interpolation transformation can be formulated as (22) and (23):   1Φ 3 2 3 2 1 ΦLR B A A A           , (22)   2 2 1 0 4 3 2 3 2 1 X X X LRX B B B A A A              (23) Interestingly, the coefficients of the transfer function of the angle and position of the pendulum have the same number in the denominator. The identification process can be determined the coefficients of the denominator and separately the coefficients of the numerator. Firstly, identification of the transfer function of pendulum angle. Chosen structure and order of the linear identification like formula (20) and (21) or (22) and (23) in the real form. The number of the unknown coefficients for the angle of the pendulum is four; it means
  • 6. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1083 that the dimension of the elements numerical characteristics is four. In a simple case, the value of nodes given be calculated by nominal distribution in the range of , where  Φ 1 1 Φ Φ 1 , 2, 1 i i i             (24) Discretization of by real interpolation method in ̅̅̅̅̅̅     1Φ Φ3 2 3 2 1 1 , 1, 1 exp * i i i i i i i M B Y i A A A                  , (25) where ∫ , for the case, the transient characteristics is given a form of function with time. For the case, the transient characteristics is given a table or massive of values with sample time .    Φ 1 1 Δi i N t i i i Y k t e t      (26) where, –sampling time, N-integrator. Here, we set replacement as:       Φ1 1 , 1, 1 exp * i i i i M F Y i           (27) From (14) and (16) we have a system of the equation like:   1Φ Φ3 2 3 2 1 1 , 1,i i i i i B F i A A A             (28) or   3 2 3 2 1 1Φ Φ1 , 1,i i i i iF A A A B i           (29)      2 3 3 2 1 1Φ Φ1 1 1 , 1,i i i i i i iF A F A A B F i             (30) Rewrite the system of equations of the unknown coefficients Matrix M1:                 2 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 2 4 4 4 4 4 1 1 1 1 1 1 1 1 1 1 1 1 1 F F F F M F F F F                                  (31) Matrix N1         3 1 1 3 2 2 3 3 3 3 4 4 1 1 1 1 1 F F N F F                      (32) Matrix X1 of the unknown coefficients
  • 7.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089 1084 3 2 1 1Φ 1 A A X A B             (33) If the matrixes M1, N1 are known, the matrix X1 can be determined as the following: 1 1 1 1X M N  (34) The identification task leads to solving a system of equation of unknown coefficients of tthe ransfer function. Identification of transfer function of pendulum’s position. We already determined the coefficient of denominator n above step, in terms of simplicity, requirement determines unknown coefficients of the numerator of the transfer function of pendulum’s position. The number of characteristics is . Transfer function of pendulum’s position has the form as (23):   2 2 1 0 4 3 2 3 2 1 X X X LRX B B B A A A              (35) Obviously, in order to reduce identification task, the coefficient has already determined in the first section, the identification task remains only determine coefficients of the numerator of the transfer function of pendulum’s position . The identification procedure of transfer function of pendulum’s position illustrates following steps. If response characteristics are given a massive of time, with the sample time , value of pendulum position at and N-number of samples, numerical characteristics are calculated by following form:     1 2 Δ  i i N t i X i i Y k t e t      (36) 3. RESULTS AND DISCUSSION The inverted pendulum has parameter as in Table 1. Table 1. Parameters of the Pendulum System Parameter Value Unit Mass of the cart (M) 0.5 Kg Mass of the pendulum (m) 0.2 Kg Length of the pendulum (l) 0.3 M Coefficient of friction for the cart (b1) 0.2 Nm-1s-1 Coefficient of friction for the pendulum (b2) 0.002 Nrad-1s-1 Mass moment of inertia of the pendulum (I) 0.006 kg/m2 Gravitation force (g) 9.81 m/s2 Force applied to the cart (F) - N Position of the cart (x) - m Pendulum angle (φ) - deg. In this example we will compare results of the nonlinear modeling, describing as the formulas (1) and (2), with a linear model, solving by the analytical approximation as formula (10) and (11) and with a RIM model, which is identified by RIM. The pendulum has a set of parameters as the Table 1. On input of the models, the rectangular input acceleration signal has 1 N in time 1 s in the Figure 2(a). In the Figure 3, obviously, the RIM model and analytical linear model highly fit with nonlinear model in the beginning of the transient process. However, in the end of the transient process, the linear model represents higher in fitness comparing with the RIM model. Moreover, Figure 4 shows that the RIM model has more accuracy at the beginning of the span time [0-7 s]. However, the fitness of RIM model is reducing with the increase of time in the range 7-70 s. From the results, the maximum error of linear is 0.24 approximation by RIM is 0.45. Generally, the linear model represents a higher accuracy model than RIM
  • 8. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1085 model when the angle deviation from the downright equilibrium position is lower than 15o . The position of the pendulum versus time is presented in Figure 5. Interestingly, the RIM model of position of pendulum illustrates higher fitness level of nonlinear model than the analytical linear model. According to the absolute errors of the analytical linear model and the RIM model presented in the Figure 6. The error of the RIM model is considerably lower than of the analytical linear model in the time range [0-10 s]. In the rest of the considered settling time, the analytical model is around 1.5 times as high as the RIM model. The maximum error of the analytical and RIM models account for 0.038 and 0.026, respectively (a) (b) Figure 2. The input acceleration signal Figure 3. Response characteristics of nonlinear – green-dot line, analytic - blue-dot-line, and identification by RIM –red line The input acceleration force with 3 N, illustrated in Figure 2(b) is active on a cart during 1 second, based on the samplings of the nonlinear model, we find the approximation model by RIM for the new set of data. For identification of the deviation angle of pendulum, the nodes considered in the range , sample time and settling time 70 s. Then making a comparison with the original nonlinear model, linear model, and RIM model, we show the results as Figure 7 and Figure 8. Overall, the Figure 7 show that the RIM model represents more accuracy than the analytical linear model, especially in the range 0-10s of the settling time. Interestingly, with a high displacement of the angle from the equilibrium point , the RIM model represents more accuracy then linear model. Maximum error of the linear model and the RIM model is 9.17 and 8.62, respectively.
  • 9.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089 1086 Figure 4. Approximation error: Linearization error-red line, approximation error by RIM - blue line Figure 5. The position of the pendulum: Nonlinear model red-line, linearization model- blue line, approximation model by RIM – green line Figure 6. The error of models: linearization model – red line, approximation model by RIM – blue line
  • 10. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1087 Figure 7. The angle of the pendulum: nonlinear model – blue dot line, linear model – violet dot line, approximation by RIM – red line Figure 8. The error of angle – linear model – red line, RIM model – blue line Position of the pendulum under acting by 3 N during 1 second presents in Figure 9. Obviously, the RIM model has a higher fitness than the analytical linear model in the beginning of the settling time [0-10 s]. In the end of the settling time the fitness of both models is approximately the same. Looking at more details in the Figure 10, the fitness of the RIM model is higher than of the analytical model in time range [0-20 s]. Figure 9. The position of the pendulum: nonlinear model –red line, linear model – blue line, RIM model – green line
  • 11.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1078 - 1089 1088 Figure 10. The error of model- linearization model –red line; approximation model by RIM – blue line However, in the further transient time, the RIM model is slightly less accuracy than the analytical linear model. The maximum error of the deviation of the pendulum’s position of both models is the same at about 0.83. According to the two careful samples, the models identified by RIM represent higher accuracy than model of analytical linear approach in cases of that the deviation of pendulum accounts for higher 45o . RIM provides adaptively tool to obtain the linear transfer function of non-linear system based on transient characteristics, which has significant advantage compare analytical methods. The method which is based on expression from physical laws of dynamic pendulum. Another advantage of RIM is that the method is conducted in real domain, which is more effective compared to other methods implement in imaginary and frequency domains. 4. CONCLUSION In this paper, the inverted pendulum system by using real interpolation method (RIM) algorithm is investigated. After the identification of the inverted pendulum system by using the real interpolation method (RIM) algorithm is analyzed, the analytical algorithm derived by Simulink Mat Lab Software. From the results, it is found that the inverted pendulum system by using real interpolation method (RIM) algorithm is simplicity, low computer source requirement, high accuracy. It can be found like the prospective solution for the inverted pendulum research direction in the near future. REFERENCES [1] Glück, Tobias, Andreas Eder, and Andreas Kugi, “Swing-up Control of a Triple Pendulum on a Cart with Experimental Validation,” Automatica 49, no. 3 (03 2013): 801-08. doi:10.1016/j.automatica.2012.12.006. [2] Aracil, J., J.á. Acosta, and F. Gordillo, “A Nonlinear Hybrid Controller for Swinging-up and Stabilizing the Furuta Pendulum,” Control Engineering Practice 21, no. 8 (08 2013): 989-93. doi:10.1016/j.conengprac.2013.04.001. [3] Åström, K.j., and K. Furuta, “Swinging up a Pendulum by Energy Control,” Automatica 36, no. 2 (02 2000): 287- 95. doi:10.1016/s0005-1098(99)00140-5. [4] Shehu, Muhammad, Mohd Ridzuan Ahmad, Auwal Shehu, and Ahmad Alhassan, “LQR, Double-PID and Pole Placement Stabilization and Tracking Control of Single Link Inverted Pendulum,” 2015 IEEE International Conference on Control System, Computing and Engineering (ICCSCE), 11 2015. doi:10.1109/iccsce.2015.7482187. [5] S., Gopikrishnan, Ameya Anil Kesarkar, and N. Selvaganesan, “Design of Fractional Controller for Cart-pendulum SIMO System,” 2012 IEEE International Conference on Advanced Communication Control and Computing Technologies (ICACCCT), 08 2012. doi:10.1109/icaccct.2012.6320764. [6] Mishra, Sunil Kumar, and Dinesh Chandra, “Stabilization and Tracking Control of Inverted Pendulum Using Fractional Order PID Controllers,” Journal of Engineering 2014 (2014): 1-9. doi:10.1155/2014/752918. [7] Prasad, Lal Bahadur, Barjeev Tyagi, and Hari Om Gupta, “Optimal Control of Nonlinear Inverted Pendulum Dynamical System with Disturbance Input Using PID Controller & LQR,” 2011 IEEE International Conference on Control System, Computing and Engineering, 11 2011. doi:10.1109/iccsce.2011.6190585. [8] Kharola, Ashwani, Pravin Patil, Suyashi Raiwani, and Deepak Rajput, “A Comparison Study for Control and Stabilisation of Inverted Pendulum on Inclined Surface (IPIS) Using PID and Fuzzy Controllers,” Perspectives in Science 8 (09 2016): 187-90. doi:10.1016/j.pisc.2016.03.016.
  • 12. Int J Elec & Comp Eng ISSN: 2088-8708  Using real interpolation method for adaptive identification of nonlinear inverted… (Phu Tran Tin) 1089 [9] Tatikonda, Ravi Chandra, Venkata Praveen Battula, and Vijay Kumar, “Control of Inverted Pendulum Using Adaptive Neuro Fuzzy Inference Structure (ANFIS),” Proceedings of 2010 IEEE International Symposium on Circuits and Systems, 05 2010. doi:10.1109/iscas.2010.5537234. [10] Elsayed, Belal A., Mohsen A. Hassan, and Saad Mekhilef, “Fuzzy Swinging-up with Sliding Mode Control for Third Order Cart-inverted Pendulum System,” International Journal of Control, Automation and Systems 13, no. 1 (12, 2014): 238-48. doi:10.1007/s12555-014-0033-4. [11] Strakos, Premysl, and Jiri Tuma, “Mathematical Modelling and Controller Design of Inverted Pendulum,” 2017 18th International Carpathian Control Conference (ICCC), 05 2017. doi:10.1109/carpathiancc.2017.7970431. [12] Peker, Fuat, and Ibrahim Kaya, “Identification and Real Time Control of an Inverted Pendulum Using PI-PD Controller,” 2017 21st International Conference on System Theory, Control and Computing (ICSTCC), 10 2017. doi:10.1109/icstcc.2017.8107130. [13] Goncharov, Nguyen Quang Dung, Than Viet Dung, Identification of object with distributed parameters, Theory and practice in the physical, mathematical and technical sciences-XXIV International Scientific and Practice Conference, London, May, 2012. doi:10.18007/gisap:tsca.v0i8.1421. [14] Tin, Phu Tran, Tran Hoang Quang Minh, Tan N. Nguyen, and Miroslav Voznak, “System Performance Analysis of Half-Duplex Relay Network over Rician Fading Channel,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control). 16, no. 1 (02, 2018): 189. doi:10.12928/telkomnika.v16i1.7491. [15] Rashid, Tarique, Sunil Kumar, Akshay Verma, Prateek Raj Gautam, and Arvind Kumar, “Pm-EEMRP: Postural Movement Based Energy Efficient Multi-hop Routing Protocol for Intra Wireless Body Sensor Network (Intra- WBSN),” TELKOMNIKA (Telecommunication, Computing, Electronics and Control) 16, no. 1 (02, 2018): 166. doi:10.12928/telkomnika.v16i1.7318. [16] A. F. Morabito, “Power Synthesis of Mask-Constrained Shaped Beams Through Maximally-Sparse Planar Arrays,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control), vol. 14, n. 4, pp. 1217-1219, 2016.