Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection

Size: px
Start display at page:

Download "Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection"

Transcription

1 Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute of Technology Guwahati, Indian Institute of Technology Bombay July 18, 2013 Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

2 Outline Problem formulation: dual/adjoint system Main results: well-posed interconnection and the regular case Ill-posed interconnection & the singular case (descriptor system) Zeros at infinity, inadmissible initial conditions Necessary and sufficient conditions for no zeros at infinity (main result) Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

3 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

4 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

5 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

6 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have impulsive initial conditions? Under what conditions will the closed-loop not have inadmissible initial conditions? al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

7 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have impulsive initial conditions? Under what conditions will the closed-loop not have inadmissible initial conditions? Well-posed interconnection of a system and its adjoint gives d dt x = Hx with H: al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

8 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have impulsive initial conditions? Under what conditions will the closed-loop not have inadmissible initial conditions? Well-posed interconnection of a system and its adjoint gives d dt x = Hx with H: a Hamiltonian matrix. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

9 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have impulsive initial conditions? Under what conditions will the closed-loop not have inadmissible initial conditions? Well-posed interconnection of a system and its adjoint gives d dt x = Hx with H: a Hamiltonian matrix. Analogue of Hamiltonian matrix for ill-posed case? al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

10 Main issues concerned in this talk When is the interconnection of a system and its adjoint (dual) a well-posed interconnection? If ill-posed, when is the interconnection autonomous? Interconnection of system and its dual: kind-of stationary trajectories. Can the interconnection cause closed-loop to have impulsive initial conditions? Under what conditions will the closed-loop not have inadmissible initial conditions? Well-posed interconnection of a system and its adjoint gives d dt x = Hx with H: a Hamiltonian matrix. Analogue of Hamiltonian matrix for ill-posed case? Link with other singular Hamiltonian pencils? Skew-Hermitian Hermitian pencil? (Mehl, Mehrmann, Meerbergen, Watkins) al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

11 Well-posed interconnection Systems B 1 and B 2 (with transfer functions G 1 & G 2 ) Call the interconnection of B 1 and B 2 d 1 u 1 y 1 G + 1 well-posed if + for each d 1, d 2 L loc 1, + y 2 u 2 there exist unique u 1, y 1, u 2 and y 2 L loc 1 G 2 such that d2 + Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

12 Well-posed interconnection Systems B 1 and B 2 (with transfer functions G 1 & G 2 ) d 1 + u 1 G 1 y 1 well-posed if + for each d 1, d 2 L loc 1, Call the interconnection of B 1 and B 2 y 2 u + 2 d2 G 2 + such that the laws (G 1, G 2 and Σ s) are satisfied. there exist unique u 1, y 1, u 2 and y 2 L loc 1 Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

13 Well-posed interconnection Systems B 1 and B 2 (with transfer functions G 1 & G 2 ) d 1 + u 1 G 1 y 1 well-posed if + for each d 1, d 2 L loc 1, Call the interconnection of B 1 and B 2 y 2 u + 2 d2 G 2 + such that the laws (G 1, G 2 and Σ s) are satisfied. there exist unique u 1, y 1, u 2 and y 2 L loc 1 Call the set of allowed trajectories of system 1 as behavior B 1 B 1 := {(u 1, y 1 ) y 1 = G 1 u 1 } Similarly, B 2 := {(y 2, u 2 ) y 2 = G 2 u 2 } Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

14 Well-posed interconnection Systems B 1 and B 2 (with transfer functions G 1 & G 2 ) d 1 + u 1 G 1 y 1 well-posed if + for each d 1, d 2 L loc 1, Call the interconnection of B 1 and B 2 y 2 u + 2 d2 G 2 + such that the laws (G 1, G 2 and Σ s) are satisfied. there exist unique u 1, y 1, u 2 and y 2 L loc 1 Call the set of allowed trajectories of system 1 as behavior B 1 Similarly, For this talk: System its behavior B 1 := {(u 1, y 1 ) y 1 = G 1 u 1 } B 2 := {(y 2, u 2 ) y 2 = G 2 u 2 } Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

