Next Article in Journal
Variational Embedding Multiscale Sample Entropy: A Tool for Complexity Analysis of Multichannel Systems
Next Article in Special Issue
Multi-Sensor Vibration Signal Based Three-Stage Fault Prediction for Rotating Mechanical Equipment
Previous Article in Journal
Entropy Effects in Intermolecular Associations of Crown-Ethers and Cyclodextrins with Amino Acids in Aqueous and in Non-Aqueous Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variance and Entropy Assignment for Continuous-Time Stochastic Nonlinear Systems

1
Engineering Research Center of the Ministry of Education (Power Grid Security Monitoring and Control Technology), Changsha University of Science and Technology, Changsha 410114, China
2
School of Engineering and The Built Environment, Edinburgh Napier University, Edinburgh EH11 4BN, UK
3
School of Automation, Central South University, Changsha 410083, China
4
Department of Computer Science, University of Bradford, Bradford BD7 1DP, UK
*
Author to whom correspondence should be addressed.
Submission received: 4 November 2021 / Revised: 9 December 2021 / Accepted: 17 December 2021 / Published: 24 December 2021

Abstract

:
This paper investigates the randomness assignment problem for a class of continuous-time stochastic nonlinear systems, where variance and entropy are employed to describe the investigated systems. In particular, the system model is formulated by a stochastic differential equation. Due to the nonlinearities of the systems, the probability density functions of the system state and system output cannot be characterised as Gaussian even if the system is subjected to Brownian motion. To deal with the non-Gaussian randomness, we present a novel backstepping-based design approach to convert the stochastic nonlinear system to a linear stochastic process, thus the variance and entropy of the system variables can be formulated analytically by the solving Fokker–Planck–Kolmogorov equation. In this way, the design parameter of the backstepping procedure can be then obtained to achieve the variance and entropy assignment. In addition, the stability of the proposed design scheme can be guaranteed and the multi-variate case is also discussed. In order to validate the design approach, the simulation results are provided to show the effectiveness of the proposed algorithm.

1. Introduction

Stochastic systems, which are usually subjected to high levels of uncertainties and randomness, have become one of the major research fields due to their presence in the real-life systems. The randomness existing in the stochastic systems will result in bad control behaviour and lead to the instability of the controlled systems. Therefore, to enhance the performance of such systems, the influence of randomness has to be well controlled. To characterise the randomness of systems, variance has always been adopted as a statistical index when the randomness can be characterised as a Gaussian process, thus showing that the variance control method is regarded as an important design for stochastic systems analysis and implementation [1,2,3].
However, for the systems that cannot meet the Gaussian assumptions due to the nonlinearities, these results from the variance control method cannot be applied directly. It has been shown that the distributions of system variables can be twisted to non-Gaussian even if the noises in nonlinear systems are Gaussian noises. Based on the probability theory, variance is only the second moment which cannot reflect the full property of the non-Gaussian randomness [4]. To address that, the probability density function (PDF) control [5] becomes the solution for the stochastic nonlinear systems as PDF contains the full stochastic properties of random variables. The existing PDF control methods mostly focus on discrete-time systems. For the continuous-time stochastic nonlinear systems, it is difficult to find a solution analytically for the Fokker–Planck–Kolmogorov (FPK) equation [6]. Even though, entropy control [7] has been developed to overcome the non-Gaussian property and high-order moment problem of variance control, the analytical formulation of entropy for continuous-time dynamic system remains challenge since the FPK equation is hard to solve. Therefore, randomness control for continuous-time stochastic nonlinear systems is still regarded as a challenging topic.
We notice that the key factor for the randomness control is the system nonlinearity as the linear stochastic system can be handled by variance control. For linear stochastic systems, the entropy value is equivalent to variance. Motivated by Brownian motion [8,9], the randomness control for the stochastic nonlinear system can be achieved if the nonlinear process is able to be converted to linear process via the control design. Furthermore, the variance and entropy assignment can be achieved analytically to reflect the randomness. Following this idea, we consider a class of continuous-time stochastic nonlinear systems in this paper, where the backstepping design is adopted to stabilise the system variables. At the same time, the investigated nonlinear systems can be converted to linear format with designed parameter. The converted linear system structure can be analysed as an Ornstein–Uhlenbeck process [6], where the associated Fokker–Planck–Kolmogorov equation is solvable. In particular, the variance and entropy can be formulated analytically. Based on the formula, the designed parameters can be further obtained to achieve the randomness assignment.
Different from the stochastic distribution control [10], the randomness is investigated in a simplified approach and the analytical solution is obtained. Furthermore, the block backstepping design [11] and covariance assignment method [12] can be applied for multi-variate systems. In practice, the randomness assignment can also be used for filtering design [13], system identification and applications [14].
The rest of the paper is organised as follows. In Section 2, the formulation is given in terms of problem description and preliminaries related to the main results. Section 3 and Section 4 are the main contents of this paper, where the backstepping-based control law is proposed to stabilise the investigated stochastic system and the design parameter is analysed for variance assignment using Fokker–Planck–Kolmogorov equation. To validate the proposed control design, a numerical example is given in Section 5, where the results show that the variance can be assigned using the proposed control algorithm. The multi-variate system extension is discussed in Section 6 and the conclusions are given in Section 7 as the last part of this paper.

2. Formulation

In this paper, we first consider a univariate continuous-time stochastic affine nonlinear system which can be modelled by the following Itô process:
d x t = f x t + v t d t + σ d W t
where W t stands for the Wiener process, f · stands for a smooth non-linear function, σ > 0 denotes a real constant, x t and v t denote the system state and control input, respectively.
We can further re-write the investigated system model (1) by introducing an integrator while the following model can be obtained,
d x t = f x t + v t d t + σ d W t d v t = u t d t
where u t denotes the new control signal while the control signal v t is the integral of u t .
Then control objective can be described as designing u t such that x t is bounded in probability sense, then the variance and entropy of x t track the given desired values.
To achieve the system design, the preliminaries [15] are recalled as follows:
Definition 1. 
For any given V x t , t C 1 , 2 R n × R + ; R + associated with the stochastic differential Equation (1), the differential operator L is defined as follows:
L V = V t + V x t f x t + V x t v t + 1 2 T r 2 V x t 2 σ 2
Lemma 1. 
Consider the stochastic non-linear system model (1) and assume that f x t is C 1 in the arguments and f 0 is bounded uniformly in t. If there exist functions V x t , t C 1 , 2 R n × R + ; R + , μ 1 · , μ 2 · K , constants c 1 > 0 , c 2 0 , and a nonnegative function V ¯ x t , t , such that
μ 1 | x t | V x t , t μ 2 | x t | L V c 1 V ¯ x t , t + c 2

3. System Stabilisation

Based on the backstepping procedure [16], the virtual input signal can be designed as
ϕ x t = f x t θ x t
where θ > 0 denotes the designed parameter for the controller.
In order to stabilise Equation (1), the virtual tracking error is formulated as
z t = v t ϕ x t = v t + f x t + θ x t
Substituting the error signal z t into Equation (32), we have
d x t = θ x t + z t d t + σ d W t
Based upon Itô’s lemma, the following result is produced.
d z t = d v t d ϕ x t = u t θ x t + z t ϕ x t x σ 2 2 2 ϕ x t x 2 d t σ ϕ x t x d W t
One Lyapunov function candidate is selected to stabilise the investigated system model (32), which is given as follows.
V t = V x + V z = 1 2 x t 2 + 1 4 z t 4
which results in
L V t = L V x + L V z
Using Definition 1, Lemma 1 and Young’s inequality, the following results can be evaluated as,
L V x = x t θ x t + z t + σ 2 2 = θ x t 2 + x t z t + σ 2 2 θ x t 2 + 1 2 x t 2 + 1 2 z t 2 + σ 2 2 = θ + 1 2 x t 2 + z t 2 + σ 2 2 θ + 1 2 x t 2 + 1 2 z t 4 + σ 2 + 1 2
and
L V z = z t 3 u t θ x t + z t ϕ x t x σ 2 2 2 ϕ x t x 2 + 3 σ 2 2 ϕ x t x 2 z t 2 z t 3 u t + θ x t ϕ x t x z t ϕ x t x σ 2 2 2 ϕ x t x 2 + 1 2 + 9 σ 4 8 ϕ x t x 4 z t 4 = z t 3 u t + θ x t ϕ x t x σ 2 2 2 ϕ x t x 2 + 9 σ 4 8 ϕ x t x 4 ϕ x t x z t 4 + 1 2
The control signal can be further developed as
u t = θ x t ϕ x t x + σ 2 2 2 ϕ x t x 2 C z t
where C stands for a designed real function.
Substituting the control signal into L V z , Equation (10) can be written as
L V t θ + 1 2 x t 2 + 1 2 z t 4 + σ 2 + 2 2 C z t 4 + 9 σ 4 8 ϕ x t x 4 ϕ x t x z t 4 = θ + 1 2 x t 2 + σ 2 + 2 2 + 1 2 C + 9 σ 4 8 ϕ x t x 4 ϕ x t x z t 4
The system state x t using u t is bounded in probability sense based on Lemma 1. Moreover, C can be further selected as follows to eliminate the z 4 nonlinear term.
C = 1 2 9 σ 4 8 ϕ x t x 4 + ϕ x t x κ
where κ 0 denotes a free designed parameter.
As a result, we have
ϕ x t x = f x t x θ
and
2 ϕ x t x 2 = 2 f x t x 2
Substituting Equations (15)–(17) into the controller design (13), the complete control scheme can then be formulated as follows:
u t = θ x t f x t x + θ σ 2 2 2 f x t x 2 1 2 9 σ 4 8 ϕ x t x 4 + ϕ x t x κ × v t + f x t + θ x t
which leads to
L V t θ + 1 2 x t 2 κ z t 4 + σ 2 + 2 2
Thus, the closed-loop system with the designed parameter θ 1 2 is bounded in probability sense.

4. Variance and Entropy Assignment

Substituting the control signal into L V z shows that the error signal z t is also bounded in the probability sense once the free parameter meets κ > 0 . Thus, the closed-loop system with the given control scheme can be further considered as the Ornstein–Uhlenbeck process.
d x t = θ x t d t + σ d W t
while the Fokker–Planck–Kolmogorov equation can be formulated as follows:
p ( x , t ) t = θ x x p ( x , t ) + σ 2 2 2 p ( x , t ) x 2
where p ( x , t ) denotes the PDF and x denotes the random variable of x t .
Note that the associated FPK equation is a linear partial differential equation; the analytical solution can be formulated analytically as follows:
p ( x , t ) = θ π σ 2 1 e 2 t θ exp θ σ 2 x x 0 e t θ 2 1 e 2 t θ )
where x 0 stands for the initial value of x t at t 0 .
Since the solution of this FPK equation is a Gaussian distribution, the mean value and variance can be obtained analytically. It implies that the presented control scheme governs the non-Gaussian PDF of x t to re-shape as a Gaussian distribution. In particular, the formula of mean value and variance can be obtained as follows:
E x t = x 0 e t θ
V a r x t = σ 2 2 θ 1 e 2 t θ
where E · and V a r · stand for the operators of the mean value and the variance value. It has been shown that both mean value and variable can be adjusted via the designed parameter turning in the presented control scheme θ . In addition, the mean value reaches zero with θ > 0 is considered as the convergence rate.
Denoting a desired variance function as r t , θ can then be further developed as a function of t such that the following equation holds.
σ 2 2 θ 1 e 2 t θ = r t
Note that this equation can be rewritten as the following form,
e 2 t θ = 2 r t σ 2 θ + 1
Using the Lambert W function [17], the equation can be solved and the solution is given as follows:
θ = W 0 t σ 2 r t e t σ 2 r t 2 t + σ 2 2 r t
where W 0 · is the Lambert W function. Note that W 0 · can be adopted when 2 r c θ σ 2 .
To ensure the W 0 · function can be used above, the stationary solution can be practically implemented if r t = r c is a real positive constant.
θ s = lim t θ = σ 2 2 r c
Note that the solution meets the condition of W 0 · where 2 r c θ s σ 2 .
We can achieve the entropy assignment for the random variable x t as the variance value is correlated to Shannon’s entropy subjected to the linear stochastic process. Based on the definition, we have
H x t = 1 2 log 2 π V a r x t + 1
where H · stands for the Shannon’s entropy.
Based on the equation above, we can establish the link between the desired variance function and the pre-specified entropy function.
r t = 1 2 π e 2 H r t 1
where H r t stands for the pre-specified desired entropy function.
Substituting Equation (30) into Equation (28), the parameter θ can be confirmed. Similar to the variance assignment, the stationary solution for entropy assignment is achieved if H r t = H c is a real constant.
θ s = π σ 2 e 1 2 H c
As Shannon’s entropy is regarded as a special case of Rényi’s entropy, the presented control scheme can be extended to various entropies. As a result, the entropy optimisation can be further achieved for non-Gaussian filtering designs [13] and performance enhancement [18].

5. Simulation

The following numerical example is designed for validating the presented scheme. Matching the investigated system model, the parameters are confirmed as follows:
d x t = x t 3 + v t d t + 0.2 d W t d v t = u t d t
To show the assignment procedure, the reference variance value of x t is given as 0.02 while it leads to the parameter selection θ = 1 . In particular, θ = 1 > 1 2 implies that the designed system should be bounded in the probability sense. To start up the backstepping design, the virtual tracking error can be described by the following dynamic model:
d z t = u t + z t x t 3 x t 3 + 1 + 0.04 3 x t d t + 0.2 3 x t 2 1 d W t
Thus, the procedure can be implemented for achieving the proposed control objective.
To show the simulation results, we first discretise the given model using with the sampling time 0.1 . Then the computational results are demonstrated in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. In particular, the trajectory of the system state x t is shown in Figure 1, where the system variable stabilisation is achieved, while the control input signal with integrator v t and the signal without integrator u t are given in Figure 2 and Figure 3, respectively. Figure 4 shows the randomness attenuation where the virtual tracking error z t reaches zero. The variance values of x t and z t are indicted in Figure 5 and Figure 6, in which the variance value of x t converges to the assigned reference value and the transient error comes from the virtual tracking error stabilisation progress. The mean values of x t and z t are also given in Figure 7 and Figure 8, respectively. Note that both mean values are approaching zero. However, Lemma 2 implies that the error of the mean value of z t still exists due to the Lyapunov theorem-based analysis, where an arbitrary small non-zero error exists. As a result, the mean value and variance value of x t may be affected by the non-zero virtual tracking error, although it is arbitrarily small.

6. Discussion

In the system design above, only the univariate system was investigated to indicate the framework via backstepping. To extend this result, the main challenge comes from the multivariate backstepping design. As a solution, the block backstepping design is an implementable solution similar to the design procedure in [11]. In particular, we can consider the following extended system model with multi-dimensional system variables.
d x ¯ t = f ¯ x ¯ t + v ¯ t d t + Σ d W ¯ t y t = C x t
where f ¯ · stands for a known smooth non-linear function f ¯ : R n × 1 R n , W ¯ t denotes n-dimensional Wiener process, Σ denotes a given as a real positive square matrix with n dimensions. x ¯ t R n and y t R 1 stand for the system state vector and system output, respectively. v ¯ t R n denotes the control input. C R n denotes a vector-valued coefficient.
Following the presented design approach, the candidate of Lyapunov functions can be re-used where the vector-value variables will be used. Since Lemma 1 holds for multivariate system, the developed result in this paper can be extended directly following the block backstepping design. Notice that the linear Ornstein–Uhlenbeck process will be in the multi-dimensional form which leads to the difficulty of solving the FPK equation, as the joint probability density function has to be involved into the multivariate case. To avoid this problem, the design parameter θ should be selected as the positive diagonal matrix. Then a set of FPK equations can be obtained where the vector state can be decomposed as single variables. Therefore, the presented parameter selection scheme can also be re-used for multivariate systems. Alternatively, the converted linear multivariate systems can be further adjusted by covariance assignment which means an additional control signal will be introduced into the system design.

7. Conclusions and Perspectives

In this paper, a new randomness assignment framework is proposed for continuous-time stochastic nonlinear systems which are described by stochastic differential equations. The core idea is to convert the nonlinear system via control design, then the randomness can be fully characterised by variance. In particular, the backstepping procedure is first used to stabilise the variables of the investigated system with designed parameters. In this way, the converted linear system can be further described by the Ornstein–Uhlenbeck process. Thus, the associated Fokker–Planck–Kolmogorov equation is then analytically solvable which leads to the formula of variance and entropy value. Then, the assignment can be achieved by selecting designed parameter of backstepping. A numerical example is given as a validation of the presented method and a discussion is further given for the multivariate systems. In the future, the covariance control theory can be further merged into the presented framework to enhance the flexibility of the assignment for the multivariate cases.

Author Contributions

Conceptualisation, Q.Z.; methodology, Q.Z.; validation, Q.Z., Y.Z. (Yiqun Zou) and X.T.; formal analysis, Q.Z.; writing—original draft preparation, Q.Z.; writing—review and editing, Y.Z. (Yuyang Zhou); project administration, X.T. and Y.Z. (Yiqun Zou); funding acquisition, X.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (NSFC) under grant number 51807010 and Natural Science Foundation of Hunan Province under grant number 2020JJ5585 and 2020JJ5777.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors would like to thank the reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, P.R.; Varaiya, P. Stochastic Systems: Estimation, Identification, and Adaptive Control; SIAM: Philadelphia, PA, USA, 2015. [Google Scholar]
  2. Hotz, A.; Skelton, R.E. Covariance control theory. Int. J. Control 1987, 46, 13–32. [Google Scholar] [CrossRef]
  3. Åström, K.J. Introduction to Stochastic Control Theory; Courier Corporation: North Chelmsford, MA, USA, 2012. [Google Scholar]
  4. Ren, M.; Zhang, Q.; Zhang, J. An introductory survey of probability density function control. Syst. Sci. Control Eng. 2019, 7, 158–170. [Google Scholar] [CrossRef]
  5. Wang, H. Bounded Dynamic Stochastic Systems: Modelling and Control; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  6. Ross, S.M.; Kelly, J.J.; Sullivan, R.J.; Perry, W.J.; Mercer, D.; Davis, R.M.; Washburn, T.D.; Sager, E.V.; Boyce, J.B.; Bristow, V.L. Stochastic Processes; Wiley: New York, NY, USA, 1996; Volume 2. [Google Scholar]
  7. Wang, H. Minimum entropy control of non-Gaussian dynamic stochastic systems. IEEE Trans. Autom. Control 2002, 47, 398–403. [Google Scholar] [CrossRef]
  8. Ayache, A.; Vehel, J.L. The generalized multifractional Brownian motion. Stat. Inference Stoch. Process. 2000, 3, 7–18. [Google Scholar] [CrossRef]
  9. El Omari, M. Mixtures of higher-order fractional Brownian motions. Commun. Stat.-Theory Methods 2021, 1–16. [Google Scholar] [CrossRef]
  10. Zhang, Q.; Wang, H. A Novel Data-based Stochastic Distribution Control for Non-Gaussian Stochastic Systems. IEEE Trans. Autom. Control 2021. [Google Scholar] [CrossRef]
  11. Zhang, Q.C.; Hu, L.; Gow, J. Output feedback stabilization for MIMO semi-linear stochastic systems with transient optimisation. Int. J. Autom. Comput. 2020, 17, 83–95. [Google Scholar] [CrossRef]
  12. Zhang, Q.; Wang, Z.; Wang, H. Parametric covariance assignment using a reduced-order closed-form covariance model. Syst. Sci. Control Eng. 2016, 4, 78–86. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, Q. Performance enhanced Kalman filter design for non-Gaussian stochastic systems with data-based minimum entropy optimisation. AIMS Electron. Electr. Eng. 2019, 3, 382–396. [Google Scholar] [CrossRef]
  14. Zhang, Q.; Sepulveda, F. Entropy-based axon-to-axon mutual interaction characterization via iterative learning identification. In Proceedings of the EMBEC & NBC 2017, Tampere, Finland, 11–15 June 2017; pp. 691–694. [Google Scholar]
  15. Liu, S.J.; Zhang, J.F.; Jiang, Z.P. Decentralized adaptive output-feedback stabilization for large-scale stochastic nonlinear systems. Automatica 2007, 43, 238–251. [Google Scholar] [CrossRef]
  16. Deng, H.; Krstic, M. Output-feedback stochastic nonlinear stabilization. IEEE Trans. Autom. Control 1999, 44, 328–333. [Google Scholar] [CrossRef] [Green Version]
  17. Mezo, I. Lambert W Function: Its Generalizations and Applications; Chapman and Hall/CRC: London, UK, 2022. [Google Scholar]
  18. Tang, X.; Zhang, Q.; Hu, L. An EKF-based performance enhancement scheme for stochastic nonlinear systems by dynamic set-point adjustment. IEEE Access 2020, 8, 62261–62272. [Google Scholar] [CrossRef]
