Next Article in Journal
Dilated Filters for Edge-Detection Algorithms
Next Article in Special Issue
Application of a Trajectory Tracking Algorithm for Underactuated Underwater Vehicles Using Quasi-Velocities
Previous Article in Journal
A Validation Method for EPID In Vivo Dosimetry Algorithms
Previous Article in Special Issue
A Mechanical Feedback Classification of Linear Mechanical Control Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control

by
Sławomir Stępień
*,† and
Paulina Superczyńska
*,†
Institute of Automatic Control and Robotics, Poznań University of Technology, Piotrowo 3a, 60-965 Poznań, Poland
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2021, 11(22), 10714; https://doi.org/10.3390/app112210714
Submission received: 24 October 2021 / Revised: 10 November 2021 / Accepted: 11 November 2021 / Published: 13 November 2021

Abstract

:
This paper presents modeling and infinite-time suboptimal control of a quadcopter device using the state-dependent Riccati equation (SDRE) method. It establishes a solution to the control problem using SDRE and proposes a new procedure for solving the problem. As a new contribution, the paper proposes a modified SDRE-based suboptimal control technique for affine nonlinear systems. The method uses a pseudolinearization of the closed-loop system employing Moore–Penrose pseudoinverse. Then, the algebraic Riccati equation (ARE), related to the feedback compensator gain, is reduced to state-independent form, and the solution can be computed only once in the whole control process. The ARE equation is applied to the problem reported in this study that provides general formulation and stability analysis. The effectiveness of the proposed control technique is demonstrated through the use of simulation results for a quadrotor device.

1. Introduction

The state-dependent Riccati equation (SDRE) method has become more popular due to the possibility of designing more effective nonlinear controllers and providing flexibility through state-dependent weighting matrices, which is well presented in [1]. This method scheme originates from the LQR method [2,3] and is also similar to the MPC (or its nonlinear version, NMPC) [4,5,6]. The main differences between SDRE and those two methods are, for example, unlike the LQR method, the SDRE sets up its control laws at every time-step. The MPC method requires a convex optimization technique, while SDRE and LQR methods utilize optimal control laws for linear systems [7,8].
The SDRE control method reformulates the nonlinear dynamics using parameterization to a linear structure with state-dependent coefficient (SDC) matrices, which are essential to solve an algebraic Riccati equation (ARE) and obtain the suboptimal control law. The coefficients of the ARE equation differ in the state space. Thus, at a given point in state space, it is essential to solve an algebraic state-dependent Riccati equation. The nonuniqueness of the state-dependent parameterization (SDP) allows us to obtain extra degrees of freedom, which may help to enhance controller performance. Minimization of the nonlinear performance index gives a quadratic-like structure [9,10,11,12].
In the last decade, this method has been used in many different areas of research, for example, satellite and spacecraft control and estimation [13,14], control of servopneumatic drives [15], control of quadrotor aircraft [16,17,18,19,20], robotics [21,22], control of non-affine systems [23], and even in medical sciences for cancer treatment allowing optimization of drug dosing control [24].
The main disadvantage of this classic SDRE strategy is that Riccati equations are state dependent. Therefore, it is necessary to solve SDRE multiple times during the control process [25,26]. Nowadays, the solution can be successfully implemented in real control systems due to developments of digital controllers. However, the solution can be found faster, reducing computational effort and introducing linearization of closed-loop form. As a new contribution, the paper presents a possible solution for suboptimal control of nonlinear systems by solving the algebraic Riccati equation only once in the whole control process. Using linearization of the closed-loop system and introducing Moore–Penrose inversion [27], it is possible to modify SDRE into ARE. For the modified SDRE method, it is possible to obtain the same solution as for classic SDRE method. The advantage of the modified case is that there is no need to solve the state-dependent Riccati equation for each state of the control process. This allows us to save space for implementation in real-time control systems and simplifies computations. The specific differences between classical and modified methods are presented in the sections below.

2. SDRE Method

Let Equation (1) be the continuous time nonlinear system
x ˙ = F ( x ) + B u ,
where x R n and u R m are the state and input vectors, respectively. F ( x ) is a continuous ( C k class) nonlinear function of x , and B R n × m is the constant control matrix. The problem consists of finding an admissible control u R m for t 0 t , which minimizes the performance index [28,29,30,31]
J ( u ) = 1 2 0 ( x T Q x + u T R u ) d t ,
where Q R n × n is a symmetric positive semidefinite, and R R m × m is a symmetric positive-definite matrix. Rewriting the nonlinear Equation (1), the state-dependent coefficient form (SDC) follows
x ˙ = A ( x ) x + B u ,
where A ( x ) R n × n .
There are many methods for SDC parameterization. Each parameterization should be true for all 0 α 1 [32]
α A I ( x ) x + ( 1 α ) A I I ( x ) x = α F ( x ) + ( 1 α ) F ( x ) = F ( x ) .
The control problem for the nonlinear continuous-time system Equation (1) can be formulated as follows. Given nonlinear functions F ( x ) R n , B R n × m , Q R n × n , and R R m × m , we attempt to find a control vector u R m for t [ t 0 , ] that allows us to control the system state vector from x 0 to x while the performance index is minimized (Equation (19)).
If the integral functions x T Q x + u T R u and A ( x ) x + B u are continuously differentiable functions of each of their arguments, then we may imply that u C [ t 0 , ] is a control that minimizes the functional J ( u ) : C [ t 0 , ] R . To solve the problem, consider the Hamiltonian
H = 1 2 ( x T Q x + u T R u ) + p T ( A ( x ) x + B u ) ,
where p R n is the co-state vector.
To be consistent with the Pontryagin minimum principle, if u C [ t 0 , ] is a control for the functional Equation (19), subject to state Equation (3), and if x denotes the corresponding state, then there exists a p C [ t 0 , ] such that
H u ( p , x , u , t ) = 0 f o r t [ t 0 , ]
and
p ˙ = H x ( p , x , u , t ) f o r t [ t 0 , ] a n d p ( ) = 0 .
Equations (6) and (7) are conditions allowing us to obtain suboptimal control while minimizing Equation (19).
It follows that any optimal input u R m and the corresponding state x R n satisfies Equation (6); that is,
H u = R u + B T p = 0 .
From Equation (8), the optimal control vector is
u = R 1 B T p ,
and subsequently, the adjoint differential equation is
p ˙ = H x = ( A ( x ) ) x x T p Q x
for t [ t 0 , ] , x ( t 0 ) = x 0 , and p ( ) = 0 , where
( A ( x ) x ) x = A ( x ) + A ( x ) x x .
Let p be a nonlinear combination of the state of the system
p = K ( x ) x ,
where K ( x ) R n × n denotes feedback compensator gain, and let x be the solution of the nonlinear state equation
x ˙ = A ( x ) x B R 1 B T K ( x ) x
for t [ t 0 , ] , x ( t 0 ) = x 0 . Then, feedback control is
u = R 1 B T K ( x ) x .
Let K ( x ) solve the state-dependent Riccati equation (SDRE), which is
K ( x ) A ( x ) + A T ( x ) K ( x ) K ( x ) B R 1 ( x ) B T K ( x ) + Q ( x ) = 0 ;
then the equation below satisfies the optimality condition [25,32]
K ˙ ( x ) + ( A ( x ) x ) x T K ( x ) = 0 .
The suboptimal control
u = R 1 B T K ( x ) x
for performance index (objective function) Equation (19) subject to dynamics Equation (1) can be found solving the state-dependent Riccati Equation (15). Often, it is difficult to find the solution of Equation (15) analytically. One approach allows us to solve the SDRE via symbolic software packages such as Matlab. However, for complex systems, the solution may become complicated, and then, it is necessary to approximate the solution. To approximate, an interpolation method or Taylor series method can be used, for instance, ref. [25].

3. Modified SDRE Method

