16. Introduction to Calculus of Variations#
Calculus of variation deals with variations - i.e. “small changes” - of functions and functionals.
The meaning of the term functional may vary on the subfield of interest. In the field of calculus of variation, a functional can be defined as a function of function, i.e. a function whose argument is another function.
Fields and applications
Fields and applications related to calculus of variations (give some examples below):
gradient-based techniques like some methods in:
optimization, either free or constrained (via Lagrange multiplier methods)
sensitivity
classical mechanics and physics in general:
analytical mechanics: Lagrangian formulation and Hamiltonian formulation of classical mechanics
…
Examples
Lagrange equations for general problem
examples:
brachistochrone for minimum time,…
catenary, i.e. static solution of wire and cables with neglibile bending stiffness
isoperimetric inequality, i.e. circle is the plane closed curve with given perimeter enclosing the largest area
sensitivity of results to parameters. Some interesting sensitivity, both in time and trasnformed domains
characteristics of a system:
equilibria
eigenvalues
…
optimal control methods
16.1. Lagrange equations#
Given the functional \(S\), with arguments a function \(q(t)\) and the independent variable \(t\),
its variation w.r.t. the function \(q(t)\) reads
where the function \(w(t)\) is arbitrary, among those satisfying the constraint of the problems: as an example here, if the function \(q(t)\) has prescribed values \(q^*\) for some values of the independent variable, \(t^*\), the variation \(w(t)\) of the function \(q(t)\) is zero there, \(w(t^*)\) so that the variated function \(q(t) + \varepsilon w(t)\) satisfies the constraint as well, i.e. \(q(t^*) + \varepsilon w(t^*) = q^*\).
Variation involves only small changes of function arguments, since these ones are the elements that can be effectively changed, while the independent variable is not.
Direct computation of the variation gives
The solution depends on the boundary conditions at the extreme points \(t_0\), \(t_1\). If the value of the function \(q(t)\) is prescribed in \(t_0\) and \(t_1\), \(q(t_0) = q_0\), \(q(t_1) = q_1\), then its variation is zero, \(w(t_0) = w(t_1) = 0\), for the reason that has been discussed above. The variation of the functional with prescribed boundary values of the argument function thus reads
having called \(w(t) =: \delta q(t)\) to stress that is the variation of function \(q(t)\). This notation - it’s just notation, it has no special properties - could be useful if the functional depends on several arguments.
Stationary conditions, \(\delta S = 0\). Stationary condition of the functional \(S\) implies that \(\delta S = 0\) for all the possible variations of the argument function, \(\forall \delta q(t)\). This condition implies that the integrand is identically zero, i.e. Lagrange equations,
Higher-order derivatives
Method 1. If the Lagrangian function \(L\) depends on higher order derivatives,
it’s possible to recast the problem defining the \(n\)-dimensional function, \(\mathbf{q}(t)\),
With some abuse of notation in \(L\), the functional \(S\) can be recasted as
todo Add constraints on components of \(\mathbf{q}(t)\)?
Repeating the computation, the variation of the functional reads
Method 2. …
16.1.1. Euler-Beltrami equation#
If the Lagrangian function \(L\) is not an explicit function of the independent variable, \(L(q'(x),q(x))\), Euler-Berltrami equation follows from the derivative of the Lagrangian,
and thus
Note 1. While Lagrange equations are a set of \(N\) equations if the functional depends on \(N\) argument functions \(q_k(t)\), \(k=1:N\), Euler-Beltrami equation is an equation only. Indeed for multiple argument functions
where Einstein’s summation notation of repeated index is used. Euler-Beltrami thus reads
Note 2. If the Lagrangian function is an explicit function of the independent variable \(x\), \(L(q'(x), q(x), x)\), it’s not hard to realize that the derivative of the Lagrangian function, along with the use of the Lagrange equation, gives
Example 16.1 (Euler-Beltrami with \(L(q'(x),q(x),x)\), Hamiltonian, energy and E.Noether)
Euler-Beltrami equation shows that if \(L(q'(x), q(x))\), thus \(L - q' \partial_{q'} L\) is constant, (or an integral of motion in dynamics). In analytical mechanics (Lagrange mechanics, Hamiltonian mechanics), Lagrangian and Hamiltonian functions of a system read
having used the common definition of the generalized momenta \(p_k := \frac{\partial L}{\partial \dot{q}^k}\). It should be immediate to realize that the Hamiltonian is just the quantity appearing in Euler-Beltrami equation (or in its “modified version” if \(\partial_t L \ne 0\)), and thus
In mechanics, if \(\partial_t L = 0\), the Hamiltonian is a constant of motion. In this case, it can be prove that the Hamiltonian is equal to the eneergy of the system.
Classical examples.
Example 16.2 (Brachistochrone)
Find the trajectory…
Elementary length: \(ds = v \, dt\)
Energy: \(E(y) = \frac{1}{2} m v^2 - m g y + C\). Setting \(E = 0\) at starting point, from rest, at \(y_0 = 0\), it implies \(C=0\); thus \(v = \sqrt{2gy}\)
\(x(s), \, y(s)\),
\(ds = \sqrt{dx^2 + dy^2} = \sqrt{1+y'^2(x)} \, dx\)
The Lagrangian doesn’t explicitly depend on \(x\), thus Euler-Beltrami equation can be used. Partial derivative of the Lagrangian function w.r.t. \(q'\) reads
and thus Euler-Beltrami equation reads
Squaring \(2 g C^2 = \frac{1}{y(1+y'^2)}\), it’s possible to write
Making the substitution \(y(x)' = \dots\)
Example 16.3 (Catenary)
Example 16.4 (Isoperimetric problem)