Next Article in Journal
Real-Time Optimization of Wastewater Treatment Plants via Constraint Adaptation
Next Article in Special Issue
Delay-Dependent Stability of Impulsive Stochastic Systems with Multiple Delays
Previous Article in Journal
Changes in Volatile Compounds during Grape Brandy Production from ‘Cabernet Sauvignon’ and ‘Syrah’ Grape Varieties
Previous Article in Special Issue
Stochastic Review Inventory Systems with Deteriorating Items; A Steady-State Non-Linear Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Stabilization for Discrete-Time System with Input Delay and Multiplicative Noise in Control Variable

School of Engineering, Qufu Normal University, Rizhao 276800, China
*
Author to whom correspondence should be addressed.
Submission received: 1 May 2022 / Revised: 12 May 2022 / Accepted: 14 May 2022 / Published: 16 May 2022
(This article belongs to the Special Issue Advances in Nonlinear and Stochastic System Control)

Abstract

:
The stabilization problems for time-delay stochastic systems with multiplicative noise in the control variable are investigated in this paper. The innovative contributions are described as follows. Since the past work on stabilization is based on some delay-dependent algebraic Riccati equation (DARE), how to numerically calculate the stabilizing solution remains an unsolved and open problem. On the one hand, an iterative algorithm for computing the unique stabilizing solution of DARE is proposed, while the convergence property is also proved. On the other hand, the concepts of critical stabilization and essential destabilization are proposed as a supplement to stochastic stabilization in terms of spectrum technique. Moreover, the Lyapunov-based necessary and sufficient conditions are developed.

1. Introduction

It is becoming increasingly clear that areas such as network control, finance, power systems, robotics, and a large range of practical engineering systems can be described as stochastic models [1,2,3,4,5]. In the actual system, due to the limited capacity between different parts of the actual system, the uncertainty of the external environment changes, etc., time delay is a common phenomenon in system modeling. For the actual system in engineering, time delay will cause some negative impacts that cannot be ignored, such as the instability of the system and the decline in system performance. In recent years, stochastic systems with time delay have been extensively studied, especially for the issues such as LQ optimal control, stochastic stabilization, observability, mixed H 2 / H control and H 2 control. One of the crucial problems is the stabilization problem, which is a scheme for obtaining closed-loop stability control.
Over the past few decades, the research on stochastic systems with input delay has attracted widespread attention; see [1,2,6,7]. The vast majority of stability studies use the method of the generalized Lyapunov equation (GLE) [8,9,10]. This approach can provide easy-to-verify stability standards, but there is nothing that can be done to deal with some important issues. Recently, in an input-delay-free stochastic system, refs. [11,12,13], a random spectral analysis method has been adopted to study the stability problem for some linear time-invariant (LTI) models. Notice that the role of spectral analysis is irreplaceable by the Lyapunov functional approach. Based on the spectrum theory, the concepts of asymptotic mean square stabilization, critical stabilization, and essential destabilization are defined, which are effective and meticulous. However, we find that there are very few stabilization results for stochastic dynamics with input delay. And due to the interference of input delay factors, the delay-dependent Lyapunov operator and operator spectrum is much more complex than the generalized Lyapunov operator and its operator spectrum. How to utilize the available information to determine the necessary and sufficient conditions for critical stabilization and essential destabilization is very difficult. As far as we know, there are no relevant research results on time-delay stochastic systems. On the other hand, since the past work on stabilization is based on the DARE, the analysis for the nonlinear DARE is challenging. How to use the convergence algorithm to numerically calculate the unique stabilizing solution of DARE deserves in-depth study. This has greatly stimulated the motivation and confidence in our scientific research of this paper.
The objective of this paper is to obtain the necessary conditions for asymptotic mean square stabilization, while at the same time we define critical stabilization and essential destabilization as a supplement to the asymptotic mean square stabilization. We propose an iterative algorithm for computing the stabilizing solution and give the proof of convergence of the algorithm. In addition, we provide a set of necessary and sufficient conditions to confirm critical stabilization and essential destabilization.
To be specific, our research layout and methodology are described as follows. Section 2 formulates the systems under consideration and also presents preliminary results which will express our results more precisely. The main results are given in Section 3, Section 4 and Section 5 in terms of methods such as convergence algorithm, recursion, delay-dependent Lyapunov equation (DLE), resizing algorithm, etc. Section 6 provides two examples to show the effectiveness of our theoretic results.
Notations: k e r ( A ) means its kernel space and let I be the identity matrix with appropriate dimension. A B represents the Kronecker product of the matrices A and B. δ k s is a Kronecker function. Let S n be the space of all n-dimensional symmetric matrices. D [ 0 , 1 ] = { z | z C , | z 0 | 1 } . e k means a sequence of real random variables defined on the complete probability space ( Ω , F , P , F t ) with F t = σ { e s , s = 0 , , t } . Define x ^ t | s = E [ x t | F s ] to be the conditional expectation x t w.r.t. the filter F s .

2. Problem Formulation and Preliminaries

Consider the following discrete time stochastic dynamics
x k + 1 = A x k + B v k d + e k C v k d ,
where e k is assumed to be a scalar random white noise satisfying E ( e k ) = 0 and E ( e k e s ) = σ 2 δ k s , d > 0 is a positive number representing time delay in control input. To condense the formula, let us define [ A , B , C | d ] representing the system (1). Note that the mean square stabilization problem of stochastic system has been fully studied in [14,15,16,17] and the study in [18,19,20,21,22,23,24,25] also extended the theoretical results to a more general level.
Remark 1.
For a stochastic system with input delay and multiplicative noise in the control variable, it is worth mentioning that this stochastic dynamic model has a wide application in engineering practice. Specifically, consider an NCS operating over a reliable communication channel, in which the control signal is assumed to suffer both packet loss and network-induced delay from the controller to the actuator. As shown in [17], the overall NCS can be modelled by
x k + 1 = A x k + η k B v k d
where the packet loss process of γ k follows the Bernoulli distribution with P ( η k = 0 ) = s [ 0 , 1 ] . When we denote e k = η k ( 1 s ) , system (2) is equivalent to
x k + 1 = A x k + ( 1 s ) B v k d + e k B v k d
which is a special case of system [ A , B , C | d ] .
Facilitated by stochastic control techniques, the objective is to find the necessary conditions for the asymptotic mean square stabilization. More importantly, explore the equivalent conditions for the critical stabilization and essential destabilization of the system (1) under consideration. To express our results more precisely, we need to introduce the following definitions and lemmas. Different from the previous work, one significant contribution is that our control law is designed as the feedback of an extended state that contains the recently available state information and part values from previous control inputs.
Here we recall some basic definitions and lemmas related to stabilization in the mean square sense, which are necessary for the exposition in this paper. For additional details and other basics, we refer the reader to [14,26,27,28].
Definition 1.
System [ A , B , C | d ] is asymptotically mean square stabilizable, if there exists a control input with v k d = K x ^ k | k d 1 , k d , such that
x k + 1 = A x k + [ B + e k C ] K x ^ k | k d 1
is asymptotically mean square stable.
In order to make the description and research more convenient, we define the following operators. For system [ A , B , C | d ] , let F K [ A , B , C | d ] be a linear operator from S n to S n defined as follows
F K [ A , B , C | d ] ( M ) = ( A + B K ) M ( A + B K ) + σ 2 A d C K M K C ( A ) d .
The adjoint operator F K * of F K [ A , B , C | d ] from S n to S n is given as
F K * ( M ) = ( A + B K ) M ( A + B K ) + σ 2 K C ( A ) d M A d C K .
The spectral set of F is represented by
ρ ( F K [ A , B , C | d ] ) = { λ , F K [ A , B , C | d ] ( M ) = λ M , M S n , M 0 } .
Moreover, for system [ A , B , C | d ] , the nonlinear Riccati operators R [ A , B , C | d ] from S n to S n is defined as
R [ A , B , C | d ] ( M ) = A M A L U 1 L + Q ,
where Q 0 , R > 0 and
L = B M A , U = B M B + R + σ 2 C Q C + σ 2 k = 0 d 1 C ( A ) k + 1 Q A k + 1 C + σ 2 C ( A ) d ( M Q ) A d C .
Lemma 1.
([29]). The equivalent conditions for the mean square stabilization of system [ A , B , C | d ] are given as follows.
(a) 
System [ A , B , C | d ] is stabilizable in the mean square sense, if and only if there admit K and P > 0 satisfying the following equation
P = Q + ( A + B K ) P ( A + B K ) + σ 2 A d C K P K C ( A ) d , Q > 0 .
(b) 
System [ A , B , C | d ] is stabilizable in the mean square sense, if and only if there exists a constant matrix K such that | ρ ( F K [ A , B , C | d ] ) | < 1 .
(c) 
The mean square stabilization of system [ A , B , C | d ] is equivalent to that of the following delay free system
z k + 1 = A z k + B u k + e k A d C u k .
(d) 
The mean square stabilizable of system [ A , B , C | d ] is equivalent to the following DARE
P = A P A L U 1 L + Q , R > 0 , Q > 0 .
admits a unique stabilizing solution P s > 0 , where L and U satisfy (9).
By Lemma 1, the stabilization of [ A , B , C | d ] is equivalent to the existence and uniqueness of stabilizing solution to DARE. It is remarkable that in contrast to the standard algebraic Riccati equation (ARE) or modified ARE, how to calculate this stabilizing solution to the DARE remains an open problem. Until now, to the best of our knowledge, there have been no firm results to these related problems.

3. The Necessary Condition of Asymptotic Mean Square Stabilization

In this section, we propose an algorithm for numerically solving the unique stabilizing solution of DARE, which is based on the asymptotic properties of the stabilizing solution. To begin with, for system [ A , B , C | d ] , let G K [ A , B , C | d ] be a linear operator from S n to S n , and it is defined as
G K [ A , B , C | d ] ( M ) = ( A + B K ) M ( A + B K ) + K R K + Q + k = 0 d 1 σ 2 K C ( A ) k + 1 Q A k + 1 C K + σ 2 K C ( A ) d ( M Q ) A d C K + σ 2 K C Q C K ,
where Q > 0 , R > 0 . Then the linear operator G K [ A , B , C | d ] can be defined in terms of (6) as
G K [ A , B , C | d ] ( M ) = F K * ( M ) + Q ^ > F K * ( M ) ,
where
Q ^ = K R K + Q + k = 0 d 1 σ 2 K C ( A ) k + 1 Q A k + 1 C K σ 2 K C ( A ) d Q A d C K + σ 2 K C Q C K > 0 .
When the gain is defined by K M = U 1 L , (8) can be written as
R [ A , B , C | d ] ( M ) = A M A L U 1 L + Q = ( A + B K M ) M ( A + B K M ) + k = 0 d 1 σ 2 K M C ( A ) k + 1 Q A k + 1 C K M + σ 2 K M C ( A ) d ( M Q ) A d C K M + σ 2 K M C Q C K M + K M R K M + Q = G K M [ A , B , C | d ] ( M ) .
By utilizing the definitions of the delay-dependent algebraic Riccati operator R [ A , B , C | d ] and the linear operator G K [ A , B , C | d ] , we have the following lemmas. Since the proof is similar to that of Lemmas 3 and 4 in [30], the details are omitted here.
Lemma 2.
Suppose d > 0 , then one has the following statements.
(a) 
For any M 0 ,
R [ A , B , C | d ] ( M ) = min K G K [ A , B , C | d ] ( M ) .
(b) 
If M 2 M 1 0 , then,
R [ A , B , C | d ] ( M 2 ) R [ A , B , C | d ] ( M 1 ) .
(c) 
For system [ A , B , C | d ] , define X k + 1 = R [ A , B , C | d ] ( X k ) and Y k + 1 = R [ A , B , C | d ] ( Y k ) , then, for any k 0 ,
X 0 Y 0 0 X k Y k 0 .
Lemma 3.
Suppose there exists K and M ¯ > 0 , such that M ¯ > F K * ( M ¯ ) . For any M 0 , there is a limit satisfying lim k F K * ( F K * ( F K * k ( M ) ) ) = 0 .
Below, by exploiting Lemmas 2 and 3, the stabilizing solution of DARE can be obtained immediately.
Theorem 1.
For P k ( k 0 ) that satisfies P k + 1 = R [ A , B , C | d ] ( P k ) , where the initial value P 0 0 is arbitrary, it is bounded and converges to the unique stabilizing solution P s > 0 of DARE (11), if the stochastic system [ A , B , C | d ] is mean square stabilizable, where Q > 0 and R > 0 .
Proof. 
Suppose system [ A , B , C | d ] is stabilizable in the mean square sense where Q > 0 , R > 0 . It follows from ( d ) of Lemma 1 that there admits a unique stabilizing solution P s > 0 for DARE, which implies that P s = R [ A , B , C | d ] ( P s ) holds.
Define the gain K P s = U s 1 L s . By (12) and (13), one has
P s = R [ A , B , C | d ] ( P s ) = G K P s [ A , B , C | d ] ( P s ) = F K P s * ( P s ) + Q ^ s > F K P s * ( P s ) ,
where
Q ^ s = K P s R K P s + Q + k = 0 d 1 σ 2 K P s C ( A ) k + 1 Q A k + 1 C K P s σ 2 K P s C ( A ) d Q A d C K P s + σ 2 K P s C Q C K P s 0 .
Lemma 3 means that there exists K P s and P s > 0 , such that P s > F K P s * ( P s ) . In this case, we have
μ P s > F K P s * ( P s ) ,
where 0 < μ < 1 , which implies that
F K P s * ( F K P s * ( F K P s * k + 1 ( P s ) ) ) < μ F K P s * ( F K P s * ( F K P s * k ( P s ) ) ) < < μ k + 1 P s .
For any P 0 0 , according to (12)–(14), one has
P k + 1 = R [ A , B , C | d ] ( P k ) = G K P k [ A , B , C | d ] ( P k ) = min K G K [ A , B , C | d ] ( P k ) G K P s [ A , B , C | d ] ( P k ) = F K P s * ( P k ) + Q ^ s .
Similarly, we have
P k F K P s * ( P k 1 ) + Q ^ s , P k 1 F K P s * ( P k 2 ) + Q ^ s , P 1 F K P s * ( P 0 ) + Q ^ s .
Via (18), the inequality (17) leads to
P k + 1 F K P s * ( P k ) + Q ^ s F K P s * ( F K P s * ( P k 1 ) + Q ^ s ) + Q ^ s = F K P s * ( F K P s * ( P k 1 ) ) + F K P s * ( Q ^ s ) + Q ^ s F K P s * ( F K P s * ( F K P s * k + 1 ( P 0 ) ) ) + F K P s * ( F K P s * ( F K P s * k ( Q ^ s ) ) ) + F K P s * ( F K P s * ( F K P s * k 1 ( Q ^ s ) ) ) + + F K P s * ( Q ^ s ) + Q ^ s .
For P 0 0 , one has P 0 δ 0 P s where δ 0 is a normal number. Similarly, for Q ^ s 0 , one has Q ^ s δ Q s P s where δ Q s is a normal number. Therefore, we have
F K P s * ( F K P s * ( F K P s * k + 1 ( P 0 ) ) ) δ 0 F K P s * ( F K P s * ( F K P s * k + 1 ( P s ) ) ) < μ k + 1 δ 0 P s .
Similarly, for 1 i k , it follows that
F K P s * ( F K P s * ( F K P s * i ( Q ^ s ) ) ) < μ i δ Q s P s .
According to the method of inequality scaling and (19)–(21), it can be seen that { P k } is bounded.
To give the stabilizing solution of (11), set X 0 and Y 0 to be 0 and P 0 + P s , respectively. Then one has X 0 P 0 Y 0 .
Suppose X 0 = 0 , by (b) in Lemma 2, it is known that
X 1 = Q X 0 = 0 ,
which implies that
X k X k 1 X 1 X 0 .
Note it has been proved that { X k } is bounded. Then the sequence { X k } must have a limit, i.e.,
lim k X k = R [ A , B , C | d ] ( lim k X k 1 ) = X ¯ ,
where X ¯ is the limit of sequence { X k } . As as result, the above equivalent can be rewritten as
X ¯ = R [ A , B , C | d ] ( X ¯ ) ,
and the unique stabilizing solution P s > 0 of (11) is equal to its positive definite solution in [30]. Thus, X ¯ = P s .
Suppose Y 0 = P 0 + P s , by (b) and (c) in Lemma 2, it is known that
Y 1 = R [ A , B , C | d ] ( P 0 + P s ) P s ,
which implies that Y k P s 0 ( k 0 ) . So, together with (17) and Lemma 3, we have
Y k P s F K P s * ( Y k 1 P s ) F K P s * ( F K P s * ( F K P s * k ( Y 0 P s ) ) ) , lim k F K P s * ( F K P s * ( F K P s * k ( Y 0 P s ) ) ) = 0 .
That is, lim k Y k = P s is established.
By item (b) in Lemma 2, one has X k P k Y k , and it is known from the Squeeze Theorem
lim k P k = P s ,
which ends this proof. □

4. Critical Stabilization

Next, we introduce another concept of stochastic stabilization which will be called critical stabilization and use DLE to give necessary and sufficient conditions for judgment.
Definition 2.
System [ A , B , C | d ] is said to be critically stabilizable, if there is a feedback control v k d = K x ^ k | k d 1 , k d such that | ρ ( F K [ A , B , C | d ] ) | D [ 0 , 1 ] .
Now, we will use the resize method to prove the following theorem to judge critical stabilization. For the sake of discussion, we introduce the following stochastic dynamics notated as [ A , B , C | d , α ]
x k + 1 = ( A α ) x k + ( B α ) v k d + e k ( C α ) v k d .
Correspondingly, for system [ A , B , C | d , α ] , define F K [ A , B , C | d , α ] be a linear operator from S n to S n satisfying
F K [ A , B , C | d , α ] ( M ) = ( A α + B α K ) M ( A α + B α K ) + σ 2 ( A α ) d C α K M K C α ( A α ) d .
The spectral set of F is represented by
ρ ( F K [ A , B , C | d , α ] ) = { λ α , F K [ A , B , C | d , α ] ( M ) = λ α M , M S n , M 0 } .
Theorem 2.
The following statements are equivalent.
(i) 
System [ A , B , C | d ] is critical stabilization.
(ii) 
There exists matrix K, K S n , such that for any α > 1 and Q > 0 , the following DLE
P = Q + ( A α + B α K ) P ( A α + B α K ) + σ 2 ( A α ) d C α K P [ ( A α ) d C α K ] ,
admits a unique positive solution P S n , P > 0 .
(iii) 
There exists matrix K, K S n , such that for any α > 1 , the following inequality
P ( A α + B α K ) P ( A α + B α K ) σ 2 ( A α ) d C α K P [ ( A α ) d C α K ] > 0 ,
admits a positive solution P S n , P > 0 .
(iv) 
There exists matrices K, K S n , such that for any α > 1 and Q > 0 , the following inequality
P = Q + ( A α + B α K ) P ( A α + B α K ) + σ 2 K C α ( A α ) d P ( A α ) d C α K ,
admits a unique solution P S n , P > 0 .
Proof. 
Interestingly, (ii)⇔(iii), what we have to do is to prove (i)⇔(ii), (iii)⇔(iv). First, we give a simple proof to show (i)⇔(ii).
Sufficiency: From the conditions, it can be seen that (25) has a unique solution P S n , P > 0 which is equivalent to the asymptotically mean square stabilization of the following system
z k + 1 α = A α z k + B α u k + e k ( A α ) d C α u k .
From Theorem 2 in [29], it can be seen that the system (28) is asymptotically mean square stabilizable, which is equivalent to that system [ A , B , C | d , α ] is stabilizable. By applying (b) in Lemma 1, it is easy to know that | ρ ( F K [ A , B , C | d , α ] ) | < 1 . Then, let α 1 , with the continuity of spectrum, it follows that | ρ ( F K [ A , B , C | d ] ) | 1 , which completes the proof of Sufficiency part.
Necessity: Using the Kronecker product theory and the definition of the spectrum yields the following formula
λ X = [ ( A + B K ) ( A + B K ) + σ 2 A d C K ( A ) d C K ] X .
Compared with
λ α X = [ σ 2 ( A α ) d C α K ( A α ) d C α K + ( A α + B α K ) ( A α + B α K ) ] X ,
where X is a column monomial matrix by all elements of X, one obtains that the relationship between λ and λ α can be inferred as λ λ α , which means that | λ α | < 1 . Combining with (a) and (b) in Lemma 1, the necessity of the proof is proved.
Next, we prove that (iii) and (iv) are equivalent. From Theorem 1 in [31], it is easy to find that (iii) and (iv) are equivalent when system [ A , B , C | d , α ] is stabilizable. □
It is worth noting that the condition of critical stabilization is weaker than that of asymptotic mean square stabilization. According to the distribution of the spectrum of the operator on the complex plane, in the next section, we will propose the concept and related properties of the essential destabilization.

5. Essential Destabilization

In this section, we focus on an interesting concept of stochastic stabilization which will be called essential destabilization with the spectrum technique. More importantly, the Lyapunov-type necessary and sufficient conditions will be proposed.
We start this section by defining essential destabilization. The following definition and lemmas are necessary to establish our main theoretical results.
Definition 3.
System [ A , B , C | d ] is called to be essentially destabilizable, if for any feedback v k d = K x ^ k | k d 1 , the closed-loop dynamic system x k + 1 = A x k + [ B + e k C ] K x ^ k | k d 1 is essentially unstable, i.e., there exists λ i ρ ( F K [ A , B , C | d ] ) satisfying | λ i | > 1 .
For any X n = ( x i j ) n × n S n , define two column vectors X n and X n ˜ , where X n is with all elements of X n and X n ˜ formed by different elements of X n . T ( n 2 , n ( n + 1 ) 2 ) denotes the transform matrix from X n ˜ to X n .
As stated in [21], T ( n 2 , n ( n + 1 ) 2 ) has the following properties
  • Matrix T ( n 2 , n ( n + 1 ) 2 ) has column full rank.
  • Matrix T ( n 2 , n ( n + 1 ) 2 ) T ( n 2 , n ( n + 1 ) 2 ) is invertible.
Besides, let us define
( T ( n 2 , n ( n + 1 ) 2 ) , A , B , C | d ) = [ T ( n 2 , n ( n + 1 ) 2 ) T ( n 2 , n ( n + 1 ) 2 ) ] 1 T ( n 2 , n ( n + 1 ) 2 ) × [ ( A + B K ) ( A + B K ) + A d C K A d C K ) ] T ( n 2 , n ( n + 1 ) 2 ) .
Lemma 4.
([31]).
ρ ( F K [ A , B , C | d ] ) = ρ ( ( T ( n 2 , n ( n + 1 ) 2 ) , A , B , C | d ) ) .
Lemma 5.
Assume that λ i 1 holds for any λ i ρ ( F K [ A , B , C | d ] ) . Then for any given symmetric matrix Q S n , the following delay-dependent equation
( A + B K ) P ( A + B K ) + σ 2 ( A d C K ) P ( A d C K ) P = Q ,
has a unique solution P S n .
Proof. 
First, Let us define the following Jordan canonical form of ( T ( n 2 , n ( n + 1 ) 2 ) , A , B , C | d )
J = J 1 J 2 J s
with Jordan blocks
J i = λ i 1 1 λ i k i × k i
where i = 1 s k i = n ( n + 1 ) 2 . Then, take the basis B = { v 1 1 v k 1 1 , v 1 2 v k 2 2 , , v 1 s v k s s } of C n ( n + 1 ) 2 and define V = [ v 1 1 v k 1 1 , v 1 2 v k 2 2 , , v 1 s v k s s ] , where v 1 1 v k 1 1 , v 1 2 v k 2 2 , ⋯, v 1 s v k s s are linear and independent eigenvectors to λ i , i = 1 , 2 , s . It follows that
V 1 ( T ( n 2 , n ( n + 1 ) 2 ) , A , B , C | d ) V = J ,
or equivalently
( T ( n 2 , n ( n + 1 ) 2 ) , A , B , C | d ) V = V J .
Then it follows that X 1 1 X k 1 1 , X 1 2 X k 2 2 , ⋯, X 1 s X k s s consititute a complete basis B of S n . Moreover, it follows that
F K [ A , B , C | d ] ( X 1 i ) = λ i X 1 i , i = 1 , 2 , , s F K [ A , B , C | d ] ( X i j i ) = λ i X j i + X j i 1 i , j i = 2 , 3 , , k i
Now, we begin to show that the symmetric solution of (30) exists and is unique. Since B is a basis of S n , each matrix Q in S n can be uniquely expressed as
Q = [ c 1 , c 2 , c 3 c n ( n + 1 ) n ( n + 1 ) 2 2 ] x 1 1 x k 1 1 x k s s
Without losing generality, it is assumed that P follows the following structure
P = [ b 1 1 , , b k 1 1 , b 1 2 , , b k 2 2 , , b 1 s b k s s ] x 1 1 x k 1 1 x k s s = i = 1 s j i = 1 k i b j i i x j i i .
Due to the fact that (31) is linear w.r.t. P, when we take (33) and (34) into (30), we obtain that
i = 1 s j i = 1 k i b j i i F K [ A , B , C | d ] ( x j i i ) i = 1 s j i = 1 k i b j i i x j i i = [ c 1 , c 2 , c 3 c n ( n + 1 ) / 2 ] x 1 1 x k 1 1 x k s s
that is
[ b 1 1 , , b k 1 1 , b 1 2 , , b k 2 2 , , b 1 s b k s s ] J [ b 1 1 , , b k 1 1 , b 1 2 , , b k 2 2 , , b 1 s b k s s ] = [ c 1 , c 2 , c n ( n + 1 ) n ( n + 1 ) 2 2 ] ,
or equivalently
P J P I = Q .
By this condition, it is known that λ i 1 . Note that J I is invertible, then P exists and is unique. □
Theorem 3.
System [ A , B , C | d ] is essentially destabilizable, if and only if for any K and W > 0 , there exists a constant α > 1 , such that the following DLE
V + W = ( A α + B α K ) V ( A α + B α K ) + σ 2 ( A α ) d C α K V K C α ( A α ) d ,
admits a unique symmetric matrix V with at least one positive eigenvalue.
Proof. 
First, let’s generalize a further relationship between λ and λ α . Utilizing Kronecker matrix theorem, (7) can be rewritten as
[ ( A + B K ) ( A + B K ) + σ 2 A d C K A d C K ] X = λ X .
Compared with
[ ( A α + B α K ) ( A α + B α K ) + σ 2 ( A α ) d C α K ( A α ) d C α K ] X = λ α X ,
one obtains the following relationship
λ i α d + 1 < λ α , i < λ i α < λ i ,
where the λ i ρ ( F k [ A , B , C | d ] ) , λ α , i ρ ( F K [ A , B , C | d , α ] ) , and X is a column monomical matrix constructed by all elements of X.
Sufficiency: We assume that the system [ A , B , C | d ] is stabilizable. Then | ρ ( F K [ A , B , C | d , α ] ) | < 1 can be derived from item (b) in Lemma 1 and the former statements.
Compared with condition, we have
( V + W ) = ( A α + B α K ) ( V ) ( A α + B α K ) + σ 2 ( A α ) d C α K ( V ) K C α ( A α ) d .
We can conclude that the matrix V is positive definite, that is, V is a negative definite matrix, which contradicts the problem condition. Therefore, [ A , B , C | d , α ] is not asymptotically mean square stabilizable.
From the former result, we can find a λ α , i { λ α , | λ α | 1 } . Denote λ α , i * to be the maximum spectrum in λ α . With the condition α > 1 and the the relationship between λ α and λ , there exists a | λ | > 1 , which implies that the system is essentially destabilizable.
Necessity: Set β = min { | λ j | , λ j ρ ( F K [ A , B , C | d ] ) , | λ j | > 1 } . Clearly, we have β > 1 , since the possible value of λ j is finitely many, there exists a α d + 1 ( 1 , β ) such that α ( 1 , β ) , for any λ i ρ ( F K [ A , B , C | d ] ) .
According to such α , conbined with (38), for any λ j ρ ( F K [ A , B , C | d ] ) , | λ j | > 1 , we have λ α , j > 1 . Moreover, for any λ i ρ ( F K [ A , B , C | d ] ) , | λ i | 1 , it is easy to get λ α , i 1 . It follows from Lemma 5 that for any given W > 0 , (36) has a unique solution V = V , which implies that V is not a negative definite matrix by item (a) in Lemma 1.
To show the existence of positive eigenvalues of V, we have the following discussions.
(i)
When k e r ( V ) = 0 , V is a column full rank matrix. Since V is not a negative definite matrix, V must have a positive eigenvalue.
(ii)
When k e r ( V ) 0 , for any non-zero x 0 k e r ( V ) , pre-multiplying x 0 and post-multiplying x 0 on both sides of (36), we have
x 0 ( A α + B α K ) V ( A α + B α K ) x 0 + x 0 σ 2 ( ( A α ) d C α K ) V ( ( A α ) d C α K ) x 0 = x 0 W x 0 .
Together with the positive definiteness of W, since V is not a negative definite matrix, We have that V must not be a zero matrix. That is, V has a positive eigenvalue. □