Let Equation (18) be the continuous time nonlinear system
x ˙ = F ( x ) + B u ,
where x R n and u R m are the state and input vectors. F ( x ) is a continuous ( C k class) nonlinear function of x , and B R n × m is the constant control matrix. The problem consists of finding an admissible control u R m for t 0 t that minimizes the performance index [28,29,30,31]
J ( u ) = 1 2 0 ( x T Q x + u T R u ) d t ,
where Q R n × n is a symmetric positive semidefinite, and R R m × m is a symmetric positive-definite matrix.
The first modification (compared to the classic approach) is to rewrite the system Equation (18) as a sum of state-independent and nonlinear state-dependent coefficient (SDC) forms Equation [25]:
x ˙ = A ( x ) x + B u = A 1 x + A 2 ( x ) x + B u ,
where
A ( x ) = A 1 + A 2 ( x )
is a sum of a constant matrix and a state-dependent matrix [12,33].
The SDC parameterization of Equation (20) can be performed the same way as for Equation (4), for all 0 α 1 [32]:
α A 1 , I + A 2 , I ( x ) x + ( 1 α ) A 1 , I I + A 2 , I I ( x ) x = α A 1 , I + ( 1 α ) A 1 , I I + α A 2 , I ( x ) + ( 1 α ) A 2 , I I ( x ) x = α F I ( x ) + ( 1 α ) F I I ( x ) = F ( x ) .
To easily formulate controllability criteria, the state-dependent controllability matrix is introduced W ( x ) [25]
W ( x ) = [ B ( A 1 + A 2 ( x ) ) B ( A 1 + A 2 ( x ) ) n 1 B ] .
If W ( x ) (state-dependent) has the full rank, then the system is controllable for all x R n [34]. As in the previous case of the SDRE method, an admissible control u ( t ) R m , that minimizes the performance index can be found for [28,29,30,31]
J ( u ) = 1 2 0 ( x T Q x + u T R u ) d t ,
where Q R n × n is a symmetric positive semi- definite matrix, and R R m × m is a symmetric positive definite matrix. The control problem for the nonlinear continuous-time system Equation (18) can be formulated as follows. Given nonlinear functions F ( x ) R n , B R n × m , Q R n × n and R R m × m , we attempt to find a control vector u R m for t [ t 0 , ] that controls the system state vector from x 0 to x while minimizing the performance index (Equation (19)).
If the integrand of the cost functions x T Q x + u T R u and A ( x ) x + B u are continuously differentiable functions of each of their arguments, then we may imply that u C [ t 0 , ] is a control that minimizes the functional J ( u ) : C [ t 0 , ] R . Note that, in this approach, two feedback compensators are defined A ( x ) = A 1 + A 2 ( x ) . To solve the problem, consider the Hamiltonian
H = 1 2 ( x T Q x + u T R u ) + p T ( ( A 1 + A 2 ( x ) ) x + B u ) ,
where p R n is the co-state vector.
If u C [ t 0 , ] is a control that minimizes Equation (19) subject to the state Equation (20), and if x denotes the corresponding state, then there exists a p C [ t 0 , ] such that
H u ( p , x , u , t ) = 0 f o r t [ t 0 , ]
and
p ˙ = H x ( p , x , u , t ) f o r t [ t 0 , ] a n d p ( ) = 0 ,
where p satisfies the condition for optimal control that minimizes Equation (19).
It follows that any optimal input u R m and the corresponding state x R n satisfy Equation (26); that is,
H u = R u + B T p = 0 .
Thus, the optimal control is
u = R 1 B T p ,
and subsequently, the adjoint differential equation is
p ˙ = H x = A 1 + ( A 2 ( x ) ) x x T p Q x ,
for t [ t 0 , ] , x ( t 0 ) = x 0 , A 2 ( x 0 ) = A 2 0 , and p ( ) = 0 , where
( A 2 ( x ) x ) x = A 2 ( x ) + A 2 ( x ) x x .
Consider p as a linear and nonlinear combination of the state of the system Equation (18), which is the second modification of the classic approach
p = K ( x ) x = K 1 x + K 2 ( x ) x ,
where K ( x ) , K 1 , K 2 ( x ) R n × n .
Consider x as the solution of the nonlinear state equation
x ˙ = A 1 x + A 2 ( x ) x B R 1 B T K 1 x B R 1 B T K 2 ( x ) x
for t [ t 0 , ] , x ( t 0 ) = x 0 , and A 2 ( x 0 ) = A 20 .
Note that, in this approach, two feedback compensators are defined.
Rearranging Equation (33), it is possible to obtain
x ˙ = ( A 1 B R 1 B T K 1 ) x + ( A 2 ( x ) B R 1 B T K 2 ( x ) ) x .
The first bracket of Equation (34) is state-independent, and the second one is state-dependent; thus, it is possible to linearize it and solve the state-dependent gain matrix- K 2 ( x ) from:
A 2 ( x ) B R 1 B T K 2 ( x ) = 0
as follows:
K 2 ( x ) = B R 1 B T + A 2 ( x ) .
Matrix B R 1 B T is singular; thus, the state-dependent matrix gain K 2 ( x ) may be computed only by the pseudoinverse operation + .
In order to perform this procedure, a Moore–Penrose pseudoinverse is used, and the matrix is unique [27,35,36].
Equating the adjoint differential Equation (30) and the differential form of Equation (32), we have
t [ K 1 + K 2 ( x ) ] x + [ K 1 + K 2 ( x ) ] x ˙ = A 1 + A 2 ( x ) + ( A 2 ( x ) ) x x T [ K 1 + K 2 ( x ) ] x Q x .
Consequently, employing Equations (22) and (29), it is possible to obtain
t [ K 1 + K 2 ( x ) ] x + [ K 1 + K 2 ( x ) ] A 1 x + [ K 1 + K 2 ( x ) ] A 2 ( x ) x [ K 1 + K 2 ( x ) ] B R 1 B T [ K 1 + K 2 ( x ) ] x = A 1 + A 2 ( x ) + ( A 2 ( x ) ) x x T [ K 1 + K 2 ( x ) ] x Q x ,
and rearranging terms in Equation (39), we find
K ˙ 2 ( x ) + ( A 2 ( x ) ) x x T [ K 1 + K 2 ( x ) ] x + [ K 1 + K 2 ( x ) ] A 1 + [ K 1 + K 2 ( x ) ] A 2 ( x ) [ K 1 + K 2 ( x ) ] B R 1 B T [ K 1 + K 2 ( x ) ] + [ A 1 + A 2 ( x ) ] T [ K 1 + K 2 ( x ) ] + Q x = 0 .
Based on linearization Equations (34)–(36), let the matrix K 1 solve Equation (39), the algebraic Riccati equation (ARE):
K 1 A 1 + A 1 T K 1 K 1 B R 1 B T K 1 + Q = 0 ;
then, the optimality condition is given by
K ˙ 2 ( x ) + A 2 ( x ) x x T K 1 + A 2 ( x ) x x T K 2 ( x ) + K 2 ( x ) A 1 + A 1 T K 2 ( x ) + K 1 A 2 ( x ) + A 2 T ( x ) K 1 + K 2 ( x ) A 2 ( x ) + A 2 T ( x ) K 2 ( x ) K 2 ( x ) B R 1 B T K 1 K 1 B R 1 B T K 2 ( x ) K 2 ( x ) B R 1 B T K 2 ( x ) = 0
for K 2 ( x ) obtained from (36). So, the suboptimal control
u = R 1 B T K 1 + K 2 ( x ) x .
For index Equation (19), subject to Equation (18), it is possible to find the solution for the state-independent Riccati Equation (39). In general, it is difficult to find the solution of Equation (39) analytically, especially for complex systems [22,25,33,37]. Symbolic software packages (such Matlab or Mathematica) are commonly used to solve such problems. However, for complex systems, the solution may become complicated, and then, it is necessary to approximate the solution. To approximate, an interpolation method or Taylor series method seems to be most suitable [15].

4. Stability Analysis

