Calculus of Variations Basics

Similar documents
Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Canonical transformations

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

PHYS 705: Classical Mechanics. Calculus of Variations II

6.3.4 Modified Euler s method of integration

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Affine transformations and convexity

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Poisson brackets and canonical transformations

2 Finite difference basics

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

More metrics on cartesian products

Quantum Mechanics I Problem set No.1

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Analytical classical dynamics

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

PHYS 705: Classical Mechanics. Canonical Transformation II

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Lecture 2: Numerical Methods for Differentiations and Integrations

Implicit Integration Henyey Method

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Modelli Clamfim Equazioni differenziali 22 settembre 2016

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Lecture 21: Numerical methods for pricing American type derivatives

Module 3: Element Properties Lecture 1: Natural Coordinates

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Lecture 12: Discrete Laplacian

Affine and Riemannian Connections

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Lecture 10 Support Vector Machines II

Bernoulli Numbers and Polynomials

Lecture Note 3. Eshelby s Inclusion II

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Integrals and Invariants of Euler-Lagrange Equations

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

Difference Equations

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

The Feynman path integral

APPENDIX A Some Linear Algebra

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Convexity preserving interpolation by splines of arbitrary degree

Uncertainty and auto-correlation in. Measurement

NUMERICAL DIFFERENTIATION

Army Ants Tunneling for Classical Simulations

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:

), it produces a response (output function g (x)

Inductance Calculation for Conductors of Arbitrary Shape

Société de Calcul Mathématique SA

Appendix B. Criterion of Riemann-Stieltjes Integrability

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

10. Canonical Transformations Michael Fowler

Review of Newtonian Mechanics

Spin-rotation coupling of the angularly accelerated rigid body

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

PETER HENR IC I TECHNICAL REPORT NO. CS 137 JULY COMPUTER SCIENCE DEPARTMENT School of Humanities and Sciences STANFORD UN IVERS ITY

Classical Mechanics Virtual Work & d Alembert s Principle

Complete subgraphs in multipartite graphs

PHYS 705: Classical Mechanics. Newtonian Mechanics

1 Matrix representations of canonical matrices

Fundamental loop-current method using virtual voltage sources technique for special cases

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Math1110 (Spring 2009) Prelim 3 - Solutions

Numerical Solution of Ordinary Differential Equations

Chapter 4 The Wave Equation

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

CALCULUS CLASSROOM CAPSULES

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Chapter 11: Simple Linear Regression and Correlation

Integrals and Invariants of

Learning Theory: Lecture Notes

PHZ 6607 Lecture Notes

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Mechanics Physics 151

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

On the symmetric character of the thermal conductivity tensor

12. The Hamilton-Jacobi Equation Michael Fowler

Three views of mechanics

k t+1 + c t A t k t, t=0

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION

Lecture 20: Noether s Theorem

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Iterative General Dynamic Model for Serial-Link Manipulators

Exercise Solutions to Real Analysis

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Transcription:

Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y 2,,y n ), 1.1) begnnng wth the case where 1.1) depends on a sngle functony and hence reduces to J [y] F x,y,y ). 1.2) We assume that all admssble curves are smooth 1, but, the end ponts of the curves for whch 1.2) s defned can move n an arbtrary way. Now, we also defne the dstance between two curvesy yx) and y y x) as ρy,y ) max y y +maxy y +ρp0,p0 )+ρp 1,P1 ), 1.3) wherep 0 andp 0 denote the left-hand end ponts of the curvesy andy, respectvely, andp 1 and P 1 denote the rght-hand end ponts. In general, the functons y and y are defned on dfferent ntervalsi and I. Thus, n order for 1.3) to make sense, we have to extendy and y onto some nterval contanng bothi andi e.g., regon[,x 1 +δx 1 ] as shown n Fg. 1.1). Now, let y and y be two neghborng curves, n the sense of dstance 1.3), and let hx) y y. Moreover, let P 0,y 0 ) andp 1 x 1,y 1 ) 1 We say the functon yx) s smooth n an nterval [,x 1 ], f both the functon as well as t s dervatves are contnuous n[,x 1 ]. 1

Fgure 1.1: Functonsy andy and ther correspondng domans denote the end ponts of the curvey, whle the end ponts of the curvey are denoted by P 0 +δ,y 0 +δy 0 ) andp 1 x 1 +δx 1,y 1 +δy 1 ). The correspondng varaton δj of the functonal J [y] s defned as the expresson whch s lnear nh,h,δ,δy 0,δx 1,δy 1, and dffers from the ncrement J J [y +h] J [y] by a quantty of the order hgher than1relatve toρy,y+h). Snce J +δx 1 +δ F x,y +h,y +h ) + [F x,y +h,y +h ) F x,y,y )] +δx 1 x 1 F x,y +h,y +h ) F x,y,y ) x0 +δ F x,y +h,y +h ), t follows by usng Taylor s theorem and lettng the symbol denote equalty except for terms of order hgher than1relatve toρy,y+h) that J [F y x,y,y )h x,y,y )h ] + F x,y,y ) xx1 δx 1 F x,y,y ) xx0 δ [ F y d ] hx)+ ) h xx1 h xx0 + F xx1 δx 1 F xx0 δ. The term nsde the round brackets of the above equaton s due to ntegraton by parts. However, t s clear from Fg.1.1 that h ) δy 0 y )δ, hx 1 ) δy 1 y x 1 )δx 1, 2