Figure 1. The trajectory of the system state x t .
Figure 1. The trajectory of the system state x t .
Entropy 24 00025 g001
Figure 2. The control input with integrator v t .
Figure 2. The control input with integrator v t .
Entropy 24 00025 g002
Figure 3. The control input without integrator u t .
Figure 3. The control input without integrator u t .
Entropy 24 00025 g003
Figure 4. The trajectory of the virtual error z t .
Figure 4. The trajectory of the virtual error z t .
Entropy 24 00025 g004
Figure 5. The variance value curve of the system state x t .
Figure 5. The variance value curve of the system state x t .
Entropy 24 00025 g005
Figure 6. The variance value curve of the virtual error z t .
Figure 6. The variance value curve of the virtual error z t .
Entropy 24 00025 g006
Figure 7. The mean value curve of the system state x t .
Figure 7. The mean value curve of the system state x t .
Entropy 24 00025 g007
Figure 8. The mean value curve of the virtual error z t .
Figure 8. The mean value curve of the virtual error z t .
Entropy 24 00025 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, X.; Zhou, Y.; Zou, Y.; Zhang, Q. Variance and Entropy Assignment for Continuous-Time Stochastic Nonlinear Systems. Entropy 2022, 24, 25. https://doi.org/10.3390/e24010025

AMA Style

Tang X, Zhou Y, Zou Y, Zhang Q. Variance and Entropy Assignment for Continuous-Time Stochastic Nonlinear Systems. Entropy. 2022; 24(1):25. https://doi.org/10.3390/e24010025

Chicago/Turabian Style

Tang, Xiafei, Yuyang Zhou, Yiqun Zou, and Qichun Zhang. 2022. "Variance and Entropy Assignment for Continuous-Time Stochastic Nonlinear Systems" Entropy 24, no. 1: 25. https://doi.org/10.3390/e24010025

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