Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Endpoint Variations and Transversality Conditions

The classical fixed-endpoint problem in the calculus of variations assumes that the initial time, final time, initial state, and final state are all prescribed. Under those assumptions, the first variation leads to the Euler–Lagrange equation, and the endpoint terms vanish because the endpoint variations are zero.

Many practical problems are less restrictive. A system may be required to reach a specified terminal state without prescribing the arrival time. Alternatively, the final time may be fixed while the final state is free. More general problems may allow both endpoint time and endpoint state to vary subject to geometric or algebraic constraints.

These cases produce additional necessary conditions known as transversality conditions. This section develops the endpoint variation identities, derives the free-final-time condition, discusses the free-final-state condition, and illustrates the results through a mechanical example.

Review of the Basic Functional

Consider the functional

J[x]=t0tfL(x(t),x˙(t),t)dt.J[x] = \int_{t_0}^{t_f} L\bigl(x(t),\dot{x}(t),t\bigr)\,\mathrm{d} t.

The variable tt is a dummy variable of integration. Therefore, after integration, J[x]J[x] is not a function of the dummy variable tt. Instead, the value of the functional depends on:

If tft_f is free, then changing tft_f can change the numerical value of the functional. Consequently, the first variation must include a term involving δtf\delta t_f.

Fixed and Free Endpoint Data

For the most restrictive scalar problem, all endpoint data are fixed:

t0=fixed,tf=fixed,x(t0)=x0,x(tf)=xf.\begin{aligned} t_0 &= \text{fixed},\\ t_f &= \text{fixed},\\ x(t_0)&=x_0,\\ x(t_f)&=x_f. \end{aligned}

The admissible variations satisfy

δx(t0)=0,δx(tf)=0.\delta x(t_0)=0, \qquad \delta x(t_f)=0.

A first extension is the free-final-time problem:

t0=fixed,x(t0)=x0,x(tf)=xf,tf is free.\begin{aligned} t_0 &= \text{fixed},\\ x(t_0)&=x_0,\\ x(t_f)&=x_f,\\ t_f & \text{ is free}. \end{aligned}

Here, the terminal state value is prescribed, but the time at which that value is reached is unknown.

Geometric Interpretation of a Free Final Time

Suppose an admissible path begins at the fixed point

(t0,x0)(t_0,x_0)

and must reach the prescribed terminal value xfx_f. If tft_f is free, the endpoint may occur anywhere along the horizontal line

x=xf.x=x_f.
When the terminal state is fixed but the terminal time is free, admissible trajectories may reach the line x=x_f at different times.

Figure 1:When the terminal state is fixed but the terminal time is free, admissible trajectories may reach the line x=xfx=x_f at different times.

Variation of an Endpoint

The endpoint of a trajectory can change in two different ways:

  1. the function itself may vary at the original endpoint time;

  2. the endpoint time may move.

These effects must be distinguished carefully.

Let the nominal final endpoint be

(tf,xf),xf=x(tf).\bigl(t_f,x_f\bigr), \qquad x_f=x(t_f).

After perturbation, let the endpoint become

(tf+δtf,  xf+δxf).\bigl(t_f+\delta t_f,\;x_f+\delta x_f\bigr).

The perturbed function is

x~(t)=x(t)+δx(t).\widetilde{x}(t) = x(t)+\delta x(t).

A first-order Taylor expansion gives

xf+δxf=x~(tf+δtf)=x(tf+δtf)+δx(tf+δtf)x(tf)+x˙(tf)δtf+δx(tf).\begin{aligned} x_f+\delta x_f &= \widetilde{x}(t_f+\delta t_f) \nonumber\\ &= x(t_f+\delta t_f) + \delta x(t_f+\delta t_f) \nonumber\\ &\approx x(t_f) + \dot{x}(t_f)\delta t_f + \delta x(t_f). \end{aligned}

Since x(tf)=xfx(t_f)=x_f, we obtain the endpoint variation identity

