Next Article in Journal
Is Green Recovery Enough? Analysing the Impacts of Post-COVID-19 Economic Packages
Next Article in Special Issue
Power Flows and Losses Calculation in Radial Networks by Representing the Network Topology in the Hierarchical Structure Form
Previous Article in Journal
Application of the Harmonic Balance Method for Spatial Harmonic Interactions Analysis in Axial Flux PM Generators
Previous Article in Special Issue
ENIGMA—A Centralised Supervisory Controller for Enhanced Onboard Electrical Energy Management with Model in the Loop Demonstration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Krasovskii Passivity and μ-Synthesis Controller Design for Quasi-Linear Affine Systems

Department of Automation, Technical University of Cluj-Napoca, Str. G. Bariţiu nr. 26-28, 400027 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2021, 14(17), 5571; https://doi.org/10.3390/en14175571
Submission received: 29 June 2021 / Revised: 15 August 2021 / Accepted: 31 August 2021 / Published: 6 September 2021
(This article belongs to the Special Issue Advanced Control Strategies for Electric Power Management)

Abstract

:
This paper presents an end-to-end method to design passivity-based controllers (PBC) for a class of input-affine nonlinear systems, named quasi-linear affine. The approach is developed using Krasovskii’s method to design a Lyapunov function for studying the asymptotic stability, and a sufficient condition to construct a storage function is given, along with a supply-rate function. The linear fractional transformation interconnection between the nonlinear system and the Krasovskii PBC (K-PBC) results in a system which manages to follow the provided input trajectory. However, given that the input and output of the closed-loop system do not have the same physical significance, a path planning is mandatory. For the path planning component, we propose a robust controller designed using the μ -synthesis mixed-sensitivity loop-shaping for the linearized system around a desired equilibrium point. As a case study, we present the proposed methodology for DC-DC converters in a unified manner, giving sufficient conditions for such systems to be Krasovskii passive in terms of Linear Matrix Inequalities (LMIs), along with the possibility to compute both the K-PBC and robust controller alike.

1. Introduction

In the recent years, passivity theory has gained renewed attention because of its advantages and practicality in modelling of multi-domain systems and constructive control techniques. Starting with the classical theory of dissipative systems of [1,2], the idea of differential passivity based on the idea of variational systems was developed in [3,4]. The same authors extended in [5] the idea to the incremental passivity concept. More recently, the idea of Krasovskii’s passivity was developed in [6]. This approach is based on the construction of a storage function using Krasovskii’s method to design a Lyapunov candidate function for stability as in [2]. Moreover, all these methods are useful for designing a simpler and more robust passivity-based controller (PBC).
However, a PBC manages to guarantee the closed-loop stability and reference tracking in terms of the command signal. But, in many cases, systems do not have the same physical significance and relevance for both input and output signals and, as such, an additional path planning procedure is necessary. For the purpose of this paper, we propose a robust controller to play this role, because robustness represents a major problem studied in Control Theory, with high impact in the literature, encompassing the sensitivity of a system with respect to both internal and external disturbances. The classical solutions for robust control problems uses H 2 and H norms as a performance measure. The most common approaches to solve such problems are Algebraic Riccati Equations (AREs) [7] and Algebraic Riccati Inequalities (ARIs) [8]. In order to solve AREs in a more numerically stable manner, a solution using Popov triplets was presented in [9]. An open-source implementation of the above method is presented in [10], along with an iterative refinement described in [11]. In order to overcome the limitations of the ARE-based solutions, consisting in the impossibility of solving singular problems, a Linear Matrix Inequality (LMI) based solution was developed in [12]. For solving robust and nonlinear control problems using LMIs, an open-source implementation is described in [13]. Moreover, an open-source toolbox for automatically modelling uncertainties as required in robust control problems using MATLAB is presented in the paper [14].
Nowadays, renewable energy sources play an increasingly important role in the energy conversion, due to the negative consequences of the global warming. DC-DC converters have been found to be an adequate solution regarding renewable energy harvesting such as for sustaining prescribed voltages for high-power loads and for maximum power point tracking, as in [15,16,17], or in DC microgrids, as in [18]. An interior point method for renewable energy management was introduced in [19], where the importance of the good performance in the DC-DC converters control is underlined. For modelling DC-DC converters, there are two major improvements: the averaging of ON and OFF state-space models from [20], obtaining a nonlinear model, along with adding parasitic components in the circuit model, such as series and parallel resistors in the switching elements and passive components alike, as in [21]. The average model greatly simplifies the computations and well approximates the physical switching system to a desired ripple tolerance, based on the PWM period. Using these ideas, a technique which adds a damping injection term to reach the passivity of the system was presented in [22], and was extended for the nonideal DC-DC boost converter in [23]. A Krasovskii passive proportional-integral (PI) type controller was designed for the ideal case of DC-DC converters in [6].
The main contributions of the current paper are: a mathematical framework for studying the property of quasi-linear input-affine nonlinear systems to be passive, a method to design a passivity-based controller for such a nonlinear systems, along with an entire software mechanism which manages to compute both PBC and robust controllers. First, we provide an extended set of sufficient conditions for a quasi-linear input-affine nonlinear system to be Krasovskii passive, along with a set of necessary and sufficient conditions for DC-DC buck, boost and single-ended primary inductor (SEPIC) converters to construct an output such that the systems are Krasovskii passive. These extended methods can be used even for nonideal converters, compared to the sufficient set of conditions presented in [6]. Also, we present a method to design a Krasovskii passivity-based controller (K-PBC) which, along with the output port variable of the DC-DC converter, manages to ensure the closed-loop stability. This means that the lower linear fractional transformation between the system and the controller is asymptotically stable and manages to track any provided input trajectory. However, the input and the output signals are of different physical significance and, as such, an input trajectory path planning is needed in order to ensure a desired set of tracking performances. As a solution to this issue, we propose a robust controller designed for a linearized system around a forced equilibrium point. All these steps are included in an open-source toolbox which extends the initial iteration presented in [14].
The paper follows the structure: Section 2 describes the mathematical tools necessary for Krasovskii passivity analysis and Krasovskii passivity-based controller (K-PBC) design, along with a basic overview of Robust Control theory; Section 3 gives an extended set of sufficient conditions for a quasi-linear input-affine nonlinear system to be Krasovskii passive and presents the proposed structure for the decentralized controller, along with relevant aspects and observations regarding the software implementation; Section 4 presents the mathematical models of the nonideal DC-DC converters, the Krasovskii’s passivity analysis of these systems and the K-PBC design manner; Section 5 provides comparative simulations of the obtained closed-loop system obtained by numerical simulation obtained using the proposed open-source toolbox; Section 6 illustrates the conclusions, along with proposed improvements and future work completions.
Notation: For x R n and Q = Q 0 R n × n , we denote by | | x | | Q : = ( x Q x ) 1 / 2 .

2. Mathematical Background

2.1. Passivity-Based Control