15 Well-posed interconnection Systems B 1 and B 2 (with transfer functions G 1 & G 2 ) d 1 + u 1 G 1 y 1 well-posed if + for each d 1, d 2 L loc 1, Call the interconnection of B 1 and B 2 y 2 u + 2 d2 G 2 + such that the laws (G 1, G 2 and Σ s) are satisfied. there exist unique u 1, y 1, u 2 and y 2 L loc 1 Call the set of allowed trajectories of system 1 as behavior B 1 Similarly, For this talk: System its behavior B 1 := {(u 1, y 1 ) y 1 = G 1 u 1 } B 2 := {(y 2, u 2 ) y 2 = G 2 u 2 } Positive or negative feedback? Assume d 1 = 0 and d 2 = 0. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

16 Well-posed interconnection Interconnection of B 1 and B 2 : Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

17 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

18 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. In this talk, B 2 is the adjoint/dual system of B 1. What is the adjoint of a system (in behavioral sense)? Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

19 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. In this talk, B 2 is the adjoint/dual system of B 1. What is the adjoint of a system (in behavioral sense)? Consider behaviors B 1 and B 2 C (R, R w ) and Σ = B 1 and B 2 are called Σ-orthogonal if [ ] I 0. 0 I w T 1 Σw 2 dt = 0 for all w 1 B 1 and w 2 B 2. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

20 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. In this talk, B 2 is the adjoint/dual system of B 1. What is the adjoint of a system (in behavioral sense)? Consider behaviors B 1 and B 2 C (R, R w ) and Σ = B 1 and B 2 are called Σ-orthogonal if [ ] I 0. 0 I w T 1 Σw 2 dt = 0 for all w 1 B 1 and w 2 B 2. Above integrals could be undefined: restrict each of B 1 and B 2 to just the compactly supported trajectories, i.e. D. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

21 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. In this talk, B 2 is the adjoint/dual system of B 1. What is the adjoint of a system (in behavioral sense)? Consider behaviors B 1 and B 2 C (R, R w ) and Σ = B 1 and B 2 are called Σ-orthogonal if [ ] I 0. 0 I w T 1 Σw 2 dt = 0 for all w 1 B 1 and w 2 B 2. Above integrals could be undefined: restrict each of B 1 and B 2 to just the compactly supported trajectories, i.e. D. (Controllability assumed throughout this talk) al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

22 Well-posed interconnection Interconnection of B 1 and B 2 : trajectories need satisfy laws of both B 1 and B 2. Allowed trajectories of the closed loop system: intersection B 1 B 2. In this talk, B 2 is the adjoint/dual system of B 1. What is the adjoint of a system (in behavioral sense)? Consider behaviors B 1 and B 2 C (R, R w ) and Σ = B 1 and B 2 are called Σ-orthogonal if [ ] I 0. 0 I w T 1 Σw 2 dt = 0 for all w 1 B 1 and w 2 B 2. Above integrals could be undefined: restrict each of B 1 and B 2 to just the compactly supported trajectories, i.e. D. (Controllability assumed throughout this talk) For a controllable system B 1 C (R, R w ), define its Σ-orthogonal complement B Σ 1 as: B Σ 1 := {v C (R, R w ) w T Σv dt = 0 for all w B 1 D}. R al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

23 Hamiltonian matrix B : ẋ = Ax + Bw 1 w 2 = Cx + Dw 1 B Σ : w 1, v 1 : inputs, w 2, v 2 : outputs. Interconnect B and B Σ : w 2 = v 1 and w 1 = v 2 ẋ A BB T BD T ż = 0 A T C T 0 C DB T I p DD T B B Σ is well-posed (I p DD T ) is nonsingular. ż = A T z C T v 1 v 2 = B T z + D T v 1 x z v 1 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

24 Hamiltonian matrix B : ẋ = Ax + Bw 1 w 2 = Cx + Dw 1 B Σ : w 1, v 1 : inputs, w 2, v 2 : outputs. Interconnect B and B Σ : w 2 = v 1 and w 1 = v 2 ẋ A BB T BD T ż = 0 A T C T 0 C DB T I p DD T ż = A T z C T v 1 v 2 = B T z + D T v 1 B B Σ is well-posed [ẋ (I p ] DD [ T ) ] is nonsingular. x Eliminating v 1, we get I 2n = H with H R ż z 2n 2n the Hamiltonian matrix: [ ] A + BD H = T (I p DD T ) 1 C BB T + BD T (I p DD T ) 1 DB T C T (I p DD T ) 1 C (A T + C T (I p DD T ) 1 DB T. ) x z v 1 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

