Next Article in Journal
An IND-CPA Analysis of a Cryptosystem Based on Bivariate Polynomial Reconstruction Problem
Next Article in Special Issue
Disturbance Attenuation Trajectory Tracking Control of Unmanned Surface Vessel Subject to Measurement Biases
Previous Article in Journal
Analysis of =P[Y<X<Z] Using Ranked Set Sampling for a Generalized Inverse Exponential Model
Previous Article in Special Issue
Optimization of Leakage Risk and Maintenance Cost for a Subsea Production System Based on Uncertain Fault Tree
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on PDF Shape Control for Nonlinear Stochastic System Using an Approximate Solution of FPK Equation

1
School of Automation, Xi’an University of Posts and Telecommunications, Xi’an 710121, China
2
School of Automation and Information Engineering, Xi’an University of Technology, Xi’an 710048, China
3
Shaanxi Key Laboratory of Surface Engineering and Remanufacturing, Xi’an University, Xi’an 710065, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(3), 303; https://doi.org/10.3390/axioms12030303
Submission received: 14 February 2023 / Revised: 13 March 2023 / Accepted: 13 March 2023 / Published: 17 March 2023
(This article belongs to the Special Issue Advances in Uncertain Optimization and Applications)

Abstract

:
In this paper, we developed a probability density function (PDF) shape control method for non-linear stochastic systems using a hybrid logistic function (HLF) as an approximate PDF of the state variable. First, the functional relationship between the hybrid logistic probability density function and the controller was established based on the Fokker–Planck–Kolmogorov (FPK) equation. Then, the optimal PDF shape controller derivation was completed using the optimization method and the inner product definition of Hilbert space. This approach is suitable for any non-linear stochastic system. To evaluate the effectiveness and performance of the proposed method, we conducted a comparison experiment with the multi-Gaussian closure (MGC) method and the exponential polynomial (EP) method. The experimental results show that, for different types of targeted PDFs (symmetric unimodal, asymmetric unimodal, bimodal, and trimodal shapes), the PDF shape controller obtained using the HLF approach can make the PDF shape of the state variable track the targeted PDF effectively. In particular, when the targeted PDF has an asymmetric or complex trimodal shape, the proposed technique has comparatively better control effects. Compared with the EP method, our method requires a much smaller number of parameters, greatly reducing the computational complexity while achieving the same control effects. This study provides another approach for controlling the PDF shape of state variables in non-linear stochastic systems, which has important research significance.

1. Introduction

In nature, stochastic phenomena are widespread, and many practical systems are not immune to their effects. Therefore, stochastic systems have been widely studied in many disciplines and engineering applications. For stochastic systems, stochastic control is essential. Stochastic control has achieved rapid development in theory in the past few decades, and its research results have been widely used in industry, economy, national defense, aerospace, and other fields, solving many practical problems. The most representative stochastic control methods are minimum variance control and linear quadratic Gaussian control, which effectively solve the mean control problem of linear systems [1,2]. However, these control methods are unable to suppress fluctuations generated from random factors, such as waves, track unevenness, voltage fluctuations, and vibrations of the vehicle engine.
To suppress these fluctuations, variance control is applied. However, the information contained in the mean and variance is limited and cannot express the complete dynamic characteristics of the system, which makes the desired goal difficult to achieve. Instead, the shape of the probability density function (PDF) is used to reflect the complete statistical characteristics of the random variable [3,4,5]. Therefore, in random systems, studies that utilize the shape of the probability density function as the control object are more comprehensive and generate more accurate results than the mean and variance control.
There are numerous PDF shape control problems, such as food processing, the emergency vehicle damping system, and the flight control system. In 1998–2009, Wang, Guo and Luan et al. obtained many relevant research achievements [6,7,8,9,10,11,12,13,14,15,16]. They mainly utilized B-sample approximation, minimum entropy control, and the linear matrix inequality (LMI) method for studying the PDF control problem. In 2012, Qian et al. addressed the complete statistical properties of discrete LQG and studied the moment control [17]. The object in these studies is the linear stochastic system.
Stochastic systems exhibit strong or weak non-linear behavior caused by various non-linear factors, so most stochastic systems are non-linear [18,19,20,21,22,23,24,25,26,27]. Therefore, the research on PDF shape control has important value and significance in non-linear systems.
In recent years, a number of studies on PDF shape control for non-linear stochastic systems have emerged [28,29,30,31,32,33,34,35,36,37,38,39,40]. Forbes et al. suggested a shape control scheme of a steady-state PDF in discrete stochastic processes to approach the desired PDF by use of Gram–Charlier expansion [28,29,30,31]. For the ARMAX (AutoRegressive Moving Average with eXogenous inputs) system with nonlinear discrete bounded parameter uncertainty, Guo et al. achieved the PDF shape control of the system output [32]. Zhu et al. put forward a method based on an equivalent non-linear system to obtain the approximate steady-state solution and designed a feedback control of non-linear stochastic systems so that the output’s PDF can track the target [33]. Qian and Zhao used the abnormal integral method to transform the non-linear system and determined the tracking control of the PDF by seeking the approximate solution [34]. Based on the relationship between the characteristic function and the probability density function, the PDF of the state response was reconstructed. In 2016, Yang et al. expressed the probability density function with an exponential function and used the FPK equation to find the derivatives of each order of the PDF [35]. In another study, the shape control problem with a PDF was transformed into a tracking optimization problem, and the optimal control gain was obtained using a particle swarm optimization algorithm [36]. In the reference [37], the superimposed PDF with a number of Gaussian PDFs was constructed according to the target PDF. This allows for the tracking of the target PDF, and the controller function is able to pass through the FPK equation. In 2014, Wang et al. proposed a PDF shape control algorithm based on the PSO algorithm, which effectively solved the PDF shape control problem for non-linear stochastic systems [38]. In 2015, a PDF shape control method was proposed using the linear least squares method for a class of non-linear stochastic systems with additive white Gaussian noise and polynomial form, overcoming the challenges of solving the FPK equation [39]. In 2016, a piecewise linear PDF shape controller was developed for non-linear stochastic systems with polynomial form [40]. This controller has fewer parameters and has recognizable advantages over polynomial controllers.
The above methods mainly solve the problem of PDF shape control for polynomial non-linear stochastic systems.
In this paper, we developed a PDF shape control method for any non-linear stochastic system based on hybrid logistic functions to approximate the target probability density function. First, the FPK equation was used to establish the functional relationship between the hybrid logistic PDF and the controller, and the optimal controller was obtained using an optimization method and the inner product of Hilbert space. Taking the probability density functions with unimodal and multimodal shapes as the target PDF, we compared the performance of our proposed approach against the MGC method and the EP method through simulation experiments. We then calculated and compared the mean, variance, and control error of the three methods.
The main innovations and improvements of the proposed approach are as follows:
(1)
The obtained PDF shape controller can be suitable for different target PDFs.
(2)
This method can be applied to any non-linear stochastic system and is no longer limited to polynomial non-linear stochastic systems.
(3)
If the target PDF is asymmetric, this method can achieve a better control effect compared to other techniques.
(4)
When obtaining the same control effect as the EP method, our approach has much fewer parameters to be solved and has a lower complexity.
The rest of the paper is arranged as follows. Section 2 provides the problem statement, and the hybrid logistic function method is introduced in Section 3. In Section 4, we provide details on the derivation of the PDF shape controller based on the hybrid logistic probability density function. Section 5 discusses sample simulations and calculation results. An analysis of results is conducted in Section 6, and the research conclusions are presented in Section 7.

