Math 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx

Similar documents
4. Calculus of Variations

1.9 C 2 inner variations

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Problem set 1: Solutions Math 207B, Winter 2016

Calculus of Variations

The Basic Functional 2 1

Calculus of Variations: The Direct Approach

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Part III. Calculus of variations

Chapter 4. Additional Variational Concepts

Solution to HW 4, Ma 1c Prac 2016

1.3 The Lemma of DuBois-Reymond

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

MATH , Calculus 2, Fall 2018

Reversing the Chain Rule. As we have seen from the Second Fundamental Theorem ( 4.3), the easiest way to evaluate an integral b

Improper Integrals. MATH 211, Calculus II. J. Robert Buchanan. Spring Department of Mathematics

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

Taylor Polynomial Inequalities

Variational Problems

Construction of Gauss Quadrature Rules

Review on Integration (Secs ) Review: Sec Origins of Calculus. Riemann Sums. New functions from old ones.

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1. 1 [(y ) 2 + yy + y 2 ] dx,

MatFys. Week 2, Nov , 2005, revised Nov. 23

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

Section 17.2 Line Integrals

Introduction to the Calculus of Variations

Homework 4. (1) If f R[a, b], show that f 3 R[a, b]. If f + (x) = max{f(x), 0}, is f + R[a, b]? Justify your answer.

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Solutions to Homework Set 3

Generalizations of the Basic Functional

(4.1) D r v(t) ω(t, v(t))

Joint distribution. Joint distribution. Marginal distributions. Joint distribution

Solutions to Assignment 1

Math 100 Review Sheet

Thin-Plate Splines. Contents

Chapter 2. Constraints, Lagrange s equations

Math 31S. Rumbos Fall Solutions to Assignment #16

21.6 Green Functions for First Order Equations

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

Theoretical foundations of Gaussian quadrature

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

1 1D heat and wave equations on a finite interval

Math 5440 Problem Set 3 Solutions

Calculus of variations with fractional derivatives and fractional integrals

Anonymous Math 361: Homework 5. x i = 1 (1 u i )

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

Math 5440 Problem Set 3 Solutions

cos 3 (x) sin(x) dx 3y + 4 dy Math 1206 Calculus Sec. 5.6: Substitution and Area Between Curves

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

Some Methods in the Calculus of Variations

INDIAN INSTITUTE OF TECHNOLOGY BOMBAY MA205 Complex Analysis Autumn 2012

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

ES 111 Mathematical Methods in the Earth Sciences Lecture Outline 1 - Thurs 28th Sept 17 Review of trigonometry and basic calculus

INTRODUCTION TO INTEGRATION

The Regulated and Riemann Integrals

The Form of Hanging Slinky

Lecture 1. Functional series. Pointwise and uniform convergence.

Orthogonal Polynomials

First variation. (one-variable problem) January 14, 2013

10 Vector Integral Calculus

The usual algebraic operations +,, (or ), on real numbers can then be extended to operations on complex numbers in a natural way: ( 2) i = 1

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Calculus 2: Integration. Differentiation. Integration

Math Fall 2006 Sample problems for the final exam: Solutions

Math Calculus with Analytic Geometry II

Math 554 Integration

Improper Integrals, and Differential Equations

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

1 The Riemann Integral

Conservation Law. Chapter Goal. 5.2 Theory

Physics 161: Black Holes: Lecture 6: 14 Jan 2011

f(a+h) f(a) x a h 0. This is the rate at which

Advanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives

Overview of Calculus I

MA 124 January 18, Derivatives are. Integrals are.

3.4 Numerical integration

Fundamental Theorem of Calculus

Best Approximation in the 2-norm

38 Riemann sums and existence of the definite integral.

Math 259 Winter Solutions to Homework #9

MATH 253 WORKSHEET 24 MORE INTEGRATION IN POLAR COORDINATES. r dr = = 4 = Here we used: (1) The half-angle formula cos 2 θ = 1 2

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

Chapter 8: Methods of Integration

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

Calculus II: Integrations and Series

Problem 1. Brachistochrone time-of-flight

Consequently, the temperature must be the same at each point in the cross section at x. Let:

221B Lecture Notes WKB Method

Abstract inner product spaces

4181H Problem Set 11 Selected Solutions. Chapter 19. n(log x) n 1 1 x x dx,

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

Transcription:

Mth 5 2006-2007) Yum-Tong Siu Lgrnge Multipliers nd Vritionl Problems with Constrints Integrl Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y] = F x,y,y )dx with y) = A nd yb) = B subject to the dditionl integrl constrint tht K[y] = G x,y,y )dx = l, where l is given constnt. Suppose we hve n extreml y = yx) for this vritionl problem. To derive necessry condition for the extreml, we embed it in 2-prmeter fmily y = yx,s,t) so tht the given extreml y = yx) corresponds to s,t) = 0, 0). The reson why we need 2-prmeter fmily insted of -prmeter fmily is tht the fmily y = ys,t) hs to stisfy the integrl constrint ) K[y] = G x,yx,s,t), x yx,s,t) dx = l for ll s,t). Consider now J[y] s function of two vribles s,t) subject to the condition K[y] = l. Since J[y] hs criticl point t s,t) = 0, 0) subject to the condition tht the point s,t) lies on the curve K[y] = l in the spce of the two rel vribles s,t), it follows tht the grdient of J[y] with respect to s,t) is proportionl to the grdient of K[y] with respect to s,t) t the point s,t) = 0, 0). Thus there exists some rel number λ such tht s J[y] = λ s K[y], t J[y] = λ t K[y] t the point s,t) = 0, 0). In other words, we hve F y ) ) b s y dx = λ G y ) ) s y dx, F y ) ) b t y dx = λ G y ) ) t y dx,

Mth 5 2006-2007) Yum-Tong Siu 2 The left-hnd side F y ) ) s y dx of the first eqution is the derivtive of J with respect to the vector y in s the spce of functions. The left-hnd side F y ) ) t y dx of the second eqution is the derivtive of J with respect to the vector y t in the spce of functions. The right-hnd side G y ) ) s y dx of the first eqution is the derivtive of K with respect to the vector y in s the spce of functions. The right-hnd side G y ) ) t y dx of the second eqution is the derivtive of K with respect to the vector t y in the spce of functions. The constnt λ is lredy determined by the first eqution which sys tht the component of the grdient of J in the direction of the vector y in s the spce of functions is equl to λ times the component of the grdient of K in the direction of the vector y in the spce of functions. The second s eqution sys tht s result the component of the grdient of J in the direction of ny other vector y in the spce of functions is equl to the t sme constnt λ times the component of the grdient of K in the direction of the vector y in the spce of functions. t As consequence we cn sy tht the full grdient of J in the spce of functions is equl to λ times the full grdient of K in the spce of functions. In other words, F y ) λ G y )) δy)dx = 0

Mth 5 2006-2007) Yum-Tong Siu 3 for ll δy with δy) ) = δy) b) = 0. Therefore we get the Euler-Lgrnge eqution ) F y = λ G y ) for some constnt λ which is known s the Lgrnge multiplier. Another wy to write it is F λg) y dx F λg) y = 0. Besides the two initil conditions y) = A nd yb) = B to determine the two constnt of integrtions for the solution of the second-order differentil eqution ), we hve n extr unknown λ which will be determined by the integrl constrint G x,y,y )dx = l. Exmple. Given > 0 nd l > 0. Find n extreml for the vritionl problem subject to y ) = y) = 0 nd K[y] = J[y] = ydx + y 2 dx = l. Solution. The Euler-Lgrnge eqution is ) y λ + y y 2 )) y λ + y dx y 2 = 0 for some Lgrnge multiplier λ. We cn rewrite it s Integrting it once, we get y + λ d dx = 0. + y 2 y x + λ = C. + y 2 Squring to remove the squre root, we get x C ) 2 = y 2 + y 2 = λ ) + y 2.

Mth 5 2006-2007) Yum-Tong Siu 4 As result, y 2 = x C ) 2 = + y 2, + y 2 = x C ) 2, + y 2 =, x C ) 2 = x C ) 2 y = ± x C λ x C ) 2 x C ) 2 x C ) 2, dx. Let x C λ = cos θ. Then cos θ sin θ)λdθ y = ± = λ sin θ + C 2. cos2 θ Eliminting θ, we get x C ) 2 + y C 2 ) 2 =. The constnts C, C 2, nd λ re to be determined by y ) = y) = 0 nd K[y] = + y 2 dx = l. Pointwise Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y,z] = F x,y,z,y,z )dx with y) = A nd yb) = B subject to the dditionl pointwise constrint tht gx,y,z) = 0. Suppose we hve n extreml y = yx),z = zx) for this vritionl problem. To derive necessry condition for the extreml, we embed it in -prmeter fmily y = yx,t),z = zx,t) which stisfy the pointwise constrint gx,yx,t),zx,t)) 0 for ll x,t so tht the given

Mth 5 2006-2007) Yum-Tong Siu 5 extreml y = yx),z = zx) corresponds to t = 0. Tking the first vrition of J, we get ) δj = F y ) δy)dx + F z ) dx F z δz) dx, where δj = dj nd δy = y nd δz = z. Differentiting t=0 t=0 t=0 dt gx, yx, t), zx, t)) 0 with respect to t yields for ll x. Solving for δz, we get Putting it into ) gives us 0 = δj = t g y δy) + δz) 0 δz = g y δy). F y g y t F z )) dx F z δy) dx. By the rbitrriness of δy, we conclude tht the integrnd vnishes nd get F y g y F z ) dx F z = 0 which cn be rewritten s F y g y = F z dx F z. Let λx) be either of the two equl sides of the bove eqution. We get F λg) y dx F λg) y = 0, F λg) z dx F λg) z = 0. Heuristiclly, we cn lso interpret the pointwise condition s -prmeter fmily of integrl constrint δ ξ ) x)gx,y,z)dx = 0 with prmeter ξ [,b], where δ ξ x) is the Dirc delt t ξ so tht we end up with one constnt multiplier λξ) for ech of such n integrl constrint.