The control method proposed in this paper is designed for a certain class of continuous time input-affine nonlinear systems. First, we consider the general form of a finite-order continuous time input-affine nonlinear system in order to describe and develop the main mathematical background:
( Σ ) : x ˙ = f ( x , u ) : = g 0 ( x ) + i = 1 m g i ( x ) u i ,
where x R n is the state vector, u = [ u 1 , , u m ] R m is the input vector and the functions g i : R n R n are of class C 1 . Consider the set of forced equilibrium points:
E = { ( x ¯ , u ¯ ) R n × R m : f ( x ¯ , u ¯ ) = 0 } .
We assume that the set E is nonempty. The notion of dissipative systems is presented according to [2] in the following definition.
Definition 1.
The system (Σ) described in Equation (1) is said to be dissipative with respect to the supply rate ω : R n × R m R if there exists a storage function S : R n R + of class C 1 such that:
S ( x ) x f ( x , u ) ω ( x , u ) , ( x , u ) R n × R m .
We assume that the supply-rate can be written as ω ( x , u ) = u h ( x ) , where h : R n R m . In this case, vectors u and y = h ( x ) are called port variables.
The critical difficulty in order to prove the passivity is characterized by a suitable choice of the storage function. A typical choice of the storage function is the total energy of the system. Although these types of storage functions can be used to prove the passivity of a system, this standard approach is not useful to design a passivity-based controller due to the port variable. A solution with a damping injection term is presented by the authors in [23]. In order to design a passivity based controller (PBC), the following extended system can be considered, as in [24]:
( Σ e ) : x ˙ = f ( x , u ) ; u ˙ = u d ,
where u d R m is the new input vector and x ˜ = [ x u ] R n + m is the new state vector. A new passivity concept was introduced by [6] using a storage function similar to the Lyapunov candidate function constructed using Krasovskii’s method presented in [2].
Definition 2.
Let h K : R n × R m R m be the function which describes the output port variable. The nonlinear system (Σ) is said to be Krasovskii passive if its extended system ( Σ e ) is dissipative with respect to the supply rate:
ω K ( x ˜ , u d ) = u d h K ( x ˜ ) ,
with a storage function:
S K ( x ˜ ) = 1 2 | | f ( x , u ) | | Q 2 , x ˜ ( x , u ) R n × R m
where Q = Q 0 .
A sufficient set of conditions for a system to be Krasovskii passive is presented in the following theorem. For simplicity, for a symmetric and positive semidefinite matrix Q R n × n we use the following shorthand notation:
Q g k ( x ) = Q g k ( x ) x + g k ( x ) x Q , k = 0 , m ¯ .
Theorem 1.
The system ( Σ ) is Krasovskii passive with the supply-rate ω K ( x ˜ , u d ) = u d h K ( x ˜ ) , where the port variable h K can be expressed as:
h K ( x , u ) = [ g 1 ( x ) g m ( x ) ] Q f ( x , u )
and with the storage function:
S K ( x ˜ ) S K ( x , u ) = 1 2 | | f ( x , u ) | | Q 2 ,
if there exists a matrix Q R n × n , Q = Q 0 , that satisfies the following condition:
Q ( x ˜ ) Q g 0 ( x ) + k = 1 m Q g k ( x ) u k 0 ,
for all x ˜ ( x , u ) R n × R m .
Proof. 
Taking the Lie derivative of S K along the vector field of ( Σ e ) :
S K ( x ˜ ) x ˜ f e ( x ˜ , u d ) = f ( x ˜ ) Q ( x ˜ ) f ( x ˜ ) + u d h K ( x ˜ ) ,
where
f e ( x ˜ , u d ) = f ( x , u ) u d ,
describes the system ( Σ e ) : x ˜ ˙ = f e ( x ˜ , u d ) .
From Q ( x ˜ ) 0 we obtain the inequality:
S K ( x ˜ ) x ˜ f e ( x ˜ , u d ) u d h K ( x , u ) = ω K ( x ˜ , u d ) ,
which completes the proof.    □
The sufficient condition indicated in the theorem is an extension of the necessary and sufficient set of conditions for a system to be differential passive, given in [4], which assert that Q g 0 ( x ) 0 and Q g k ( x ) = 0 , for each x R n . These necessary and sufficient conditions were extended in [6] to be sufficient conditions to Krasovskii passivity, but this approach can also lead to a port variable unusable for controller design, as it will be shown in Section 4.
Next, we want to prove that the classical interconnection between two Krasovskii passive systems presented in Figure 1 forms also a Krasovskii passive systems. Consider two such systems ( Σ 1 ) and ( Σ 2 ) given by:
( Σ j ) : x ˙ j = f ( x j , u j ) , j { 1 , 2 } ,
having the extended systems ( Σ e , 1 ) and ( Σ e , 2 ) with the states x ˜ j ( x j , u j ) , which are Krasovskii passive with respect to the supply rates ω K , j ( x ˜ j , u d , j ) = u d , j h K , j ( x ˜ j ) and the storage functions S K , j ( x ˜ j ) . For the systems ( Σ e , j ) , the interconnection considered in Figure 1 can be written as:
( Σ i ) : u d , 1 u d , 2 = 0 1 1 0 h K , 1 h K , 2 + e d 1 e d 2 ,
where e 1 , e 2 R m are external inputs and e d , j = e ˙ j are their derivatives. The extended system state can be expressed as
x ˜ = [ x ˜ 1 + [ 0 n e 1 ] x ˜ 2 + [ 0 n e 2 ] ] .
Lemma 1
([6]). The resulting system ( Σ i ) is Krasovskii passive with respect to the supply-rate
ω K , i ( x ˜ , e d 1 , e d 2 ) = e d 1 h K , 1 ( x ˜ 1 ) + e d 2 h K , 2 ( x ˜ 2 )
and the storage function
S K , i ( x ˜ ) = S K , 1 ( x ˜ 1 ) + S K , 2 ( x ˜ 2 ) .
Proof. 
From the assumptions of the lemma it results that:
S K , j ( x ˜ j ) x ˜ j f e , j ( x ˜ j , u d , j ) u d , j h K , j ( x ˜ j ) , j { 1 , 2 } .
Consider the Lie derivative of S K , i along the vector fields of ( Σ e , 1 ) and ( Σ e , 2 ) and using the previous relations:
S K , i ( x ˜ ) x ˜ f e , 1 ( x ˜ 1 , u d 1 ) f e , 2 ( x ˜ 2 , u d 2 ) u d , 1 h K , 1 ( x ˜ 1 ) + u d , 2 h K , 2 ( x ˜ 2 ) =
= e d , 1 h K , 1 ( x ˜ 1 ) + e d , 2 h K , 2 ( x ˜ 2 ) = ω K , i ( x ˜ , e d 1 , e d 2 ) ,
which concludes the proof.    □
Now, we show a procedure to design a Krasovskii passivity-based controller (K-PBC). Let us define:
( Σ c ) : y c ( t ) = x ˙ c ( t ) = K 1 ( K 2 x c u c ( t ) ) f c ( x c , u c ) ,
where x c R n c is the state vector, u c and y c are the input vector and output vector, respectively, while the matrices K 1 , K 2 R n c × n c must be symmetric and negative definite, along with symmetric and positive definite, respectively.
Lemma 2
([6]). The controller ( Σ c ) is Krasovskii passive with respect to the supply-rate ω K , c ( y c , u c ) = y c u c and with the storage function S K , c = 1 2 | | x c | | K 2 2 .
Proof. 
Consider the Lie derivative of the S K , c along the vector field of ( Σ c ) :
S K , c ( x c ) x c f c ( x c , u c ) = x ˙ c K 2 x c = x ˙ c ( K 1 1 x ˙ c + u c ) .
Using K 1 < 0 and y c = x ˙ c , we obtain:
S K , c ( x c ) x c f c ( x c , u c ) x ˙ c u c = ω K , c ( y c , u c ) ,
and the proof is done.    □
Now, consider the interconnection of the systems ( Σ e ) and ( Σ c ) given by:
u d u c = 0 1 1 0 h K y c + 0 u 🟉 ,
where u 🟉 is the external input. The following theorem gives the sufficient condition for which the interconnection described above is Krasovskii passive.
Theorem 2
([2]). Consider the extended system ( Σ e ) which satisfies the conditions of Theorem 1 for a symmetric and positive definite matrix Q R n × n with at least one equilibrium isolated point ( x ¯ , u ¯ ) . Also, consider the controller given by ( Σ c ) with the states x c = u 🟉 u and the interconnection given by Equation (13). The closed-loop system is dissipative with respect to the supply rate ω o ( u d , u 🟉 ) = u d u 🟉 with the storage function:
S o ( x , u , x c ) = 1 2 | | f ( x , u ) | | Q 2 + 1 2 | | x c | | K 2 2 .
Proof. 
Consider x ˜ o = [ x ˜ x c ] . The Lie derivative of the storage function S o ( x ˜ o ) along the closed-loop system trajectory respects the inequality from Lemma 1:
S o ( x ˜ o ) x ˜ o f e ( x ˜ , u d ) f c ( x c , u c ) ( u 🟉 ) · y c .
   □
All these results will be used to develop a unified approach to design K-PBC for a class of input-affine nonlinear systems in Section 3, which will be particularized to the case of DC-DC converters in Section 4 in a unified manner.

2.2. Robust Control

The proposed K-PBC will manage to ensure that the closed-loop system follows the input trajectory u 🟉 , but in most applications, the desired trajectory is with regards to the output signal, of different physical significance compared to the input signal. Therefore, this command signal needs to be computed using a second controller. For the purpose of this paper, we will use the Robust Control Framework in order to compute the input trajectories. This framework is designed for the class of linear and time invariant systems (LTI). The generalized plant used in the Robust Control Framework considering two types of uncertainties is the following:
P Δ : x ˙ ( t ) i n e v ( t ) z ( t ) y ( t ) = A i n e C v C z C y | B d D v d D z d D y d B w D v w D z w D y w B u D v u D z u D y u x ( t ) i n e d ( t ) w ( t ) u ( t ) ,
where the vectors presented are: the disturbance input d R n d , the exogenous input w R n w , the control input u R n u , the disturbance output v R n v , the performance output z R n z and the measurement output y R n y . The disturbance vector d encompasses two types of uncertainties: parametric and unstructured. These uncertainties are modelled using the following set:
Δ = Δ = diag δ 1 I n 1 , , δ s I n s , Δ 1 , , Δ f | δ k R , Δ j R m j × m j , k = 1 , s ¯ , j = 1 , f ¯ .
The entire structure is presented in Figure 2. As observed in the figure, the generalized plant has an upper linear fractional transformation (ULFT) with the uncertainty block Δ and a lower linear fractional transformation (LLFT) with the controller K. The goal of the Robust Controller is to ensure both robust stability and robust performance. In order to deal with uncertainties, the mathematical framework used for this is represented by the structural singular value:
μ Δ ( LLFT ( P , K ) ) = sup ω R + 1 min Δ Δ { σ ¯ ( Δ ) | det ( I LLFT ( P , K ) ( j ω ) Δ ) = 0 } .
However, the structural singular values are hard to be explicitly computed, therefore an upper bound can be used [25]:
μ Δ ( LLFT ( P , K ) ( j ω ) ) inf D D σ ¯ ( D · LLFT ( P , K ) ( j ω ) · D 1 ) ,
where the set D is defined as:
D = diag D 1 , , D s , d 1 I m 1 , , d f I m f | D k = D k R n k × n k , d j > 0 , k = 1 , s ¯ , j = 1 , f ¯ .
Using the main loop theorem, the robust stability and robust performance can be guaranteed by the small gain principle:
μ Δ ( LLFT ( P , K ) ) < 1 ,
which means that the robust control problem can be formulated in terms of optimization:
inf K stab . sup ω R + μ Δ ( LLFT ( P , K ) ( j ω ) ) ,
and, using Equation (18), the following quasi-convex problem results:
inf K stab . sup ω R + inf D D σ ¯ D ( j ω ) · LLFT ( P , K ) ( j ω ) · ( D ( j ω ) ) 1 .
As detailed in [14,26], the μ -synthesis control problem can be solved using the so-called DK iteration procedure. This starts with an initial guess for D (usually the unitary matrix), then the controller step, consisting in solving a H control problem, and a scaling step, consisting in solving several Parrot problems in a prescribed set of frequencies followed by fitting a nonminimum-phase system using the obtained solutions, are executed in a loop sequence until a stopping criterion is fulfilled.
The robust control problem used in this paper is the mixed sensitivity loop-shaping. In this procedure, the LTI plant G is augmented using three weighting functions, one for each specific function: sensitivity S, complementary sensitivity T and control effort K S . This controller design technique provides a good trade-off between stability and performance, and can be expressed as:
inf K stab . sup ω R + inf D D σ ¯ D ( j ω ) · LLFT ( P , K ) ( j ω ) · ( D ( j ω ) ) 1 , such that W S S W T T W R R T < 1 ,
where W S , W T and W K S are dynamic weighting functions for the closed-loop functions mentioned above, selected in order to penalise certain frequencies in different amounts.

3. Krasoviskii Passivity-Based Controller Design with μ -Synthesis Path Planning

3.1. Proposed Method

For the purpose of this paper, we consider the following nonlinear input-affine systems in order to develop the proposed controller design technique. First, some initial considerations should be made.
Definition 3.
A nonlinear input-affine system is calledquasi-linear affineif
( Σ a ) : x ˙ = f ( x , u ) = g 0 ( x ) + k = 1 m g k ( x ) u k ,
where the nonlinear functions are affine, i.e., g k ( x ) = A k x + b k , for each k = 0 , m ¯ , with A k R n × n and b k R n .
Considering that:
g k ( x ) x = A k , k = 0 , m ¯ ,
and using Theorem 1, the necessary and sufficient condition for a system ( Σ a ) to be Krasovskii passive can be written as:
Q A 0 + A 0 Q + i = 1 m Q A i + A i Q u i 0 , u = [ u 1 , u 2 , , u m ] R m .
But, due to the physical constraints of the process, the command vector is bounded: u [ u 1 l , u 1 u ] × [ u 2 l , u 2 u ] × × [ u m l , u m u ] . Given that the space spanned by the values of the command signal describes a polytope, Equation (26) can be formulated as a reunion of LMI problems, one for each vertex of the polytope. Using this remark, the necessary and sufficient conditions for a quasi-linear affine input-affine system ( Σ a ) to be Krasovskii passive are presented in the following theorem.
Theorem 3.
The system ( Σ a ) is Krasovskii passive with the supply-rate ω K ( x , u ) = u h K ( x ) , where the port variable h K can be expressed as:
h K ( x , u ) = [ x A 1 + b 1 x A m + b m ] · Q · x ˙ ,
and with the storage function:
S K ( x ) = 1 2 x ˙ Q 2
if and only if there exists a symmetrical matrix Q = Q 0 R n × n which satisfies the following conditions:
Q A 0 + A 0 Q + i = 1 m Q A i + A i Q e i u i l + Q A i + A i Q ( 1 e i ) u i u 0 ,
for each binary word e = e 1 e 2 e m Z 2 m .
Proof. 
Taking the Lie derivative of S K along the vector field of ( Σ a ) :
S K ( x ) x f ( x , u ) = x ˙ · Q A 0 + A 0 Q + i = 1 m Q A i + A i Q u i · x ˙ + u h K ( x ) .
Now, in order to impose the desired behaviour, we need the inequality:
S K ( x ) x f ( x , u ) u h K ( x )
to be true for each state trajectory x R n and for each input trajectory u R m , which is equivalent with the following LMIs:
Q A 0 + A 0 Q + i = 1 m Q A i + A i Q u i 0
to have a common solution for each command signal u [ u 1 l , u 1 u ] × [ u 2 l , u 2 u ] × × [ u m l , u m u ] , which is a polytope problem and can be solved in its vertices only, and the proof is complete.    □
After the output port-variable is computed as in the previous theorem, a K-PBC, having Equation (12), is computed such that the closed-loop system is Krasosvkii passive. This property ensures asymptotic stability. As such, the closed-loop system manages to follow the input trajectory u 🟉 . However, in order to give a trajectory according to the desired reference for the output y 🟉 , another component, i.e., trajectory path planning, is needed. For this component we use a robust controller K rob , obtaining the cascade structure presented in Figure 3.
The robust controller which computes the input trajectory for the K-PBC is designed for a linearized model of the system around a forced equilibrium point x ¯ ˜ = x ¯ u ¯ . The LTI model of the linearized plant of the input-affine quasi-linear system ( Σ a ) can be computed as:
A = f ( x , u ) x x ˜ = x ¯ ˜ = A 0 + k = 1 m A i u ¯ i ; B = f ( x , u ) u x ˜ = x ¯ ˜ = k = 1 m A i x ¯ + b i ,
followed by a similar extension to the output matrices C and D by replacing the state equation function f with the output function, usually denoted by h. Now, the linearized plant can be represented as:
( Σ a , l i n ) : Δ x ˙ ( t ) = A Δ x ( t ) + B Δ u ( t ) ; Δ y ( t ) = C Δ x ( t ) + D Δ u ( t ) ,
and now all the procedures described in the paper [14] can be applied in order to obtain the uncertainty block Δ , the augmented plant and, thus, to prepare the augmented plant for synthesizing the robust controller. Figure 4 encompasses all these steps. After the mixed sensitivity loop-shaping μ -synthesis control Equation (23) is solved, the order reduction for the resulting high-order controller is applied. All these details are underlined in the previously mentioned paper and are summarily described in next subsection, where the software integration of the new features is presented.

3.2. Software Implementation

