Order Preserving Properties of Vehicle Dynamics with Respect to the Driver s Input
|
|
- Alexis Daisy Carter
- 6 years ago
- Views:
Transcription
1 Order Preserving Properties of Vehicle Dynamics with Respect to the Driver s Input Mojtaba Forghani and Domitilla Del Vecchio Massachusetts Institute of Technology September 19, Introduction In this report, we consider longitudinal dynamics of a vehicle and prove that its longitudinal position and velocity are order preserving with respect to the driver s (or automatic controller s) input. We also prove that the position is strictly order preserving with respect to the constant term of the acceleration replaced with a linear function of the position and the velocity. This estimation of the acceleration is used for the cases in which no direct measurement of the acceleration is available. 2 System Model Before introducing the model that we are considering, we define the strict and non-strict order preserving properties. Definition 1. (Order preserving) For all w, z R n we have that w z (w < z) if and only if w i z i (w i < z i ) for all i {1, 2,..., n}, in which w i denotes the ith component of w. We denote the piecewise continuous signal on U by S(U) : R + U. For U R m we define the partial order (strict partial order) by component-wise ordering for all times, that is, for all w, z S(U) we have that w z (w < z) provided w(t) z(t) (w(t) < z(t)) for all t R +. The map f : P Q is order preserving (strict order preserving) provided if for x, y P we have x y (x < y), then f(x) f(y) (f(x) < f(y)). Definition 2. (Continuous system) A continuous system is a tuple Σ = (X, U,, O, f, h), with state x X R n, control input u U R m, disturbance input d R q, output y O X, vector field in the form of f : X U X, and output map h : X O. We denote the flow of a system Σ at time t R + by φ(t, x, u, d), with initial condition x X, control input signal u S(U) and disturbance input signal d S( ). We also denote the ith component of the flow by φ i (t, x, u, d). Definition 3. (Control input/output order preserving) A continuous system Σ = (X, U,, O, f, h) is called input/output order preserving (strict input/output order preserving) with respect to the control input, if the map h(φ(t, x,, d)) : U O, for any fixed t, x and d, is order preserving (strict order preserving). 1
2 Definition 4. (Disturbance input/output order preserving) A continuous system Σ = (X, U,, O, f, h) is called input/output order preserving (strict input/output order preserving) with respect to the disturbance input, if the map h(φ(t, x, u, )) : O, for any fixed t, x and u, is order preserving (strict order preserving). We denote the position and the velocity of vehicle by p and v, respectively. We consider two systems Σ 1 and Σ 2. For Σ 1 we assume that the input is the control input u and there is no disturbance. Therefore it is reasonable to represent the flow of Σ 1 by φ 1 (t, x 1, u), where x 1 is the initial state of the system. The deceleration due to the road load (rolling resistance) and the slope of the road are represented by a r and a s, respectively, and the drag coefficient is denoted by D. We also impose a condition that the speed of the vehicle must be non-negative. We use the superscript T to denote the transpose of a vector or matrix, e.g., A T represents the transpose of matrix A. The dynamics of Σ 1 is as the following: x 1 X R 2, where x 1 = (p, v) T, (1) u U R, with U = {u u [u m, u M ]}, (2) where u m is the minimal control input, and u M is the maximal control input, where f 1 : X U X, with ẋ 1 = f 1 (x 1, u), (3) f 1 (x 1, u) = The function f 1 (x 1, u) is also in the following form: [ f 1 (x 1 v, u) = u Dv 2 a r a s { f 1 (x 1, u) if v > if v. (4) ]. (5) Since we cannot measure the acceleration for system Σ 2, we assume the total acceleration of the vehicle to be in the form of a linear function of the position and the velocity. Moreover we assume the input term (the constant term of the linear function), be a disturbance input. Therefore it is reasonable to represent the flow of Σ 2 by φ 2 (t, x 2, d), where x 2 is the initial state of the system. The dynamics of system Σ 2 is in the following form: where x 2 X R 2, where x 2 = (p, v) T, (6) d R, (7) f 2 : X X, with ẋ 2 = f 2 (x 2, d), (8) f 2 (x 2, d) = The function f 2 (x 2, d) is also in the following form: [ f 2 (x 2 v, d) = ap + bv + d { f 2 (x 2, d) if v > if v. (9) ]. (1) The term ap + bv + d is the total acceleration of the vehicle in Σ 2, where a and b are known constants. Since in Σ 2, unlike Σ 1, we cannot measure the acceleration directly, this model can be a good approximation to the acceleration of the vehicle. 2
3 3 Order Preserving Properties We prove that the flows of the position and the velocity in Σ 1, φ 1 1 (t, x1, u) and φ 1 2 (t, x1, u), respectively, have order preserving property with respect to the control input signal u. Proposition 1. The flows φ 1 1 (t, x1, u) and φ 1 2 (t, x1, u) are order preserving with respect to the control input signal u. Proof. If we consider two different input signals u 1 and u 2, such that u 1 > u 2, then for the velocity of Σ 1 at time t corresponding to these two input signals, with the same initial conditions p 1 () = p 2 () = p() and v 1 () = v 2 () = v(), we have v 1 (t) = u 1 (t) Dv 1 (t) 2 a r a s and v 2 (t) = u 2 (t) Dv2 2(t) a r a s, if both v 1 (t) > and v 2 (t) >. Let the function g(t) := v 1 (t) v 2 (t). At an arbitrary time t we have ġ(t) = v 1 (t) v 2 (t) = (u 1 (t) u 2 (t)) D ( v 2 1(t) v 2 2(t) ). (11) Note that since we have chosen the same initial conditions, we have g() = v 1 () v 2 () =. Because of the continuity of flow of the system with respect to time, if order in state v is not preserved, we must have a time t R + such that g(t ) =, since otherwise for all t R +, either g(t) < or g(t) >. Therefore we can define t := min{t R + g(t) = }. Since ġ( + ) = u 1 ( + ) u 2 ( + ) >, ġ(t ) = u 1 (t ) u 2 (t ) >, and g() = g(t ) =, for the interval t (, t ) we have and similarly ġ( + g(h) g() g(h) ) = lim = lim h + h h + h > since h > : h = h 1 (, t ) s.t. g(h 1 ) >, (12) ġ(t g(t ) g(t + h) g(t + h) ) = lim h t (t = lim > + h) h h since h < : h = h 2 (, t ) s.t. g(h 2 ) <, (13) and because of the continuity of the flow with respect to time, there is a t [h 1, h 2 ] such that g(t) =, which is in contradiction with the initial assumption that t := min{t R + g(t) = }. Therefore there is no such t, and for all t R + such that v 1 (t) > and v 2 (t) > we have either g(t) = v 1 (t) v 2 (t) > or g(t) = v 1 (t) v 2 (t) <. From (12) we conclude that the former is true. We had assumed initially that v 1 (t) > and v 2 (t) >. For a case that for some t R + we have v 1 (t ) = and v 2 (t ) =, we let t 1 := min{t R + v 1 (t) = } and t 2 = min{t R + v 2 (t) = }. Because of the non-negativity of v 1 (t), v 2 (t) and v 1 (t) v 2 (t), we must have t 2 t 1. If an arbitrary time t (, t 2 ), then g(t) = v 1 (t) v 2 (t) > ; if t [ t 2, t 1 ), then v 1 (t) v 2 (t) = v 1 (t) > ; and if t [ t 1, ), then g(t) = v 1 (t) v 2 (t) =. Therefore, in any case the order of the flow of the velocity is preserved with respect to the control input signal. Since p 1 () = p 2 () = p(), then based on equation (5), g(s)ds, which implies that the order preserving property of the flow of p is also satisfied with respect to the control input signal. In Proposition 2 we prove that the position in Σ 2, φ 2 1 (t, x2, d), is strictly order preserving with respect to d. 3
4 Proposition 2. input d. The flow φ 2 1 (t, x2, d) is strictly order preserving with respect to the disturbance Proof. Let x 2 := (p, v ) T be the initial condition, where v >. Since d is not a function of time, then by differentiating v from equations (8) and (1), we have that the velocity in Σ 2, for v(t) >, satisfies the following differential equation: v b v av = where v() = v and v() = ap + bv + d. (14) The above differential equation has the solution in the form v(t) = k 1 e λ 1t + k 2 e λ 2t where, λ 1 =.5(b + b 2 + 4a) and λ 2 =.5(b b 2 + 4a). (15) Since complex and real values of λ 1 and λ 2 reveal different behaviors for v(t), we consider different possible cases and analyze the behaviors of v(t) and p(t) with respect to d for each of them. We divide the problem into three different cases; (1): b 2 + 4a >, (2): b 2 + 4a < and (3): b 2 + 4a =. For each case we consider two input signals d 1 = d 1 and d 2 = d 2 such that d 1 > d 2 and determine the relationship between v 1 (t) and v 2 (t) and then between p 1 (t) and p 2 (t), the velocity and the position at time t corresponding to d 1 and d 2, respectively. Case (1): If b 2 + 4a >, then λ 1 and λ 2 in (15) are real numbers. The solution of (14) then takes the form 1 ( ) v(t) = (v (λ 2 b) ap d) e λ1t (v (λ 1 b) ap d) e λ 2t. (16) If we replace d in equation (16) with d 1 and d 2 in order to obtain their corresponding velocities at time t, represented by v 1 (t) and v 2 (t), respectively, we have v 1 (t) v 2 (t) = d1 d 2 ( e λ 2t e λ 1t ). (17) Equation (17) can become zero only when t =. Therefore because of the continuity of flow of the system with respect to time, for all t R +, either v 1 (t) v 2 (t) > or v 1 (t) v 2 (t) <. To determine which of these two cases holds, we note that in general for any x R {} we have that if x >, then e x 1 > and if x <, then e x 1 <. These two statements together imply that ex 1 x >. Since in Case (1) and we are considering t R +, then ( )t. Therefore we can replace x with ( )t in order to obtain e (λ 2 λ 1 )t 1 ( )t > te λ 1t [ e (λ 2 λ 1 )t 1 ( )t ] > eλ 2t e λ 1t > d 1 d 2 ( ) e λ2t e λ 1t > v 1 (t) > v 2 (t), (18) where we have used the facts that t R +, e λ1t > and d 1 d 2 >. By integrating both sides of (18) to determine the position at time t, we obtain v 1 (s)ds > 4 v 2 (s)ds
5 p + v 1 (s)ds > p + v 2 (s)ds p 1 (t) > p 2 (t). (19) Case (2): If b 2 + 4a <, then λ 1 and λ 2 in (15) are complex numbers. The solution of (14) then takes the form ( ) v(t) = e αt ap + (b α)v + d sin βt + v cos βt, β with α =.5b, and β =.5 (b 2 + 4a). (2) If we replace d in equation (2) with d 1 and d 2 in order to obtain their corresponding velocities at time t, represented by v 1 (t) and v 2 (t), respectively, we have v 1 (t) v 2 (t) = d1 d 2 e αt sin βt. (21) β We observe that in Case (2), unlike Case (1), we cannot guarantee that for all t R +, v 1 (t) v 2 (t). Note that v 1 (t) v 2 (t) = for all t such that sin βt = or alternatively, βt = kπ, for all k Z. The smallest t R + that satisfies sin βt = is t = π β. The velocity at time t corresponding to d 1 and d 2, based on equation (2), is given by v i (t ) = v e αt cos βπ β = v e αt < for i {1, 2}. (22) Since for all t R +, v(t), we must have v i (t ) =, or in other words, for all t [, t ] we have either v 1 (t) v 2 (t) or v 1 (t) v 2 (t). In order to determine which case holds, we note that for all t [, t ] we have β > and sin βt 1. Therefore in any case, for all t [, t ] we have. Also e αt (d 1 d 2 ) >. These two statements along with (21) imply sin βt β that v 1 (t) v 2 (t). Since in (2) we have v 2 () = v > and v 2 (t ) = v e αt <, then because of the continuity of flow of the system with respect to time, there is a t (, t ) such that t := min{t (, t ) v 2 (t) = }. Then we have for all t (, t), v 1 (t) v 2 (t) >. For a t (, t), we have for a t [ t, t ) we have + and for a t [t, ), we have t t (v 1 (s) v 2 (s))ds + (v 1 (s) v 2 (s))ds > ; (23) t (v 1 (s) v 2 (s))ds > (v 1 (s) v 2 (s))ds p 1 (t) p 2 (t) > ; (24) (v 1 (s) v 2 (s))ds + (v 1 (s) v 2 (s))ds = t (v 1 (s) v 2 (s))ds + > p 1 (t) p 2 (t) >. (25) 5
6 Case (3): If b 2 + 4a =, then λ 1 = λ 2 = λ, which is also a real number. The solution of (14) then takes the form v(t) = e λt [v + (ap + (b λ)v + d) t], (26) and for v 1 (t) v 2 (t) we have which implies v 1 (t) v 2 (t) = t(d 1 d 2 )e λt >, (27) (v 1 (s) v 2 (s))ds >. (28) We had assumed initially that v 1 (t) > and v 2 (t) >. In general, we may have a time t such that v 1 (t ) = and v 2 (t ) =. In this case, because of the continuity of flow of the system with respect to time, there are times t 1 and t 2 such that t 1 = sup{t (, t ) v 1 (t) > } and t 2 = sup{t (, t ) v 2 (t) > }. Since we have proved through Cases (1)-(3) that as long as v 1 (t) > and v 2 (t) > we have v 1 (t) v 2 (t) >, then t 2 t 1. For an arbitrary time τ (, t 2 ) Cases (1)-(3) imply that p 1 (τ) p 2 (τ) > ; if τ [ t 2, t 1 ), then and if τ [ t 1, ), then p 1 (τ) p 2 (τ) = and the proof is complete. p 1 (τ) p 2 (τ) = t 2 t 2 τ (v 1 (s) v 2 (s))ds + (v 1 (s) )ds > ; t 2 (29) t 1 τ (v 1 (s) v 2 (s))ds + t 2 (v 1 (s) )ds + ( )ds >, t 1 (3) 6
Design of Driver-Assist Systems Under Probabilistic Safety Specifications Near Stop Signs
Design of Driver-Assist Systems Under Probabilistic Safety Specifications Near Stop Signs The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationPUTNAM PROBLEMS DIFFERENTIAL EQUATIONS. First Order Equations. p(x)dx)) = q(x) exp(
PUTNAM PROBLEMS DIFFERENTIAL EQUATIONS First Order Equations 1. Linear y + p(x)y = q(x) Muliply through by the integrating factor exp( p(x)) to obtain (y exp( p(x))) = q(x) exp( p(x)). 2. Separation of
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationSection 4.2: The Mean Value Theorem
Section 4.2: The Mean Value Theorem Before we continue with the problem of describing graphs using calculus we shall briefly pause to examine some interesting applications of the derivative. In previous
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More informationWhat will you learn?
Section 2.2 Basic Differentiation Rules & Rates of Change Calc What will you learn? Find the derivative using the Constant Rule Find the derivative using the Power Rule Find the derivative using the Constant
More informationFormal Groups. Niki Myrto Mavraki
Formal Groups Niki Myrto Mavraki Contents 1. Introduction 1 2. Some preliminaries 2 3. Formal Groups (1 dimensional) 2 4. Groups associated to formal groups 9 5. The Invariant Differential 11 6. The Formal
More informationA Separation Principle for a Class of Hybrid Automata on a Partial Order
2009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 10-12, 2009 ThC12.4 A Separation Principle for a Class of Hybrid Automata on a Partial Order Domitilla Del Vecchio, Michael
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationSecond Order Differential Equations Lecture 6
Second Order Differential Equations Lecture 6 Dibyajyoti Deb 6.1. Outline of Lecture Repeated Roots; Reduction of Order Nonhomogeneous Equations; Method of Undetermined Coefficients Variation of Parameters
More informationStudent Study Session. Theorems
Students should be able to apply and have a geometric understanding of the following: Intermediate Value Theorem Mean Value Theorem for derivatives Extreme Value Theorem Name Formal Statement Restatement
More informationLecture 5: Lyapunov Functions and Storage Functions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.243j (Fall 2003): DYNAMICS OF NONLINEAR SYSTEMS by A. Megretski Lecture 5: Lyapunov Functions and Storage
More informationLecture Notes in Mathematics. Arkansas Tech University Department of Mathematics
Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics Introductory Notes in Ordinary Differential Equations for Physical Sciences and Engineering Marcel B. Finan c All Rights
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationAP Physics Free Response Practice Kinematics ANSWERS 1982B1 2
AP Physics Free Response Practice Kinematics ANSWERS 198B1 a. For the first seconds, while acceleration is constant, d = ½ at Substituting the given values d = 10 meters, t = seconds gives a = 5 m/s b.
More informationMath 106 Answers to Exam 3a Fall 2015
Math 6 Answers to Exam 3a Fall 5.. Consider the curve given parametrically by x(t) = cos(t), y(t) = (t 3 ) 3, for t from π to π. (a) (6 points) Find all the points (x, y) where the graph has either a vertical
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationPredicting the future with Newton s Second Law
Predicting the future with Newton s Second Law To represent the motion of an object (ignoring rotations for now), we need three functions x(t), y(t), and z(t), which describe the spatial coordinates of
More informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationObserver design for a general class of triangular systems
1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals
More informationL2 gains and system approximation quality 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility
More informationINHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS
INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS Definitions and a general fact If A is an n n matrix and f(t) is some given vector function, then the system of differential equations () x (t) Ax(t)
More informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More informationMA5510 Ordinary Differential Equation, Fall, 2014
MA551 Ordinary Differential Equation, Fall, 214 9/3/214 Linear systems of ordinary differential equations: Consider the first order linear equation for some constant a. The solution is given by such that
More informationLecture 2. 1D motion with Constant Acceleration. Vertical Motion.
Lecture 2 1D motion with Constant Acceleration. Vertical Motion. Types of motion Trajectory is the line drawn to track the position of an abject in coordinates space (no time axis). y 1D motion: Trajectory
More informationOHSx XM521 Multivariable Differential Calculus: Homework Solutions 14.1
OHSx XM5 Multivariable Differential Calculus: Homework Solutions 4. (8) Describe the graph of the equation. r = i + tj + (t )k. Solution: Let y(t) = t, so that z(t) = t = y. In the yz-plane, this is just
More informationTheory of Ordinary Differential Equations
Theory of Ordinary Differential Equations Existence, Uniqueness and Stability Jishan Hu and Wei-Ping Li Department of Mathematics The Hong Kong University of Science and Technology ii Copyright c 24 by
More informationCHAPTER 7. Connectedness
CHAPTER 7 Connectedness 7.1. Connected topological spaces Definition 7.1. A topological space (X, T X ) is said to be connected if there is no continuous surjection f : X {0, 1} where the two point set
More informationNonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems
Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationDifferential Equations: Homework 2
Differential Equations: Homework Alvin Lin January 08 - May 08 Section.3 Exercise The direction field for provided x 0. dx = 4x y is shown. Verify that the straight lines y = ±x are solution curves, y
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationCS Tutorial 5 - Differential Geometry I - Surfaces
236861 Numerical Geometry of Images Tutorial 5 Differential Geometry II Surfaces c 2009 Parameterized surfaces A parameterized surface X : U R 2 R 3 a differentiable map 1 X from an open set U R 2 to R
More information2 Statement of the problem and assumptions
Mathematical Notes, 25, vol. 78, no. 4, pp. 466 48. Existence Theorem for Optimal Control Problems on an Infinite Time Interval A.V. Dmitruk and N.V. Kuz kina We consider an optimal control problem on
More informationFind the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x
Assignment 5 Name Find the indicated derivative. ) Find y(4) if y = sin x. ) A) y(4) = cos x B) y(4) = sin x y(4) = - cos x y(4) = - sin x ) y = (csc x + cot x)(csc x - cot x) ) A) y = 0 B) y = y = - csc
More informationSection , #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation.
Section.3.5.3, #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation dq = 1 4 (1 + sin(t) ) + Q, Q(0) = 50. (1) 100 (a) The differential equation given
More informationParticle Motion Problems
Particle Motion Problems Particle motion problems deal with particles that are moving along the x or y axis. Thus, we are speaking of horizontal or vertical movement. The position, velocity, or acceleration
More informationINTRODUCTION TO ALGEBRAIC TOPOLOGY. (1) Let k < j 1 and 0 j n, where 1 n. We want to show that e j n e k n 1 = e k n e j 1
INTRODUCTION TO ALGEBRAIC TOPOLOGY Exercise 7, solutions 1) Let k < j 1 0 j n, where 1 n. We want to show that e j n e k n 1 = e k n e j 1 n 1. Recall that the map e j n : n 1 n is defined by e j nt 0,...,
More informationMA Ordinary Differential Equations and Control Semester 1 Lecture Notes on ODEs and Laplace Transform (Parts I and II)
MA222 - Ordinary Differential Equations and Control Semester Lecture Notes on ODEs and Laplace Transform (Parts I and II Dr J D Evans October 22 Special thanks (in alphabetical order are due to Kenneth
More informationNotes: DERIVATIVES. Velocity and Other Rates of Change
Notes: DERIVATIVES Velocity and Oter Rates of Cange I. Average Rate of Cange A.) Def.- Te average rate of cange of f(x) on te interval [a, b] is f( b) f( a) b a secant ( ) ( ) m troug a, f ( a ) and b,
More informationk=1 = e λ λ k λk 1 k! (k + 1)! = λ. k=0
11. Poisson process. Poissonian White Boise. Telegraphic signal 11.1. Poisson process. First, recall a definition of Poisson random variable π. A random random variable π, valued in the set {, 1,, }, is
More information( ) 4 and 20, find the value. v c is equal to this average CALCULUS WORKSHEET 1 ON PARTICLE MOTION
CALCULUS WORKSHEET 1 ON PARTICLE MOTION Work these on notebook paper. Use your calculator only on part (f) of problems 1. Do not use your calculator on the other problems. Write your justifications in
More informationTeam-Exercises for DGC 100 Modelica Course
Team-Exercises for DGC 100 Modelica Course Hubertus Tummescheit United Technologies Research Center, East Hartford, CT 06108. November 4, 2003 Abstract This document is a preliminary version and is going
More informationf(s)ds, i.e. to write down the general solution in
1. First order ODEs. Simplest examples. Existence and uniqueness theorem Example 1.1. Consider the equation (1.1) x (t) = 1 This is an equation because there is an unknown: the function x(t). This is a
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationPassivity-based Stabilization of Non-Compact Sets
Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained
More informationMATH 308 Differential Equations
MATH 308 Differential Equations Summer, 2014, SET 6 JoungDong Kim Set 6: Section 3.3, 3.4, 3.5, 3.6 Section 3.3 Complex Roots of the Characteristic Equation Recall that a second order ODE with constant
More information4.1 Areas and Distances. The Area Problem: GOAL: Find the area of the region S that lies under the curve y = f(x) from a to b.
4.1 Areas and Distances The Area Problem: GOAL: Find the area of the region S that lies under the curve y = f(x) from a to b. Easier Problems: Find the area of a rectangle with length l and width w. Find
More informationAnswer Key. Calculus I Math 141 Fall 2003 Professor Ben Richert. Exam 2
Answer Key Calculus I Math 141 Fall 2003 Professor Ben Richert Exam 2 November 18, 2003 Please do all your work in this booklet and show all the steps. Calculators and note-cards are not allowed. Problem
More informationStatic Problem Set 2 Solutions
Static Problem Set Solutions Jonathan Kreamer July, 0 Question (i) Let g, h be two concave functions. Is f = g + h a concave function? Prove it. Yes. Proof: Consider any two points x, x and α [0, ]. Let
More informationCalculus of Variations. Final Examination
Université Paris-Saclay M AMS and Optimization January 18th, 018 Calculus of Variations Final Examination Duration : 3h ; all kind of paper documents (notes, books...) are authorized. The total score of
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
More informationNonlinear Systems and Elements of Control
Nonlinear Systems and Elements of Control Rafal Goebel April 3, 3 Brief review of basic differential equations A differential equation is an equation involving a function and one or more of the function
More informationUNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems
UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems Robert M. Freund February 2016 c 2016 Massachusetts Institute of Technology. All rights reserved. 1 1 Introduction
More informationRoot Locus Design Example #3
Root Locus Design Example #3 A. Introduction The system represents a linear model for vertical motion of an underwater vehicle at zero forward speed. The vehicle is assumed to have zero pitch and roll
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationKatznelson Problems. Prakash Balachandran Duke University. June 19, 2009
Katznelson Problems Prakash Balachandran Duke University June 9, 9 Chapter. Compute the Fourier coefficients of the following functions (defined by their values on [ π, π)): f(t) { t < t π (t) g(t) { t
More informationRelative Controllability of Fractional Dynamical Systems with Multiple Delays in Control
Chapter 4 Relative Controllability of Fractional Dynamical Systems with Multiple Delays in Control 4.1 Introduction A mathematical model for the dynamical systems with delayed controls denoting time varying
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationA note on linear differential equations with periodic coefficients.
A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.
More informationUnit #5 Applications of the Derivative Part II Homework Packet
Unit #5 Applications of the Derivative Part II Homework Packet 1. For which of the following functions is the Extreme Value Theorem NOT APPLICABLE on the interval [a, b]? Give a reason for your answer.
More informationLinear Algebra Lecture Notes-I
Linear Algebra Lecture Notes-I Vikas Bist Department of Mathematics Panjab University, Chandigarh-6004 email: bistvikas@gmail.com Last revised on February 9, 208 This text is based on the lectures delivered
More informationLeaving Cert Differentiation
Leaving Cert Differentiation Types of Differentiation 1. From First Principles 2. Using the Rules From First Principles You will be told when to use this, the question will say differentiate with respect
More informationComments on integral variants of ISS 1
Systems & Control Letters 34 (1998) 93 1 Comments on integral variants of ISS 1 Eduardo D. Sontag Department of Mathematics, Rutgers University, Piscataway, NJ 8854-819, USA Received 2 June 1997; received
More informationMath : Solutions to Assignment 10
Math -3: Solutions to Assignment. There are two tanks. The first tank initially has gallons of pure water. The second tank initially has 8 gallons of a water/salt solution with oz of salt. Both tanks drain
More informationAverage and Instantaneous Velocity. p(a) p(b) Average Velocity on a < t < b =, where p(t) is the position a b
Particle Motion Problems Particle motion problems deal with particles that are moving along the x or y axis. Thus, we are speaking of horizontal of vertical movement. The position, velocity or acceleration
More informationRiesz Representation Theorems
Chapter 6 Riesz Representation Theorems 6.1 Dual Spaces Definition 6.1.1. Let V and W be vector spaces over R. We let L(V, W ) = {T : V W T is linear}. The space L(V, R) is denoted by V and elements of
More informationXinfu Chen. Topics in Differential Equations. Department of mathematics, University of Pittsburgh Pittsburgh, PA 15260, USA
Xinfu Chen Topics in Differential Equations Department of mathematics, University of Pittsburgh Pittsburgh, PA 1526, USA June 26, 26 ii Preface This notes is for math 3926, topics in differential equations,
More informationFormulas that must be memorized:
Formulas that must be memorized: Position, Velocity, Acceleration Speed is increasing when v(t) and a(t) have the same signs. Speed is decreasing when v(t) and a(t) have different signs. Section I: Limits
More information(4.1) D r v(t) ω(t, v(t))
1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution
More informationLecture 2. 1D motion with Constant Acceleration. Vertical Motion.
Lecture 2 1D motion with Constant Acceleration. Vertical Motion. Types of motion Trajectory is the line drawn to track the position of an abject in coordinates space (no time axis). y 1D motion: Trajectory
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationProblem 1: Ship Path-Following Control System (35%)
Problem 1: Ship Path-Following Control System (35%) Consider the kinematic equations: Figure 1: NTNU s research vessel, R/V Gunnerus, and Nomoto model: T ṙ + r = Kδ (1) with T = 22.0 s and K = 0.1 s 1.
More information1 y 2 dy = (2 + t)dt 1 y. = 2t + t2 2 + C 1 2t + t 2 /2 + C. 1 t 2 /2 + 2t 1. e y y = 2t 2. e y dy = 2t 2 dt. e y = 2 3 t3 + C. y = ln( 2 3 t3 + C).
Math 53 First Midterm Page. Solve each of the following initial value problems. (a) y = y + ty, y() = 3. The equation is separable : y = y ( + t). Thus y = y dy = ( + t)dt y = t + t + C t + t / + C. For
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationMath 1310 Lab 10. (Sections )
Math 131 Lab 1. (Sections 5.1-5.3) Name/Unid: Lab section: 1. (Properties of the integral) Use the properties of the integral in section 5.2 for answering the following question. (a) Knowing that 2 2f(x)
More informationMA304 Differential Geometry
MA304 Differential Geoetry Hoework 4 solutions Spring 018 6% of the final ark 1. The paraeterised curve αt = t cosh t for t R is called the catenary. Find the curvature of αt. Solution. Fro hoework question
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More information2nd-Order Linear Equations
4 2nd-Order Linear Equations 4.1 Linear Independence of Functions In linear algebra the notion of linear independence arises frequently in the context of vector spaces. If V is a vector space over the
More informationPHYS1100 Practice problem set, Chapter 2: 6, 10, 13, 17, 22, 24, 25, 34, 42, 50, 55, 65, 71, 82
PHYS00 Practice problem set, Chapter : 6, 0, 3, 7,, 4, 5, 34, 4, 50, 55, 65, 7, 8.6. Model: We will consider Larry to be a particle. Solve: Since Larry s speed is constant, we can use the following equation
More informationNONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality
M athematical Inequalities & Applications [2407] First Galley Proofs NONLINEAR DIFFERENTIAL INEQUALITY N. S. HOANG AND A. G. RAMM Abstract. A nonlinear differential inequality is formulated in the paper.
More informationRecitation 1. Gradients and Directional Derivatives. Brett Bernstein. CDS at NYU. January 21, 2018
Gradients and Directional Derivatives Brett Bernstein CDS at NYU January 21, 2018 Brett Bernstein (CDS at NYU) Recitation 1 January 21, 2018 1 / 23 Initial Question Intro Question Question We are given
More informationGENERALIZATION OF GRONWALL S INEQUALITY AND ITS APPLICATIONS IN FUNCTIONAL DIFFERENTIAL EQUATIONS
Communications in Applied Analysis 19 (215), 679 688 GENERALIZATION OF GRONWALL S INEQUALITY AND ITS APPLICATIONS IN FUNCTIONAL DIFFERENTIAL EQUATIONS TINGXIU WANG Department of Mathematics, Texas A&M
More informationHOMEWORK FOR , FALL 2007 ASSIGNMENT 1 SOLUTIONS. T = sup T(x) Ax A := sup
HOMEWORK FOR 18.101, FALL 007 ASSIGNMENT 1 SOLUTIONS (1) Given a linear map T : R m R n Define the operator norm of T as follows: T(x) T := sup x 0 x Similarly, if A is a matrix, define the operator norm
More information2019 Spring MATH2060A Mathematical Analysis II 1
2019 Spring MATH2060A Mathematical Analysis II 1 Notes 1. CONVEX FUNCTIONS First we define what a convex function is. Let f be a function on an interval I. For x < y in I, the straight line connecting
More information2.4 Differences Between Linear and Nonlinear Equations 75
.4 Differences Between Linear and Nonlinear Equations 75 fying regions of the ty-plane where solutions exhibit interesting features that merit more detailed analytical or numerical investigation. Graphical
More informationarxiv: v1 [math.na] 28 Jan 2009
The Dynamical Systems Method for solving nonlinear equations with monotone operators arxiv:0901.4377v1 [math.na] 28 Jan 2009 N. S. Hoang and A. G. Ramm Mathematics Department, Kansas State University,
More informationSafety control of piece-wise continuous order preserving systems
Safety control of piece-wise continuous order preserving systems Reza Ghaemi and Domitilla Del Vecchio Abstract This paper is concerned with safety control of systems with imperfect state information and
More informationChapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. D) u = - x15 cos (x15) + C
AP Calculus AB Exam Review Differential Equations and Mathematical Modelling MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the general solution
More informationDynamic interpretation of eigenvectors
EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation
More information4.3 Linear, Homogeneous Equations with Constant Coefficients. Jiwen He
4.3 Exercises Math 3331 Differential Equations 4.3 Linear, Homogeneous Equations with Constant Coefficients Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math3331
More information