2. Problem Statement

Consider a single-degree-of-freedom non-linear system in random environments:
d x d t = y x +   w t x t 0 =   x 0
where x t R is the state variable describing the system; w t R is Gaussian white noise, E w t = 0 , and E w t w t + τ = 2 π S 0 δ τ , where E · is the expectation, S 0 is the power spectral density of Gaussian white noise, and δ · denotes the (Dirac) delta function; y x R is a non-linear function of x; and x 0 is the initial state determined by the system at the initial time t 0 .
In this study, we aimed to determine a state feedback function u x so that the probability density function p x of the state variable can track the shape of the expected probability density function p d x under the regulation of u x when the system tends to be stable. Here, p d x is determined in advance and must conform to the distribution of a real random variable.
Therefore, the system (1) can be rewritten after applying the controller u x as:
d x d t = y x   + u x + w t x t 0 = x 0

3. PDF Shape Control Based on Hybrid Logistic Function Method

The PDF of the state variable corresponds to the steady-state solution of the FPK equation. However, the solution for the FPK equation is very difficult, and only under certain conditions can the exact solution be obtained. Therefore, we suggest a probability density function as the approximate solution of the FPK equation, which would be used to approximate a variety of target probability density functions p d x . In this paper, the hybrid logistic probability density function was used as the approximate solution of the FPK equation to obtain the corresponding controller so that the PDF of the state variable tracks the target probability density function p d x .
Logistic distribution is a continuous distribution defined by its location and scale parameters. The logistic distribution has a similar shape to the Gaussian distribution but with a longer tail. If the target PDF has a longer tail and higher crest, the hybrid logistic probability density function is more suitable as the approximate PDF of the state variable than the MGC method.
The logistic distribution function is:
F x = 1 + e x μ σ 1 , < x < +
where μ   < μ < + is the location parameter and σ   σ > 0 is the scale parameter.
The corresponding probability density function is:
f x = 1 σ e x μ σ 1 + e x μ σ 2
The following results are obtained by hybridizing the probability density functions of N logistic functions:
p x = i = 1 N w i 1 σ i e x μ i σ i 1 + e x μ i σ i 2
where w i is the weight coefficient of the ith logistic function, such that i = 1 N w i = 1 ; μ i is the location parameter of the ith logistic function; and σ i is the scale parameter of the ith logistic function.
In order for the hybrid logistic probability density function p x to track the target probability density function well, we constructed the following optimization problem:
min e = x d x D p d x p x 2 d x s . t . i = 1 N w i = 1 σ i > 0 < μ i < +
where e is the residual sum of squares of the target PDF and tracking PDF (i.e., the hybrid logistic probability density function), x D is the maximum value of the state variable x, and x d is the minimum value of the state variable x. Meanwhile, e is used to calculate the control error.

4. Solving the PDF Shape Controller