25 Hamiltonian matrix B : ẋ = Ax + Bw 1 w 2 = Cx + Dw 1 B Σ : w 1, v 1 : inputs, w 2, v 2 : outputs. Interconnect B and B Σ : w 2 = v 1 and w 1 = v 2 ẋ A BB T BD T ż = 0 A T C T 0 C DB T I p DD T ż = A T z C T v 1 v 2 = B T z + D T v 1 B B Σ is well-posed [ẋ (I p ] DD [ T ) ] is nonsingular. x Eliminating v 1, we get I 2n = H with H R ż z 2n 2n the Hamiltonian matrix: [ ] A + BD H = T (I p DD T ) 1 C BB T + BD T (I p DD T ) 1 DB T C T (I p DD T ) 1 C (A T + C T (I p DD T ) 1 DB T. ) With small variations (due to the Σ-matrix), x z v 1 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

26 Hamiltonian matrix B : ẋ = Ax + Bw 1 w 2 = Cx + Dw 1 B Σ : w 1, v 1 : inputs, w 2, v 2 : outputs. Interconnect B and B Σ : w 2 = v 1 and w 1 = v 2 ẋ A BB T BD T ż = 0 A T C T 0 C DB T I p DD T ż = A T z C T v 1 v 2 = B T z + D T v 1 B B Σ is well-posed [ẋ (I p ] DD [ T ) ] is nonsingular. x Eliminating v 1, we get I 2n = H with H R ż z 2n 2n the Hamiltonian matrix: [ ] A + BD H = T (I p DD T ) 1 C BB T + BD T (I p DD T ) 1 DB T C T (I p DD T ) 1 C (A T + C T (I p DD T ) 1 DB T. ) With small variations (due to the Σ-matrix), H arises in LQ, LQG, H 2 and H -control problems. For LQ control, above trajectories are stationary : Willems, CDC-92. x z v 1 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

27 Hamiltonian matrix B : ẋ = Ax + Bw 1 w 2 = Cx + Dw 1 B Σ : w 1, v 1 : inputs, w 2, v 2 : outputs. Interconnect B and B Σ : w 2 = v 1 and w 1 = v 2 ẋ A BB T BD T ż = 0 A T C T 0 C DB T I p DD T ż = A T z C T v 1 v 2 = B T z + D T v 1 B B Σ is well-posed [ẋ (I p ] DD [ T ) ] is nonsingular. x Eliminating v 1, we get I 2n = H with H R ż z 2n 2n the Hamiltonian matrix: [ ] A + BD H = T (I p DD T ) 1 C BB T + BD T (I p DD T ) 1 DB T C T (I p DD T ) 1 C (A T + C T (I p DD T ) 1 DB T. ) With small variations (due to the Σ-matrix), H arises in LQ, LQG, H 2 and H -control problems. For LQ control, above trajectories are stationary : Willems, CDC-92. Nonsingularity of (I p DD T ) required for Riccati (in)equality. al,kalaimani,belur (IIT Bombay/Guwahati) T ECC / 14 x z v 1

28 Extreme case: D = I We focus on the case when (I p DD T ) is singular: the ill-posed case. To understand better the ill-posed case, assume G is square and D = I. Also assume B has full column rank Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

29 Extreme case: D = I We focus on the case when (I p DD T ) is singular: the ill-posed case. To understand better the ill-posed case, assume G is square and D = I. Also assume B has full column rank (otherwise, can show, B B Σ is non-autonomous). Under what conditions is B B Σ autonomous? Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

30 Ill-posed interconnection B B Σ Assume D = I and B is full column rank. Define [ ] [ ] A BB T Ã := 0 A T, B B := C T, C := [ C B T ] (1) Theorem Consider the interconnection of the behaviors B and B Σ. Then the following are equivalent. 1 The interconnected system is autonomous. 2 Ce Ãt B is nonsingular for some t. 3 ker ( C B) ker ( CÃ B) ker ( CÃ2n 1 B) = 0 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

31 Inadmissible initial conditions For a descriptor system E d dtx = Ax, with det (se A) 0, some initial conditions x(0 ) could result in impulsive solutions x(t). Call these initial conditions inadmissible. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

32 Inadmissible initial conditions For a descriptor system E d dtx = Ax, with det (se A) 0, some initial conditions x(0 ) could result in impulsive solutions x(t). Call these initial conditions inadmissible. No inadmissible initial conditions (se A) has no zeros at infinity. (Like a polynomial matrix P (s) can have finite zeros, P (s) can also have zeros at s =.) Zeros at infinity generalized eigenvalue at. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