6. Simulation

In this section, we will verify the validity of the developed theoretic results, including Theorem 1 and Theorem 3, through two illustrative examples.
Example 1.
Consider system [ A , B , C | d ] , with d = 1 , σ 2 = 1 , Q = I , R = I and A, B and C respectively meet the values as shown in Table 1.
First, we need to verify the mean square stabilization for each system in Table 1. Similar to Theorem 2 in [29], we have that system [ A , B , C | d ] is stabilizable if and only if there exist matrices H and Z > 0 satisfying
Z * * A Z + B H Z * σ 2 A d C H 0 Z < 0 ,
where * represents the corresponding transpose part. By using the LMI toolbox in MATLAB, we can obtain the corresponding values of H, Z > 0 for three different stochastic systems as shown in Table 1, which illustrates that each system is stabilizable.
It is remarkable that for DARE with Q > 0 and R > 0 , Theorem 1 defines a numerically iterative algorithm for obtaining the unique stabilizing solution, i.e., for any initial value P 0 0 , the matrix sequence satisfying P k + 1 = R [ A , B , C | d ] ( P k ) can converge to the unique solution. In Table 2, the value of P 0 is assumed to be positive semi-definite. Then, for any state dimension case, we can choose any positive semi-definite P 0 0 as an initial value. Perform iterative calculation P k + 1 = R [ A , B , C | d ] ( P k ) until P k 0 + 1 = P k 0 . In this case, we have P k 0 + 1 = P k 0 = P s . Besides, to illustrate the effectiveness of the developed algorithm, we provide the iteration N in Table 3. Based on the simulations, the convergence rate is fast.
By introducing the obtained lim k P k into Equation (11), we can verify that the limit is the stabilizing solution of DARE.
Example 2.
For system [ A , B , C | d ] , we take the following matrix coefficients into account
A = 1 5 1 3 , B = 1 0 2 0.7 , C = 1 2 1 1 , d = 1 .
In this case, one obtains that the control law follows
v k 1 = K x ^ k | k 2 = K ( A x k 1 + B u k 2 ) .
To begin with, we first check the given system is not stabilizable in the mean square sense by utilizing inequality (39). Then, to study the essential destabilization, we choose different feedback gain K i as shown in Table 4. It follows from Theorem 3, the following inequality holds for each K i , which can be solved by LMI in Matlab with
V ( A α + B α K ) V ( A α + B α K ) σ 2 ( A α ) d C α K V K C α ( A α ) d > 0 .
For the given matrix M, there exist α = 2 and a unique symmetric matrix, such that (36) holds. Utilizing the conventional calculation method, one obtains that V has at least one positive eigenvalue. Figure 1 shows the essential destabilization of the considered system.

7. Conclusions

We study the stochastic stabilization problems with the co-existence of input delay and multiplicative noise in the control variable. First, as the necessary condition of asymptotic mean square stabilization, we derive an iterative algorithm to solve the DARE. Then, the concepts of critical stabilization and essential destabilization are defined by operator spectrum theory. Utilizing DLE, the necessary and sufficient conditions are developed for dynamic models under consideration. However, how to utilize the spectrum theory and Lyapunov technology to study the stochastic stabilization for a more general time-delay stochastic system with multiplicative noises remains an open question, which defines a meaningful research direction. On the other hand, applying online Reinforcement Learning (RL) technology to solve DARE with partial system information is another promising research direction.

Author Contributions

Conceptualization, C.T.; methodology, C.T.; software, B.Z.; validation, C.T. and Z.C.; formal analysis, J.D. and M.X.; investigation, J.D. and M.X.; writing—original draft preparation, J.D. and M.X.; writing—review and editing, C.T.; supervision, C.T.; project administration, C.T.; funding acquisition, C.T. and Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant number 62173206; the Natural Science Foundation of Shandong Province grant number ZR2019QF005, ZR2021QF075, ZR2021ZD13; China Postdoctoral Science Foundation grant number 2021M691849.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank anonymous reviewers for constructive suggestions to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, L.; Xie, X. State feedback stabilization for stochastic feedforward nonlinear systems with time-varying delay. Automatica 2013, 49, 936–942. [Google Scholar] [CrossRef]
  2. Zhang, H.; Li, L.; Xu, J.; Fu, M.Y. Linear quadratic regulation and stabilization of discrete-time systems with delay and multiplicative noise. IEEE Trans. Autom. Control 2015, 60, 2599–2613. [Google Scholar] [CrossRef]
  3. Mesbah, A. Stochastic model predictive control: An overview and perspectives for future research. IEEE Control Syst. Mag. 2016, 36, 30–44. [Google Scholar]
  4. Jiang, X.; Zhao, D. Event-triggered fault detection for nonlinear discrete-time switched stochastic systems: A convex function method. Sci. China Inform. Sci. 2021, 64, 200204. [Google Scholar] [CrossRef]
  5. Zhang, T.; Deng, F.; Sun, Y.; Shi, P. Fault estimation and fault-tolerant control for linear discrete time-varying stochastic systems. Sci. China Inform. Sci. 2021, 64, 200201. [Google Scholar] [CrossRef]
  6. Wang, Z.; Ho, D.W.C.; Liu, Y.; Liu, X. Robust H∞ control for a class of nonlinear discrete time-delay stochastic systems with missing measurements. Automatica 2009, 45, 684–691. [Google Scholar] [CrossRef] [Green Version]
  7. Li, X.; Guay, M.; Huang, B.; Fisher, D.G. Delay-Dependent Robust H∞ Control of Uncertain Linear Systems with Input Delay. In Proceedings of the American Control Conference, San Diego, CA, USA, 2–4 June 1999. [Google Scholar]
  8. Zhao, C.-R.; Zhang, K.; Xie, X.-J. Output feedback stabilization of stochastic feedforward nonlinear systems with input and state delay. Int. J. Robust Nonlinear Control 2016, 26, 1422–1436. [Google Scholar] [CrossRef]
  9. Krstic, M. Delay Compensation for Nonlinear, Adaptive, and PDE Systems; Birkhäuser Boston: Boston, MA, USA, 2009. [Google Scholar]
  10. Cacace, F.; Germani, A.; Manes, C. Exponential stabilization of linear systems with time-varying delayed state feedback via partial spectrum assignment. Syst. Control Lett. 2014, 69, 47–52. [Google Scholar] [CrossRef]
  11. Hou, T.; Zhang, W.; Ma, H. Conditions for essential instability and essential destabilization of linear stochastic systems. In Proceedings of the World Congress on Intelligent Control and Automation, Jinan, China, 7–9 July 2010; pp. 1770–1775. [Google Scholar]
  12. Hou, T.; Ma, H. A small gain theorem for discrete-time stochastic periodic systems. In Proceedings of the IEEE International Conference on Control and Automation, Kathmandu, Nepal, 1–3 June 2016; pp. 282–287. [Google Scholar]
  13. Hou, T.; Zhang, W.; Chen, B.S. Regional stability and stabilizability of linear stochastic systems: Discrete-time case. In Proceedings of the IEEE International Conference on Control and Automation, Xiamen, China, 9–11 June 2010; pp. 1660–1665. [Google Scholar]
  14. Peaucelle, D.; Arzelier, D.; Bachelier, O.; Bernussou, J. A new robust D-stability condition for real convex polytopic uncertainty. Syst. Control Lett. 2000, 40, 21–30. [Google Scholar] [CrossRef]
  15. Dragan, V.; Morozan, T. Observability and detectability of a class of discrete-time stochastic linear systems. IMA J. Math. Control Inform. 2006, 23, 371–394. [Google Scholar] [CrossRef]
  16. Tan, C.; Zhang, H.; Wong, W.S. Delay-dependent algebraic Riccati equation to stabilization of networked control systems: Continuous-time case. IEEE Trans. Cybern. 2018, 48, 2783–2794. [Google Scholar] [CrossRef] [PubMed]
  17. Tan, C.; Zhang, H. Necessary and sufficient stabilizing conditions for networked control systems with simultaneous transmission delay and packet dropout. IEEE Trans. Autom. Control 2017, 62, 4011–4016. [Google Scholar] [CrossRef]
  18. Zhang, H.; Xie, L. Linear quadratic regulation for discrete-time systems with multiplicative noises and input delays. In Proceedings of the American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; pp. 1718–1723. [Google Scholar]
  19. Sipahi, R.; Niculescu, S.-I.; Abdallah, C.T.; Michiels, W.; Gu, K. Stability and Stabilization of Systems with Time Delay. IEEE Control Syst. Mag. 2011, 31, 38–65. [Google Scholar]
  20. Zhang, W.; Zhang, H.; Chen, B.S. Generalized Lyapunov Equation Approach to State-Dependent Stochastic Stabilization/Detectability Criterion. IEEE Trans. Autom. Control 2008, 53, 1630–1642. [Google Scholar] [CrossRef]
  21. Hou, T.; Zhang, W.; Ma, H. Essential instability and essential destabilisation of linear stochastic systems. IET Control Theory Appl. 2011, 5, 334–340. [Google Scholar] [CrossRef]
  22. Freiling, G.; Jank, G.; Abou-Kandil, H. Generalized Riccati difference and differential equations. Linear Algebra Its Appl. 1996, 241, 291–303. [Google Scholar] [CrossRef] [Green Version]
  23. Ran, A.C.M.; Vreugdenhil, R. Existence and comparison theorems for algebraic Riccati equations for continuous- and discrete-time systems. Linear Algebra Its Appl. 1988, 99, 63–83. [Google Scholar] [CrossRef] [Green Version]
  24. Huang, Y.; Zhang, W.; Zhang, H. Infinite horizon linear quadratic optimal control for discrete-time stochastic systems. Asian J. Control 2008, 10, 608–615. [Google Scholar] [CrossRef]
  25. Zhang, W.; Chen, B.S. On stabilizability and exact observability of stochastic systems with their applications. Automatica 2004, 40, 87–94. [Google Scholar] [CrossRef]
  26. Rami, M.A.; Zhou, X.Y. Linear Matrix Inequalities, Riccati Equations, and Indefinite Stochastic Linear Quadratic Controls. IEEE Trans. Autom. Control 2000, 45, 1131–1143. [Google Scholar] [CrossRef]
  27. Feng, Y.; Anderson, B.D.O. An iterative algorithm to solve state-perturbed stochastic algebraic Riccati equations in LQ zero-sum games. Syst. Control Lett. 2010, 59, 50–56. [Google Scholar] [CrossRef]
  28. Rami, M.A.; Chen, X.; Moore, J.B.; Zhou, X.Y. Solvability and asymptotic behavior of generalized Riccati equations arising in indefinite stochastic LQ controls. IEEE Trans. Autom. Control 2002, 46, 428–440. [Google Scholar] [CrossRef] [Green Version]
  29. Tan, C.; Yang, L.; Zhang, F.; Zhang, Z.; Wong, W.S. Stabilization of discrete time stochastic system with input delay and control dependent noise. Syst. Control Lett. 2019, 123, 62–68. [Google Scholar] [CrossRef]
  30. Tan, C.; Wong, W.S.; Zhang, H. On delay-dependent algebraic Riccati equation. IET Control Theory Appl. 2017, 11, 2506–2513. [Google Scholar] [CrossRef]
  31. Hou, T.; Zhang, W.; Chen, B.S. Study on general stability and stabilizability of linear discrete-time stochastic systems. Asian J. Control 2011, 13, 977–987. [Google Scholar] [CrossRef]
