Calculus of Variations
|
|
- Melvyn Mills
- 5 years ago
- Views:
Transcription
1 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t Figure by MIT OCW.
2 Spr Calculus of Variations Goal: Develop alternative approach to solve general optimization problems for continuous systems variational calculus Formal approach will provie new insights for constraine solutions, an a more irect path to the solution for other problems. Main issue Cost metric is a function of the functions x(t) an u(t). min J = h(x(t f )) + g(x(t), u(t), t)) subject to ẋ = f (x, u, t) x( ), = given m(x(t f ), t f ) = 0 Call J (x(t), u(t)) a functional. Nee to investigate how to fin the optimal values of a functional. For a function, we foun the graient, an set it to zero to fin the stationary points, an then investigate the higher orer erivatives to etermine if it is a maximum or minimum. Will investigate something similar for functionals.
3 Spr Maximum an minimum: function A function f (x) has a local minimum at x if f (x) f (x ) for all amissible x in x x ɛ Minimum can occur at (i) stationary point, (ii) at a bounary, or (iii) a point of iscontinuous erivative. If we only consier stationary points of the ifferentiable function f (x), then this conition is equivalent to requiring that the ifferential of f satisfy: f f = x = 0 x for all small x, which gives the same necessary conition as before f = 0 x Note that this efinition use norms to compare two vectors. Can o the same thing with functions: istance between two functions where = x 2 (t) x 1 (t) 1. x(t) 0 for all x(t), an x(t) = 0 only if x(t) = 0 for all t in the interval of efinition. 2. ax(t) = a x(t) for all real scalars a. 3. x 1 (t) + x 2 ( t) x 1 (t) + x 2 (t) Common functional norm: x(t) 2 = ( ) 1/2 x(t) T x(t)
4 Spr Maximum an minimum: functional A functional J(x(t)) has a local minimum at x (t) if J(x(t)) J(x (t)) for all amissible x(t) in x(t) x (t) ɛ Now efine something equivalent to the ifferential of a function calle a variation of a functional. An increment of a functional ΔJ(x(t), δx(t)) = J(x(t) + δx(t)) J(x(t)) A variation of the functional is a linear approximation of this increment: ΔJ(x(t), δx(t)) = δj(x(t), δx(t)) + H.O.T. i.e. δj(x(t), δx(t)) is linear in δx(t). f(t) Slope equals f' (t 1 ) f f t t t 1 Differential f versus increment f shown for a function, but the same ifference hols for a functional. Figure by MIT OCW.
5 Spr x(t) x* x*+ αδx (1) x*- αδx(1) αδx (1) αδx (1) t f t Figure by MIT OCW. Figure 2: Visualization of perturbations to function x(t) by δx(t) it is a potential change in the value of x over the entire time perio of interest. Typically require that if x(t) is in some class (i.e., continuous), that x(t) + δx(t) is also in that class. Funamental theorem of the calculus of variations Let x be a function of t in the class Ω, an J(x) be a ifferentiable functional of x. Assume that the functions in Ω are not constraine by any bounaries. If x is an extremal function, then the variation of J must vanish on x, i.e. for all amissible δx, δj(x(t), δx(t)) = 0 Proof is in Kirk, page 121, but it is relatively straightforwar. How o we compute the variation? If J(x(t)) = f(x(t)) where f has cts first an secon erivatives with respect to x, then { } f(x(t)) δj(x(t), δx) = (t) δx + f(x(t f ))δt f f(x( ))δ x tf = f x (x(t))δx + f(x(t f ))δt f f(x( ))δ
6 Spr 2006 Variation Examples: Scalar For more general problems, consier the cost evaluate on a scalar function x(t) J(x(t)) = tf g(x(t), ẋ(t), t) with, t f an the curve enpoints fixe. Then δj(x(t), δx) = [ g x (x(t), ẋ(t), t)δx + g ẋ (x(t), ẋ(t), t)δẋ] Note that δẋ = δx so δx an δẋ are not inepenent. Integrate by parts: uv = uv vu with u = g ẋ an v = δ ẋ to get: δj(x(t), δx) t = f gx (x(t), ẋ(t), t)δx + [g ẋ (x(t), ẋ(t), t)δx] t0 gẋ(x(t), ẋ(t), t)δx t 0 t f = g x (x(t), ẋ(t), t) g ẋ (x(t), ẋ(t), t) δx(t) t + f [gẋ (x(t), ẋ(t), t)δx] t0 But since x( ) an x(t f ) are given, then δx( ) = δx(t f ) = 0, yieling δj(x(t), δx) = g x (x(t), ẋ(t), t) g ẋ (x(t), ẋ(t), t) δx(t)
7 Spr Recall that we nee δj = 0 for all amissible δx(t), which are arbitrary within (, t f ), then the (first orer) necessary conition for a maximum or minimum is that ( ) g(x(t), ẋ(t), t) g(x(t), ẋ(t), t) = 0 x ẋ which is calle the Euler equation. Example: Fin the curve that gives the shortest istance between 2 points in a plane (x 0, y 0 ) an (x f, y f ). Cost function sum of ifferential arc lengths: x f x f J = s = (x) 2 + (y) 2 x 0 x 0 x f ( ) 2 y = 1 + x x 0 x Take y as the epenent variable, an x as the inepenent one y x ẏ New form of the cost: x f x f J = 1 + ẏ 2 x g(ẏ)x x 0 x 0 Take partials: g/ y = 0, an ( ) ( ) ( ) g g ẏ ẏ ÿ = = ÿ = = 0 x ẏ ẏ ẏ x ẏ (1 + y 2 ) 1/2 (1 + y 2 ) 3/2 which implies that ÿ = 0 Most general curve with ÿ = 0 is a line y = c 1 x + c 2
8 Spr 2006 Vector Functions Can generalize the problem by incluing several (N) functions x i (t) an possibly free enpoints J (x(t)) = g(x(t), ẋ(t), t) with, t f, x( ) fixe. Then (rop the arguments for brevity) δj (x(t), δx) = [ g x δx(t) + g ẋ δẋ(t)] Integrate by parts to get: δj (x(t), δx) = g x g ẋ δx(t) + g ẋ (x(t f ), ẋ(t f ), t f )δx(t f ) The requirement then is that for t (, t f ), x(t) must satisfy g g = 0 x ẋ where x( ) = x 0 which are the given N bounary conitions, an the remaining N more BC follow from: x(t f ) = x f if x f is given as fixe, If x(t f ) are free, then g(x(t), ẋ(t), t) = 0 ẋ(t f ) Note that we coul also have a mixture, where parts of x(t f ) are given as fixe, an other parts are free just use the rules above on each component of x i (t f )
9 Spr 2006 Free Terminal Time Now consier a slight variation: the goal is to minimize J (x(t)) = g(x(t), ẋ(t), t) with, x( ) fixe, t f free, an various constraints on x(t f ) Compute the variation of the functional where we consier 2 caniate solutions: x(t), which we consier to be a perturbation of the optimal x (t) (that we nee to fin) δj (x (t), δx) = [ g x δx(t) + g ẋ δẋ(t)] + g(x (t f ), ẋ (t f ), t f )δt f Integrate by parts to get: δj (x t f (t), δx) = g x g ẋ δx(t) + g ẋ (x (t f ), ẋ (t f ), t f )δx(t f ) + g(x (t f ), ẋ (t f ), t f )δt f Looks stanar so, but we have to be careful how we hanle the terminal conitions
10 Spr x (t) δx(t f ) x f x x* δx f x 0 t f + δt f Comparison of possible changes to function at en time when t f is free. t f t By efinition, δx(t f ) is the ifference between two amissible functions at time t f (in this case the optimal solution x an another caniate x. But in this case, we must also account for possible changes to δt f. Define δx f as being the ifference between the ens of the two possible functions total possible change in the final state: δx f δx(t f ) + ẋ (t f )δt f so δx(t f ) = δx f in general. Figure by MIT OCW. Substitute to get δj (x t f (t), δx) = g x g ẋ δx(t) + g ẋ (x (t f ), ẋ (t f ), t f )δx f + [g(x (t f ), ẋ (t f ), t f ) g ẋ (x (t f ), ẋ (t f ), t f )ẋ (t f )] δt f
11 Spr Inepenent of the terminal constraint, the conitions on the solution x (t) to be an extremal for this case are that it satisfy the Euler equations gẋ (x (t), g x (x (t), ẋ (t), t) x (t), t) = 0 Now consier the aitional constraints on the iniviual elements of x (t f ) an t f to fin the other bounary conitions Type of terminal constraints etermines how we treat δx f an δt f Unrelate Relate by a simple function x(t f ) = Θ(t f ) Specifie by a more complex constraint m(x(t f ), t f ) = 0 If t f an x(t f ) are free but unrelate, then δx f an δt f are inepenent an arbitrary their coefficients must both be zero. g x (x (t), ẋ (t), t) g ẋ(x (t), ẋ (t), t) = 0 g(x (t f ), ẋ (t f ), t f ) g ẋ(x (t f ), ẋ (t f ), t f ) ẋ (t f ) = 0 g ẋ(x (t f ), ẋ (t f ), t f ) = 0 Which makes it clear that this is a two point bounary value problem, as we have conitions at both an t f
12 Spr If t f an x(t f ) are free but relate as x(t f ) = Θ(t f ), then Θ δx f = (t f )δt f Substitute an collect terms gives [ Θ δj = g g ẋ (x (t f ), x g ẋ δx + x (t f ), t f ) (t f ) ] + g(x (t f ), ẋ (t f ), t f ) g ẋ (x (t f ), ẋ (t f ), t f )ẋ (t f ) δt f Set coefficient of δt f to zero (it is arbitrary) full conitions g ẋ(x (t f ), ẋ (t f ), t f ) g x (x (t), ẋ (t), t) g ẋ(x (t), ẋ (t), t) = 0 Θ (t f ) ẋ (t f ) + g(x (t f ), ẋ (t f ), t f ) = 0 Last conition calle the transversality conition
13 DETERMINATION OF BOUNDARY-VALUE RELATIONSHIPS Problem Description Substitution Bounary Conitions Remarks 1 x(t f ), t f both specifie (Problem 1) δx f = δx(t f ) = 0 δt f = 0 x*( ) = x 0 x*(t f ) = x f 2n equations to etermine 2n constants of integration 2 x(t f ) free; t f specifie (Problem 2) δx f = δx(t f ) δt f = 0 x*( ) = x 0 g x (x*(t f ), x*(t.. f ),t f ) = 0 2n equations to etermine 2n constants of integration 3 t f free; x(t f ) specifie (Problem 3) δx f = 0 x*( ) = x 0 x*(t f ) = x f. g(x*(t f ), x*(t f ),t f ) Ṭ g x (x*(t.. f), x*(t f ),t f ) x*(t f ) = 0 (2n + 1) equations to etermine 2n constants of integration an t f 4 t f, x(t f ) free an inepenent (Problem 4) x*( ) = x 0 g x (x*(t.. f), x*(t f ),t f ) = 0. g(x*(t f ), x*(t f ),t f ) = 0 (2n + 1) equations to etermine 2n constants of integration an t f 5 t f, x(t f ) free but relate by x(t f ) = θ(t f ) (Problem 4) θ δx f = δt (t f) f x*( ) = x 0 x*(t f ) = θ(t f ). g(x*(t f ), x*(t f ),t f ) T g x (x*(t.. f), x*(t f ),t f ) θ (t f) - x*(t f ) = 0 (2n + 1) equations to etermine 2n constants of integration an t f θ θ 1 θ enotes the n x 1 column vector 2... θ 7 T Figure by MIT OCW. Aapte from: figure 4 in Kirk, Donal E. Optimal Control Theory: An Introuction. New York, NY: Dover, ISBN: Figure 4: Possible terminal constraints To hanle the thir type of terminal conition, we must aress the solution of constraine problems.
14 Spr 2006 Example: Fin the shortest curve from the origin to a specifie line. Goal: minimize the cost functional (See page 5 6) J = 1 + ẋ 2 (t) Since g(x, ẋ, t) is only a function of ẋ, the Euler equation reuces to ẋ (t) = 0 [1 + ẋ (t) 2 ] 1/2 given that = 0, x(0) = 0, an t f an x(t f ) are free, but x(t f ) must line on the line θ(t) = 5t + 15 which after ifferentiating an simplifying, gives ẍ (t) = 0 from which we know that the answer is a straight line but since x(0) = 0, then c 0 = 0 x (t) = c 1 t + c 0 The transversality conition gives ẋ (t f ) 2 [ 5 ẋ (t f )] + [1 + ẋ (t f ) ] 1/2 = 0 [1 + ẋ (t f ) 2 ] 1/2 that simplifies to [ẋ (t f )] [ 5 ẋ (t f )] + [1 + ẋ (t f ) 2 ] = 5ẋ (t f ) + 1 = 0 so that ẋ (t f ) = c 1 = 1/5 Not a surprise, as this gives the slope of a line orthogonal to the constraint line. To fin final time: x(t f ) = 5t f + 15 = t f /5 which gives t f 2.88
15 Spr 2006 Example: Ha the terminal constraint been a bit more challenging, such as 1 Θ Θ(t) = ([t 5] 2 1) = t 5 2 Then the transversality conition gives ẋ (t f ) [t f 5 ẋ (t f )] + [1 + ẋ (t f ) 2 ] 1/2 = 0 [1 + ẋ (t f ) 2 ] 1/2 [ẋ (t f )] [t f 5 ẋ (t f )] + [1 + ẋ (t f ) 2 ] = 0 c 1 [t f 5] + 1 = 0 Now look at x (t) an Θ(t) at t f t f 1 x (t = ([t f 5] 2 f ) = 1) (tf 5) 2 which gives t f = 3, c 1 = 1/2 an x (t f ) = t/ x t Figure 5: Quaratic terminal constraint.
16 Spr 2006 Constraine Solutions Now consier variations of the basic problem that inclue constraints. For example, if the goal is to fin the extremal function x that minimizes tf J(x(t), t) = g(x(t), ẋ(t), t) subject to the constraint that a given set of n ifferential equations be satisfie f (x(t), ẋ(t), t) = 0 where we assume that x R n+m (take t f an x(t f ) to be fixe) As with the basic optimization problems in Lecture 2, we procee by augmenting the cost with the constraints using Lagrange multipliers Since the constraints must be satisfie at all time, these multipliers are also assume to be functions of time. { } J a (x(t), t) = g(x, ẋ, t) + p(t) T f (x, ẋ, t) Which oes not change the cost if the constraints are satisfie. Take variation of augmente functional consiering perturbations to both x(t) an p(t) δj(x(t), δx(t), p(t), δp(t)) { } = g x + p(t) T f x δx(t) + g ẋ + p(t) T f ẋ δẋ(t) + f T δp(t) As before, integrate by parts to get: δj(x(t), δx(t), p(t), δp(t)) ({ } ) = g x + p(t) T f x g ẋ + p(t) T f ẋ δx(t) + f T δp(t)
17 Spr To simplify things a bit, efine g a (x(t), ẋ(t), t) g(x(t), ẋ(t), t) + p(t) T f (x(t), ẋ(t), t) On the extremal, the variation must be zero, but since δx(t) an δp(t) can be arbitrary, can only occur if g a (x(t), ẋ(t), t) x g a (x(t), ẋ(t), t) ẋ = 0 f (x(t), ẋ(t), t) = 0 which are obviously a generalize version of the Euler equations that we saw before.
18 Spr General Terminal Conitions Now re consier the case on 5 10 with very general terminal constraints, with t f free, an x(t f ) given by the conition: m(x(t f ), t f ) = 0 Constraine optimization, so as before, augment the cost functional J (x(t), t) = h(x(t f ), t f ) + g(x(t), ẋ(t), t) with the constraint using Lagrange multipliers: J a (x(t), ν, t) = h(x(t f ), t f )+ν T m(x(t f ), t f )+ g(x(t), ẋ(t), t) Consiering changes to x(t), t f, x(t f ) an ν, the variation for J a is δj a = h x (t f )δx f + h tf δt f + m T (t f )δν + ν T m x (t f )δx f +ν T m tf (t f )δt f + [g x δx + g ẋ δẋ] + g(t f )δt f = h x (t f ) + ν T m x (t f ) δx f + h tf + ν T m tf (t f ) + g(t f ) δt f T +m (t f )δν + g x g ẋ δx + g ẋ (t f )δx(t f ) Now use that δx f = δx(t f ) + ẋ(t f )δt f as before to get δj a = h [ x (t f ) + ν T m x (t f ) + g ẋ (t f ) δx f ] + h tf + ν T m tf (t f ) + g(t f ) g ẋ (t f )ẋ(t f ) δt f + m T (t f )δν + g x g ẋ δx
19 Spr Looks like a bit of a mess, but we can clean it up a bit using to get w(x(t f ), ν, t f ) = h(x(t f ), t f ) + ν T m(x(t f ), t f ) δj a = [w [ x (t f ) + g ẋ (t f )] δx f ] + w tf + g(t f ) g ẋ (t f )ẋ(t f ) δt f + m T (t f )δν + g x g ẋ δx Given the extra egrees of freeom in the multipliers, we can treat all of the variations as inepenent all coefficients must be zero to achieve δj a = 0 So the necessary conitions are g x g ẋ = 0 (im n) (2) l x (t f ) + g ẋ (t f ) = 0 (im n) (3) l tf + g(t f ) g ẋ (t f )ẋ(t f ) = 0 (im 1) (4) with x( ) = x 0 (there are n of these) an m(x(t f ), t f ) = 0 (m of these), (3) an (4) From the solution of Euler s equation, there are 2n constants of integration to fin for x(t), we must also fin ν (m of those), an 2n + m + 1 variables to fin t f Some special cases: If t f is fixe, h = h(x(t f )), m m(x(t f )) an we lose the last conition (which is the transversality conition (4)) others remain unchange If t f is fixe, x(t f ) free, then there is no m, no ν an w reuces to h. Kirk consiers several other type of constraints.
Principles of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationLecture 27: Generalized Coordinates and Lagrange s Equations of Motion
Lecture 27: Generalize Coorinates an Lagrange s Equations of Motion Calculating T an V in terms of generalize coorinates. Example: Penulum attache to a movable support 6 Cartesian Coorinates: (X, Y, Z)
More informationExam 2 Review Solutions
Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon
More informationLecture 6: Calculus. In Song Kim. September 7, 2011
Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationMath 2163, Practice Exam II, Solution
Math 63, Practice Exam II, Solution. (a) f =< f s, f t >=< s e t, s e t >, an v v = , so D v f(, ) =< ()e, e > =< 4, 4 > = 4. (b) f =< xy 3, 3x y 4y 3 > an v =< cos π, sin π >=, so
More informationCalculus of variations - Lecture 11
Calculus of variations - Lecture 11 1 Introuction It is easiest to formulate the problem with a specific example. The classical problem of the brachistochrone (1696 Johann Bernoulli) is the search to fin
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationMath 115 Section 018 Course Note
Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of
More informationMath 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v
Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationSYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is
SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationFirst Order Linear Differential Equations
LECTURE 6 First Orer Linear Differential Equations A linear first orer orinary ifferential equation is a ifferential equation of the form ( a(xy + b(xy = c(x. Here y represents the unknown function, y
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationFall 2016: Calculus I Final
Answer the questions in the spaces provie on the question sheets. If you run out of room for an answer, continue on the back of the page. NO calculators or other electronic evices, books or notes are allowe
More informationCalculus Class Notes for the Combined Calculus and Physics Course Semester I
Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average
More informationd dx [xn ] = nx n 1. (1) dy dx = 4x4 1 = 4x 3. Theorem 1.3 (Derivative of a constant function). If f(x) = k and k is a constant, then f (x) = 0.
Calculus refresher Disclaimer: I claim no original content on this ocument, which is mostly a summary-rewrite of what any stanar college calculus book offers. (Here I ve use Calculus by Dennis Zill.) I
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationMathematical Review Problems
Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationApplications of the Wronskian to ordinary linear differential equations
Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More information1 Lecture 20: Implicit differentiation
Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationChapter 3 Definitions and Theorems
Chapter 3 Definitions an Theorems (from 3.1) Definition of Tangent Line with slope of m If f is efine on an open interval containing c an the limit Δy lim Δx 0 Δx = lim f (c + Δx) f (c) = m Δx 0 Δx exists,
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationSOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE. Bing Ye Wu
ARCHIVUM MATHEMATICUM (BRNO Tomus 46 (21, 177 184 SOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE Bing Ye Wu Abstract. In this paper we stuy the geometry of Minkowski plane an obtain some results. We focus
More informationWUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationCalculus and optimization
Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function
More informationTopic 2.3: The Geometry of Derivatives of Vector Functions
BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
More informationPure Further Mathematics 1. Revision Notes
Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationLecture 2 - First order linear PDEs and PDEs from physics
18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationIMPLICIT DIFFERENTIATION
IMPLICIT DIFFERENTIATION CALCULUS 3 INU0115/515 (MATHS 2) Dr Arian Jannetta MIMA CMath FRAS Implicit Differentiation 1/ 11 Arian Jannetta Explicit an implicit functions Explicit functions An explicit function
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationLecture 3 Notes. Dan Sheldon. September 17, 2012
Lecture 3 Notes Dan Shelon September 17, 2012 0 Errata Section 4, Equation (2): yn 2 shoul be x2 N. Fixe 9/17/12 Section 5.3, Example 3: shoul rea w 0 = 0, w 1 = 1. Fixe 9/17/12. 1 Review: Linear Regression
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationDefine each term or concept.
Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationArmenian Transformation Equations For Relativity
Armenian Transformation Equations For Relativity Robert Nazaryan an Haik Nazaryan 00th Anniversary of the Special Relativity Theory This Research one in Armenia 968-988, Translate from the Armenian Manuscript
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
More informationImplicit Differentiation. Lecture 16.
Implicit Differentiation. Lecture 16. We are use to working only with functions that are efine explicitly. That is, ones like f(x) = 5x 3 + 7x x 2 + 1 or s(t) = e t5 3, in which the function is escribe
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationMath Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis
Math 231 - Chapter 2 Essentials of Calculus by James Stewart Prepare by Jason Gais Chapter 2 - Derivatives 21 - Derivatives an Rates of Change Definition A tangent to a curve is a line that intersects
More informationKramers Relation. Douglas H. Laurence. Department of Physical Sciences, Broward College, Davie, FL 33314
Kramers Relation Douglas H. Laurence Department of Physical Sciences, Browar College, Davie, FL 333 Introuction Kramers relation, name after the Dutch physicist Hans Kramers, is a relationship between
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationA new approach to explicit MPC using self-optimizing control
28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur
More informationmodel considered before, but the prey obey logistic growth in the absence of predators. In
5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationWell-posedness of hyperbolic Initial Boundary Value Problems
Well-poseness of hyperbolic Initial Bounary Value Problems Jean-François Coulombel CNRS & Université Lille 1 Laboratoire e mathématiques Paul Painlevé Cité scientifique 59655 VILLENEUVE D ASCQ CEDEX, France
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More information