On perturbations in the leading coefficient matrix of time-varying index-1 DAEs

Similar documents
Bohl exponent for time-varying linear differential-algebraic equations

Stability analysis of differential algebraic systems

Nonlinear Control Systems

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is

Linear Differential Equations. Problems

6 Linear Equation. 6.1 Equation with constant coefficients

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

First and Second Order Differential Equations Lecture 4

Model reduction for linear systems by balancing

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

Chapter 2: First Order DE 2.4 Linear vs. Nonlinear DEs

LINEAR EQUATIONS OF HIGHER ORDER. EXAMPLES. General framework

Calculating the domain of attraction: Zubov s method and extensions

LMI Methods in Optimal and Robust Control

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

Nonuniform in time state estimation of dynamic systems

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

2. Metric Spaces. 2.1 Definitions etc.

Approximation of spectral intervals and associated leading directions for linear DAEs via smooth SVDs*

Ordinary Differential Equations II

Differential-algebraic equations. Control and Numerics V

Chapter III. Stability of Linear Systems

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

Exponential stability of families of linear delay systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

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).

1. Stochastic Processes and filtrations

= 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 +

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Presenter: Noriyoshi Fukaya

Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations

Jordan Canonical Form

Econ 508B: Lecture 5

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).

Problem List MATH 5173 Spring, 2014

Chap. 1. Some Differential Geometric Tools

Stabilization of persistently excited linear systems

Differential equations

Continuous Distributions

Controller Design for Robust Output Regulation of Regular Linear Systems

On a general definition of transition waves and their properties

Algebraic Properties of Solutions of Linear Systems

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Math Ordinary Differential Equations

Generalized transition waves and their properties

Backward Stochastic Differential Equations with Infinite Time Horizon

ODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

6.241 Dynamic Systems and Control

ON THE STABILITY OF THE QUASI-LINEAR IMPLICIT EQUATIONS IN HILBERT SPACES

Deterministic Dynamic Programming

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

8 Periodic Linear Di erential Equations - Floquet Theory

Weak Singular Optimal Controls Subject to State-Control Equality Constraints

EDP with strong anisotropy : transport, heat, waves equations

Nonlinear Control Lecture 5: Stability Analysis II

Sequential State Estimation (Crassidas and Junkins, Chapter 5)

Section , #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation.

(2) Let f(x) = a 2 x if x<2, 4 2x 2 ifx 2. (b) Find the lim f(x). (c) Find all values of a that make f continuous at 2. Justify your answer.

Equations with Separated Variables on Time Scales

Last lecture: Recurrence relations and differential equations. The solution to the differential equation dx

Asymptotic Stability by Linearization

DISCRETE GRONWALL LEMMA AND APPLICATIONS

Before you begin read these instructions carefully.

Math 273, Final Exam Solutions

Advanced Mechatronics Engineering

NOTES ON CALCULUS OF VARIATIONS. September 13, 2012

1+t 2 (l) y = 2xy 3 (m) x = 2tx + 1 (n) x = 2tx + t (o) y = 1 + y (p) y = ty (q) y =

4. Conditional risk measures and their robust representation

{ } is an asymptotic sequence.

Some results in Floquet theory, with application to periodic epidemic models

Functions of a Complex Variable (S1) Lecture 11. VII. Integral Transforms. Integral transforms from application of complex calculus

Impulse free solutions for switched differential algebraic equations

A note on linear differential equations with periodic coefficients.

Periodic Linear Systems

Journal of Computational and Applied Mathematics

11.2 Basic First-order System Methods

On Controllability of Linear Systems 1

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

Functional Analysis. Martin Brokate. 1 Normed Spaces 2. 2 Hilbert Spaces The Principle of Uniform Boundedness 32

Observability and state estimation

Numerical Treatment of Unstructured. Differential-Algebraic Equations. with Arbitrary Index

Course 212: Academic Year Section 1: Metric Spaces

Solution of Additional Exercises for Chapter 4

Stability of Stochastic Differential Equations

Chap. 3. Controlled Systems, Controllability

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

A Bang-Bang Principle of Time Optimal Internal Controls of the Heat Equation

Problem 1 In each of the following problems find the general solution of the given differential

EML Spring 2012

Lecture 19 Observability and state estimation

Ordinary Differential Equation Theory

Nonlinear Systems and Control Lecture # 19 Perturbed Systems & Input-to-State Stability

STAT/MATH 395 A - PROBABILITY II UW Winter Quarter Moment functions. x r p X (x) (1) E[X r ] = x r f X (x) dx (2) (x E[X]) r p X (x) (3)

Nonlinear Systems Theory

1 Lyapunov theory of stability

PREPRINT 2009:10. Continuous-discrete optimal control problems KJELL HOLMÅKER

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

Transcription:

On perturbations in the leading coefficient matrix of time-varying index-1 DAEs Institute of Mathematics, Ilmenau University of Technology Darmstadt, March 27, 2012 Institute of Mathematics, Ilmenau University of Technology Page 1 / 11

Perturbations in the leading coefficient of DAEs E(t)ẋ = A(t)x, E, A C(R 0 ; R n n ) Institute of Mathematics, Ilmenau University of Technology Page 2 / 11

Perturbations in the leading coefficient of DAEs E(t)ẋ = A(t)x, E, A C(R 0 ; R n n ) [März 1991]: (E, A) is index-1 Q C 1 : Q(t) 2 = Q(t) im Q(t) = ker E(t), and D C 0 : Eẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x P = I Q Institute of Mathematics, Ilmenau University of Technology Page 2 / 11

Perturbations in the leading coefficient of DAEs E(t)ẋ = A(t)x, E, A C(R 0 ; R n n ) [März 1991]: (E, A) is index-1 Q C 1 : Q(t) 2 = Q(t) im Q(t) = ker E(t), and D C 0 : Eẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x P = I Q ẋ 1 = x 2 0 = x 1 ẋ1 = D 11 x 1 x 2 = D 21 x 1 Institute of Mathematics, Ilmenau University of Technology Page 2 / 11

Eẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x Institute of Mathematics, Ilmenau University of Technology Page 3 / 11

Eẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x Eẋ = Ax d dt x = P x + Qx = (I + QD)P x, (P x) = ( P + P D)P x Institute of Mathematics, Ilmenau University of Technology Page 3 / 11

Eẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x Eẋ = Ax d dt x = P x + Qx = (I + QD)P x, (P x) = ( P + P D)P x ẏ = ( P + P D)y, y(t 0 ) = P (t 0 )x(t 0 ) uniqueness = y(t) = P (t)x(t) Institute of Mathematics, Ilmenau University of Technology Page 3 / 11

crucial: E = EP Eẋ = Ax EP ẋ = Ax d dt (P x) = ( P + P D)P x, Qx = QDP x Institute of Mathematics, Ilmenau University of Technology Page 4 / 11

crucial: E = EP Eẋ = (A + A )x EP ẋ = (A + A )x d dt (P x) = ( P + P D)P x + P G A x, Qx = QDP x + QG A x G 1 Institute of Mathematics, Ilmenau University of Technology Page 4 / 11

crucial: E = EP i.g. (E + E )ẋ = Ax (E + E )P ẋ = Ax Institute of Mathematics, Ilmenau University of Technology Page 4 / 11

crucial: E = EP i.g. (E + E )ẋ = Ax (E + E )P ẋ = Ax (A): E C(R 0 ; R n n ) s.t. (E + E, A) is index-1 and ker E(t) = ker(e(t) + E (t)) Institute of Mathematics, Ilmenau University of Technology Page 4 / 11

crucial: E = EP i.g. (E + E )ẋ = Ax (E + E )P ẋ = Ax (A): E C(R 0 ; R n n ) s.t. (E + E, A) is index-1 and ker E(t) = ker(e(t) + E (t)) (E + E )ẋ = Ax d dt (P x) = ( P + P D)P x + P G P x, Qx = QDP x + QG P x Institute of Mathematics, Ilmenau University of Technology Page 4 / 11

[ ] ) [ ] ( ) In1 0 (ẋ1 A11 (t) A = 12 (t) x1 0 0 A 21 (t) A 22 (t) ẋ 2 x 2 Institute of Mathematics, Ilmenau University of Technology Page 5 / 11

[ ] ) [ ] ( ) In1 0 (ẋ1 A11 (t) A = 12 (t) x1 (t) 0 A 21 (t) A 22 (t) ẋ 2 0 = (A 21 A 11 )(t)x 1 + (A 22 A 12 )(t)x 2 x 2 Institute of Mathematics, Ilmenau University of Technology Page 5 / 11

[ ] ) [ ] ( ) In1 0 (ẋ1 A11 (t) A = 12 (t) x1 (t) 0 A 21 (t) A 22 (t) ẋ 2 0 = (A 21 A 11 )(t)x 1 + (A 22 A 12 )(t)x 2 [ ] ) 1 0 (ẋ1 = 0 0 ẋ 2 [ 1 0 0 1 ] ( ) x1 x 2 x 2 x 1 (t) = e t x 0 1 x 2 (t) = 0 exp. stable Institute of Mathematics, Ilmenau University of Technology Page 5 / 11

[ ] ) [ ] ( ) In1 0 (ẋ1 A11 (t) A = 12 (t) x1 (t) 0 A 21 (t) A 22 (t) ẋ 2 0 = (A 21 A 11 )(t)x 1 + (A 22 A 12 )(t)x 2 [ ] ) 1 0 (ẋ1 = 0 ε ẋ 2 [ 1 0 0 1 ] ( ) x1 x 2 x 2 x 1 (t) = e t x 0 1 x 2 (t) = e t/ε x 0 2 not exp. stable Institute of Mathematics, Ilmenau University of Technology Page 5 / 11

Bohl exponent and perturbation operator E(t) d dt Φ(t, t 0) = A(t)Φ(t, t 0 ), P (t 0 )(Φ(t 0, t 0 ) I) = 0. (E, A) exp. stab. µ, M > 0 t t 0 : Φ(t, t 0 ) Me µ(t t 0) Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 6 / 11

Bohl exponent and perturbation operator E(t) d dt Φ(t, t 0) = A(t)Φ(t, t 0 ), P (t 0 )(Φ(t 0, t 0 ) I) = 0. (E, A) exp. stab. µ, M > 0 t t 0 : Φ(t, t 0 ) Me µ(t t 0) k B (E, A) = inf ρ R N ρ > 0 t s : Φ(t, s) N ρ e ρ(t s) } Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 6 / 11

Bohl exponent and perturbation operator E(t) d dt Φ(t, t 0) = A(t)Φ(t, t 0 ), P (t 0 )(Φ(t 0, t 0 ) I) = 0. (E, A) exp. stab. µ, M > 0 t t 0 : Φ(t, t 0 ) Me µ(t t 0) k B (E, A) = inf ρ R N ρ > 0 t s : Φ(t, s) N ρ e ρ(t s) } Rem.: k B (E, A) < 0 (E, A) exp. stable Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 6 / 11

Bohl exponent and perturbation operator E(t) d dt Φ(t, t 0) = A(t)Φ(t, t 0 ), P (t 0 )(Φ(t 0, t 0 ) I) = 0. (E, A) exp. stab. µ, M > 0 t t 0 : Φ(t, t 0 ) Me µ(t t 0) k B (E, A) = inf ρ R N ρ > 0 t s : Φ(t, s) N ρ e ρ(t s) } Rem.: k B (E, A) < 0 (E, A) exp. stable L t0 : L 2 ([t 0, ); R n ) L 2 ([t 0, ); R n ), f( ) x( ), x solves E(t)ẋ = A(t)x + f(t), P (t 0 )x(t 0 ) = 0 Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 6 / 11

Bohl exponent and perturbation operator E(t) d dt Φ(t, t 0) = A(t)Φ(t, t 0 ), P (t 0 )(Φ(t 0, t 0 ) I) = 0. (E, A) exp. stab. µ, M > 0 t t 0 : Φ(t, t 0 ) Me µ(t t 0) k B (E, A) = inf ρ R N ρ > 0 t s : Φ(t, s) N ρ e ρ(t s) } Rem.: k B (E, A) < 0 (E, A) exp. stable L t0 : L 2 ([t 0, ); R n ) L 2 ([t 0, ); R n ), f( ) x( ), x solves E(t)ẋ = A(t)x + f(t), P (t 0 )x(t 0 ) = 0 Lemma [Du et al. 2006]: (E, A) exp. stable and (BC) hold = L t0 is linear bd. operator and t 0 L t0 is mon. nonincreasing Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 6 / 11

Theorem (Robustness of Bohl exponent) (E, A) index-1, Q bounded, given ε > 0: E satisfies (A), E suff. small = k B (E + E, A) k B (E, A) + ε Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 7 / 11

Theorem (Robustness via perturbation operator) (E, A) index-1 and exp. stable, (BC) hold, E satisfies (A): } κ i = κ i (E, A, Q), i = 1, 2, 3, α := min lim L t t 0 0 1, κ 3, lim E [t0, ) < t 0 = (E + E, A) is exponentially stable α κ 1 + κ 2 α Perturbation theory Institute of Mathematics, Ilmenau University of Technology Page 8 / 11

Stability radius P := [ E, A ] B(R 0 ; R n 2n ) inf [ E, A ] P (E + E, A + A ) is index-1, ker E(t) = ker(e(t) + E (t)) S := (E, A) C(R 0 ; R n n ) } 2 (E, A) is exponentially stable, [ r(e, A) := [ E, A ] E, A ] P or (E + E, A + A ) S }, } Stability radius Institute of Mathematics, Ilmenau University of Technology Page 9 / 11

Proposition (Properties of the stability radius) r(e, A) = 0 (E, A) S r(α(e, A)) = r(αe, αa) = α r(e, A) for all α 0 V(t) time-varying subspace of R n with constant dimension, K V := [E, A] B(R 0 ; R n 2n ) (E, A) is index-1, ker E(t) = V(t) = K V [E, A] r(e, A) is continuous }, Stability radius Institute of Mathematics, Ilmenau University of Technology Page 10 / 11

Theorem (Lower bound for the stability radius) (E, A) index-1 and exp. stable, (BC) holds = } κ i = κ i (E, A, Q), i = 1, 2, 3, α := min lim L t t 0 0 1, κ 3, α r(e, A) κ 1 + κ 2 α Stability radius Institute of Mathematics, Ilmenau University of Technology Page 11 / 11

Theorem (Lower bound for the stability radius) (E, A) index-1 and exp. stable, (BC) holds = } κ i = κ i (E, A, Q), i = 1, 2, 3, α := min lim L t t 0 0 1, κ 3, α r(e, A) κ 1 + κ 2 α Cor.: V(t) time-varying subspace of R n with constant dimension, S V := [E, A] B(R 0 ; R n 2n ) (E, A) is index-1 and exp. stable, ker E(t) = V(t) for all t R 0 } = S V is open in K V Stability radius Institute of Mathematics, Ilmenau University of Technology Page 11 / 11

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback