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:
where
stands for the Wiener process,
stands for a smooth non-linear function,
denotes a real constant,
and
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,
where
denotes the new control signal while the control signal
is the integral of
.
Then control objective can be described as designing such that is bounded in probability sense, then the variance and entropy of track the given desired values.
To achieve the system design, the preliminaries [
15] are recalled as follows:
Definition 1. For any given associated with the stochastic differential Equation (1), the differential operator is defined as follows: Lemma 1. Consider the stochastic non-linear system model (1) and assume that is in the arguments and is bounded uniformly in t. If there exist functions , , , constants , , and a nonnegative function , such that 3. System Stabilisation
Based on the backstepping procedure [
16], the virtual input signal can be designed as
where
denotes the designed parameter for the controller.
In order to stabilise Equation (1), the virtual tracking error is formulated as
Substituting the error signal
into Equation (32), we have
Based upon Itô’s lemma, the following result is produced.
One Lyapunov function candidate is selected to stabilise the investigated system model (32), which is given as follows.
which results in
Using Definition 1, Lemma 1 and Young’s inequality, the following results can be evaluated as,
and
The control signal can be further developed as
where
C stands for a designed real function.
Substituting the control signal into
, Equation (10) can be written as
The system state
using
is bounded in probability sense based on Lemma 1. Moreover,
C can be further selected as follows to eliminate the
nonlinear term.
where
denotes a free designed parameter.
Substituting Equations (15)–(17) into the controller design (13), the complete control scheme can then be formulated as follows:
which leads to
Thus, the closed-loop system with the designed parameter
is bounded in probability sense.
4. Variance and Entropy Assignment
Substituting the control signal into
shows that the error signal
is also bounded in the probability sense once the free parameter meets
. Thus, the closed-loop system with the given control scheme can be further considered as the Ornstein–Uhlenbeck process.
while the Fokker–Planck–Kolmogorov equation can be formulated as follows:
where
denotes the PDF and
x denotes the random variable of
.
Note that the associated FPK equation is a linear partial differential equation; the analytical solution can be formulated analytically as follows:
where
stands for the initial value of
at
.
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
to re-shape as a Gaussian distribution. In particular, the formula of mean value and variance can be obtained as follows:
where
and
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
is considered as the convergence rate.
Denoting a desired variance function as
,
can then be further developed as a function of
t such that the following equation holds.
Note that this equation can be rewritten as the following form,
Using the Lambert W function [
17], the equation can be solved and the solution is given as follows:
where
is the Lambert W function. Note that
can be adopted when
.
To ensure the
function can be used above, the stationary solution can be practically implemented if
is a real positive constant.
Note that the solution meets the condition of
where
.
We can achieve the entropy assignment for the random variable
as the variance value is correlated to Shannon’s entropy subjected to the linear stochastic process. Based on the definition, we have
where
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.
where
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
is a real constant.
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:
To show the assignment procedure, the reference variance value of
is given as
while it leads to the parameter selection
. In particular,
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:
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
. 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
is shown in
Figure 1, where the system variable stabilisation is achieved, while the control input signal with integrator
and the signal without integrator
are given in
Figure 2 and
Figure 3, respectively.
Figure 4 shows the randomness attenuation where the virtual tracking error
reaches zero. The variance values of
and
are indicted in
Figure 5 and
Figure 6, in which the variance value of
converges to the assigned reference value and the transient error comes from the virtual tracking error stabilisation progress. The mean values of
and
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
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
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.
where
stands for a known smooth non-linear function
,
denotes n-dimensional Wiener process, Σ denotes a given as a real positive square matrix with n dimensions.
and
stand for the system state vector and system output, respectively.
denotes the control input.
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.