On Irregular Linear Quadratic Control: Deterministic Case

Size: px
Start display at page:

Download "On Irregular Linear Quadratic Control: Deterministic Case"

Transcription

1 1 On Irregular Linear Quadratic Control: Deterministic Case Huanshui Zhang and Juanjuan Xu arxiv: v2 math.oc 1 Mar 218 Abstract This paper is concerned with fundamental optimal linear quadratic (LQ control with irregular Riccati equation where the controller could not be solved from the equilibrium condition like the classical LQ control. The problem remains challenging because it is not clear what tool should be applied instead to design the controller although the problem has been investigated since 7 s last century Different from the standard regular Riccati equation based LQ-control theory, a new approach of two-layer optimization is proposed. With the approach, we show the controller entries of irregular-lq controller are derived from two equilibrium conditions in two different layers, which is different from the classical regular LQ control where all the controller entries were obtained from equilibrium condition in one layer based on regular Riccati equation. It is shown that the irregular LQ control is totally different from the regular one due to the irregular controller must guarantee the state terminal constraint ofp 1(Tx(T. The presented results also clarify the differences of open-loop control from closed-loop control in the aspects of solvability and controller design. Examples are presented to show the effectiveness of the proposed approach. The presented results contain fundamental contribution to the LQ-control theory. The stochastic case will be reported in another paper. Index Terms Irregular LQ control; Riccati equation; Openloop solvability; Closed-loop solvability. I. INTRODUCTION Linear quadratic (LQ control is one of the most fundamental and widely used tools in modern engineering. Two approaches including maximum principle 6 and dynamic programming 7 have been extensively studied since 195 s. See 8, 9 and references therein. In the literature it is typically assumed that the cost function has a positive definite weighting matrix for the control term, and a positive semidefinite weighting matrix for the state term 14. In this case, the Riccati equation is used to characterize the solvability of the underlying control problem considered a special case of singular LQ control problem where the weighting matrix of control is zero. In 21, the LQ control with semi-positive weighting matrix of control term was studied. It is noted that only the sub-optimal controller is given although it shows that a constructed sequence of cost function converges to the optimal one. 12, 16 found the well-posedness of the stochastic LQ problem where the cost weighting matrices for the state and the control are allowed to be indefinite. A new type of Riccati equation, called the general Riccati equation, is proposed to explicitly construct the optimal feedback control. In very recently, 19 studied irregular LQ problem with indefinite weighting matrices, focusing on the closed-loop and open-loop solvability for the stochastic control. 19 showed that the solvability obtained in previous work 16 is only applicable to the case of closed-loop control, and the openloop solvability given in 16 is incorrect by presenting a counterexample. In particular, the open-loop control was obtained by constructing a minimizing sequence and the closed-loop solution was given in the feedback form of the state where the feedback gain is required to be in the space L 2 (t,t;r m n. In this paper, following-up the previous works 16, 19, we will study the irregular LQ control problem and and try to clarify the following questions: Question 1: The solvability conditions for the irregular- LQ problem? Question 2: What is the difference of irregular-lq from the regular-lq? Question 3: How to design the irregular controller in open-loop and closed-loop form? By proposing so-called two-layer optimization approach, Question 1 and 2 are to be answered in Section IV, while Question 3 will be answered in Section V. The remainder of the paper is organized as follows. Section II presents the studied problem and some preliminaries. The optimal solution is shown in Section III when Riccati equation is regular. The results for irregular LQ control problem are given in Section IV and V respectively. Two numerical examples are given in Section VI. Some concluding remarks are given in Section VII. The following notations will be used throughout this paper: R n denotes the family of n dimensional vectors. x means the transpose of x. A symmetric matrix M > ( means strictly positive definite (positive semi-definite. Range(M represents the range of the matrix M. M means the Moore- Penrose inverse of the matrix M. We also introduce the following set: This work is supported by the National Natural Science Foundation of China ( , This version has been submitted to IEEE Trans. on Automatic Control for publication H. Zhang and J. Xu are with School of Control Science and Engineering, Shandong University, Jinan, Shandong, P.R.China hszhang@sdu.edu.cn, jnxujuanjuan@163.com L 2 (t,t;r m n {K(t R m n, K(t 2 dt < }, t Ut,T {u(t R m, t u(t 2 dt < }.

2 2 II. DEFINITIONS AND PRELIMINARIES Consider the continuous-time system: ẋ(t A(tx(t+B(tu(t, x(t x, (1 where x R n is the state, u R m is the control input. The matrices A(t, B(t are time-varying matrices with appropriate dimension. x is the given initial value. The cost function is given by J(t,x ;u t x (tq(tx(t+u (tr(tu(tdt +x (THx(T. (2 where Q(t, R(t, H are symmetric matrices with appropriate dimension. Problem (LQ: For any given initial pair (t,x, find a u such that J (t,x ;u min u( J(t,x ;u. V(t,x. To guarantee the convexity of Problem (LQ, it is assumed that R(t. We now introduce the definitions of solvability. Definition 1. Problem (LQ is said to be finite at initial pair (t,x,t R n, if V(t,x >. (3 Problem is said to be finite at t,t if (3 holds for all x R n. Problem is said to be finite if (3 holds for all (t,x,t R n. Definition 2. An element u is called an open-loop optimal control of Problem (LQ for the initial pair (t,x,t R n, if J (t,x ;u J(t,x ;u, u(. (4 If an open-loop optimal control (uniquely exists for (t,x,t R n, Problem (LQ is said to be (uniquely open-loop solvable at (t,x,t R n. Problem (LQ is said to be (uniquely open-loop solvable at t,t if for the given t, (4 holds for all x R n. Problem (LQ is said to be (uniquely open-loop solvable if it is open-loop solvable at all (t,x,t R n. Definition 3. An element u ( K( x ( Ut,T is called a closed-loop optimal control of Problem (LQ on t,t if J (t,x ;K( x ( J(t,x ;u, x R n, u( Ut,T. (5 wherex ( is the solution to the following closed-loop system: ẋ (t A(t+B(tK(tx (t, x(t x. (6 If a closed-loop optimal controller exists on t,t, Problem (LQ is said to be (uniquely closed-loop solvable on t,t. Problem (LQ is said to be (uniquely closed-loop solvable if it is (uniquely closed-loop solvable on any t,t. Remark 1. Noting that the closed-loop solvability is only assumed that K( x ( Ut,T which allows the case that K( L 2 (t,t;r m n. We firstly recall the Moore-Penrose inverse of a matrix which will be used this paper repeatedly. From 23, for a given matrix M R n m, there exists a unique matrix in R m n denoted by M such that MM M M, M MM M, (MM MM, (M M M M. The matrix M is called the Moore-Penrose inverse of M. The following lemma is from 22. Lemma 1. Let matrices L, M and N be given with appropriate size. Then, LXM N has a solution X if and only if LL NMM N. Moreover, the solution of LXM N can be expressed as X L NM +Y L LYMM, where Y is a matrix with appropriate size. In particular, let M I, we have LX N has a solution if and only if LL N N. This is also equivalent to Range(N Range(L where Range(N is the range of N. We now present the the maximum principle. Lemma 2. Problem (LQ is solvable if and only if the following forward and backward differential equations (FBDEs admit a solution: ẋ(t A(tx(t + B(tu(t, (7 ṗ(t A (tp(t+q(tx(t, (8 R(tu(t+B (tp(t. (9 with x(t x and p(t Hx(T. Proof. The result follows from the maximum principle 11. So we omit it. From Lemma 2, we note that the key to solve Problem (LQ is to obtain the solution to (7-(9. To this end, without loss of generality, we assume that p(t P(tx(t+Θ(t (1 where P, Θ are to be determined with terminal values P(T H and Θ(T. Taking derivative to (1 yields ṗ(t P(tx(t+P(t A(tx(t+B(tu(t + Θ(t. (11 From (8 and (1, it is obtained that ṗ(t A (tp(tx(t+a (tθ(t+q(tx(t. (12 With a comparison of (11 and (12, it follows that P(tx(t+P(tA(tx(t +P(tB(tu(t + Θ(t+A (tp(tx(t+a (tθ(t +Q(tx(t. (13 Using (1, the equilibrium condition (9 is reformulated as R(tu(t+B (tp(t Υ (tu(t+γ (tx(t+b (tθ(t, (14

3 3 where Υ (t R(t and Γ (t B (tp(t. There exists two cases: Range Γ (t Range Υ (t and Range Γ (t Range Υ (t. III. OPTIMAL CONTROL WITH REGULAR RICCATI EQUATION In this section, we consider the case: Range Γ (t Range Υ (t, (15 In this case, the following results are obtained directly from 16. Lemma 3. If (15 holds in (14 where P(t is given by the following Riccati equation: P(t+A (tp(t+p(ta(t +Q(t Γ (tυ (tγ (t, (16 with P(T H, then Problem (LQ is solvable and the optimal solution is given by u(t Υ (tγ (tx(t+ I Υ (tυ (t z(t, (17 where R (t represents the Moore-Penrose inverse ofr(t and z(t is an arbitrary vector with compatible dimension. In this case, Θ(t holds in (1. In the case of Range Γ (t Range Υ (t, the Riccati equation (16 is called regular. IV. SOLUTION TO IRREGULAR LQ CONTROL This section focuses on the LQ control problem when (15 does not hold, i.e., Range Γ (t Range Υ (t. (18 In this case the problem is termed as the irregular LQ control problem (IR-LQ, and the controller u(t can not be solved from the equilibrium condition (14. From (18, it is clear that Υ (t is not invertible. ( We assume that rank(υ (t m (t < m. Thus rank I Υ (tυ (t m m (t >. It is not difficult to know that there is an elementary row transformation matrix T (t such that T (t I Υ (tυ (t Υ T (t, (19 where Υ T (t R m m(t m is full row rank. Further denote C (t Γ (t I Υ (tυ (t T 1 (t, B (t B(t( I Υ (tυ (t T 1 (t, (2 and A (t A(t B(tΥ (tγ (t, D (t B(tΥ (tb (t. Lemma 4. Under the condition (18, Problem (IR-LQ is solvable if and only if there exists u 1 (t R m m(t such that C (tx(t+b (tθ(t,θ(t (21 holds, where u 1 (t, x(t,θ(t satisfy the FBSDEs ẋ(t A (tx(t+d (tθ(t+b (tu 1 (t, (22 Θ(t A (tθ(t+c (tu 1(t. (23 Proof. Necessity. If the problem is solvable, then from Lemma 2, (14 must hold. Thus, (14 can be equivalently written as u(t Υ (Γ (t (tx(t+b (tθ(t + I Υ (tυ (t z(t, (24 where z(t is an arbitrary vector with compatible dimension such that following equality hold Denote (I Υ (tυ (t Γ (tx(t+b (tθ(t.(25 T (t I Υ (tυ (t z(t u 1 (t, (26 where u 1 (t Υ T (tz(t R m m(t. Now we rewrite (25 as (21. First, note that I Υ (tυ (I (t Υ (tυ (I (t Υ (tυ (t (I Υ (tυ T (t (t T 1 (t (I Υ (tυ (t Υ T (t ( T 1 (t (I Υ (tυ (t, (27 where (19 has been used in the derivation of the last equality. By using (2, (25 can be written as Υ T (t C (tx(t+b (tθ(t. (28 Note that Υ T (t is full column rank, (28 is rewritten as (21 directly. By substituting (24 into (13 and using (16, it yields that P(tx(t+P(tA(tx(t+A (tp(tx(t +A (tθ(t+q(tx(t Γ (tυ (Γ (t (tx(t +B (tθ(t +Γ (t I Υ (tυ (t z(t Θ(t+ ( A (t Γ (tυ (tb (t Θ(t +Γ (t I Υ (tυ (t z(t. (29

4 4 2 In view of the fact that I Υ (tυ (t I Υ (tυ (t, it is obtained that Γ (t I Υ (tυ (t z(t Γ (I (t Υ (tυ (t T 1 (tt (t I Υ (tυ (t z(t Γ (t I Υ (tυ (t T 1 (t u 1 (t C (t u 1 (t C (tu 1(t. (3 Thus, from (29 we have Θ(t A (tθ(t+c (tu 1(t, (31 this implies that the dynamic of Θ is as (23. By substituting ( (24 into (7 and combining with the fact that B(t I Υ (tυ (t z(t B (tu 1 (t which can be obtained similar to (3, one has the dynamic (22 of the state. Sufficiency We now show Problem (IR-LQ is solvable if there exists u 1 (t to achieve (21. In fact, if (21 is true then (24 and (25 can be rewritten as (14. Further, by taking reverse procedures to (1-(14, it is easily verified that p(t P(tx(t+Θ(t, where x, Θ satisfy (21-(23, are the solutions to (7-(9. Thus, Problem (IR-LQ is solvable according to Lemma 2. The proof is now completed. Remark 2. Applying (26 we have I Υ (tυ (t z(t T 1 (t u 1 (t G (tu 1 (t, (32 where G (t T 1 (t. Thus (24 can be rewritten as u(t Υ (Γ (t (tx(t+b (tθ(t +G (tu 1 (t. (33 Observation 1. The solution to FBSDEs (7-(9 is homogeneous, i.e., Θ(t, if and only if the Riccati equation (16 is regular. Proof. Necessity Assume that the Riccati equation (16 is not regular, then (I Υ (tυ (t Γ (t, (34 which implies C (t using (3. Thus according to Lemma 4, Θ(t satisfies (23 which is not equal to zero. This is a contradiction. Sufficiency The necessity follows from Theorem 3. The proof is now completed. To solve forward and backward differential equations (FB- DEs (21-(23, Denote that Θ 1 (t Θ(t P 1 (tx(t, (35 where P 1 (t obeys a Riccati-like equation associated with (21-(23 as P 1 (t+p 1 (ta (t+a (tp 1 (t+p 1 (td (tp 1 (t while Γ 1(tΥ 1 (tγ 1(t, P 1 (T, (36 Υ 1 (t, Γ 1 (t C (t+b (tp 1(t, From (35 it is clear that Θ 1 (T P 1 (Tx(T. Remark 3. It is interesting to know that P(t+P 1 (t obeys the same Riccati equation as P(t P(t+ P 1 (t+ P(t+P 1 (t A(t+A (t P(t B(tR +P 1 (t P(t+P 1 (t (tb (t P(t +P 1 (t +Q(t, with terminal value P(T+P 1 (T. Theorem 1. Problem (Deterministic IR-LQ is solvable if and only if there exists P 1 (t of (36 with terminal value P 1 (T such that Γ 1 (t, i.e., C (t+b (tp 1(t, t T, (37 and there exists u 1 (t to achieve Θ(T P 1 (Tx(T, (38 where x(t obeys ẋ(t A (t+d (tp 1 (t x(t+b (tu 1 (t, (39 with initial value x(t x. In this case, the controller u(t is given by (33. Proof. Sufficiency By taking derivative to P 1 (tx(t, it yields that dp 1 (tx(t dt P 1 (tx(t+p 1 (ta (t+d (tp 1 (t x(t+p 1 (tb (tu 1 (t A (tp 1 (tx(t Γ 1(tΥ 1 (tγ 1(tx(t +P 1 (tb (tu 1 (t A (tp 1(tx(t C (tu 1(t, (4 where (36 with Υ 1 (t and (37 has been used in the derivation of the last equality. On the other hand, using again (37, we have C (tx(t+b (tp 1 (tx(t. (41 Making comparison between (21, (22, (23 and (39, (4, (41, it is obtained that (21-(23 is solvable with Θ(t P 1 (tx(t if P 1 (Tx(T. Following Lemma 4, Problem (LQ is solvable. Necessity The proof will be completed by applying Reduction to Absurdity. Now we assume that (37 does not hold

5 5 for any terminal value P 1 (T. With similar lines to Θ(t, we have ( Θ 1 (t A (t+p 1(tD (t Θ 1 (t+ C (t +P 1 (tb (t u 1 (t, (42 By using (35, (21 and (22 are respectively rewritten as C (t+b 1P 1 (t x(t+b Θ 1 (t, (43 ẋ(t A (t+d (tp 1 (t x(t+d (tθ 1 (t +B (tu 1 (t, (44 Since (37 does not hold, i.e., C (t+b 1 P 1(t, it is easy to know that Θ 1 (t from (43. So, the control u 1 (t can not be determined at this stage. In order to solve u 1 (t, it is necessary to solve FBDEs (42-(44, we further denote that Θ 2 (t Θ 1 (t P 2 (tx(t, (45 where P 2 (t, similar to P 1 (t, obeys a Riccati-like equation associated with (42-(44 as P ( 2 (t+p 2 (t A (t+d (tp 1 (t + A (t +P 1 (td (t P 2 (t+p 2 (td (tp 2 (t, P 2 (T. (46 Now we will show that Θ 2 (t. In fact, by combining (36 and (46 we know that P 1 (t + P 2 (t satisfies the following Riccati-like equation. P 1 (t+ P 2 (t+a (tp 1 (t+p 2 (t+p 1 (t +P 2 (ta (t+p 1 (t+p 2 (td (tp 1 (t+p 2 (t. (47 It is noted that (47 is the same as (36. On the other hand, substituting (45 into (43 yields C (t+b 1 P 1 (t+p 2 (t x(t+b Θ 2(t, (48 Recall the assumption that (37 does not hold for any terminal value P 1 (T( and (47 is the same as (36, we know that C (t + B 1 P 1 (t + P 2 (t for any P 1 (T + P 2 (T. Thus Θ 2 (t from (48. Similarly, we can prove Θ i+1 (t Θ i (t P i+1 (tx(t for i 2,3, if (37 does not hold for any terminal value P 1 (T. This implies that FBDEs (21-(23 are unsolvable, and thus the Problem (Deterministic IR-LQ is unsolvable which leads to contradiction. In other words, if the problem is solvable, there will exist P 1 (t of (36 with terminal value P 1 (T such that (37 is satisfied, i.e., Θ 1 (t for t T. Now we show that there exists u 1 (t to achieve (38. Otherwise, Θ(T P 1 (Tx(T+Θ 1 (T P 1 (Tx(T, thus the problem is unsolvable from Lemma 4. The proof is now completed. Remark 4. Now the Question 1 in introduction has been answered in Theorem 1. Also, the Question 2 is answered at the same time, i.e., a irregular controller must guarantee the state terminal constraint of P 1 (Tx(T which is the essential difference from the regular LQ. Remark 5. Using Theorem 1, it is easy to know that the terminal value P 1 (T of (36 should be chosen such that C (T+B (TP 1(T. V. OPEN-LOOP AND CLOSED-LOOP SOLUTION Theorem 2. Problem (Deterministic IR-LQ is open-loop solvable if and only if there exists P 1 (t of (36 such that Γ 1 (t, and Range P 1 (t Range G 1 t,t, (49 where the Gramian matrix G 1 t,t is defined by G 1 t,t while P 2 (t,s satisfies t P 2 (t,sc (sc (sp 2 (t,sds, (5 P 2 (t,s A (tp 2 (t,s P 2 (t,t I. In this case, the open-loop solution can be given by (33 while u 1 (t is given by u 1 (t C (tp 2 (t,tg 1 t,tp 1 (t x. (51 Proof. Sufficiency : If (49 holds, i.e., Range P 1 (t Range G 1 t,t, then from the proof of Theorem 1 we have that Θ(t P 1 (tx(t ( t T. Also, if (49 holds, then for any x, there exists ζ such that where P 1 (t x G 1 t,tζ, Then it is easy to verify that ζ G 1 t,tp 1 (t x. u 1 (t C (tp 2(t,tζ C (tp 2(t,tG 1 t,tp 1 (t x, is the open-loop controller to achieve Θ(T P 1 (Tx(T. In fact, from (23, we have Θ(t t P 2 (t,sc (su 1 (sds+c, (52 where C is arbitrary constant. It is clear that Θ(T C from (52. On the other hand, by substituting u 1 (s of (51 into the above equation, we have Θ(t P 2 (t,sc (sc (sp 2(t,sds ζ +C t P 1 (t x +C, it follows that C from the fact Θ(t P 1 (t x(t, which implies that Θ(T. Thus the proof of sufficiency is now completed based on Theorem 1. Necessity If the control problem is solvable, it follows from Theorem 1 that there exists P 1 (t such that Γ 1 (t.

6 6 We now prove that (49 holds. Otherwise, Range P 1 (t Range G 1 t,t. Then there exists a nonzero vector ρ such that It is further obtained that ρ P 1 (t ρ, ρ G 1 t,tρ. ρ G 1 t,tρ ρ P 2 (t,sc (sc (sp 2(t,sds ρ t This implies that Thus, we further have ρ t C (sp 2(t,sρ 2 ds. C (sp 2(t,sρ. t P 2 (t,sc (su 1(sds. Let x ρ. From Θ(t P 1 (t x, one has Θ(t P 1 (t ρ. Combining with (52, it gives that ρ P 1 (t ρ ρ P 2 (t,sc (su 1 (sds ρ. t This is a contradiction. Thus (49 holds. The proof is now completed. We further consider the closed-loop solution. To this end, denote P 1 P 1 (t 11 (t P12 1 (t P12 1 (t P 1, 22 (t Assume that P 1 (t is singular, then there exists T 1 (t such that ˆP(t T 1(tP 1 (tt 1 (t, where ˆP(t is invertible. Let and T 1(tT 1 (t T 1(tB (tt 1 (t T1 (t T 2 (t B1 (t B 2 (t T 1(t Â1 (t A (t+d (tp 1 (t T 1 (t  2 (t. (53 Theorem 3. Assume that Γ 1 (t holds and there exists K(t such that T 1 (t+â1(t+b 1 (tt 1 (tk(tt 1(t I t T, (54 then Problem (Deterministic IR-LQ is closed-loop solvable. In this case, the closed-loop solution is given by (4 where where K(t satisfies that (54. u 1 (t K(tx(t, Proof. Denote y(t T 1(tx(t y1 (t y 2 (t, (55 then using (22 and the feedback control z(t K(tx(t, we have ẏ(t T 1 (tx(t+t 1 (A (t (tx(t+d (tθ(t +B (tk(tx(t ( T 1(tT 1 (t+t 1(t A (t+d (tp 1 (t +B (tk(t T 1 (t 1 (tx(t ( T1 (t Â1 (t B1 (t + + T 2 (t  2 (t B 2 (t K(tT 1 (t y(t T1 (t+â1(t+b 1 (tt 1 (tk(tt 1(t T 2 (t+â2(t+b 2 (tt 1(tK(tT 1 (t T By applying (54, we obtain that I ẏ1 (t t T ẏ 2 (t This yields that ẏ 1 (t I t T y 1(t. y1 (t y 2 (t By solving the above equation, it holds that This further implies that Since y 1 (t T t T t y 1 (t,. T 1 (t y(t. y 1 (T. (56 T 1 (TP 1(TT 1 (TT 1 (Tx(T ˆP(T y(t ˆP(Ty 1 (T, combining with the invertibility of ˆP(t, we have thaty 1 (T. On the other hand, note (55, y(t T 1 (Tx(T, then y 2 (T P 1 (Tx(T P 1 (TT 1 (T y 2 (T T 1 (TT 1 (TP 1(TT 1 (T ˆP(T T 1 (T The proof is now completed. y 2 (T y 2 (T

7 7 Remark 6. When P 1 (t is nonsingular, we have the following results. Problem (IR-LQ is closed-loop solvable if and only if Γ 1 (t holds and there exists K(t such that A (t+d (tp 1 (t+b (tk(t 1 I. (57 t T In this case, the closed-loop solution is given by (4 where where K(t satisfies that (57. u 1 (t K(tx(t, VI. NUMERICAL EXAMPLES Example 1: Consider the system and the cost function J(t,x ;u ẋ(t 1 1 u(t, x(t x, t ( 1 u (t u(t dt+x (Tx(T. The solution to Riccati equation (16 with P(T 1 is 1 P(t 1 t T 1. Noting that Υ (t and 1 Γ (t P(t, the LQ problem is irregular. Following 1 Theorem 1, the solvability of the irregular problem is equivalent to that (37 and (38 hold. By solving equation (36 with P 1 (T 1, it yields that P 1 (t 1 t T 1. This further implies that Γ 1 (t C (t+b (tp 1 (t P(t+P 1 (t, that is (37 holds. We then apply Theorem 2 and 3 to find the open-loop and closed-loop controller to achieve (38. Firstly, following Theorem 2 and using (4 and (51, the open-loop solution is u(t where u 1 (t u 1 (t C (tp 2 (t,tg 1 t,tp 1 (t x x(t T t. Based on Theorem 3, the closed-loop solution is u(t where u u 1 (t 1 (t x(t t T. Example 2: Consider the system and the cost function J(t,x ;u ẋ(t 1 1 u(t, x(t x, t ( 1 x (tx(t+u (t +2x (Tx(T. The Riccati equation (16 is in fact u(t dt P(t P 2 (t+1, (58 with terminal value P(T 2. The solution is P(t 3+e 2(t T 3 e 2(t T. Thus A (t P(t,B (t 1,D (t 1,C (t P(t. The Riccati equation of P 1 (t is as P 1 (t 2P 1 (tp(t P1 2 (t. (59 Based on Theorem 1, Problem (IR-LQ is solvable only if there exists an terminal value P 1 (T such P(t+P 1 (t. However, this is impossible. In fact, by substituting P(t P 1 (t into (59, it yields that P 1 (t+p 2 1 (t. Together with (58, it is obtained that P(t+ P 1 (t+p 2 1(t P 2 (t+1. It is obvious that P(t+P 1 (t 1 e2(t T 1+e 2(t T which is not zero for t < T. Thus, following Theorem 1, the problem (IR-LQ is unsolvable. VII. CONCLUSIONS By proposing the so-called two-layer optimization approach, i.e., solving the FBDEs, the irregular LQ control problem has been converted into the controllability of a specified system. Then we have presented the solution to the open-loop and closed-loop control in this paper. In more specified, the closed-loop controller is designed by solving two Riccati equations and one linear equation, while the openloop solvability condition (necessary and sufficient is checked with the Range of two matrices and the open-loop controller is designed with a Gramian matrix and two Riccati equations. The presented results in this paper have clarified the essential difference of a irregular LQ control from a regular one, i.e., a irregular controller must guarantee the state terminal constraint of P 1 (Tx(T. We have also shown that the difference of a closed-loop control from a open-loop control, i.e., a closed-loop control is to chose the controller in the current state to make P 1 (Tx(T, while the open-loop control is to seek the controller in the form of initial state to reach of P 1 (Tx(T. REFERENCES 1 W. Wonham. On a matrix Riccati equation of stochastic control. SIAM Journal on Control and Optimization, 6(4, , D. L. Kleinman, Optimal stationary contorl of linear systems with control dependent noise. IEEE Trans. Autom. Control, 14(6, pp , W. L. Levison, D. L. Kleinman and S. Baron. A model for human controller remnant. Bolt, Beranek, and Newman, Inc., Cambridge, Mass., Tech. Reppt. pp. 1731, J. Baillieul and P. Antsaklis, Control and communication challenges in networked real-time systems. Proceeding of the IEEE. vol. 95, no. 1, pp. 9-28, J. M. Bismut. Linear quadratic optimal stochastic control with random coefficient. SIAM Journal on Control and Optimization, 14(3, pp , L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, E. F. Mishchenko. The mathematical theory of optimal process. English translation. Interscience, Bellman, Richard. The theory of dynamic programming. Bulletin of the American Mathematical Society, 6(6, , B. D. O. Anderson, J. B. Moore. Optimal control: linear quadratic methods. Englewood Cliffs, NJ: Prentice Hall, F. L. Lewis, D. L. Vrabie, V. L. Syrmos Optimal control. John Wiley & Sons, Inc., B. Oksendal. Stochastic Differential Equations: An Introduction with Applications. Springer, Berlin, H. Zhang, L. Lin, J. Xu, M. Fu, Linear quadratic regulation and stabilization of discrete-time Systems with delay and multiplicative noise, IEEE Transactions on Automatic Control, 6(1, , S. Chen, X. Li, X. Zhou, Stochastic linear quadratic regulators with indefinite control weight costs, SIAM Journal on Control and Optimization, 36(5, , 1998.

8 13 M. Ait Rami, J. B. Moore, and X. Y. Zhou, Indefinite stochastic linear quadratic control and generalized differential Riccati equation, SIAM Journal of Control and Optimization, 4, , M. Davis, Linear Estimation and Stochastic Control. Chapman and Hall, London, W. Wonham, On the separation theorem of stochastic control, SIAM Journal on Control, 6, , A. E. B. Lim and X. Y. Zhou, Stochastic optimal LQR control with integral quadratic constraints and indefinite control weights, IEEE Transactions on Automatic Control, 44, , X. Y. Zhou, D. Li, Continuous-time mean-variance portfoli selection: A stochastic LQ framework, Appl. Math. Optim., 42, 19-33, M. Kohlmann, X. Y. Zhou, Relationship between backward stochastic differential equations and stochastic controls: A linear-quadratic approach, SIAM Journal of Control and Optimization, 38, 1392C147, J. Sun, X. Li, J. Yong, Open-loop and closed-loop solvabilities for stochastic linear quadratic optimal control problems, SIAM Journal of Control and Optimization, 54(5, , Y-C Ho, Linear stochastic singular control problems, Journal of Optimization Theory and Applications, 9(1, 24-31, H-F. Chen, Unified controls applicable to general case under quadratic index, Acta Mathematicae Applicatae Sinica, 5(1, 45-52, M. Ait Rami, X. Chen, X. Y. Zhou, Discrete-time indefinite LQ control with state and control dependent noise, Journal of Global Optimization, 23, , R. Penrose, A generalized inverse of matrices, Mathematical Proceedings of the Cambridge Philosophical Society, 52, 17-19,

Consensus Control of Multi-agent Systems with Optimal Performance

Consensus Control of Multi-agent Systems with Optimal Performance 1 Consensus Control of Multi-agent Systems with Optimal Performance Juanjuan Xu, Huanshui Zhang arxiv:183.941v1 [math.oc 6 Mar 18 Abstract The consensus control with optimal cost remains major challenging

More information

Generalized Riccati Equations Arising in Stochastic Games

Generalized Riccati Equations Arising in Stochastic Games Generalized Riccati Equations Arising in Stochastic Games Michael McAsey a a Department of Mathematics Bradley University Peoria IL 61625 USA mcasey@bradley.edu Libin Mou b b Department of Mathematics

More information

Research Article Indefinite LQ Control for Discrete-Time Stochastic Systems via Semidefinite Programming

Research Article Indefinite LQ Control for Discrete-Time Stochastic Systems via Semidefinite Programming Mathematical Problems in Engineering Volume 2012, Article ID 674087, 14 pages doi:10.1155/2012/674087 Research Article Indefinite LQ Control for Discrete-Time Stochastic Systems via Semidefinite Programming

More information

Yuefen Chen & Yuanguo Zhu

Yuefen Chen & Yuanguo Zhu Indefinite LQ optimal control with equality constraint for discrete-time uncertain systems Yuefen Chen & Yuanguo Zhu Japan Journal of Industrial and Applied Mathematics ISSN 0916-7005 Volume 33 Number

More information

Chapter 2 Event-Triggered Sampling

Chapter 2 Event-Triggered Sampling Chapter Event-Triggered Sampling In this chapter, some general ideas and basic results on event-triggered sampling are introduced. The process considered is described by a first-order stochastic differential

More information

Non-zero Sum Stochastic Differential Games of Fully Coupled Forward-Backward Stochastic Systems

Non-zero Sum Stochastic Differential Games of Fully Coupled Forward-Backward Stochastic Systems Non-zero Sum Stochastic Differential Games of Fully Coupled Forward-Backward Stochastic Systems Maoning ang 1 Qingxin Meng 1 Yongzheng Sun 2 arxiv:11.236v1 math.oc 12 Oct 21 Abstract In this paper, an

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

More information

USE OF SEMIDEFINITE PROGRAMMING FOR SOLVING THE LQR PROBLEM SUBJECT TO RECTANGULAR DESCRIPTOR SYSTEMS

USE OF SEMIDEFINITE PROGRAMMING FOR SOLVING THE LQR PROBLEM SUBJECT TO RECTANGULAR DESCRIPTOR SYSTEMS Int. J. Appl. Math. Comput. Sci. 21 Vol. 2 No. 4 655 664 DOI: 1.2478/v16-1-48-9 USE OF SEMIDEFINITE PROGRAMMING FOR SOLVING THE LQR PROBLEM SUBJECT TO RECTANGULAR DESCRIPTOR SYSTEMS MUHAFZAN Department

More information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the

More information

On Solvability of an Indefinite Riccati Equation

On Solvability of an Indefinite Riccati Equation Working Paper Series National Centre of Competence in Research Financial Valuation and Risk Management Working Paper No. 883 On Solvability of an Indefinite Riccati Equation Kai Du First version: December

More information

A new robust delay-dependent stability criterion for a class of uncertain systems with delay

A new robust delay-dependent stability criterion for a class of uncertain systems with delay A new robust delay-dependent stability criterion for a class of uncertain systems with delay Fei Hao Long Wang and Tianguang Chu Abstract A new robust delay-dependent stability criterion for a class of

More information

Discrete-Time H Gaussian Filter

Discrete-Time H Gaussian Filter Proceedings of the 17th World Congress The International Federation of Automatic Control Discrete-Time H Gaussian Filter Ali Tahmasebi and Xiang Chen Department of Electrical and Computer Engineering,

More information

On linear quadratic optimal control of linear time-varying singular systems

On linear quadratic optimal control of linear time-varying singular systems On linear quadratic optimal control of linear time-varying singular systems Chi-Jo Wang Department of Electrical Engineering Southern Taiwan University of Technology 1 Nan-Tai Street, Yungkung, Tainan

More information

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,

More information

Tilburg University. An equivalence result in linear-quadratic theory van den Broek, W.A.; Engwerda, Jacob; Schumacher, Hans. Published in: Automatica

Tilburg University. An equivalence result in linear-quadratic theory van den Broek, W.A.; Engwerda, Jacob; Schumacher, Hans. Published in: Automatica Tilburg University An equivalence result in linear-quadratic theory van den Broek, W.A.; Engwerda, Jacob; Schumacher, Hans Published in: Automatica Publication date: 2003 Link to publication Citation for

More information

An Observation on the Positive Real Lemma

An Observation on the Positive Real Lemma Journal of Mathematical Analysis and Applications 255, 48 49 (21) doi:1.16/jmaa.2.7241, available online at http://www.idealibrary.com on An Observation on the Positive Real Lemma Luciano Pandolfi Dipartimento

More information

L 2 -induced Gains of Switched Systems and Classes of Switching Signals

L 2 -induced Gains of Switched Systems and Classes of Switching Signals L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit

More information

A New Algorithm for Solving Cross Coupled Algebraic Riccati Equations of Singularly Perturbed Nash Games

A New Algorithm for Solving Cross Coupled Algebraic Riccati Equations of Singularly Perturbed Nash Games A New Algorithm for Solving Cross Coupled Algebraic Riccati Equations of Singularly Perturbed Nash Games Hiroaki Mukaidani Hua Xu and Koichi Mizukami Faculty of Information Sciences Hiroshima City University

More information

Linear Quadratic Zero-Sum Two-Person Differential Games

Linear Quadratic Zero-Sum Two-Person Differential Games Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,

More information

UNIQUENESS OF SOLUTIONS TO MATRIX EQUATIONS ON TIME SCALES

UNIQUENESS OF SOLUTIONS TO MATRIX EQUATIONS ON TIME SCALES Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 50, pp. 1 13. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIQUENESS OF

More information

On the linear quadratic optimization problems and associated Riccati equations for systems modeled by Ito linear differential equations

On the linear quadratic optimization problems and associated Riccati equations for systems modeled by Ito linear differential equations Innovativity in Modeling and Analytics Journal of Research vol. 1, 216, pp.13-33, http://imajor.info/ c imajor ISSN XXXX-XXXX On the linear quadratic optimization problems and associated Riccati equations

More information

BESSEL MATRIX DIFFERENTIAL EQUATIONS: EXPLICIT SOLUTIONS OF INITIAL AND TWO-POINT BOUNDARY VALUE PROBLEMS

BESSEL MATRIX DIFFERENTIAL EQUATIONS: EXPLICIT SOLUTIONS OF INITIAL AND TWO-POINT BOUNDARY VALUE PROBLEMS APPLICATIONES MATHEMATICAE 22,1 (1993), pp. 11 23 E. NAVARRO, R. COMPANY and L. JÓDAR (Valencia) BESSEL MATRIX DIFFERENTIAL EQUATIONS: EXPLICIT SOLUTIONS OF INITIAL AND TWO-POINT BOUNDARY VALUE PROBLEMS

More information

Performance assessment of MIMO systems under partial information

Performance assessment of MIMO systems under partial information Performance assessment of MIMO systems under partial information H Xia P Majecki A Ordys M Grimble Abstract Minimum variance (MV) can characterize the most fundamental performance limitation of a system,

More information

LINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC COMPENSATION. Received October 2010; revised March 2011

LINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC COMPENSATION. Received October 2010; revised March 2011 International Journal of Innovative Computing, Information and Control ICIC International c 22 ISSN 349-498 Volume 8, Number 5(B), May 22 pp. 3743 3754 LINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC

More information

H 1 optimisation. Is hoped that the practical advantages of receding horizon control might be combined with the robustness advantages of H 1 control.

H 1 optimisation. Is hoped that the practical advantages of receding horizon control might be combined with the robustness advantages of H 1 control. A game theoretic approach to moving horizon control Sanjay Lall and Keith Glover Abstract A control law is constructed for a linear time varying system by solving a two player zero sum dierential game

More information

Optimal Polynomial Control for Discrete-Time Systems

Optimal Polynomial Control for Discrete-Time Systems 1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should

More information

Observer design for a general class of triangular systems

Observer design for a general class of triangular systems 1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals

More information

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery

Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Linear Quadratic Gausssian Control Design with Loop Transfer Recovery Leonid Freidovich Department of Mathematics Michigan State University MI 48824, USA e-mail:leonid@math.msu.edu http://www.math.msu.edu/

More information

Exam 2, Solution 4. Jordan Paschke

Exam 2, Solution 4. Jordan Paschke Exam 2, Solution 4 Jordan Paschke Problem 4. Suppose that y 1 (t) = t and y 2 (t) = t 3 are solutions to the equation 3t 3 y t 2 y 6ty + 6y = 0, t > 0. ( ) Find the general solution of the above equation.

More information

M. S. N. Murty, B. V. Appa Rao, and G. Suresh Kumar

M. S. N. Murty, B. V. Appa Rao, and G. Suresh Kumar Bull. Korean Math. Soc. 43 (2006), No. 1, pp. 149 159 CONTROLLABILITY, OBSERVABILITY, AND REALIZABILITY OF MATRIX LYAPUNOV SYSTEMS M. S. N. Murty, B. V. Appa Rao, and G. Suresh Kumar Abstract. This paper

More information

A note on linear differential equations with periodic coefficients.

A note on linear differential equations with periodic coefficients. A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.

More information

Time-Invariant Linear Quadratic Regulators!

Time-Invariant Linear Quadratic Regulators! Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting

More information

Hamilton-Jacobi-Bellman Equation Feb 25, 2008

Hamilton-Jacobi-Bellman Equation Feb 25, 2008 Hamilton-Jacobi-Bellman Equation Feb 25, 2008 What is it? The Hamilton-Jacobi-Bellman (HJB) equation is the continuous-time analog to the discrete deterministic dynamic programming algorithm Discrete VS

More information

Input-output finite-time stabilization for a class of hybrid systems

Input-output finite-time stabilization for a class of hybrid systems Input-output finite-time stabilization for a class of hybrid systems Francesco Amato 1 Gianmaria De Tommasi 2 1 Università degli Studi Magna Græcia di Catanzaro, Catanzaro, Italy, 2 Università degli Studi

More information

Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field

Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department

More information

The generalised continuous algebraic Riccati equation and impulse-free continuous-time LQ optimal control

The generalised continuous algebraic Riccati equation and impulse-free continuous-time LQ optimal control The generalised continuous algebraic iccati equation and impulse-free continuous-time LQ optimal control Augusto Ferrante Lorenzo Ntogramatzidis arxiv:135.5312v1 math.ds 23 May 213 Dipartimento di Ingegneria

More information

= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt +

= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt + ECE 553, Spring 8 Posted: May nd, 8 Problem Set #7 Solution Solutions: 1. The optimal controller is still the one given in the solution to the Problem 6 in Homework #5: u (x, t) = p(t)x k(t), t. The minimum

More information

On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method

On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method Ahmet Taha Koru Akın Delibaşı and Hitay Özbay Abstract In this paper we present a quasi-convex minimization method

More information

Research Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities

Research Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,

More information

Z i Q ij Z j. J(x, φ; U) = X T φ(t ) 2 h + where h k k, H(t) k k and R(t) r r are nonnegative definite matrices (R(t) is uniformly in t nonsingular).

Z i Q ij Z j. J(x, φ; U) = X T φ(t ) 2 h + where h k k, H(t) k k and R(t) r r are nonnegative definite matrices (R(t) is uniformly in t nonsingular). 2. LINEAR QUADRATIC DETERMINISTIC PROBLEM Notations: For a vector Z, Z = Z, Z is the Euclidean norm here Z, Z = i Z2 i is the inner product; For a vector Z and nonnegative definite matrix Q, Z Q = Z, QZ

More information

Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015

Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015 Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting

More information

The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation

The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse

More information

WHAT IS THE MINIMUM FUNCTION OBSERVER ORDER

WHAT IS THE MINIMUM FUNCTION OBSERVER ORDER WHAT IS THE MINIMUM FUNCTION OBSERVER ORDER Chia-Chi Tsui 743 Clove Road, NY 10310, USA, ctsui@ny.devry.edu poles (or the eigenvalues of F) of Keywords:function observer order, much lower & lowest possible

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May

More information

A conjecture on sustained oscillations for a closed-loop heat equation

A conjecture on sustained oscillations for a closed-loop heat equation A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping

More information

OPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28

OPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28 OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from

More information

Unbiased minimum variance estimation for systems with unknown exogenous inputs

Unbiased minimum variance estimation for systems with unknown exogenous inputs Unbiased minimum variance estimation for systems with unknown exogenous inputs Mohamed Darouach, Michel Zasadzinski To cite this version: Mohamed Darouach, Michel Zasadzinski. Unbiased minimum variance

More information

Positive Markov Jump Linear Systems (PMJLS) with applications

Positive Markov Jump Linear Systems (PMJLS) with applications Positive Markov Jump Linear Systems (PMJLS) with applications P. Bolzern, P. Colaneri DEIB, Politecnico di Milano - Italy December 12, 2015 Summary Positive Markov Jump Linear Systems Mean stability Input-output

More information

OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN

OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN Dynamic Systems and Applications 16 (2007) 393-406 OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN College of Mathematics and Computer

More information

Module 05 Introduction to Optimal Control

Module 05 Introduction to Optimal Control Module 05 Introduction to Optimal Control Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html October 8, 2015

More information

Sliding Mode Regulator as Solution to Optimal Control Problem for Nonlinear Polynomial Systems

Sliding Mode Regulator as Solution to Optimal Control Problem for Nonlinear Polynomial Systems 29 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -2, 29 WeA3.5 Sliding Mode Regulator as Solution to Optimal Control Problem for Nonlinear Polynomial Systems Michael Basin

More information

EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS

EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS Int. J. Appl. Math. Comput. Sci., 1, Vol.11, No.4, 957 964 EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS Tadeusz KACZOREK The notions of externally and internally positive time-varying

More information

An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game

An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game Marzieh Khakestari (Corresponding author) Institute For Mathematical Research, Universiti Putra Malaysia, 43400

More information

Noncausal Optimal Tracking of Linear Switched Systems

Noncausal Optimal Tracking of Linear Switched Systems Noncausal Optimal Tracking of Linear Switched Systems Gou Nakura Osaka University, Department of Engineering 2-1, Yamadaoka, Suita, Osaka, 565-0871, Japan nakura@watt.mech.eng.osaka-u.ac.jp Abstract. In

More information

(q 1)t. Control theory lends itself well such unification, as the structure and behavior of discrete control

(q 1)t. Control theory lends itself well such unification, as the structure and behavior of discrete control My general research area is the study of differential and difference equations. Currently I am working in an emerging field in dynamical systems. I would describe my work as a cross between the theoretical

More information

CONVERGENCE PROOF FOR RECURSIVE SOLUTION OF LINEAR-QUADRATIC NASH GAMES FOR QUASI-SINGULARLY PERTURBED SYSTEMS. S. Koskie, D. Skataric and B.

CONVERGENCE PROOF FOR RECURSIVE SOLUTION OF LINEAR-QUADRATIC NASH GAMES FOR QUASI-SINGULARLY PERTURBED SYSTEMS. S. Koskie, D. Skataric and B. To appear in Dynamics of Continuous, Discrete and Impulsive Systems http:monotone.uwaterloo.ca/ journal CONVERGENCE PROOF FOR RECURSIVE SOLUTION OF LINEAR-QUADRATIC NASH GAMES FOR QUASI-SINGULARLY PERTURBED

More information

Constrained controllability of semilinear systems with delayed controls

Constrained controllability of semilinear systems with delayed controls BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Constrained controllability of semilinear systems with delayed controls J. KLAMKA Institute of Control Engineering, Silesian

More information

Deterministic Dynamic Programming

Deterministic Dynamic Programming Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0

More information

Fixed-Order Robust H Filter Design for Markovian Jump Systems With Uncertain Switching Probabilities

Fixed-Order Robust H Filter Design for Markovian Jump Systems With Uncertain Switching Probabilities IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 4, APRIL 2006 1421 Fixed-Order Robust H Filter Design for Markovian Jump Systems With Uncertain Switching Probabilities Junlin Xiong and James Lam,

More information

Linear-Quadratic Optimal Control: Full-State Feedback

Linear-Quadratic Optimal Control: Full-State Feedback Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually

More information

Linear Quadratic Optimal Control Topics

Linear Quadratic Optimal Control Topics Linear Quadratic Optimal Control Topics Finite time LQR problem for time varying systems Open loop solution via Lagrange multiplier Closed loop solution Dynamic programming (DP) principle Cost-to-go function

More information

Complex Laplacians and Applications in Multi-Agent Systems

Complex Laplacians and Applications in Multi-Agent Systems 1 Complex Laplacians and Applications in Multi-Agent Systems Jiu-Gang Dong, and Li Qiu, Fellow, IEEE arxiv:1406.186v [math.oc] 14 Apr 015 Abstract Complex-valued Laplacians have been shown to be powerful

More information

Research Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field

Research Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field Complexity, Article ID 6235649, 9 pages https://doi.org/10.1155/2018/6235649 Research Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field Jinwang Liu, Dongmei

More information

Chapter 2 Optimal Control Problem

Chapter 2 Optimal Control Problem Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter

More information

A LaSalle version of Matrosov theorem

A LaSalle version of Matrosov theorem 5th IEEE Conference on Decision Control European Control Conference (CDC-ECC) Orlo, FL, USA, December -5, A LaSalle version of Matrosov theorem Alessro Astolfi Laurent Praly Abstract A weak version of

More information

A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems

A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems 53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,

More information

H Synchronization of Chaotic Systems via Delayed Feedback Control

H Synchronization of Chaotic Systems via Delayed Feedback Control International Journal of Automation and Computing 7(2), May 21, 23-235 DOI: 1.17/s11633-1-23-4 H Synchronization of Chaotic Systems via Delayed Feedback Control Li Sheng 1, 2 Hui-Zhong Yang 1 1 Institute

More information

Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control

Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control Outline Background Preliminaries Consensus Numerical simulations Conclusions Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control Email: lzhx@nankai.edu.cn, chenzq@nankai.edu.cn

More information

Balancing of Lossless and Passive Systems

Balancing of Lossless and Passive Systems Balancing of Lossless and Passive Systems Arjan van der Schaft Abstract Different balancing techniques are applied to lossless nonlinear systems, with open-loop balancing applied to their scattering representation.

More information

MODERN CONTROL DESIGN

MODERN CONTROL DESIGN CHAPTER 8 MODERN CONTROL DESIGN The classical design techniques of Chapters 6 and 7 are based on the root-locus and frequency response that utilize only the plant output for feedback with a dynamic controller

More information

Mean-field SDE driven by a fractional BM. A related stochastic control problem

Mean-field SDE driven by a fractional BM. A related stochastic control problem Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,

More information

Interval solutions for interval algebraic equations

Interval solutions for interval algebraic equations Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya

More information

Time-optimal control of a 3-level quantum system and its generalization to an n-level system

Time-optimal control of a 3-level quantum system and its generalization to an n-level system Proceedings of the 7 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 7 Time-optimal control of a 3-level quantum system and its generalization to an n-level

More information

Chapter III. Stability of Linear Systems

Chapter III. Stability of Linear Systems 1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,

More information

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,

More information

DIGITAL STABILIZATION OF LINEAR CONTINUOUS-TIME PERIODIC PROCESSES WITH PURE DELAY

DIGITAL STABILIZATION OF LINEAR CONTINUOUS-TIME PERIODIC PROCESSES WITH PURE DELAY DIGITAL STABILIZATION OF LINEAR CONTINUOUS-TIME PERIODIC PROCESSES WITH PURE DELAY B.P. Lampe E.N. Rosenwasser University of Rostock, Department of Computer Science and Electrical Engineering, D-18051

More information

Research Article Constrained Solutions of a System of Matrix Equations

Research Article Constrained Solutions of a System of Matrix Equations Journal of Applied Mathematics Volume 2012, Article ID 471573, 19 pages doi:10.1155/2012/471573 Research Article Constrained Solutions of a System of Matrix Equations Qing-Wen Wang 1 and Juan Yu 1, 2 1

More information

Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate

Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate www.scichina.com info.scichina.com www.springerlin.com Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate WEI Chen & CHEN ZongJi School of Automation

More information

The best generalised inverse of the linear operator in normed linear space

The best generalised inverse of the linear operator in normed linear space Linear Algebra and its Applications 420 (2007) 9 19 www.elsevier.com/locate/laa The best generalised inverse of the linear operator in normed linear space Ping Liu, Yu-wen Wang School of Mathematics and

More information

On Design of Reduced-Order H Filters for Discrete-Time Systems from Incomplete Measurements

On Design of Reduced-Order H Filters for Discrete-Time Systems from Incomplete Measurements Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 2008 On Design of Reduced-Order H Filters for Discrete-Time Systems from Incomplete Measurements Shaosheng Zhou

More information

A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1

A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,

More information

Exam in Systems Engineering/Process Control

Exam in Systems Engineering/Process Control Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total

More information

ELA THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES

ELA THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES Volume 22, pp. 480-489, May 20 THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES XUZHOU CHEN AND JUN JI Abstract. In this paper, we study the Moore-Penrose inverse

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

An infinite dimensional pursuit evasion differential game with finite number of players

An infinite dimensional pursuit evasion differential game with finite number of players An infinite dimensional pursuit evasion differential game with finite number of players arxiv:156.919v1 [math.oc] 2 Jun 215 Mehdi Salimi a Massimiliano Ferrara b a Center for Dynamics, Department of Mathematics,

More information

Advanced Mechatronics Engineering

Advanced Mechatronics Engineering Advanced Mechatronics Engineering German University in Cairo 21 December, 2013 Outline Necessary conditions for optimal input Example Linear regulator problem Example Necessary conditions for optimal input

More information

Optimal control of nonstructured nonlinear descriptor systems

Optimal control of nonstructured nonlinear descriptor systems Optimal control of nonstructured nonlinear descriptor systems TU Berlin DFG Research Center Institut für Mathematik MATHEON Workshop Elgersburg 19.02.07 joint work with Peter Kunkel Overview Applications

More information

Stability of interval positive continuous-time linear systems

Stability of interval positive continuous-time linear systems BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 66, No. 1, 2018 DOI: 10.24425/119056 Stability of interval positive continuous-time linear systems T. KACZOREK Białystok University of

More information

On Solving Large Algebraic. Riccati Matrix Equations

On Solving Large Algebraic. Riccati Matrix Equations International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University

More information

STABILIZATION OF LINEAR SYSTEMS VIA DELAYED STATE FEEDBACK CONTROLLER. El-Kébir Boukas. N. K. M Sirdi. Received December 2007; accepted February 2008

STABILIZATION OF LINEAR SYSTEMS VIA DELAYED STATE FEEDBACK CONTROLLER. El-Kébir Boukas. N. K. M Sirdi. Received December 2007; accepted February 2008 ICIC Express Letters ICIC International c 28 ISSN 1881-83X Volume 2, Number 1, March 28 pp. 1 6 STABILIZATION OF LINEAR SYSTEMS VIA DELAYED STATE FEEDBACK CONTROLLER El-Kébir Boukas Department of Mechanical

More information

Gramians based model reduction for hybrid switched systems

Gramians based model reduction for hybrid switched systems Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics

More information

Periodic Linear Systems

Periodic Linear Systems Periodic Linear Systems Lecture 17 Math 634 10/8/99 We now consider ẋ = A(t)x (1) when A is a continuous periodic n n matrix function of t; i.e., when there is a constant T>0 such that A(t + T )=A(t) for

More information

Nonplanar Periodic Solutions for Restricted N+1-body Problems

Nonplanar Periodic Solutions for Restricted N+1-body Problems Non-planar Periodic Solutions for Spatial Restricted N+1-body Problems Xiaoxiao Zhao, Fengying Li and Shiqing Zhang Sichuan University New Perspectives on the N-body Problem, BIRS, 2013 Construction Content

More information

Robotics. Control Theory. Marc Toussaint U Stuttgart

Robotics. Control Theory. Marc Toussaint U Stuttgart Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,

More information

Delay-dependent Stability Analysis for Markovian Jump Systems with Interval Time-varying-delays

Delay-dependent Stability Analysis for Markovian Jump Systems with Interval Time-varying-delays International Journal of Automation and Computing 7(2), May 2010, 224-229 DOI: 10.1007/s11633-010-0224-2 Delay-dependent Stability Analysis for Markovian Jump Systems with Interval Time-varying-delays

More information

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state

More information

On the projection onto a finitely generated cone

On the projection onto a finitely generated cone Acta Cybernetica 00 (0000) 1 15. On the projection onto a finitely generated cone Miklós Ujvári Abstract In the paper we study the properties of the projection onto a finitely generated cone. We show for

More information

MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem

MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear

More information

Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems

Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems Linear Riccati Dynamics, Constant Feedback, and Controllability in Linear Quadratic Control Problems July 2001 Revised: December 2005 Ronald J. Balvers Douglas W. Mitchell Department of Economics Department

More information

Research Article Mean Square Stability of Impulsive Stochastic Differential Systems

Research Article Mean Square Stability of Impulsive Stochastic Differential Systems International Differential Equations Volume 011, Article ID 613695, 13 pages doi:10.1155/011/613695 Research Article Mean Square Stability of Impulsive Stochastic Differential Systems Shujie Yang, Bao

More information