Fundamental lemma of the calculus of variations

Short description: Initial result in using test functions to find extremum

In mathematics, specifically in the calculus of variations, a variation δf of a function f can be concentrated on an arbitrarily small interval, but not a single point. Accordingly, the necessary condition of extremum (functional derivative equal zero) appears in a weak formulation (variational form) integrated with an arbitrary function δf. The fundamental lemma of the calculus of variations is typically used to transform this weak formulation into the strong formulation (differential equation), free of the integration with arbitrary function. The proof usually exploits the possibility to choose δf concentrated on an interval on which f keeps sign (positive or negative). Several versions of the lemma are in use. Basic versions are easy to formulate and prove. More powerful versions are used when needed.

Basic version

If a continuous function [math]\displaystyle{ f }[/math] on an open interval [math]\displaystyle{ (a,b) }[/math] satisfies the equality

[math]\displaystyle{ \int_a^b f(x)h(x)\,\mathrm{d}x = 0 }[/math]

for all compactly supported smooth functions [math]\displaystyle{ h }[/math] on [math]\displaystyle{ (a,b) }[/math], then [math]\displaystyle{ f }[/math] is identically zero.

Here "smooth" may be interpreted as "infinitely differentiable", but often is interpreted as "twice continuously differentiable" or "continuously differentiable" or even just "continuous", since these weaker statements may be strong enough for a given task. "Compactly supported" means "vanishes outside [math]\displaystyle{ [c,d] }[/math] for some [math]\displaystyle{ c }[/math], [math]\displaystyle{ d }[/math] such that [math]\displaystyle{ a\lt c\lt d\lt b }[/math]"; but often a weaker statement suffices, assuming only that [math]\displaystyle{ h }[/math] (or [math]\displaystyle{ h }[/math] and a number of its derivatives) vanishes at the endpoints [math]\displaystyle{ a }[/math], [math]\displaystyle{ b }[/math]; in this case the closed interval [math]\displaystyle{ [a,b] }[/math] is used.

Version for two given functions

If a pair of continuous functions f, g on an interval (a,b) satisfies the equality

[math]\displaystyle{ \int_a^b ( f(x) \, h(x) + g(x) \, h'(x) ) \, \mathrm{d}x = 0 }[/math]

for all compactly supported smooth functions h on (a,b), then g is differentiable, and g' = f  everywhere.

The special case for g = 0 is just the basic version.

Here is the special case for f = 0 (often sufficient).

If a continuous function g on an interval (a,b) satisfies the equality

[math]\displaystyle{ \int_a^b g(x) \, h'(x) \, \mathrm{d}x = 0 }[/math]

for all smooth functions h on (a,b) such that [math]\displaystyle{ h(a)=h(b)=0 }[/math], then g is constant.

If, in addition, continuous differentiability of g is assumed, then integration by parts reduces both statements to the basic version; this case is attributed to Joseph-Louis Lagrange, while the proof of differentiability of g is due to Paul du Bois-Reymond.

Versions for discontinuous functions

The given functions (f, g) may be discontinuous, provided that they are locally integrable (on the given interval). In this case, Lebesgue integration is meant, the conclusions hold almost everywhere (thus, in all continuity points), and differentiability of g is interpreted as local absolute continuity (rather than continuous differentiability). Sometimes the given functions are assumed to be piecewise continuous, in which case Riemann integration suffices, and the conclusions are stated everywhere except the finite set of discontinuity points.

Higher derivatives

If a tuple of continuous functions [math]\displaystyle{ f_0,f_1,\dots,f_n }[/math] on an interval (a,b) satisfies the equality

[math]\displaystyle{ \int_a^b ( f_0(x) \, h(x) + f_1(x) \, h'(x) + \dots + f_n(x) \, h^{(n)}(x) ) \, \mathrm{d}x = 0 }[/math]

for all compactly supported smooth functions h on (a,b), then there exist continuously differentiable functions [math]\displaystyle{ u_0,u_1,\dots,u_{n-1} }[/math] on (a,b) such that

[math]\displaystyle{ \begin{align} f_0 &= u'_0, \\ f_1 &= u_0 + u'_1, \\ f_2 &= u_1 + u'_2 \\ \vdots \\ f_{n-1} &= u_{n-2} + u'_{n-1}, \\ f_n &= u_{n-1} \end{align} }[/math]

everywhere.

This necessary condition is also sufficient, since the integrand becomes [math]\displaystyle{ (u_0 h)' + (u_1 h')' + \dots + (u_{n-1} h^{(n-1)})'. }[/math]

The case n = 1 is just the version for two given functions, since [math]\displaystyle{ f=f_0=u'_0 }[/math] and [math]\displaystyle{ f_1=u_0, }[/math] thus, [math]\displaystyle{ f_0-f'_1 = 0. }[/math]

In contrast, the case n=2 does not lead to the relation [math]\displaystyle{ f_0 - f'_1 + f''_2 = 0, }[/math] since the function [math]\displaystyle{ f_2 = u_1 }[/math] need not be differentiable twice. The sufficient condition [math]\displaystyle{ f_0 - f'_1 + f''_2 = 0 }[/math] is not necessary. Rather, the necessary and sufficient condition may be written as [math]\displaystyle{ f_0 - (f_1 - f'_2)' = 0 }[/math] for n=2, [math]\displaystyle{ f_0 - (f_1 - (f_2-f'_3)')' = 0 }[/math] for n=3, and so on; in general, the brackets cannot be opened because of non-differentiability.

Vector-valued functions

Generalization to vector-valued functions [math]\displaystyle{ (a,b)\to\mathbb{R}^d }[/math] is straightforward; one applies the results for scalar functions to each coordinate separately, or treats the vector-valued case from the beginning.

Multivariable functions

If a continuous multivariable function f on an open set [math]\displaystyle{ \Omega\subset\mathbb{R}^d }[/math] satisfies the equality

[math]\displaystyle{ \int_\Omega f(x) \, h(x) \, \mathrm{d}x = 0 }[/math]

for all compactly supported smooth functions h on Ω, then f is identically zero.

Similarly to the basic version, one may consider a continuous function f on the closure of Ω, assuming that h vanishes on the boundary of Ω (rather than compactly supported).

Here is a version for discontinuous multivariable functions.

Let [math]\displaystyle{ \Omega\subset\mathbb{R}^d }[/math] be an open set, and [math]\displaystyle{ f\in L^2(\Omega) }[/math] satisfy the equality

[math]\displaystyle{ \int_\Omega f(x) \, h(x) \, \mathrm{d}x = 0 }[/math]

for all compactly supported smooth functions h on Ω. Then f=0 (in L², that is, almost everywhere).

Applications

This lemma is used to prove that extrema of the functional

[math]\displaystyle{ J[y] = \int_{x_0}^{x_1} L(t,y(t),\dot y(t)) \, \mathrm{d}t }[/math]

are weak solutions [math]\displaystyle{ y:[x_0,x_1]\to V }[/math] (for an appropriate vector space [math]\displaystyle{ V }[/math]) of the Euler–Lagrange equation

[math]\displaystyle{ {\partial L(t,y(t),\dot y(t)) \over \partial y} = {\mathrm{d} \over \mathrm{d}t} {\partial L(t,y(t),\dot y(t)) \over \partial \dot y} . }[/math]

The Euler–Lagrange equation plays a prominent role in classical mechanics and differential geometry.

References

• Jost, Jürgen; Li-Jost, Xianqing (1998), Calculus of variations, Cambridge University 
• Gelfand, I.M.; Fomin, S.V. (1963), Calculus of variations, Prentice-Hall  (transl. from Russian).
• Hestenes, Magnus R. (1966), Calculus of variations and optimal control theory, John Wiley 
• Giaquinta, Mariano; Hildebrandt, Stefan (1996), Calculus of Variations I, Springer