Lecture 10: Euler s Equations for Multivariable

Similar documents
Poisson brackets and canonical transformations

1 Matrix representations of canonical matrices

Section 8.3 Polar Form of Complex Numbers

Mechanics Physics 151

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

Integrals and Invariants of Euler-Lagrange Equations

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

PHYS 705: Classical Mechanics. Calculus of Variations II

Economics 101. Lecture 4 - Equilibrium and Efficiency

Lecture 4. Instructor: Haipeng Luo

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

Mechanics Physics 151

Integrals and Invariants of

Week 5: Neural Networks

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

Expectation propagation

12. The Hamilton-Jacobi Equation Michael Fowler

From Biot-Savart Law to Divergence of B (1)

Chapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Chapter Newton s Method

NUMERICAL DIFFERENTIATION

10.34 Fall 2015 Metropolis Monte Carlo Algorithm

coordinates. Then, the position vectors are described by

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

Section 3.6 Complex Zeros

Canonical transformations

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Spring Force and Power

( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.

Société de Calcul Mathématique SA

Lecture 20: Noether s Theorem

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

4DVAR, according to the name, is a four-dimensional variational method.

10. Canonical Transformations Michael Fowler

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

Linear Feature Engineering 11

Mean Field / Variational Approximations

PHYS 705: Classical Mechanics. Canonical Transformation II

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

Unit 5: Quadratic Equations & Functions

Group Analysis of Ordinary Differential Equations of the Order n>2

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

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

One Dimension Again. Chapter Fourteen

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

A Quantum Gauss-Bonnet Theorem

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Error Probability for M Signals

8.6 The Complex Number System

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

Solutions to Problem Set 6

Maximum Likelihood Estimation

8.323 Relativistic Quantum Field Theory I

Lecture Notes on Linear Regression

Radar Trackers. Study Guide. All chapters, problems, examples and page numbers refer to Applied Optimal Estimation, A. Gelb, Ed.

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

More metrics on cartesian products

Lecture 10 Support Vector Machines II

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3 Basic boundary value problems for analytic function in the upper half plane

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD

Lecture 10 Support Vector Machines. Oct

Lecture 12: Discrete Laplacian

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Lecture 19. Endogenous Regressors and Instrumental Variables

Restricted divisor sums

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

Module 3: Element Properties Lecture 1: Natural Coordinates

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

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

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

14 The Postulates of Quantum mechanics

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Estimation: Part 2. Chapter GREG estimation

Some modelling aspects for the Matlab implementation of MMA

Generalized Linear Methods

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Work is the change in energy of a system (neglecting heat transfer). To examine what could

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Feature Selection: Part 1

Assortment Optimization under MNL

Differentiating Gaussian Processes

A new Approach for Solving Linear Ordinary Differential Equations

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

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Problem Solving in Math (Math 43900) Fall 2013

Linear Approximation with Regularization and Moving Least Squares

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Lecture 10. Reading: Notes and Brennan Chapter 5

Solutions Homework 4 March 5, 2018

Which Separator? Spring 1

Calculus of Variations Basics

Appendix B. The Finite Difference Scheme

Transcription:

Lecture 0: Euler s Equatons for Multvarable Problems Let s say we re tryng to mnmze an ntegral of the form: {,,,,,, ; } J f y y y y y y d We can start by wrtng each of the y s as we dd before: y (, ) ( 0) ( ) α y + αη Remember that y (0) s the functon we want J f { y, y, y, y,, y, y ; } d

Applyng the chan rule, we fnd: J f y f y f y f y d d α + y α y α + + y y f y + + + f y d y y Whch means we have ntegrals of the type encountered n the one-varable problem J So, for to be 0, each ntegral must be zero Ths means we have to satsfy Euler equatons, of the form: f d f y d y 0

Euler s Equaton wth Constrants In some cases, we want to restrct the set of paths one can take between ponts and, and fnd the mnmum of J subect to ths constrant For eample, we mght want to fnd the shortest dstance between two ponts on a surface wth a fed shape If ths happens, our dervaton of the many-varable Euler Equaton s altered. We stll have: J f y f y f y f y d d α + y α y α + + y y f y f y + + + d 0 when α 0 y y but now the y are not ndependent, so we can t assume that each ntegral s zero ndvdually

Usng ntegraton by parts on each ntegral gves: J f y d f y f y d f y d d α y α d y α + y d y y + + + f d f y d y d y y f d f y f d f d d y d y + y d y f d f y + + + d y d y

In general, each constrant can be epressed as a functon relatng the varables to each other: g y ; 0 { } and there may be an arbtrary number m of such constrant equatons (well, not qute arbtrary m must be less than, the number of varables n the problem) For any one of the constrant equatons we have: g { y } ; 0 y dg dα η dα y α y Ths last relaton holds even f we multply by an arbtrary functon of : λ ( ) η dα 0 y 0

We can sum over all m constrants, and ntegrate from to, and stll get zero: m λ ( ) y ow we can add ths rather eotc-lookng form of the number 0 to our equaton for J : ηdα m J f d f λ ( ) α + + y d y y But snce the varables are related, we stll can t assume that each term n [] s zero ndvdually 0 η d

However, we can always redefne the varables of the problem such that the frst -m of them are ndependent The remanng m varables are related by the m constrant equatons For the -m ndependent varables, we must have the terms n [] be zero ndvdually: m f d f + + λ ( ) 0 for < m y d y y Snce we ve sad nothng so far about the value of the functons λ (), we are free to choose ones that are convenent. In partcular, we choose them such that: m f d f + + λ ( ) 0 for > m y d y y

Combnng the epressons on the prevous sldes, we arrve at Euler s Equatons when constrants are appled: m f d f + + λ ( ) y d y y 0 for any You may be concerned that there are +m varables n the above equatons, but only equatons avalable But remember that the m equatons of constrant are also avalable So we can solve for everythng ncludng the λ s! Keep n mnd that we made a specal choce for the λ s durng the dervaton. When we apply these deas to physcs we ll see what they represent

Eample: Geodesc on a Sphere Let s say we want to fnd the shortest dstance between two ponts, but we have to move on the surface of a sphere between them That means we need to mnmze the ntegral: + + z ( ) ( ) ( ) J ds d + dy + dz y where only paths that satsfy the followng constrant are consdered: d ( ) g, y, z + y + z R 0 f

Snce f depends on two varables n ths problem, we have the followng Euler equatons: f d f λ 0 y d y y f d f λ 0 z d yz z The frst equaton gves: d y 0 y 0 d λ + y + z d z 0 λz 0 d + y + z

ote that we can solve each of these for λ: d y d z λ y d z d + y + z + y + z y y y + y z z 3/ + y + z ( + y + z ) z z z + y z y 3/ + y + z ( + y + z ) ( ) + + z y y z y y y z z ( ) yz + y + z z z y z z + zy y + zy y zy zy z zy z z yz + yz y + yz z yz z yy z z

zy + zy z zy z z yz + yz y yy z z ote that we haven t used any nformaton from the equaton of constrant yet. We can do so f we take a dervatve of g: dg + yy + zz 0 d yy + zz Pluggng ths back n above gves: ( ) ( ) zy + z y zz + yy yz + y z zz + yy zy z y yz y z One can verfy that a plane passng through the center of the sphere ( A + By z ) s the soluton to the above equaton

Here s the verfcaton: z A + By; z A + By ; z By [ ] [ ] A + By y A + By y Byy By y Ay + Byy Ay By y Byy By y ote that the tet does the same problem (Eample 6.4) wthout usng λ In ths case that s easer, snce we re not nterestng n what the functon λ s But sometmes for physcsts, λ s qute mportant we ll see why n a few lectures