The software framework developed and implemented in [14] contains two main original functionalities, which are also extended in the present paper: the former regarding model-in-the-loop (MiL) simulations for a wide range of system configurations in a unified manner, while the latter allows automatic computation of uncertain plant sets near an operating point as required for robust control synthesis, in the sense of Figure 4.
Starting from said framework, an extension of the class diagram previously-published is illustrated in Figure 5. The building block for all provided system configurations starts from the abstract class System, which provides an interface for finite order, possibly time-varying dynamical systems, with a state and output equation, respectively:
x ˙ ( t ) = f x ( t ) , u ( t ) , t ; y ( t ) = h x ( t ) , u ( t ) , t ,
with the vector-valued maps f : R n + m × R + R n and h : R n + m × R + R p representing multivariable Lipschitz functions. In general, the input u ( t ) size is denoted by m, state signal x ( t ) has dimension n, while the output signal y ( t ) has size p. The system initial conditions are denoted x ( 0 ) = x 0 R n . Using this interface, linear and time-invariant systems were supported, along with affine LTI systems and system interconnections, such as series, parallel, LLFT and ULFT, respectively. Another means of branching the provided System class for switching systems and other hybrid-driven behaviour is through the abstract class HybridSystem, which also provides the previously-mentioned interconnections. The abstract class UncertainPlantFactory provides means for specifying and instantiating System objects based on a family of uncertain plant dynamics. It provides two interface methods to return such systems through getNominalPlant and getRandomPlant. Based on the considered nominal plant, by convention, it also automatically generates uncertainty models using the standard models considered in the literature, such as input/output multiplicative, additive and inverse additive uncertainties and so on. Additionally, the class ExtendedDcDcConverter implements the input-affine systems of Equation (4), extended with the port variable h K ( x , u ) from Equation (6). In a similar manner, the class KrasovskiiPassivityBasedController implements the K-PBC system as described in Equation (12). An intermediary class which provides an interface between the plant ( Σ e ) with K-PBC K K P B C in Equation (13) and the robust controller K r o b is implemented in KPBCWrapper, as in the left part of Figure 3.
The K-PBC and μ -synthesis controllers are decoupled, can be used individually to control the quasi-linear affine plant, but are used to gather the benefits of both approaches, as the K-PBC assures asymptotic stability of the nonlinear plant, while the μ -synthesis controller is used for robust stability and performance around the operating point. A final LLFT connection, implemented by means of class LLFTConnection, as in Figure 3, gives the proposed closed-loop system.
For the MiL component of the toolbox, closed-loop simulations will be conducted through the high-level interface of the abstract class System which links all system state and output equation definitions with the ordinary differential equation ode solver framework, inherited through the many subclasses from Figure 5 and, also mentioned in the previous paragraph. An auxiliary functionality was created in order to specify various input signals, such as combinations of steps, ramps, sine waves and interpolated lookup tables, on desired input channels, as necessary for MISO and MIMO systems, based on function handles, cell arrays and flexible low-level signal definitions.
The software structure and workflow for the automatic uncertainty set computation for a family of plants can be described in the following manner: the toolbox user must inherit the abstract UncertainPlantFactory class and adapt it with means of returning a nominal plant System object through the method getNominalPlant(), along with means of returning a random plant System object from the considered uncertainty set through the method getRandomPlant(). This interface allows encapsulating the main required functionalities of the class, without restricting the plant structure and definition, allowing the support for parametric uncertainties, unstructured uncertainties, or both. After this definition, a Monte Carlo simulation is conducted based on the nominal and randomized instances of the nonlinear plant, followed by a linearization procedure around each particular instance’s adapted operating point using the System routine linearize, resulting in a family of models of Equation (31). Finally, an uncertainty model Δ of given structure is fitted, as needed for the μ -synthesis procedure. All these steps are computed automatically based on a set of specifications encompassing the desired operating point, the relevant frequency domain of the plant, number of randomized instances, uncertainty type (additive, inverse additive, input multiplicative, inverse input multiplicative and so on), order and structure of the transfer functions involved in the optimization.
The μ -synthesis controller can then be computed using the results obtained in the previous paragraph, after augmenting the uncertain LTI plant family using the well-known closed-loop-shaping method, by weighting the sensitivity, complementary sensitivity and control effort functions, respectively, based on the optimization Equation (23). Additionally, an order reduction step may be performed on the resulting controller.

4. Case Study: DC-DC Converter Control—An Unified Approach

A relevant class of quasi-linear affine systems with respect to the command signal is that of DC-DC converter circuits. In the first subsection, we will provide a mathematical framework applicable to all DC-DC converter topologies, followed by three illustrations of the framework’s flexibility on buck, boost and SEPIC converters alike.

4.1. Mathematical Models

A nonideal DC-DC converter has a set of switching elements S i , i 2 , with at least one of them being a transistor and, otherwise, represented by other transistors with synchronized or complementary switching, or directly diodes. Among the advantages of using multiple transistors, one can recall their smaller voltage drops, along with ease of manufacturing. The components of such a converter are:
  • L i : converter inductors, i = 1 , k 1 ¯ ;
  • r L i : resistances associated with the inductors, i = 1 , k 1 ¯ ;
  • C j : converter capacitors, j = 1 , k 2 ¯ ;
  • r C j : resistances associated with the output capacitors, j = 1 , k 2 ¯ ;
  • R: variable output load resistor;
  • E: source voltage;
  • r D S i : parasitic resistances associated with the ON state of the switching elements;
  • V F i : constant voltage drops associated with the ON state of the switching elements;
  • μ : normalized duty cycle of the main switching device, i.e., μ 0 , 1 .
The previously-mentioned components and notations are illustrated in Figure 6 for the three classic circuit topologies studied in this section.
The mathematical modelling of the converter circuits is made using Kirchhoff’s laws. As such, the state variable vector will be defined as:
x ( t ) = x L ( t ) x C ( t ) = i L 1 ( t ) i L k 1 ( t ) | u C 1 ( t ) u C k 2 ( t ) ,
storing the inductor currents i L i , i = 1 , k 1 ¯ and capacitor voltages u C j , j = 1 , k 2 ¯ alike. Having the structure of a commutation system. dependent on at least one switching element, two different LTI state-space models are implied. As such, the ON and OFF states of the transistor provide the state-space models:
x ˙ ( t ) = f O N ( x ( t ) , u ( t ) ) ; x ˙ ( t ) = f O F F ( x ( t ) , u ( t ) ) .
Using the duty cycle command signal in a convex combination of the two previous state-space models, a practical nonlinear approximation of the switching system follows:
x ˙ ( t ) = μ ( t ) · f O N x ( t ) , u ( t ) + ( 1 μ ( t ) ) · f O F F x ( t ) , u ( t ) .
Replacing the vector maps f O N and f O F F with their actual expressions returns the averaged state-space model, which well approximates the physical switching system to a desired voltage and current ripple tolerance, based on the PWM frequency.
The average model of each DC-DC converter circuit can be written as:
( Σ ) : x ˙ = g 0 ( x ) + g 1 ( x ) μ f ( x , μ ) ,
where the functions g i : R n R n can be expressed as g i ( x ) = A i x + b i . Next, we will provide the matrices which describe three topologies: Buck, Boost and SEPIC.

4.2. Buck Converter

The electrical scheme of nonideal DC-DC buck converter is presented in Figure 6, with the electrical components described above. The vector of state variables is:
x ( t ) = x L ( t ) x C ( t ) = i L ( t ) u C ( t ) .
The input-affine quasi-linear model of the system can be written as:
( Σ b u c k ) : x ˙ = A 0 x + b 0 + A 1 x μ + b 1 μ f b u c k ( x , μ ) ,
with the components being:
A 0 = 1 L ( r L + R | | r c r D S ) 1 L R R + r C 1 C R R + r C 1 C ( R + r C ) ;
A 1 = 2 L r D S 0 0 0 ; b 0 = V F L 0 ; b 1 = E L 0 ,
with r C | | R denoting the parallel connection of the resistors r C and R.
The system ( Σ b u c k ) is Krasovskii passive if the following LMI problem has a feasible solution Q = Q 0 :
Q A 0 + A 0 Q 0 ; Q A 0 + A 0 Q + Q A 1 + A 1 Q 0 .
However, the matrix A 0 depends on the value of the resistance R, which is an exogenous input for the system. Therefore, the fist LMI condition can be replaced with:
A 0 = 1 L ( r L r D S 1 r C ) 1 L 1 C 0 A 01 + 1 r C + R r C 2 L r C L r C C 1 C A 02 .
Now, considering a lower and an upper bound for this exogenous input, we can express the passivity of the buck converter using the following claim.
Claim 1.
The DC-DC buck converter system ( Σ b u c k ) is Krasovskii passive for each R [ R m i n , R m a x ] if and only of the following LMI problem has a solution Q = Q 0 :
Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q + Q A 1 + A 1 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q + Q A 1 + A 1 Q 0 .
The supply-rate is ω K ( x , μ ) = μ · h K ( x ) , having the output port-variable:
h K ( x ) = x A 1 + b 1 · Q · x ˙ ,
with the storage function given by:
S K ( x ) = x ˙ Q 2 .
Remark 1.
The passivity analysis using the classical storage function:
S ( x ) = 1 2 L i L 2 + C u C 2
does not lead to a possible PBC, because the Lie derivative of this storage function along the state trajectory x follows:
S ( x ) x f b u c k ( x , μ ) μ E V F i L ,
which means that the system is passive with respect to the voltage source E, voltage drop V F and the inductor current i L , where the voltage source and the voltage drop cannot be controlled.
Remark 2.
The Krasovskii passivity analysis of the nonideal DC-DC buck converter performed using the techniques presented in [6] also leads to the impossibility of constructing a PBC, because by solving the system:
Q A 0 + A 0 Q 0 ; Q A 1 + A 1 Q = 0 , Q = 0 0 0 q , q > 0 ,
having the output port-variable h K ( x ) 0 , thus, a K-PBC cannot be constructed.