For a stochastic non-linear system, the FPK equation can fully describe the dynamic evolution law of the probability density function γ x , t of the state x t in time and space. Its specific form is:
γ x , t t = γ x , t Φ x x + 1 2 2 γ x , t 2 π S 0 x 2
where S 0 is the power spectral density of Gaussian white noise and Φ x is the controlled function of the state variable.
When the system response tends to be stable, the PDF of the steady-state response is independent of the initial conditions and time and satisfies the following steady-state FPK equation:
0 = d γ x Φ x d x + 1 2 d 2 γ x 2 π S 0 d x 2
Based on the control system (2), Φ x in Equations (7) and (8) should meet the relationship:
Φ x = y x + u x
Let π S 0 = H . Since S 0 is a constant, H is also a constant. When the hybrid logistic probability density function p x is regarded as an approximate solution of the FPK equation, Equation (7) can be transformed into:
p x , t t = p x , t Φ x x + H 2 p x , t x 2
Therefore, the steady-state FPK equation satisfying the hybrid logistic probability density function p x is:
0 = d p x Φ x d x + H d 2 p x d x 2
In Hilbert space, we usually define the inner product as:
q x , g x R n q x g x d x q , g H
Thus, we multiply both sides of Equation (11) by the test functions h l l = 1 L : h l x l l = 1 , 2 , , L , and integrate both sides over the entire region simultaneously to obtain the following equation:
d p x Φ x d x + H d 2 p x d x 2 , h l x = R d p x Φ x d x + H d 2 p x d x 2 h l x d x = 0
The introduction of the test function h l l = 1 L weakly solves the above integral so that the hybrid logistic function p ( x ) acts as the target probability density function p d ( x ) in the weak sense and the error between them vanishes.
We choose h l x l l = 1 , 2 , , L , and obtain the following equation set from (13):
R d p x Φ x d x + H d 2 p x d x 2 x d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x 2 d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x L d x = 0
We solve the differential Equation (14) and obtain all of the general solutions Φ x . Then, the controller function u ( x ) is obtained from Equation (9) as:
u x = Φ x y x
Since the controller function structure is unknown, the structure of the function Φ x is not uniquely determined, which, in turn, increases the difficulty in solving the equation. Therefore, it is not easy to obtain the general solution of Φ x directly. In summary, we will make assumptions on the controller function structure so as to complete solving the controller.
If the non-linear function of the system is polynomial, the expression is:
y x = a 0 + a 1 x + + a n x n = Σ n j = 0 a j x j
We adopt common non-linear controllers given by the expression:
u x = b 0 + b 1 x + + b m x m = Σ m k = 0 b k x k
Therefore, we substitute the non-linear function (16) and the controller function (17) into Equation (9) and obtain the following expression:
Φ ( x ) = a 0 + a 1 x + + a n x n + b 0 + b 1 x + + b m x m
The order of the controller u ( x ) is determined by the control requirements of the system. For ease of calculation, we assume that m n . Then, we organize:
Φ x = j = 0 m a j + b j x j m = n Φ x = j = 0 n a j + b j x j + b n + 1 x n + 1 + + b m x m m > n
where a j is known and b j is unknown.
We take Equation (19) and the best approximation hybrid logistic function p ( x ) into Equation (14) and choose h l = x l to solve for the unknowns in Equation (19). Then, the final solution is obtained for the controller u ( x ) :
u x = b 0 ^ + b 1 ^ x + b 2 ^ x 2 + + b m ^ x m = j = 0 m b j ^ x j
where b j ^ is the optimal solution of the controller parameters
If the system function y ( x ) is in other forms, we need to make reasonable assumptions about the controller function according to the design requirements or the structure of the system equations, and then solve the controller in the same way. Therefore, this method is not limited to polynomial systems, but is applicable to any non-linear stochastic system.

5. Simulations

The proposed algorithm for the design of probability density function shape controller for nonlinear stochastic systems based on hybrid logistic functions is as follows:
Step 1. Specify N in the hybrid logistic function according to the shape of the target probability density function;
Step 2. Solve optimization problem (6), obtain the values of parameters w i , μ i , and σ i , and calculate e;
Step 3. Judge whether e meets the precision requirement. If so, go to Step 4; otherwise, let N = N + 1 and return to Step 2;
Step 4. If the system has a polynomial structure, design the controller u ( x ) according to Equation (17) and then bring it into Equation (19); otherwise, design the controller u ( x ) according to the system control requirements, and then bring it into (9), thus determining the structure of Φ x .
Step 5. Calculate Φ x and choose a suitable value of L in h l x l l = 1 , 2 , , L . Determine the selection of L by the number of unknown parameters in the controller u ( x ) . Then, substitute Φ x and h l x l into (14) together with the best approximate hybrid logistic function to obtain the controller parameters.
We consider two single-degree-of-freedom stochastic system models with Gaussian white noise. The first model is the polynomial system model and the second model is the non-polynomial system model.
For the polynomial system model:
d x d t = x   x 2 0.5 x 3 + w t x t 0 = x 0
For the non-polynomial system model:
d x d t = sin ( x ) + w t x t 0 = x 0
For computational simplicity, the spectral density of w t is taken as S 0 = 1 π , and x 0 = 1 is taken to be the initial state of the system.
To evaluate the performance of the PDF shape control method based on the hybrid logistic probability density function, we conducted simulation experiments. We used the probability density function with three kinds (unimodal, bimodal, and trimodal) as the target PDFs and used MGC, EP, and HLF methods to track the target PDFs for different shapes.
(1) Case 1
The target PDF Gaussian probability density function has a unimodal symmetric shape and is given by the mathematical expression:
p d x = 1 2 π σ e x μ 2 2 σ 2
where σ = 0.5 and μ = 1 .
The target PDF is first approximated using a hybrid logistic probability density function. Here, we used N = 1 and applied the non-linear least squares method to obtain the hybrid logistic PDF p x using the formula:
p x = 1 0.307 e x 1 0.307 1 + e x 1 0.307 2
In the the polynomial system model:
From (21), we know that:
y x = x x 2 0.5 x 3
Thus, we assume that the functional expression for the controller is:
u x = b 0 + b 1 x + b 2 x 2 + b 3 x 3
Substitute (25) and (26) into (19), generating:
Φ x = b 0 + ( b 1 + 1 ) x + ( b 2 1 ) x 2 + ( b 3 0.5 ) x 3
Substituting (24) and (27) into (14) and selecting h l x l l = 1 , 2 , 3 , 4 generates:
R d p x Φ x d x + H d 2 p x d x 2 x d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x 2 d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x 3 d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x 4 d x = 0
Solving the above equation yields:
b 0 ^ = 3.39911 , b 1 ^ = 3.17513 , b 2 ^ = 0.835433 , b 3 ^ = 1.11182
Therefore, the optimal controller is:
u x = 3.39911 3.17513 x 0.835433 x 2 + 1.11182 x 3
In the non-polynomial system model:
From (22), we know:
y ( x ) = s i n ( x )
Thus, we assume that the functional expression for the controller is:
u x = b 0 + b 1 x
Then, we substitute the above equation into (9), obtaining:
Φ ( x ) = s i n ( x ) + b 0 + b 1 x
Substituting (24) and (32) into (14) and selecting h l x l l = 1 , 2 generates:
R d p x Φ x d x + H d 2 p x d x 2 x d x = 0 R d p x Φ x d x + H d 2 p x d x 2 x 2 d x = 0
Solving the above equation yields:
b 0 ^ = 2.94417 , b 1 ^ = 3.66833
Then, the optimal controller is:
u x = 2.94417 3.66833 x
When the target PDF has a symmetric unimodal shape, the simulation results are as shown in Figure 1. Figure 1a is the PDF shape control effect of the state variable for the three methods, and Figure 1b is the tracking error.
(2) Case 2
The target PDF beta probability density function has an asymmetric unimodal shape. The mathematical expression is:
p d ( x ) = 1 B 2 , 11 x 1 x 10 , x 0 , 1
where B is the beta function for normalization purposes, such that the integral is 1.
In the hybrid logistic probability density function, taking N = 2 , the optimum is determined by the expression:
p x = 0.5303 × 1 0.02275 e x + 0.03148 0.02275 1 + e x + 0.03148 0.02275 2 + 1.5303 × 1 0.08086 e x 0.07954 0.08086 1 + e x 0.07954 0.08086 2
In the the polynomial system model:
Select h l x l l = 1 , 2 , 3 , 4 , combine Equations (27), (28) and (36), and solve b j ^ in the controller:
b 0 ^ = 43.4223 , b 1 ^ = 570.292 , b 2 ^ = 1776.84 , b 3 ^ = 1633.25
The optimal controller is:
u x = 43.4223 570.292 x + 1776.84 x 2 1633.25 x 3
In the non-polynomial system model:
Select h l x l l = 1 , 2 , combine Equations (32), (33) and (36), and solve b j ^ in the controller:
b 0 ^ = 14.8129 , b 1 ^ = 95.432
The optimal controller is:
u x = 14.8129 95.432 x
When the target PDF has an asymmetric unimodal shape, the control effect and the tracking error are as shown in Figure 2.
(3) Case 3
In Cases 1 and 2, the target PDF only behaves unimodal, whereas the actual system is generally required to be in a multimodal shape. For example, the following target PDF is bimodal in shape:
p d x = 1.1284 x 2 e x 2
We take N = 4 in the hybrid logistic probability density function and obtain the optimal solution as:
p x = 0.5794 × 1 0.3427 e x + 1.0038 0.3427 1 + e x + 1.0038 0.3427 2 + 0.1396 × 1 0.2660 e x 1.4310 0.2660 1 + e x 1.4310 0.2660 2 + 0.4075 × 1 0.2879 e x 0.8979 0.2879 1 + e x 0.8979 0.2879 2 0.1265 × 1 0.2171 e x + 0.0155 0.2171 1 + e x + 0.0155 0.2171 2
In the the polynomial system model:
Select h l x l l = 1 , 2 , 3 , 4 , combine Equations (27), (28) and (40), and solve b j ^ in the controller:
b 0 ^ = 0.145343 , b 1 ^ = 0.834639 , b 2 ^ = 0.88514 , b 3 ^ = 0.200824
The optimal controller is:
u x = 0.145343 0.834639 x + 0.88514 x 2 0.200824 x 3
In the non-polynomial system model:
Select h l x l l = 1 , 2 , combine Equations (32), (33) and (40), and solve b j ^ in the controller:
b 0 ^ = 0.0163366 , b 1 ^ = 1.27504
The optimal controller is:
u x = 0.0163366 1.27504 x
When the target PDF has a symmetric bimodal shape, the simulation results are as shown in Figure 3, including the control effect and the error of the three methods on controlling the PDF shape of the state variable.
(4) Case 4
When we choose the target PDF to have a complex trimodal shape:
p d x = 0.1659 e x 2 5 x 4 3 + 0.2551 e 7 x x 2 2 49 2
We take N = 4 in the hybrid logistic probability density function and obtain the optimal as:
p x = 0.6391 × 1 0.6145 e x + 7 0.6145 1 + e x + 7 0.6145 2 + 0.1449 × 1 0.1996 e x + 0.5892 0.1996 1 + e x + 0.5892 0.1996 2 + 1 0.1950 e x 0.3703 0.1950 1 + e x 0.3703 0.1950 2 0.7840 × 1 0.1757 e x 0.3619 0.1757 1 + e x 0.3619 0.1757 2
In the the polynomial system model:
Select h l x l l = 1 , 2 , 3 , 4 , combine Equations (27), (28) and (44), and solve b j ^ in the controller:
b 0 ^ = 0.0117224 , b 1 ^ = 2.33374 , b 2 ^ = 0.563835 , b 3 ^ = 0.465683
The optimal controller is:
u x = 0.0117224 2.33374 x + 0.563835 x 2 0.465683 x 3
In the non-polynomial system model:
Select h l x l l = 1 , 2 , combine Equations (32), (33) and (40), and solve b j ^ in the controller:
b 0 ^ = 0.496586 , b 1 ^ = 0.164573
The optimal controller is:
u x = 0.496586 0.164573 x
When the target PDF has a complex trimodal shape, the simulation results are as shown in Figure 4.
For the above four cases, in order to more clearly analyze and compare the PDF shape control performance of the three methods, we provide the cumulative distribution function (CDF) curves of the PDF from the three methods and the target PDF, as shown in Figure 5a–d.
Table 1 and Table 2 shows the calculated mean and variance using the different methods for the different target PDFs, while Table 3 summarizes the control errors.

6. Result Analysis

As presented in the simulation results in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5, the hybrid logistic probability density function method can better control the PDF shape of the state variables for the non-linear stochastic systems regardless of the target PDF being unimodal, bimodal, and complex trimodal- shaped.
When the target PDF has a Gaussian distribution, the MGC method and EP method have a better control effect than our method. However, the PDF shape obtained from our proposed approach matches the target PDF shape. In Figure 2a, our method outperforms the other two methods when the target PDF has a beta distribution. The EP method came second, and the MGC method performed the worst. In Figure 3a, when the target PDF is bimodal in shape and an equal number of parameters is used in the three methods, the control effect of the EP method is significantly inferior compared to the other two, with the MGC method performing slightly better than our method. However, if the number of parameters of the EP method is increased to 21, its control effect would be similar to the other two. When the target PDF has a complex trimodal shape, our method outperforms the others. The PDF shape of the state variable obtained by our proposed approach has the same trimodal shape as the target PDF, whereas the other two only showed bimodal forms. For this scenario, the MGC method and the EP method fail to achieve PDF shape control.
We analyze Figure 1b, Figure 2b, Figure 3b and Figure 4b and find that, when the target PDF has a Gaussian distribution, the tracking error of the EP method and MGC method is smaller than that of the HLF method. However, when the target has a beta distribution, the tracking error value of our method is significantly better than that of the EP method and the MGC method, and its error fluctuation is also much smaller. For the tracking error of a bimodal target PDF, the MGC method is the smallest, the HLF method is the second, and the EP method is the largest. If the target is a complex multimodal shape, the EP method has the largest fluctuation range in terms of the tracking error, and is therefore inferior to the MGC method and our method. Overall, when the target PDF is asymmetric or has a more complex trimodal shape, the resulting PDF from the HLF method tracks the target PDF significantly better than the other two methods.
The CDF curves in Figure 5a show that the EP method and the MGC method have slightly better control effects than our method, but the differences between the three methods are small. The CDF curves in Figure 5b,d also demonstrate that our proposed method has a more obvious advantage over the other two methods in terms of PDF shape performances. The CDF curves in Figure 5c show that the MGC method works best, our method is second, and the EP method is last.
Analyzing the mean and variance of the PDFs in Table 1 and Table 2, in the first case, the mean values calculated from three methods and the target PDF are exactly the same. The variance values obtained by the MGC method and the EP method are equal to that of the target PDF, while there exists an error in the variance calculated by our method. This is because the expression of the state PDF in the two former methods is of exponential form, which is the same as the target PDF. In the second case, the mean and variance calculated by our method are closest to the mean and variance of the target PDF, which is consistent with the previous analysis. In the third case, the mean and variance computed by the MGC method are closer to the target PDF than those calculated by the other two approaches. In the fourth case, the mean and variance calculated by our proposed method are closest to the target PDF, with the MGC method being the next most effective and the EP method being the least effective. In contrast, the PDF obtained by our approach is trimodal, which is closest to the shape of the target PDF.
From the control errors in Table 3, the EP method and the MGC method in the first case both have an error of 0, whereas the control error of our method is relatively larger. In the second and fourth cases, our method’s control error is significantly smaller than the other two approaches. In the third case, the MGC method has the smallest control error, followed by our proposed technique, and the EP method has the largest control error. However, the degrees of the control errors from our method are all no more than 10 2 for all four cases.
Based on the above analysis, the MGC method is more suitable for cases when the target PDF is symmetrical (the first and third cases). The EP method would require more parameters to be determined in order to obtain good control results. For the first scenario, the EP method required 3 parameters, whereas the other two approaches only needed 2. For the third scenario, when the same control effect was obtained, the EP method required 21 parameters, whereas the other two only required 11. Our method can achieve better control effects for various shapes of the target PDF. The focus of our control is the PDF shape of the state variable. In Figure 1, Figure 2, Figure 3 and Figure 4, although our method may not always yield the smallest errors, it ensures that the PDF shape of the state variable is consistent with the target PDF, thereby achieving the purpose of PDF shape control.

7. Conclusions

For non-linear stochastic systems, the state variables have random uncertainties, and the PDF shape of random variables completely characterizes the statistical properties of uncertain systems. In this paper, an innovative approach for controlling the PDF shape of non-linear stochastic systems was proposed. The method uses the hybrid logistic probability density function as an approximate target-tracking PDF of the state variable. After introducing the definition of inner product of Hilbert space, the hybrid logistic function is multiplied with the test function and integrated over the whole space, and then the optimal controller is obtained by solving the equation. This approach is applicable to arbitrary non-linear stochastic systems.
The simulation results show that the HLF method can achieve the purpose of PDF shape control regardless of if the target PDF shape is unimodal, bimodal, or trimodal. In particular, when the target PDF has an asymmetric distribution or complex trimodal shape, the proposed approach presents noticeable advantages over other techniques. Moreover, the method requires fewer parameters than the EP method when achieving the same control effect. For solving the controller, the solution is more difficult when the controller order is higher. Overall, the controller solution of polynomial systems is simpler than the non-polynomial system.
Although the proposed method can achieve good PDF control effects, the primary challenge is determining the optimization algorithm to obtain the global optimal solution for the parameters. In addition, the applicability of the method should be further evaluated for systems suffering from non-Gaussian white noise. These considerations will be explored and addressed in our subsequent studies.

Author Contributions

L.W.: project administration, methodology, supervision, funding acquisition; K.Z.: conceptualization, data curation, formal analysis, investigation, writing—original draft, writing—review and editing; F.Q.: validation, funding acquisition; X.Z.: funding acquisition. 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 (61903298, 62073259), Natural Science Foundation of Shaanxi Provincial Department of Education (20JS130), and Xi’an Science and Technology Planning Project (21XJZZ0060).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PDFprobability density function
HLFhybrid logistic function
FPKFokker–Planck–Kolmogorov
MGCmulti-Gaussian closure
EPexponential polynomial
CDFcumulative distribution function