Asymptotic stability of the closed-loop system ensures the possibility of controlling the states from the initial values to the final ones. The controlled system with the SDRE compensator-based feedback is locally asymptotically stable [25]. Let r > 0 be the largest radius, such that B r ( 0 ) Ω . Assuming that the system is stabilizable at x = 0 , it is possible to use optimal control theory to define matrix K 1 , such that all eigenvalues of ( A 1 B K 1 ) have negative real parts. There exist β > 0 , such that R e ( p ) < β for all eigenvalues p of ( A 1 B K 1 ) , having the system
x ˙ = ( A 1 B R 1 B T K 1 ) x + ( A 2 ( x ) B R 1 B T K 2 ( x ) ) x ,
where ( A 1 + A 2 ( x ) ) x and ( A 1 x + A 2 ( x ) x ) x are continuous in x for | | x | | < r , where r > 0 is the largest radius around the origin x = 0 .
Let
g ( x ) = A 2 ( x ) B K 2 ( x ) ,
and h ( x ) = g ( x ) x ; then, the system (34) is described by
x ˙ = ( A 1 B K 1 ) x + h ( x ) .
Evaluation of h ( x ) shows that it is an almost linear system:
lim x 0 h ( x ) x = 0 ,
from the inequality
h ( x ) g ( x ) x = A 2 ( x ) B K 2 ( x ) x .
Thus, g ( x ) 0 when x 0 ; then, h ( x ) satisfies the condition Equation (46). The theoretical considerations of almost linear systems leads to the conclusion that, if eigenvalues ( A 1 B K 1 ) have negative real parts, h ( x ) is continuous around the origin, and the condition of Equation (46) holds, then x = 0 is asymptotically stable [38].
Let δ > 0 ; then, there exists a η ^ ( 0 , r ) , such that h ( x ) δ x whenever x η ^ . Let x ( t 0 ) = x 0 B η ^ ( 0 ) ; assuming that f is continuous, the solution exists and remains in B η ^ ( 0 ) for t [ t 0 , ] , then
x ( t ) = e ( A 1 B K 1 ) t x ( t 0 ) + t 0 t e ( A 1 B K 1 ) ( t s ) h ( x ( s ) ) d s .
Considering the norm and the assumption | | x ( t ) | | < η ^ , it is possible to obtain
x ( t ) e ( A 1 B K 1 ) t x ( t 0 ) + δ t 0 t e ( A 1 B K 1 ) ( t s ) x ( s ) d s .
There exists a positive constant G and β , such that
e ( A 1 B K 1 ) t G e β t
and
x ( t ) G e β t x ( t 0 ) + δ G t 0 t e β ( t s ) x ( s ) d s .
Referring to the Gronwall inequality and after multiplication by e β t , we obtain
x ( t ) G x ( t 0 ) e ( β δ G ) t
for all t > 0 , such that x ( t ) < η ^ . δ is chosen from δ ( 0 , β G ) ; then, with respect to η ^ having γ ( 0 , η ^ ) , it is possible to obtain η = min { η ^ , γ G } . For all x 0 B η , it is possible to obtain x ( t ) < γ η ^ for all t > 0 , and x = 0 is stable.
The condition of Equation (52) is satisfied for all t > 0 if x 0 B η , x = 0 is asymptotically stable, since β δ G > 0 [25].
For the proposed method, a global uniform stability can also be proved using the Lyapunov function. The proof can be obtained more easily than for classic SDRE, because feedback compensator gains are state-independent, and for a closed-loop system, there is no problem with generalization and definition of an attraction region.
Assuming that SDC parameterization is stabilizable and detectable for all x ( t ) , then the closed-loop matrix A C L is symmetric for all x ( t ) , and the solution of SDRE is asymptotically stable when t [ t 0 , ] . Let the control law be in the form
u = R 1 B T ( K 1 + K 2 ( x ) ) x ,
and K 1 satisfies ARE
A 1 T K 1 + K 1 A 1 K 1 B R 1 B T K 1 + Q = 0 ,
while K 2 ( x ) is taken from
K 2 ( x ) = [ B R 1 B T ] + A 2 ( x ) .
Then, the system in closed-loop form
x ˙ = ( A 1 B R 1 B T K 1 ) x + ( A 2 ( x ) B R 1 B T K 2 ( x ) ) x
is globally asymptotically stable if K 1 is the unique symmetric positive-defined matrix.
Let us define the as Lyapunov function
V ( x ) = x T K 1 x
and
k 1 x 2 V ( x ) k 2 x 2
for x U , where U in R n has the origin, and constants k 1 and k 2 are obtained from
k 1 = inf σ i ( K 1 )
and
k 2 = sup σ i ( K 1 )
for i = 1 , 2 , . . . , n .
Using Equations (54) and (56), it is possible to obtain
V ˙ ( x ) = x T K 1 x ˙ + x ˙ T K 1 x .
Taking Equation (56), assuming that K 2 ( x ) assures linearization of Equations (45) and (46), and that the solution of ARE, Equation (54), is always symmetric, then
V ˙ ( x ) = x T A 1 T K 1 + K 1 A 1 2 K 1 B R 1 B T K 1 x .
Then, taking into consideration (54), it is possible to obtain
V ˙ ( x ) = x T K 1 B R 1 B T K 1 Q x .
Let us assume that
V ˙ ( x ) k 3 x 2 ,
where
k 3 = inf σ i ( K 1 B R 1 B T K 1 + Q )
for x U and i = 1 , 2 , . . . , n .
Generally, in an SDRE control problem with infinite time horizon, the existence of constants k 1 , k 2 , k 3 do not assure global stability. The closed-loop matrix is state-dependent but can be generalized to global stability by defining the region of attraction [39].
In the proposed method, the existence of constants k 1 , k 2 , k 3 proposed in Equations (59), (60), and (65) assures global stability, because the matrix K 1 is state-independent, so there is no reason to define the region of attraction for the Riccati Equation (54).