4.3. Boost Converter

The electrical circuit of the nonideal DC-DC boost topology is presented in Figure 6, with the electrical components described above. The state variable vector is:
x ( t ) = x L ( t ) x C ( t ) = i L ( t ) u C ( t ) .
The input-affine nonlinear model of the system is:
( Σ b o o s t ) : x ˙ = A 0 x + b 0 + A 1 x μ + b 1 μ ,
where the components are:
A 0 = 1 L ( r L + R | | r c + r D S ) 1 L R R + r C 1 C R R + r C 1 C ( R + r C ) ;
A 1 = 1 L r C | | R R L ( r C + R ) R C ( r C + R ) 0 ; b 0 = E V F L 0 ; b 1 = 0 0 ,
with r C | | R denoting the parallel connection of the resistors r C and R.
The system ( Σ b o o s t ) is Krasovskii passive if the following LMI problem has a feasible solution Q = Q 0 :
Q A 0 + A 0 Q 0 ; Q A 0 + A 0 Q + Q A 1 + A 1 Q 0 .
However, the matrices A 0 and A 1 depend on the value of the resistance R, which is an exogenous input for the system. Therefore, these matrices can be rewritten as:
A 0 = 1 L ( r L + r D S + r C ) 1 L ; 1 C 0 A 01 + 1 r C + R r C 2 L r C L r C C 1 C A 02 ;
A 1 = R r C + R r C L 1 L 1 C 0 A 10 .
Now, considering a lower and an upper bound for this exogenous input, we can express the passivity of the boost converter using the following claim.
Claim 2.
The DC-DC boost converter system ( Σ b o o s t ) is Krasovskii passive for each R [ R m i n , R m a x ] if and only of the following LMI problem has a solution Q = Q 0 :
Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q + R m i n r C + R m a x Q A 10 + A 10 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q + R m a x r C + R m i n Q A 10 + A 10 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q + R m i n r C + R m a x Q A 10 + A 10 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q + R m a x r C + R m i n Q A 10 + A 10 Q 0 .
The supply-rate is ω K ( x , μ ) = μ · h K ( x ) , having the output port-variable:
h K ( x ) = x A 1 + b 1 · Q · x ˙ ,
with the storage function given by:
S K ( x ) = x ˙ Q 2 .

4.4. SEPIC Converter

The electrical schematic of the nonideal DC-DC SEPIC topology is provided in Figure 6, with the electrical components as described above. The state variable vector is:
x ( t ) = x L ( t ) x C ( t ) = u C i n ( t ) i L 1 ( t ) u C 1 ( t ) i L 2 ( t ) u C 2 ( t ) .
The input-affine nonlinear model of the system is:
( Σ s e p i c ) : x ˙ = A 0 x + b 0 + A 1 x μ + b 1 μ ,
where the components are:
A 0 = 1 r C i n C i n 0 0 0 0 1 L 1 r C i n + r L 1 + r C 1 + r D S 1 + r C 2 L 1 1 L 1 r D S + r C 2 L 1 1 L 1 0 1 C 1 0 0 0 0 r D S 2 + r C 2 L 2 0 r D S 2 + r L 2 + r C 2 L 2 1 L 2 0 R C 2 ( R + r C 2 ) 0 R C 2 ( R + r C 2 ) 1 C 2 ( R + r C 2 ) ; b 0 = E r C i n C i n V F L 1 0 V F L 2 0 ;
A 1 = 0 0 0 0 0 0 r C 1 + r C 2 L 1 1 L 1 r C 2 L 1 1 L 1 0 1 C 1 0 1 C 1 0 0 r C 2 + r D S 1 + r D S 2 L 2 1 L 2 r C 2 r C 1 + r D S 2 r D S 1 L 2 1 L 2 0 R C 2 ( R + r C 2 ) 0 R C 2 ( R + r C 2 ) 0 ; b 1 = 0 1 L 1 ( V F 1 V F 2 ) 0 1 L 2 ( V F 1 V F 2 ) 0 ,
with r C 2 | | R being the parallel connection of the resistors r C and R.
The system ( Σ s e p i c ) is Krasovskii passive if the following LMI problem has a feasible solution Q = Q 0 :
Q A 0 + A 0 Q 0 ; Q A 0 + A 0 Q + Q A 1 + A 1 Q 0 .
However, the matrices A 0 and A 1 depend on the value of the resistance R, which is an exogenous input for the system. Therefore, these matrices can be rewritten as:
A 0 = 1 r C i n C i n 0 0 0 0 1 L 1 r C i n + r L 1 + r C 1 + r D S 1 + r C 2 L 1 1 L 1 r D S 2 + r C 2 L 1 1 L 1 0 1 C 1 0 0 0 0 r D S 2 + r C 2 L 2 0 r D S + r L 2 + r C 2 L 2 1 L 2 0 1 C 2 0 1 C 2 0 A 01 + 1 r C 2 + R 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 r C 2 C 2 0 r C 2 C 2 1 C 2 A 02 ;
A 1 = 0 0 0 0 0 0 r C 1 + r C 2 L 1 1 L 1 r C 2 L 1 1 L 1 0 1 C 1 0 1 C 1 0 0 r C 2 + r D S 1 + r D S 2 L 2 1 L 2 r C 2 r C 1 + r D S 2 r D S 1 L 2 1 L 2 0 0 0 0 0 A 10 + R r C + R 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 C 2 0 1 C 2 0 A 11 .
Now, considering a lower and an upper bound for this exogenous input, we can express the passivity of the SEPIC converter using the following claim.
Claim 3.
The DC-DC SEPIC converter system ( Σ b o o s t ) is Krasovskii passive for each R [ R m i n , R m a x ] if and only of the following LMI problem has a solution Q = Q 0 :
Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q + Q A 10 + A 10 Q + R m i n r C + R m a x Q A 11 + A 11 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m i n Q A 02 + A 02 Q + Q A 10 + A 10 Q + R m a x r C + R m i n Q A 11 + A 11 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q + Q A 10 + A 10 Q + R m i n r C + R m a x Q A 11 + A 11 Q 0 ; Q A 01 + A 01 Q + 1 r C + R m a x Q A 02 + A 02 Q + Q A 10 + A 10 Q + R m a x r C + R m i n Q A 11 + A 11 Q 0 .
The supply-rate is ω K ( x , μ ) = μ · h K ( x ) , having the output port-variable:
h K ( x ) = x A 1 + b 1 · Q · x ˙ ,
with the storage function given by:
S K ( x ) = x ˙ Q 2 .

5. Numerical Results

In this section we present detailed design and analysis steps of the proposed method and toolbox workflow only for the single-ended primary-inductor converter (SEPIC), for brevity. The nominal values of the SEPIC converter parameters used for this set of numerical simulations are presented in Table 1, along with their tolerances.
Firstly, we proceed to obtain the robust controller used for path planning, using the methodology presented in [14]. As such, we need to linearize the quasi-linear input-affine nonlinear system ( Σ s e p i c ) from Equation (57) according to Equation (30). This high-power configuration of the converter is targeted for renewable energy sources, with the desired input/state/output operating specifications being: output voltage y 0 ( t ) = u R ( t ) = 400 [V], with a nominal external voltage source E 0 = 300 [V] and load input resistance R 0 = 80 [ Ω ] . An initial guess for the state vector equilibrium point values was x ˜ 0 = [ 300 , 10 , 300 , 10 , 400 ] , followed by μ ˜ 0 = 0.55 for the duty cycle command. Through numeric computation, the nominal nonlinear plant equilibrium point follows:
u 0 , x 0 , y 0 = 300 , 80 , 0.5788 , 300 , 6.8711 , 297.722 , 5 , 400 , 400 .
After calling the System.linearize() routine for the LTI plant model, the uncertainty model needs to be computed. For modelling the uncertainty of the electrical components and external influences, an input multiplicative model was considered, i.e., G ( s ) = G n ( s ) 1 + Δ ( s ) W u n c ( s ) , with | | Δ | | 1 . This model has been numerically computed in an automatic manner from the interface of input u 3 ( t ) μ ( t ) to output y 1 ( t ) u R ( t ) , with the additional tolerances E 290 , 310 [V] and R 75 , 85 [ Ω ] , based on N = 1000 Monte Carlo simulations. A comprehensive frequency range for this SEPIC configuration is ω ̲ = 10 2 , ω ¯ = 10 8 , with 300 equally distributed samples in logarithmic domain. For the particle swarm optimization algorithm used for computing the uncertainty bounds, a successful set of hyperparameters consists in the swarm size of 1500, initial swarm span of 10 3 , minimum neighbors fraction of 0.85 , and inertia range of 0.15 , 1.15 , for a complex pole and zero pair transfer function, resulting in:
W u n c ( s ) = K s 2 + 2 ζ z ω n z s + ω n z 2 s 2 + 2 ζ p ω n p s + ω n p 2 = 0.53201 · s 2 + 1289 s + 4.623 e + 5 s 2 + 158.6 s + 6.286 e + 7 .
This procedure and its validation is portrayed in Figure 7.
In minimal form, the linearized SEPIC plant family is spanned by fourth-order stable and proper systems, with four zeros, three of which are of nonminimum phase, from the control path to the load resistor voltage output. The nominal MISO transfer matrix model from u ( t ) Δ R ( t ) , Δ E ( t ) , Δ μ ( t ) T to y ( t ) Δ u R ( t ) is:
G n ( s ) = s + 8.003 e + 5 α ( s ) 5.925 e + 7 ( s 2 + 162.5 s + 7.202 e + 7 ) 0.021779 ( s + 346.5 ) ( s 2 + 179.7 s + 5.637 e + 7 ) 4.1368 ( s + 2.304 e + 4 ) ( s 2 717.4 s + 5.145 e + 7 ) ,
where α ( s ) = ( s 2 + 2673 s + 3.749 e + 7 ) ( s 2 + 1339 s + 6.493 e + 7 ) .
Besides the classical scope of obtaining good transient and steady-state responses, another main purpose for the controller synthesis was to penalize near-resonance control effort, as the obtained controllers would otherwise be difficult to implement in practice, requiring small sampling periods, i.e., T s < 50 [ μ s]. As such, the sensitivity, complementary sensitivity and control effort loop-shaping filters are:
W S ( s ) = 0.5 s + 200 s + 2 , W T ( s ) = s 2 + 4000 s + 4 e + 6 1 e 4 s 2 + 56.57 s + 8 e + 6 , W K S ( s ) = 100 s + 346.5 s + 346.5 ,
with an imposed sensitivity bandwidth of at least ω B = 200 [rad/s], admissible steady-state error of maximum A = 10 2 , allowed sensitivity peak at M = 2 . Lower bandwidth was imposed for damping the frequency response resonance and, also, to move away from the physical limitations of the multiple nonminimum phase zeros. The complementary sensitivity bandwidth must be faster than ω B T = 2000 [rad/s], forcing a high frequency attenuation of at least A T = 10 4 with a roll-off order n = 2 , with an allowed peak of M T = 2 . The dynamic control effort weighting function imposes all frequencies with magnitude less than 1, with a further penalization of maximum 0.1 above ω 1000 [rad/s], due to the resonant peak centered at ω 9000 [rad/s].
From the direct application of the μ -synthesis procedure for the augmented plant with uncertainties, a controller of order n = 19 is obtained. This procedure is followed by a balanced order reduction operation, with the smallest-order controller which manages to ensure all imposed performance specifications, with a peak value μ Δ ( LLFT ( P , K ) ) 0.901513 < 1 , being of order n = 3 :
K r o b S E P I C = 1.997 3.257 3.706 0.5324 3.257 2294 6474 0.434 3.702 6474 1.539 e + 04 0.4941 i n e 0.5324 0.434 0.4941 0 .
The robust path-planning design of the weighting filters, along with frequency response closed-loop performance plots with the reduced-order controller for the uncertain SEPIC converter plant set are illustrated in Figure 8 and Figure 9. The closed-loop control system presents itself with large stability margins, with a phase margin γ k 80 [ ] and gain margin of m k d B 19 [dB]. Additionally, as specified by the parameters n = 2 and A T = 10 4 of W T ( s ) , the closed-loop control system practically nullifies stochastic sensor noise signals spanning from ω B T > 2000 [rad/s], using an initial roll-off of 40 [dB/dec], followed by an attenuation of at least four orders of magnitude. In practice, the attenuation continues with a further 20 [dB/dec] roll-off afterwards.
From the frequency response of Figure 8, the bandwidth is measured at ω B 390 [rad/s], resulting in an equivalent rise time less than t r 2.5 [ms], with a negligible steady-state error ε s s 10 2 × ( y s s y 0 ) = 0.2 [V], where y s s represents the steady-state value, relative to the desired equilibrium output y 0 = 400 [V]. Additionally, the system presents no overshoot in the vicinity of the operating point, as it was designed to behave similar to a first-order low-pass filter.
The SEPIC converter is a highly-nonlinear plant with respect to the command signal μ ( t ) , and there are use cases where the operating point may be dynamic, along with being generally affected by a diverse set of perturbations. Besides the previously-mentioned advantages of the already designed robust controller, it will be used as an auxiliary path planning component along with the K-PBC controller used to guarantee asymptotic stability for the entire domain of operation for the converter.
In order to design the K-PBC controller, we need to construct the output port-variable h K ( x , u ) such that the SEPIC converter is Krasovskii passive. The LMI problem presented in Claim 3, having the bounds of the load resistance R [ 10 , 1000 ] [ Ω ] , has been solved using the method described in [13], with a possible solution from the feasibility cone:
Q = 1 e 03 · 0.000714 0 0 0 0 0 0.514 0 0 0 0 0 0.00094 0 0 0 0 0 0.342 0 0 0 0 0 0.000714 ,
from which the output port-variable have the form presented in Claim 3, as well.
After this extra output is designed, the K-PBC having the form:
( Σ K P B C ) : x ˙ c = K 1 ( K 2 ( x c μ 🟉 ) h K ( x ) ) y c = x c ,
having two inputs: the output port-variable h K ( x ) of the extended ( Σ sepic ) and the desired input trajectory μ 🟉 , while the output is the actual value of the duty-cycle μ ( t ) y c ( t ) . The parameters of the controller (Equation (70)) are K 1 = 3 × 10 5 and K 2 = 10 8 .
Figure 9 and Figure 10 illustrate time domain simulations of the SEPIC converter in closed-loop configuration using the proposed control structure of cascaded K-PBC with μ -synthesis path planning K r o b at various operating points, as in Figure 3, along with different types of disturbances applied, in comparison to using the same robust controller-only approach, using the structure from Figure 4. The former compares the response of the proposed closed-loop system by comparison with the robust controller only counterpart, while the latter figure illustrates the robust stability and robust performance of the proposed method for a set of 50 Monte Carlo simulations sampled from the tolerance set of Table 1.
Both experiments start from null initial conditions applied to all state variables, in order to emphasize the behaviour to sudden jumps in the dynamics, followed by a succession of disturbances applied at the converter inputs, such as:
  • a sequence of load resistance steps:
    -
    from R ( t 1 = 0.035 ) = 80 [ Ω ] to R ( t 1 + = 0.035 + ) = 160 [ Ω ];
    -
    from R ( t 2 = 0.065 + ) = 160 [ Ω ] to R ( t 2 + = 0.065 + ) = 60 [ Ω ];
    -
    from R ( t 3 = 0.115 + ) = 60 [ Ω ] to R ( t 3 + = 0.115 + ) = 260 [ Ω ].
  • a ramp disturbance on the external source voltage from E ( t < 0.090 ) 300 [V], starting at t 4 = 0.090 [s] with an increasing slope of 500 and saturated at 50 [V];
  • a sequence of output reference steps:
    -
    from u R * ( t 5 = 0.025 ) = 400 [V] to u R * ( t 5 + = 0.025 + ) = 550 [V];
    -
    from u R * ( t 6 = 0.050 ) = 550 [V] to u R * ( t 6 + = 0.050 + ) = 250 [V];
    -
    from u R * ( t 7 = 0.075 ) = 250 [V] to u R * ( t 7 + = 0.075 + ) = 550 [V];
    -
    from u R * ( t 8 = 0.135 ) = 550 [V] to u R * ( t 8 + = 0.135 + ) = 450 [V].
As noticeable in Figure 9, the proposed method not only tracks the desired voltage reference at all operating points, but it also considerably improves transients caused by changes in the disturbance signals compared to using K r o b only, such as for the moments t = 0, t 1 , t 2 , t 3 , with smaller overshoots and more damped oscillations. A small compromise is that it adds insignificant overshoots when reference changes occur, such as at time t 6 . The steady-state performance is not affected by the addition of the K-PBC.
In addition to the previous results, Figure 10 shows the robustness of the method when subjected to the parametric uncertainties inherent to the SEPIC converter circuit, in which all the dynamic and steady-state performance indicators remain fundamentally unchanged even after ± 20 % variations of the circuit’s main component values.

6. Conclusions

The current paper presents a mathematical framework which allows to construct an output vector such that an input-affine nonlinear system is Krasovskii passive. Moreover, Theorem 3, which is a convex particularization of Theorem 1, without any dependency on the state variable x , provides the set of necessary and sufficient conditions for a quasi-linear affine system to be Krasovskii passive in terms of LMIs. To illustrate the flexibility and usefulness of the method, as a case study, a unified treatment for DC-DC converters is presented. After this interface is set for such a system, a method to construct a K-PBC is presented. The LLFT interconnection between the nonlinear system and K-PBC ensures the asymptotic stability, which means that the closed-loop system manages to follow the given input trajectory. However, the input trajectory has the same physical significance as the input of the nonlinear system, which, in general, differs from that of the output. As such, another component, called path planner, is mandatory in order to obtained the desired tracking performance. The proposed method includes a dynamical path planning as a robust controller computed for the linearized plant around the desired equilibrium point by solving the mixed-sensitivity μ -synthesis loop-shaping problem. In brief, the decentralized controller manages to ensure the asymptotic stability due to the K-PBC component, while the robust controller manages to compute the input trajectory such that the closed-loop system fulfills the robust performance.
Although Section 4 presents the possibility to compute the output port-variable such that various DC-DC converter topologies are Krasovskii passive, in Section 5 we present the SEPIC converter as a case study, due to its highly-nonlinear behaviour. As shown in the previous section, there are some important improvements if we compare the results obtained with the proposed method against the results obtained with the robust controller only. Moreover, the robustness of the proposed method has been proved using Monte Carlo simulations in time domain.
In this paper we managed to develop the mathematical background, which was successfully implemented in a second version of our toolbox, initially presented in [14]. However, in this current iteration we assume that all signals required for the output port-variable construction are available. But, in practice, an estimator will be necessary. As such, one possible extension for practical implementation could consist in (1) considering a high-gain observer for quasi-linear affine systems. Additionally, (2) there are three degrees of freedom for the K-PBC controller—the matrices Q, K 1 and K 2 —which may also be introduced in the optimization problem when computing the robust controller. Finally, (3) the implementation of the proposed method on microcontroller units will be studied, by also taking into consideration quantization effects.

Author Contributions

Conceptualization, V.M. and M.Ş.; methodology, V.M.; software, M.Ş. and V.M.; validation, P.D. and M.Ş.; formal analysis, P.D.; investigation, V.M. and M.Ş.; resources, M.Ş. and V.M.; data curation, M.Ş. and V.M.; writing—original draft preparation, V.M. and M.Ş.; writing—review and editing, M.Ş., V.M. and P.D.; visualization, M.Ş.; supervision, P.D.; project administration, V.M.; funding acquisition, P.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This paper was financially supported by the Project “Entrepreneurial Competences and Excellence Research in Doctoral and Postdoctoral Programs—ANTREDOC”, project co-funded by the European Social Fund financing agreement no. 56437/24.07.2019.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AREAlgebraic Riccati Equation
ARIAlgebraic Riccati Inequality
DC-DCDirect Current to Direct Current
K-PBCKrasovskii Passivity-Based Controller
LLFTLower Linear Fractional Transformation
LMILinear Matrix Inequality
LTILinear Time-Invariant
MiLModel-in-the-Loop
MIMOMultiple Input and Multiple Output
MISOMultiple Input and Single Output
PBCPassivity-Based Controller
PIDProportional Integral Derivative
PWMPulse-Width Modulation
SEPICSingle-Ended Primary Inductor Converter
SIMOSingle Input and Multiple Output
SISOSingle Input and Single Output
ULFTUpper Linear Fractional Transformation

References

  1. Willems, J.C. Dissipative dynamical systems part II: Linear systems with quadratic supply rates. Arch. Ration. Mech. Anal. 1972, 45, 352–393. [Google Scholar] [CrossRef]
  2. Khalil, H. Nonlinear Systems; Prentice Hall: Prentice, NJ, USA, 1996. [Google Scholar]
  3. Forni, F.; Sepulchre, R. On differential dissipative dynamical systems. IFAC Proc. Vol. 2013, 46, 15–20. [Google Scholar] [CrossRef] [Green Version]
  4. van der Schaft, A.J. On differential passivity. IFAC Proc. Vol. 2013, 46, 21–25. [Google Scholar] [CrossRef] [Green Version]
  5. Forni, D.; Sepulchre, R.; van der Schaft, A.J. On differential passivity of physical systems. In Proceedings of the 52nd IEEE Conference on Decision and Control, Firenze, Italy, 10–13 December 2013; pp. 6580–6585. [Google Scholar]
  6. Kosaraju, K.C.; Kawano, Y.; Scherpen, J.M.A. Krasovskii’s Passivity. In Proceedings of the 11th IFAC Symposium on Nonlinear Control Systems NOLCOS, Vienna, Austria, 4–6 September 2019; Volume 52, pp. 466–471. [Google Scholar]
  7. Doyle, J.; Glover, K.; Khargonekar, P.; Francis, B.A. State-space solutions to standard H 2 and H control problems. IEEE Trans. Autom. Control 1989, 34, 831–847. [Google Scholar] [CrossRef]
  8. Liu, K.Z.; Yao, Y. Robust Control—Theory and Applications; John Wiley & Sons: Singapore, 2016. [Google Scholar]
  9. Ionescu, V.; Oară, C.; Weiss, M. Generalized Riccati Theory and Robust Control—A Popov Function Approach; John Wiley & Sons: Chichester, UK, 1999. [Google Scholar]
  10. Şuşcă, M. Solving Algebraic Riccati Equations Using Proper Deflating Subspaces for H 2 / H Synthesis. Master’s Thesis, Technical University of Cluj-Napoca, Cluj-Napoca, Romania, 2019. [Google Scholar]
  11. Şuşcă, M.; Mihaly, V.; Stănese, M.; Dobra, P. Iterative Refinement Procedure for Solutions to Algebraic Riccati Equations. In Proceedings of the 2020 IEEE International Conference on Automation, Quality and Testing, Robotics (AQTR), Cluj-Napoca, Romania, 21–23 May 2020. [Google Scholar]
  12. Gahinet, P.; Apkarian, P. A linear matrix inequality approach to H control. Int. J. Robust Nonlinear Control 1994, 4, 421–448. [Google Scholar] [CrossRef]
  13. Mihaly, V. General Purpose Linear Matrix Inequality Solver with Applications in Robust and Nonlinear Control. Master’s Thesis, Technical University of Cluj-Napoca, Cluj-Napoca, Romania, 2020. [Google Scholar]
  14. Şuşcă, M.; Mihaly, V.; Stănese, M.; Morar, D.; Dobra, P. Unified CACSD Toolbox for Hybrid Simulation and Robust Controller Synthesis with Applications in DC-to-DC Power Converter Control. Mathematics 2021, 9, 731. [Google Scholar] [CrossRef]
  15. Iovine, A.; Mazenc, F. Bounded Control for DC/DC Converters: Application to Renewable Sources. In Proceedings of the 57th IEEE Conference on Decision and Control (CDC), Miami Beach, FL, USA, 17–19 December 2018; pp. 3415–3420. [Google Scholar]
  16. Cavallo, A.; Canciello, G.; Russo, A. Buck-Boost Converter Control for Constant Power Loads in Aeronautical Applications. In Proceedings of the 57th IEEE Conference on Decision and Control (CDC), Miami Beach, FL, USA, 17–19 December 2018; pp. 6741–6747. [Google Scholar]
  17. Cavallo, A.; Canciello, G.; Guida, B.; Kulsangcharoen, P.; Yeoh, S.S.; Rashed, M.; Bozhko, S. Multi-Objective Supervisory Control for DC/DC Converters in Advanced Aeronautic Applications. Energies 2018, 11, 3216. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, M.; Tang, F.; Wu, X.; Niu, J.; Zhang, Y.; Wang, J. A Nonlinear Control Strategy for DC-DC Converter with Unknown Constant Power Load Using Damping and Interconnection Injecting. Energies 2021, 14, 3031. [Google Scholar] [CrossRef]
  19. Mihaly, V.; Stănese, M.; Şuşcă, M.; Dobra, P. Interior Point Methods for Renewable Energy Management. In Proceedings of the 2020 IEEE International Conference on Automation, Quality and Testing, Robotics (AQTR), Cluj-Napoca, Romania, 21–23 May 2020. [Google Scholar]
  20. Middlebrook, R.D.; Cuk, S. A General Unified Approach to Modeling Switching-Converter Power Stages. Int. J. Electron. 1976, 42, 73–86. [Google Scholar]
  21. Czarkowski, D.; Kazimierczuk, M.K. Energy-conservation approach to modeling PWM DC-DC converters. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 1059–1063. [Google Scholar] [CrossRef]
  22. Sira-Ramirez, H.; Perez-Moreno, R.A.; Ortega, R.; Garcia-Esteban, M. Passivity-Based Controllers for the Stabilization of DC-to-DC Power Converters. Automatica 1997, 33, 499–513. [Google Scholar] [CrossRef]
  23. Mihaly, V.; Şuşcă, M.; Dobra, P. Passivity-Based Controller for Nonideal DC-to-DC Boost Converter. In Proceedings of the 22nd International Conference on Control Systems and Computer Science (CSCS), Bucharest, Romania, 28–30 May 2019; pp. 30–35. [Google Scholar]
  24. van der Schaft, A.J. Observability and controllability for smooth nonlinear systems. SIAM J. Control Optim. 1982, 20, 338–354. [Google Scholar] [CrossRef]
  25. Packard, A.; Doyle, J.; Balas, G. Linear Multivariable Robust Control with a μ Perspective. J. Dyn. Syst. Meas. Control 1993, 115, 426–438. [Google Scholar] [CrossRef]
  26. Mihaly, V.; Şuşcă, M.; Morar, D.; Stănese, M.; Dobra, P. μ-Synthesis for Fractional-Order Robust Controllers. Mathematics 2021, 9, 911. [Google Scholar] [CrossRef]
