Preprint
Article

Control of Three-Phase Two-Level Inverters: A Stochastic Lpv Model Approach

Submitted:

22 October 2024

Posted:

22 October 2024

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
This paper aims to propose a stochastic linear parameter varying (LPV) model approach to design a state feedback controller for three-phase two-level inverters. To deal with the parameter changes, stochastic noise and delays faced by the inverter, it is modeled as a stochastic LPV system with time delay. The stability analysis and control synthesis are conducted for the LPV system. With parameter-dependent Lyapunov functionals, the sufficient stability condition for asymptotical mean-square stability is obtained. In addition, the slack matrix technique is employed to improve the feasibility and reduce the conservatism of the conditions. The obtained theoretical results are applied to the three-phase two-level inverter, whose currents are treated as state variables and controlled to reach the equilibrium point. The simulation results validate the effectiveness of the proposed theories and demonstrate the advantages of using the slack matrix.
Keywords: 
Subject: 
Engineering  -   Electrical and Electronic Engineering

1. Introduction

Three-phase two-level inverters are the most widely used power converters in power electronics and motor drives. The main components of the inverters are controllable power switches such as IGBTs, MOSFETs, GTOs, etc. [1]. By implementing appropriate control strategies, they can output active and reactive power as required, meanwhile operate with desirable performance in terms of good current quality, high power efficiency, and strong robustness against disturbances. In real applications, the inverter faces uncertain conditions including time-varying grid voltage, filter capacitance and inductance changes, measurement delays, and environment noises. To effectively control the inverter, all these uncertainties should be taken into account [2,3].
Linear parameter-varying (LPV) systems are used to describe systems with time-varying parameters. Since the seminal work of Shamma in 1991 [4], the LPV control and filtering approach has attracted a lot of attention from the control field. This is because 1 ) it is effective to cope with nonlinearities and time-varying dynamics of the system, and 2 ) the maturely developed LTI control methodologies, such as sensitivity shaping and modeling tools, can be extended to LPV systems. During the past decades, several advances in theoretical and practical aspects have been done. Among others, theoretical works are related to analysis and synthesis issues [5,6], robust filtering problems [7,8], fault detection and isolation [9], etc. Practical engineering works are also popular, for instance, an LPV static output feedback control method is proposed to improve lateral stability and driving comfort of narrow tilting vehicles [10], LPV modeling method is adopted to ensure large signal stability for multi-mode buck-boost converters [11], the uncertain robot system is controlled by an LPV controller with sliding mode optimization [12]. Time delays exist in many practical systems, such as network controlled systems, electronic networks, hydraulic systems, and chemical processes, which are mostly nonlinear systems. Time delays may cause instability and oscillation in these systems and severely degrade performance [13,14]. To target this problem, research on time-delayed LPV system starts coming on stage in recent years, see [15,16].
Stochastic systems are used to model systems with stochastic perturbations. This approach is used in several fields and engineering applications, such as mechanical systems, economic systems, etc. Over the past decades, the study of stabilization and filtering problems for stochastic systems has received much attention, and many valuable results have been obtained. For stabilization problems, the asymptotic stability of semiswitched stochastic systems has been investigated in [17], the problem of input-to-state stability for nonlinear systems with stochastic impulses is systematically studied [18]. For filtering problems, the recursive filtering problem for a class of uncertain stochastic systems with amplify and forward relays has been investigated in [19], the problem of distributed robust filtering for switched stochastic time-delay systems with fading measurements over sensor networks is addressed in [20].
This paper aims to propose a stochastic LPV model approach to deal with the uncertain conditions and time delay faced by a three-phase two-level inverter, where the grid voltage, filter capacitance, and inductance are treated as time-varying parameters, and the environment noises as stochastic disturbances. Sufficient conditions for stability analysis and controller synthesis are obtained for the stochastic LPV system. The theoretical results obtained are applied to the inverter, whose currents are treated as state variables and controlled to reach the equilibrium point. The contributions of this paper are of three folds: 1) the LPV approach is used to model the three-phase two-level inverter, which is a new practical method; 2) comparing with previous theoretical works, this paper further considers the stochastic perturbation and parameter-varying time delay, based on which the stochastic LPV model is established and the stability conditions are obtained; 3) the slack matrix technique is adopted to improve the feasibility of the stability conditions, which facilitates the controller design process. The remaining of this work includes system description and inverter modeling in Section II. Main theoretical analysis is presented in Section III. Simulations to validate the proposal are provided in Secition IV. Finally, conclusioms are addressed in Section V.

2. System Description and Inverter Modeling

  • Figure 1 is a two-level three-phase inverter, which is normally used as grid-connecting power converter (the case in this paper) or motor drive. In ideal case, the inverter can be modelled in synchronous reference frame as follows [21]:
L d i d d t = u d V i n r i d ω L i q v d ,
L d i q d t = u q V i n r i q + ω L i d v q ,
where V i n is the input dc power source, ω is the frequency of grid voltage, L is the filtering inductor and r is its equivalent series resistance, and the dq variables are transformed from the three-phase abc variables, i.e., v d , v q are grid voltages transformed from v a , v b and v c , i d ,   i q are grid currents transformed from i a , i b and i c , and u d ,   u q are switching functions generated by controller, which is to be used to generate the switching signals S a , S b and S c , though PWM modulator. To provide required power to grid, the dq currents must be controlled to desired value. Taking into account the inverter parameter variations, stochastic grid perturbations, and control/communication delays, the inverter model can be formulated as a general time-delayed stochastically perturbed LPV system:
d x ( t ) = [ A ( ρ ( t ) ) x ( t ) + A h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) ] d t + B ( ρ ( t ) ) u ( t ) d t + B v ( ρ ( t ) ) v ( t ) d t + B ω ( ρ ( t ) ) x ( t ) d ω ( t ) ,
y ( t ) = C ( ρ ( t ) ) x ( t ) + C h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) + D ( ρ ( t ) ) v ( t ) ,
x ( θ ) = ϕ ( θ ) , θ [ h ( ρ ( 0 ) ) , 0 ] ,
where x ( t ) = [ i d ; i q ] is the state vector, A ( ρ ( t ) ) is the system matrix related to L, r, ω ; u ( t ) = [ u d V i n v d ; u q V i n v q ] is the control input, B ( ρ ( t ) ) is the control matrix related to L; v ( t ) R q belonging to L 2 [ 0 , ) , is either a disturbance input or a reference signal; ω ( t ) is a one-dimensional (1-D) Brownian motion satisfying E { d ω ( t ) } = 0 and E { d ω 2 ( t ) } = d t ; ρ ( t ) R s is vector-valued parameter evolving continuously over time and its range is limited to a compact subset; h ( ρ ( t ) ) is the parameter varying delay and satisfies:
0 h ( ρ ( t ) ) H < , | h ˙ ( ρ ( t ) ) | σ < ,
where H and σ are constant scalars; x ( θ ) is the initial data function given in the time interval [ H , 0 ] . The value of ρ ( t ) is unknown, but can be measured in real time.
Disregarding the control input u ( t ) , the autonomous system of (2) is formulated as:
d x ( t ) = [ A ( ρ ( t ) ) x ( t ) + A h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) ] d t + B v ( ρ ( t ) ) v ( t ) d t + B ω ( ρ ( t ) ) x ( t ) d ω ( t ) ,
y ( t ) = C ( ρ ( t ) ) x ( t ) + C h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) + D ( ρ ( t ) ) v ( t ) ,
x ( θ ) = ϕ ( θ ) , θ [ h ( ρ ( 0 ) ) , 0 ] .
On the other hand, to carry out stability analysis and synthesis, the exogenous disturbance “ v ( t ) ” in (2) is considered to be zero, yielding following system:
d x ( t ) = A ( ρ ( t ) ) x ( t ) + A h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) + B ( ρ ( t ) ) u ( t ) d t + B ω ( ρ ( t ) ) x ( t ) d ω ( t ) ,
x ( θ ) = ϕ ( θ ) , θ [ h ( ρ ( 0 ) ) , 0 ] .
Therefore the autonomous system of (4) can be formulated as
d x ( t ) = [ A ( ρ ( t ) ) x ( t ) + A h ( ρ ( t ) ) x ( t h ( ρ ( t ) ) ) ] d t + B ω ( ρ ( t ) ) x ( t ) d ω ( t ) ,
x ( θ ) = ϕ ( θ ) , θ [ h ( ρ ( 0 ) ) , 0 ] .
For simplicity, ρ will be used instead of ρ ( t ) in the rest of this paper. Before going further, some definitions and lemmas should be given, which are important for deriving main results.
Definition 1: For any initial state x ( 0 ) R n , the time-delayed LPV stochastic system (3) is said to be robustly stable with disturbance attenuation index γ if for all v ( t ) L 2 [ 0 , ) ] and all parameter trajectories, it holds that
y ( t ) 2 < γ v ( t ) 2 .
Also, the system is asymptotically mean square stable according to Definition 1 in [22].

3. System Analysis and Controller Design

3.1. Stability Analysis and Synthesis

