SlideShare a Scribd company logo
International Journal of Electrical and Computer Engineering (IJECE)
Vol. 9, No. 2, April 2019, pp. 1067~1077
ISSN: 2088-8708, DOI: 10.11591/ijece.v9i2.pp1067-1077  1067
Journal homepage: http://iaescore.com/journals/index.php/IJECE
Multi parametric model predictive control based on laguerre
model for permanent magnet linear synchronous motors
Nguyen Hong Quang1
, Nguyen Phung2
, Dao Phuong Na3
, Nguyen Thanh Binh4
1
Departement of Automation, Thai Nguyen University of Technology, Viet Nam
2
Institude for Control Engineering and Automation, Hanoi University of Science and Technology, Viet Nam
3
Departement of Automatic Control, Hanoi University of Science and Technology, Viet Nam
4
University of Ulsan, South Korea
Article Info ABSTRACT
Article history:
Received May 26 2018
Revised Sep 13, 2018
Accepted Oct 10, 2018
The permanent magnet linear motors are widely used in various industrial
applications due to its advantages in comparisons with rotary motors such as
mechanical durability and directly creating linear motions without gears or
belts. The main difficulties of its control design are that the control
performances include the tracking of position and velocity as well as
guarantee limitations of the voltage control and its variation. In this work, a
cascade control strategy including an inner and an outer loop is applied to
synchronous linear motor. Particularly, an offline MPC controller based on
MPP method and Laguerre model was proposed for inner loop and the outer
controller was designed with the aid of nonlinear damping method. The
numerical simulation was implemented to validate performance of the
proposed controller under voltage input constraints.
Keywords:
Model predictive control
Multi parametric programming
Permanent magnet linear
synchronous motor
Polysolenoid linear motor Copyright © 2019 Institute of Advanced Engineering and Science.
All rights reserved.
Corresponding Author:
Nguyen Hong Quang,
Department of Automation,
Thai Nguyen University of Technology,
666, 3/2 Street, Tich Luong ward, Thai Nguyen city, Viet Nam.
Email: quang.nguyenhong@tnut.edu.vn
1. INTRODUCTION
The Permanent Magnet Linear Synchronous Motors (PMLSM) is extensively used in various
industries due to the ability to directly create linear motion without gears or belts. Although the mediate
mechanical actuators are eliminated, the system become weak robustness, in that external impact such as
frictional force, end – effect, changed load and non-sine of flux cause damage to control performances.
Generally, principle operation of synchronous linear motor is similar to permanent magnet synchronous
motor; however their physical construction is different in [1]-[5], and various applications such as CNC
Lathe [6], sliding door [7].
In recent years, there has been many researches for control problem of permanent magnet linear
motor. An adaptive fuzzy neural network in [8] was proposed to control the permanent magnet linear
synchronous motor. The authors in [9] presented a control design to regulate velocity based on PI – self
tuning combining with appropriate estimation technique at slow velocity zone, but if load is changed, PI –
self tuning controller will be not efficient. In order to overcome changed load, model reference control
method based on Lyapunov stability theory employed in [10]. Additionally, the backstepping technique in
[11], was applied to reduce influence of frictional force and controller is designed based on appropriate
frictional estimated model. In [12], the advantage of that the sliding mode control applied in Linear Motor is
that real position value tracks set point. However, the disadvantages of this method is that sliding surface is
complicated and chattering problem occurred. It is clear that the previous researches do not mention position,
velocity and current constraints. To solve this problem, the MPC approach in [13] was proposed as a single
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077
1068
controller for speed control. The authors in [14] built a new mathematic model and use optimal control
approach to result in linear quadratic regulation (LQR). However, the considered model did not include
disturbance load as well as friction force. In addition, the implementation of this MPC controller on a
microcontroller is very difficult because of calculation burden.
In this paper, we apply cascade control strategy to synchronous linear motor including an inner and
an outer loop. The offline MPC controller based on MPP method in [15] was proposed for inner loop to make
motor current to follow the reference signal from the outer controller. We modify optimization problem in
the MPC controller by using a Laguerre Model approach in [13] to reduce the number of optimal variables.
The major advantage of our MPC controller lies in the ability to solve constraints problem and reducing
amount of calculation because the optimal problem is offline solved. The outer controller was designed based
on nonlinear damping method in [16] to guarantee the error between real and reference velocity converge to
arbitrary small value.
2. PRIMARY RESULTS
2.1. Laguerre orthogonal polynomials
As represented in [13] Laguerre polynomials are defined as follows:
 
11
2
1
( ) (1 )
(1 )
i
i i
z a
l z a
az



 

(1)
where a is a positive constant, 0 1a  and 0,1,...i  . The application of Laguerre polynomials is mainly
in the area of system identification, in which the discrete-time impulse response of a dynamic system is
represented by a Laguerre model (see Wahlberg, 1991). In this work, based on Wang 2009, we obtain the
main result: Suppose that the impulse response of a stable system is  H k , then with a given number of
terms N ,  H k is written as:
       1 1 2 2 ... N NH k c l k c l k c l k   
(2)
with 1 2, ,..., Nc c c are the coefficients to be determined from the system data. The discrete-time Laguerre
functions are orthonormal functions, and with these orthonormal properties, the coefficients of the Laguerre
network are defined by the following relation:
     
0
, 0,1,...,i i
k
c H k l k i N


  
(3)
2.2. Model predictive control based on Laguerre function
In this section, we consider this discrete time linear system described by:
( 1) ( ) ( )m m m mx k A x k B u k  
(4)
In which, mx is vector of state variable and ( )u k is the input at the time k .
Convert the Equation (1) as follows:
( 1) ( ) ( )m m m mx k A x k B u k     
(5)
Where: ( 1) ( 1) ( ), ( ) ( ) ( 1)m m mx k x k x k u k u k u k        
From (4), (5) and by letting ( 1) [ ( 1) ( 1)]T
m mx k x k x k     , the Equation (5) becomes:
( 1) ( ) ( )x k Ax k B u k    (6)
Int J Elec & Comp Eng ISSN: 2088-8708 
Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang)
1069
With: ,
mm n n
m n n m
BA
A B
A I B


   
    
   
Predictive model of (6) at ik as follows:
   
( 1| ) ( | ) ( | )
0,1,...,
|
x k i k Ax k i k B u k i k
i N
x k k x k
       


  (7)
N is prediction horizon. From (7), the sake of designing the MPC controller is finding the sequence of input
signal ( )u k minimizing this under cost function:
1 0
( | ) ( | ) ( ) ( )
p pN N
T T
j j
J x k j k Qx k j k u k j R u k j
 
        
(8)
Where ,Q R are positive definite matries.
Denote that:
( 1| ) ( 2 | ) ( | )
( | ) ( 1| ) ( 1| )
TT T T
TT T T
X x k k x k k x k N k
U u k k u k k u k N k
     
        
(9)
Putting (7) into a form admitting variable (9) as;
0
ˆ ˆX Ax BU  (10)
Where:
2
3 2
1 2
0 0 0 0
0 0 0
ˆ ˆ, 0 0
0
N N N
A B
A AB B
A BA A B AB B
A A B A B AB B 
   
   
   
    
   
   
      
and 0 ( )x x k is initial state of prediction horizon. Substituting (10) into the cost function (8):
0 0 0
1 ˆ ˆ ˆ
2
T T T T
J U HU x FU x A QAx
 
   
  (11)
Where:
 ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ2 , 2 , ,T T T
Q R
Q R
H B QB R F A QB Q R
P R
        
               
    
   
        