Figure 1. Simulation of lim k E | | x k | | 2 with K 1 .
Figure 1. Simulation of lim k E | | x k | | 2 with K 1 .
Processes 10 00989 g001
Table 1. Three sets of values of real system parameters in stochastic system [ A , B , C | d ] .
Table 1. Three sets of values of real system parameters in stochastic system [ A , B , C | d ] .
ABC
1 1 3 0 ; 0 1 5 1 1 5 ; 0 1 3 1 5 1 3 ; 0 1 5
2 1 1 5 ; 1 1 7 1 1 4 ; 2 1 4 1 3 1 5 ; 0 2
3 1 1 4 ; 1 1 7 2 1 3 ; 1 1 2 1 1 2 ; 1 3 1 4
Table 2. Three sets of simulation results of LMI toolbox.
Table 2. Three sets of simulation results of LMI toolbox.
HZ P 0
1 174.5693 52.8967 11.6459 276.5355 532.7641 0.3806 0.3806 531.4675 1 2 2 4
2 44.0803 6.0853 150.2109 29.4270 109.7799 1.6274 1.6274 115.1478 4 6 6 9
3 69.6452 1.8577 113.4397 18.3662 117.1685 2.5141 2.5141 117.3061 1 0 0 1
Table 3. Results of three sets of simulation series.
Table 3. Results of three sets of simulation series.
lim x P k N
1 1.0304 0.9834 0.9834 1.0120 > 011
2 1.6470 0.8792 0.8792 1.0231 > 027
3 2.6428 0.8405 0.8405 1.0700 > 019
Table 4. Results of three kinds of feedback design for the system.
Table 4. Results of three kinds of feedback design for the system.
K 10 8   W 10 6   V 10 5   λ V
1 1 0 0 1 0.3936 0.0582 0.0582 1.8208 6.5492 6.8033 6.8033 6.8989 0.8200 135.00
2 5 0 0 5 0.6375 0.2810 0.2810 1.8007 0.5420 0.48857 0.48857 0.476510 0.0398 9.8566
3 0.5 0 0 0.5 0.6467 0.5164 0.5164 2.5623 6.0390 22.4000 22.4000 31.8010 692.0000 447.5900
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, C.; Di, J.; Xiang, M.; Chen, Z.; Zhu, B. Stochastic Stabilization for Discrete-Time System with Input Delay and Multiplicative Noise in Control Variable. Processes 2022, 10, 989. https://0-doi-org.brum.beds.ac.uk/10.3390/pr10050989

AMA Style

Tan C, Di J, Xiang M, Chen Z, Zhu B. Stochastic Stabilization for Discrete-Time System with Input Delay and Multiplicative Noise in Control Variable. Processes. 2022; 10(5):989. https://0-doi-org.brum.beds.ac.uk/10.3390/pr10050989

Chicago/Turabian Style

Tan, Cheng, Jianying Di, Mingyue Xiang, Ziran Chen, and Binlian Zhu. 2022. "Stochastic Stabilization for Discrete-Time System with Input Delay and Multiplicative Noise in Control Variable" Processes 10, no. 5: 989. https://0-doi-org.brum.beds.ac.uk/10.3390/pr10050989

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