δxf=δx(tf)+x˙(tf)δtf.\boxed{ \delta x_f = \delta x(t_f) + \dot{x}(t_f)\delta t_f }.

Similarly, at the initial endpoint,

δx0=δx(t0)+x˙(t0)δt0.\boxed{ \delta x_0 = \delta x(t_0) + \dot{x}(t_0)\delta t_0 }.

Endpoint Parameters Versus the Dummy Variable

The symbols t0t_0 and tft_f denote endpoint parameters. They are real numbers defining the integration interval. They are not the same object as the dummy integration variable tt.

Thus, while the dummy variable is integrated out, the endpoint parameters may remain as arguments of the resulting functional:

J=J[x(),t0,tf].J=J[x(\cdot),t_0,t_f].

If tft_f is free, then

δtf\delta t_f

is an admissible independent endpoint variation.

Variation of an Integral with a Moving Upper Limit

Consider

J[x,tf]=t0tfL(x,x˙,t)dt,J[x,t_f] = \int_{t_0}^{t_f} L(x,\dot{x},t)\,\mathrm{d} t,

where t0t_0 is fixed but tft_f may vary.

The variation has two contributions:

δJ=t0tfδLdt+L(x(tf),x˙(tf),tf)δtf.\delta J = \int_{t_0}^{t_f}\delta L\,\mathrm{d} t + L\bigl(x(t_f),\dot{x}(t_f),t_f\bigr)\delta t_f.

The second term follows from the fundamental theorem of calculus:

tft0tfL(x,x˙,t)dt=L(x(tf),x˙(tf),tf).\frac{\partial}{\partial t_f} \int_{t_0}^{t_f} L(x,\dot{x},t)\,\mathrm{d} t = L\bigl(x(t_f),\dot{x}(t_f),t_f\bigr).

First Variation for a Free Final Time

To first order,

δL=Lxδx+Lx˙δx˙.\delta L = \frac{\partial L}{\partial x}\delta x + \frac{\partial L}{\partial\dot{x}}\delta\dot{x}.

Therefore,

δJ=t0tf[Lxδx+Lx˙δx˙]dt+L(tf)δtf,\begin{aligned} \delta J &= \int_{t_0}^{t_f} \left[ \frac{\partial L}{\partial x}\delta x + \frac{\partial L}{\partial\dot{x}}\delta\dot{x} \right]\,\mathrm{d} t \nonumber\\ &\quad + L(t_f)\delta t_f, \end{aligned}

where

L(tf):=L(x(tf),x˙(tf),tf).L(t_f) := L\bigl(x(t_f),\dot{x}(t_f),t_f\bigr).

Because

δx˙=ddtδx,\delta\dot{x} = \frac{\,\mathrm{d}}{\,\mathrm{d} t}\delta x,

integration by parts gives

δJ=t0tf[Lxddt(Lx˙)]δxdt+Lx˙δxt0tf+L(tf)δtf.\begin{aligned} \delta J &= \int_{t_0}^{t_f} \left[ \frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) \right] \delta x\,\mathrm{d} t \nonumber\\ &\quad + \left. \frac{\partial L}{\partial\dot{x}} \delta x \right|_{t_0}^{t_f} + L(t_f)\delta t_f. \end{aligned}

The initial endpoint is fixed, so

δx(t0)=0.\delta x(t_0)=0.

Hence,

δJ=t0tf[Lxddt(Lx˙)]δxdt+Lx˙tfδx(tf)+L(tf)δtf.\begin{aligned} \delta J &= \int_{t_0}^{t_f} \left[ \frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) \right] \delta x\,\mathrm{d} t \nonumber\\ &\quad + \left. \frac{\partial L}{\partial\dot{x}} \right|_{t_f} \delta x(t_f) + L(t_f)\delta t_f. \end{aligned}

Using the Fixed Terminal State

The terminal state value is fixed:

xf=constant.x_f=\text{constant}.

Therefore,

δxf=0.\delta x_f=0.