,
In this work, we regard the sequence of signals      | , 1| ,..., 1|u k k u k k u k N k      can be
approximated by the by a discrete Laguerre polynomial functions. In other words, a set of Laguerre functions
     1 2, ,..., Ml k l k l k are used to capture the input ( | )iu k i k  in prediction horizon:
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077
1070
1
( | ) ( )
M
i j j
j
u k i k c l i

    (12)
Where jc are the coefficients which depend only on the initial of prediction horizon. ( )jl i are the orthogonal
Laguerre functions having form as following:
1
1
1
1 1 1
2
2 2 2
2
2
( ) ( ).... ( )
( ) ( ).... ( )
( ) ( ).... ( )
T
M
M
n n n
M
c k c k c k
c k c k c k
c k c k c k

 
 
 
  
 
 
 
,
1 2
1 2
1 2
(0) (0) (0).
(1) (1) (1).
( 1) ( 1) ( 1)
T
M
M
M
l l l
l l l
L
l N l N l N
 
 
 
 
 
   
The equation (12) is rewrited to:
( ) ( )T
iu k k L k    (13)
Substuting (13) model (4), we can obtain:
1 1
1 1
0 0
( | ) ( ) ( ) ( ) ( ) ( ) ( )
j j
j j i j j i T j T
i i i
i i
x k j k A x k A B u k i A x k A BL i A x k j  
 
   
 
         
Where
1
1
0
( ) ( )
j
T j i T
i
j A BL i

 

 
Substituting into the cost function we obtain:
   
1 0
1 0 1 1
( ) ( ) ( ) ( ) ( ( ) ) ( ( ) )
( ) ( ) ( ) ( ) 2 ( ) ( ) ( ) ( ) ( )
p p
p p p p
N N
Tj T j T T T T
i i i i
j j
N N N N
T T T j T T m m
i i i i i
j j j j
J A x k j Q A x k j L k j R L k j
j Q j L k j RL k j j QA x k x k A QA x k
     
     
 
   
     
   
           
   
 
   
Let:
1 0
( ) ( ) ( ) ( )
p pN N
T T
j j
j Q j L j RL j 
 
    ;
1
( )
pN
j
j
j QA 

 
We set the optimization problem as:
 *
argmin 2 ( )T T
ix k       (14)
 Unconstrained MPC controller
In this case, (7) was considered without any constraints on inputs and state variables:
2 2 ( )T
i
J
x k 


  

Which is equivalent:
1
* 1
1 0 1
( ) ( ) ( ) ( ) ( ) ( ) ( )
p p pN N N
T T j
i i
j j j
x k j Q j L j RL j j QA x k    


  
 
        
 
  
Int J Elec & Comp Eng ISSN: 2088-8708 
Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang)
1071
So that, the optimal input signal at time ik :
*
( ) (0)T
iu k L  
 Constrained MPC controller
We suppose that (14) has constraints as follows:
min max
min max
( )
( )m m m
u u k u
x x k x
 
 
(15)
We consider prediction model (7) and constraints (15) at time ik :
min max
min max
min max
( 1) ( | ) ( 1)
( ) ( )
( 1| )
i i i i
m m i m m i
i i
m m
u u k u k k u u k
x x k x x k
x k k
x x
      
    
     
   
(16)
(15), (16) deduce to:
min max( )iU u k U     (17)
Furthermore, by using ( ) (0)T
iu k L   , rewriting (17):
M  (18)
Thus, the optimization problem (14):
min mi
. .
n 2 ( )T T
iJ x k
s t M
   
 
  
 (19)
Since J is the quadratic function, optimal issue in (19) can be solved by Quadratic Programming (QP) to
obtain the solution *
 , the control signal is calculated by:
*
( ) (0)T
iu k L   (20)
2.3. Multi parametric programming
In this section, we remind the results of multi parametric programming method (MPP) in [11]. The
basic idea is that the space of parameter is separated into critical regions in that each critical regions, the
solution of optimization problem is in the same form. Consider the quadratic cost function:
1
min
2
. .
T
z
V z Hz
s t Gz W S
  
  
 
 
(21)
with s
z  is the vector of optimization variables. n
  is the vector of parameters, and the matrices:
s s
H 
 , q s
G 
 , q
W  , q n
S 
 . Let D be a polytopic set of parameters. In multi parametric
programming, we consider finding the solution of optimizing V on D . By using first order Karush-Kuhn-
Tucker (KKT) optimality conditions we obtain:
0,T q
Hz G    (22)
  0, 1,2,...,i i i
i G z W S i q     (23)
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077
1072
0  (24)
where i is the row index. From (22) we have:
1 T
z GH 
  (25)
Substituting into (23) we arrive at the conditions as follow:
 1
0 , 1,2,...,i iT i i
i iG H G W S i q  
     (26)
Let  and  are the Lagrange multiples corresponding to the inactive constraints and active constraints
respectively. With the inactive constraints we have 0  . With the active constraints we receive:
   
1
1 T
GH W SG 


   (27)
Where G ,W , S are the matrices corresponding to the active constraints. There exist  
1
1 T
GH G


because
the rows of G is linear independence. Substituting  from (24) into (21) we obtain:
   
1
1 1
( ) T T
z HG GG W Sx H 

 
  (28)
And from (20),  in (24) must be satisfied:
   
1
1
0T
GH W SG 


   (29)
The results (25) and (26) are the basic of multi parametric programming method in this case. Based on the
above results, the main steps of the off-line mp-QP solver are outlined in the following algorithm [11]:
Step 1: Defining the current region be the whole space D of the vector of parameter .
Step 2: Choose vector 0 in the current region.
Step 3: With 0  , find the optimal solution  0 0,z  by QP method.
Step 4: Define the active and inactive constraints in case of  0 0,z  , and then build those matrices: ,,G W S .
Step 5: Find  , z 
from (11) and (12).
Step 6: Characterize the 0CR of x from (13) in which the optimal solution is in (12).
Step 7: Redefine the current region be the 0D CR and go to step 2.
Step 8: When all regions have been explored, exit.
3. MAIN RESULTS
3.1. MPC controller for current sub-system
As mentioned in [17], current loop model of permanent magnet linear synchronous motor is
described by:
2
2 2
sqsd s sd
sd sq
sd sd sd
sq p sqs sd
sq sd
sq sq sq sq
Ldi R U
i v i
dt L L L
di UR L
i v i v
dt L L L L


 
 
  
     
 

               
(30)
Int J Elec & Comp Eng ISSN: 2088-8708 
Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang)
1073
Denote that d
sdi , d
sqi are desired outputs and ,d d
sd sd sd sq sq sqe i i e i i    . Substituting into (30):
2 2
2 2
sq sqd dsd s sd s
sd sq sd sq
sd sd sd sd sd
d
sq sq sd sd pds sd s
sq sd sq
sq sq sq sq sq
L Lde R U R
e v e i v i
dt L L L L L
de U L iR L R
e v e i v
dt L L L L L
 
 
 
 
     
           
    

    
               
(31)
Define:
212 ( )0( )
( ) , ( ) , ( )
12 20( ) ( )
sqd dsd ssqs
sd sq
sd sd sdsdsd sd
m m d
s sd sq sd sd pds
sq
sqsq sq sq sq sq
LU RLR i v t iv t
L L LLL L
A t B t u t
R L U L iRv t i v tLL L L L L


 
 
                         
                         
 
 
 
   
 
The current error model (31) is rewriten as:
( ) ( ) ( )m m
dx
A t x t B u t
dt
  (32)
