Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 7679165, 9 pages

http://dx.doi.org/10.1155/2016/7679165

## Optimal Preview Control for a Class of Linear Continuous Stochastic Control Systems in the Infinite Horizon

^{1}School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China^{2}Leeds Sustainability Institute, Leeds Beckett University, Leeds LS2 9EN, UK

Received 1 July 2016; Accepted 26 September 2016

Academic Editor: Andrzej Swierniak

Copyright © 2016 Jiang Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper discusses the optimal preview control problem for a class of linear continuous stochastic control systems in the infinite horizon, based on the augmented error system method. Firstly, an assistant system is designed and the state equation is translated to the assistant system. Then, an integrator is introduced to construct a stochastic augmented error system. As a result, the tracking problem is converted to a regulation problem. Secondly, the optimal regulator is solved based on dynamic programming principle for the stochastic system, and the optimal preview controller of the original system is obtained. Compared with the finite horizon, we simplify the performance index. We also study the stability of the stochastic augmented error system and design the observer for the original stochastic system. Finally, the simulation example shows the effectiveness of the conclusion in this paper.

#### 1. Introduction

Future reference signals or disturbance signals are known in certain circumstances. All the known future information can be utilized by preview control theory to improve the performance of the dynamic system. The original idea of preview control theory is that in order to minimize the error between the reference signals and the controlled terms, not only the past and present information but also the future information should be concentrated on [1–5]. An augmented error system is constructed while designing the optimal preview controller for the discrete system. And also a group of identical equations of future reference signals is added to the augmented error system [6–9]. Since the continuous system is of infinite dimensions, the method in dealing with the discrete condition is no longer useful. In [10], the augmented error system was constructed by differentiating the state equation on both sides, combining with the error equation. According to the maximum principle, the optimal preview controller was obtained by solving a differential equation on reference signals backward in time. This method was extended to systems with previewable disturbance signals in [11] and to singular continuous systems in [12]. The application of preview control method to continuous systems with time delay is studied in [13].

According to automatic control system theory, the controlled systems can be regarded as falling into two categories: deterministic systems and stochastic systems. Stochastic systems are a collection of dynamic systems that contain internal stochastic parameters, external stochastic disturbances, or observation noises [14]. A typical stochastic system is the stochastic differential equation. It is a class of differential equations driven by one or more stochastic processes. Therefore, the solution of a stochastic differential equation is also a stochastic process. Up to the present, a phenomenon such as stock market volatility or thermal motion in a physical system is usually described by a stochastic differential equation. Typically, a white noise stochastic variable represented by the differential form of Brownian motion or Weiner process is contained in a stochastic differential equation.

In this paper, preview control theory is applied to a class of linear continuous stochastic control systems in the infinite horizon. In the finite horizon [15], an assistant system is designed, and the state equation is translated to the assistant system. Then an integrator is introduced to construct a stochastic augmented error system. As a result, the tracking problem is converted to a regulation problem. The optimal regulator is obtained based on the dynamic programming principle for stochastic systems, which means the optimal preview controller of the original stochastic system is gained. For such a system, when the integrator approaches zero at infinity, the error also approaches zero. Therefore, this property can be utilized to simplify the performance index in the infinite horizon. Due to the fact that the relative terms of reference signals are included in the stochastic augmented error system, the conclusion in the infinite horizon cannot be directly employed when solving the optimal regulation problem in this paper. Firstly, the corresponding optimal regulation problem is solved in the finite horizon with the new performance index and then the time was set to approach infinity. The stability of the stochastic augmented error system is studied and the sufficient and necessary criteria which guarantee that there exists a unique semipositive definite solution to the corresponding Riccati equation have been met. Also an observer for the original stochastic system is designed. The introduction of the integrator can eliminate the static error and the simulation example shows the effectiveness of the conclusion in this paper [16].

*Notation*. The notations are standard. is a complete probability space. Adaptive procedure is algebra generated by , and denotes Brownian motion of dimensions. denotes matrix. denotes a positive definite (semidefinite) matrix. denotes the unit matrix. The symbol denotes the symmetric terms in a symmetric matrix. denotes the expectation of process with initial time and state . denotes the trace of a matrix.