33 Inadmissible initial conditions For a descriptor system E d dtx = Ax, with det (se A) 0, some initial conditions x(0 ) could result in impulsive solutions x(t). Call these initial conditions inadmissible. No inadmissible initial conditions (se A) has no zeros at infinity. (Like a polynomial matrix P (s) can have finite zeros, P (s) can also have zeros at s =.) Zeros at infinity generalized eigenvalue at. When does an autonomous singular Hamiltonian system have inadmissible initial conditions? al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

34 Main result [ A BB T Ã := 0 A T ] [, B := B C T ], C := [ C B T ] (2) Theorem Assume the singular Hamiltonian system is autonomous. The following are equivalent: 1 There are no inadmissible initial conditions. 2 Ce Ãt B is nonsingular at t = 0. 3 ker ( C B) = 0. 4 rank (CB B T C T ) = p al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

35 Main result [ A BB T Ã := 0 A T ] [, B := B C T ], C := [ C B T ] (2) Theorem Assume the singular Hamiltonian system is autonomous. The following are equivalent: 1 There are no inadmissible initial conditions. 2 Ce Ãt B is nonsingular at t = 0. 3 ker ( C B) = 0. 4 rank (CB B T C T ) = p Recall our theorem for the interconnection of B and B Σ. The following are equivalent. 1 The interconnected system is autonomous. 2 Ce Ãt B is nonsingular for some t. 3 ker ( C B) ker ( CÃ B) ker ( CÃ2n 1 B) = 0 al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

36 Relation with all pass MIMO transfer matrix det (CB (CB) T ) 0 is opposite to requirement for being all-pass. G(s) is all-pass I G( jω) T G(jω) = 0 for each ω R. Assumption D = I and G is all-pass CB B T C T = 0, CAB + (CAB) T B T C T CB = 0, Notice that CB is the first moment of G(s) about s =. Thus a necessary condition on the first moment for G to be all-pass is: Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

37 Relation with all pass MIMO transfer matrix det (CB (CB) T ) 0 is opposite to requirement for being all-pass. G(s) is all-pass I G( jω) T G(jω) = 0 for each ω R. Assumption D = I and G is all-pass CB B T C T = 0, CAB + (CAB) T B T C T CB = 0, Notice that CB is the first moment of G(s) about s =. Thus a necessary condition on the first moment for G to be all-pass is: the skew-symmetric part of CB is zero. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

38 Relation with all pass MIMO transfer matrix det (CB (CB) T ) 0 is opposite to requirement for being all-pass. G(s) is all-pass I G( jω) T G(jω) = 0 for each ω R. Assumption D = I and G is all-pass CB B T C T = 0, CAB + (CAB) T B T C T CB = 0, Notice that CB is the first moment of G(s) about s =. Thus a necessary condition on the first moment for G to be all-pass is: the skew-symmetric part of CB is zero. In fact, the interconnection of B and B Σ is autonomous no all-pass subsystem. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

39 Relation with all pass MIMO transfer matrix det (CB (CB) T ) 0 is opposite to requirement for being all-pass. G(s) is all-pass I G( jω) T G(jω) = 0 for each ω R. Assumption D = I and G is all-pass CB B T C T = 0, CAB + (CAB) T B T C T CB = 0, Notice that CB is the first moment of G(s) about s =. Thus a necessary condition on the first moment for G to be all-pass is: the skew-symmetric part of CB is zero. In fact, the interconnection of B and B Σ is autonomous no all-pass subsystem. On the other hand, a singular Hamiltonian system (assumed autonomous) has no inadmissible initial conditions if and only if CB B T C T is nonsingular. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

40 Relation with all pass MIMO transfer matrix det (CB (CB) T ) 0 is opposite to requirement for being all-pass. G(s) is all-pass I G( jω) T G(jω) = 0 for each ω R. Assumption D = I and G is all-pass CB B T C T = 0, CAB + (CAB) T B T C T CB = 0, Notice that CB is the first moment of G(s) about s =. Thus a necessary condition on the first moment for G to be all-pass is: the skew-symmetric part of CB is zero. In fact, the interconnection of B and B Σ is autonomous no all-pass subsystem. On the other hand, a singular Hamiltonian system (assumed autonomous) has no inadmissible initial conditions if and only if CB B T C T is nonsingular. Odd number of inputs (and autonomous, ill-posed) there exist inadmissible initial conditions. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