Obtaining discrete time model from (32) by using ZOH method:
( 1) ( ) ( ) ( )m mx k A k x k B u k   (33)
Where: ( ) sAT
mA k e and
0
sT
m m s mB B dt T B 
Then, we can apply MPC controller designed in 2.1 to current loop and the optimal control signal is
presented in (20).
3.2. Control design for outter loop
Let us consider model of PMLSM proposed in [18]:
 
 2
c
p sq sd
p
sq sd sq
dv p
F F
dt m
F i L L i i
dx
v
dt




 


     


 

(34)
Without loss of generality, we choose desired current in d - axis on d-q coordinate: 0d
sdi  .
By letting 1 2,d dx x x x v x    with dx is desired position, we obtain the model in state space as:
1 2
2 d
x x
x au d x


  
(35)
where
2
p
p
a



 , squ i and cF
d
m
  .
Lemma 1: As is presented in [7], by using controller (36) and disturbance observer (37), the state variables of
system (35) converges to a ball centered at the origin.
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077
1074
 1 2 2 1
1 ˆtanh x tanh(x ) du L L d x
a
      
(36)
0 2
0 0 0 2
ˆ( )
( )d
d t K x
K K au x K x

 
  

    
(37)
Proof:
Let du au x  , we obtain:
 
2
0 0 0 2
0 2
1 2
2 0
ˆ( )
x x
x u d K x
K K u K x
d t K x

 



   

   

  (38)
where ˆd d d  . Choose Lyapunov candidate function as follows:
   
22 2
1 2
2
2 2 1
1 1 1
tanh tanh
2 2 2
V L x x L x d      (39)
Differentiating both side of (39) along solution of (38) with respect to t:
        
    
2
1 1 2 2 1
2
2 1 2 2 1
3 2
2
2
tanh 1 tanh tanh
1 tanh tanh
V L x x x L x u d dd
L x x L x
        
    
(40)
Substituting (36) into (40):
      
   
22
1 1 1 2 2 1
2
2 2
1 2 2 1 2 2
3
2
1
2
1
2 tanh 1 tanh tanh
tanh tanh
4
V L x x L x L x
d
k x L x L x L x dd
k
      
           
(41)
From (38):
   20 0 0 0 2 0 0 2 0
ˆd K x K K u K x K u d K x K d            
Then, refer to:
0d d K d  (42)
Rewriting (42) with the aid of (41):
        
22 2
1 1 1 2 1 2 2 1
2
0
1
3
2
1
tanh 1 tanh tanh
4
V L x x L L k x L x K d dd
k
 
              
  (43)
Int J Elec & Comp Eng ISSN: 2088-8708 
Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang)
1075
And using the relation:
2
2
2
24
d
dd k d
k
  , we finally get:
        
2
22 2
1 1 1 2 1 2 2 1 0 2
1 2
3 2
2
1
tanh 1 tanh tanh
4 4
d
V L x x L L k x L x K k d
k k
 
               
  (44)
By selecting control parameter constant 1 2 0, ,L L K such that:
1 2 1 0 2
1
1 2
1
, , 0, 0
4
L L k K k k k
k
      (45)
And assume that d   , since 2k can be chosen arbitrarily large and from (59), Lemma 1 is proved.
4. RESULTS AND ANALYSIS
Base on the above conclusions, the simulation model of PMLSM and controller are constructed in Matlab
environment. The parameters of PMLSM is given:
Parameter Value
Number of Pole 2
Pole step 72mm
Rotor mass 3.5kg
Phase coil Resistance 3.1
d-axis inductance
q- axis inductance
Flux
4.1mH
4.1mH
0.8Wb
Figure 1 and Figure 2 describe the responses of PMLSM in cases unconstrainted and constrainted
MPC controller. In first case, Figure 1. display that actual trajectory’s motor tracking designed trajectory very
fast, but it requires the large started voltage. This is a reason to we must constraint the voltage input. Figure 2.
shows responses of motor when the input voltage squ is limited by max 90squ V .
Considering the desired trajectory of motors is expressed by: ( )d
x t t , we obtain following
efficients:
Figure 1. Actual trajectory and control signal with unconstrainted MPC
 ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077