#### 2. Problem Statement

Consider the following stochastic control system: where is the state vector of dimensions, is the output vector of dimensions, is the adaptive input of dimensions, and is the Brownian motion of adaptive of dimensions. , , and are constant matrices. .

First of all, the following assumptions are introduced.

*Assumption 1. *Suppose the matrices pair is stabilizable.

*Assumption 2. *Suppose the matrices pair is detectable and the matrix is of full row rank.

By denoting the reference signals of dimensions with , another assumption about is needed.

*Assumption 3. *Suppose the reference signal is a piecewise continuously differentiable function defined on . Moreover, () is available at any moment and when , ; that is, . is the preview length of the reference signal.

*Remark 4. *Since there are few effects on the system when the reference signals exceed the previewable range, it is always assumed that the signals are usually constant [11].

The tracking error signal is defined as the difference between the output vector and the reference signal ; that is,

In this paper, an optimal controller will be designed to make the output in (1) tracking the reference signal as accurately as possible without static error.

Then, an assistant system is defined. The following assumption is needed.

*Assumption 5. *Suppose

With Assumption 5, the following lemma holds.

Lemma 6 (see [15]). *There exist matrices pairs which solve the following:if and only if Assumption 5 holds, where and .*

*Definition 7 (see [17]). *If matrix and the rank of is equal to , then has a left inverse: such that .

*Remark 8. *If , there exists unique solution of matrices pair , while if , according to Definition 7, any left inverse of matrix can be applied to calculate the solution of (4) and the solutions are infinite. In this paper only the square full-rank case is used.

Let the state vector and input vector be denoted by and , respectively. Let and . According to Lemma 6, at the moment , for any time , meetsSince and according to (4), substituting these two equations into (5) yieldsthat is,DefineSubtracting (7) by (1) from both sides, the following could be obtained:Now an integrator is introduced. Integrating the tracking error signal in the interval of to yieldsObviously . Let , and then holds. Define and the differential form of becomes

*Remark 9. *According to Assumption 3, the term in (11) is also previewable while . In particular, while , .

In order to make the output tracking the reference signal as accurately as possible, the performance index is designed aswhere and are matrices of proper dimensions.

*Remark 10. *Compared with the finite horizon, can also be obtained in the infinite horizon while . Therefore, the term in the performance index of [15] can be replaced by in (12). Furthermore, since there exists unique semipositive definite solution of the corresponding Riccati equation when , no other terms need to be included in the performance index, which simplifies the calculation.

To solve the optimal control problem of (9) with the performance index (12), the augmented error system method can be employed.

#### 3. Construction of the Stochastic Augmented Error System

Combining (9) and (11), the following holds:whereReferring to the output equation of (9) and the reference signal, the output of (13) can be defined aswhere . Joining (13) and (15) together, the following hold:

Equation (16) is the needed stochastic augmented error system.

#### 4. Design of the Optimal Controller for the Stochastic System

Denoting the performance index (12) with the related variables in (16) yields

Due to the existence of the term in (16), the conclusion of the optimal regulation problem in dealing with the infinite horizon cannot be directly employed. Therefore, similar to the condition in deterministic systems [11], the optimal regulator of the stochastic augmented error system (16) can be obtained through the following three steps.

Firstly, (17) is revised to

Secondly, the finite-time horizon optimal regulation problem with the performance indexis solved, where is the terminal time.

Compared with the performance index in [15], the two terms and are removed in (19). Therefore, according to the dynamic programming principle for stochastic systems [16, 18], the following lemma about the finite horizon optimal control problem with the terminal time can be obtained by letting in Theorem 1 in [15].

Lemma 11 (see [15]). *If Assumptions 3 and 5 hold, the optimal controller of (16), with the performance index (19), iswhere satisfies the Riccati differential equationWith the terminal condition and the differential equation is obtained:with the terminal condition .*

Thirdly, Theorem 12 is received by letting *.*