41 Example G(s) := s+1 s+2 with input u and output y. State space realization : (A, B, C, D) = ( 2, 1, 1, 1). Dual system := s 1 s 2. State space realization : (A, B, C, D) = (2, 1, 1, 1). The interconnection of G and its dual: x x d z = z dt y y Above Es A has a zero at infinity. The differential equation in just x and z has initial conditions that have impulsive solutions. From our result too: CB B T C T = 0. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

42 Conclusion For the Σ we considered: B B Σ is an all-pass subsystem (possibly autonomous) In the context of Σ-dissipativity, this is set of stationary trajectories with respect to w T Σw. The stationary trajectories are interconnection of B and B Σ. Singular descriptor system if and only if ill-posed interconnection. Assuming D = I (the zeroth moment of G), no inadmissible initial conditions G s first moment has its skew-symmetric part nonsingular. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

43 Conclusion For the Σ we considered: B B Σ is an all-pass subsystem (possibly autonomous) In the context of Σ-dissipativity, this is set of stationary trajectories with respect to w T Σw. The stationary trajectories are interconnection of B and B Σ. Singular descriptor system if and only if ill-posed interconnection. Assuming D = I (the zeroth moment of G), no inadmissible initial conditions G s first moment has its skew-symmetric part nonsingular. Note: square MIMO all-pass G has first moment symmetric. No obvious way to have E as skew-hamitonian and H as Hamiltonian (as considered by Mehl, Mehrmann, et al.) But: for our pencil (E, H) also, generalized eigenvalues occur in quadruplets (λ, λ, λ, λ): like Mehl, Mehrmann and others. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

44 Questions, thank you! Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

45 Inadmissible initial condition Consider an autonomous system P ( d dt )w(t) = 0, with P (ξ) Rw w [ξ] nonsingular. A vector w(0) R zw is said to be an inadmissible initial condition vector if the corresponding solution w(t) contains the Dirac impulse δ(t) and/or its distributional derivatives. There exist no inadmissible initial conditions for P ( d dt )w = 0 P has no zeros at infinity. Pal,Kalaimani,Belur (IIT Bombay/Guwahati) ECC / 14

46 Stationary trajectories Consider a behavior B L w cont and a symmetric nonsingular matrix Σ R w w. A trajectory w B is Σ-stationary if Assume (wlog) that Σ := w T Σv dt = 0 for all v B D. [ ] Im 0 0 I p Definition a Given a controllable behavior B L w cont and Σ R w w, the Σ-orthogonal complement of B, denoted by B Σ is the set of all the trajectories v L loc 1 (R, R w ) such that vt Σw dt = 0 for all w B D. a J.C. Willems and H.L. Trentelman, On quadratic differential forms, SIAM Journal on Control and Optimization, The set of Σ-stationary trajectories is equal to B B Σ. al,kalaimani,belur (IIT Bombay/Guwahati) ECC / 14

The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications

The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department

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

Robust Control 2 Controllability, Observability & Transfer Functions

Robust Control 2 Controllability, Observability & Transfer Functions Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable

More information

The norms can also be characterized in terms of Riccati inequalities.

The norms can also be characterized in terms of Riccati inequalities. 9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements

More information

Denis ARZELIER arzelier

Denis ARZELIER   arzelier COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15

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

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

Linear Hamiltonian systems

Linear Hamiltonian systems Linear Hamiltonian systems P. Rapisarda H.L. Trentelman Abstract We study linear Hamiltonian systems using bilinear and quadratic differential forms. Such a representation-free approach allows to use the

More information

Zeros and zero dynamics

Zeros and zero dynamics CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)

More information

Stabilization and Passivity-Based Control

Stabilization and Passivity-Based Control DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive

More information

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and

More information

Control Systems. Frequency domain analysis. L. Lanari

Control Systems. Frequency domain analysis. L. Lanari Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic

More information

1 Linear Algebra Problems

1 Linear Algebra Problems Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and

More information

Dissipativity analysis of SISO systems using Nyquist-Plot-Compatible (NPC) supply rates

Dissipativity analysis of SISO systems using Nyquist-Plot-Compatible (NPC) supply rates Dissipativity analysis of SISO systems using Nyquist-Plot-Compatible (NPC) supply rates N. Santosh Kumar Debasattam Pal Abstract In this paper we deal with a special class of quadratic supply rates, dissipativity

More information