Figure 1. Passivity interconnection.
Figure 1. Passivity interconnection.
Energies 14 05571 g001
Figure 2. Generalized plant with uncertainties.
Figure 2. Generalized plant with uncertainties.
Energies 14 05571 g002
Figure 3. Proposed closed-loop control structure, comprised in extended affine plant Σ e , Krasovskii passivity-based controller K K P B C for nonlinear asymptotic stability and loop-shaping controller K r o b for robust performance.
Figure 3. Proposed closed-loop control structure, comprised in extended affine plant Σ e , Krasovskii passivity-based controller K K P B C for nonlinear asymptotic stability and loop-shaping controller K r o b for robust performance.
Energies 14 05571 g003
Figure 4. Robust controller synthesis diagram for the plant G, with the uncertainty set Δ , linearized at the operating point ( u 0 , x 0 , y 0 ) ; by convention, control inputs Δ u c of the linearized plant G are indexed after disturbance inputs Δ d .
Figure 4. Robust controller synthesis diagram for the plant G, with the uncertainty set Δ , linearized at the operating point ( u 0 , x 0 , y 0 ) ; by convention, control inputs Δ u c of the linearized plant G are indexed after disturbance inputs Δ d .
Energies 14 05571 g004
Figure 5. Class diagram for general-purpose system implementations, i.e., nonlinear, LTI, linearized, hybrid, along with the uncertain plant factory class, system interconnections, extended Krasovskii passivity framework and associated functionalities.
Figure 5. Class diagram for general-purpose system implementations, i.e., nonlinear, LTI, linearized, hybrid, along with the uncertain plant factory class, system interconnections, extended Krasovskii passivity framework and associated functionalities.
Energies 14 05571 g005
Figure 6. Nonideal DC-DC converters electrical circuits [14].
Figure 6. Nonideal DC-DC converters electrical circuits [14].
Energies 14 05571 g006
Figure 7. Computation procedure of the SEPIC input multiplicative uncertainty set, along with open loop responses for the equilibrium point having E = 300 [V], R = 80 [ Ω ] , U R = 400 [V]: magnitude fitting procedure and cost functional values, along with set of Monte Carlo simulations by sampling plant models starting from the initial nonlinear uncertainty set.
Figure 7. Computation procedure of the SEPIC input multiplicative uncertainty set, along with open loop responses for the equilibrium point having E = 300 [V], R = 80 [ Ω ] , U R = 400 [V]: magnitude fitting procedure and cost functional values, along with set of Monte Carlo simulations by sampling plant models starting from the initial nonlinear uncertainty set.
Energies 14 05571 g007
Figure 8. Open-loop SEPIC converter models, along with the S, T, K S functions using the reduced-order robust controller, which ensures imposed specifications, provides high stability margins and guarantees closed-loop response practically similar to a first-order low-pass filter’s response.
Figure 8. Open-loop SEPIC converter models, along with the S, T, K S functions using the reduced-order robust controller, which ensures imposed specifications, provides high stability margins and guarantees closed-loop response practically similar to a first-order low-pass filter’s response.
Energies 14 05571 g008
Figure 9. SEPIC averaged state-space closed-loop simulations using the robust controller K r o b with the structure from Figure 4 vs. the proposed control structure from Figure 3: subplot 1—time-varying disturbance inputs E ( t ) and R ( t ) ; subplot 2—output voltages u R ( t ) ; subplot 3—inductor currents i L 1 ( t ) and i L 2 ( t ) ; subplot 4—command signals μ ( t ) .
Figure 9. SEPIC averaged state-space closed-loop simulations using the robust controller K r o b with the structure from Figure 4 vs. the proposed control structure from Figure 3: subplot 1—time-varying disturbance inputs E ( t ) and R ( t ) ; subplot 2—output voltages u R ( t ) ; subplot 3—inductor currents i L 1 ( t ) and i L 2 ( t ) ; subplot 4—command signals μ ( t ) .
Energies 14 05571 g009
Figure 10. SEPIC-averaged state-space closed-loop simulations using the proposed control structure from Figure 3, illustrating 50 Monte Carlo simulations by sampling plant models from the tolerance set of Table 1: subplot 1—time-varying disturbance inputs E ( t ) and R ( t ) ; subplot 2—output voltages u R ( t ) ; subplot 3—command signals μ ( t ) .
Figure 10. SEPIC-averaged state-space closed-loop simulations using the proposed control structure from Figure 3, illustrating 50 Monte Carlo simulations by sampling plant models from the tolerance set of Table 1: subplot 1—time-varying disturbance inputs E ( t ) and R ( t ) ; subplot 2—output voltages u R ( t ) ; subplot 3—command signals μ ( t ) .
Energies 14 05571 g010
Table 1. SEPIC topology parameters, values, along with corresponding tolerances.
Table 1. SEPIC topology parameters, values, along with corresponding tolerances.
Param.Val.Tol.Param.Val.Tol.
L 1 2.57 [mH] ± 20 % L 2 1.71 [mH] ± 20 %
r L 1 130 [m Ω ] ± 10 % r L 2 110 [m Ω ] ± 10 %
r D S 1 0.01 [ Ω ] ± 10 % r D S 2 80 [m Ω ] ± 10 %
C 1 4.7 [ μ F] ± 20 % C 2 3.57 [ μ F] ± 20 %
r C 1 270 [m Ω ] ± 10 % r C 2 350 [m Ω ] ± 10 %
C i n 3.57 [ μ F] ± 20 % r C i n 270 [m Ω ] ± 10 %
V F 1 0.2 [V] ± 10 % V F 2 0.62 [V] ± 10 %
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mihaly, V.; Şuşcă, M.; Dobra, P. Krasovskii Passivity and μ-Synthesis Controller Design for Quasi-Linear Affine Systems. Energies 2021, 14, 5571. https://doi.org/10.3390/en14175571

AMA Style

Mihaly V, Şuşcă M, Dobra P. Krasovskii Passivity and μ-Synthesis Controller Design for Quasi-Linear Affine Systems. Energies. 2021; 14(17):5571. https://doi.org/10.3390/en14175571

Chicago/Turabian Style

Mihaly, Vlad, Mircea Şuşcă, and Petru Dobra. 2021. "Krasovskii Passivity and μ-Synthesis Controller Design for Quasi-Linear Affine Systems" Energies 14, no. 17: 5571. https://doi.org/10.3390/en14175571

APA Style

Mihaly, V., Şuşcă, M., & Dobra, P. (2021). Krasovskii Passivity and μ-Synthesis Controller Design for Quasi-Linear Affine Systems. Energies, 14(17), 5571. https://doi.org/10.3390/en14175571

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