Theorem 12. *If Assumptions 1 to 5 hold, the optimal controller of (16) with the performance index (17) can be expressed aswhere is the unique semipositive definite solution of Riccati differential equation*

*Proof. *In terms of the performance index (17), if is stabilizable and is detectable, according to the known conclusion, there exists a unique semipositive definite solution of Riccati algebra equation And when , As a result, the solving of in (22) is also simplified.

Denote . Equation (22) and its terminal condition can be expressed asBased on linear system theory [19], the fundamental solution matrix of the homogeneous system corresponding to (27) is with the propertyTherefore, the solution of (27) isand the backward form isWith the terminal condition , the following will be obtained:According to the property of the reference signals in Assumption 3, can be expressed asSubstituting (33) into (20) yieldswhere .

The conclusion is that when , the following theorem can be obtained.

Theorem 13. *If Assumptions 1 to 5 hold, the optimal controller of (16) with the performance index (17) can be expressed aswhere is the unique semipositive definite solution of Riccati differential equation**Then, let the matrix be separated intowhere and are matrices of and . According to (8), (35), , and , the following will be obtained:where satisfies Riccati differential equation (24).*

Due to the fact that , , and , the following corollary is obtained.

Corollary 14. *If Assumptions 1 to 5 hold, the optimal preview controller of (1) with the performance index (12) can be expressed as*

From (39) it can be seen that there are four parts contained in the preview controller of (1): the first part is the integral of tracking error term , stemmed from the inducing of the integrator. The second part is the state feedback . The third one is the reference preview compensation. The last one is the previewable complement of the assistant system.

#### 5. Stability of Closed-Loop System

With being an external disturbance, the stability and detectability of system (16) are studied only when . Here, (16) becomes the deterministic systemAs a result, the sufficient and necessary criteria which guarantee that there exists a unique semipositive definite solution to the Riccati equation (24) can be gained by the known conclusion, as follows.

Lemma 15 (see [11]). *Suppose the following conditions are satisfied: *(1)*The matrices and are both positive definite.*(2)*The matrix is of full row rank.*(3)*The matrices pair is stabilizable.*(4)*The matrices pair is detectable.**Then there exists unique semipositive definite solution in the Riccati equation (24), and there exists the asymptotically stable coefficient matrix in the closed-loop system (27).*

*Therefore, the following theorem can be received.*

*Theorem 16. If and are both positive definite matrices and Assumptions 1 to 5 hold, the full regulation is achieved in the closed-loop system:namely,Furthermore, the following will be obtained:*

*6. State Observer*

*If the state vector in (1) cannot be measured directly, the optimal preview controller (39) cannot be realized. In order to solve this problem, the state observer could be designed.*

*Subtracting (45) from (1) on both sides yieldswhere .*

*So, if is stable, the following will hold:which means that the state vector in the observer equation (45) approximates the state vector in (1) when . Based on linear system theory [19], it is known that if is detectable, is stable. Then the following will be obtained.*

*Theorem 17. If the conditions in Theorem 16 hold, the closed-loop system can be received:where is stable and (48) achieve the complete regulation.*

*7. Numerical Simulation*

*Example 1. *Consider the stochastic control systemwhere the coefficient matrices arerespectively. Therefore, the coefficient matrices in (16) areand Brownian motion satisfiesLet the initial state of (1) bewhere and are mutually independent. The preview lengths of the reference signals are , respectively. The weight matrices of the performance index are and the solution of Riccati equationiswhich is a positive definite matrix. The coefficient matrix in (25) isWhen the reference signal isFigures 1–3 can be obtained.

Figure 1 shows the tracking controller with nonpreview. Figure 2 shows the tracking controller with preview. Figure 3 shows the tracking controller with preview. Comparing the above three figures, it can be seen that based on the preview controller, the output signals can track the reference signals much faster and with less tracking error.

When the reference signal isFigures 4–6 can be obtained.

Figure 4 shows the tracking controller with nonpreview. Figure 5 shows the tracking controller with preview. Figure 6 shows the tracking controller with preview. Comparing the above three figures, it can be seen that based on the preview controller, the output signals can track the reference signals much faster and with less tracking error.