where has the same meanng as before, and hence [ δj F y d ] hx)+ δy xx 1 x + F y )δx xx 1 x, where we defne δx xx δx and, δy xx δy 0,1). Next, we return to the more general functonal 1.1), whch depends onnfunctonsy 1,y 2,,y n. From the analogy, varaton for the general functonal s gven as δj F y d ] } ) xx F y h x) + δy xx 1 1 + F F xx y 0 y δx 1.4) where, as before, we defne δx xxj δx j and, δy xx δy j j 0,1). x 1.2 Varaton n terms of the Canoncal Varables We now wrte an even more concse formula for the varaton 1.4), at the same tme ntroducng some new mportant deas. Let and suppose that the Jacoban p 1,,n), p 1,,p n ) y 1,,y n ) det Fy y k s nonzero so that the mappng s nvertble). Then we can solve the equatons fory 1,,y n as functons of the varables x,y 1,,y n,p 1,,p n. 1.5) Next we can express the functonf x,y 1,y 2,,y n,y 1,y 2,,y n) n terms of a new functon Hx,y 1,,y n,p 1,,p n ) related tof by the formula H F + y F + p y, 1.6) where y are regarded as the functons of the varables. The functon H s called the Hamltonan functon) correspondng to the functonal J [y 1,,y n ]. In ths way, we can make a local transformaton from the varables x,y 1,y 2,,y n,y 1,y 2,,y n,f to the new quanttes x,y 1,,y n,p 1,,p n,h called the canoncal varables. In terms of the canoncal varables, varaton can be wrtten n the form δj F y dp ] } h x) + p δy Hδx) xx 1 x. 1.7) 3

1.3 Euler-Lagrange Equatons & Correspondng Constrants The necessary condton for the dfferentable functonal J [y 1,y 2,,y n ] to have an extremum for[ŷ 1,ŷ 2,,ŷ n ] s that ts varaton vansh for[ŷ 1,ŷ 2,,ŷ n ],.e., that δj [h] 0 for[ŷ 1,ŷ 2,,ŷ n ] and all admssbleh. So, from 1.4) and 1.7), the necessary condtonδj [h] 0 mples F y d ] } F y h x) + δy xx 1 x + F y dp ] } h x) + F ) y δx xx 1 x p δy Hδx) xx 1 x 0. 1.8) Suppose the functonal J [y 1,y 2,,y n ] has an extremum n certan class of admssble curves) for some curve y ŷ jonng the ponts P 0,y 0 1,,y0 n) andp1,y 1 1,,y1 n). Then, snce J [y 1,y 2,,y n ] has an extremum for ŷ compared to all admssble curves, t certanly has an extremum forŷ compared to all curves wth fxed 2 end pontsp 0 andp 1. Therefore, F y d ] } F y h x) Because of the ndependence, the above equaton smply reduces to F y dp ] } h x) 0. F y d ] } F y h x) F y dp ] } h x) 0 1.9) for allvalues. Lemma 1. If αx) s contnuous n[a,b], and f b a αx)hx) 0 for every contnuous functonhx) such thatha) hb) 0, thenαx) 0 for allxn[a,b]. 2 When we say fxed, t obvously meansh ) h ) 0,δx j 0 andδy j 0. 4

So, from the above lemma and 1.9), F y d ) F y F y dp ) 0, 1,,n). 1.10) The above equatons known as Euler-Lagrange equatons represent a system of n second-order 3 dfferental equatons and ts soluton can be obtaned from 2n arbtrary constrants gven as ) δy + F F xx y 0 y δx xx0 p δy xx0 Hδx xx0 0, 1,,n), 1.11) δy + F xx 1 ) y δx xx1 p δy xx1 Hδx xx1 0, 1,,n). 1.12) The above equaton s from 1.8) and due to the ndependence ofδ andδx 1. Before concludng, a few mportant relatons between F and H are gven below. These relatons can be derved from 1.6). H x F x, H y F y and H p y, 1,,n) 1.13) Substtutng the above relatons nto 1.10) gves2n frst order dfferental equatons dy H p and dp H y, 1,,n). 1.14) Once agan, 2n arbtrary constrants requred to solve the above equatons are gven by 1.11) and 1.12). The materal presented here s from the book: [1] Gelfand, I. M., Fomn, S. V., Calculus of varatons, Prentce-Hall, Inc., 1963 3 or2n frst order dfferental equatons n the case of canoncal varables, whch are gven by 1.14). 5