Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2016, Article ID 7679165, 9 pages http://dx.doi.org/10.1155/2016/7679165
Research Article Optimal Preview Control for a Class of Linear Continuous Stochastic Control Systems in the Infinite Horizon Jiang Wu,1 Fucheng Liao,1 and Jiamei Deng2 1
School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China Leeds Sustainability Institute, Leeds Beckett University, Leeds LS2 9EN, UK
2
Correspondence should be addressed to Fucheng Liao;
[email protected] 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. 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
2
Mathematical Problems in Engineering
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. π > 0 (π β₯ 0) denotes a positive definite (semidefinite) matrix. πΌ denotes the unit matrix. The symbol β denotes the symmetric terms in a symmetric matrix. πΈπ‘0 ,π₯0 denotes the expectation of process (π‘, π₯) with initial time π‘0 and state π₯0 . tr(β
) denotes the trace of a matrix.
2. Problem Statement
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, π (π‘) = π¦ (π‘) β π (π‘) .
(2)
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 π π ] = π + π, (full row rank) . πΆ 0
rank [
(3)
With Assumption 5, the following lemma holds. Lemma 6 (see [15]). There exist matrices pairs (Ξ, πΎ) which solve the following: πΞ + ππΎ = 0, πΆΞ = πΌ,
(4)
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 π΅π΄ = πΌ.
Consider the following stochastic control system: ππ₯ (π‘) = [ππ₯ (π‘) + ππ’ (π‘)] ππ‘ + π (π‘) ππ΅π‘ , π¦ (π‘) = πΆπ₯ (π‘) ,
Moreover, π(π) (π‘ β€ π β€ π‘ + ππ ) is available at any moment Μ = 0. ππ is π‘ and when π > π‘ + ππ , π(π) = π(π‘ + ππ ); that is, π(π) the preview length of the reference signal.
(1) π₯ (π‘0 ) = π₯0 , π‘ β [π‘0 , β) ,
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 (π‘ β₯ π‘0 ).
Remark 8. If π = π, there exists unique solution of matrices pair (Ξ, πΎ), while if π < π, according to Definition 7, any π left inverse of matrix [ π πΆ 0 ] 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 π‘ (π β€ π‘ β€ π‘π ), π₯β (π‘) meets ππ₯β (π‘) = ΞπΜ (π‘ + ππ ) ππ‘, πΆπ₯β (π‘) = πΆΞπ (π‘ + ππ ) .
(5)
Since [πΞ + ππΎ]π(π‘ + ππ ) = 0 and πΆΞπ(π‘ + ππ ) = π(π‘ + ππ ) according to (4), substituting these two equations into (5) yields ππ₯β (π‘) = {[πΞ + ππΎ] π (π‘ + ππ ) + ΞπΜ (π‘ + ππ )} ππ‘, π (π‘ + ππ ) = πΆπ₯β (π‘) ;
(6)
Mathematical Problems in Engineering
3
3. Construction of the Stochastic Augmented Error System
that is, β
β
β
ππ₯ (π‘) = {[ππ₯ (π‘) + ππ’ (π‘)] + ΞπΜ (π‘ + ππ )} ππ‘, π (π‘ + ππ ) = πΆπ₯β (π‘) .
(7)
Combining (9) and (11), the following holds: Μ (π‘) + πΜ Μ π’ (π‘) + π·Μ Μ π (π‘) + ΞΜπΜ (π‘ + ππ )] ππ‘ ππ§ (π‘) = [ππ§
Define β
Μ (π‘) = π₯ (π‘) β π₯ (π‘) , π₯ Μ (π‘) = π’ (π‘) β π’β (π‘) , π’ β
Μ (π‘) = π¦ (π‘) β πΆπ₯ (π‘) , π¦
+ π (π‘) ππ΅π‘ , (8)
where
], 0 π
0 Μ = [ ], π π
(9)
Μ=[ π·
Μ (π‘0 ) = π₯0 β π₯β (π‘0 ) , π‘ β [π‘0 , β) . π₯ Now an integrator is introduced. Integrating the tracking error signal π(π‘) in the interval of π‘0 to π‘ yields π‘
π‘0
π‘0
π (π‘) = β« π (π ) ππ = β« [π¦ (π ) β π (π )] ππ .
(10)
π‘
Obviously π(π‘0 ) = 0. Let πβ (π‘) = β«π‘ [πΆπ₯β (π ) β π(π + ππ )]ππ , and 0 then πβ (π‘) = 0 holds. Define πΜ(π‘) = π(π‘) β πβ (π‘) and the differential form of πΜ(π‘) becomes πΜπ (π‘) = π (π‘) ππ‘ = [πΆΜ π₯ (π‘) β Μπ (π‘)] ππ‘.
(11)
Remark 9. According to Assumption 3, the term Μπ(π‘) in (11) is also previewable while Μπ(π) (π‘ β€ π β€ π‘ + ππ ). In particular, while π β₯ π‘ + ππ , Μπ(π) = 0. In order to make the output π¦(π‘) tracking the reference signal π(π‘) as accurately as possible, the performance index π½ is designed as
ΞΜ = [
where ππΜ = ππΜπ > 0 and π
= π
π > 0 are matrices of proper dimensions. Remark 10. Compared with the finite horizon, limπ‘ββ π(π‘) = 0 can also be obtained in the infinite horizon while limπ‘ββ π(π‘) = 0. 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 ππΜ > 0, 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.
βπΌ 0
(14)
],
0
]. βΞ
Referring to the output equation of (9) and the reference signal, the output of (13) can be defined as Μ (π‘) , πΜ (π‘) = πΆπ§
(15)
Μ = [πΌ 0]. Joining (13) and (15) together, the following where πΆ hold: Μ (π‘) + πΜ Μ π’ (π‘) + π·Μ Μ π (π‘) + ΞΜπΜ (π‘ + ππ )] ππ‘ ππ§ (π‘) = [ππ§ + π (π‘) ππ΅π‘ , Μ (π‘) , πΜ (π‘) = πΆπ§
(16) πΜ (π‘0 ) ] = π§0 . π§ (π‘0 ) = [ Μ 0 (π‘0 ) π₯
Μ (β
)) π½ (π‘0 , π₯0 ; π’ (12) 1 β Μ π (π‘) π
Μ π’ (π‘)] ππ‘} , = πΈπ‘0 ,π₯0 { β« [Μππ (π‘) ππΜ πΜ (π‘) + π’ 2 π‘0
],
Μ=[ π
πΜ π₯ (π‘) = [πΜ π₯ (π‘) + πΜ π’ (π‘) + (βΞ) πΜ (π‘ + ππ )] ππ‘
π‘
Μ (π‘) π₯
0 πΆ
Subtracting (7) by (1) from both sides, the following could be obtained:
Μ (π‘) = πΆΜ π¦ π₯ (π‘) ,
πΜ (π‘)
π§ (π‘) = [
Μπ (π‘) = π (π‘) β π (π‘ + ππ ) .
+ π (π‘) ππ΅π‘ ,
(13)
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 1 Μπ½ (π‘0 , π§0 ; π’ Μ (β
)) = πΈπ‘0 ,π§0 { 2 β
Μ π§ (π‘) + π’ Μπ ππΜ πΆ) Μ π (π‘) π
Μ β
β« [π§ (π‘) (πΆ π’ (π‘)] ππ‘} .
(17)
π
π‘0
Μ π(π‘) + ΞΜπ(π‘ Μ + ππ ) in Due to the existence of the term π·Μ (16), the conclusion of the optimal regulation problem in
4
Mathematical Problems in Engineering
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 Μπ½ = lim πΈπ‘0 ,π§0 { 1 π‘π ββ 2 π‘π
(18)
π
Μ ππΜ πΆπ§ Μ (π‘) + π’ Μ (π‘) π
Μ π’ (π‘)] ππ‘} . β
β« [π§ (π‘) πΆ π
π
Μ ΓΜ π is the unique semipositive definite solution of where π β π
π Riccati differential equation
Μ ππ + πΆ Μ = 0. Μ β πππ
Μ β1 π Μπ ππΜ πΆ Μ π π + ππ π
Μ π) Μ is Proof. In terms of the performance index (17), if (π, Μ π) Μ is detectable, according to the stabilizable and (ππΜ1/2 πΆ, known conclusion, there exists a unique semipositive definite solution π of Riccati algebra equation
π‘0
Secondly, the finite-time horizon optimal regulation problem with the performance index ΜΜπ½ = πΈπ‘0 ,π§0 { 1 2 π‘π
π
π‘0
is 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 ππ = ππ(π‘π ) = 0 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), is π
Μ {π (π‘) π§ (π‘) + π (π‘)} , Μ (π‘) = βπ
β1 π π’ 1 (π+π)Γ(π+π)
where π(π‘) β₯ 0 β π
equation
(20)
π
lim π (π‘) = π.
π
As a result, the solving of π(π‘) in (22) is also simplified. Μ π π. Equation (22) and its termΜ ππ
Μ β1 π Denote ππ = πβ inal condition can be expressed as Μ π (π‘) + πΞΜπΜ (π‘ + ππ )) = 0, πΜ (π‘) + πππ π (π‘) + (ππ·Μ π (π‘π ) = 0.
Μ π ) π (π‘) Μ π β π (π‘) ππ
Μ β1 π πΜ (π‘) + (π Μ π (π‘) + π (π‘) ΞΜπΜ (π‘ + ππ )) = 0, + (π (π‘) π·Μ
π‘
Μ (π‘, π‘0 ) = exp [β β« ππ ππ‘] = exp [(π‘0 β π‘) ππ ] , Ξ¦ π π π‘0
Μ (π‘, π‘0 ) , Μ β1 (π‘0 , π‘) = Ξ¦ Ξ¦
π‘0 β€ π β€ π‘. Therefore, the solution of (27) is Μ (π‘, π‘0 ) π (π‘0 ) π (π‘) = Ξ¦ π‘
Μ (π‘, π ) (ππ·Μ Μ π (π‘) + πΞΜπΜ (π‘ + ππ )) ππ , ββ« Ξ¦
(30)
Μ (π‘π , π‘) π (π‘) π (π‘π ) = Ξ¦
π
Μ ππ§ (π‘) Μ (π‘) = βπ
β1 π π’ π
Μ π π) ] Μ β ππ
Μ β1 π exp [(π‘0 β π‘) (π
Μ π (π ) + πΞΜπΜ (π + ππ )) ππ , β
(ππ·Μ
(29)
and the backward form is
Theorem 12. If Assumptions 1 to 5 hold, the optimal controller Μ (π‘) of (16) with the performance index (17) can be expressed π’ as
π‘
Μ (π‘, π‘0 ) = Ξ¦ Μ (π‘, π‘π ) Ξ¦ Μ (π‘π , π‘0 ) , Ξ¦
π‘0
Thirdly, Theorem 12 is received by letting π‘π β β.
π‘+ππ
(28)
with the property
(22)
with the terminal condition π(π‘π ) = 0.
Μ β« β π
β1 π
(27)
Based on linear system theory [19], the fundamental solution matrix of the homogeneous system corresponding to (27) is
(21)
With the terminal condition π(π‘π ) = 0 and π(π‘) the differential equation is obtained:
(26)
π‘ββ
π
Μ Μ ππΜ πΆ. +πΆ
(25)
And when π‘ β β,
satisfies the Riccati differential
Μ ππ (π‘) Μ π (π‘) + π (π‘) π Μ β π (π‘) ππ
Μ β1 π βπΜ (π‘) = π
π
Μ ππ + πΆ Μ = 0. Μ π π + ππ Μ β πππ
Μ β1 π Μπ ππΜ πΆ π
(19)
Μ ππΜ πΆπ§ Μ (π‘) + π’ Μ π (π‘) π
Μ π’ (π‘)] ππ‘} β
β« [π§π (π‘) πΆ
(24)
(23)
π‘π
Μ (π‘π , π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ . ββ« Ξ¦
(31)
π‘
With the terminal condition π(π‘π ) = 0, the following will be obtained: π (π‘) π‘π
Μ (π‘π , π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ . Μ β1 (π‘π , π‘) β« Ξ¦ =Ξ¦ π‘
(32)
Mathematical Problems in Engineering
5
According to the property of the reference signals in Assumption 3, π(π‘) can be expressed as Μ β1 (π‘π , π‘) π (π‘) = Ξ¦ β
β«
π‘π β§(π‘+ππ )
π‘
π‘π β§(π‘+ππ )
π‘
=β«
Corollary 14. If Assumptions 1 to 5 hold, the optimal preview controller π’(π‘) of (1) with the performance index (12) can be expressed as
Μ (π‘π , π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ Ξ¦
π
(33) π
Μ β« β π
β1 π
Μ (π‘π , π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ Ξ¦
π‘π β§(π‘+ππ )
π‘
π
π
π‘π β§(π‘+ππ )
π‘
Μ (π‘, π ) Ξ¦
(34)
where π‘π β§ (π‘ + ππ ) = min{π‘π , π‘ + ππ }. 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 π’ as Μ (π‘) π’ π
π
(35)
π‘+ππ
Μ (π‘, π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ , Ξ¦
π‘
Μ ΓΜ π where π β π
π is the unique semipositive definite solution of Riccati differential equation π
β1 Μ π
π
Μ = 0. Μ β πππ
Μ π π+πΆ Μ ππΜ πΆ Μ π + ππ π
(36)
β1 Μ π
π
π‘+ππ
π‘
β1 Μ π
β
Μ (π‘, π ) ππ·Μ Μ π (π ) ππ + π’β (π‘) Ξ¦ β1 Μ π
With π(π‘)ππ΅π‘ being an external disturbance, the stability and detectability of system (16) are studied only when π(π‘) = 0. Here, (16) becomes the deterministic system Μ (π‘) + πΜ Μ π’ (π‘) + π·Μ Μ π (π‘) + ΞΜπΜ (π‘ + ππ )] ππ‘. ππ§ (π‘) = [ππ§
(40)
As 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:
π π ] Ξ =[ πΆ 0 ] [
(41)
is of full row rank. (3) The matrices pair (π, π) is stabilizable.
π’ (π‘) = βπ
π π1 π (π‘) β π
π π2 π₯ (π‘) Μ β« β π
β1 π
5. Stability of Closed-Loop System
(2) The matrix (37)
Μ Γ π and π Μ Γ π. According to where π1 and π2 are matrices of π π‘+ππ πΜ(π‘) Μ Μ Μ + ππ )ππ = 0, the (8), (35), π§(π‘) = [ π₯Μ(π‘) ], and β«π‘ Ξ¦(π‘, π )πΞπ(π following will be obtained: β1 Μ π
From (39) it can be seen that there are four parts contained in the preview controller of (1): the first part is the inteΜ π π1 π(π‘), stemmed from the gral of tracking error term βπ
β1 π inducing of the integrator. The second part is the state feedΜ π π2 π₯(π‘). The third one βπ
β1 π Μ π β«π‘+ππ Ξ¦(π‘, Μ back βπ
β1 π π‘ Μ π )ππ·Μπ(π )ππ is the reference preview compensation. The last Μ π π2 Ξ]π(π‘ + ππ ) is the previewable complement one [πΎ + π
β1 π of the assistant system.
(1) The matrices ππΜ and π
are both positive definite.
Then, let the matrix π be separated into π = [π1 π2 ] ,
(39)
π
Μ π (π ) + πΞΜπΜ (π + ππ )) ππ , β
(ππ·Μ
Μ ππ§ (π‘) = βπ
β1 π
Μ (π‘, π ) ππ·Μ Μ π (π ) ππ Ξ¦
Μ π2 Ξ] π (π‘ + ππ ) . + [πΎ + π
β1 π
Μ (π‘, π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ . Ξ¦
Μ ππ§ (π‘) β π
β1 π Μ β« Μ (π‘) = βπ
β1 π π’ 1
π‘+ππ
π‘
Substituting (33) into (20) yields
Μ β« β π
β1 π
π
Μ π1 π (π‘) β π
β1 π Μ π2 π₯ (π‘) π’ (π‘) = βπ
β1 π
Μ (π‘, π‘π ) =Ξ¦ β
β«
Due to the fact that π₯β (π‘) = Ξπβ (π‘ + ππ ), π’β (π‘) = πΎπβ (π‘ + ππ ), and πβ (π‘) = 0, the following corollary is obtained.
(4) The matrices pair (πΆ, π) is detectable. (38)
β
+ π
π π2 π (π‘) + π
π π2 π₯ (π‘) , Μ ΓΜ π where π β₯ 0 β π
π satisfies Riccati differential equation (24).
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.
6
Mathematical Problems in Engineering
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: ππ§ (π‘)
Theorem 17. If the conditions in Theorem 16 hold, the closedloop system can be received: ππ₯ (π‘) = [ππ₯ (π‘) + ππ’ (π‘)] ππ‘ + π (π‘) ππ΅π‘ , πΜ π₯ (π‘) = [πΜ π₯ (π‘) + ππ’ (π‘) + πΏ (π¦ (π‘) β πΆΜ π₯ (π‘))] ππ‘
Μ (π‘) + πΜ Μ π’ (π‘) + π·Μ Μ π (π‘) + ΞΜπΜ (π‘ + ππ )] ππ‘ = [ππ§
+ π (π‘) ππ΅π‘ , π¦ (π‘) = πΆπ₯ (π‘) ,
+ π (π‘) ππ΅π‘ , Μ (π‘) π’
π‘
π
Μ π1 β« (πΆΜ π₯ (π ) β π (π )) ππ π’ (π‘) = βπ
β1 π π‘0
β1 Μ π
= βπ
π ππ§ (π‘) π
Μ β« β π
β1 π
π‘+ππ
π‘
(42)
(48)
π
Μ π2 π₯ Μ (π‘) β π
β1 π
Μ (π‘, π ) (ππ·Μ Μ π (π ) + πΞΜπΜ (π + ππ )) ππ , Ξ¦
π
Μ β« β π
β1 π
Μ (π‘) πΜ (π‘) = πΆπ§
π‘+ππ
π‘
Μ (π‘, π ) ππ·Μ Μ π (π ) ππ Ξ¦
π
Μ π2 Ξ] π (π‘ + ππ ) , + [πΎ + π
β1 π
πΜ (π‘0 )
] = π§0 ; π§ (π‘0 ) = [ Μ 0 (π‘0 ) π₯
where (π β πΏπΆ) is stable and (48) achieve the complete regulation.
namely, lim πΜ (π‘) = 0.
π‘ββ
7. Numerical Simulation (43)
Furthermore, the following will be obtained: lim π (π‘) = 0.
π‘ββ
Example 1. Consider the stochastic control system π[
(44)
π₯1 (π‘) π₯2 (π‘)
] = [[
0
1
β1 β1
π₯1 (π‘)
][
0 ] + [ ] π’ (π‘)] ππ‘ 1 π₯2 (π‘)
+ 0.1 ππ΅π‘ , π₯1 (π‘)
6. State Observer
π¦ (π‘) = [1 0] [
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 πΜ π₯ (π‘) = [πΜ π₯ (π‘) + ππ’ (π‘) + πΏ (π¦ (π‘) β πΆΜ π₯ (π‘))] ππ‘ + π (π‘) ππ΅π‘ ,
π₯ (π‘0 ) = π₯0 , π‘ β [π‘0 , β)
(45)
π‘ββ
where the coefficient matrices are 0
1
], β1 β1
0 π = [ ], 1
(50)
πΆ = [1 0] , (46)
Μ (π‘). where π₯(π‘) = π₯(π‘) β π₯ So, if (π β πΏπΆ) is stable, the following will hold: lim π₯ (π‘) = 0,
],
π=[
could be designed. Subtracting (45) from (1) on both sides yields ππ₯ (π‘) = [(π β πΏπΆ) π₯ (π‘)] ππ‘,
π₯2 (π‘)
(49)
π (π‘) = [
0.1 0.1
],
respectively. Therefore, the coefficient matrices in (16) are (47)
Μ (π‘) in the observer which means that the state vector π₯ 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.
0 1 [ Μ = [0 0 π
0
1] ],
[0 β1 β1] 0 [0] Μ π = [ ], [1]
7
1.2
1.2
1
1
0.8
0.8
0.6
0.6 y&r
y&r
Mathematical Problems in Engineering
0.4
0.4
0.2
0.2
0
0
β0.2
0
5
10
15
20 t (s)
25
30
35
β0.2
40
0
Reference signal: r(t) lr = 0 s, output: y(t)
5
10
15
20 t (s)
25
30
35
40
Reference signal: r(t) lr = 0.5 s, output: y(t)
Figure 1: Closed-loop output response with nonpreview and reference signal π(π‘).
Figure 2: Closed-loop output response with ππ = 0.5 s preview and reference signal.
β1 0 [ 0 β1] Μ π·=[ ] [0
1.2 1
0]
0.8
(51) and Brownian motion π΅π‘ satisfies π΅π‘ βΌ π (0, 0.01) .
(52)
y&r
0.6 0.4
Let the initial state of (1) be 0.2
0 π₯ (0) = [ ] , 0
(53)
where π₯1 and π₯2 are mutually independent. The preview lengths of the reference signals are πππ = 0, 0.5, 0.75 s (π = 1, 2, 3), respectively. The weight matrices of the performance index are ππΜ = 1, (54) π
= 1, and the solution of Riccati equation Μ π π + ππ Μ ππ + πΆ Μ=0 Μ β πππ
Μ β1 π Μπ ππΜ πΆ π
0 β0.2
0
5
10
15
20 t (s)
25
30
35
40
Reference signal: r(t) lr = 0.75 s, output: y(t)
Figure 3: Closed-loop output response with ππ = 0.75 s preview and reference signal π(π‘).
(55) When the reference signal is
is 2.1533 1.8184 1.000 [1.8184 1.9157 1.153] π=[ ], [ 1.000
(56)
1.153 0.818]
which is a positive definite matrix. The coefficient matrix ππ in (25) is 0
Μ ππ = [ Μ β ππ
Μ β1 π ππ = π [0
1
0
0
1
] ].
[β1 β2.153 β1.818]
(57)
0, { { { { π (π‘) = {0.2 (π‘ β 5) , { { { {1,
0 β€ π‘ β€ 5, 5 < π‘ β€ 10,
(58)
π‘ > 10,
Figures 1β3 can be obtained. Figure 1 shows the tracking controller with nonpreview. Figure 2 shows the tracking controller with ππ = 0.5 s preview. Figure 3 shows the tracking controller with ππ = 0.75 s preview. Comparing the above three figures, it can be seen that
8
Mathematical Problems in Engineering 0.2
0.15
0.15
0.1
0.1
0.05 0 y&r
y&r
0.05 0 β0.05
β0.1
β0.1
β0.15
β0.15 β0.2
0
5
10
15
20 t (s)
25
30
35
40
Reference signal: r(t) lr = 0 s, output: y(t)
0.2
0
5
10
15
20 t (s)
25
30
35
40
Figure 6: Closed-loop output response with ππ = 0.75 s preview and reference signal π(π‘).
Figure 6 shows the tracking controller with ππ = 0.75 s 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.
0.15 0.1
8. Conclusion
0.05 y&r
β0.2
Reference signal: r(t) lr = 0.75 s, output: y(t)
Figure 4: Closed-loop output response with nonpreview and reference signal π(π‘).
0 β0.05 β0.1 β0.15
β0.05
0
5
10
15
20
25
30
35
40
t (s) Reference signal: r(t) lr = 0.5 s, output: y(t)
This paper has studied the optimal preview control problem for a class of continuous stochastic system in the infinite horizon. By introducing the integrator and the assistant system, the stochastic augmented error system is constructed. Compared with the finite horizon, the performance index is simplified. The stability of the stochastic augmented error system is also studied and the observer for the original stochastic system is designed. Finally, the simulation example shows the effectiveness of the conclusion in this paper.
Competing Interests
Figure 5: Closed-loop output response with ππ = 0.5 s preview and reference signal.
The authors declare that they have no competing interests.
Acknowledgments based on the preview controller, the output signals can track the reference signals much faster and with less tracking error. When the reference signal is 0, 0 < π‘ β€ 5, { { { { π (π‘) = {0.1 sin 0.5 (π‘ β 5) , 5 < π‘ β€ 5 + 5π, { { { π‘ > 5 + 5π, {0,
(59)
Figures 4β6 can be obtained. Figure 4 shows the tracking controller with nonpreview. Figure 5 shows the tracking controller with ππ = 0.5 s preview.
This work was supported by the National Natural Science Foundation of China (Grant no. 61174209).
References [1] M. Tomizuka, βOptimal continuous finite preview problem,β IEEE Transactions on Automatic Control, vol. 20, no. 3, pp. 362β 365, 1975. [2] Y. Xu, F. Liao, L. Cui, and J. Wu, βPreview control for a class of linear continuous time-varying systems,β International Journal of Wavelets, Multiresolution and Information Processing, vol. 11, no. 2, Article ID 1350018, 1350018, 14 pages, 2013. [3] F. Liao, M. Cao, Z. Hu, and P. An, βDesign of an optimal preview controller for linear discrete time causal descriptor systems,β
Mathematical Problems in Engineering
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]
[14] [15]
[16]
[17]
[18] [19]
International Journal of Control, vol. 85, no. 10, pp. 1616β1624, 2012. F. Liao, K. Takaba, T. Katayama, and J. Katsuura, βDesign of an optimal preview servomechanism for discrete-time systems in a multirate setting,β Dynamics of Continuous, Discrete & Impulsive Systems. Series B. Applications & Algorithms, vol. 10, no. 5, pp. 727β744, 2003. T. Katayama, T. Ohki, T. Inoue, and T. Kato, βDesign of an optimal controller for a discrete-time system subject to previewable demand,β International Journal of Control, vol. 41, no. 3, pp. 677β699, 1985. F. Liao, Y. Guo, and Y. Y. Tang, βDesign of an optimal preview controller for linear time-varying discrete systems in a multirate setting,β International Journal of Wavelets, Multiresolution and Information Processing, vol. 13, no. 6, Article ID 1550050, 1550050, 19 pages, 2015. T. Tsuchiya and T. Egami, Digital Preview and Predictive Control, Translated by Fucheng Liao, Beijing Science and Technology Press, Beijing, China, 1994. M. Cao and F. Liao, βDesign of an optimal preview controller for linear discrete-time descriptor systems with state delay,β International Journal of Systems Science, vol. 46, no. 5, pp. 932β943, 2015. Z. Y. Zhen, Z. S. Wang, and D. B. Wang, βInformation fusion estimation based preview control for a discrete linear system,β Acta Automatica Sinica, vol. 36, no. 2, pp. 347β352, 2010. T. Katayama and T. Hirono, βDesign of an optimal servomechanism with preview action and its dual problem,β International Journal of Control, vol. 45, no. 2, pp. 407β420, 1987. F. Liao, Y. Y. Tang, H. Liu, and Y. Wang, βDesign of an optimal preview controller for continuous-time systems,β International Journal of Wavelets, Multiresolution and Information Processing, vol. 9, no. 4, pp. 655β673, 2011. F. Liao, Z. Ren, M. Tomizuka, and J. Wu, βPreview control for impulse-free continuous-time descriptor systems,β International Journal of Control, vol. 88, no. 6, pp. 1142β1149, 2015. F. Liao and H. Xu, βApplication of the preview control method to the optimal tracking control problem for continuous-time systems with time-delay,β Mathematical Problems in Engineering, vol. 2015, Article ID 423580, 8 pages, 2015. Λ om, Introduction to Stochastic Control Theory, Translated A. AstrΒ¨ by Yuhuan Pan, Science Press, Beijing, China, 1983. J. Wu, F. Liao, and M. Tomizuka, βOptimal preview control for a linear continuous-time stochastic control system in finite-time horizon,β International Journal of Systems Science, vol. 48, no. 1, pp. 129β137, 2017. B. Oksendal, Stochastic Differential Equations: An Introduction with Applications, Springer, New York, NY, USA, 6th edition, 2005. T. H. Cormen, C. E. Leiserson, R. Rivest, and C. Stein, Introduction to Algorithms, MIT Press and McGraw-Hill, Boston, Mass, USA, 2nd edition, 2001. J. Yong and X. Y. Zhou, Stochastic Control: Hamiltonian Systems and HJB Equations, vol. 43, Springer, New York, NY, USA, 1999. D. Zheng, Linear System Theory, Tsinghua University Press, Beijing, China, 2002.
9
Advances in
Operations Research Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Advances in
Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Applied Mathematics
Algebra
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Probability and Statistics Volume 2014
The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Differential Equations Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Submit your manuscripts at http://www.hindawi.com International Journal of
Advances in
Combinatorics Hindawi Publishing Corporation http://www.hindawi.com
Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of Mathematics and Mathematical Sciences
Mathematical Problems in Engineering
Journal of
Mathematics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Discrete Mathematics
Journal of
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Discrete Dynamics in Nature and Society
Journal of
Function Spaces Hindawi Publishing Corporation http://www.hindawi.com
Abstract and Applied Analysis
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Journal of
Stochastic Analysis
Optimization
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014