Proposition 1: If there exists a family of parameter-dependent continuous differentiable symmetric positive matrices P ( ρ ) and Q ( ρ ) , such that for all the parameter trajectories, it satisfies
Π 11 P ( ρ ) A h ( ρ ) B ω T ( ρ ) P ( ρ ) Π 12 0 P ( ρ ) < 0 ,
where Π 11 = i = 1 s τ i P ( ρ ) ρ i + sym P ( ρ ) A ( ρ ) + Q ( ρ ) , Π 12 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) , then system (5) is asymptotically mean square stable.
Proof. Choose Lyapunov-Krasovskii functional as:
V ( x , ρ ) = x T ( t ) P ( ρ ) x ( t ) + t h ( ρ ) t x T ( s ) Q ( s ) x ( s ) d s .
According to Itô’s formula [23],
d V ( x , ρ ) = L V ( x , ρ ) d t + 2 x T ( t ) P ( ρ ) B ω ( ρ ) x ( t ) d ω ( t ) ,
where
L V ( x , ρ ) = x T ( t ) P ˙ ( ρ ) x ( t ) + 2 x T ( t ) P ( ρ ) A ( ρ ) x ( t ) + 2 x T ( t ) P ( ρ ) A h ( ρ ) x ( t h ( ρ ) ) + x T ( t ) B ω T ( ρ ) P ( ρ ) B ω ( ρ ) x ( t ) + x T ( t ) Q ( ρ ) x ( t ) 1 i = 1 s τ i h ( ρ ) ρ i · x T ( t h ( ρ ) ) Q ( t h ( ρ ) ) x ( t h ( ρ ) ) = x ( t ) x ( t h ( ρ ) ) T Π x ( t ) x ( t h ( ρ ) ) ,
with Π = Π 11 P ( ρ ) A h ( ρ ) 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) , and Π 11 = i = 1 s τ i P ( ρ ) ρ i + sym P ( ρ ) A ( ρ ) + B ω T ( ρ ) P ( ρ ) B ω ( ρ ) + Q ( ρ ) . Thus, according to [24], system (5) is asymptotically mean square stable if Π < 0 , which is equal to the inequality (7) with Schur complement transformation. Proof is completed.
Proposition 2: If there exists a family of parameter-dependent continuous differentiable symmetric positive matrices R ( ρ ) , Q ¯ ( ρ ) and matrices F ( ρ ) , such that for all the parameter trajectories, it holds that
Π 11 A h ( ρ ) R ( t h ( ρ ) ) R ( ρ ) B ω T ( ρ ) Π 22 0 R ( ρ ) < 0 ,
where
Π 11 = sym A ( ρ ) R ( ρ ) + B ( ρ ) F ( ρ ) i = 1 s τ i R ( ρ ) ρ i + Q ¯ ( ρ ) , Π 22 = 1 i = 1 s τ i h ( ρ ) ρ i Q ¯ ( t h ( ρ ) ) ,
then system (4) is mean square stable under u ( t ) = K ( ρ ) x ( t ) with K ( ρ ) = F ( ρ ) R 1 ( ρ ) .
Proof. Consider system (4) with u ( t ) = K ( ρ ) x ( t ) , K ( ρ ) is a parameter-dependent state feedback controller.
Substitute A ( ρ ) in (7) with A ( ρ ) + B ( ρ ) K ( ρ ) , and perform a congruence transformation to it with diag { P 1 ( ρ ) , P 1 ( t h ( ρ ) ) , P 1 ( ρ ) } , it is obtained that
Π ¯ 11 A h ( ρ ) P 1 ( t h ( ρ ) ) P 1 ( ρ ) B ω T ( ρ ) Π ¯ 12 0 P 1 ( ρ ) < 0 ,
where
Π ¯ 11 = P 1 ( ρ ) i = 1 s τ i P ( ρ ) ρ i P 1 ( ρ ) + P 1 ( ρ ) Q ( ρ ) P 1 ( ρ ) + sym { ( A ( ρ ) + B ( ρ ) K ( ρ ) ) P 1 ( ρ ) } , Π ¯ 12 = 1 i = 1 s τ i h ( ρ ) ρ i · P 1 ( t h ( ρ ) ) Q ( t h ( ρ ) ) P 1 ( t h ( ρ ) ) .
Let matrix R ( ρ ) = P 1 ( ρ ) , Q ¯ ( ρ ) = P 1 ( ρ ) Q ( ρ ) P 1 ( ρ ) , F ( ρ ) = K ( ρ ) P 1 ( ρ ) , then it is obtained that (10), and K ( ρ ) = F ( ρ ) R 1 ( ρ ) . Proof is completed.
Remark 1: τ i denotes the varying rate of parameter i which is assumed to be measurable in real time. If τ i is not measurable, but its bound is known a prior, i.e. | τ i | υ i , then i = 1 s τ i R ( ρ ) ρ i can be approximated by i = 1 s ± υ i R ( ρ ) ρ i to obtain a new stability condition. i = 1 s ± ( · ) represents the sum of every combination of + ( · ) and ( · ) , which contains a total of 2 s combinations.
Proposition 3: If there exists a family of parameter-dependent continuous differentiable symmetric positive matrices P ( ρ ) , Q ( ρ ) and matrix W, such that for all the parameter trajectories, it satisfies that
W + W T P ρ + W T A ρ W T A h ρ 0 W T Π 22 0 B ω T ρ 0 Π 33 0 0 P 1 ρ 0 P ρ < 0 ,
where Π 22 = P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) , Π 33 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) , then system (5) is asymptotically mean square stable.
Proof. The inequality (12) can be written as
0 P ρ 0 0 0 Π 22 0 B ω T ρ 0 Π 33 0 0 P 1 ρ 0 P ρ + I A T ρ A h T ρ 0 I W I 0 0 0 0 T + < 0 .
The null spaces of [ I A T ( ρ ) A h T ( ρ ) 0 I ] and [ I 0 0 0 0 ] are
A ( ρ ) A h ( ρ ) 0 I I 0 0 0 0 I 0 0 0 0 I 0 0 0 0 I and 0 0 0 0 I 0 0 0 0 I 0 0 0 0 I 0 0 0 0 I .
According to Projection Lemma [25], following inequalities are obtained:
Π ¯ 11 P ( ρ ) A h ( ρ ) B ω T ( ρ ) P ( ρ ) Π ¯ 12 0 0 P 1 ( ρ ) 0 P ( ρ ) < 0 ,
i = 1 s τ i P ρ ρ i + Q ρ P ρ 0 B ω T ρ 0 1 i = 1 s τ i h ρ ρ i Q t h ρ 0 0 P 1 ρ 0 P ρ < 0 ,
where
Π ¯ 11 = i = 1 s τ i P ( ρ ) ρ i + sym P ( ρ ) A ( ρ ) + Q ( ρ ) P ( ρ ) , Π ¯ 12 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) .
The inequality (14) equals to the inequality (7) with Schur complement transformation, which indicates (12) can ensure the system to be asymptotically mean square stable. Proof is completed.
Remark 2: With the introduction of a new additional matrix W, the matrix P ( ρ ) in the Lyapunov function is separated from the system matrix. This extra degree of freedom reduces the usually strong interrelations between plant data and Lyapunov variables, hence improves the solution feasibility and reduces the conservatism.
Theorem 1: If there exists a family of parameter-dependent continuous differentiable symmetric positive matrices X ( ρ ) , Y ( ρ ) , P ¯ ( ρ ) and matrices R ( ρ ) , V, such that for all the parameter trajectories, it holds that
V + V T Π 12 A h ρ V 0 V Π 22 0 V T B ω T ρ 0 Π 33 0 0 P ¯ ρ 0 X ρ < 0 ,
where
Π 12 = X ( ρ ) + A ( ρ ) V + B ( ρ ) R ( ρ ) , Π 22 = X ( ρ ) + i = 1 s τ i X ( ρ ) ρ i + Y ( ρ ) , Π 33 = 1 i = 1 s τ i h ( ρ ) ρ i Y ( t h ( ρ ) ) ,
then system (4) is mean square stable under u = K ( ρ ) x ( t ) with K ( ρ ) = R ( ρ ) V 1 ( ρ ) .
Proof. Consider system (4) with u ( t ) = K ( ρ ) x ( t ) , substituting A ( ρ ) in (12) with A ( ρ ) + B ( ρ ) K ( ρ ) , and perform a congruence transformation with diag { W 1 , W 1 , W 1 , I , W 1 } , it is obtained that
Π ¯ 11 Π ¯ 12 A h ( ρ ) W 1 0 W 1 Π ¯ 22 0 W T B ω T ( ρ ) 0 Π ¯ 33 0 0 P 1 ( ρ ) 0 W T P ( ρ ) W 1 < 0 ,
where
Π ¯ 11 = ( W 1 + W T ) , Π ¯ 12 = W T P ( ρ ) W T + ( A ( ρ ) + B ( ρ ) K ( ρ ) ) W 1 , Π ¯ 22 = W T P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) W 1 , Π ¯ 33 = 1 i = 1 s τ i h ( ρ ) ρ i W T Q ( t h ( ρ ) ) W 1 .
Let matrix X ( ρ ) = W T P ( ρ ) W 1 , Y ( ρ ) = W T Q ( ρ ) W 1 , V = W 1 , P ¯ ( ρ ) = P 1 ( ρ ) , (16) is obtained, and K ( ρ ) = R ( ρ ) V 1 . Proof is completed.

3.2. Performance Analysis and Synthesis

Proposition 4: The time-delayed LPV stochastic system (3) is said to be robustly stable with disturbance attenuation γ if there exists a family of parameter-dependent continuous differentiable symmetric positive matrices R ( ρ ) and Q ( ρ ) , such that for all the parameter trajectories, it holds that
Π 11 P ( ρ ) A h ( ρ ) P ( ρ ) B v ( ρ ) C T ( ρ ) B ω T ( ρ ) P ( ρ ) Π 22 0 C h T ( ρ ) 0 γ 2 I D T ( ρ ) 0 I 0 P ( ρ ) < 0 .
where
Π 11 = P ( ρ ) A ( ρ ) + A T ( ρ ) P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) , Π 22 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) .
Proof. Assume that the system has zero initial state, i.e., x ( t ) = 0 when t [ h ( ρ ) , 0 ] , then according to Itô’s formula, it can be obtained that
E V ( x ( t ) , t ) = E 0 t L V ( x ( s ) , s ) d s .
Let
J ( t ) = E 0 t y T ( s ) y ( s ) γ 2 v T ( s ) v ( s ) d s = E 0 t x ( s ) x ( s h ( ρ ) ) v ( s ) T Π ¯ x ( s ) x ( s h ( ρ ) ) v ( s ) d s ,
where
Π ¯ = Π ¯ 11 P ( ρ ) A h ( ρ ) P ( ρ ) B v ( ρ ) Π ¯ 22 0 γ 2 I + C T ( ρ ) C h T ( ρ ) D T ( ρ ) C T ( ρ ) C h T ( ρ ) D T ( ρ ) T ,
with
Π ¯ 11 = P ( ρ ) A ( ρ ) + A T ( ρ ) P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) + B ω T ( ρ ) P ( ρ ) B ω ( ρ ) , Π ¯ 22 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) .
With Schur complement, condition (18) ensures Π ¯ < 0 , thus J ( t ) < 0 , meaning y T ( t ) y ( t ) < γ 2 v T ( t ) v ( t ) . Thus system (3) is robustly stable in the sense of definition 2. Proof is completed.
Proposition 5: The time-delayed LPV stochastic system (2) is said to be robustly stabilized by state feedback controller K ( ρ ) with disturbance attenuation γ if there exists a family of parameter-dependent continuous differentiable symmetric positive matrices R ( ρ ) , Q ¯ ( ρ ) and matrices F ( ρ ) , such that for all the parameter trajectories, it holds that
Π 11 A h R ( t h ( ρ ) ) B v ( ρ ) R ( ρ ) C T R ( ρ ) B ω T ( ρ ) Π 22 0 R ( t h ( ρ ) ) C h T 0 γ 2 I D T 0 I 0 R ( ρ ) < 0 ,
where
Π 11 = sym A ( ρ ) R ( ρ ) + B ( ρ ) F ( ρ ) i = 1 s τ i R ( ρ ) ρ i + Q ¯ ( ρ ) , Π 22 = 1 i = 1 s τ i h ( ρ ) ρ i Q ¯ ( t h ( ρ ) ) .
and K ( ρ ) = F ( ρ ) R 1 ( ρ ) . Proof. Consider system (2) with u ( t ) = K ( ρ ) x ( t ) , substituting A ( ρ ) in (18) with A ( ρ ) + B ( ρ ) K ( ρ ) , and perform a congruence transformation with diag { P 1 ( ρ ) , P 1 ( t h ( ρ ) ) , I , I , P 1 ( ρ ) } , it is obtained that
Π ¯ 11 A h P 1 ( t h ( ρ ) ) B v ( ρ ) P 1 ( ρ ) C T P 1 ( ρ ) B ω T ( ρ ) Π ¯ 22 0 P 1 ( t h ( ρ ) ) C h T 0 γ 2 I D T 0 I 0 P 1 ( ρ ) < 0 ,
where
Π ¯ 11 = sym ( A ( ρ ) + B ( ρ ) K ( ρ ) ) P 1 ( ρ ) + P 1 ( ρ ) Q ( ρ ) + i = 1 s τ i P ( ρ ) ρ i P 1 ( ρ ) , Π ¯ 22 = 1 i = 1 s τ i h ( ρ ) ρ i · P 1 ( t h ( ρ ) ) Q ( t h ( ρ ) ) P 1 ( t h ( ρ ) ) .
Let matrix R ( ρ ) = P 1 ( ρ ) , Q ¯ ( ρ ) = P 1 ( ρ ) Q ( ρ ) P 1 ( ρ ) , F ( ρ ) = K ( ρ ) P 1 ( ρ ) , then (21) is obtained, and K ( ρ ) = F ( ρ ) R 1 ( ρ ) . Proof is completed.
Proposition 6: The time-delayed LPV stochastic system (3) is said to be robustly stable with disturbance attenuation γ if there exists a family of parameter-dependent continuous differentiable symmetric positive matrices P ( ρ ) , Q ( ρ ) and matrix W, such that for all the parameter trajectories, it holds that
Π 11 Π 12 W T A h ( ρ ) W T B v ( ρ ) 0 0 W T Π 22 0 0 C T ( ρ ) B ω T ( ρ ) 0 Π 33 0 C h T ( ρ ) 0 0 γ 2 I D T ( ρ ) 0 0 I 0 0 P 1 ( ρ ) 0 P ( ρ ) < 0 ,
where Π 11 = ( W + W T ) , Π 12 = P ( ρ ) + W T A ( ρ ) , Π 22 = P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) , Π 33 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) .
Proof. The inequality (23) can be written as:
0 P ( ρ ) 0 0 0 0 0 Π 22 0 0 C T ( ρ ) B ω T ( ρ ) 0 Π 33 0 C h T ( ρ ) 0 0 γ 2 I D T ( ρ ) 0 0 I 0 0 P 1 ( ρ ) 0 P ( ρ ) + I A T ( ρ ) A h T ( ρ ) B v T ( ρ ) 0 0 I W ( ρ ) I 0 0 0 0 0 0 T + I A T ( ρ ) A h T ( ρ ) B v T ( ρ ) 0 0 I T W T ( ρ ) I 0 0 0 0 0 0 < 0 .
The null spaces of
[ I A ( ρ ) A h ( ρ ) B v ( ρ ) 0 0 I ] and [ I 0 0 0 0 0 0 ] are
A ( ρ ) A h ( ρ ) B v ( ρ ) 0 0 I I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I , 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I 0 0 0 0 0 0 I respectively. According to Projection Lemma, (24) is obtained,
Π ¯ 11 P ( ρ ) A h ( ρ ) P ( ρ ) B v ( ρ ) C T ( ρ ) B ω T ( ρ ) P ( ρ ) Π ¯ 12 0 C h T ( ρ ) 0 0 γ 2 I D T ( ρ ) 0 0 I 0 0 P 1 ( ρ ) 0 P ( ρ ) < 0 ,
with
Π ¯ 11 = i = 1 s τ i P ( ρ ) ρ i + P ( ρ ) A ( ρ ) + A T ( ρ ) P ( ρ ) + Q ( ρ ) P ( ρ ) , Π ¯ 12 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) ,
and
Π ˜ 11 0 0 C T ( ρ ) B ω T ( ρ ) 0 Π ˜ 22 0 C h T ( ρ ) 0 0 γ 2 I D T ( ρ ) 0 0 I 0 0 P 1 ( ρ ) 0 P ( ρ ) < 0 ,
with Π ˜ 11 = i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) P ( ρ ) and Π ˜ 22 = 1 i = 1 s τ i h ( ρ ) ρ i Q ( t h ( ρ ) ) .
The inequality (24) equals to the inequality (18) with Schur complement transformation, which indicates (23) can ensure the system to be robustly stable. Proof is completed.
Theorem 2: The time-delayed LPV stochastic system (2) is said to be robustly stabilized by state feedback controller K ( ρ ) with disturbance attenuation γ if there exists a family of parameter-dependent continuous differentiable symmetric positive matrices X ( ρ ) , Y ( ρ ) , P ¯ ( ρ ) and matrices R ( ρ ) , V, such that for all the parameter trajectories, it holds that
Π 11 Π 12 Π 13 B v ( ρ ) 0 0 V Π 22 0 0 V T C T ( ρ ) V T B ω T ( ρ ) 0 Π 33 0 V T C h T ( ρ ) 0 0 γ 2 I V T D T ( ρ ) 0 0 I 0 0 P ¯ ( ρ ) 0 X ( ρ ) < 0
where
Π 11 = ( V + V T ) , Π 12 = X ( ρ ) + A ( ρ ) V + B ( ρ ) R ( ρ ) , Π 13 = A h ( ρ ) V , Π 22 = X ( ρ ) + i = 1 s τ i X ( ρ ) ρ i + Y ( ρ ) , Π 33 = 1 i = 1 s τ i h ( ρ ) ρ i Y ( t h ( ρ ) ) .
and K ( ρ ) = R ( ρ ) V 1 ( ρ ) .
Proof. Consider system (2) with u ( t ) = K ( ρ ) x ( t ) , substituting A ( ρ ) in (23) with A ( ρ ) + B ( ρ ) K ( ρ ) , and perform a congruence transformation with diag { W 1 , W 1 , W 1 , I , I , I , W 1 } , it is obtained that
Π ¯ 11 Π ¯ 12 Π ¯ 13 B v ( ρ ) W 1 0 0 W 1 Π ¯ 22 0 0 W T C T ( ρ ) 0 0 Π ¯ 33 0 W T C h T ( ρ ) 0 0 γ 2 I W T D T ( ρ ) 0 0 I 0 0 P 1 ( ρ ) 0 Π ¯ 77 < 0 ,
where
Π ¯ 11 = ( W 1 + W T ) , Π ¯ 12 = W T P ( ρ ) W T + ( A ( ρ ) + B ( ρ ) K ( ρ ) ) W 1 , Π ¯ 13 = A h ( ρ ) W 1 , Π ¯ 22 = W T P ( ρ ) + i = 1 s τ i P ( ρ ) ρ i + Q ( ρ ) W 1 , Π ¯ 33 = 1 i = 1 s τ i h ( ρ ) ρ i W T Q ( t h ( ρ ) ) W 1 , Π ¯ 77 = W T P ( ρ ) W 1 .
Let matrix X ( ρ ) = W T P ( ρ ) W 1 , Y ( ρ ) = W T Q ( ρ ) W 1 , V = W 1 , P ¯ ( ρ ) = P 1 ( ρ ) , (26) is obtained, and K ( ρ ) = R ( ρ ) V 1 . Proof is completed.