5. Numerical Simulation

To present a numerical example, the mentioned quadcopter model was chosen. Typical SDRE control technique was compared with its modified version.
Considering a quadcopter (presented in Figure 1) as a vehicle with four independent drives and with an electric power system located at its center of gravity, it is important to define coordinate systems in order to determine the position. A model device has six degrees of freedom, where vertical movement in the global coordinate system and three Euler angles are controlled parameters. There is a need to assume a body fixed frame and a symmetrical structure of the model with the origin in the center of mass, where every drive is independently controlled and aerodynamic effects are negligible. The mathematical model is taken from [26].
As mentioned, this model needs to be described using Euler’s angles: roll angle ϕ π , π , pitch angle θ π 2 , π 2 , and yaw angle ψ π , π . These angles are a sequence of three rotations
R x ( ϕ ) = 1 0 0 0 cos ( ϕ ) sin ( ϕ ) 0 sin ( ϕ ) cos ( ϕ ) , R y ( θ ) = cos ( θ ) 0 sin ( θ ) 0 1 0 sin ( θ ) 0 cos ( θ ) , R z ( ψ ) = cos ( ψ ) sin ( ψ ) 0 sin ( ψ ) cos ( ψ ) 0 0 0 1 .
This leads to the rotation matrix R , which describes rotations of the local system in terms of the global system.
R z y x ( ψ , θ , ϕ ) = R z ( ψ ) R y ( θ ) R x ( ϕ )
The representation of the actual model in the simulations requires us to determine differential equations that describe all relations in the model. The orientation of the vehicle is described by two vectors, r T = ( x , y , z ) and Ω T = ( ψ , ϕ , θ ) . The quadrotor’s acceleration is described by the differential equation
r ¨ = g 0 0 1 + R b m i = 1 4 ω i 2 0 0 1 ,
where:
g—the gravitational acceleration;
R —the rotation matrix;
b—the thrust factor;
ω i —the speed of rotor i = 1 , 2 , 3 , 4 .
The second differential equation is presented in
I Ω ¨ = Ω ˙ × I Ω ˙ i = 1 4 J R Ω ˙ × 0 0 1 ω i + τ .
Vector τ is defined as
τ = l b ( ω 4 2 ω 2 2 ) l b ( ω 3 2 ω 1 2 ) d ( ω 2 2 + ω 4 2 ω 1 2 ω 3 2 ) ,
where:
l—the length the of quadrotor’s arm;
d—the drag factor.
The quadrotor’s model has four control inputs (vertical movement in the global coordinate system and three Euler angles). It is assumed that the input values, which represent forces, are directly proportional to the squared angular velocity of the motor.
u 1 = b ( ω 1 2 + ω 2 2 + ω 3 2 + ω 4 2 ) u 2 = b ( ω 4 2 ω 2 2 ) u 3 = b ( ω 3 2 ω 1 2 ) u 4 = d ( ω 2 2 + ω 4 2 ω 1 2 ω 3 2 ) .
Let us create a new variable that describes relative motor speed
ω d = ω 2 + ω 4 ω 1 ω 3 .
Using Equations (68)–(72), a dynamic model is presented as a system of six equations
x ¨ = ( cos ( ϕ ) sin ( θ ) cos ( ψ ) + sin ( ϕ ) sin ( ψ ) ) u 1 m y ¨ = ( cos ( ϕ ) sin ( θ ) sin ( ψ ) sin ( ϕ ) cos ( ψ ) ) u 1 m z ¨ = g + ( cos ( ϕ ) cos ( θ ) ) u 1 m ϕ ¨ = θ ˙ ψ ˙ ( I y I z ) I x J R I x θ ˙ ω d + l I x u 2 θ ¨ = ϕ ˙ ψ ˙ ( I z I x ) I y + J R I y ϕ ˙ ω d + l I y u 3 ψ ¨ = ϕ ˙ θ ˙ ( I x I y ) I z + l I z u 4 .
It is possible to describe the model in the state-space form, where
u = u 1 u 2 u 3 u 4 T
is the control vector, and x R 12 is the state vector
x = x x ˙ y y ˙ z z ˙ ϕ ϕ ˙ θ θ ˙ ψ ψ ˙ T = x 1 x 2 x 12 T .
Combining Equations (73) and (75), we obtain
x ˙ = x 2 ( cos x 7 sin x 9 cos x 11 + sin x 7 sin x 11 ) u 1 m x 4 ( cos x 7 sin x 9 sin x 11 sin x 7 cos x 11 ) u 1 m x 6 g + ( cos x 7 cos x 9 ) u 1 m x 8 x 12 x 10 I 1 J R I x x 10 ω d + l I x u 2 x 10 x 12 x 8 I 2 + J R I y x 8 ω d + l I y u 3 x 12 x 10 x 8 I 3 + l I z u 4 ,
where I 1 = ( I y I z ) I x , I 2 = ( I z I x ) I y , and I 3 = ( I x I y ) I z .
Considering the attitude of the control, and applying SDRE method, the state vector
x I = x 7 x 8 x 9 x 10 x 11 x 12 T ,
and the control vector
u I = u 2 u 3 u 4 T .
After factorization, it is possible to obtain
x ˙ I = 0 1 0 0 0 0 0 0 0 x 12 I 1 0 0 0 0 0 1 0 0 0 x 12 I 2 0 0 0 0 0 0 0 0 0 1 0 0 0 x 8 I 3 0 0 x I + 0 0 0 l I x 0 0 0 0 0 0 l I y 0 0 0 0 0 0 l I z u I .
The above state-space model has SDC structure, where A ( x I ) is state-dependent matrix, and B is constant, as assumed in Equation (18).
Applying the proposed methodology and replacing A by a sum of two matrices as presented in Equation (20) (for proposed method), it is possible to obtain
x ˙ I = 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 x I + 0 0 0 0 0 0 0 0 0 x 12 I 1 0 0 0 0 0 0 0 0 0 x 12 I 2 0 0 0 0 0 0 0 0 0 0 0 0 0 x 8 I 3 0 0 x I + 0 0 0 l I x 0 0 0 0 0 0 l I y 0 0 0 0 0 0 l I z u I .
To present the influence on the system, weight matrices were changed during this experiment. There were three different configurations:
Q = 1 · I 6 × 6 , R = 100 · I 3 × 3 , Q = 1 · I 6 × 6 , R = 1 · I 3 × 3 , Q = 100 · I 6 × 6 , R = 1 · I 3 × 3 .
The values of parameters used in modeling of the object are presented below in Table 1 [40].
Initial state and velocities were set as
x 0 = π / 4 π / 2 3 π / 4 0 0 0 rad rad rad rad / s rad / s rad / s .
Simulation time was 5 s and time step was 1 × 10 3 s .
The aim of the control process is stabilization at “zero” point.
The comparison of the execution time for both methods (Table 2) shows that the proposed method reduces computational effort, which can be seen in the shortened time of computation below. The computations were executed on Windows 7 Professional, 64-bit Intel® Core™ i5-5200M CPU @ 2.40GHz.
Simulation results for configuration R 1 , Q 100 are presented. Differences between classical and modified SDRE method for next two configurations are unnoticeable.
In every case of Q and R matrices configuration after a short transition phase, all angles are stabilized at a required point (Figure 2, Figure 3 and Figure 4). The waveforms for both methods are convergent. Changing the value of Q and R matrices affects only the speed of regulation. Small differences are due to numerical errors and pseudoinversion inaccuracy.

6. Conclusions

The modeling and control of a quadcopter using the SDRE method and modified SDRE-based proposition were presented in this paper. A new method for solving a suboptimal control problem is established, and a procedure for solving the problem was presented and illustrated with a numerical example. As shown, the modified SDRE method gives the same solution as the classic SDRE method, but in the modified case, there is no need to solve a state-dependent Riccati equation for each state of the control process. This allows us to save space for implementation in real-time control systems and simplifies computations.

Author Contributions

Writing—original draft preparation, S.S. and P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to pay our gratitude and our respects to our long term scientific supervisor and colleague Krzysztof Kozłowski who passed away in 2021.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cloutier, J.R. State-dependent Riccati equation techniques: An overview. In Proceedings of the 1997 American Control Conference (Cat. No. 97CH36041), Albuquerque, NM, USA, 6 June 1997; Volume 2, pp. 932–936. [Google Scholar]
  2. Cohen, M.R.; Abdulrahim, K.; Forbes, J.R. Finite-Horizon LQR Control of Quadrotors on SE_2(3). IEEE Robot. Autom. Lett. 2020, 5, 5748–5755. [Google Scholar] [CrossRef]
  3. Okyere, E.; Bousbaine, A.; Poyi, G.T.; Joseph, A.K.; Andrade, J.M. LQR controller design for quad-rotor helicopters. J. Eng. 2019, 2019, 4003–4007. [Google Scholar] [CrossRef]
  4. Gros, S.; Quirynen, R.; Diehl, M. Aircraft control based on fast non-linear MPC & multiple-shooting. In Proceedings of the 2012 IEEE 51st IEEE Conference on Decision and Control (CDC), Maui, HI, USA, 10–13 December 2012; pp. 1142–1147. [Google Scholar]
  5. Richter, S.; Jones, C.N.; Morari, M. Real-time input-constrained MPC using fast gradient methods. In Proceedings of the 48h IEEE Conference on Decision and Control (CDC) Held Jointly with 2009 28th Chinese Control Conference, Shanghai, China, 15–18 December 2009; pp. 7387–7393. [Google Scholar]
  6. Fnadi, M.; Alexandre dit Sandretto, J. Experimental Validation of a Guaranteed Nonlinear Model Predictive Control. Algorithms 2021, 14, 248. [Google Scholar] [CrossRef]
  7. Lin, L.G. Nonlinear Control Systems: A ‘State-Dependent (Differential) Riccati Equation’ Approach. Ph.D. Thesis, Faculty of Engineering Science, KU Lueven and NCTU, Leuven, Belgium, 2014. [Google Scholar]
  8. Fnadi, M.; Du, W.; Plumet, F.; Benamar, F. Constrained Model Predictive Control for dynamic path tracking of a bi-steerable rover on slippery grounds. Control Eng. Pract. 2021, 107, 104693. [Google Scholar] [CrossRef]
  9. Mracek, C.P.; Cloutier, J.R. Control designs for the nonlinear benchmark problem via the state-dependent Riccati equation method. Int. J. Robust Nonlinear Control 1998, 8, 401–433. [Google Scholar] [CrossRef]
  10. Çimen, T. State-dependent Riccati equation (SDRE) control: A survey. IFAC Proc. Vol. 2008, 41, 3761–3775. [Google Scholar] [CrossRef] [Green Version]
  11. Kaczorek, T. Minimum energy control of positive continuous- time linear systems with bounded inputs. Int. J. Appl. Math. Comput. Sci. 2013, 23, 725–730. [Google Scholar] [CrossRef] [Green Version]
  12. Bernat, J.; Stepien, S.J.; Stranz, A.; Superczynska, P. Infinite-time linear–quadratic optimal control of the BLDC motor exploiting a nonlinear finite element model. COMPEL- Int. J. Comput. Math. Electr. Electron. Eng. 2017, 36, 633–648. [Google Scholar] [CrossRef]
  13. Stansbery, D.T.; Cloutier, J.R. Position and attitude control of a spacecraft using the state-dependent Riccati equation technique. In Proceedings of the 2000 American Control Conference, ACC (IEEE Cat. No. 00CH36334), Chicago, IL, USA, 28–30 June 2000; Volume 3, pp. 1867–1871. [Google Scholar]
  14. Razzaghi, P.; Al Khatib, E.; Bakhtiari, S. Sliding mode and SDRE control laws on a tethered satellite system to de-orbit space debris. Adv. Space Res. 2019, 64, 18–27. [Google Scholar] [CrossRef]
  15. Weickgenannt, M.; Zimmert, N.; Klumpp, S.; Sawodny, O. Application of SDRE control to servopneumatic drives. In Proceedings of the 2010 IEEE International Conference on Control Applications, Yokohama, Japan, 8–10 September 2010; pp. 1725–1730. [Google Scholar]
  16. Jing-Liang, S.; Chun-Sheng, L.; Ke, L.; Hao-Ming, S. Optimal robust control for attitude of quad-rotor aircraft based on sdre. In Proceedings of the 2015 34th Chinese Control Conference (CCC), Hangzhou, China, 28–30 July 2015; pp. 2333–2337. [Google Scholar]
  17. Guo, R.; Wu, A.; Lang, Z.; Zhang, X. A nonlinear attitude control method for an unmanned helicopter. In Proceedings of the 2010 2nd International Asia Conference on Informatics in Control, Automation and Robotics (CAR 2010), Wuhan, China, 6–7 March 2010; Volume 1, pp. 166–169. [Google Scholar]
  18. Guo, R.; Chen, J. Sdre attitude control with global asymptotic stability for an unmanned helicopter. In Proceedings of the 2011 International Conference on Transportation, Mechanical, and Electrical Engineering (TMEE), Changchun, China, 16–18 December 2011; pp. 1044–1049. [Google Scholar]
  19. dos Santos, G.P.; Balthazar, J.M.; Janzen, F.C.; Rocha, R.T.; Nabarrete, A.; Tusset, A.M. Nonlinear dynamics and SDRE control applied to a high-performance aircraft in a longitudinal flight considering atmospheric turbulence in flight. J. Sound Vib. 2018, 436, 273–285. [Google Scholar] [CrossRef]
  20. Chipofya, M.; Lee, D.J. Position and altitude control of a quadcopter using state-dependent Riccati equation (SDRE) control. In Proceedings of the 2017 17th International Conference on Control, Automation and Systems (ICCAS), Jeju, Korea, 18–21 October 2017; pp. 1242–1244. [Google Scholar]
  21. Shihabudheen, K.; Thankachan, J.; Vasista, C. SDRE control of flexible beam manipulator. In Proceedings of the 2015 International Conference on Soft Computing Techniques and Implementations (ICSCTI), Faridabad, India, 8–10 October 2015; pp. 148–154. [Google Scholar]
  22. Nekoo, S.R. Nonlinear closed loop optimal control: A modified state-dependent Riccati equation. ISA Trans. 2013, 52, 285–290. [Google Scholar] [CrossRef] [PubMed]
  23. Roudkenary, K.A.; Khaloozadeh, H.; Sedigh, A.K. SDRE control of non-affine systems. In Proceedings of the 2016 4th International Conference on Control, Instrumentation, and Automation (ICCIA), Qazvin, Iran, 27–28 January 2016; pp. 239–244. [Google Scholar]
  24. Itik, M.; Salamci, M.U.; Banks, S.P. SDRE optimal control of drug administration in cancer treatment. Turk. J. Electr. Eng. Comput. Sci. 2010, 18, 715–730. [Google Scholar]
  25. Banks, H.; Lewis, B.; Tran, H.T. Nonlinear feedback controllers and compensators: A state-dependent Riccati equation approach. Comput. Optim. Appl. 2007, 37, 177–218. [Google Scholar] [CrossRef] [Green Version]
  26. Voos, H. Nonlinear state-dependent Riccati equation control of a quadrotor UAV. In Proceedings of the 2006 IEEE Conference on Computer Aided Control System Design, 2006 IEEE International Conference on Control Applications, 2006 IEEE International Symposium on Intelligent Control, Munich, Germany, 4–6 October 2006; pp. 2547–2552. [Google Scholar]
  27. Barata, J.C.A.; Hussein, M.S. The Moore–Penrose pseudoinverse: A tutorial review of the theory. Braz. J. Phys. 2012, 42, 146–165. [Google Scholar] [CrossRef] [Green Version]
  28. Knobloch, H.W.; Kwakernaak, H. Lineare Kontrolltheorie; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  29. Kwakernaak, H.; Sivan, R. Linear Optimal Control Systems; Wiley: New York, NY, USA, 1972; Volume 1. [Google Scholar]
  30. Anderson, B.D.; Moore, J.B. Optimal Control: Linear Quadratic Methods; Courier Corporation: Tokyo, Japan, 2007. [Google Scholar]
  31. Wernli, A.; Cook, G. Suboptimal control for the nonlinear quadratic regulator problem. Automatica 1975, 11, 75–84. [Google Scholar] [CrossRef]
  32. Cloutier, J.R.; D’Souza, C.N.; Mracek, C.P. Nonlinear regulation and nonlinear H control via the state-dependent Riccati equation technique: Part 1, theory. In Proceedings of the First International Conference on Nonlinear Problems in Aviation and Aerospace, Daytona Beach, FL, USA, 9–11 May 1996; Embry-Riddle Aeronautical Univ. Press: Daytona Beach, FL, USA, 1996; pp. 117–130. [Google Scholar]
  33. Cloutier, J.R.; Stansbery, D.T.; Sznaier, M. On the recoverability of nonlinear state feedback laws by extended linearization control techniques. In Proceedings of the 1999 American Control Conference (Cat. No. 99CH36251), San Diego, CA, USA, 2–4 June 1999; Volume 3, pp. 1515–1519. [Google Scholar]
  34. Liang, Y.W.; Lin, L.G. Analysis of SDC matrices for successfully implementing the SDRE scheme. Automatica 2013, 49, 3120–3124. [Google Scholar] [CrossRef]
  35. Meyer, C.D., Jr. Generalized inversion of modified matrices. SIAM J. Appl. Math. 1973, 24, 315–323. [Google Scholar] [CrossRef]
  36. Kaczorek, T.; Dzieliński, A.; Dabrowski, W.; Łopatka, R. Podstawy Teorii Sterowania; Wydawnictwo WNT: Warsaw, Poland, 2014. [Google Scholar]
  37. Erdem, E.; Alleyne, A.G. Globally stabilizing second order nonlinear systems by SDRE control. In Proceedings of the 1999 American Control Conference (Cat. No. 99CH36251), San Diego, CA, USA, 2–4 June 1999; Volume 4, pp. 2501–2505. [Google Scholar]
  38. Brauer, F.; Nohel, J.A. The Qualitative Theory of Ordinary Differential Equations: An Introduction; Courier Corporation: Tokyo, Japan, 2012. [Google Scholar]
  39. Heydari, A.; Balakrishnan, S. Closed-form solution to finite-horizon suboptimal control of nonlinear systems. Int. J. Robust Nonlinear Control 2015, 25, 2687–2704. [Google Scholar] [CrossRef]
  40. Zhao, W.; Go, T.H. Quadcopter formation flight control combining MPC and robust feedback linearization. J. Frankl. Inst. 2014, 351, 1335–1355. [Google Scholar] [CrossRef]
Figure 1. Structure of a quadrotor aircraft.
Figure 1. Structure of a quadrotor aircraft.
Applsci 11 10714 g001
Figure 2. Time plot of angles.
Figure 2. Time plot of angles.
Applsci 11 10714 g002
Figure 3. Time plot of velocities.
Figure 3. Time plot of velocities.
Applsci 11 10714 g003
Figure 4. Time plot of control.
Figure 4. Time plot of control.
Applsci 11 10714 g004
Table 1. Parameters used in simulation.
Table 1. Parameters used in simulation.
ParameterValueUnit
g9.81 m s 2
m0.5kg
l0.3m
I x 0.0081kgm 2
I y 0.0081kgm 2
I z 0.0162kgm 2
J R 0.01N
Table 2. Simulation time.
Table 2. Simulation time.
Comparison of Simulation Time
classical SDRE
R 1 , Q 100 11.276407 s
R 1 , Q 1 11.978995 s
R 100 , Q 1 11.224440 s
modified SDRE
R 1 , Q 100 2.893052 s
R 1 , Q 1 2.433471 s
R 100 , Q 1 2.590412 s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stępień, S.; Superczyńska, P. Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control. Appl. Sci. 2021, 11, 10714. https://doi.org/10.3390/app112210714

AMA Style

Stępień S, Superczyńska P. Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control. Applied Sciences. 2021; 11(22):10714. https://doi.org/10.3390/app112210714

Chicago/Turabian Style

Stępień, Sławomir, and Paulina Superczyńska. 2021. "Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control" Applied Sciences 11, no. 22: 10714. https://doi.org/10.3390/app112210714

APA Style

Stępień, S., & Superczyńska, P. (2021). Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control. Applied Sciences, 11(22), 10714. https://doi.org/10.3390/app112210714

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop