Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas

Similar documents
Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton Jacobi equations

Math Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Chapter 2. First Order Scalar Equations

Math 334 Fall 2011 Homework 11 Solutions

Homogenization of random Hamilton Jacobi Bellman Equations

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Second Order Linear Differential Equations

Matlab and Python programming: how to get started

MA 366 Review - Test # 1

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

An Introduction to Malliavin calculus and its applications

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Chapter 6. Systems of First Order Linear Differential Equations

Announcements: Warm-up Exercise:

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Differential Equations

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Lecture 20: Riccati Equations and Least Squares Feedback Control

Predator - Prey Model Trajectories and the nonlinear conservation law

10. State Space Methods

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

EXERCISES FOR SECTION 1.5

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

ME 391 Mechanical Engineering Analysis

Two Coupled Oscillators / Normal Modes

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Solutions from Chapter 9.1 and 9.2

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

MA Study Guide #1

Math Final Exam Solutions

Math 2214 Solution Test 1B Fall 2017

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

An random variable is a quantity that assumes different values with certain probabilities.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

1 Solutions to selected problems

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

arxiv: v1 [math.ca] 15 Nov 2016

2. Nonlinear Conservation Law Equations

Some Basic Information about M-S-D Systems

Elementary Differential Equations and Boundary Value Problems

INSTANTANEOUS VELOCITY

Math 10B: Mock Mid II. April 13, 2016

Example on p. 157

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

Lecture 10: The Poincaré Inequality in Euclidean space

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Solutions to Assignment 1

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

Second Order Linear Differential Equations

Let us start with a two dimensional case. We consider a vector ( x,

EECE 301 Signals & Systems Prof. Mark Fowler

Notes for Lecture 17-18

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Maxwell s Equations and Electromagnetic Waves

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

MATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

15. Vector Valued Functions

Chapter 7: Solving Trig Equations

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension

System of Linear Differential Equations

INDEX. Transient analysis 1 Initial Conditions 1

t 2 B F x,t n dsdt t u x,t dxdt

Boundedness and Stability of Solutions of Some Nonlinear Differential Equations of the Third-Order.

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Chapter 3 Boundary Value Problem

Fractional Method of Characteristics for Fractional Partial Differential Equations

Kinematics in two dimensions

A remark on the H -calculus

ON THE DEGREES OF RATIONAL KNOTS

Math 2214 Solution Test 1A Spring 2016

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Lie Derivatives operator vector field flow push back Lie derivative of

Math 116 Second Midterm March 21, 2016

Lecture 9: September 25

Linear Dynamic Models

Exam 1 Solutions. 1 Question 1. February 10, Part (A) 1.2 Part (B) To find equilibrium solutions, set P (t) = C = dp

Math 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm

From Complex Fourier Series to Fourier Transforms

Homework sheet Exercises done during the lecture of March 12, 2014

) were both constant and we brought them from under the integral.

Testing What You Know Now

Transcription:

Mah 527 Lecure 6: Hamilon-Jacobi Equaion: Explici Formulas Sep. 23, 2 Mehod of characerisics. We r o appl he mehod of characerisics o he Hamilon-Jacobi equaion: u +Hx, Du = in R n, u = g on R n =. 2 To avoid confusion, we use he following noaion: x, u z, x u p p =. 3 Du p Then we can re-wrie he equaion o where The characerisics ODEs hen are = x = D p F = ẋ D p H ż = D p F p = D p H ṗ ṗ Fx, z, p = 4 Fx, z, p 6 p + Hx, p. 5, 6 = p = D z F p D x F = p p = p + D p H p =D p H p Hp, x, 7. 8 D x H Now we r o solve he characerisic ODEs. Firs noice ha, since =, we can simpl use as he parameer insead of s. Thus he equaions become ẋ = D p H, 9 ż = D p H p Hx, p, N ṗ = D x H, p = p =. 2 I is clear he all we need o do is o solve he firs 3 equaions. Losing a bi rigor, we assume for now onl H is differeniable and sricl convex. We also assume H grows super-linearl a infini: Hx, p lim =+, 3 p ր p Now ake Legendre ransform: The z equaion hen becomes where q saisfies Therefore he soluion u is given b ux=ux + where x and x are relaed b Lx, v 6 sup p R n v p Hx, p. 4 ż = Lx, v 5 q =D p Hx, p. 6 Lxτ, vτ dτ. 7 ẋ = D p H =v, x =x. 8

2 Mah 527 Lecure 6: Hamilon-Jacobi Equaion: Explici Formulas To furher simplif he ssem, we noice ha implies which implies ha q, x minimizes wih x, x fixed. To see his, wrie Lx, v=v px, v Hx, px, v, and compue ẋ = D p H, ṗ = D x H 9 d d D vl+d x L = 2 Lxτ, vτ dτ 2 D v L= p + v D q p D p H D v p = p, 22 D x L=v D x p D x H D p H D x p = D x H, 23 where we have used v =D p H. Now he equaion ṗ = D x H gives wha we wan. Thus we see ha he Hamilon-Jacobi equaion can be solved as soon as we find ou he rajecories x and v. Below we will see ha in a special case, his can indeed be done in some sense. The Hopf-Lax formula. This special case is when H is independen of x, ha is H = HDu. The characerisic equaions can hen be furher simplified o ẋ = D p H, 24 ż = D p H p Hp=Lv, 25 ṗ = D x H =, 26 p n+ = p n+ N =. 27 We see ha p is a consan vecor along he characerisic curve, and as a consequence ẋ =D p H is a consan vecor, and herefore he characerisics x are sraigh lines. Furhermore we know ha he veloci q = ẋ is consan. Thus if x= and x=x, we mus have As a consequence d d z =Lv=L v =. 28 z=z+l = g+l. 29 Now he onl problem is ha is no known. Now hink of g as no merel an iniial funcion, bu as an inermediae record. In oher words, insead of saring a =, imagine our ssem sars from =, sa,. We consider all possible rajecories emanaing from some poin a =, passing a =, and finall reach ime a x. Think of g as he record of work done from = o =. Obviousl he correc rajecor should be he one ha is he minimizer among hem all. Remark. Noe ha he above explanaion means ha he rajecor ma no sa C as i crosses =. This should be expeced. Because in general he given g canno be he resul of a dnamical ssem wih H = HDu, ha is free paricle. The H producing g has o be dependen on x or even u. Thus here is no surprise ha his sudden change leads o a sudden change of direcion in he rajecor v. Following his idea, we reach he following Hopf-Lax formula: u, x=z= inf L + g. 3 R n Remark 2. I can be shown ha L grows superlinearl a infini. As a consequence, if we assume g o be Lipschiz coninuous, hen he infimum is acuall a minimum.

Sep. 23, 2 3 Remark 3. Noe ha convex funcions are coninuous. The proof can go roughl as follows. Firs one can show ha f he convex funcion is bounded, le he bound be denoed M. Then using he definiion of convexi we have, for an fixed x,, Leing α we see ha u + α u +α ux u u +2αM. 3 On he oher hand, for an x n x we have, b convexi This gives limsupux n ux. 32 x n x ux 2 [ux n+u2x x n ]. 33 ux 2 liminf[ux n +u2x x n ]. 34 x n x Coninui hen follows. One can in fac prove ha an convex funcion is Lipschiz coninuous, see e.g. B. Dacorogna Direc Mehods in he Calculus of Variaions, 2nd ed., Springer, 28, 2.3. Soluion of he H-J equaion. Now we show ha he Hopf-Lax formula ux, = inf L R n indeed solves he Hamilon-Jacobi equaion, albei onl almos everwhere. + g. 35 Remark 4. I is eas o see ha in general one canno expec he exisence of classical soluions due o possible inersecions of characerisics. There are hree hings o show.. u = g on R n =, 2. u, Du exis almos everwhere, 3. u + HDu= a.e. We show hem one b one.. u = g on R n =. Recall he formula: Taking = x we have u, x=min L u, x gx+l On he oher hand, we compue u, x = min L = gx+min gx max = gx max z = gx max + g. 36 limsupu, x gx. 37 ց + g + g gx L Lipg x L Lipg z Lz max w z Lz w B Lipg z = gx max w B Lipg Hw. 38

4 Mah 527 Lecure 6: Hamilon-Jacobi Equaion: Explici Formulas As H is coninuous, we have Thus ends he proof. 2. u, Du exis almos everwhere. I suffices o show ha u is Lipschiz wih respec o x and o. liminfu, x gx. 39 ց u is Lipschiz w.r.. x. We esimae u, xˆ u, x. Choose such ha u, x=l + g. 4 Then xˆ z u, xˆ u, x = min L + gz L g. 4 Taking z =xˆ x+ such ha xˆ z = we have u, xˆ u, x gxˆ x+ g Lipg xˆ x. 42 Similarl we can show u, x u, xˆ Lipg xˆ x. 43 The Lipschiz coninui of u hen follows. u is Lipschiz w.r... This follows from he following proper of he Hopf-Lax formula: u, x= min sl +us,. 44 R n s Tha his should hold is inuiivel ver clear following our derivaion of he formula. For a proof see Evans p. 26. Using his formula, we see ha esimaing u, x us, x is no differen han esimaing ux, gx. Thus a similar argumen as in Sep. gives u, x us, x C s. 45 3. u + HDu= a.e. Fix an q R n, we compue x+hq u + h, x+hq = min h L + u, h hlq+u, x. 46 This implies for all q R n. Therefore and u + q Du Lq u Du q Lq 47 u max Du q Lq = HDu 48 q u + HDu. 49 For he oher direcion ha is u + HDu, we onl need o find one q such ha or more specificall where is in he direcion of q. u + q Du Lq 5 u, x us, s Lq 5

Sep. 23, 2 5 As u is a minimum, o ge u, x us, somehing, we ge rid of he minimum in u, x. Take z such ha x z u, x=l + gz. 52 Now ha q = Then we have This gives As we ge and finishes he proof. is alread chosen, has o be on he line segmen connecing x and z. Thus we ake s = h, = s x + s z. 53 = z s = q 54 [ ] z u, x us, L + gz s L + gz s = sl. 55 u, x us, s u + x z u + Du 56 x z Du L 57