4. Simulation Results

For following simulation, the grid and inverter parameters are chosen as follows: v d = 400 V, v d = 0 V, V i n = 750 V, L = 0.002 H, r = 0.004 Ω , ω = 314 rad/s. These values are chosen to comply with European low voltage grid, which has a voltage of 400V, a frequency of 50Hz, with small line resistance and inductance. The matrices of the LPV inverter model (2) are as follows:
A ( ρ ) = r L + 0.1 ρ 1 ( t ) ω ω r L + 0.1 ρ 1 ( t ) , A h ( ρ ) = r L + 0.1 ρ 2 ( t ) ω ω r L + 0.1 ρ 2 ( t ) , B ( ρ ) = 1 L + 0.1 ρ 1 ( t ) 1 L + 0.1 ρ 1 ( t ) , B υ ( ρ ) = 0.2 + 0.1 ρ 1 ( t ) 0.1 + 0.1 ρ 2 ( t ) , B ω ( ρ ) = 1.4 + 0.1 ρ 1 ( t ) 1.0 + 0.2 ρ 1 ( t ) 0.7 0.9 ,
C ( ρ ) = 1.2 + 0.2 ρ 1 ( t ) + 0.1 ρ 2 ( t ) 0.8 T , C h ( ρ ) = 0.6 + 0.1 ρ 1 ( t ) + 0.1 ρ 2 ( t ) 0.4 , D ( ρ ) = 1 + 0.1 ρ 1 ( t ) + 0.1 ρ 2 ( t ) , h ( ρ ) = 0.01 ρ 2 ( t ) ,
where ρ 1 ( t ) = sin ( t ) and ρ 2 ( t ) = | cos 5 t | are time-varying parameters, satisfying ρ 1 ( t ) [ 1 , 1 ] , ρ 2 ( t ) [ 0 , 1 ] , τ 1 [ 1 , 1 ] , and τ 2 [ 5 , 5 ] ; h ( ρ ) is the parameter varying delay with a bounded varying rate. The initial states are set as x ( t ) = [ 1.0 , 1.0 ] T , where t [ 1 , 0 ] .
To deal with the problem of infinite LMIs to be solved during controller synthesis, gridding technique and basis functions should be adopted. The basis functions are chosen as f 1 ( ρ ) = 1 , f 2 ( ρ ) = ρ 1 ( t ) , f 3 ( ρ ) = ρ 2 ( t ) , leading to R ( ρ ) = j = 1 3 f j ( ρ ) R j , F ( ρ ) = j = 1 3 f j ( ρ ) F j , therefore the controller synthesis problem becomes finding matrices R j ( j = 1 , 2 , 3 ) and F j ( j = 1 , 2 , 3 ) to satisfy (10), (16), (21), and (26).

4.1. Stabilization

The state feedback control K ( ρ ) is designed to ensure that system (4) is asymptotically mean square stable. By solving the conditions in (10) in Proposition (2), it is obtained that,
R 1 = 10 26 0.1402 0.0557 0.0557 0.1030 ,   R 2 = 10 28 0.3142 0.2191 0.2191 0.1663 , R 3 = 10 28 0.5181 0 . 5006 0 . 5006 0.5802 , F 1 = 10 26 0.1251 0.1250 , F 2 = 10 28 0.1031 0.0825 , F 3 = 10 28 0.4335 0.4369 . K 1 = 0.5222 0.9310 ,   K 2 = 0.2204 0.7866 , K 3 = 9.3997 8.8628 .
Figure 2-(a) shows the state response of one arbitrarily chosen set, which means a random set of currents can be controlled to a certain reference. Figure 2-(b) shows the state response of ten arbitrarily chosen sets, which shows that the states of all the sets converge to zero. As can be seen, the stochastic phenomenon can influence the trajectories of currents, but does not affect the stability. Therefore, the time-delayed LPV inverter system is asymptotically stabilized under the stochastic perturbations.

4.2. Stabilization with Slack Matrix

Using slack matrix technique and solve the conditions in (16) in Theorem 1, it is obtained that:
V = 10 28 0.2004 0.1983 ; 0.2192 0.2194 , R 1 = 10 26 0.2288 0.2284 , R 2 = 10 28 0.1545 0.1342 , R 3 = 10 30 0.2733 0.2741 . K 1 = 28.4143 78.3968 ,   K 2 = 8.8494 7 . 3862 , K 3 = 0.0018 0.0141 .
Similar as in the previous subsection, Figure 2-(c) shows the state response of one arbitrarily chosen set, and Figure 2-(d) shows the state response of ten arbitrarily chosen sets. As can be seen, with the controller obtained by adopting slack matrix, the currents converge to the equilibrium with faster speed, compared to the controller obtained in previous subsection, which demonstrates the advantage of slack matrix technique.