References

  1. Sain, M. Control of linear systems according to the minimal variance criterion: A new approach to the disturbance problem. IEEE Trans. Autom. Control 1966, 11, 118–122. [Google Scholar] [CrossRef]
  2. Anderson, B.; Moore, J. Optimal Control, Linear Quadratic Methods; Prentice-Hall: Englewood Cliffs, NJ, USA, 1990. [Google Scholar]
  3. Paola, M.; Ricciardi, G.; Vasta, M. A method for the probabilistic analysis of non-linear systems. Probabilistic Eng. Mech. 1995, 10, 1–10. [Google Scholar] [CrossRef]
  4. Karny, M. Towards fully probabilistic control design. Automatica 1996, 32, 1719–1722. [Google Scholar] [CrossRef]
  5. Li, G.; Wang, Y.; Zeng, Y.; He, W. A new maximum entropy method for estimation of multimodal probability density function. Appl. Math. Model. 2022, 102, 137–152. [Google Scholar] [CrossRef]
  6. Wang, H. Robust control of the output probability density functions for multivariable stochastic systems. In Proceedings of the 37th IEEE Conference on Decision and Control, Tampa, FL, USA, 18 December 1998; pp. 1305–1310. [Google Scholar]
  7. Wang, H. Robust control of the output probability density functions for multivariable stochastic systems with guaranteed stability. IEEE Trans. Autom. Control 1999, 44, 2103–2107. [Google Scholar] [CrossRef]
  8. Wang, H. Control for bounded pseudo ARMAX stochastic systems via linear B-spline approximations. IEEE Conf. Decis. Control 2000, 4, 3369–3374. [Google Scholar]
  9. Wang, H.; Zhang, J. Bounded stochastic distributions control for pseudo-ARMAX stochastic systems. IEEE Trans. Autom. Control 2001, 46, 486–490. [Google Scholar] [CrossRef]
  10. Wang, H. Minimum entropy control of non-Gaussian dynamic stochastic systems. IEEE Trans. Autom. Control 2002, 47, 398–403. [Google Scholar] [CrossRef]
  11. Wang, H.; Sun, X. Neural network based probability density function shape control for unknown stochastic systems. In Proceedings of the 2004 IEEE International Symposium on Intelligent Control, Taipei, Taiwan, 2–4 September 2004; pp. 120–125. [Google Scholar]
  12. Wang, H.; Zhang, J.; Yue, H. Multi-step Predictive Control of a PDF-shaping Problem. Acta Autom. Sin. 2005, 31, 274–279. [Google Scholar]
  13. Guo, L.; Wang, H. PID controller design for output PDFs of stochastic systems using linear matrix inequalities. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2005, 35, 65–71. [Google Scholar]
  14. Luan, X.; Liu, F. Robust Resilient Optimal Tracking Control for Output Probability Design Function. Control Eng. China 2008, 15, 493–496. [Google Scholar]
  15. Luan, X.; Liu, F. Finite time stabilization for output probability density function of stochastic system. Control Decis. 2009, 24, 1161–1166. [Google Scholar]
  16. Wang, H. Afshar Puya. ILC-Based Fixed-Structure Controller Design for Output PDF Shaping in Stochastic systems Using LMI Techniques. IEEE Trans. Autom. 2009, 54, 760–773. [Google Scholar] [CrossRef]
  17. Qian, F.; Gao, J.; Li, D. Complete statistical characterization of discrete-time LQG and cumulant control. IEEE Trans. Autom. Control 2012, 57, 2110–2115. [Google Scholar] [CrossRef]
  18. Mao, J.; Xiang, Z.; Zhai, G. Sampled-data output feedback stabilization for a class of switched stochastic nonlinear systems. Int. J. Robust Nonlinear Control 2019, 29, 2844–2861. [Google Scholar] [CrossRef]
  19. Kim, H.; Guo, P.; Zhu, M.; Liu, P. Simultaneous input and state estimation for stochastic nonlinear systems with additive unknown inputs. Automatica 2020, 111, 108588. [Google Scholar] [CrossRef]
  20. Zhang, T.; Deng, F.; Shi, P. Event-triggered H∞ filtering for nonlinear discrete-time stochastic systems with application to vehicle roll stability control. Int. J. Robust Nonlinear Control 2020, 30, 8430–8448. [Google Scholar] [CrossRef]
  21. Chen, J.; Liu, S.; Li, J.; Xie, J. Adaptive neural tracking control for stochastic nonlinear multi-agent periodic time-varying systems. Appl. Math. Model. 2022, 102, 228–242. [Google Scholar] [CrossRef]
  22. Yuan, T.; Yang, H.; Ivanov, I.G. Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives. Mathematics 2021, 9, 2856. [Google Scholar] [CrossRef]
  23. Drăgan, V.; Ivanov, I.G.; Popa, I.L. A Game—Theoretic Model for a Stochastic Linear Quadratic Tracking Problem. Axioms 2023, 12, 76. [Google Scholar] [CrossRef]
  24. Hu, Z.; Mu, X. Event-Triggered Impulsive Control for Nonlinear Stochastic Systems. IEEE Trans. Cybern. 2022, 52, 7805–7813. [Google Scholar] [CrossRef] [PubMed]
  25. Ji, R.; Liang, Y.; Xu, L.; Wei, Z. Estimation of dual-mode nonlinear stochastic systems with unknown parameters. Int. J. Robust Nonlinear Control 2022, 32, 9258–9274. [Google Scholar] [CrossRef]
  26. Yin, J.; Khoo, S.; Yu, X. Domain stabilization in probability in a fixed time for nonlinear stochastic systems via feedback control. Int. J. Robust Nonlinear Control 2023, 33, 1738–1749. [Google Scholar] [CrossRef]
  27. Gashi, B.; Hua, H. Optimal regulators for a class of nonlinear stochastic systems. Int. J. Robust Nonlinear Control 2023, 96, 136–146. [Google Scholar] [CrossRef]
  28. Forbes, M.; Guay, M.; Forbes, J. Control design for first-order processes: Shaping the probability density of the process state. J. Process Control 2004, 14, 399–410. [Google Scholar] [CrossRef]
  29. Forbes, M.; Forbes, J.; Guay, M. Control design for discrete-time stochastic non-linear processes with a non-quadratic performance objective. In Proceedings of the 42nd IEEE International Conference on Decision and Control, Maui, HI, USA, 9–12 December 2003; pp. 4243–4248. [Google Scholar]
  30. Forbes, M. Performance Characterization and Regulatory Feedback Control Design for Time-Invariant Discrete-Time Stochastic Processes. Ph.D. Thesis, International University of Alberta, Edmonton, AB, Canada, 2003. [Google Scholar]
  31. Michael, G.; Martin, G.; Forbes, J. Probabilistic control design for continuous time stochastic non-linear systems: A PDF shaping approach. In Proceedings of the 2004 IEEE International Symposium on Intelligent Control, Taipei, Taiwan, 2–4 September 2004; pp. 132–136. [Google Scholar]
  32. Guo, L.; Wang, H.; Wang, A. Optimal probability density function control for NARMAX stochastic Systems. Automatica 2008, 44, 1904–1911. [Google Scholar] [CrossRef]
  33. Zhu, C.; Yang, L. Target Control Design for Stationary Probability Density Function of Nonlinear Stochastic System. Acta Autom. Sinca 2012, 38, 197–205. [Google Scholar] [CrossRef]
  34. Zhao, Y.; Huang, J.; Qin, X.; Qian, F. The study of stochastic non-linear state variable density function control method. J. Xi’an Univ. Technol. 2013, 29, 144–148. [Google Scholar]
  35. Yang, H.; Qian, F.; Huang, J.; Gao, S. The complete statistical characterization control for a class of stochastic systems. Control Theory Appl. 2016, 33, 669–675. [Google Scholar]
  36. Yang, H.; Qian, F.; Gao, S.; Jiang, T. The shape control of probability density function for a class of stochastic systems. Syst. Eng. Theory Pract. 2016, 36, 2424–2431. [Google Scholar]
  37. Yang, H.; Fu, Y.; Gao, S.; Qian, F. PDF control of non-linear stochastic systems based on MGC method. Control Decis. 2019, 34, 1463–1468. [Google Scholar]
  38. Wang, L.; Qian, F.; Liu, J. Shape Control on Probability Density Function in Stochastic Systems. J. Syst. Eng. Electron. 2014, 25, 144–149. [Google Scholar] [CrossRef]
  39. Wang, L.; Qian, F.; Liu, J. The PDF Shape Control of the State Variable for a Class of Stochastic Systems. Int. J. Syst. Sci. 2015, 46, 2231–2239. [Google Scholar] [CrossRef]
  40. Wang, L.; Qian, F. Shaping PDF of the state variable based on piecewise Linear control for non-linear stochastic systems. Sci. China Inf. Sci. 2016, 59, 106–116. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Case 1: unimodal target PDF (Gaussian distribution). (a) Control effect. (b) Tracking error.