1076
Figure 2. Actual trajectory and control signal with constrainted MPC
5. CONCLUSION
In this work, we apply cascade control strategy for polysolenoid linear motor to saparate the motor
into outer and inner loop. The offline MPC controller based on MPP method was proposed for inner loop to
make motor current to follow the reference signal from the outer controller. Optimization problem in the
MPC controller was modified by using a Laguerre Model approach to reduce the number of optimal
variables. The major advantage of our MPC controller lies in the ability to solve constraints problem and
reducing amount of calculation because the optimal problem is offline solved. The outer controller was
designed based on nonlinear damping method to guarantee the error between real and reference velocity
converge to arbitrary small value.
ACKNOWLEDGEMENTS
This research was supported by Research Foundation funded by Thainguyen University of
Technology.
REFERENCES
[1] Ng. Ph. Quang, J.-A. Dittrich; “Vector Control of Three-Phase AC Machines – System Development in the
Practice,” 2nd Edition 2015, Springer Berlin Heidelberg.
[2] Jacek F. Gieras, Zbigniew J. Piech, Bronislaw Tomczuk, “Linear Synchronous Motors Transportation and
Automation Systems,” 2nd Edition. CRC press, 2011.
[3] I.Boldea; Linear Electric Machines, Drives, and MAGLEVs Handbook, CRC press, 2013.
[4] Daniel Ausderau, Polysolenoid – Linearantrieb mit genutetem Stator, Zurich. PhD Thessis, 2004.
[5] Huilai Li, Xiaomin Li, Zhiyuan Li, “Performance Assessment and Comparison of Two Types Linear Motors for
Electromagnetic Catapult,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control), 2014 Apr.,
12(4): 2506-2515
[6] Yang Zeqing, Liu Libing, Wangzuojie, Chen Yingshu, Xiao Quanyang, “Static and Dynamic Characteristic
Simulation of Feed System Driven by Linear Motor in High Speed Computer Numerical Control Lathe,”
TELKOMNIKA (Telecommunication, Computing, Electronics and Control), 2013 July, 11(7): 3673-3683.
[7] Aymen Lachheb, Jalel Khediri, Lilia El Amraoui, “Performances Analysis of a Linear Motor for Sliding Door
Application,” International Journal of Power Electronics and Drive System (IJPEDS), 2017 Sep., 8(3): 1139-1146.
[8] Faa – Jeng Lin, Po – Hung Shen, “Adaptive fuzzy-neural-network control for a DSP-based permanent magnet
linear synchronous motor servo drive,” IEEE Transactions on Fuzzy Systems, 2006, pp. 481 - 495.
[9] Jul – Ki Seok, Jong – Kun Lee, Dong – Choon Lee, “Sensorless Speed Control of Nonsalient Permanent Magnet
Synchronous Motor Using Rotor – Position – Tracking PI Controller,” IEEE Transactions on Industrial
Electronics, Vol. 53, No. 2, pp.399 – 405, 2006.
[10] Yuan – Rui Chen, Jie Wu, Nobert Cheung (2004), “Lyapunov’s Stability Theory – Based Model Reference
Adaptive Control for Permanent Magnet Linear Motor Drives,” Proc of Power Electronics Systems and
Application, 2004, pp. 260 – 266
Int J Elec & Comp Eng ISSN: 2088-8708 
Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang)
1077
[11] Chin – I Huang, Li – Chen Fu, “Adaptive Back stepping Speed/Position Control with Friction Compensation for
Linear Induction Motor,” in Proceeding of the 41st
IEEE Conference on Decision and Control, 2002, USA, pp. 474
– 479.
[12] Gerardo Tapia, Arantxa Tapia (2007), “Sliding – Mode Control for Linear Permanent – Magnet motor Position
Tracking,” Proc of the IFAC World Congress.
[13] Liuping Wang; “Model Predictive Control System Design and Implementation Using MATLAB,” Springer-Verlag
London Limited, 2009.
[14] Hossein Komijani, Saeed Masoumi Kazraji, Ehsan Baneshi, Milad Janghorban Lariche, “Modeling and State
Feedback Controller Design of Tubular Linear Permanent Magnet Synchronous Motor”, International Journal of
Power Electronics and Drive System (IJPEDS), 2016 Dec., 7(4): 1410-1419.
[15] Alexandra Grancharova, Arne Hohansen, Explicit “Nonlinear Model Predictive Control – Theory and
Applications”. Springer, 2012, ch. 1.
[16] H. K. Khalil. Nonlinear Systems. 3rd ed. Upper Saddle River, NJ:Prentice-Hall, 2002
[17] Quang H. Nguyen, Nam P. Dao, Ty T. Nguyen, Hung M. Nguyen, Hien N. Nguyen, Tan D. Vu, “Flatness Based
Control Structure for Polysolenoid Permanent Stimulation Linear Motors,” SSRG International Journal of
Electrical and Electronics Engineering, Volume-3 Issue-12, 2016, pp 31-37.
[18] Quang H. Nguyen, Nam P. Dao, Hung M. Nguyen, Hien N. Nguyen, Ty T. Nguyen, Chi P. Nguyen, “Design an
Exact Linearization Controller for Permanent Stimulation Synchronous Linear Motor Polysolenoid,” SSRG
International Journal of Electrical and Electronics Engineering, Volume-4 Issue-1, 2017, pp. 7-12.
BIOGRAPHIES OF AUTHORS
Nguyen Hong Quang received the B.S degree in electrical engineering from Thai Nguyen
University of technology (TNUT), Vietnam, in 2007, the Master’s degree in control engineering
and automation from Hanoi University of Science and Technology (HUST), Viet Nam, in 2012.
He is currently with TNUT as Ph.D student and Lecturer. His Research Interests are Electrical
Drive Systems, Adaptive Dynamic Programming Control, Robust Nonlinear Model Predictive
Control.
Nguyen Phung Quang received his Dipl.-Ing. (Uni.), Dr.-Ing. and Dr.-Ing. habil. degrees from
TU Dresden, Germany in 1975, 1991 and 1994 respectively. Prior to his return to Vietnam, he had
worked in Germany industry for many years, contributed to create inverters REFU 402 Vectovar,
RD500 (REFU Elektronik); Simovert 6SE42, Master Drive MC (Siemens). From 1996 to 1998, he
served as lecturer of TU Dresden where he was conferred as Privatdozent in 1997. He joined
Hanoi University of Science and Technology in 1999, as lecturer up to now. He is currently a
professor of HUST and honorary professor of TU Dresden. He was author/co-author of more than
170 journal and conference papers; 8 books with three among them was written in German and one
in English entitled “Vector Control of Three-Phase AC Machines – System Development in the
Practice” published by Springer in 2008, and 2nd edition in June 2015. . His Research Interests are
Electrical Drive Systems, Motion Control, Robotic Control, Vector Control of Electrical Machines,
Wind and Solar Power Systems, Digital Control Systems, Modeling and Simulation.
Dao Phuong Nam obtained Doctor degree on January – 2013 at Hanoi University of Science and
Technology (Vietnam). Currently, he holds the position as lecturer at Hanoi University of Science
and Technology, Vietnam. His research interests include control of robotic systems and
robust/adaptive, optimal control He is author/co-author of more than 70 papers (Journals,
Conferences,…)
Nguyen Thanh Binh received the B.S and M.S degree in control engineering and automation from
Hanoi University of Science and Technology (HUST), Viet Nam, in 2014 and 2017. Currently, he
holds the position as lecturer at Thuyloi University, Vietnam. Since 2018, He is currently with
University of Ulsan, Ulsan, Korea as PhD student. His research interests include Control of land,
air, underwater vehicles, and Robust Nonlinear Model Predictive Control.

More Related Content

PDF
Iaetsd estimation of damping torque for small-signal
PDF
Explicit model predictive control of fast dynamic system
PDF
B010341317
PDF
Application of neuro fuzzy controller in torque ripple minimization of vec
PDF
Image Based Visual Servoing for Omnidirectional Wheeled Mobile Robots in Volt...
PDF
Article 1
PDF
THE RESEARCH OF QUANTUM PHASE ESTIMATION ALGORITHM
PDF
922214 e002013
Iaetsd estimation of damping torque for small-signal
Explicit model predictive control of fast dynamic system
B010341317
Application of neuro fuzzy controller in torque ripple minimization of vec
Image Based Visual Servoing for Omnidirectional Wheeled Mobile Robots in Volt...
Article 1
THE RESEARCH OF QUANTUM PHASE ESTIMATION ALGORITHM
922214 e002013

What's hot (18)

PDF
OPTIMAL TRAJECTORY OF ROBOT MANIPULATOR FOR ENERGY MINIMIZATION WITH QUARTIC ...
PDF
Inverted Pendulum Control: A Brief Overview
PDF
PID Controller Design for a Real Time Ball and Beam System – A Double Integra...
PDF
Ball and beam
PDF
On tracking control problem for polysolenoid motor model predictive approach
PPT
12EE62R11_Final Presentation
PDF
A Comparison Study between Inferred State-Space and Neural Network Based Syst...
PDF
Adaptive pi based on direct synthesis nishant
PDF
Controller design of inverted pendulum using pole placement and lqr
PDF
Single Machine Power Network Load Frequency Control Using T-S Fuzzy Based Rob...
PDF
The Neural Network-Combined Optimal Control System of Induction Motor
PDF
Nonlinear Control of an Active Magnetic Bearing with Output Constraint
PDF
Iaetsd modelling and controller design of cart inverted pendulum system using...
PDF
STABILIZATION AT UPRIGHT EQUILIBRIUM POSITION OF A DOUBLE INVERTED PENDULUM W...
PDF
Study and Development of an Energy Saving Mechanical System
PDF
MODELLING AND SIMULATION OF INVERTED PENDULUM USING INTERNAL MODEL CONTROL
PDF
Multistage condition-monitoring-system-of-aircraft-gas-turbine-engine
DOCX
OPTIMAL TRAJECTORY OF ROBOT MANIPULATOR FOR ENERGY MINIMIZATION WITH QUARTIC ...
Inverted Pendulum Control: A Brief Overview
PID Controller Design for a Real Time Ball and Beam System – A Double Integra...
Ball and beam
On tracking control problem for polysolenoid motor model predictive approach
12EE62R11_Final Presentation
A Comparison Study between Inferred State-Space and Neural Network Based Syst...
Adaptive pi based on direct synthesis nishant
Controller design of inverted pendulum using pole placement and lqr
Single Machine Power Network Load Frequency Control Using T-S Fuzzy Based Rob...
The Neural Network-Combined Optimal Control System of Induction Motor
Nonlinear Control of an Active Magnetic Bearing with Output Constraint
Iaetsd modelling and controller design of cart inverted pendulum system using...
STABILIZATION AT UPRIGHT EQUILIBRIUM POSITION OF A DOUBLE INVERTED PENDULUM W...
Study and Development of an Energy Saving Mechanical System
MODELLING AND SIMULATION OF INVERTED PENDULUM USING INTERNAL MODEL CONTROL
Multistage condition-monitoring-system-of-aircraft-gas-turbine-engine
Ad