4.3. Disturbance Attenuation

Let’s consider in the system a disturbance signal v ( t ) = 1 / ( 1 + 0.5 t 2 ) . Then, the state feedback control K ( ρ ) is designed to ensure the system (2) is stabilized and achieves some disturbance attenuation index γ . By solving the conditions in (21) in Theorem 2, it is obtained that,
R 1 = 0.0846 0.0846 0.0846 0.0847 ,   R 2 = 0.0033 0.0035 0 . 0035 0.0037 , R 3 = 0.0111 0.0111 0.0111 0.0111 ,   F 1 = 0.3018 0.3017 , F 2 = 0.0023 0.0019 ,   F 3 = 0.0031 0.0031 . K 1 = 3.8245 0.2580 ,   K 2 = 65.2790 62.5179 , K 3 = 1.3953 1.6769 .
With this controller, the currents converge to equilibrium and the system is stabilized under stochastic perturbations, as shown in Figure 3-(a) and Figure 3-(b), and the achieved disturbance attenuation index is γ is 21.4113 .

4.4. Disturbance Attenuation with Slack Matrix

Similar to Section 4-C, let’s consider in the system a disturbance signal v ( t ) = 1 / ( 1 + 0.5 t 2 ) . By solving the conditions in (26) in Proposition (5), it is obtained that,
V = 0.1136 0.1150 ; 0.1167 0.1185 , R 1 = 0.1501 0.1500 ,   R 2 = 0.0075 0.0076 , R 3 = 0.0012 0.0012 , K 1 = 5.8506 4.4128 ,   K 2 = 0.0117 0..0758 , K 3 = 0.0369 0.0260 .
Applying the controller to the inverter system, as shown in Figure 3-(c) and Figure 3-(d), the currents converge at faster speed, comparing to the controller obtained without adopting slack matrix technique. Furthermore, the achieved disturbance attenuation index γ is 12.0578 , which is smaller than the previous method. These results again demonstrate the advantage of slack matrix technique.

5. Conclusions

A stochastic LPV approach has been proposed to control the three-phase two-level inverter. The time-delayed stochastic LPV system is first established, for which the analysis and synthesis of the system are carried out. With the parameter-dependent Lyapunov functionals, sufficient conditions have been proposed for stability analysis and control synthesis in terms of parameter-dependent LMIs. The slack matrix approach has been used to derive a new set of stability LMIs, which improves the feasibility and reduces the conservatism of the solution. The theoretical results are verified on the inverter, whose currents are controlled to reach the equilibrium point. Simulation results validate the effectiveness of the proposed theories and demonstrate the advantages of adopting a slack matrix, in terms of faster state response and lower disturbance attenuation index.

Author Contributions

Conceptualization, W. Luo and L. Wu; methodology, W. Luo; validation, W. Luo, R. Zhang and J. Zhang; formal analysis, W. Luo; investigation, W. Luo; data curation, R. Zhang and J. Zhang; writing—original draft preparation, W. Luo, R. Zhang and J. Zhang; writing—review and editing, W. Luo and S. Vazquez; supervision, L. G. Franquelo. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (62003114, 62033005, 62320106001), the China Postdoctoral Science Foundation (2020M681097), the Heilongjiang Postdoctoral Fund (LBH-Z20134), the Natural Science Foundation of Heilongjiang Province (ZD2021F001), the Fundamental Research Funds for the Central Universities (Grant No. HIT.NSFJG202207).