Figure 1. Case 1: unimodal target PDF (Gaussian distribution). (a) Control effect. (b) Tracking error.
Axioms 12 00303 g001
Figure 2. Case 2: unimodal target PDF (beta distribution). (a) Control effect. (b) Tracking error.
Figure 2. Case 2: unimodal target PDF (beta distribution). (a) Control effect. (b) Tracking error.
Axioms 12 00303 g002
Figure 3. Case 3: Bimodal target PDF. (a) Control effect. (b) Tracking error.
Figure 3. Case 3: Bimodal target PDF. (a) Control effect. (b) Tracking error.
Axioms 12 00303 g003
Figure 4. Case 4: Complex trimodal target PDF. (a) Control effect. (b) Tracking error.
Figure 4. Case 4: Complex trimodal target PDF. (a) Control effect. (b) Tracking error.
Axioms 12 00303 g004
Figure 5. CDFs of the three methods. (a) Case 1. (b) Case 2. (c) Case 3. (d) Case 4.
Figure 5. CDFs of the three methods. (a) Case 1. (b) Case 2. (c) Case 3. (d) Case 4.
Axioms 12 00303 g005
Table 1. Means of the PDFs based on the three methods.
Table 1. Means of the PDFs based on the three methods.
Case 1Case 2Case 3Case 4
Target PDF 1.0 0.1538 0 4.4761
PDF-EP 1.0 0.1458 0.0058 0
PDF-MGC 1.0 0.1417 0 4.3029
PDF-HLF 1.0 0.1583 0.0140 4.4725
Table 2. Variances of the PDFs based on the three methods.
Table 2. Variances of the PDFs based on the three methods.
Case 1Case 2Case 3Case 4
Target PDF 0.25 0.0093 1.5 12.0680
PDF-EP 0.25 0.0080 1.3934 0.1315
PDF-MGC 0.25 0.0079 1.5000 15.1631
PDF-HLF 0.31 0.0105 1.5458 12.2556
Table 3. Control errors of the three methods.
Table 3. Control errors of the three methods.
Case 1Case 2Case 3Case 4
PDF-EP0 2.0933 0.5310 11.5344
PDF-MGC0 2.8212 8.5008 × 10 6 0.3236
PDF-HLF 0.01395 0.0446419 9.0914 × 10 4 0.0058
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.

Share and Cite

MDPI and ACS Style

Wang, L.; Zhang, K.; Qian, F.; Zhang, X. Research on PDF Shape Control for Nonlinear Stochastic System Using an Approximate Solution of FPK Equation. Axioms 2023, 12, 303. https://doi.org/10.3390/axioms12030303

AMA Style

Wang L, Zhang K, Qian F, Zhang X. Research on PDF Shape Control for Nonlinear Stochastic System Using an Approximate Solution of FPK Equation. Axioms. 2023; 12(3):303. https://doi.org/10.3390/axioms12030303

Chicago/Turabian Style

Wang, Lingzhi, Kun Zhang, Fucai Qian, and Xiaoli Zhang. 2023. "Research on PDF Shape Control for Nonlinear Stochastic System Using an Approximate Solution of FPK Equation" Axioms 12, no. 3: 303. https://doi.org/10.3390/axioms12030303

APA Style

Wang, L., Zhang, K., Qian, F., & Zhang, X. (2023). Research on PDF Shape Control for Nonlinear Stochastic System Using an Approximate Solution of FPK Equation. Axioms, 12(3), 303. https://doi.org/10.3390/axioms12030303

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