Similar to Multi parametric model predictive control based on laguerre model for permanent magnet linear synchronous motors (20)

PDF
Min Max Model Predictive Control for Polysolenoid Linear Motor
PDF
Iaetsd modelling of one link flexible arm manipulator using
PDF
An Adaptive Internal Model for Load Frequency Control using Extreme Learning ...
PDF
Fuzzy Gain-Scheduling Proportional–Integral Control for Improving the Speed B...
PDF
Sliding mode control-based system for the two-link robot arm
PDF
Back-Stepping Control of Free-Floating Space Robot based on Adaptive Neural N...
PDF
Direct torque control and dynamic performance of induction motor using fract...
PDF
Explicit model predictive control of fast dynamic system
PDF
Design of Quadratic Optimal Regulator for DC Motor
PDF
17.pmsm speed sensor less direct torque control based on ekf
PDF
Asymptotic features of Hessian Matrix in Receding Horizon Model Predictive Co...
PDF
Fractional-order sliding mode controller for the two-link robot arm
PDF
Simulation, bifurcation, and stability analysis of a SEPIC converter control...
PDF
Modelling design and control of an electromechanical mass lifting system usin...
PDF
The Control Structure for DC Motor based on the Flatness Control
PDF
Fuzzy logic based direct torque control of induction motor with space vector ...
PDF
OPTIMAL PID CONTROLLER DESIGN FOR SPEED CONTROL OF A SEPARATELY EXCITED DC MO...
PDF
EFFECT OF TWO EXOSYSTEM STRUCTURES ON OUTPUT REGULATION OF THE RTAC SYSTEM
PDF
Fuzzy logic based direct torque control of induction motor with space vector ...
PDF
Optimal PID Controller Design for Speed Control of a Separately Excited DC Mo...
Min Max Model Predictive Control for Polysolenoid Linear Motor
Iaetsd modelling of one link flexible arm manipulator using
An Adaptive Internal Model for Load Frequency Control using Extreme Learning ...
Fuzzy Gain-Scheduling Proportional–Integral Control for Improving the Speed B...
Sliding mode control-based system for the two-link robot arm
Back-Stepping Control of Free-Floating Space Robot based on Adaptive Neural N...
Direct torque control and dynamic performance of induction motor using fract...
Explicit model predictive control of fast dynamic system
Design of Quadratic Optimal Regulator for DC Motor
17.pmsm speed sensor less direct torque control based on ekf
Asymptotic features of Hessian Matrix in Receding Horizon Model Predictive Co...
Fractional-order sliding mode controller for the two-link robot arm
Simulation, bifurcation, and stability analysis of a SEPIC converter control...
Modelling design and control of an electromechanical mass lifting system usin...
The Control Structure for DC Motor based on the Flatness Control
Fuzzy logic based direct torque control of induction motor with space vector ...
OPTIMAL PID CONTROLLER DESIGN FOR SPEED CONTROL OF A SEPARATELY EXCITED DC MO...
EFFECT OF TWO EXOSYSTEM STRUCTURES ON OUTPUT REGULATION OF THE RTAC SYSTEM
Fuzzy logic based direct torque control of induction motor with space vector ...
Optimal PID Controller Design for Speed Control of a Separately Excited DC Mo...
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
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
PPTX
Geodesy 1.pptx...............................................
PDF
Model Code of Practice - Construction Work - 21102022 .pdf
PPTX
Internet of Things (IOT) - A guide to understanding
PPT
Mechanical Engineering MATERIALS Selection
PPTX
web development for engineering and engineering
PPTX
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
PPTX
Sustainable Sites - Green Building Construction
PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PPTX
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
PPTX
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
PDF
composite construction of structures.pdf
PDF
July 2025 - Top 10 Read Articles in International Journal of Software Enginee...
PDF
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PPT
Project quality management in manufacturing
PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PPTX
additive manufacturing of ss316l using mig welding
PDF
R24 SURVEYING LAB MANUAL for civil enggi
KTU 2019 -S7-MCN 401 MODULE 2-VINAY.pptx
Geodesy 1.pptx...............................................
Model Code of Practice - Construction Work - 21102022 .pdf
Internet of Things (IOT) - A guide to understanding
Mechanical Engineering MATERIALS Selection
web development for engineering and engineering
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
Sustainable Sites - Green Building Construction
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
IOT PPTs Week 10 Lecture Material.pptx of NPTEL Smart Cities contd
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
composite construction of structures.pdf
July 2025 - Top 10 Read Articles in International Journal of Software Enginee...
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Project quality management in manufacturing
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
additive manufacturing of ss316l using mig welding
R24 SURVEYING LAB MANUAL for civil enggi

Multi parametric model predictive control based on laguerre model for permanent magnet linear synchronous motors

  • 1. International Journal of Electrical and Computer Engineering (IJECE) Vol. 9, No. 2, April 2019, pp. 1067~1077 ISSN: 2088-8708, DOI: 10.11591/ijece.v9i2.pp1067-1077  1067 Journal homepage: http://iaescore.com/journals/index.php/IJECE Multi parametric model predictive control based on laguerre model for permanent magnet linear synchronous motors Nguyen Hong Quang1 , Nguyen Phung2 , Dao Phuong Na3 , Nguyen Thanh Binh4 1 Departement of Automation, Thai Nguyen University of Technology, Viet Nam 2 Institude for Control Engineering and Automation, Hanoi University of Science and Technology, Viet Nam 3 Departement of Automatic Control, Hanoi University of Science and Technology, Viet Nam 4 University of Ulsan, South Korea Article Info ABSTRACT Article history: Received May 26 2018 Revised Sep 13, 2018 Accepted Oct 10, 2018 The permanent magnet linear motors are widely used in various industrial applications due to its advantages in comparisons with rotary motors such as mechanical durability and directly creating linear motions without gears or belts. The main difficulties of its control design are that the control performances include the tracking of position and velocity as well as guarantee limitations of the voltage control and its variation. In this work, a cascade control strategy including an inner and an outer loop is applied to synchronous linear motor. Particularly, an offline MPC controller based on MPP method and Laguerre model was proposed for inner loop and the outer controller was designed with the aid of nonlinear damping method. The numerical simulation was implemented to validate performance of the proposed controller under voltage input constraints. Keywords: Model predictive control Multi parametric programming Permanent magnet linear synchronous motor Polysolenoid linear motor Copyright © 2019 Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Nguyen Hong Quang, Department of Automation, Thai Nguyen University of Technology, 666, 3/2 Street, Tich Luong ward, Thai Nguyen city, Viet Nam. Email: quang.nguyenhong@tnut.edu.vn 1. INTRODUCTION The Permanent Magnet Linear Synchronous Motors (PMLSM) is extensively used in various industries due to the ability to directly create linear motion without gears or belts. Although the mediate mechanical actuators are eliminated, the system become weak robustness, in that external impact such as frictional force, end – effect, changed load and non-sine of flux cause damage to control performances. Generally, principle operation of synchronous linear motor is similar to permanent magnet synchronous motor; however their physical construction is different in [1]-[5], and various applications such as CNC Lathe [6], sliding door [7]. In recent years, there has been many researches for control problem of permanent magnet linear motor. An adaptive fuzzy neural network in [8] was proposed to control the permanent magnet linear synchronous motor. The authors in [9] presented a control design to regulate velocity based on PI – self tuning combining with appropriate estimation technique at slow velocity zone, but if load is changed, PI – self tuning controller will be not efficient. In order to overcome changed load, model reference control method based on Lyapunov stability theory employed in [10]. Additionally, the backstepping technique in [11], was applied to reduce influence of frictional force and controller is designed based on appropriate frictional estimated model. In [12], the advantage of that the sliding mode control applied in Linear Motor is that real position value tracks set point. However, the disadvantages of this method is that sliding surface is complicated and chattering problem occurred. It is clear that the previous researches do not mention position, velocity and current constraints. To solve this problem, the MPC approach in [13] was proposed as a single
  • 2.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077 1068 controller for speed control. The authors in [14] built a new mathematic model and use optimal control approach to result in linear quadratic regulation (LQR). However, the considered model did not include disturbance load as well as friction force. In addition, the implementation of this MPC controller on a microcontroller is very difficult because of calculation burden. In this paper, we apply cascade control strategy to synchronous linear motor including an inner and an outer loop. The offline MPC controller based on MPP method in [15] was proposed for inner loop to make motor current to follow the reference signal from the outer controller. We modify optimization problem in the MPC controller by using a Laguerre Model approach in [13] to reduce the number of optimal variables. The major advantage of our MPC controller lies in the ability to solve constraints problem and reducing amount of calculation because the optimal problem is offline solved. The outer controller was designed based on nonlinear damping method in [16] to guarantee the error between real and reference velocity converge to arbitrary small value. 2. PRIMARY RESULTS 2.1. Laguerre orthogonal polynomials As represented in [13] Laguerre polynomials are defined as follows:   11 2 1 ( ) (1 ) (1 ) i i i z a l z a az       (1) where a is a positive constant, 0 1a  and 0,1,...i  . The application of Laguerre polynomials is mainly in the area of system identification, in which the discrete-time impulse response of a dynamic system is represented by a Laguerre model (see Wahlberg, 1991). In this work, based on Wang 2009, we obtain the main result: Suppose that the impulse response of a stable system is  H k , then with a given number of terms N ,  H k is written as:        1 1 2 2 ... N NH k c l k c l k c l k    (2) with 1 2, ,..., Nc c c are the coefficients to be determined from the system data. The discrete-time Laguerre functions are orthonormal functions, and with these orthonormal properties, the coefficients of the Laguerre network are defined by the following relation:       0 , 0,1,...,i i k c H k l k i N      (3) 2.2. Model predictive control based on Laguerre function In this section, we consider this discrete time linear system described by: ( 1) ( ) ( )m m m mx k A x k B u k   (4) In which, mx is vector of state variable and ( )u k is the input at the time k . Convert the Equation (1) as follows: ( 1) ( ) ( )m m m mx k A x k B u k      (5) Where: ( 1) ( 1) ( ), ( ) ( ) ( 1)m m mx k x k x k u k u k u k         From (4), (5) and by letting ( 1) [ ( 1) ( 1)]T m mx k x k x k     , the Equation (5) becomes: ( 1) ( ) ( )x k Ax k B u k    (6)
  • 3. Int J Elec & Comp Eng ISSN: 2088-8708  Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang) 1069 With: , mm n n m n n m BA A B A I B                Predictive model of (6) at ik as follows:     ( 1| ) ( | ) ( | ) 0,1,..., | x k i k Ax k i k B u k i k i N x k k x k             (7) N is prediction horizon. From (7), the sake of designing the MPC controller is finding the sequence of input signal ( )u k minimizing this under cost function: 1 0 ( | ) ( | ) ( ) ( ) p pN N T T j j J x k j k Qx k j k u k j R u k j            (8) Where ,Q R are positive definite matries. Denote that: ( 1| ) ( 2 | ) ( | ) ( | ) ( 1| ) ( 1| ) TT T T TT T T X x k k x k k x k N k U u k k u k k u k N k                (9) Putting (7) into a form admitting variable (9) as; 0 ˆ ˆX Ax BU  (10) Where: 2 3 2 1 2 0 0 0 0 0 0 0 ˆ ˆ, 0 0 0 N N N A B A AB B A BA A B AB B A A B A B AB B                                  and 0 ( )x x k is initial state of prediction horizon. Substituting (10) into the cost function (8): 0 0 0 1 ˆ ˆ ˆ 2 T T T T J U HU x FU x A QAx         (11) Where:  ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ2 , 2 , ,T T T Q R Q R H B QB R F A QB Q R P R                                            , In this work, we regard the sequence of signals      | , 1| ,..., 1|u k k u k k u k N k      can be approximated by the by a discrete Laguerre polynomial functions. In other words, a set of Laguerre functions      1 2, ,..., Ml k l k l k are used to capture the input ( | )iu k i k  in prediction horizon:
  • 4.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077 1070 1 ( | ) ( ) M i j j j u k i k c l i      (12) Where jc are the coefficients which depend only on the initial of prediction horizon. ( )jl i are the orthogonal Laguerre functions having form as following: 1 1 1 1 1 1 2 2 2 2 2 2 ( ) ( ).... ( ) ( ) ( ).... ( ) ( ) ( ).... ( ) T M M n n n M c k c k c k c k c k c k c k c k c k                 , 1 2 1 2 1 2 (0) (0) (0). (1) (1) (1). ( 1) ( 1) ( 1) T M M M l l l l l l L l N l N l N               The equation (12) is rewrited to: ( ) ( )T iu k k L k    (13) Substuting (13) model (4), we can obtain: 1 1 1 1 0 0 ( | ) ( ) ( ) ( ) ( ) ( ) ( ) j j j j i j j i T j T i i i i i x k j k A x k A B u k i A x k A BL i A x k j                     Where 1 1 0 ( ) ( ) j T j i T i j A BL i       Substituting into the cost function we obtain:     1 0 1 0 1 1 ( ) ( ) ( ) ( ) ( ( ) ) ( ( ) ) ( ) ( ) ( ) ( ) 2 ( ) ( ) ( ) ( ) ( ) p p p p p p N N Tj T j T T T T i i i i j j N N N N T T T j T T m m i i i i i j j j j J A x k j Q A x k j L k j R L k j j Q j L k j RL k j j QA x k x k A QA x k                                                   Let: 1 0 ( ) ( ) ( ) ( ) p pN N T T j j j Q j L j RL j        ; 1 ( ) pN j j j QA     We set the optimization problem as:  * argmin 2 ( )T T ix k       (14)  Unconstrained MPC controller In this case, (7) was considered without any constraints on inputs and state variables: 2 2 ( )T i J x k        Which is equivalent: 1 * 1 1 0 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) p p pN N N T T j i i j j j x k j Q j L j RL j j QA x k                         
  • 5. Int J Elec & Comp Eng ISSN: 2088-8708  Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang) 1071 So that, the optimal input signal at time ik : * ( ) (0)T iu k L    Constrained MPC controller We suppose that (14) has constraints as follows: min max min max ( ) ( )m m m u u k u x x k x     (15) We consider prediction model (7) and constraints (15) at time ik : min max min max min max ( 1) ( | ) ( 1) ( ) ( ) ( 1| ) i i i i m m i m m i i i m m u u k u k k u u k x x k x x k x k k x x                       (16) (15), (16) deduce to: min max( )iU u k U     (17) Furthermore, by using ( ) (0)T iu k L   , rewriting (17): M  (18) Thus, the optimization problem (14): min mi . . n 2 ( )T T iJ x k s t M           (19) Since J is the quadratic function, optimal issue in (19) can be solved by Quadratic Programming (QP) to obtain the solution *  , the control signal is calculated by: * ( ) (0)T iu k L   (20) 2.3. Multi parametric programming In this section, we remind the results of multi parametric programming method (MPP) in [11]. The basic idea is that the space of parameter is separated into critical regions in that each critical regions, the solution of optimization problem is in the same form. Consider the quadratic cost function: 1 min 2 . . T z V z Hz s t Gz W S           (21) with s z  is the vector of optimization variables. n   is the vector of parameters, and the matrices: s s H   , q s G   , q W  , q n S   . Let D be a polytopic set of parameters. In multi parametric programming, we consider finding the solution of optimizing V on D . By using first order Karush-Kuhn- Tucker (KKT) optimality conditions we obtain: 0,T q Hz G    (22)   0, 1,2,...,i i i i G z W S i q     (23)
  • 6.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077 1072 0  (24) where i is the row index. From (22) we have: 1 T z GH    (25) Substituting into (23) we arrive at the conditions as follow:  1 0 , 1,2,...,i iT i i i iG H G W S i q        (26) Let  and  are the Lagrange multiples corresponding to the inactive constraints and active constraints respectively. With the inactive constraints we have 0  . With the active constraints we receive:     1 1 T GH W SG       (27) Where G ,W , S are the matrices corresponding to the active constraints. There exist   1 1 T GH G   because the rows of G is linear independence. Substituting  from (24) into (21) we obtain:     1 1 1 ( ) T T z HG GG W Sx H       (28) And from (20),  in (24) must be satisfied:     1 1 0T GH W SG       (29) The results (25) and (26) are the basic of multi parametric programming method in this case. Based on the above results, the main steps of the off-line mp-QP solver are outlined in the following algorithm [11]: Step 1: Defining the current region be the whole space D of the vector of parameter . Step 2: Choose vector 0 in the current region. Step 3: With 0  , find the optimal solution  0 0,z  by QP method. Step 4: Define the active and inactive constraints in case of  0 0,z  , and then build those matrices: ,,G W S . Step 5: Find  , z  from (11) and (12). Step 6: Characterize the 0CR of x from (13) in which the optimal solution is in (12). Step 7: Redefine the current region be the 0D CR and go to step 2. Step 8: When all regions have been explored, exit. 3. MAIN RESULTS 3.1. MPC controller for current sub-system As mentioned in [17], current loop model of permanent magnet linear synchronous motor is described by: 2 2 2 sqsd s sd sd sq sd sd sd sq p sqs sd sq sd sq sq sq sq Ldi R U i v i dt L L L di UR L i v i v dt L L L L                                   (30)
  • 7. Int J Elec & Comp Eng ISSN: 2088-8708  Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang) 1073 Denote that d sdi , d sqi are desired outputs and ,d d sd sd sd sq sq sqe i i e i i    . Substituting into (30): 2 2 2 2 sq sqd dsd s sd s sd sq sd sq sd sd sd sd sd d sq sq sd sd pds sd s sq sd sq sq sq sq sq sq L Lde R U R e v e i v i dt L L L L L de U L iR L R e v e i v dt L L L L L                                                      (31) Define: 212 ( )0( ) ( ) , ( ) , ( ) 12 20( ) ( ) sqd dsd ssqs sd sq sd sd sdsdsd sd m m d s sd sq sd sd pds sq sqsq sq sq sq sq LU RLR i v t iv t L L LLL L A t B t u t R L U L iRv t i v tLL L L L L                                                                       The current error model (31) is rewriten as: ( ) ( ) ( )m m dx A t x t B u t dt   (32) Obtaining discrete time model from (32) by using ZOH method: ( 1) ( ) ( ) ( )m mx k A k x k B u k   (33) Where: ( ) sAT mA k e and 0 sT m m s mB B dt T B  Then, we can apply MPC controller designed in 2.1 to current loop and the optimal control signal is presented in (20). 3.2. Control design for outter loop Let us consider model of PMLSM proposed in [18]:    2 c p sq sd p sq sd sq dv p F F dt m F i L L i i dx v dt                    (34) Without loss of generality, we choose desired current in d - axis on d-q coordinate: 0d sdi  . By letting 1 2,d dx x x x v x    with dx is desired position, we obtain the model in state space as: 1 2 2 d x x x au d x      (35) where 2 p p a     , squ i and cF d m   . Lemma 1: As is presented in [7], by using controller (36) and disturbance observer (37), the state variables of system (35) converges to a ball centered at the origin.
  • 8.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077 1074  1 2 2 1 1 ˆtanh x tanh(x ) du L L d x a        (36) 0 2 0 0 0 2 ˆ( ) ( )d d t K x K K au x K x             (37) Proof: Let du au x  , we obtain:   2 0 0 0 2 0 2 1 2 2 0 ˆ( ) x x x u d K x K K u K x d t K x                   (38) where ˆd d d  . Choose Lyapunov candidate function as follows:     22 2 1 2 2 2 2 1 1 1 1 tanh tanh 2 2 2 V L x x L x d      (39) Differentiating both side of (39) along solution of (38) with respect to t:               2 1 1 2 2 1 2 2 1 2 2 1 3 2 2 2 tanh 1 tanh tanh 1 tanh tanh V L x x x L x u d dd L x x L x               (40) Substituting (36) into (40):            22 1 1 1 2 2 1 2 2 2 1 2 2 1 2 2 3 2 1 2 1 2 tanh 1 tanh tanh tanh tanh 4 V L x x L x L x d k x L x L x L x dd k                    (41) From (38):    20 0 0 0 2 0 0 2 0 ˆd K x K K u K x K u d K x K d             Then, refer to: 0d d K d  (42) Rewriting (42) with the aid of (41):          22 2 1 1 1 2 1 2 2 1 2 0 1 3 2 1 tanh 1 tanh tanh 4 V L x x L L k x L x K d dd k                    (43)
  • 9. Int J Elec & Comp Eng ISSN: 2088-8708  Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang) 1075 And using the relation: 2 2 2 24 d dd k d k   , we finally get:          2 22 2 1 1 1 2 1 2 2 1 0 2 1 2 3 2 2 1 tanh 1 tanh tanh 4 4 d V L x x L L k x L x K k d k k                     (44) By selecting control parameter constant 1 2 0, ,L L K such that: 1 2 1 0 2 1 1 2 1 , , 0, 0 4 L L k K k k k k       (45) And assume that d   , since 2k can be chosen arbitrarily large and from (59), Lemma 1 is proved. 4. RESULTS AND ANALYSIS Base on the above conclusions, the simulation model of PMLSM and controller are constructed in Matlab environment. The parameters of PMLSM is given: Parameter Value Number of Pole 2 Pole step 72mm Rotor mass 3.5kg Phase coil Resistance 3.1 d-axis inductance q- axis inductance Flux 4.1mH 4.1mH 0.8Wb Figure 1 and Figure 2 describe the responses of PMLSM in cases unconstrainted and constrainted MPC controller. In first case, Figure 1. display that actual trajectory’s motor tracking designed trajectory very fast, but it requires the large started voltage. This is a reason to we must constraint the voltage input. Figure 2. shows responses of motor when the input voltage squ is limited by max 90squ V . Considering the desired trajectory of motors is expressed by: ( )d x t t , we obtain following efficients: Figure 1. Actual trajectory and control signal with unconstrainted MPC
  • 10.  ISSN: 2088-8708 Int J Elec & Comp Eng, Vol. 9, No. 2, April 2019 : 1067 - 1077 1076 Figure 2. Actual trajectory and control signal with constrainted MPC 5. CONCLUSION In this work, we apply cascade control strategy for polysolenoid linear motor to saparate the motor into outer and inner loop. The offline MPC controller based on MPP method was proposed for inner loop to make motor current to follow the reference signal from the outer controller. Optimization problem in the MPC controller was modified by using a Laguerre Model approach to reduce the number of optimal variables. The major advantage of our MPC controller lies in the ability to solve constraints problem and reducing amount of calculation because the optimal problem is offline solved. The outer controller was designed based on nonlinear damping method to guarantee the error between real and reference velocity converge to arbitrary small value. ACKNOWLEDGEMENTS This research was supported by Research Foundation funded by Thainguyen University of Technology. REFERENCES [1] Ng. Ph. Quang, J.-A. Dittrich; “Vector Control of Three-Phase AC Machines – System Development in the Practice,” 2nd Edition 2015, Springer Berlin Heidelberg. [2] Jacek F. Gieras, Zbigniew J. Piech, Bronislaw Tomczuk, “Linear Synchronous Motors Transportation and Automation Systems,” 2nd Edition. CRC press, 2011. [3] I.Boldea; Linear Electric Machines, Drives, and MAGLEVs Handbook, CRC press, 2013. [4] Daniel Ausderau, Polysolenoid – Linearantrieb mit genutetem Stator, Zurich. PhD Thessis, 2004. [5] Huilai Li, Xiaomin Li, Zhiyuan Li, “Performance Assessment and Comparison of Two Types Linear Motors for Electromagnetic Catapult,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control), 2014 Apr., 12(4): 2506-2515 [6] Yang Zeqing, Liu Libing, Wangzuojie, Chen Yingshu, Xiao Quanyang, “Static and Dynamic Characteristic Simulation of Feed System Driven by Linear Motor in High Speed Computer Numerical Control Lathe,” TELKOMNIKA (Telecommunication, Computing, Electronics and Control), 2013 July, 11(7): 3673-3683. [7] Aymen Lachheb, Jalel Khediri, Lilia El Amraoui, “Performances Analysis of a Linear Motor for Sliding Door Application,” International Journal of Power Electronics and Drive System (IJPEDS), 2017 Sep., 8(3): 1139-1146. [8] Faa – Jeng Lin, Po – Hung Shen, “Adaptive fuzzy-neural-network control for a DSP-based permanent magnet linear synchronous motor servo drive,” IEEE Transactions on Fuzzy Systems, 2006, pp. 481 - 495. [9] Jul – Ki Seok, Jong – Kun Lee, Dong – Choon Lee, “Sensorless Speed Control of Nonsalient Permanent Magnet Synchronous Motor Using Rotor – Position – Tracking PI Controller,” IEEE Transactions on Industrial Electronics, Vol. 53, No. 2, pp.399 – 405, 2006. [10] Yuan – Rui Chen, Jie Wu, Nobert Cheung (2004), “Lyapunov’s Stability Theory – Based Model Reference Adaptive Control for Permanent Magnet Linear Motor Drives,” Proc of Power Electronics Systems and Application, 2004, pp. 260 – 266
  • 11. Int J Elec & Comp Eng ISSN: 2088-8708  Multi parametric model predictive control based on laguerre model… (Nguyen Hong Quang) 1077 [11] Chin – I Huang, Li – Chen Fu, “Adaptive Back stepping Speed/Position Control with Friction Compensation for Linear Induction Motor,” in Proceeding of the 41st IEEE Conference on Decision and Control, 2002, USA, pp. 474 – 479. [12] Gerardo Tapia, Arantxa Tapia (2007), “Sliding – Mode Control for Linear Permanent – Magnet motor Position Tracking,” Proc of the IFAC World Congress. [13] Liuping Wang; “Model Predictive Control System Design and Implementation Using MATLAB,” Springer-Verlag London Limited, 2009. [14] Hossein Komijani, Saeed Masoumi Kazraji, Ehsan Baneshi, Milad Janghorban Lariche, “Modeling and State Feedback Controller Design of Tubular Linear Permanent Magnet Synchronous Motor”, International Journal of Power Electronics and Drive System (IJPEDS), 2016 Dec., 7(4): 1410-1419. [15] Alexandra Grancharova, Arne Hohansen, Explicit “Nonlinear Model Predictive Control – Theory and Applications”. Springer, 2012, ch. 1. [16] H. K. Khalil. Nonlinear Systems. 3rd ed. Upper Saddle River, NJ:Prentice-Hall, 2002 [17] Quang H. Nguyen, Nam P. Dao, Ty T. Nguyen, Hung M. Nguyen, Hien N. Nguyen, Tan D. Vu, “Flatness Based Control Structure for Polysolenoid Permanent Stimulation Linear Motors,” SSRG International Journal of Electrical and Electronics Engineering, Volume-3 Issue-12, 2016, pp 31-37. [18] Quang H. Nguyen, Nam P. Dao, Hung M. Nguyen, Hien N. Nguyen, Ty T. Nguyen, Chi P. Nguyen, “Design an Exact Linearization Controller for Permanent Stimulation Synchronous Linear Motor Polysolenoid,” SSRG International Journal of Electrical and Electronics Engineering, Volume-4 Issue-1, 2017, pp. 7-12. BIOGRAPHIES OF AUTHORS Nguyen Hong Quang received the B.S degree in electrical engineering from Thai Nguyen University of technology (TNUT), Vietnam, in 2007, the Master’s degree in control engineering and automation from Hanoi University of Science and Technology (HUST), Viet Nam, in 2012. He is currently with TNUT as Ph.D student and Lecturer. His Research Interests are Electrical Drive Systems, Adaptive Dynamic Programming Control, Robust Nonlinear Model Predictive Control. Nguyen Phung Quang received his Dipl.-Ing. (Uni.), Dr.-Ing. and Dr.-Ing. habil. degrees from TU Dresden, Germany in 1975, 1991 and 1994 respectively. Prior to his return to Vietnam, he had worked in Germany industry for many years, contributed to create inverters REFU 402 Vectovar, RD500 (REFU Elektronik); Simovert 6SE42, Master Drive MC (Siemens). From 1996 to 1998, he served as lecturer of TU Dresden where he was conferred as Privatdozent in 1997. He joined Hanoi University of Science and Technology in 1999, as lecturer up to now. He is currently a professor of HUST and honorary professor of TU Dresden. He was author/co-author of more than 170 journal and conference papers; 8 books with three among them was written in German and one in English entitled “Vector Control of Three-Phase AC Machines – System Development in the Practice” published by Springer in 2008, and 2nd edition in June 2015. . His Research Interests are Electrical Drive Systems, Motion Control, Robotic Control, Vector Control of Electrical Machines, Wind and Solar Power Systems, Digital Control Systems, Modeling and Simulation. Dao Phuong Nam obtained Doctor degree on January – 2013 at Hanoi University of Science and Technology (Vietnam). Currently, he holds the position as lecturer at Hanoi University of Science and Technology, Vietnam. His research interests include control of robotic systems and robust/adaptive, optimal control He is author/co-author of more than 70 papers (Journals, Conferences,…) Nguyen Thanh Binh received the B.S and M.S degree in control engineering and automation from Hanoi University of Science and Technology (HUST), Viet Nam, in 2014 and 2017. Currently, he holds the position as lecturer at Thuyloi University, Vietnam. Since 2018, He is currently with University of Ulsan, Ulsan, Korea as PhD student. His research interests include Control of land, air, underwater vehicles, and Robust Nonlinear Model Predictive Control.