Using the endpoint identity (11),

0=δx(tf)+x˙(tf)δtf,0 = \delta x(t_f) + \dot{x}(t_f)\delta t_f,

so

δx(tf)=x˙(tf)δtf.\boxed{ \delta x(t_f) = - \dot{x}(t_f)\delta t_f }.

Substituting into (25) gives

δJ=t0tf[Lxddt(Lx˙)]δxdt+[Lx˙Lx˙]tfδtf.\begin{aligned} \delta J &= \int_{t_0}^{t_f} \left[ \frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) \right] \delta x\,\mathrm{d} t \nonumber\\ &\quad + \left[ L - \dot{x} \frac{\partial L}{\partial\dot{x}} \right]_{t_f} \delta t_f. \end{aligned}

Now the variation has been separated into two independent parts:

Necessary Conditions for the Free-Final-Time Problem

Because δx(t)\delta x(t) and δtf\delta t_f vary independently, both coefficients must vanish.

The interior condition is the Euler–Lagrange equation:

Lxddt(Lx˙)=0.\boxed{ \frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) = 0 }.

The free-final-time transversality condition is

[Lx˙Lx˙]tf=0.\boxed{ \left[ L - \dot{x} \frac{\partial L}{\partial\dot{x}} \right]_{t_f} = 0 }.

Interpretation Through the Hamiltonian-Like Quantity

Define

p=Lx˙.p = \frac{\partial L}{\partial\dot{x}}.

Then define the quantity

H=px˙L.H = p\dot{x}-L.

The free-final-time condition (32) is equivalent to

H(tf)=0.\boxed{ H(t_f)=0 }.

This relation anticipates the Hamiltonian transversality conditions that appear later in optimal control.

Why Independent Variations Matter

Before integration by parts, the first variation contains both δx\delta x and δx˙\delta\dot{x}. These quantities are not independent because

δx˙=ddtδx.\delta\dot{x} = \frac{\,\mathrm{d}}{\,\mathrm{d} t}\delta x.

Integration by parts rewrites the expression entirely in terms of δx\delta x and endpoint variations.

In the free-final-time problem, the endpoint relation

δx(tf)=x˙(tf)δtf\delta x(t_f) = -\dot{x}(t_f)\delta t_f

further removes the dependence between the endpoint state variation and the endpoint time variation.

Only after these dependencies have been eliminated can the coefficients of independent variations be set to zero separately.

Free Final State with Fixed Final Time

Consider a different problem:

t0=fixed,tf=fixed,x(t0)=x0,x(tf) is free.\begin{aligned} t_0&=\text{fixed},\\ t_f&=\text{fixed},\\ x(t_0)&=x_0,\\ x(t_f)&\text{ is free}. \end{aligned}

Since tft_f is fixed,

δtf=0.\delta t_f=0.

The endpoint identity becomes

δxf=δx(tf).\delta x_f = \delta x(t_f).

The first variation is

δJ=t0tf[Lxddt(Lx˙)]δxdt+Lx˙tfδx(tf).\begin{aligned} \delta J &= \int_{t_0}^{t_f} \left[ \frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) \right] \delta x\,\mathrm{d} t \nonumber\\ &\quad + \left. \frac{\partial L}{\partial\dot{x}} \right|_{t_f} \delta x(t_f). \end{aligned}

Because δx(tf)\delta x(t_f) is arbitrary, the natural boundary condition is

Lx˙tf=0.\boxed{ \left. \frac{\partial L}{\partial\dot{x}} \right|_{t_f} = 0 }.

Thus, a free endpoint value generates an additional boundary condition.

Summary of Common Endpoint Cases

Final endpoint dataVariation relationAdditional condition
xfx_f fixed, tft_f fixedδx(tf)=0\delta x(t_f)=0None beyond prescribed boundary data
xfx_f fixed, tft_f freeδx(tf)=x˙(tf)δtf\delta x(t_f)=-\dot{x}(t_f)\delta t_f[Lx˙Lx˙]tf=0\left[L-\dot{x}L_{\dot{x}}\right]_{t_f}=0
xfx_f free, tft_f fixedδtf=0\delta t_f=0, δx(tf)\delta x(t_f) arbitraryLx˙tf=0\left.L_{\dot{x}}\right\rvert_{t_f}=0

Necessary Versus Sufficient Conditions

The Euler–Lagrange and transversality equations are first-order necessary conditions. They identify stationary candidates.

They do not by themselves determine whether a candidate is:

A condition is necessary if every optimizer must satisfy it. A condition is sufficient if satisfying it guarantees optimality. The first-order conditions developed here are generally necessary but not sufficient.

Worked Example: Mass–Spring System

Consider a unit-mass, unit-stiffness spring system. Let the Lagrangian be

L(x,x˙)=TV=12x˙212x2.L(x,\dot{x}) = T-V = \frac{1}{2}\dot{x}^2 - \frac{1}{2}x^2.

The action functional is

J[x]=t0tf(12x˙212x2)dt.J[x] = \int_{t_0}^{t_f} \left( \frac{1}{2}\dot{x}^2 - \frac{1}{2}x^2 \right) \,\mathrm{d} t.

The required derivatives are

Lx=x,\frac{\partial L}{\partial x} = -x,

and

Lx˙=x˙.\frac{\partial L}{\partial\dot{x}} = \dot{x}.

Therefore,

ddt(Lx˙)=x¨.\frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) = \ddot{x}.

The Euler–Lagrange equation becomes

xx¨=0,-x-\ddot{x}=0,

or

x¨+x=0.\boxed{ \ddot{x}+x=0 }.

The general solution is

x(t)=Acost+Bsint.x(t) = A\cos t + B\sin t.

For prescribed endpoint values

x(t0)=x0,x(tf)=xf,x(t_0)=x_0, \qquad x(t_f)=x_f,

the constants AA and BB are determined from

Acost0+Bsint0=x0,Acostf+Bsintf=xf.\begin{aligned} A\cos t_0+B\sin t_0&=x_0,\\ A\cos t_f+B\sin t_f&=x_f. \end{aligned}

In matrix form,

[cost0sint0costfsintf][AB]=[x0xf].\begin{bmatrix} \cos t_0 & \sin t_0\\ \cos t_f & \sin t_f \end{bmatrix} \begin{bmatrix} A\\ B \end{bmatrix} = \begin{bmatrix} x_0\\ x_f \end{bmatrix}.

Provided the matrix is nonsingular,

sin(tft0)0,\sin(t_f-t_0)\neq 0,

the constants are uniquely determined.

General Mass–Spring Parameters

For mass mm and spring stiffness kk,

L=12mx˙212kx2.L = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2.

Then

Lx=kx,\frac{\partial L}{\partial x} = -kx,

and

Lx˙=mx˙.\frac{\partial L}{\partial\dot{x}} = m\dot{x}.

The Euler–Lagrange equation gives

kxmx¨=0,-kx-m\ddot{x}=0,

or

mx¨+kx=0.\boxed{ m\ddot{x}+kx=0 }.

This is exactly the familiar equation of motion for an undamped mass–spring oscillator.

Free-Final-Time Condition for the Mass–Spring Lagrangian

For

L=12mx˙212kx2,L = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2,

we have

Lx˙=mx˙.\frac{\partial L}{\partial\dot{x}} = m\dot{x}.

The free-final-time condition is

0=[Lx˙Lx˙]tf=[12mx˙212kx2mx˙2]tf=[12mx˙212kx2]tf.\begin{aligned} 0 &= \left[ L - \dot{x} \frac{\partial L}{\partial\dot{x}} \right]_{t_f} \nonumber\\ &= \left[ \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2 - m\dot{x}^2 \right]_{t_f} \nonumber\\ &= \left[ -\frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2 \right]_{t_f}. \end{aligned}

Thus,

mx˙2(tf)+kx2(tf)=0.m\dot{x}^2(t_f)+kx^2(t_f)=0.

For positive mm and kk, this implies

x(tf)=0,x˙(tf)=0.x(t_f)=0, \qquad \dot{x}(t_f)=0.

This shows that a free-final-time condition can strongly restrict the terminal behavior. It also illustrates why endpoint conditions must be interpreted together with the full problem statement.

Extension to Vector-Valued Functions

Let

x(t)Rn\boldsymbol{x}(t)\in\mathbb{R}^n

and consider

J[x]=t0tfL(x,x˙,t)dt.J[\boldsymbol{x}] = \int_{t_0}^{t_f} L(\boldsymbol{x},\dot{\boldsymbol{x}},t)\,\mathrm{d} t.

The endpoint variation identity becomes

δxf=δx(tf)+x˙(tf)δtf.\delta\boldsymbol{x}_f = \delta\boldsymbol{x}(t_f) + \dot{\boldsymbol{x}}(t_f)\delta t_f.

The Euler–Lagrange equations are

Lxddt(Lx˙)=0.\frac{\partial L}{\partial\boldsymbol{x}} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{\boldsymbol{x}}} \right) = \boldsymbol{0}.

For fixed terminal state and free terminal time, the transversality condition becomes

[Lx˙Lx˙]tf=0.\boxed{ \left[ L - \dot{\boldsymbol{x}}^\top \frac{\partial L}{\partial\dot{\boldsymbol{x}}} \right]_{t_f} = 0 }.

General Endpoint Curves

A still more general endpoint may be constrained to lie on a curve

ψ(xf,tf)=0.\psi\bigl(x_f,t_f\bigr)=0.

Then

ψxδxf+ψtδtf=0.\psi_x\,\delta x_f + \psi_t\,\delta t_f = 0.

This relation links δxf\delta x_f and δtf\delta t_f. Substitution into the general boundary variation produces a transversality condition associated with the endpoint curve.

This geometric viewpoint is fundamental in advanced optimal control, where terminal states may be constrained to manifolds rather than prescribed exactly.

Chapter Summary

The main results are:

  1. The functional

J[x]=t0tfL(x,x˙,t)dtJ[x] = \int_{t_0}^{t_f} L(x,\dot{x},t)\,\mathrm{d} t

depends on the entire path and may also depend on free endpoint parameters.

  1. The dummy integration variable tt is not the same object as the endpoint parameters t0t_0 and tft_f.

  2. The total final-state variation satisfies

δxf=δx(tf)+x˙(tf)δtf.\delta x_f = \delta x(t_f) + \dot{x}(t_f)\delta t_f.
  1. If xfx_f is fixed and tft_f is free, then

δx(tf)=x˙(tf)δtf.\delta x(t_f) = -\dot{x}(t_f)\delta t_f.
  1. A moving upper integration limit contributes

L(tf)δtfL(t_f)\delta t_f

to the first variation.

  1. The interior stationarity condition remains the Euler–Lagrange equation:

Lxddt(Lx˙)=0.\frac{\partial L}{\partial x} - \frac{\,\mathrm{d}}{\,\mathrm{d} t} \left( \frac{\partial L}{\partial\dot{x}} \right) = 0.
  1. For fixed xfx_f and free tft_f, the additional transversality condition is

[Lx˙Lx˙]tf=0.\left[ L - \dot{x} \frac{\partial L}{\partial\dot{x}} \right]_{t_f} = 0.
  1. For free xfx_f and fixed tft_f, the natural boundary condition is

Lx˙tf=0.\left. \frac{\partial L}{\partial\dot{x}} \right|_{t_f} = 0.
  1. These conditions are necessary, not generally sufficient.

  2. For a mass–spring Lagrangian, the Euler–Lagrange equation reproduces the physical equation of motion.

Connection. The transversality calculation becomes more useful when organized as natural boundary conditions and embedded in the standard structure of an optimal-control problem.