Iterative methods to compute center and center-stable manifolds with application to the optimal output regulation problem

Iterative methods to compute center and center-stable manifolds with application to the optimal output regulation problem Iterative methods to compute center and center-stable manifolds with application to the optimal output regulation problem Noboru Sakamoto, Branislav Rehak N.S.: Nagoya University, Department of Aerospace

More information

Citation for published version (APA): Minh, H. B. (2009). Model Reduction in a Behavioral Framework s.n.

Citation for published version (APA): Minh, H. B. (2009). Model Reduction in a Behavioral Framework s.n. University of Groningen Model Reduction in a Behavioral Framework Minh, Ha Binh IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please

More information

Problem Set 5 Solutions 1

Problem Set 5 Solutions 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel

More information

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear

More information

Controllability, Observability, Full State Feedback, Observer Based Control

Controllability, Observability, Full State Feedback, Observer Based Control Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time

More information

H 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University

H 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process

More information

Model reduction for linear systems by balancing

Model reduction for linear systems by balancing Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,

More information

Lecture 15: H Control Synthesis

Lecture 15: H Control Synthesis c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control

More information

Chap 4. State-Space Solutions and

Chap 4. State-Space Solutions and Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations

More information

POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK

POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK 1 POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK H.L. Trentelman, Member, IEEE, R. Zavala Yoe, and C. Praagman Abstract In this paper we will establish polynomial algorithms

More information

Problem Set 4 Solution 1

Problem Set 4 Solution 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 4 Solution Problem 4. For the SISO feedback

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

Intro. Computer Control Systems: F8

Intro. Computer Control Systems: F8 Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2

More information

Multivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006

Multivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006 Multivariable Control Lecture 3 Description of Linear Time Invariant Systems John T. Wen September 7, 26 Outline Mathematical description of LTI Systems Ref: 3.1-3.4 of text September 7, 26Copyrighted

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

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control

Chapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design

More information

CME 345: MODEL REDUCTION

CME 345: MODEL REDUCTION CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of

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

SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL

SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic

More information

Approximation of the Linearized Boussinesq Equations

Approximation of the Linearized Boussinesq Equations Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan

More information

CDS Solutions to the Midterm Exam

CDS Solutions to the Midterm Exam CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2

More information

Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach

Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach Bhaskar Ramasubramanian, Swanand R Khare and Madhu N Belur Abstract This paper formulates the problem

More information

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin

LQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system

More information

Stability, Pole Placement, Observers and Stabilization

Stability, Pole Placement, Observers and Stabilization Stability, Pole Placement, Observers and Stabilization 1 1, The Netherlands DISC Course Mathematical Models of Systems Outline 1 Stability of autonomous systems 2 The pole placement problem 3 Stabilization

More information

SYSTEMTEORI - ÖVNING Stability of linear systems Exercise 3.1 (LTI system). Consider the following matrix:

SYSTEMTEORI - ÖVNING Stability of linear systems Exercise 3.1 (LTI system). Consider the following matrix: SYSTEMTEORI - ÖVNING 3 1. Stability of linear systems Exercise 3.1 (LTI system. Consider the following matrix: ( A = 2 1 Use the Lyapunov method to determine if A is a stability matrix: a: in continuous

More information

DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS

DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS DISTANCE BETWEEN BEHAVIORS AND RATIONAL REPRESENTATIONS H.L. TRENTELMAN AND S.V. GOTTIMUKKALA Abstract. In this paper we study notions of distance between behaviors of linear differential systems. We introduce

More information

Robust Multivariable Control

Robust Multivariable Control Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms

More information

Introduction to Nonlinear Control Lecture # 4 Passivity

Introduction to Nonlinear Control Lecture # 4 Passivity p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive

More information

Optimization Based Output Feedback Control Design in Descriptor Systems

Optimization Based Output Feedback Control Design in Descriptor Systems Trabalho apresentado no XXXVII CNMAC, S.J. dos Campos - SP, 017. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics Optimization Based Output Feedback Control Design in

More information

LQR and H 2 control problems

LQR and H 2 control problems LQR and H 2 control problems Domenico Prattichizzo DII, University of Siena, Italy MTNS - July 5-9, 2010 D. Prattichizzo (Siena, Italy) MTNS - July 5-9, 2010 1 / 23 Geometric Control Theory for Linear

More information

Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control

Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control Electronic Journal of Linear Algebra Volume 34 Volume 34 08) Article 39 08 Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control Christian Mehl Technische Universitaet

More information

Suppose that we have a specific single stage dynamic system governed by the following equation:

Suppose that we have a specific single stage dynamic system governed by the following equation: Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where

More information

18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in

18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in 806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state

More information

Linear System Theory

Linear System Theory Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

Problem Set (T) If A is an m n matrix, B is an n p matrix and D is a p s matrix, then show

Problem Set (T) If A is an m n matrix, B is an n p matrix and D is a p s matrix, then show MTH 0: Linear Algebra Department of Mathematics and Statistics Indian Institute of Technology - Kanpur Problem Set Problems marked (T) are for discussions in Tutorial sessions (T) If A is an m n matrix,

More information

3 Gramians and Balanced Realizations

3 Gramians and Balanced Realizations 3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations

More information

Riccati Equations and Inequalities in Robust Control

Riccati Equations and Inequalities in Robust Control Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem

More information

Math 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:

Math 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations: Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =

More information

Introduction to Modern Control MT 2016

Introduction to Modern Control MT 2016 CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear

More information

EL2520 Control Theory and Practice

EL2520 Control Theory and Practice So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal

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

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael

More information

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA Kent State University Department of Mathematical Sciences Compiled and Maintained by Donald L. White Version: August 29, 2017 CONTENTS LINEAR ALGEBRA AND

More information

THE STABLE EMBEDDING PROBLEM

THE STABLE EMBEDDING PROBLEM THE STABLE EMBEDDING PROBLEM R. Zavala Yoé C. Praagman H.L. Trentelman Department of Econometrics, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands Research Institute for Mathematics

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1

More information

Network Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems

Network Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems Preprints of the 19th World Congress he International Federation of Automatic Control Network Reconstruction from Intrinsic Noise: Non-Minimum-Phase Systems David Hayden, Ye Yuan Jorge Goncalves Department

More information

OPTIMAL CONTROL SYSTEMS

OPTIMAL CONTROL SYSTEMS SYSTEMS MIN-MAX Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University OUTLINE MIN-MAX CONTROL Problem Definition HJB Equation Example GAME THEORY Differential Games Isaacs

More information

On Positive Real Lemma for Non-minimal Realization Systems

On Positive Real Lemma for Non-minimal Realization Systems Proceedings of the 17th World Congress The International Federation of Automatic Control On Positive Real Lemma for Non-minimal Realization Systems Sadaaki Kunimatsu Kim Sang-Hoon Takao Fujii Mitsuaki

More information

Lecture notes: Applied linear algebra Part 1. Version 2

Lecture notes: Applied linear algebra Part 1. Version 2 Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and

More information

Nonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México

Nonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October

More information

Problem Set 2 Solutions 1

Problem Set 2 Solutions 1 Massachusetts Institute of echnology Department of Electrical Engineering and Computer Science 6.45: MULIVARIABLE CONROL SYSEMS by A. Megretski Problem Set Solutions Problem. For each of the statements

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems . AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...

More information

Port-Hamiltonian Realizations of Linear Time Invariant Systems

Port-Hamiltonian Realizations of Linear Time Invariant Systems Port-Hamiltonian Realizations of Linear Time Invariant Systems Christopher Beattie and Volker Mehrmann and Hongguo Xu December 16 215 Abstract The question when a general linear time invariant control

More information

Topic # Feedback Control Systems

Topic # Feedback Control Systems Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control

More information

EL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)

EL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1) EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the

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

SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER

SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER Exercise 54 Consider the system: ẍ aẋ bx u where u is the input and x the output signal (a): Determine a state space realization (b): Is the

More information

Linear Algebra Practice Problems

Linear Algebra Practice Problems Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,

More information

Realization theory for systems biology

Realization theory for systems biology Realization theory for systems biology Mihály Petreczky CNRS Ecole Central Lille, France February 3, 2016 Outline Control theory and its role in biology. Realization problem: an abstract formulation. Realization

More information

Passivity Indices for Symmetrically Interconnected Distributed Systems

Passivity Indices for Symmetrically Interconnected Distributed Systems 9th Mediterranean Conference on Control and Automation Aquis Corfu Holiday Palace, Corfu, Greece June 0-3, 0 TuAT Passivity Indices for Symmetrically Interconnected Distributed Systems Po Wu and Panos

More information

Stabilization, Pole Placement, and Regular Implementability

Stabilization, Pole Placement, and Regular Implementability IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 5, MAY 2002 735 Stabilization, Pole Placement, and Regular Implementability Madhu N. Belur and H. L. Trentelman, Senior Member, IEEE Abstract In this

More information

Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach

Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach Generic degree structure of the minimal polynomial nullspace basis: a block Toeplitz matrix approach Bhaskar Ramasubramanian 1, Swanand R. Khare and Madhu N. Belur 3 December 17, 014 1 Department of Electrical

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear

More information

16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1

16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform

More information

Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications

Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications L. Faybusovich T. Mouktonglang Department of Mathematics, University of Notre Dame, Notre Dame, IN

More information

Some New Results on Linear Quadratic Regulator Design for Lossless Systems

Some New Results on Linear Quadratic Regulator Design for Lossless Systems Some New Results on Linear Quadratic Regulator Design for Lossless Systems Luigi Fortuna, Giovanni Muscato Maria Gabriella Xibilia Dipartimento Elettrico Elettronico e Sistemistico Universitá degli Studi

More information

Semidefinite Programming Duality and Linear Time-invariant Systems

Semidefinite Programming Duality and Linear Time-invariant Systems Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,

More information

MIT Final Exam Solutions, Spring 2017

MIT Final Exam Solutions, Spring 2017 MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of

More information

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution

More information

A BEHAVIORAL APPROACH TO THE LQ OPTIMAL CONTROL PROBLEM

A BEHAVIORAL APPROACH TO THE LQ OPTIMAL CONTROL PROBLEM A BEHAVIORAL APPROACH TO THE LQ OPTIMAL CONTROL PROBLEM GIANFRANCO PARLANGELI AND MARIA ELENA VALCHER Abstract This paper aims at providing new insights into the linear quadratic optimization problem within

More information

Singular-value-like decomposition for complex matrix triples

Singular-value-like decomposition for complex matrix triples Singular-value-like decomposition for complex matrix triples Christian Mehl Volker Mehrmann Hongguo Xu December 17 27 Dedicated to William B Gragg on the occasion of his 7th birthday Abstract The classical

More information

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10) Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must

More information

Lecture 4 Continuous time linear quadratic regulator

Lecture 4 Continuous time linear quadratic regulator EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon

More information

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation

More information

On at systems behaviors and observable image representations

On at systems behaviors and observable image representations Available online at www.sciencedirect.com Systems & Control Letters 51 (2004) 51 55 www.elsevier.com/locate/sysconle On at systems behaviors and observable image representations H.L. Trentelman Institute

More information

Maxwell Equations Waves Helmholtz Wave Equation Dirac Delta Sources Power Matrices Complex Analysis Optimization. EM & Math - Basics

Maxwell Equations Waves Helmholtz Wave Equation Dirac Delta Sources Power Matrices Complex Analysis Optimization. EM & Math - Basics EM & Math - S. R. Zinka srinivasa_zinka@daiict.ac.in October 16, 2014 Outline 1 Maxwell Equations 2 Waves 3 Helmholtz Wave Equation 4 Dirac Delta Sources 5 Power 6 Matrices 7 Complex Analysis 8 Optimization

More information

Eigenvalue perturbation theory of classes of structured matrices under generic structured rank one perturbations

Eigenvalue perturbation theory of classes of structured matrices under generic structured rank one perturbations Eigenvalue perturbation theory of classes of structured matrices under generic structured rank one perturbations, VU Amsterdam and North West University, South Africa joint work with: Leonard Batzke, Christian

More information

Control Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation

Control Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation Lecture 9: State Feedback and s [IFAC PB Ch 9] State Feedback s Disturbance Estimation & Integral Action Control Design Many factors to consider, for example: Attenuation of load disturbances Reduction

More information

Control Systems. Laplace domain analysis

Control Systems. Laplace domain analysis Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output

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

Intro. Computer Control Systems: F9

Intro. Computer Control Systems: F9 Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se

More information

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response .. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........

More information

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 15 Review of Matrix Theory III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Matrix An m n matrix is a rectangular or square array of

More information

APPLIED LINEAR ALGEBRA

APPLIED LINEAR ALGEBRA APPLIED LINEAR ALGEBRA Giorgio Picci November 24, 2015 1 Contents 1 LINEAR VECTOR SPACES AND LINEAR MAPS 10 1.1 Linear Maps and Matrices........................... 11 1.2 Inverse of a Linear Map............................

More information