References

  1. W. Luo, S. Stynski, A. Chub, L. G. Franquelo, M. Malinowski and D. Vinnikov, “Utility-Scale Energy Storage Systems: A Comprehensive Review of Their Applications, Challenges, and Future Directions,” in IEEE Industrial Electronics Magazine, vol. 15, no. 4, pp. 17–27, 2021. [CrossRef]
  2. Y. Yin, S. Vazquez, A. Marquez, J. Liu, J. I. Leon, L. Wu, and L. G. Franquelo, “Observer-based sliding-mode control for grid-connected power converters under unbalanced grid conditions,” IEEE Transactions on Industrial Electronics, vol. 69, no. 1, pp. 517–527, 2022. [CrossRef]
  3. E. Zafra, J. Granado, V. B. Lecuyer, S. Vazquez, A. M. Alcaide, J. I. Leon, and L. G. Franquelo, “Computationally Efficient Sphere Decoding Algorithm Based on Artificial Neural Networks for Long-Horizon FCS-MPC,” IEEE Transactions on Industrial Electronics, vol. 71, no. 1, pp. 39–48, 2024. [CrossRef]
  4. J. S. Shamma and M. Athans, “Guaranteed properties of gain scheduled control for linear parameter-varying plants,” Automatica, vol. 27, no. 3, pp. 559–564, 1991. [CrossRef]
  5. C. Briat, “Stability analysis and state-feedback control of LPV systems with piecewise constant parameters subject to spontaneous poissonian jumps,” IEEE Control Systems Letters, vol. 2, no. 2, pp. 230–235, 2018. [CrossRef]
  6. L. O. de Araújo, P. C. Pellanda, J. F. Galdino, and A. M. Simões, “Haar-based stability analysis of LPV systems,” IEEE Transactions on Automatic Control, vol. 60, no. 1, pp. 192–198, 2015. [CrossRef]
  7. M. S. de Oliveira and R. L. Pereira, “LMI-based filter design conditions for discrete-time LPV systems with bounded parameter variation,” IEEE Transactions on Automatic Control, vol. 66, no. 2, pp. 910–915, 2021. [CrossRef]
  8. A. P. Pandey and M. C. de Oliveira, “On the necessity of LMI-based design conditions for discrete time LPV filters,” IEEE Transactions on Automatic Control, vol. 63, no. 9, pp. 3187–3188, 2018. [CrossRef]
  9. J. Tan, S. Olaru, F. Xu, and X. Wang, “Toward a convex design framework for online active fault diagnosis of LPV systems,” IEEE Transactions on Automatic Control, vol. 67, no. 8, pp. 4154–4161, 2022. [CrossRef]
  10. A.-T. Nguyen, P. Chevrel, and F. Claveau, “LPV static output feedback for constrained direct tilt control of narrow tilting vehicles,” IEEE Transactions on Control Systems Technology, vol. 28, no. 2, pp. 661–670, 2020. [CrossRef]
  11. J. Ma, M. Zhu, Y. Li, and X. Cai, “Dynamic analysis of multimode buck–boost converter: An LPV system model point of view,” IEEE Transactions on Power Electronics, vol. 36, no. 7, pp. 8539–8551, 2021. [CrossRef]
  12. G. P. Incremona, A. Ferrara, and V. I. Utkin, “Sliding mode optimization in robot dynamics with LPV controller design,” IEEE Control Systems Letters, vol. 6, pp. 1760–1765, 2022. [CrossRef]
  13. Y. Shen, Z.-G. Wu, P. Shi, Z. Shu, and H. R. Karimi, “H control of Markov jump time-delay systems under asynchronous controller and quantizer,” Automatica, vol. 99, pp. 352–360, 2019. [CrossRef]
  14. L. Wu, X. Su, and P. Shi, “Sliding mode control with bounded L2 gain performance of Markovian jump singular time-delay systems,” Automatica, vol. 48, no. 8, pp. 1929–1933, 2012. [CrossRef]
  15. X. Yang, S. Yin, and O. Kaynak, “Robust identification of LPV time-delay system with randomly missing measurements,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, vol. 48, no. 12, pp. 2198–2208, 2018. [CrossRef]
  16. X. Yang and S. Yin, “Robust global identification and output estimation for LPV dual-rate systems subjected to random output time-delays,” IEEE Transactions on Industrial Informatics, vol. 13, no. 6, pp. 2876–2885, 2017. [CrossRef]
  17. Q. Wang, Baoand Zhu, “Stability analysis of semi-Markov switched stochastic systems,” Automatica, vol. 94, pp. 72–80, 2018. [CrossRef]
  18. Y. Tang, X. Wu, P. Shi, and F. Qian, “Input-to-state stability for nonlinear systems with stochastic impulses,” Automatica, vol. 113, p. 108766, 2020. [CrossRef]
  19. H. Tan, B. Shen, P. Kaixiang, and H. Liu, “Robust recursive filtering for uncertain stochastic systems with amplify-and-forward relays,” International Journal of Systems Science, vol. 51, no. 7, pp. 1188–1199, 2020. [CrossRef]
  20. Y. Chen, Z. Wang, Y. Yuan, and P. Date, “Distributed H filtering for switched stochastic delayed systems over sensor networks with fading measurements,” IEEE Transactions on Cybernetics, vol. 50, no. 1, pp. 2–14, 2020. [CrossRef]
  21. W. Luo, S. Vazquez, J. Liu, F. Gordillo, L. G. Franquelo, and L. Wu, “Control system design of a three-phase active front end using a sliding-mode observer,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, vol. 52, no. 2, pp. 739–748, 2022. [CrossRef]
  22. W. Luo, S. Li, S. Vazquez, J. Du, L. Wu, and L. G. Franquelo, “Grid-connected inverter control via linear parameter-varying system approach,” in IECON 2022 – 48th Annual Conference of the IEEE Industrial Electronics Society, 2022, pp. 1–6. [CrossRef]
  23. X. Mao, “Stability of stochastic differential equations with Markovian switching,” Stochastic Processes and Their Applications, vol. 79, no. 1, pp. 45–67, 1999. [CrossRef]
  24. S. Xu and T. Chen, “Robust H control for uncertain stochastic systems with state delay,” IEEE Transactions on Automatic Control, vol. 47, no. 12, pp. 2089–2094, 2002. [CrossRef]
  25. P. Gahinet and P. Apkarian, “A linear matrix inequality approach to H control,” International Journal of Robust and Nonlinear Control, vol. 4, no. 4, pp. 421–448, 1994. [CrossRef]
Figure 1. Two-level three-phase inverter.
Figure 1. Two-level three-phase inverter.
Preprints 121974 g001
Figure 2. Stabilization: (a) one set of state response without slack matrix; (b) ten sets of state response without slack matrix; (c) one set of state response with slack matrix; (d) ten sets of state response with slack matrix.
Figure 2. Stabilization: (a) one set of state response without slack matrix; (b) ten sets of state response without slack matrix; (c) one set of state response with slack matrix; (d) ten sets of state response with slack matrix.
Preprints 121974 g002
Figure 3. Disturbance attenuation: (a) one set of state response without slack matrix; (b) ten sets of state response without slack matrix; (c) one set of state response with slack matrix; (d) ten sets of state response with slack matrix.
Figure 3. Disturbance attenuation: (a) one set of state response without slack matrix; (b) ten sets of state response without slack matrix; (c) one set of state response with slack matrix; (d) ten sets of state response with slack matrix.
Preprints 121974 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Alerts
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated