Linear Matrix Inequality-Based Robust Model Predictive Speed Control for a Permanent Magnetic Synchronous Motor with a Disturbance Observer

: In this paper, a linear matrix inequality (LMI)-based robust model predictive speed control (RMPSC) with a disturbance observer (DOB) is proposed to guarantee stability and control performance against the parameter uncertainty and disturbance of a permanent magnetic synchronous motor (PMSM). All external torques applied to the PMSM are defined as disturbance, estimated by the DOB


Introduction
A permanent magnet synchronous motor (PMSM) is a type of electric motor whose utilization is increasing due to its various advantages [1,2] such as high efficiency, compact design, precise control, high torque-to-inertia ratio, and high power density.Therefore, PMSMs are used in a variety of fields [3][4][5][6] including CNC machines and robotics, which traditionally require precise control, as well as factory automation equipment, new and renewable energy production, space-and defense-related fields, and, recently, even for electric power steering (EPS) in the automobile industry.However, a PMSM is a timevarying system with nonlinearity that is mechanically coupled with other systems and operates with complex dynamics, causing many difficulties [7][8][9] in control system design.
To efficiently control such complex PMSM systems, control methods such as modified PI control [10][11][12][13][14], sliding mode control (SMC) [15][16][17][18], model predictive control (MPC) [19][20][21], and deep learning (DL)-based control [22][23][24] have been proposed.In [11][12][13][14], modified PI control-based methods were implemented.A vector-controlled PMSM drive using a hybrid fuzzy PI speed controller [11] that combines the advantages of a PI controller and a fuzzy controller was proposed.It is characterized by being a combination of two outputs to which a weight for a predetermined speed error or a weight for the actual speed error is assigned.In [12], a novel adaptive fuzzy proportional integral (AFPI) control method was introduce to address the system's uncertain saturation, uncertain disturbances, and system delay.An online optimization adaptive fuzzy proportional integral (AFPI) control method was applied using the designed adaptive fuzzy tuner.On the other hand, Ref. [13] proposed a mixed PI controller consisting of angle compensation, complex coefficient, and proportional integrator controllers that provide decoupling performance for the delay time of digital control.In addition, a neural network-based PID controller is presented in [14].Four control modules were combined to enable optimal gain adjustment of the online PID controller while considering nonlinearity.However, the PI control method generally has the limitation that it can only satisfy control performance within a certain range.
Meanwhile, several studies have been carried out to apply sliding mode control (SMC) [15][16][17][18], a method capable of effectively regulating nonlinear systems, in PMSMs.Simple model information is required to construct the SMC method, but it is used in many fields because it shows excellent control performance against disturbances.An extended state observer (ESO)-based anti-disturbance sliding-mode speed controller [15] was proposed.This method alleviates the chattering of the SMC and improves the response characteristics of the speed control, allowing the sliding mode surface to be reached quickly.A hybrid SMC combining sliding-mode-and deadbeat-based predictive current control was proposed in [16], and the results show robustness to variation in load and machine parameters by an integrated high-order sliding-mode observer.In [17], a singular-perturbationbased second-order sliding-mode controller that separates the slow-speed section and high-speed section was proposed.Depending on the speed of the PMSM, it was determined whether or not to apply a differentiator, and this resulted in a chattering reduction effect.There is also research on SMCs applying fuzzy control [18], considering disturbance inputs that represent PMSM system nonlinearity or unmodeled uncertainty.
As high-performance computing devices become more popular and cheaper, research combining model predictive control (MPC) with a PMSM [19][20][21] is also being actively conducted.In general, MPC can be divided into a finite control set (FCS) MPC and a continuous control set (CCS) MPC according to the control input space.FCS-MPC performs optimization based on a limited number of switching states and is widely used because it has a lower computational load and is easier to implement compared to CCS-MPC.A direct-torque control method [19], consisting of eight voltage space vectors of a two-level inverter has been presented, and torque ripple was reduced by using Lyapunov's theory to determine the weights of the cost function.In [20], a generalized multiple-vector-based MPC for PMSM drives, which is integrated into one frame, was studied to reduce the amount of computation compared to existing MPCs.To reduce the amount of computation as in [20], a combination of deadbeat and MPC [21] was proposed to optimize the switchingvector selection procedure.
Recently, research on applying deep learning to PMSMs [22][23][24] has been conducted.To improve the performance of PMSMs in the electric aircraft field, a novel active disturbance rejection control based on deep reinforcement learning [22] was proposed.It trains a neural network using the twin-delayed deep deterministic (TDDD) policy gradient algorithm and performs optimization of the DRL model.Ref. [23] proposed a multi-scenario parameter optimization method that applied deep reinforcement learning to the optimization process of finding the parameters of the active disturbance rejection controller.In [24], a deep neural network (DNN)-based long-horizon MPC was proposed to achieve offset-free torque tracking.It highlights the practical advantages of utilizing MPC approximation in real-time through a DNN, specifically emphasizing the elimination of the need for solver development in real-time implementation.
This paper proposes a linear matrix inequality (LMI)-based robust model predictive speed control (RMPSC) with a disturbance observer (DOB) to regulate the speed of a PMSM that is robust in the face of unknown disturbances and parameter variations.All external torques connected with the PMSM are defined as disturbances and estimated by the DOB.In [7], a disturbance observer was designed based on the assumption that disturbances change slowly, and in [25] one was designed under the assumption that disturbances decrease.Thus, the estimated performance is limited for rapid changes in disturbances, and no specific method was provided to obtain disturbance observer gain.However, in order to ensure the estimated performance of a disturbance, the proposed DOB derives multiple LMIs based on the Lyapunov function and formalizes them as an optimization problem to obtain optimal disturbance gain.And the proposed RMPSC method with a DOB regulates the speed of the PMSM considering the uncertainty of the parameters.In addition, we propose a gain-tuning procedure to find the optimal gains of observer and controller.This proposed method is expected to be suitable for stable operation in environments where the PMSM's parameters change due to long-term operation or torque changes rapidly [26,27].
The contribution of this study is to design a DOB to be used for prediction and control, and to demonstrate that the RMPSC method can stably control the speed of a PMSM with parameter uncertainty within a certain range due to manufacturing errors or deterioration caused by aging.And an optimization-based systematic-tuning method is proposed to find the optimal observer and controller gains.Additionally, we propose an offline optimization method that ensures a low computational load to enable real-time processing despite CCS-MPC.
The remainder of this paper is organized as follows.The system description of the PMSM is provided in Section 2. In Section 3, the disturbance observer and state observer are presented as the optimization methods.Section 4 explains the proposed RMPSC method, including parameter uncertainty and disturbance.The tuning procedure and simulation results are presented in Section 5. Finally, the conclusion is summarized in Section 6.

System Description
In many papers, the mathematical model of a PMSM is converted from a three-phase fixed-coordinate system to a two-axis rotation coordinate system for convenience of analysis and control [15].The generalized nonlinear system dynamics of a PMSM can be represented in the synchronous reference frame (SRF) as follows: where V d , V q , I d , I q , R s , Φ d , Φ q , and ω m are the stator voltage, stator current, stator resistance, stator flux linkage, and rotor speed, respectively.The PMSM stator flux linkage is expressed as: where L d and I q are the stator inductance and λ r is the flux linkage of a permanent magnet, respectively.
The electromagnetic torque equation of the PMSM can be expressed as: where P is pole pairs of a PMSM.Since a surface mounted type of PMSM (SPMSM) is used in this study, the value of L d and L q are the same.Therefore, the dynamics of electromagnetic torque are rewritten as: The mechanical motion of the PMSM is expressed as: ) where θ m , F and τ m are the rotor position, viscosity coefficient, and mechanical torque, and δ is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance d.Substituting ( 6) in (7) yields where K Energies 2024, 17, 869 4 of 17 where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: where () =   , () = [ ] , () = [] , and () =  are the state vector, control input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space model in the continuous time of (13) can be converted to the following linear state-space model in discrete time using step time  : where

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into contact with an external system, the speed of the PMSM may change due to external mechanical torque.Therefore, if the external mechanical torque applied to the PMSM is predictable, the speed regulation performance of the PMSM can be improved.Meanwhile, the mechanical encoder of the PMSM has many limitations in installation and use due to wear and tear, lifespan, sensor sensitivity, and cost issues depending on the usage environment.Therefore, estimation of the PMSM's speed, which is the role of the encoder, is needed.where  =  .
To eliminate the interaction between  ,  , and  , the reference c ally assumed to be zero [28].Therefore, the only control input that regu the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented space model as follows: where () =   , () = [ ] , () = [] , and () =  are the trol input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space mod ous time of (13) can be converted to the following linear state-space mod using step time  : where  ≔   ,  ≔    ,  ≔    ,  ≔

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes in external system, the speed of the PMSM may change due to external m Therefore, if the external mechanical torque applied to the PMSM is pred regulation performance of the PMSM can be improved.Meanwhile, th coder of the PMSM has many limitations in installation and use due lifespan, sensor sensitivity, and cost issues depending on the usage env fore, estimation of the PMSM's speed, which is the role of the encoder, i −τ m + δ.The above Equations ( 1), ( 2) and ( 9) can be represented as follows: where ω e = Pω m .
To eliminate the interaction between ω e , I d , and I q , the reference current I d is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage V q .Equations ( 11) and ( 12) can be presented as a linear state-space model as follows: .
where x(t) = I q ω e T , u(t) = [V q ], ϵ(t) = [ϵ], and y(t) = I q are the state vector, control input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space model in the continuous time of (13) can be converted to the following linear state-space model in discrete time using step time T s : x where where  =  .
To eliminate the interaction between  ,  , and  , the reference current  is us ally assumed to be zero [28].Therefore, the only control input that regulates the speed the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear stat space model as follows: () = () + () +   () () = () (13 where  =  . To eliminate the interaction between  ,  , and  , the reference curren ally assumed to be zero [28].Therefore, the only control input that regulates the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a space model as follows: where  =  .
To eliminate the interaction between  ,  , and  , th ally assumed to be zero [28].Therefore, the only control inp the PMSM is stator voltage  .Equations (11)  where  =  .
To eliminate the interaction between ally assumed to be zero [28].Therefore, t the PMSM is stator voltage  .Equation space model as follows:

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into contact with an external system, the speed of the PMSM may change due to external mechanical torque.Therefore, if the external mechanical torque applied to the PMSM is predictable, the speed regulation performance of the PMSM can be improved.Meanwhile, the mechanical encoder of the PMSM has many limitations in installation and use due to wear and tear, lifespan, sensor sensitivity, and cost issues depending on the usage environment.Therefore, estimation of the PMSM's speed, which is the role of the encoder, is needed.
In this section, we propose the disturbance observer and the state observer design method to estimate the disturbance ϵ and the state x, which are denoted by ε and x, respectively.
Let us assume that the disturbance ϵ(k) is a constant: Energies 2024, 17, 869 5 of 17 It can be observed as where N 2 is disturbance observer gain.And the full state observer of x is proposed as follows using state observer gain N 1 and the observed disturbance ε.
Subtracting the state observer (17) from state ( 14) and the disturbance observer ( 16) from disturbance (15) gives the following augmented equation:

z(k)
Energies 2024, 17, 869 4 of 1 where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque and  is the unknown external torque, respectively.Since the mechanical torque and un known external torque cannot be measured in the system, these values are defined as dis turbance .Substituting ( 6) in ( 7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usu ally assumed to be zero [28].Therefore, the only control input that regulates the speed o the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear state space model as follows: where () =   , () = [ ] , () = [] , and () =  are the state vector, con trol input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space model in the continu ous time of ( 13) can be converted to the following linear state-space model in discrete tim using step time  : where  ≔   ,  ≔    ,  ≔    ,  ≔

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into contact with a external system, the speed of the PMSM may change due to external mechanical torque Therefore, if the external mechanical torque applied to the PMSM is predictable, the spee regulation performance of the PMSM can be improved.Meanwhile, the mechanical en coder of the PMSM has many limitations in installation and use due to wear and tea lifespan, sensor sensitivity, and cost issues depending on the usage environment.There fore, estimation of the PMSM's speed, which is the role of the encoder, is needed.
Energies 2024, 17, 869  =    where  ,  and  are the rotor position, viscosity coefficient, and mec and  is the unknown external torque, respectively.Since the mechanical known external torque cannot be measured in the system, these values are turbance .
Substituting ( 6) in (7) yields The above Equations ( 1), (2), and ( 9) can be represented as follows: To eliminate the interaction between  ,  , and  , the reference cur ally assumed to be zero [28].Therefore, the only control input that regulat the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as space model as follows: , and () =  are the sta trol input, and measurement output, respectively.
Utilizing the Euler approximation [29], the linear state-space model ous time of ( 13) can be converted to the following linear state-space model i using step time  :

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into c external system, the speed of the PMSM may change due to external mec Therefore, if the external mechanical torque applied to the PMSM is predict regulation performance of the PMSM can be improved.Meanwhile, the m coder of the PMSM has many limitations in installation and use due to lifespan, sensor sensitivity, and cost issues depending on the usage enviro fore, estimation of the PMSM's speed, which is the role of the encoder, is n To eliminate the interaction between  ,  , and  ally assumed to be zero [28].Therefore, the only control the PMSM is stator voltage  .Equations ( 11) and ( 12) c space model as follows: , and () = trol input, and measurement output, respectively.
Utilizing the Euler approximation [29], the linear s ous time of (13) can be converted to the following linear s using step time  :

Disturbance and State Observer Design
As the mechanical rotating body connected to the P external system, the speed of the PMSM may change du Therefore, if the external mechanical torque applied to th regulation performance of the PMSM can be improved coder of the PMSM has many limitations in installation lifespan, sensor sensitivity, and cost issues depending o fore, estimation of the PMSM's speed, which is the role o To eliminate the interaction between  ally assumed to be zero [28].Therefore, the the PMSM is stator voltage  .Equations (1 space model as follows:

and measurement output, respect
Utilizing the Euler approximation [29], ous time of (13) can be converted to the follow using step time  :

Disturbance and State Observer Design
As the mechanical rotating body conne external system, the speed of the PMSM ma Therefore, if the external mechanical torque a regulation performance of the PMSM can b coder of the PMSM has many limitations i lifespan, sensor sensitivity, and cost issues d fore, estimation of the PMSM's speed, which where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: , and () =  are the state vector, control input, and measurement output, respectively.
Utilizing the Euler approximation [29], the linear state-space model in the continuous time of (13) can be converted to the following linear state-space model in discrete time using step time  : where  ≔   ,  ≔    ,  ≔    ,  ≔

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into contact with an external system, the speed of the PMSM may change due to external mechanical torque.Therefore, if the external mechanical torque applied to the PMSM is predictable, the speed regulation performance of the PMSM can be improved.Meanwhile, the mechanical encoder of the PMSM has many limitations in installation and use due to wear and tear, lifespan, sensor sensitivity, and cost issues depending on the usage environment.There- To eliminate the interaction between  ,  , and  , the reference ally assumed to be zero [28].Therefore, the only control input that regu the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented space model as follows: , and () =  are the trol input, and measurement output, respectively.
Utilizing the Euler approximation [29], the linear state-space mod ous time of (13) can be converted to the following linear state-space mod using step time  :

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes in external system, the speed of the PMSM may change due to external m Therefore, if the external mechanical torque applied to the PMSM is pred regulation performance of the PMSM can be improved.Meanwhile, t coder of the PMSM has many limitations in installation and use due lifespan, sensor sensitivity, and cost issues depending on the usage en Suppose the Lyapunov function, where P is a positive definite matrix, is defined as follows: In Equation ( 19) above, the disturbance and state observer gain N is determined so that V(k) decreases [30] and it can be expressed as follows: where △V(k + 1) Energies 2024, 17, 869 4 of 17 where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: , and () =  are the state vector, control input, and measurement output, respectively.
Utilizing the Euler approximation [29], the linear state-space model in the continuous time of (13) can be converted to the following linear state-space model in discrete time using step time  : where 20) is satisfied under the condition that z(k) is non-zero, and the following equation is fulfilled: [A e − NC e ] T P[A e − NC e ] − P < 0 (21) If there exists a P o satisfying P > P o > 0, and Equation (22), it is clear that Equation (21) holds.
In (20), in order to reduce the Lyapunov function V(k) as quickly as possible at each step, the value of [A e − NC e ] T P[A e − NC e ] must become smaller.Therefore, to control the decay rate of the Lyapunov function, the value α is suggested as follows: where 0 < α < 1.
In Equation (23), it is evident that a decrease in the value of alpha α leads to a faster reduction in the Lyapunov function V(k).And by applying Schur supplementary lemma [28], it can be converted as an LMI as follows.

P o
[ where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .

PN.
Energies 2024, 17, 869 6 of 17 Assuming that the parameters of a PMSM have uncertainty relative to their nominal values due to manufacturing errors or aging degradation, the parameters containing uncertainty can be expressed as values within the following range: where R o s , L o q , J o and λ o r are nominal values and µ ≥ 1.The matrices A, B and B d containing parameter R o s , L o q , J o and λ o r belong to the polyhedral uncertainty set: Therefore, considering parameter uncertainty, Equation ( 24) is represented as the following LMIs for (A i , B i , B d,i ) ∈ Ω: where A e,i Energies 2024, 17, 869 4 of 17 where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
The above Equations ( 1), (2), and ( 9) can be represented as follows: where  =  .To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: where () =   , () = [ ] , () = [] , and () =  are the state vector, control input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space model in the continuous time of ( 13) can be converted to the following linear state-space model in discrete time using step time  : where  ≔   ,  ≔    ,  ≔    ,  ≔

Disturbance and State Observer Design
The observer gain N, which allows for the fastest convergence of the estimation errors of the disturbance observer and the state observer, can be obtained through the following optimization problem: Through the optimization problem in (28), the decision matrices P, P o , and Y are calculated, and as a result, the gain matrices of the disturbance observer and state observer are determined as follows.minimize P,P o ,Y α subject to (23) and (27). ( The decision matrix P, P o , and Y are calculated via the optimization problem in (28) and the gain matrix of N is determined by

Robust Model Predictive Speed Controller (RMPSC) Design
The design of a robust model predictive speed controller (RMPSC) to regulate the speed of a PMSM under parameter uncertainty and disturbance is presented in this section.At each step, the reference speed ω * e is assumed to be known and optimal control gain is determined by optimization of the problem of cost function.

Steady-State Condition
Since the proposed RMPSC method is based on a cost function, defining this cost function requires calculating the steady state x 0 of the PMSM from the reference input.Therefore, substituting the reference speed ω * e into Equation ( 14) to obtain x 0 is as follows: where x 0 = [I o q ω o e ] T and u 0 = [V o q ] are the steady state and steady state control input, respectively.In (30), the steady-state speed ω o e of the PMSM is the same as the reference speed value ω * e .Since the disturbance torque ϵ is unmeasurable, it will be replaced by ε.Thus, Equation ( 30) is expressed as: The steady-state currents I o q of x 0 and steady-state control input V o q of u 0 are determined using the reference speed of PMSM ω * e :

Cost Function of the RMPSC
The cost function J(k) can be defined as follows based on the difference between the steady state x 0 and the steady-state control input u 0 : where  =  .

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into c external system, the speed of the PMSM may change due to external mech Therefore, if the external mechanical torque applied to the PMSM is predict regulation performance of the PMSM can be improved.Meanwhile, the m coder of the PMSM has many limitations in installation and use due to w lifespan, sensor sensitivity, and cost issues depending on the usage enviro fore, estimation of the PMSM's speed, which is the role of the encoder, is n where ψ(k + 1) Energies 2024, 17, 869 4 of 17 where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: () = () + () +   () () = () where () =   , () = [ ] , () = [] , and () =  are the state vector, control input, and measurement output, respectively.Utilizing the Euler approximation [29], the linear state-space model in the continuous time of ( 13) can be converted to the following linear state-space model in discrete time using step time  : where  ≔   ,  ≔    ,  ≔    ,  ≔

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM comes into contact with an external system, the speed of the PMSM may change due to external mechanical torque.Therefore, if the external mechanical torque applied to the PMSM is predictable, the speed regulation performance of the PMSM can be improved.Meanwhile, the mechanical encoder of the PMSM has many limitations in installation and use due to wear and tear, lifespan, sensor sensitivity, and cost issues depending on the usage environment.There- where  =  .

Disturbance and State Observer Design
As the mechanical rotating body connected to the PMSM co external system, the speed of the PMSM may change due to ext Therefore, if the external mechanical torque applied to the PMSM regulation performance of the PMSM can be improved.Meanw coder of the PMSM has many limitations in installation and u lifespan, sensor sensitivity, and cost issues depending on the us where  =  .
To eliminate the interaction between  ally assumed to be zero [28].Therefore, the o the PMSM is stator voltage  .Equations (11 space model as follows: where trol input, and measurement output, respecti Utilizing the Euler approximation [29], ous time of ( 13) can be converted to the follow using step time  : where  ≔   ,  ≔    ,  ≔

Disturbance and State Observer Design
As the mechanical rotating body connec external system, the speed of the PMSM may Therefore, if the external mechanical torque a regulation performance of the PMSM can be coder of the PMSM has many limitations in lifespan, sensor sensitivity, and cost issues d x T Hx, and H > 0, R > 0.Then, we propose the control law as follows: Using Equations ( 17) and (31), ψ(k + 1) can be expressed as follows: Assuming that x(k) converges to x(k), Equation ( 35) can be derived using the control law (34).
By substituting (34) and ( 35) into (33), the following cost function J(k) can be obtained: In order to achieve the control objective, the cost function J(k) should decrease at each step, so the following equation, including the controller gain k, should be satisfied: where △J(k) where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
To eliminate the interaction between  ,  , and  , the reference current  is usually assumed to be zero [28].Therefore, the only control input that regulates the speed of the PMSM is stator voltage  .Equations ( 11) and ( 12) can be presented as a linear statespace model as follows: 1 holds.Multiplying the matrix Q on the left and right side of (39), respectively, we obtain where where  ,  and  are the rotor position, viscosity coefficient, and mechanical torque, and  is the unknown external torque, respectively.Since the mechanical torque and unknown external torque cannot be measured in the system, these values are defined as disturbance .Substituting ( 6) in (7) yields where  ≔  and disturbance  ≔ − + .
where  ,  and  are the rotor position, viscosity coefficient, and mechanical torqu and  is the unknown external torque, respectively.Since the mechanical torque and u known external torque cannot be measured in the system, these values are defined as di turbance .Substituting ( 6) in (7)  where  =  .
KQ. Additionally, by applying the method used in previous the Equations ( 21) and (22) to Equation (40), we can derive the following equation with where  ≔   ,  ≔  .Additionally, by applying the method used in previo Equations ( 21) and (22) to Equation (40), we can derive the following equation with the range of 0 <   < .

𝐐 𝐨 < β𝐐
where 0 < β < 1.By applying the Schur supplementary lemma [30] previously u the observer design, this can be expressed as the following LMI: Considering the parameter uncertainty in ( 26), (43) must be satisfied f  ,  ,  , ∈  using the following LMIs: Having the minimum β value means ensuring that J(k) converges to 0 the f Therefore, the value of the control gain matrix  can be determined through the f ing optimization problem: minimize , , β subject to (42) and (44).
The optimization problem (45) can be solved using the same procedure as i using the design parameter .Here, the minimum value of β secures that  conv fastest to  and  is determined using the following decision matrices  and :  =   Consequently, the DOB and RMPCS proposed in Sections 3 and 4 of this pap illustrated as shown in Figure 1.

Simulation Result
The simulation was performed with several cases and the results are presented to demonstrate the effectiveness of the proposed RMPSC method in this section.The nominal value of the PMSM parameters used in the simulation are listed in Table 1.As per the characteristics of the surface mounted PMSM used in this simulation, the stator inductance L d and L q values are the same.In this study, parameter uncertainties are taken into account to ensure robust control of PMSM speed against manufacturing errors or aging degradation.In (25), the range of parameter uncertainties such as stator inductance L d , stator resistance R s , rotational inertia J, and flux linkage of permanent magnets λ r are expressed with uncertainty ratio µ(≥ 1).As mentioned in Equation (25), we consider uncertainties for four parameters, so the total number of corner points of uncertainty Ω is 16.Based on uncertainty ratio µ, the disturbance and state observer gain and controller gain are determined via the optimization problem with multiple LMIs in ( 27) and (44).In other words, the proposed RMPSC method is able to optimize control performance using uncertainty ratio µ.
A large value of µ means that the parameter values have a large range compared to the rated value.This can result in slow response characteristics because the performance of the RMPSC method must be guaranteed over parameters containing a wide range of uncertainties.Therefore, it is important to select the optimal µ value for speed control of the PMSM under various operating environments.In this paper, in order to select the optimal µ value, we compare the behavior of the PMSM according to the µ value when the reference speed value of the PMSM, which operates at an initial value of 1200 rpm, is changed to 1500 rpm at 0.5 s.In addition, the condition under which the initial disturbance value changes from 1.0 Nm to 2.0 Nm at 0.7 s is also included.The parameters of the PMSM used in this procedure are listed in Table 1.A total of five µ values were applied, and the results are shown in Figure 2.
An overshoot, undershoot, integral of squared (ISE), and integral of time-weighted squared error (ISET) were used to present the target value tracking characteristics of the PMSM according to µ as a performance index, and the results are shown in Table 2.It can be confirmed through the performance index that as the value of µ increases, the speed control performance for the reference value decreases.Considering the performance of the controller along with the parameter uncertainty, the value of µ is determined to be 1.3, and this was used in the simulations.The process of selecting the value of µ is presented as a tuning procedure of determining µ, as shown in Figure 3.
timal  value, we compare the behavior of the PMSM according to the  value when the reference speed value of the PMSM, which operates at an initial value of 1200 rpm, is changed to 1500 rpm at 0.5 s.In addition, the condition under which the initial disturbance value changes from 1.0 Nm to 2.0 Nm at 0.7 s is also included.The parameters of the PMSM used in this procedure are listed in Table 1.A total of five  values were applied, and the results are shown in Figure 2.An overshoot, undershoot, integral of squared (ISE), and integral of time-weighted squared error (ISET) were used to present the target value tracking characteristics of the PMSM according to  as a performance index, and the results are shown in Table 2.It can be confirmed through the performance index that as the value of  increases, the speed control performance for the reference value decreases.Considering the performance of the controller along with the parameter uncertainty, the value of  is determined to be 1.3, and this was used in the simulations.The process of selecting the value of  is presented as a tuning procedure of determining , as shown in Figure 3.

Simulation Result
To validate the proposed DOB and RMPSC, simulations were carried out under three cases.The simulation model containing the proposed method was built via MAT-LAB/Simulink with the optimal uncertainty µ values determined in Section 5.1.The performance of the previous method and proposed method were compared with the PI decoupling method [10] and sliding-mode method [17], respectively, and the main parameters of each controller, including the proposed RMPSC method, are shown in Table 3.The parameters of the PMSM used in the simulation are listed in Table 1.In the first case, the reference speed ω * e of the PMSM is changed, the disturbance ϵ applied to the PMSM is fixed to a constant value, and the PMSM has nominal parameters.Figure 4 shows the comparison results of the proposed RMPSC method and comparative methods for speed control of PMSM when the reference speed ω * e of the PMSM is increased.The initial state of the PMSM is a rotation speed of 1500 rpm and a disturbance of 1.5 Nm. Figure 4 shows the simulation results for the proposed method and the comparative method when the reference speed of the PMSM is increased to 1800 rpm at 0.6 s.The speed results of the PMSM applying the proposed RMPSC and comparative methods are shown in Figure 4a. Figure 4b-d show the three-phase current of the PMSM, the electromagnetic torque value, and the disturbance ε estimated by the DOB when the proposed RMPSC method is applied.It can be seen in Figure 4a that the speed results applying the proposed RMPSC method reach the reference speed ω * e faster than other comparative methods.

Simulation Result
To validate the proposed DOB and RMPSC, simulations were carried out under three cases.The simulation model containing the proposed method was built via MATLAB/Simulink with the optimal uncertainty  values determined in Section 5.1.The performance of the previous method and proposed method were compared with the PI decoupling method [10] and sliding-mode method [17], respectively, and the main parameters of each controller, including the proposed RMPSC method, are shown in Table 3.The parameters of the PMSM used in the simulation are listed in Table 1.

Controller
Description Value PI decoupling [10] P gain/I gain 7/600 SMC [17] k /k 1000/200 SMC [17] P gain/I gain In the first case, the reference speed  * of the PMSM is changed, the disturbance  applied to the PMSM is fixed to a constant value, and the PMSM has nominal parameters.Figure 4 shows the comparison results of the proposed RMPSC method and comparative methods for speed control of PMSM when the reference speed  * of the PMSM is increased.The initial state of the PMSM is a rotation speed of 1500 rpm and a disturbance of 1.5 Nm. Figure 4 shows the simulation results for the proposed method and the comparative method when the reference speed of the PMSM is increased to 1800 rpm at 0.6 s.The speed results of the PMSM applying the proposed RMPSC and comparative methods are shown in Figure 4a. Figure 4b-d show the three-phase current of the PMSM, the electromagnetic torque value, and the disturbance ϵ estimated by the DOB when the proposed RMPSC method is applied.It can be seen in Figure 4a that the speed results applying the proposed RMPSC method reach the reference speed  * faster than other comparative methods.In the second case, the disturbance  applied to the PMSM is changed, the reference speed  * of the PMSM is kept at a constant value, and the PMSM has nominal parameters.Figure 5 shows the simulation results of the proposed RMPSC and a comparative method for speed control of the PMSM when the disturbance  is increased.Simulations were performed when the disturbance increases to 2.5 Nm in 0.6 s. under the same initial conditions as in the first case.The speed control results, applying the proposed method and the comparative methods, can be seen in Figure 5a.The three-phase current and electromagnetic torque values of the PMSM, and the disturbance ϵ estimated by the DOB are the results of applying the proposed RMPSC method and can be confirmed in Figure 5bd.When an additional disturbance is applied, it can be seen in Figure 5a that the proposed RMPSC method has a small speed dip and quickly recovers to the steady state  * compared to the other comparative methods.
(a) In the second case, the disturbance ϵ applied to the PMSM is changed, the reference speed ω * e of the PMSM is kept at a constant value, and the PMSM has nominal parameters.Figure 5 shows the simulation results of the proposed RMPSC and a comparative method for speed control of the PMSM when the disturbance ϵ is increased.Simulations were performed when the disturbance increases to 2.5 Nm in 0.6 s. under the same initial conditions as in the first case.The speed control results, applying the proposed method and the comparative methods, can be seen in Figure 5a.The three-phase current and electromagnetic torque values of the PMSM, and the disturbance ε estimated by the DOB are the results of applying the proposed RMPSC method and can be confirmed in Figure 5b-d.When an additional disturbance is applied, it can be seen in Figure 5a that the proposed RMPSC method has a small speed dip and quickly recovers to the steady state ω * e compared to the other comparative methods.In the second case, the disturbance  applied to the PMSM is changed, the reference speed  * of the PMSM is kept at a constant value, and the PMSM has nominal parameters.Figure 5 shows the simulation results of the proposed RMPSC and a comparative method for speed control of the PMSM when the disturbance  is increased.Simulations were performed when the disturbance increases to 2.5 Nm in 0.6 s. under the same initial conditions as in the first case.The speed control results, applying the proposed method and the comparative methods, can be seen in Figure 5a.The three-phase current and electromagnetic torque values of the PMSM, and the disturbance ϵ estimated by the DOB are the results of applying the proposed RMPSC method and can be confirmed in Figure 5bd.When an additional disturbance is applied, it can be seen in Figure 5a that the proposed RMPSC method has a small speed dip and quickly recovers to the steady state  * compared to the other comparative methods.
(a) In the last case, a simulation was carried out with parameter uncertainty added to the previous two cases.The initial conditions were the same as the previous cases, the disturbance  and reference speed  * changed at each specific times, and the parameters were set to have a 20% error compared to the nominal parameters.The simulation started from the initial conditions, and the reference speed  * was increased to 1800 rpm at 0.5 s, and the disturbance  was changed to 3.0 Nm at 0.7 s, and the results are shown in Figure 6.The speed control results adapted with the proposed and comparative methods are shown in Figure 6a, and the three-phase current and electromagnetic torque values and the estimated disturbance can be seen in Figure 6b-d, as the results of applying the proposed RMPSC method.Figure 6a shows that the proposed RMPSC method has better speed control performance compared to other comparative methods for reference speed changes and disturbances under parameter uncertainty.In the last case, a simulation was carried out with parameter uncertainty added to the previous two cases.The initial conditions were the same as the previous cases, the disturbance ϵ and reference speed ω * e changed at each specific times, and the parameters were set to have a 20% error compared to the nominal parameters.The simulation started from the initial conditions, and the reference speed ω * e was increased to 1800 rpm at 0.5 s, and the disturbance ϵ was changed to 3.0 Nm at 0.7 s, and the results are shown in Figure 6.The speed control results adapted with the proposed and comparative methods are shown in Figure 6a, and the three-phase current and electromagnetic torque values and the estimated disturbance can be seen in Figure 6b-d, as the results of applying the proposed RMPSC method.Figure 6a shows that the proposed RMPSC method has better speed control performance compared to other comparative methods for reference speed changes and disturbances under parameter uncertainty.To analyze the speed control performance of the proposed RMPSC method, based on simulation results for the previous three cases, the performance index is expressed as shown in Table 4.The performance index confirms that the proposed RMPSC method provides good speed control performance for a PMSM against comparative methods in three simulation cases.To analyze the speed control performance of the proposed RMPSC method, based on simulation results for the previous three cases, the performance index is expressed as shown in Table 4.The performance index confirms that the proposed RMPSC method provides good speed control performance for a PMSM against comparative methods in three simulation cases.

Conclusions
This study proposes a robust model predictive speed control (RMPSC) method for the speed control of a PMSM that is robust to disturbance and parameter uncertainty.External mechanical torque, which affects the PMSM, was considered a disturbance, and this value was estimated through a disturbance observer.In addition, uncertainty in parameters was considered to ensure stable speed control performance of the PMSM under manufacturing error or aging degradation affecting the PMSM conditions.The disturbance and state observer were designed to satisfy Lyapunov stability and a cost function was constructed based on the steady state of the PMSM at the reference speed.Both the observer and controller design method were transformed into LMI-based optimization problems, and a tuning procedure of the disturbance observer, state observer, and RMPSC controller was proposed so that optimal gains could be obtained systematically.To verify the performance of the proposed methods, simulations were carried out under predefined test conditions.
The simulation results showed the effectiveness of the proposed RMPSC method in providing speed control performance to the PMSM even under conditions where the target speed changes and external disturbances are applied to the PMSM.In particular, the fast response characteristics and speed control performance were ensured by quickly estimating external disturbances through a disturbance observer, reflecting them in the RMPSC controller, and controlling the speed of the PMSM through state feedback.Furthermore, considering the uncertainty of parameters, speed control performance for the PMSM was secured even in an environment where parameters are not nominal values and optimal control performance was achieved through the proposed optimization-based systematic tuning procedure to find the optimal observer and controller gains.

Figure 1 .
Figure 1.Block diagram of the proposed RMPSC method with a disturbance observer.

Figure 1 .
Figure 1.Block diagram of the proposed RMPSC method with a disturbance observer.

Figure 2 .
Figure 2. Simulation results of the dynamic response of the PMSM when changing the reference speed value ω * e and disturbance ϵ according to µ values: (a) Speed results of the PMSM according to the proposed RMPSC, (b) measured disturbance ϵ and estimated disturbance ε using the DOB.

Figure 2 .
Figure 2. Simulation results of the dynamic response of the PMSM when changing the reference speed value  * and disturbance  according to  values: (a) Speed results of the PMSM according to the proposed RMPSC, (b) measured disturbance  and estimated disturbance  using the DOB.

Figure 3 .
Figure 3. Tuning procedure to obtain the optimal gains by determining .Figure 3. Tuning procedure to obtain the optimal gains by determining µ.

Figure 3 .
Figure 3. Tuning procedure to obtain the optimal gains by determining .Figure 3. Tuning procedure to obtain the optimal gains by determining µ.

Figure 4 .
Figure 4. Simulation results of the dynamic response of the PMSM when changing the reference speed value  * : (a) Speed results of the PMSM applying the PI-decoupling method, sliding mode control method, and proposed RMPSC method, respectively, (b) the three-phase current of the PMSM under the proposed RMPSC method, (c) the electromagnetic torque  of the PMSM under the proposed RMPSC method, (d) the measured disturbance  and estimated disturbance  according to the DOB.

Figure 4 .
Figure 4. Simulation results of the dynamic response of the PMSM when changing the reference speed value ω * e : (a) Speed results of the PMSM applying the PI-decoupling method, sliding mode control method, and proposed RMPSC method, respectively, (b) the three-phase current of the PMSM under the proposed RMPSC method, (c) the electromagnetic torque τ e of the PMSM under the proposed RMPSC method, (d) the measured disturbance ϵ and estimated disturbance ε according to the DOB.

Figure 4 .
Figure 4. Simulation results of the dynamic response of the PMSM when changing the reference speed value  * : (a) Speed results of the PMSM applying the PI-decoupling method, sliding mode control method, and proposed RMPSC method, respectively, (b) the three-phase current of the PMSM under the proposed RMPSC method, (c) the electromagnetic torque  of the PMSM under the proposed RMPSC method, (d) the measured disturbance  and estimated disturbance  according to the DOB.

Figure 5 .
Figure 5. Simulation results of the dynamic response of the PMSM when disturbance ϵ changes in the PMSM: (a) speed results of the PMSM applying the PI-decoupling method, sliding-mode control method, and proposed RMPSC method, respectively, (b) the three-phase current of the PMSM under proposed the RMPSC method, (c) the electromagnetic torque τ e of the PMSM under the proposed RMPSC method, (d) the measured disturbance ϵ and estimated disturbance ε according to the DOB.

Figure 6 .
Figure 6.Simulation results of the dynamic response of the PMSM when disturbance  changes in the PMSM with parameter uncertainty: (a) speed results of the PMSM applying the PI-decoupling method, sliding-mode control method, and the proposed RMPSC method, respectively, (b) the three-phase current of the PMSM under the proposed RMPSC method, (c) the electromagnetic torque  of the PMSM under the proposed RMPSC method, (d) the measured disturbance  and estimated disturbance  according to the DOB.

Figure 6 .
Figure 6.Simulation results of the dynamic response of the PMSM when disturbance ϵ changes in the PMSM with parameter uncertainty: (a) speed results of the PMSM applying the PI-decoupling method, sliding-mode control method, and the proposed RMPSC method, respectively, (b) the threephase current of the PMSM under the proposed RMPSC method, (c) the electromagnetic torque τ e of the PMSM under the proposed RMPSC method, (d) the measured disturbance ϵ and estimated disturbance ε according to the DOB.

3 2
Pλ r and disturbance ϵ  and  are the rotor position, viscosity coefficient, and m and  is the unknown external torque, respectively.Since the mechanic known external torque cannot be measured in the system, these values a turbance .
and  are the rotor position, viscosity co and  is the unknown external torque, respectively.Sin known external torque cannot be measured in the system turbance .
and  are the rotor position and  is the unknown external torque, resp known external torque cannot be measured turbance .
and  are the rotor position, viscosity coefficient, and m and  is the unknown external torque, respectively.Since the mechani known external torque cannot be measured in the system, these values turbance .

Table 2 .
Tuning result of the proposed RMPSC method: PMSM speed control.

Table 2 .
Tuning result of the proposed RMPSC method: PMSM speed control.

Table 3 .
Main parameters of each controller.

Table 3 .
Main parameters of each controller.

Table 4 .
Performance index comparison results for PMSM speed control.