Sampled-data implementation of derivative-dependent control using artificial delays

Size: px
Start display at page:

Download "Sampled-data implementation of derivative-dependent control using artificial delays"

Transcription

1 1 Sampled-data implementation of derivative-dependent control using artificial delays Anton Selivanov and Emilia Fridman, Senior Member, IEEE arxiv: v1 [mat.oc] 17 Jan 2018 Abstract We study a sampled-data implementation of linear controllers tat depend on te output and its derivatives. First, we consider an LTI system of relative degree r 2 tat can be stabilized using r 1 output derivatives. Ten, we consider PID control of a second order system. In bot cases, te Euler approximation is used for te derivatives giving rise to a delayed sampled-data controller. Given a derivative-dependent controller tat stabilizes te system, we sow ow to coose te parameters of te delayed sampled-data controller tat preserves te stability under fast enoug sampling. Te maximum sampling period is obtained from LMIs tat are derived using te Taylor s expansion of te delayed terms wit te remainders compensated by appropriate Lyapunov-Krasovskii functionals. Finally, we introduce te event-triggering mecanism tat may reduce te amount of sampled control signals used for stabilization. I. INTRODUCTION Control laws tat depend on output derivatives are used to stabilize systems wit relative degrees greater tan one. To estimate te derivatives, wic can ardly be measured directly, one can use te Euler approximation ẏ (y(t) y(t τ))/τ. Tis replaces te derivative-dependent control wit te delaydependent one [2] [5]. It as been sown in [6] tat suc approximation preserves te stability if τ > 0 is small enoug. Similarly, te output derivative in PID controller can be replaced by its Euler approximation. Te resulting controller was studied in [7] and [8] using te frequency domain analysis. In tis paper, we study sampled-data implementation of te delay-dependent controllers. For double-integrators, tis as been done in [9] using complete Lyapunov-Krasovskii functionals wit a Wirtinger-based term and in [10] via impulsive system representation and looped-functionals. Bot metods lead to complicated linear matrix inequalities (LMIs) containing many decision variables. In tis paper, we obtain simpler LMIs for more general systems and prove teir feasibility for small enoug sampling periods. A simple Lyapunov-based metod for delay-induced stabilization was proposed in [11], [12]. Te key idea is to use te Taylor s expansion of te delayed terms wit te remainders in te integral form tat are compensated by appropriate terms in te Lyapunov-Krasovskii functional. Tis leads to simple LMIs feasible for small delays if te derivative-dependent controller stabilizes te system. In tis paper, we study sampled-data implementation of two types of derivative-dependent controllers. In Section II, we consider an LTI system of relative degree r 2 tat can be Te autors are wit te Scool of Electrical Engineering, Tel Aviv University, Israel ( antonselivanov@gmail.com; emilia@eng.tau.ac.il). Supported by Israel Science Foundation (grant No. 1128/1). Preliminary results ave been presented in [1]. stabilized using r 1 output derivatives. In Section III, we consider PID control of a second order system. In bot cases, te Euler approximation is used for te derivatives giving rise to a delayed sampled-data controller. Assuming tat te derivativedependent controller exponentially stabilizes te system wit a decay rate α > 0, we sow ow to coose te parameters of its sampled-data implementation tat exponentially stabilizes te system wit any decay rate α < α if te sampling period is small enoug. Te maximum sampling period is obtained from LMIs tat are derived using te ideas of [11], [12]. Finally, we introduce te event-triggering mecanism tat may reduce te amount of sampled control signals used for stabilization [13] [17]. In te preliminary paper [1], we studied delayed sampled-data control for systems wit relative degree two. Notations: N 0 = N {0}, 1 r = [1, 1,..., 1] T R r, I l R l l is te identity matrix, stands for te Kronecker product, x = max{n N n x} for x R, col{a 1,..., a r } denotes te column vector composed from te vectors a 1,..., a r. For p R, f() = O( p ) if tere exist positive M and 0 suc tat f() M p for (0, 0 ). Auxiliary lemmas Lemma 1 (Exponential Wirtinger inequality [18]): Let f H 1 (a, b) be suc tat f(a) = 0 or f(b) = 0. Ten b a e2αt f T (t)w f(t) dt e 2 α (b a) (b a) 2 π 2 b a e2αt f T (t)w f(t) dt for any α R and 0 W R n n. Lemma 2 (Jensen s inequality [19]): Let ρ: [a, b] [0, ) and f : [a, b] R n be suc tat te integration concerned is well-defined. Ten for any 0 < Q R n n, [ b a ρ(s)f(s) ds ] T Q [ b a ρ(s)f(s) ds ] b a ρ(s) ds b a ρ(s)f T (s)qf(s) ds. II. DERIVATIVE-DEPENDENT CONTROL USING DISCRETE-TIME MEASUREMENTS Consider te LTI system ẋ(t) = Ax(t) + Bu(t), y(t) = Cx(t), wit relative degree r 2, i.e., x R n, u R m, y R l (1) CA i B = 0, i = 0, 1,..., r 2, CA r 1 B 0. (2) Relative degree is ow many times te output y(t) needs to be differentiated before te input u(t) appears explicitly. In particular, (2) implies y (i) = CA i x, i = 0, 1,..., r 1. (3)

2 2 To prove (3), note tat it is trivial for i = 0 and, if it as been proved for i < r 1, it olds for i + 1: y (i+1) = ( y (i)) (3) = (CA i x) = CA i [Ax + Bu] (2) = CA i+1 x. For LTI systems wit relative degree r, it is common to look for a stabilizing controller of te form u(t) = K 0 y(t) + K 1 ẏ(t) K r 1 y (r 1) (t) () wit K i R m l for i = 0,..., r 1. Remark 1: Te control law () essentially reduces te system s relative degree from r 2 to r = 1. Indeed, te transfer matrix of (1) as te form W (s) = βrsn r + +β n s n +α 1s n 1 + +α n wit β r = CA r 1 B 0. Taking u(t) = K 0u 0 (t)+k 1 u 0 (t)+ + K r 1u (r 1) 0 (t), one as ỹ(s) = (βrsn r + +β n)(k r 1 sr 1 + +K 0 ) s n +α 1s n 1 + +α n ũ 0 (s), were ỹ and ũ 0 are te Laplace transforms of y and u 0. If β r K r 1 0, te latter system as relative degree one. If it can be stabilized by u 0 = Ky ten (1) can be stabilized by () wit K i = K i K. Te controller () depends on te output derivatives, wic are ard to measure directly. Instead, te derivatives can be approximated by te finite-differences ẏ(t) y(t) y(t τ1) τ 1, ÿ(t) 1 τ 1 ( y(t) y(t τ1) τ 1 Tis leads to te delay-dependent control ) y(t τ1) y(t τ2) (τ 2 τ 1),... u(t) = K 0 y(t) + K 1 y(t τ 1 ) + + K r 1 y(t τ r 1 ), (5) were te gains K 0,..., K r 1 depend on te delays 0 < τ 1 < < τ r 1. If (1) can be stabilized by te derivative-dependent control (), ten it can be stabilized by te delayed control (5) wit small enoug delays [6]. In tis paper, we study te sampled-data implementation of (5): u(t) = K 0 y(t k )+ r 1 i=1 K iy(t k q i ), t [t k, t k+1 ), (6) were > 0 is a sampling period, t k = k, k N 0, are sampling instants, 0 < q 1 < < q r 1, q i N, are discretetime delays, and y(t) = 0 if t < 0. In te next section, we prove tat if (1) can be stabilized by te derivative-dependent controller (), ten it can be stabilized by te delayed sampled-data controller (6) wit a small enoug sampling period. Moreover, we sow ow to coose appropriate sampling period, controller gains K 0,..., K r 1, and discrete-time delays q 1,..., q r 1. A. Stability conditions Introduce te errors due to sampling δ 0 (t) = y(t k ) y(t), δ i (t) = y(t k q i ) y(t q i ), t [t k, t k+1 ), k N 0, were i = 1,..., r 1. Following [12], we employ Taylor s expansion wit te remainder in te integral form: y(t q i ) = r 1 y (j) (t) j=0 j! ( q i ) j + κ i (t), were κ i (t) = ( 1)r t (r 1)! t q (s t + q i i) r 1 y (r) (s) ds. Combining tese representations wit (3), we rewrite (6) as u = [K 0, K]M Cx + K 0 δ 0 + Kδ + Kκ, (7) were δ = col{δ 1,..., δ r 1 }, κ = col{κ 1,..., κ r 1 }, M =. I l ( q I l q 1I 1 ) 2 l 2! I l I l q 2I l ( q 2 ) 2 2! I l. I l q r 1I l ( q r 1 ) 2 2! I l. ( q 1 ) r 1 I (r 1)! l ( q 2 ) r 1 I (r 1)! l.... ( q r 1 ) r 1 I (r 1)! l C CA K = [ ] K 1 K 2 K r 1, C = Te closed-loop system (1), (6) takes te form. CA r 1 ẋ = Dx + BK 0 δ 0 + BKδ + BKκ, D = A + B[K 0, K]M C., Using (3), te closed-loop system (1), () can be written as Coosing. (8) (9) ẋ = Dx, D = A + B[ K0,..., K r 1 ] C. (10) [K 0, K 1,..., K r 1 ] = [ K 0,..., K r 1 ]M 1, (11) we obtain D = D. (Te Vandermonde-type matrix M is invertible, since te delays q i are different.) If (1), () is stable, D must be Hurwitz and (9) will be stable for zero δ0, δ, κ. Te following teorem provides LMIs guaranteeing tat δ 0, δ, and κ do not destroy te stability of (9). Teorem 1: Consider te LTI system (1) subject to (2). (i) For given sampling period > 0, discrete-time delays 0 < q 1 <... < q r 1, controller gains K 0,..., K r 1 R m l, and decay rate α > 0, let tere exist positivedefinite matrices P R n n, W 0, W i, R i R m m (i = 1,..., r 1) suc tat 1 Φ 0, were Φ = {Φ ij } is te symmetric matrix composed from Φ 11 = P D + D T P + 2αP + 2 e 2α r 1 i=0 (K ica) T W i (K i CA), Φ 12 = 1 T r P B, Φ 13 = 1 T r 1 P B, Φ 1 = (CA r 1 D) T H, Φ 2 = 1 r (CA r 1 B) T H, Φ 22 = diag{w { 0, e 2αq1 W 1,..., e 2αqr 1 W r 1 }, Φ 33 = (r!) 2 e diag 2αq 1 (q 1) R r 1,..., e 2αq r 1 (q r 1) R r r 1 }, Φ 3 = 1 r 1 (CA r 1 B) T H, Φ = H wit H = r 1 i=1 (q i) r K T i R ik i and D defined in (9). Ten te delayed sampled-data controller (6) exponentially stabilizes te system (1) wit te decay rate α. (ii) Let tere exist K0,..., K r 1 R m l suc tat te derivative-dependent controller () exponentially stabilizes (1) wit a decay rate α. Ten, te delayed sampleddata controller (6) wit K 0,..., K r 1 given by (11) wit 1 MATLAB codes for solving te LMIs are available at ttps://gitub.com/antonselivanov/tac18

3 3 Fig. 1. Event-triggering wit respect to te control signal M from (8) and q i = i 1 r 1 (i = 1,..., r 1) 2 exponentially stabilizes (1) wit any given decay rate α < α if te sampling period > 0 is small enoug. Proof is given in Appendix A. Remark 2: Teorem 1(ii) explicitly defines te controller parameters K 0, K i, q i (i = 1,..., r 1), wic depend on. To find appropriate sampling period, one sould reduce until te LMIs from (i) start to be feasible. B. Event-triggered control Event-triggered control allows to reduce te number of signals transmitted troug a communication network [13] [17]. Te idea is to transmit te signal only wen it canges a lot. Te event-triggering mecanism for measurements was implemented in [1] for te system (1), (6) wit relative degree r = 2. Here, we consider te system wit r 2 and introduce te event-triggering for control signals, since te output eventtriggering leads to complicated conditions (see Remark 3). Consider te system (Fig. 1) ẋ(t) = Ax(t) + Bû k, t [t k, t k+1 ), k N 0, y(t) = Cx(t), x R n, u R m, y R l, (12) were û k = u(t k ) if u(t k ) from (6) was transmitted and û k = û k 1 oterwise. Te signal u(t k ) is transmitted if its relative cange since te last transmission is large enoug, namely, if (u(t k ) û k 1 ) T Ω(u(t k ) û k 1 ) > σu T (t k )Ωu(t k ), (13) were σ [0, 1) and 0 < Ω R m m are event-triggering parameters. Tus, û 0 = u(t 0 ) and { u(tk ), (13) is true, û k = (1) û k 1, (13) is false. Teorem 2: Consider te system (12) subject to (2). For given sampling period > 0, discrete-time delays 0 < q 1 <... < q r 1, controller gains K 0,..., K r 1 R m l, eventtriggering tresold σ [0, 1), and decay rate α > 0, let tere exist positive-definite matrices P R n n, Ω, W 0, W i, R i R m m (i = 1,..., r 1) suc tat 3 Φ e 0, were Φ e = Φ P B σ([k 0, K]M C) T Ω 0 σ1 r Ω 0 σ1 r 1 Ω HCA r 1 B 0 Ω 0 0 σω 2 Note tat q i = i 1 r 1 N for small > 0 3 MATLAB codes for solving te LMIs are available at ttps://gitub.com/antonselivanov/tac18 wit M, K, C defined in (8) and Φ, H given in Teorem 1. Ten te event-triggered controller (6), (13), (1) exponentially stabilizes te system (12) wit te decay rate α. Proof is given in Appendix B. Remark 3: Te event-triggering mecanism (13), (1) is constructed wit respect to te control signal. Tis allows to reduce te workload of a controller-to-actuator network. To compensate te event-triggering error, we add (29) to V, wic leads to two additional block-columns and block-rows in te LMI (confer Φ of Teorem 1 and Φ e of Teorem 2). One can study te event-triggering mecanism wit respect to te measurements by replacing y(t k ), y(t k q i ) wit ŷ k = y(t k ) + e k, ŷ k qi = y(t k q i ) + e k qi in (6). Tis may reduce te workload of a sensor-to-controller network but would require to add expressions similar to (29) to V for eac error e k, e k q1,..., e k qr 1. Tis would lead to more complicated LMIs wit two additional block-columns and block-rows for eac error. We study te event-triggering mecanism wit respect to te control for simplicity. Remark : Taking Ω = ωi wit large ω > 0, one can sow tat Φ e 0 and Φ 0 are equivalent for σ = 0. Tis appens since te event-triggered control (6), (13), (1) wit σ = 0 degenerates into periodic sampled-data control (6). Terefore, an appropriate σ can be found by increasing its value from zero wile preserving te feasibility of te LMIs from Teorem 2. C. Example Consider te triple integrator... y = u, wic can be presented in te form (1) wit [ ] [ ] 00 A= 0 0 1, B =, C = [ ]. (15) Tese parameters satisfy (2) wit r = 3. Te derivativedependent control () wit K 0 = 2 10, K1 = 0.06, K2 = 0.32 stabilizes te system (1), (15). Te LMIs of Teorem 1 are feasible for = 0.0, q 1 = 30, q 2 = 60, α = 10 3, K , K , K , were K i are calculated using (11). Terefore, te delayed sampled-data controller (6) also stabilizes te system (1), (15). Consider now te system (12), (15). Te LMIs of Teorem 2 are feasible for = 0.02, σ = wit te same control gains K 0, K1, K2, delays q 1, q 2, and decay rate α. Tus, te event-triggered control (6), (13), (1) stabilizes te system (12), (15). Performing numerical simulations for 10 randomly cosen initial conditions x(0) 1, we find tat te eventtriggered control (6), (13), (1) requires to transmit on average 55.6 control signals during 100 seconds. Te amount of transmissions for te sampled-data control (6) is given by + 1 = Tus, te event-triggering mecanism reduces te workload of te controller-to-actuator network by almost 80% preserving te decay rate α. Note tat σ > 0 leads to a smaller sampling period. Terefore, te event-triggering mecanism requires to transmit more measurements troug sensor-to-controller network. However, te total workload of bot networks is reduced by over 37%. 100

4 III. EVENT-TRIGGERED PID CONTROL Consider te scalar system and te PID controller ÿ(t) + a 1 ẏ(t) + a 2 y(t) = bu(t) (16) u(t) = k p y(t) + k i y(s) ds + k d ẏ(t). (17) Here, we study sampled-data implementation of te PID controller (17) tat is obtained using te approximations 0 y(s) ds k 0 y(s) ds k 1 j=0 y(t j), ẏ(t) ẏ(t k ) y(t k) y(t k q ) t [t k, t k+1 ), q, were > 0 is a sampling period, t k = k, k N 0, are sampling instants, q N is a discrete-time delay, and y(t k q ) = 0 for k < q. Substituting tese approximations into (17), we obtain te sampled-data controller u(t) = k p y(t k ) + k i k 1 j=0 y(t j) + k d y(t k q ), t [t k, t k+1 ), k N 0, wit y(t k q ) = 0 for k < q and 0 (18) k p = k p + k d q, k i = k i, k d = k d q. (19) Similarly to Section II-B, we introduce te event-triggering mecanism to reduce te amount of transmitted control signals. Namely, we consider te system ÿ(t)+a 1 ẏ(t)+a 2 y(t) = bû k, t [t k, t k+1 ), k N 0, (20) were û k is te event-triggered control: û 0 = u(t 0 ), { u(tk ), if (22) is true, û k = û k 1, if (22) is false, wit u(t) from (18) and te event-triggering condition (21) (u(t k ) û k 1 ) 2 > σu 2 (t k ). (22) Here, σ [0, 1) is te event-triggering tresold. Remark 5: We consider te event-triggering mecanism wit respect to te control signal, since te event-triggering wit respect to te measurements ŷ k = y(t k ) + e k leads to an accumulating error in te integral term: k 0 y(s) ds k 1 j=0 ŷj = k 1 j=0 y(t j) + k 1 j=0 e j. A. Stability conditions To study te stability of (20) under te event-triggered PID control (18), (21), (22), we rewrite te closed-loop system in te state space. Let x 1 = y, x 2 = ẏ, and x 3 (t) = (t t k )y(t k ) + k 1 j=0 y(t j), t [t k, t k+1 ). Introduce te errors due to sampling v(t) = x(t k ) x(t) δ(t) = y(t k q ) y(t q), t [t k, t k+1 ), k N 0. Using Taylor s expansion for y(t q) wit te remainder in te integral form, we ave y(t k q ) = y(t q) + δ(t) = y(t) ẏ(t)q + κ(t) + δ(t), were κ(t) = (s t + q)ÿ(s) ds. t q Using tese representations in (18), we obtain u(t k ) = k p x 1 (t k ) + k i x 3 (t k ) + k d y(t k q ) = [k p + k d, qk d, k i ]x + [k p, 0, k i ]v + k d (κ + δ). (23) Introduce te event-triggering error e k = û k u(t k ). Ten te system (20) under te event-triggered PID control (21), (22), (23) can be presented as ẋ = Ax + A v v + Bk d (κ + δ) + Be k, y = Cx, for t [t k, t k+1 ), k N 0, were A = a 2 + b(k p + k d ) a 1 qbk d bk i, A v = bk p 0 bk i, B = 0 b, C = [ ] (2) (25) Note tat te integral term in (18) requires to introduce te error due to sampling v tat appears in (2) but was absent in (9). Te analysis of v is te key difference between Teorem 2 and te next result. Teorem 3: Consider te system (20). (i) For given sampling period > 0, discrete-time delay q > 0, controller gains k p, k i, k d, event-triggering tresold σ [0, 1), and decay rate α > 0, let tere exist positivedefinite matrices P, S R 3 3 and nonnegative scalars W, R, ω suc tat Ψ 0, were Ψ = {Ψ ij } is te symmetric matrix composed from Ψ 11 = P A + A T P + 2αP + [ ] W k 2 d 2 e 2α, Ψ 12 = P A v, Ψ13 = Ψ 1 = Ψ 15 = P B, Ψ 16 = k p+k d k p qk d ωσ, Ψ 26 = 0 ωσ, k i Ψ 17 = A T G, Ψ 22 = S, Ψ 27 = A T v G, Ψ 36 = Ψ 6 = ωσ, Ψ 37 = Ψ 7 = Ψ 57 = B T G, Ψ 33 = W π2 e 2αq, Ψ = R (q) e 2αq, 2 Ψ 55 = ω, Ψ 66 = ωσ, ] Ψ 77 = G, G = 2 e 2α S + Rkd 2(q)2 [ wit A, A v, B, C given in (25). Ten, te event-triggered PID controller (18), (21), (22) exponentially stabilizes te system (20) wit te decay rate α. (ii) Let tere exist k p, k i, k d suc tat te PID controller (17) exponentially stabilizes te system (16) wit a decay rate α. Ten, te event-triggered PID controller (18), (21), (22) wit k p, k i, k d given by (11) and q = 1 2 exponentially stabilizes te system (20) wit any given decay rate α < α if te sampling period > 0 and te event-triggering tresold σ [0, 1) are small enoug. Proof is given in Appendix C. MATLAB codes for solving te LMIs are available at ttps://gitub.com/antonselivanov/tac18 k i

5 5 Remark 6: Te event-triggered control (18), (21), (22) wit σ = 0 degenerates into sampled-data control (18). Terefore, Teorem 3 wit σ = 0 gives te stability conditions for te system (16) under te sampled-data PID control (18). Remark 7: Appropriate values of and σ can be found in a manner similar to Remarks 2 and. B. Example Following [8], we consider (16) wit a 1 = 8., a 2 = 0, b = Te system is not asymptotically stable if u = 0. Te PID controller (17) wit k p = 10, k i = 0, k d = 0.65 exponentially stabilizes it wit te decay rate α 10.. Teorem 3 wit σ = 0 (see Remark 6) guarantees tat te sampled-data PID controller (18) can acieve any decay rate α < α if te sampling period > 0 is small enoug. Since α is on te verge of stability, α close to α requires to use small. Tus, for α = 10.3, te LMIs of Teorem 3 are feasible wit 10 7, q = 272, and k p, k i, k d given by (19). To avoid small sampling period, we take α = 5. For σ = 0, α = 5 and eac q = 1, 2, 3,... we find te maximum sampling period > 0 suc tat te LMIs of Teorem 3 are feasible. Te largest corresponds to α = 5, σ = 0, q = 7, = , k p 29.76, k i = 0, k d 19.76, were k p, k i, k d are calculated using (19). Remark 6 implies tat te sampled-data PID controller (18) stabilizes (16). Teorem 3 remains feasible for α = 5, σ = , q = 7, = 10 3, k p 33.21, k i = 0, k d 23.21, were k p, k i, k d are calculated using (19). Tus, te eventtriggered PID control (18), (21), (22) exponentially stabilizes (20). Performing numerical simulations in a manner described in Section II-C, we find tat te event-triggered PID control requires to transmit on average 628. control signals during 10 seconds. Te sampled-data controller (18) requires = 2128 transmissions. Tus, te eventtriggering mecanism reduces te workload of te controllerto-actuator network by more tan 70%. Te total workload of bot networks is reduced by more tan 26%. wit (i) Consider te functional APPENDIX A PROOF OF THEOREM 1 V = V 0 + V δ0 + V δ + V κ (26) V 0 = x T P x, V δ0 = 2 e 2α t t k e 2α(t s) ẏ T (s)k0 T W 0 K 0 ẏ(s) ds t k e 2α(t s) [y(s) y(t k )] T K0 T W 0 K 0 [y(s) y(t k )] ds, V δ = 2 e 2α r 1 r 1 i=1 i=1 t k q e 2α(t s) i ẏ T (s)ki T W ik i ẏ(s) ds qi t k q e 2α(t s) i [y(s) y(t k q i )] T Ki T W ik i [y(s) y(t k q i )] ds, V κ = r 1 i=1 t q e 2α(t s) i (s t + q i ) r (y (r) (s)) T Ki T R ik i y (r) (s) ds. Te term V κ, introduced in [12], compensates Taylor s remainders κ i, wile V δ0 and V δ, introduced in [9], compensate te sampling errors δ 0 and δ. Te Wirtinger inequality (Lemma 1) implies V δ0 0 and V δ 0. Using (9) and (3), we obtain V 0 + 2αV 0 =2x T P [Dx+BK 0 δ 0 +BKδ+BKκ]+2αx T P x, V δ0 + 2αV δ0 = 2 e 2α x T (K 0 CA) T W 0 (K 0 CA)x δt 0 K0 T W 0 K 0 δ 0, V δ + 2αV δ = 2 e 2α r 1 i=1 xt (K i CA) T W i (K i CA)x r 1 i=1 e 2αqi δ T i KT i W ik i δ i. Using y (r) = CA r 1 ẋ (wic follows from (3)) and Jensen s inequality (Lemma 2) wit ρ(s) = (s t + q i ) r 1, we ave V κ + 2αV κ = r 1 i=1 (q i) r (y (r) (t)) T Ki T R ik i y (r) (t) r 1 i=1 r t q e 2α(t s) i (s t + q i ) r 1 (y (r) (s)) T Ki T R ik i y (r) (s) ds r 1 i=1 (q i) r ẋ T (A r 1 ) T C T Ki T R ik i CA r 1 ẋ r 1 (r!) 2 i=1 (q i) e 2αqi κ T r i KT i R ik i κ i. Summing up, we obtain were V + 2αV ϕ T Φϕ + ẋ T (CA r 1 ) T H(CA r 1 )ẋ, (27) ϕ = col{x, K 0 δ 0,..., K r 1 δ r 1, K 1 κ 1,..., K r 1 κ r 1 } (28) and Φ is obtained from Φ by removing te last block-column and block-row. Substituting (9) for ẋ and applying te Scur complement, we find tat Φ 0 guarantees V 2αV. Since V (t k ) V (t k ), te latter implies exponential stability of te system (9) and, terefore, (1), (6). (ii) Since q i = O( 1 r 1 ), [12, Lemma 2.1] guaranties M 1 = O( 1 r 1 ), wic implies K i = O( 1 r 1 ) for i = 0,..., r 1. Since D = D (wit D defined in (10)) and (1), () is exponentially stable wit te decay rate α, tere exists P > 0 suc tat P D + D T P + 2αP < 0 for any α < α. Coose W 0 = O( 1 r ), W i = O( 1 r ), and R i = O( 1 1 r ) for i = 1,..., r 1. Applying te Scur complement to Φ 0, we obtain P D + D T P + 2αP + O( 1 r ) < 0, wic olds for small > 0. Tus, (i) guarantees (ii). APPENDIX B PROOF OF THEOREM 2 Denote e k = û k u(t k ), k N 0. Te event-triggering mecanism (13), (1) guarantees 0 σu T (t k )Ωu(t k ) e T k Ωe k. (29) Substituting û k = u(t k )+e k into (12) and using (7), we obtain (cf. (9)) ẋ = Dx + BK 0 δ 0 + BKδ + BKκ + Be k, t [t k, t k+1 ) (30) wit D given in (9). Consider V from (26). Calculations similar to tose from te proof of Teorem 1 lead to (cf. (27)) V + 2αV (29) V + 2αV + σu T (t k )Ωu(t k ) e T k Ωe k ϕ T e Φ e ϕ e + ẋ T (CA r 1 ) T H(CA r 1 )ẋ + σu T (t k )Ωu(t k ),

6 6 were ϕ e = col{ϕ, e k } (wit ϕ from (28)) and Φ e is obtained from Φ e by removing te blocks Φ ij wit i {, 6} or j {, 6}. Substituting (30) for ẋ and (7) for u(t k ) and applying te Scur complement, we find tat Φ e 0 guarantees V 2αV. Since V (t k ) V (t k ), te latter implies exponential stability of te system (30) and, terefore, (12) under te controller (6), (13), (1). wit (i) Consider te functional V 0 = x T P x, V v = 2 e 2α APPENDIX C PROOF OF THEOREM 3 V = V 0 + V v + V δ + V κ t k e 2α(t s) ẋ T (s)sẋ(s) ds t k e 2α(t s) v T (s)sv(s) ds, V δ = W kd 22 e 2α t t k q e 2α(t s) (ẏ(s)) 2 ds W kd 2 π2 V κ = Rk 2 d q t k q e 2α(t s) [y(s) y(t k q)] 2 ds, t q e 2α(t s) (s t + q) 2 (ÿ(s)) 2 ds. t Te Wirtinger inequality (Lemma 1) implies V v 0 and V δ 0. Using te representation (2), we obtain V 0 + 2αV 0 = 2x T P [Ax(t) + A v v(t) + Bk d (κ + δ) + Be k ] +2αx T P x, V v + 2αV v = 2 e 2α ẋ T (t)sẋ(t) vt (t)sv(t), V δ + 2αV δ = W kd 22 e 2α (x 2 (t)) 2 W kd 2 π2 e 2αq δ 2 (t). Using Jensen s inequality (Lemma 2) wit ρ(s) = (s t+q i ), we obtain V κ + 2αV κ = Rk 2 d (q)2 (ÿ(t)) 2 Rkd 22 t q e 2α(t s) (s t + q)(ÿ(s)) 2 ds Rkd 2(q)2 (ẋ 2 (t)) 2 Rkd 2 (q) e 2αq κ 2 (t). 2 For ω 0, te event-triggering rule (21), (22) guarantees Tus, we ave 0 ωσu 2 (t k ) ωe 2 k. V + 2αV V + 2αV + [ωσu 2 (t k ) ωe 2 k ] ψ T Ψψ + ẋ T (t)gẋ(t) + ωσu 2 (t k ), were ψ = col{x, v/, k d δ, k d κ, e k } and Ψ is obtained from Ψ by removing te last two block-columns and block-rows. Substituting (2) for ẋ and (23) for u(t k ) and applying te Scur complement, we find tat Ψ 0 guarantees V 2αV. Since V (t k ) V (t k ), te latter implies exponential stability of te system (2) and, terefore, (18), (20) (22). (ii) Te closed-loop system (16), (17) is equivalent to ẋ = Āx wit y Ā = a 2 + b k p a 1 + b k d b k i, x = ẏ y(s) ds Since q = O( 1 ), relations (19) imply k p = O( 1 ), k d = O( 1 ). Since (16), (17) is exponentially stable wit te decay 0 rate α and (19) implies A = Ā, tere exists P > 0 suc tat P A+A T P +2αP < 0 for any α < α. Coose S = O( 1 ), W = O( 1 ), R = O( ), and ω = O( 1 ). Applying te Scur complement to Ψ 0, we obtain P A + A T P + 2αP + O( ) + σf < 0 wit some F independent of σ. Te latter olds for small > 0 and σ 0. Tus, (i) guarantees (ii). REFERENCES [1] A. Selivanov and E. Fridman, Simple conditions for sampled-data stabilization by using artificial delay, in 20t IFAC World Congress, 2017, pp [2] A. Ilcmann and C. J. Sangwin, Output feedback stabilisation of minimum pase systems by delays, Systems & Control Letters, vol. 52, no. 3-, pp , 200. [3] S. I. Niculescu and W. Miciels, Stabilizing a cain of integrators using multiple delays, IEEE Transactions on Automatic Control, vol. 9, no. 5, pp , 200. [] I. Karafyllis, Robust global stabilization by means of discrete-delay output feedback, Systems & Control Letters, vol. 57, no. 12, pp , [5] A. Ramírez, R. Sipai, S. Mondié, and R. R. Garrido, An Analytical Approac to Tuning of Delay-Based Controllers for LTI-SISO Systems, SIAM Journal on Control and Optimization, vol. 55, no. 1, pp , [6] M. Frenc, A. Ilcmann, and M. Mueller, Robust stabilization by linear output delay feedback, SIAM Journal on Control and Optimization, vol. 8, no., pp , [7] Y. Okuyama, Discretized PID Control and Robust Stabilization for Continuous Plants, in IFAC Proceedings Volumes, vol. 1, no. 2, 2008, pp [8] A. Ramírez, S. Mondié, R. Garrido, and R. Sipai, Design of Proportional-Integral-Retarded (PIR) Controllers for Second-Order LTI Systems, IEEE Transactions on Automatic Control, vol. 61, no. 6, pp , [9] K. Liu and E. Fridman, Wirtinger s inequality and Lyapunov-based sampled-data stabilization, Automatica, vol. 8, no. 1, pp , [10] A. Seuret and C. Briat, Stability analysis of uncertain sampled-data systems wit incremental delay using looped-functionals, Automatica, vol. 55, pp , [11] E. Fridman and L. Saiket, Delay-induced stability of vector secondorder systems via simple Lyapunov functionals, Automatica, vol. 7, pp , [12], Stabilization by using artificial delays: An LMI approac, Automatica, vol. 81, pp , [13] K. J. Åström and B. Bernardsson, Comparison of Periodic and Event Based Sampling for First-Order Stocastic Systems, in 1t IFAC World Congress, 1999, pp [1] P. Tabuada, Event-Triggered Real-Time Sceduling of Stabilizing Control Tasks, IEEE Transactions on Automatic Control, vol. 52, no. 9, pp , [15] W. P. M. H. Heemels, K. H. Joansson, and P. Tabuada, An introduction to event-triggered and self-triggered control, in IEEE Conference on Decision and Control, 2012, pp [16] W. P. M. H. Heemels, M. C. F. Donkers, and A. R. Teel, Periodic Event- Triggered Control for Linear Systems, IEEE Transactions on Automatic Control, vol. 58, no., pp , [17] A. Selivanov and E. Fridman, Event-Triggered H Control: A Switcing Approac, IEEE Transactions on Automatic Control, vol. 61, no. 10, pp , [18], Observer-based input-to-state stabilization of networked control systems wit large uncertain delays, Automatica, vol. 7, pp , [19] O. Solomon and E. Fridman, New stability conditions for systems wit distributed delays, Automatica, vol. 9, no. 11, pp , 2013.

Copyright 2012 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future

Copyright 2012 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future Copyrigt 212 IEEE. Personal use of tis material is permitted. Permission from IEEE must be obtained for all oter uses, in any current or future media, including reprinting/republising tis material for

More information

Stability and Stabilizability of Linear Parameter Dependent System with Time Delay

Stability and Stabilizability of Linear Parameter Dependent System with Time Delay Proceedings of te International MultiConference of Engineers and Computer Scientists 2008 Vol II IMECS 2008, 19-21 Marc, 2008, Hong Kong Stability and Stabilizability of Linear Parameter Dependent System

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

NETWORKED control systems (NCS), that are comprised

NETWORKED control systems (NCS), that are comprised 1 Event-Triggered H Control: a Switching Approach Anton Selivanov and Emilia Fridman, Senior Member, IEEE Abstract Event-triggered approach to networked control systems is used to reduce the workload of

More information

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter

More information

Copyright c 2008 Kevin Long

Copyright c 2008 Kevin Long Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula

More information

Zuerst erschienen in: IEEE Trans. on Autom. Control 44 (1999), S DOI: /

Zuerst erschienen in: IEEE Trans. on Autom. Control 44 (1999), S DOI: / Ilcmann, Acim ; Townley, Stuart: Adaptive Sampling Control of Hig-Gain Stabilizable Systems Zuerst erscienen in: IEEE Trans. on Autom. Control 44 (999), S.96-966 DOI:.9/9.793786 IEEE TRASACTIOS O AUTOMATIC

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

Exam 1 Review Solutions

Exam 1 Review Solutions Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),

More information

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t)) Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous

More information

The Verlet Algorithm for Molecular Dynamics Simulations

The Verlet Algorithm for Molecular Dynamics Simulations Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical

More information

Event-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems

Event-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems Event-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems Pavankumar Tallapragada Nikhil Chopra Department of Mechanical Engineering, University of Maryland, College Park, 2742 MD,

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

Networked Control With Stochastic Scheduling

Networked Control With Stochastic Scheduling IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 60, NO., NOVEMBER 205 307 Networked Control With Stochastic Scheduling Kun Liu, Emilia Fridman, and Karl Henrik Johansson Abstract This note develops the time-delay

More information

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1 Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.

More information

Finite Difference Methods Assignments

Finite Difference Methods Assignments Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation

More information

Lyapunov characterization of input-to-state stability for semilinear control systems over Banach spaces

Lyapunov characterization of input-to-state stability for semilinear control systems over Banach spaces Lyapunov caracterization of input-to-state stability for semilinear control systems over Banac spaces Andrii Mironcenko a, Fabian Wirt a a Faculty of Computer Science and Matematics, University of Passau,

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

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

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

arxiv: v1 [math.oc] 18 May 2018

arxiv: v1 [math.oc] 18 May 2018 Derivative-Free Optimization Algoritms based on Non-Commutative Maps * Jan Feiling,, Amelie Zeller, and Cristian Ebenbauer arxiv:805.0748v [mat.oc] 8 May 08 Institute for Systems Teory and Automatic Control,

More information

Complete Quadratic Lyapunov functionals using Bessel-Legendre inequality

Complete Quadratic Lyapunov functionals using Bessel-Legendre inequality Complete Quadratic Lyapunov functionals using Bessel-Legendre inequality Alexandre Seuret Frédéric Gouaisbaut To cite tis version: Alexandre Seuret Frédéric Gouaisbaut Complete Quadratic Lyapunov functionals

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

Math 312 Lecture Notes Modeling

Math 312 Lecture Notes Modeling Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a

More information

HOMEWORK HELP 2 FOR MATH 151

HOMEWORK HELP 2 FOR MATH 151 HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,

More information

Continuity and Differentiability of the Trigonometric Functions

Continuity and Differentiability of the Trigonometric Functions [Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te

More information

Analysis of an SEIR Epidemic Model with a General Feedback Vaccination Law

Analysis of an SEIR Epidemic Model with a General Feedback Vaccination Law Proceedings of te World Congress on Engineering 5 Vol WCE 5 July - 3 5 London U.K. Analysis of an ER Epidemic Model wit a General Feedback Vaccination Law M. De la en. Alonso-Quesada A.beas and R. Nistal

More information

Output Stabilization of Time-Varying Input Delay System using Interval Observer Technique

Output Stabilization of Time-Varying Input Delay System using Interval Observer Technique Output Stabilization of Time-Varying Input Delay System using Interval Observer Technique Andrey Polyakov a, Denis Efimov a, Wilfrid Perruquetti a,b and Jean-Pierre Richard a,b a - NON-A, INRIA Lille Nord-Europe

More information

7.1 Using Antiderivatives to find Area

7.1 Using Antiderivatives to find Area 7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 = Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff

More information

Quasiperiodic phenomena in the Van der Pol - Mathieu equation

Quasiperiodic phenomena in the Van der Pol - Mathieu equation Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van

More information

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error. Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error

More information

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a? Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c) Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

Global Output Feedback Stabilization of a Class of Upper-Triangular Nonlinear Systems

Global Output Feedback Stabilization of a Class of Upper-Triangular Nonlinear Systems 9 American Control Conference Hyatt Regency Riverfront St Louis MO USA June - 9 FrA4 Global Output Feedbac Stabilization of a Class of Upper-Triangular Nonlinear Systems Cunjiang Qian Abstract Tis paper

More information

MAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points

MAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x

More information

5.74 Introductory Quantum Mechanics II

5.74 Introductory Quantum Mechanics II MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT

More information

The total error in numerical differentiation

The total error in numerical differentiation AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

Delay-dependent stability and stabilization of neutral time-delay systems

Delay-dependent stability and stabilization of neutral time-delay systems INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL Int. J. Robust Nonlinear Control 2009; 19:1364 1375 Published online 6 October 2008 in Wiley InterScience (www.interscience.wiley.com)..1384 Delay-dependent

More information

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Input/output delay approach to robust sampled-data H control

Input/output delay approach to robust sampled-data H control Systems & Control Letters 54 (5) 71 8 www.elsevier.com/locate/sysconle Input/output delay approach to robust sampled-data H control E. Fridman, U. Shaked, V. Suplin Department of Electrical Engineering-Systems,

More information

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XL 2002 FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS by Anna Baranowska Zdzis law Kamont Abstract. Classical

More information

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014 Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.

More information

Generic maximum nullity of a graph

Generic maximum nullity of a graph Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n

More information

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015 Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes

More information

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!

Precalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here! Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.

More information

2.8 The Derivative as a Function

2.8 The Derivative as a Function .8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

Digital Filter Structures

Digital Filter Structures Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical

More information

Mass Lumping for Constant Density Acoustics

Mass Lumping for Constant Density Acoustics Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial

More information

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t). . Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

Eects of small delays on stability of singularly perturbed systems

Eects of small delays on stability of singularly perturbed systems Automatica 38 (2002) 897 902 www.elsevier.com/locate/automatica Technical Communique Eects of small delays on stability of singularly perturbed systems Emilia Fridman Department of Electrical Engineering

More information

4.2 - Richardson Extrapolation

4.2 - Richardson Extrapolation . - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence

More information

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to

More information

OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL

OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL Gang FENG Department of Mechanical and Biomedical Engineering City University of Hong Kong July 25, 2014 Department

More information

The Derivative as a Function

The Derivative as a Function Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )

More information

Stability of linear systems with general sawtooth delay

Stability of linear systems with general sawtooth delay IMA Journal of Matematical Control an Information Page of 8 oi:0.093/imamci/nq03 Stability of linear systems wit general sawtoot elay KUN LIU, VLADIMIR SUPLIN AND EMILIA FRIDMAN Department of Electrical

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

IN this paper, we control a robotic sensor network to estimate

IN this paper, we control a robotic sensor network to estimate IEEE TRANSACTIONS ON AUTOMATIC CONTROL, JUNE 2017 1 Distributed Active State Estimation wit User-Specified Accuracy Carles Freundlic, Student Member, IEEE, Soomin Lee, Member, IEEE, and Micael M. Zavlanos,

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems

New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems Systems & Control Letters 43 (21 39 319 www.elsevier.com/locate/sysconle New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems E. Fridman Department of Electrical

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

Continuous Stochastic Processes

Continuous Stochastic Processes Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling

More information

Complexity of Decoding Positive-Rate Reed-Solomon Codes

Complexity of Decoding Positive-Rate Reed-Solomon Codes Complexity of Decoding Positive-Rate Reed-Solomon Codes Qi Ceng 1 and Daqing Wan 1 Scool of Computer Science Te University of Oklaoma Norman, OK73019 Email: qceng@cs.ou.edu Department of Matematics University

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

The Priestley-Chao Estimator

The Priestley-Chao Estimator Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are

More information

Name: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).

Name: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ). Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate

More information

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,

More information

1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2

1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2 MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are

More information

2.11 That s So Derivative

2.11 That s So Derivative 2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point

More information

Robotic manipulation project

Robotic manipulation project Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

ROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF SINGULAR SYSTEMS

ROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF SINGULAR SYSTEMS INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 5 Number 3-4 Pages 480 487 c 2009 Institute for Scientific Computing and Information ROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF

More information

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Chapters 19 & 20 Heat and the First Law of Thermodynamics Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,

More information

MAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016

MAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016 MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0

More information

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of

More information

A method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps

A method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.

More information

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001 Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x

More information

ETNA Kent State University

ETNA Kent State University Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS

More information

IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang

IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS Shumei Mu Tianguang Chu and Long Wang Intelligent Control Laboratory Center for Systems and Control Department of Mechanics

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Exercises for numerical differentiation. Øyvind Ryan

Exercises for numerical differentiation. Øyvind Ryan Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can

More information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019 ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 6. Differential Calculus 6.. Differentiation from First Principles. In tis capter, we will introduce

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

160 Chapter 3: Differentiation

160 Chapter 3: Differentiation 3. Differentiation Rules 159 3. Differentiation Rules Tis section introuces a few rules tat allow us to ifferentiate a great variety of functions. By proving tese rules ere, we can ifferentiate functions

More information

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous

More information

A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES

A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties

More information