optimization—theory and applications || examples and exercises on classical problems

43
CHAPTER 3 Examples and Exercises on Classical Problems 3.1 An Introductory Example Let us briefly consider, in terms of the integral 1 = (l + X'2)1/2 dt, the question of the path of minimum length between two fixed points, or between a fixed point and a given curve in the tx-plane. Here fo = (1 + X'2)1/2 depends on x' only, and satisfies condition (S) of Section 2.8. Any extremal must satisfy Euler's equation in the reduced form (2.2.10), fo,,' = C, or x'(1 + X,2)-1/2 = c. Here fo,,' = x'(1 + X,2)-1/2 is a strictly increasing function of x' with range ( -1,1). Thus -1 < C < 1, and there is one and only one value x' = m, depending on c, such that x'(1 + X,2)-1/2 = c. From (2.6.iii) we know that any optimal solution is of class C 2 and therefore an extremal. Thus an optimal solution must be a segment x' = m, x(t) = mt + b, m, b constants. (a) Fixed end points. For both end points fixed, 1 = (t1o Xl)' 2 = (t2' X2), tl < t2, then x(t) = Xl + m(t - td, tl t t2, m = (X2 - Xl)(t2 - t1)-1, corresponding to the seg- ment s = 12 of length L = «t2 - tl)2 + (X2 - Xl)2)1/2. That this is optimal can be seen as follows. Since fo"',,, = (1 + x') - 3/2 > 0, we have, by Taylor's formula, for any x' and where m is some number between m and x'. Thus, for any AC arc C:x = x(t), tl t t2, x(t 1) = X10 x(t 2) = X2' we have and equality holds if and only if x'(t) = m in [tl' t2] (a.e.), that is, if C is the segment s = 12. Another way to show that s = 12 is the optimal path, is to imbed s = 12 in the field of extremals made up of all straight lines parallel to s and covering the whole tx-plane. Here fo"',,, > 0; hence E(x', X') > 0 for all X' ¥- x', and s = 12 is the path of minimum length between 1 and 2 by Weierstrass's statement (2.1l.iv), so I[C 12 ] I[s] 116 L. Cesari, Optimization—Theory and Applications © Springer-Verlag New York Inc. 1983

Upload: lamberto

Post on 08-Dec-2016

214 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

CHAPTER 3

Examples and Exercises on Classical Problems

3.1 An Introductory Example

Let us briefly consider, in terms of the integral 1 = J:~ (l + X'2)1/2 dt, the question of the path of minimum length between two fixed points, or between a fixed point and a given curve in the tx-plane. Here fo = (1 + X'2)1/2 depends on x' only, and satisfies condition (S) of Section 2.8. Any extremal must satisfy Euler's equation in the reduced form (2.2.10), fo,,' = C, or x'(1 + X,2)-1/2 = c. Here fo,,' = x'(1 + X,2)-1/2 is a strictly increasing function of x' with range ( -1,1). Thus -1 < C < 1, and there is one and only one value x' = m, depending on c, such that x'(1 + X,2)-1/2 = c. From (2.6.iii) we know that any optimal solution is of class C2 and therefore an extremal. Thus an optimal solution must be a segment x' = m, x(t) = mt + b, m, b constants.

(a) Fixed end points. For both end points fixed, 1 = (t1o Xl)' 2 = (t2' X2), tl < t2, then x(t) = Xl + m(t - td, tl ~ t ~ t2, m = (X2 - Xl)(t2 - t1)-1, corresponding to the seg­ment s = 12 of length L = «t2 - tl)2 + (X2 - Xl)2)1/2. That this is optimal can be seen as follows. Since fo"',,, = (1 + x') - 3/2 > 0, we have, by Taylor's formula,

for any x' and where m is some number between m and x'. Thus, for any AC arc C:x =

x(t), tl ~ t ~ t2, x(t1) = X10 x(t2) = X2' we have

and equality holds if and only if x'(t) = m in [tl' t2] (a.e.), that is, if C is the segment s = 12.

Another way to show that s = 12 is the optimal path, is to imbed s = 12 in the field of extremals made up of all straight lines parallel to s and covering the whole tx-plane. Here fo"',,, > 0; hence E(x', X') > 0 for all X' ¥- x', and s = 12 is the path of minimum length between 1 and 2 by Weierstrass's statement (2.1l.iv), so I[C12] ~ I[s]

116 L. Cesari, Optimization—Theory and Applications© Springer-Verlag New York Inc. 1983

Page 2: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.2 Geodesics 117

for all C12 :x = x(t), t1 :::;; t:::;; t2, x AC, joining 1 and 2, with equality holding only if C12 = s.

(b) Variable end point. In the case of the first end point 1 = (t1,X1) being fixed, and the second end point 2 = (t2,X2) being on any arc r:x = g(t), t':::;; t:::;; t", of class C1, then either t' < t 2 < t", or t 2 = t', or t 2 = t". If t' < t 2 < t", then B' is the linear manifold (dt2, dX2 = g'(t2) dt2), dt2 arbitrary (the tangent line to r at (t2, g(t2)))' Here fox' = x'(1 + X'2)-1/2,fo - X}ox' = (1 + X'2)-1/2, and (2.2.9) yields Ll = dt2 [1 + x'(t2)g'(t2)] =

0, dt2 arbitrary; thus X'(t2)g'(t2) = -1. The segment s = 12 is therefore orthogonal to rat 2.

If t2 = t', then B' is the cone dX2 = g'(t2)dt2 with dt2 ~ 0, and (2.2.i) yields Ll = dt2 [1 + x'(t2)g'(t2)] ~ 0 with dt2 ~ 0, or X'(t2)g'(t2) ~ -1. Thus the angie IX between the oriented tangents T to s = 12 and T' to r is 0 :::;; IX :::;; n12.

If t2 = t", then B' is the cone dX2 = g'(t2)dt2 with dt2 :::;; 0, and (2.2.i) yields Ll = dt2 [1 + x'(t2)g'(t2)] ~ 0 with dt2 :::;; 0, or X'(t2)g'(t2) :::;; -1. The angle IX is nl2 s IX S n.

I

",;: /PI T

~ 2 I ........... T' I

I

:1,(' I , I T

Remark. Critical examples concerning the same length integral I above have been mentioned in Section 1.6, Example 3.

3.2 Geodesics

A. The Equation of the Geodesics

Let us briefly consider the problem of joining two given points 1 and 2 on a surface S by an arc C lying on the surface and having the shortest possible length. Any such arc is called a geodesic on S.

Let us assume that S is given in parametric form

x = x(u, v), y = y(u, v), z = z(u, v), (u,v) E D,

as a map of class C1 from a fixed connected domain in D in the uv-plane. Then on any curve C of class C1 on S the arc length parameter s has differential ds given by

(dS)2 = (dX)2 + (dy)2 + (dZ)2

= P(u, V)(dU)2 + 2Q(u, v) du dv + R(u, V)(dV)2,

Page 3: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

118 Chapter 3 Examples and Exercises on Classical Problems

If we think of the given points 1 and 2 on S as the images of certain points on D, say briefly 1 = (Ui> Vi), 2 = (U2, V2), U i < U2, and C as the image of some path Co in D, say of the form V = v(u), Ui :;;; U:;;; U2, V(Ui) = Vi> V(U2) = V2, (u, v(u)) E D, then we have the problem of the minimum of the functional

I = f.N2 (P(u, v) + 2Q(u, v)v' + R(u, V)V'2)1/2 du,

N,

with v' = dv/du. Here Jo(u, v, v') = (P + 2Qv' + RV'2)1/2, and the Euler equation (2.2.12), Jov = (d/du)Jov" is

(3.2.1) 2- i(Pv + 2Qvv' + RvV'2)(p + 2Qv' + RV'2)-1/2

= (d/du)(Q + Rv')(P + 2Qv' + RV'2)-1/2.

If P, Q, R depend on u only, then the reduced Euler equation (2.2.10), or Jov' = C i ,

becomes

If, in addition, Q = 0, that is, the lines u = constant and v = constant on S are orthog­onal, then we have the equation

(3.2.2)

from which we derive (R2 - CiR)v,2 = cip, and then

(3.2.3)

where P and R are functions of u only.

Remark. The existence of at least one parametric path curve of minimum Jordan length joining two given points of S follows easily from general considerations (Cf. Section 14.1A. For the purpose of framing the problem in the present nonparametric discussion we have assumed that a curve of minimum length on S has a nonparametric representation v = vo(u), Ui :;;; U :;;; U2, and we compare its length only with the length of nonparametric curves on S having the same end points. Existence theorem for nonparametric integrals including the integral I are discussed in Section 14. For a discussion ofthe problem only in terms of parametric integrals see L. Tonelli [I].

B. Geodesics on a Sphere

By using polar coordinates, the sphere S of radius a > 0 has equations

(3.2.4) x = a sin u cos v, y = a sin u sin v, z=acosu,

where 0:;;; u:;;; 1t,0 :;;; V :;;; 21t represent colatitude and longitude respectively (and u = 0 and u = 1t represent the poles). Here

P = x; + y; + z; = a2,

Q = XNXv + Y.Yv + Z.Zv = 0,

R = x; + Y; + z; = a2 sin2 u

Page 4: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.2 Geodesics

depend only on u. Thus, for curves C:v = v(u), Ul ::$; U::$; U2, on S we have

I = a S."2 (1 + V'2 sin2 U)1/2 du, ",

119

fo(u, v') c.epends only on u and v', and the Euler equation fov' = Cb i.e. (3.2.2), becomes v' sin2 u = C 1(1 + V'2 sin2 U)1/2, or

(3.2.5)

Equation (3.2.3) now is

v = C1 f (sin4 u - ci sin2 U)-1/2 du;

hence, v = -arcsin(,.-l cot u) + C2 , or sin(C2 - v) = ,.-1 cot u, and finally

(3.2.6) (sin C2)(a sin u cos v) - (cos C2)(a sin u sin v) - ,.-l(a cos u) = 0,

where,. is a constant with ,.2 = Cl 2 - 1. In other words, any optimal path is on a plane Ax + By + Cz = 0 through (0,0,0). Thus, any geodesic on a sphere S is an arc of a great circle.

A few details on the integration above should be mentioned. First we have taken points 1 = (Ul' vd, 2 = (U2' V2) not at antipodes on S, and we have chosen the polar representation in such a way that 0 < Ul < U2 < n, -n12 < Vl, V2 < nl2 (or even Vl = V2 = 0). Thus, we seek optimal paths oftheform C:v = V(U),Ul::$; U::$; U2,V(Ul) = Vl'

V(U2) = V2' Now we note that C1 = 0 in equation (3.2.5) corresponds to the solutions v' = 0, or v = constant, the arcs of meridians on S. For C1 #- 0, then we must have IC11 < 1, and actually, if q = min [sin Ub sin U2], then sin U ~ q > 0 as U describes [Ub U2]. From equation (3.2.5) we derive now IC11::$; q < 1, and Cl 2 ~ q-2 > 1,,.2 ~ q-2 _ 1 > O. Moreover, forul < U < U2, we have sin U > q > 0, 1 + cotl U = sin- 2 U < q-2, cot2 U < q-2 - 1 ::$; ,.2, and ,.-llcot ul < 1. Thus, arcsin(,.-l cot u) is defined in (Ub U2) with values in (0, n), and

(dldu) arcsin(,.-l cot u) = -(1 _,.-2 cot2 U)-1/2,.-1 csc2 U

= _(,.2 + 1 - sin- 2 U)-1/2 sin- 2 U

= - C1(sin4 U - Ci sin2 U)-1/2,

from which (3.2.6) follows. We have proved that any geodesic between two points 1 and 2 on S not at antipodes is the smaller of the two arcs of the great circle which joins them. This statement can be extended to points 1 and 2 at antipodes on S by the remark that any subarc also must be a geodesic, but now there are infinitely-many such arcs of great circles joining 1 and 2.

Since we know from existence theorems (cf. Section 14.1A) that paths of minimum length between any two points on S exist, we have also proved that the arcs of great circles mentioned above are the geodesics on S.

Note that, without the use of existence theorems, one can obtain the same results from sufficient conditions and the theory of fields. For instance, for points 1 and 2 not at antipodes, it is not restrictive to assume 1 = (Ul,O), 2 = (U2' 0), 0 < Ul < U2 < n. We take C 12 :v = 0, Ul ::$; U::$; U2, and then C12 is imbedded in the field of extremals v = constant (the family of meridians). Finally, Cll is optimal, as follows from Section 2.11, sincefov'v' = (1 + v'2 sin2u)-3/2sin2u > 0 and E(u, v', V') > 0 for all V' #- v'.

Page 5: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

120 Chapter 3 Examples and Exercises on Classical Problems

3.3 Exercises

1. Study the geodesics on a surface of revolution y2 + Z2 = [g(x)J2, or in parametric form

(3.3.1) x=u, y = g(u) cos v, z = g(u) sin v.

2. Show that the family of geodesics of the paraboloid of revolution

x=u, y = UI /2 cos V, z = UI /2 sin v has the form

u - C2 = u(l + 4C2) cos- 2 {v + 2C log[k(2(u - C2)1/2 + (4u + 1)1/2)]}.

where C and k are arbitrary constants. 3. Prove that any geodesic on one nappe of the right circular cone

(3.3.2)

has the following property: If the nappe is cut from the vertex along a generator and the surface of the cone is made to lie flat on a plane surface, then the geodesic becomes a straight line. Hint: use the following representation of the cone:

in terms of the parameters r and e. Then rand e become usual polar coordinates on the flattened surface of the cone.

4. Prove the property analogous to the one in Exercise 3 for geodesics on a right circular cylinder.

5. Prove the same for an arbitrary cylindrical surface. 6. Show that for integrals 1 = J:~ fo(t,x')dt, with fo of class C2 in x' and fox'x' ~ 0, the

function 1(8) = J:~ fo(t, X') dt is of class C2 in 8 if X = x(t) + 811(t), where I1(tI) = I1(t2) = 0, and x, 11 are continuous with sectionally continuous derivative. Moreover, 1;~) = J:~ fox'x,(t, x')112 dt. Finally, show that if x satisfies the Euler equation for 1, then x is an absolute minimum for 1.

7. Show that the function fo for the geodesics on a surface of revolution satisfies the property of Exercise 6.

3.4 Fermat's Principle

Fermat expressed the principle that the time elapsed in the passage of light from a source A to an observer B should be the minimum possible. Thus, in an optical homog­eneous material where the velocity is constant, the light path is the segment AB.

(a) If A and B are in two different media, each optically homogeneous with light velocities VI and V2 respectively, and they are separated by a plane n, then the Fermat principle implies the usual law of refraction, or Snell's law. Indeed, we can choose xyz coordinates so that n is the plane y = 0, and A and B are on the plane z = 0, say A = (Xl> YI), B = (X2' Y2), with Zl = Z2 = 0, Xl < x2, YI < 0, Y2 > 0; then the light path

Page 6: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.4 Fermat's Principle 121

ACB must be in the same plane and made up of two segments AC and CB, C = (x, 0) again with z = O. Then for the time T from A to B we have

T = vll[(x - xd2 +.yi]I/2 + V21[(X2 - X)2 + ynl/2,

and the minimum of T is attained for x satisfying

dTjdx = vll[(x - xd2 + yi]-1/2(X - XI) - V21[(X2 - X)2 + yn- I /2(X2 - x) = 0

or

(3.4.1) sin CPI sin CP2

where CPI is the angle between the normal to the interface y = 0 and the path AC, and CP2 is the correspondent angie for CB. This is the law of refraction.

" A

~------------~r-------~X

B

B

y

(b) If we have now a set of N contiguous parallel faced homogeneous media in which the light velocities are Vb V2, ... , VN , and ifC I , C2, ... , CN - I are the points at which the light path crosses the separation planes y = hi, Y = h2 , ••• , y = hN - I , hi < hz < ... < hN- I, then the same principle will require that the light path ACICz ··· CN_IB be a path of minimum time, and thus a polygonal line with vertices C I , Cz, ... , CN-I. Consequently, each subarc ACICz, CI CZC3 ' ••• , CN-ZCN-IB must be an optimal path between A and Cz, C I and C3 , ••• , CN- Z and B respectively, and from the above we derive that the path ACICz ... CN_IB is contained in a plane orthogonal to the interfaces and satisfies

(3.4.2)

where CPI, cpz, ... , CPN are the angles of the consecutive segments with the normal to the interfaces. In other words, vi I sin CPi = K, i = 1, ... , N, K a constant.

(c) Let us assume now that the medium is such that the light velocity V = v(y) is a continuous positive function. Then the Fermat principle stating that the light path

Page 7: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

122 Chapter 3 Examples and Exercises on Classical Problems

from A to B is a path of minimum time, is reducible here to the minimization of the integral

where y' = dy/dx and we think of I as depending on the admissible trajectory y = y(x), Xl ::;; X ::;; X2' with y(XI) = YI' y(X2) = Y2' A = (Xl> yd, B = (X2' Y2), Xl < X2' Then to = v- l (y)(1 + y,2)1/2, and the DuBois-Reymond equation (2.2.14), or to - Y'!Oy' = C, reduces here, after simplification to

(3.4.3)

r---------------------~x

A

B

y

If c/> denotes the angle of the tangent t to the path at P = (x, y(x» with the normal to the planes y = constant, we have sin c/> =(1 + y'2)-1/2, and hence

(3.4.4) V-I sin c/>(x) = c. This extends the law ofrefraction, as stated by (3.4.1) and (3.4.2) to the continuous case. For instance, for v = a, then (3.4.3) yields y' = m and y(x) = mx + c, a, m, C constants. For instance, for v = ay-I/2, then (3.4.3) reduces to y' = (Ca)-l(y - C2a2)1/2, and then X = 2Ca(y - C2a2)1/2 + D, y > Ca, a, C, D constants, C, a > O.

Note that the problem ofthe minimum time of descent (the classical brachistochrone problem of Section 1.6, Example 1) is a particular case of the problem above. Indeed, we have only to take v = (y - 1X)1/2, Y 2! IX, whence Y2 > YI 2! IX, and then the functional above reduces to the one we have mentioned in Section 1.6, Example 1 (see Section 3.12 for more details). Actually, John Bernoulli in 1696 solved the classical brachistochrone problem as a problem of optics. The existence of the minimum for integrals of the type above is proved in general in Section 14.3, and in detail for the brachistochrone problem in Section 14.4, Example 3.

Exercises

1. Write the integral which must be extremized, according to Fermat's principle, ifthe light paths are not restricted to plane curves, and with v = v(x, y, z). Let X be the independent variable.

2. Describe the plane paths of light in the two dimensional media in which the light velocities are given respectively by (a) v = ay; (b) v = a/y; (c) v = ayl/2, where a > 0, y>O.

Page 8: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.5 The Ramsay Model of Economic Growth 123

3.5 The Ramsay Model of Economic Growth

We consider here an economy in which a single homogeneous good is produced with the aid of a capital K(t) which may depend on t, and in which the total output yet) is either consumed or invested. Thus, if C(t) is the total consumption, we have yet) = C(t) + K'(t). It is assumed that there is no deterioration or depreciation of capital, and that the production yet) is a known function Y = 'l'(K) ofthe capital. We shall require C ~ 0, 'l'(K) ~ 0, while K' may be positive or negative: since there is no deterioration or depreciation, the capital can be consumed.

Since the economic objective of any planning concerns the standard of living, we assume that a utility function U( C) is known which measures the instantaneous well­being of the economy, and we assume that in any planning we should try to maximize the global utility

w = 1"2 U(C(t»dt Jtl

in a finite interval oftime [t 10 t 2] (though we need not exclude t 2 = + 00, infinite horizon). Here C(t) = yet) - K'(t), and Y = 'l'(K). Thus, we have the problem of the calculus of variations concerning the maximum of

w = 1"2 U('l'(K(t» - K'(t»dt. Jtl

Here we assume that U( C) is a smooth, positive, nondecreasing function of C for C ~ 0, and two obvious constraints are here that K(t) ~ ° and K' S; 'l'(K(t». However, neither of the two extreme cases K = ° and K' = 'l'(K) should be taken into consid­eration. Let us see whether an optimal solution can be visualized in the interior of the domain, that is, with K(t) > 0, K(t l ) = KI > ° initial capital, and K'(t) < 'l'(K(t». In this sense, we have a free problem of the calculus of variations with

fo(K,K') = U('l'(K» - K'), fOK' = - U'('l'(K) - K'),

and the DuBois-Reymond equation fo - K'!OK' = c becomes

(3.5.1) U('l'(K) - K') + K'U'('l'(K) - K') = c,

where c is a constant. Here we consider only one particular choice for the function U(C), namely, U(C) =

U* - rt(C - C*)2 for ° S; C S; C*, rt > 0, and U(C) = U* for C ~ C*. (Here it can be said that the utility function U "saturates" at C = C*.) For the case t10 t2 finite, the value of U* is not relevant, and thus we may take U* = 0. Also, we assume that the economy is well below the point of saturation: ° < C(t) < C* for tl S; t S; t2. We shall choose for the production function 'l'(K) a linear one, Y = 'l'(K) = PK, P a positive constant. With these conventions and choices, equation (3.5.1) with C = 'l'(K) - K' = PK - K' and with y = c/rt becomes

- rt(PK - K' - C*)2 - 2rtK'(PK - K' - C*) = c, or

(3.5.2) K,2 - (PK - C*)2 = y.

In the phase plane (K, K'), equation (3.5.2) represents a family of hyperbolas with center K' = 0, K = K* = C* /P, for y oF 0, and corresponding asymptotes K' = ± (PK - C*).

Page 9: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

124 Chapter 3 Examples and Exercises on Classical Problems

K'

---r----~-i-r--~--1_~~------.K

K'

_.:..-C_~

All these curves are traveled in the sense indicated in the picture. The hyperbolas above and below the asymptotes correspond to positive constants y > 0; those on the right and left, to y < O.

Now let us consider specific values Ko < Kh say Ko < Kl ~ K*. There are infinitely many such curves along which the economy could move from Ko to Kl (e.g., the bold arc on trajectory 1 in the second illustration). These differ one from the other, however, in the time T required to go from Ko to K 1• The higher curves have larger values of K' all the way and hence reach Kl in a shorter time than the lower curves. We can label these curves with the value of T, the time required to go from Ko to K 1• Thus, the specification of T determines the value of y. A small T corresponds to a large positive value of y. As T increases, y decreases and finally becomes negative. For T large, y < 0, the curves take the value K twice, at two different times Tl and T2 (e.g., the bold arc on trajectory 3 in the second illustration): K first rises above Kl and then decreases, returning to K 1• Note that, in general, K is positive at t = Th that is, at t = Tl the economy is not consuming all it produces. If Kl = K*, then K' -+ 0 as K -+ K*, and

Page 10: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.6 Two Isoperimetric Problems 125

this occurs for y = ° (bold arc on the trajectory 2). In this situation, K' = fJ(K* - K) and by integration K(t) = K* - D exp( - fJt), D a constant, and this shows that K -+ K* as t -+ + 00. For y > 0, the point of equilibrium K' = ° is reached by an optimal path. Note that at this point dK'/dK = (dK'/dt)(dt/dK) = K"/K', or K" = 0, an inflection point on the trajectories K = K(t) depicted in the third illustration. Five arcs in the phase plane (K, K') have been labeled 1 to 5, and the corresponding arcs of trajectories K(t) are depicted in the third illustration.

K

5

The results above have been derived directly by an analysis of the differential equation (3.5.2) in the phase plane (K, K'). However, the same equation can be integrated formally. Thus, for y = fJK - C* and y = 0, we have the equation y' = ± fJy and the arc of trajectory 2 above is given by

K(t) = K* - (K* - K 1) exp( - fJ(t - t1)),

with K(t1) = Kl and K(+oo) = K*. For y> 0, we take y = K - C*/fJ = K - K*, Y = fJ2 L 2 and (3.5.2) becomes y,2 = fJ2(L 2 + y2). Hence, y = - L sinh fJ(t* - t) and the arc of trajectory 1 above with K(t1) = K 1, K(t*) = K* is given by

K(t) = K* - L sinh fJ(t* - t),

where L = yl/2 p-l and t* is determined by L -l(K* - K 1) = sinh fJ(t* - t 1).

Exercise

1. Discuss Ramsay's model with 'l'(K) = fJK and U(C) = - Uoe-·c, IX constant.

3.6 Two Isoperimetric Problems

1. The classical isoperimetric problem. We are to find the curves C in R2 of given distint end points A, B, of given length L, and such that the area between C and the chord AB is maximum. We may well assume A = (thO), B = (t2' 0), t2 - tl = I> 0,

Page 11: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

126 Chapter 3 Examples and Exercises on Classical Problems

and we want to maximize

I[y] = f.12 ydt

II

with

J[y] = f.12(1 + y'2)dt = L. II

We have the classical isoperimetric problem withfo = y,jl = (1 + y'2)1/2. By Section 4.8 the optimal solutions (if any) are among those satisfying the Euler equation and the other necessary conditions corresponding to F = fo + Afl = Y + A(1 + y'2)1/2, where A is an undetermined constant. Since F does not depend on t, the DuBois-Reymond equation F - y'Fy' = C yields

Hence,

and finally

(y - C)(1 + y'2)1/2 + A = o.

±(y - C)(A2 - (y - C)2)-1/2 dy = dt,

+W - (y - C)2)1/2 = t - D,

(y _ C)2 + (t _ D)2 = A2.

These are all the circles in R2 (for A oF 0) and all the straight lines y = C (for A = 0). The passage trough A, B implies that D = 2- I(t l + t2)' If 1:s; L:S; 1tl/2 then there is one and only one arc y = y(t), tl :s; t:s; t2, from such loci of length L. For L < 1 and L > 1tl/2 the problem has no solutions (in nonparametric form). This problem will be resumed in Section 6.8.

2. The shape of a hanging rope. We are to find the slope of a heavy rope or chain, of length 1 extended between two fixed points A and B. Since in the rest position the center of gravity takes the lowest possible position, the problem is reduced to that of finding the minimum of the static moment I with respect to the t-axis, which is directed horizontally. Thus, if 1 = (tl' YI), 2 = (t2' Y2), tl < t2, and if C:y = y(t), tl :s; t:s; t2, y(tl ) = Yl> y(t2) = Y2' denotes any AC real function, or curve joining 1 to 2, we have the problem of minimizing

with J = /,,2 (1 + y'2)1/2 dt = 1.

JII

We have here a classical isoperimetric problem withfo = y(1 + y'2)1/2,j1 = (1 + y'2)1/2. By Section 4.8 optimal solutions are among those satisfying the Euler equation cor­responding to F = fo + Afl = (y + ),)(1 + y'2)1/2, where A is an undetermined constant. Since F does not depend on t, the DuBois-Reymond equation F - y'Fy' = CI yields

y + A = CI(1 + y'2)1/2.

By taking y'= sinh v, wehaveY+A=CI cosh v,dy/dv=C I sinh v,and dt=(sinh V)-l dy=

CI dv. Hence

t=C I V+C2,

and by eliminating v we obtain

y + A = CI cosh v,

y + A = CI cosh«t - C2)/Cd,

a family of catenaries. The undetermined constants A, Cl> C2 should be determined

Page 12: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.7 More Examples of Classical Problems 127

in terms of t1, Yl' t2, Y2, I. A necessary condition for compatibility is of course that I;;:: «t2 - tl)2 + (Y2 - Yl)2)1/2. This problem will be resumed in Section 6.8.

3. The classical isoperimetric problem in parametric form. Weare to find the closed curves C in R2 of given length L which enclose the maximum area. In other words, we want to maximize

I[x,y] = 2- 1 f:(y'x - x'y)dt

with

J[ x, y] = f: (X,2 + y'2)1/2 dt = L,

in the class of all path curves C:x = x(t), y = y(t), 0 ~ t ~ a, with x(O) = x(a), y(O) = y(a). (Here x, yare AC with X,2 + y,2 =F O. If we take for t the arc length parameter s, then X,2 + y,2 = 1 a.e., and a = L). Both I and J are parametric integrals as mentioned in Remark 8 of Section 2.2 since fo = 2- 1(y'x - x'y) and fl = (X,2 + y'2)1/2 do not depend on the independent variable and are positively homogeneous of degree one in x', y'. By Section 4.8 the optimal solutions (if any) are among those satisfying the Euler equations and the other necessary conditions for F = fo + Afl = r 1(y'x - x'y) + A(X,2 + y'2)1/2 for some constant A. Here we have

Fx'y = -2-1,Fy'x = 2-1, and the function F* relative to F defined in Remark 8 of Section 2.2 is

F* = (y,)-2Fx'x' = -(x'y')-lFx'y' = (x')-2Fy'y' = A(X,2 + y,2)-3/2.

If we take s, the arc length parameter, as independent variable, we have X,2 + y,2 = 1, and the Weierstrass form of the Euler equations is now l/r = - A -1, a constant. Thus any optimal solution has constant curvature, and is therefore a circle in R2.

3.7 More Examples of Classical Problems

A. Fixed End Point Problems

We briefly. consider here a few problems for which we seek the minimum in the class of all trajectories with fixed end points.

1. J:~ X,2 dt, J:~ (1 + X,2) dt. Here fo depends on x' o~y, and also satisfies condition (S). Any optimal solution must satisfy Euler's equation in the reduced form fox' = c, or x' = C, c, C constants. Since fox'x' = 2 > 0, by Section 2.6 we know that any optimal solution is of class C2 and an extremal. (From Section 2.2 we know that no comer is possible). Thus any optimal arc is a segment. For both end points fixed, 1 = (t 1> X d, 2 = (t2,X2), the absolute minimum is known to exist (Section 2.20, Example 1 and Exercise 1). Thus it must be given by the segment s = 12.

2. J:~(tx' + x'2)dt. Here fo does not depend on x and satisfies condition (S) of Section 2.7. Hence any optimal solution satisfies the Euler equation fox' = C, or t + 2x = C. Since fox'x' = 2 > 0, by Section 2.6 any optimal solution is of class C2

and an extremal. By integration we have x(t) = -4- 1(t - a)2 + b, a, b constants. For both end points fixed, 1 = (t1>xd, 2 = (t2,X2), tl < t2, the absolute minimum is

Page 13: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

128 Chapter 3 Examples and Exercises on Classical Problems

known to exist (Section 2.20C, Exercise 2). Thus, it must be given by an arc of one of the above parabolas joining 1 and 2. For instance for 1 = (0,1),2 = (1,0), the absolute minimum is given by x(t) = 16 -1(25 - (2t + 3)Z), 0 :0:; t :0:; 1.

3. s:~ (1 + (1 + tZ)x'l)dt. Here fo does not depend on x and satisfies condition (S) of Section 2.7. Hence, any optimal solution satisfies the Euler equation fox' = C, or 2(1 + tZ)x' = C. Since fox'x' = 2(1 + tZ) > 0, by Section 2.6 any optimal solution is of class CZ and an extremal. By integration x(t) = a arctant t + b, a, b constants. For both end points fixed, 1 = (tj, XI), 2 = (tz, Xz), tl < t2 , the absolute minimum is known to exist (Section 2.20C, Exercise 3). Thus, it must be given by such an arc joining 1 to 2.

4. J:~ t-Ix'z dt, 0 < tl < tz, x(t l ) = XI' X(t2) = x2 • Here fo does not depend on X and satisfies condition (S) of Section 2.7. Hence, any optimal solution satisfies the reduced Euler equation fox' = C, or 2t- Ix' = C. Since fox'x' = 2t- 1 > 0 for tl < t < tz, any optimal solution is of class COO and an extremal. Then x(t) = atZ + b, a, b constants. For both end points fixed, 1 = (tl,X I), 2 = (tz, xz), 0 < tl < tz, the absolute minimum is known to exist (Section 2.20C, Exercise 4). Thus, it must be given by such an arc joining 1 to 2,

5. J:; tl/Zx'z dt, 0 = tl :0:; tz, x(t l ) = Xl> x(tz) = x 2. Here j~ does not depend on X and satisfies (S) of Section 2.7. Hence, any optimal solution satisfies the reduced Euler equation fox' = C, or 2tl /1 X' = C. Since fox'x' = 2tl /z > 0 for t > 0, by Section 2.6 any optimal solution x(t), 0 :0:; t:o:; t2 , must be of class Cl (actually COO) for 0 < t:o:; t2 and actually an extremal in such an half open interval. Then, x(t) is such an extremal in the closed interval [0, t2l Thus, x(t) = at l / z + b, a, b 0:0:; tl < tz, the absolute minimum is known to exist (Section 2.20B, Example 4), Thus, it must be given by such an arc joining 1 to 2.

6. S:; tx'l dt, 0 < tl < tl , x(t d = XI' x(tl ) = xz' The extremals are segments X = C and arcs of X = a + b log t, a, b constants. From Section 1.5, Example 4, we know that this problem with fixed end points 1 = (0,1), 2 = (1,0) has no absolute minimum. Actually, no such arc can join 1 and 2. On the other hand, for the same problem with fixed end points 1 = (t1'XI), 2 = (t2,X 2), 0 < t1 < tz, any optimal arc is of class Coo and an extremal. From Section 2.20C, Exercise 5, we know that an optimal solution exists, and therefore it is given by such an arc joining 1 and 2.

7, S:; x'2(l + X')2 dt. The extremals are segments. There may be corner points, since fox' = 2x' + 6X'2 + 4X'3 is not monotone, See Section 3,11 for more details,

8, J:; XX'2 dt, X ~ O. The extremals are solutions of the DuBois-Reymond equa­tion (2,2.14), fo - X}ox' = C, and are segments X = C and arcs of X = aCt - b)2i3, See Section 3.10 for more details.

9. J:; (Xl + X'Z)I/Z dt. The extremals are solutions of the DuBois-Reymond equation fo - X}ox' = C, or C2X'2 = Xl(Xl - C2), and are arcs of X sin(t - a) = b, a, b constants. For fixed end points 1 = (0,1),2 = (nI2, 1) we can take a = -nI4, b = 2-1/2, and the extremal E 12: x(t) = r 1/2 csc(t + nI4), 0 :0:; t :0:; n12, joins 1 and 2, The family of extre­mals x(t) = b csc(t + nI4), 0 :0:; t :0:; n12, covers simply the strip R = [0, n12] x R as b ranges over R. We have a Weierstrass field in R. Since j~x'x' = 2 > 0, by (2.1 Liv) the extremal E12 gives an absolute minimum, Cf. Section 1.5, Example 5, where the same problem with fixed end points 1 = (0,0), 2 = (1,1) was shown to have no absolute minimum. Indeed no arc X = b csc(t - a) can join 1 and 2.

10. J:~ (x - 1X)-1/2(1 + X'Z)1/2 dt (brachistochrone: minimum time of descent), The extremals are arcs of cycloids. See Section 3.12 for details,

11. S:~ x(l + X'Z)I/2 dt, x ~ 0 (surface of revolution with minimum area). The extre­mals are arcs of a catenary, See Section 3.13 for details,

Page 14: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.7 More Examples of Classical Problems 129

12. J~ (x2 + (t2 - 2t)x') dt, both end points fixed 1 = (0,0),2 = (3,1), in the rectangle A = [0 ~ t ~ 3, 0 ~ x ~ 1]. Here n = 1, fo = x2 + (t2 - 2t)x' is linear in x', and in the terms of Remark 4 of Section 2.12, we have F = x2, G = t2 - 2t, w(t,x) = Gt - Fx = 2t - 2 - 2x. Thus, w = 0 on the line r:x = t -1, and w < 0 above r, and w > 0 below r. The minimum of the integral in A is given by the polygonal line 1342 made up of the segments 13 and 42 on the boundary of A and of the singular arcs 3, and I min = l The polygonal lines 152 and 162 are parametric curves in A on which I = 3, and represent the maximum of I in the class of the parametric curves in A. The same integral has no maximum in the class of the AC curves x = x(t), x(O) = 0, x(3) = 1, since for these curves I can only take values lower than 3 and as close to 3 as we want, but this value 3 cannot be attained.

x

2

4

3 6

13. Show that the integral J~ (x2 + 2tx') dt with constraints x(O) = 0, x(3) = 1, (t, x(t» E A = [0 ~ t ~ 3, 0 ~ x ~ 1], has no minimum and no maximum in the class of the AC curves x = x(t), 0 ~ t ~ 3, under the constraints.

14. J:l [(2- 1X2 + 3- 1X3) + (tx2 + 2t2)x'] dt, both end points fixed 1 = (-1, -1), 2 = (1,1), in the class of the AC curves x = x(t), -1 ~ t ~ 1, in the region A = [(t, x)1 It I + Ixl ~ 2]. Show that the integral has a minimum in this class but no maximum. Compute I min•

15. Determine the maximum and the minimum of I = g (- x - x2 - t2x') dt in the class of all AC curves x = x(t), 0 ~ t ~ 1, joining BI = [t = 0, 0 ~ x ~ 1] to 2 = (2, ~), within the rectangle A = [O::S; t::s; 2,0::s; x::s; 1]. Ans.: Imax = -35/12, Imin = -4.

B. Computation of the Functions E and H

1. J:~ x,2dt. Here E(x',X') = X,2 - X,2 - 2x'(X' - x') = (X' - X,)2 ~ 0, H(x', A) = X,2 + AX', and H as a function of x' has a minimum at x' = -A/2.

2. J:~ x,4dt. Here E(x', Xl) = (X' - X')2(X'2 + 2X'x' + 3X'2)~0, H(x', A) = X,4 + AX', and H as a function of x' has a minimum at x' = ( - A/4)1/3.

3. J:~(X2 + X'2)1/2. Here E(x, x', X') = (x2 + X,2)-1/2[(X2 + X'2)1/2(X2 + X'2)1/2_

x2 - x' X')] ~ 0, H(x, x', A) = (x2 + X'2)1/2 + Ax', and H as a function of x' has a mini­mum at x' = -Alxl(1 - A2)-1/2 for IAI < 1.

C. Variable End Point Problems

1. J:~(1 + x'2)dt, first end point 1 = (tl,XI) fixed, second end point 2 = (t2,X2) on the straight line B2 = [(t,X2)it ~ t l , X2 fixed, X2 # Xl]. Then from Subsection A, Exer­cise 1, any optimal arc is a segment. Now fox' = 2x', fo - x' fox' = 1 - X,2, and the

Page 15: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

130 Chapter 3 Examples and Exercises on Classical Problems

transversality relation in (2.2.i) yields (1 - X,2) dtz = 0, dtz arbitrary. Hence, x'(tz) = ± 1. We see that the segment s = 12 from 1 to Bz at an angle n/4 with B2 is the only element satisfying the necessary conditions. From Section 2.2OC, Exercise 6, we know that the minimum exists. Thus, the minimum is given by the indicated segment s = 12.

2. J:~(1 + x,z)dt, first end point 1 = (0,1), second end point on the segment Bz =

[(t, 0) 13 :::;; t :::;; 4]. Obviously, no segment s = 12 can hit B2 in an interior point at an angle n/4, so, from the above, the second end point in the optimal path is not in the interior of B. Thus s = 12 must be the segment from 1 = (0,1) to either (3,0), or (4,0). For 2 = (3,0) we have x' = -t; hence LI = (1 - x,2)dt2 = (8/9)dt2 ~ 0 for all dt2 ~ 0, as required by the transversality relation in (2.2.i). The point (4,0) does not satisfy the analogous relation. Since from Section 2.2OC, Exercise 6, we know that the minimum exists, it must be given by s = 12 with 2 = (3,0).

3. J:~ (l + x,2)dt, fixed first end point 1 = (0,0), second end point 2 = (t2' x2) on the locus B2 = [(t,x)ltx = 1, t > 0, x > 0]. As in no. 1 the minimum is given by a segment s = 12, or x = mt, and by the transversality relation (1 - x,2)dt2 + 2x'dxz = 0, where dt2, dX2 are computed on Bz, while x' is the slope of the trajectory at 2. The equations xdt + tdx = 0 and (1 - x,2)dt + 2x'dx = 0 yield 2xx' - t(l - X,2) = 0 where x' = x/to Thus, 3x2/t = t with tx = 1, and finally t4 = 3. Since from Section 2.2OC, Exercise 7, we know that the minimum exists, it must be given by the segment s = 12, 1 = (0,0), 2 = (3 1/4,3 -114), or x = 3 -l/Zt, 0 :::;; t :::;; 31/4.

4. J:~(tx' + x'2)dt, first end point 1 = (0, Xl) fixed, second end point on the straight line B2 = [(t2' x) 1 t2 > 0 fixed, x arbitrary]. Here B' is the linear manifold B' = [dt2 = 0, dX2 arbitrary]. Then (2.2.i) yields (t + 2x')dx2 = 0, and x' = -t/2. In other words, the angle IX between the locus 12 and B2 at 2 depends on t2. From Subsection A, Example 2 above, x(t) = -4 -l(t - a)2 + b. Then a and b are determined by the system of equations _4- la2 + b = Xl' -2- 1(t2 - a) = -2- l t2. Thus, a = 0, b = Xl. Since we know from Section 2.2OC, Exercise 2, that the absolute minimum exists, it must be given by x(t) = _4- l t 2 + Xl' 0:::;; t:::;; tz.

5. J:~ (1 + (1 + t2)X,2) dt, first end point 1 = (0, Xl), Xl > 0, second end point on the half straight line B2 = [(t, 0) 1 t ~ 0]. Then (2.2.i) yields [1 - (1 + t2)X'2] dt2 = 0, dt2 arbitrary, or x' = ±(1 + t2}-1/2 at the point 2 = (t2' 0). From Subsection A, Example 3 above, x(t) = a arctan t + b. Thus, b = Xl> and tz, a are determined by the system of equations a arctan t2 + Xl = O,a(1 + t~)-l = ±(1 + ti)-l/z. Weknowthattheabsolute minimum exists. Thus, it must be given by the corresponding arc. For instance, for Xl = 2n/3, we have t2 = 31/2, a = -2, and x(t) = x(t) = -2 arctan t + 2n/3, 0:::;; t:::;; 31/2, is optimal.

6. J:~ t- lx,2dt, first end point (1,0), second end point on the locus B2 = [(t,x)lx = (13/3) + (t - 2)z, t ~ 1]. From Subsection A above, Example 4, we obtained the parabolas X = atZ + b. Thus b = - a for the passage through (1,0), and x(t) = a(t2 - 1), 1 :::;; t:::;; t2· The transversality relation in (2.2.i) yields ( - t -1 X,2) dtz + (2C 1 X') dX2 = 0 with dX2 = 2(t - 2)dt2; hence, x' = 4(t - 2), where x' = 2at. Now a and t2 are determined by the system of equations a(t~ -1) = (13/3)+ (t2 - 2)Z, 2atz = 4(t2 - 2). Thus t2 = -4(a - 2)-1, and by computation we obtain a third degree equation in a, namely 3a3 + 13a2 - 88a + 52 = 0, with three real roots a = t, a = 3.17891, a = - 8.17891. Each ofthe correspond­ing parabolas x(t) = a(t2 - 1) has two real intersections with the parabola B2. For instance, for a = t, x(t) = t(t2 - 1) has two intersections (3,16/3) and (9,160/3) with B2, and at these points x'(3) = 4, x'(9) = 12. Only for t2 = 3 is the equation x'(t2) =

4(t2 - 2) satisfied with the common value 4. For analogous reasons we must discard all the other cases. The absolute minimum is known to exist (Section 2.20C, Exercise 4

Page 16: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.8 Miscellaneous Exercises 131

and remark below). It must be given by x(t) = 1(t2 - 1), 0 ~ t ~ 3. Thus, the only candidate for optimality is E 12 :x = 1(t2 - 1), 0 ~ t ~ 3, with I[E 12] = 8, t2 = 3.

As noticed in Section 2.20, Exercise 4, the problem above has certainly an absolute minimum if we limit 2 = (t2 , x 2) to be on any arc x = (13/3) + (t - 2)2, 1 ~ t ~ N, N finite. Now t2 can be bounded by the following argument. Any minimum better than E12, if any, must be of the form E:x = xo(t) = a(t2 - 1), 1 ~ t ~ t2, with I[xo] ~ 8, a(t~ - 1) E B2 , hence, 2a2t~ ~ 8, and a(t - 1) = (13/3) + (t2 - 2)2. As t2 ..... 00 we derive a ..... 0 from the first relation and a ..... 1 from the second one. Thus, we must have t2 ~ N for some constant N, and an optimal solution exists by (2.20.i). This is given by E12 •

7. J:~(X2 + X'2) 112 dt first end point fixed 1 = (t1,X1), second end point 2 = (t2,X2) on the straight line B2: t = t2, t2 fixed, t2 > t1. We have already computed the extremals in the Subsection A above, Example 9. Now fox' = x'(x2 + X,2)-112, B' is the linear manifold dt2 = 0, dX2 arbitrary, and the transversality relation in (2.2.i) yields x~ dX2 = 0, and X'(t2) = O. The optimal curve C = 12 must be orthogonal to B2 at 2.

8. H~(x - a:)-112(1 + X'2)1/2d~, first end point 1 = (~1,X1)' fixed, second end point 2 = (~2,X2) on the parabola r:~ = ~o - x 2, ~o fixed, ~1 < ~o - xi, Xl > 0 (brachisto­chrone: minimum time of descent from 1 to the parabola n. (Cf. Section 1.5, Example 1, and Subsection A above, Example 10). The extremals are arcs of cycloids. Here B' is the linear manifold B' = [d~2,dx2 = (2X2)d~2' d~2 arbitrary], and the transversality relation in (2.2.i) yields d~2 + x'dx2 = 0, or (dx2/de 2)x'(e2) = -1. We see that the cy­cloid C12 must be orthogonal to r at 2.

3.8 Miscellaneous Exercises

1. (a) Regarding the left-hand member ofthe obvious inequality

112 [g(t) + Gh(t)] 2 dt ~ 0 Jll

as a quadratic function of the parameter G, prove the Schwartz inequality

(3.8.1)

where equality holds if and only if g(t) = Ah(t) for some constant A. (b) Given x(t1) = Xl> x(t2) = X2, and p(t) a known function, p(t) > 0, use (3.8.1) to prove that the absolute minimum ofJ = J:! p2X'2 dt is m = (X2 - X1)2[J:! p-2 dt]-1, and that this minimum is attained if and only if x' = Ap-2 and then A = (X2 - X1)[J:! p- 2dt]-1. (c) Show that x' = Ap-2 is a first integral of the Euler equation associated with the integral J of part (b). Verify that m is the actual value of J when x' = Ap-2.

2. (a) Show that if x(t) satisfies the Euler equation associated with J = J:~ (p2X'2 + q2x 2)dt, where p and q are known functions of t, then J has the value [p2 XX']:!. (b) Show that if x satisfies the Euler equation associated with J, ifl1(t) is an arbitrary continuous function with sectionally continuous derivative, and ifl1(td = I1(t2) = 0, then

Page 17: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

132 Chapter 3 Examples and Exercises on Classical Problems

(c) Show that by replacing x in J with x + 1'/ (thus x and x + 1'/ have the same end values), the value of J is increased by the nonnegative amount

r'z (p21'/,2 + q21'/2)dt. JI! 3. Investigate the extrema of the following functionals:

(a) l[ x] = Sf x'(l + t2x') dt, x(l) = -1, x(2) = 1. (b) l[x] = J~4(4x2 - X,2 + 8x)dt, x(O) = -1, x(n/4) = O.

(c) l[x] = H(t2x,2 + 12x2)dt, x(l) = 1, x(2) = 8.

(d) l[x] = J6(X,2 + x2 + 2xe2t )dt, x(O) = 1, x(l) = (t)e2.

4. Prove that y = bx/a is a weak minimum but not a strong minimum of the functional

l[x] = f: x,3dt,

where x(O) = 0, x(a) = b, a > 0, b > O. 5. Prove that the extrema of the functional

s: 1'/(t, x)(l + X,2)1/2 dt

are always strong minima if 1'/(t, x) ;::: 0 for all t and x. 6. Show that for the same integral as in Exercise 5, if the first end point is fixed and

the second end point is on a curve r, then the transversality relation implies the orthogonality of the optimal path with r.

9. rlz x'(l + t2x')dt. JI!

14. rtz [(1 + X,2)1/2 - tx] dt. JI!

15. rtz (X,2 - 1)(x,2 - 4) dt. Jt!

3.9 The Integral 1= f(X'2 - x2) dt

Here fo = X,2 - x2 is independent of t, fox = - 2x, fox' = 2x', fox'x' = 2 > 0, and any extremum for 1 is a minimum. The Euler equation is x = -(d/dt)x', hence x" + x = 0, and the extremals are of the form x = Cl sin t + C2 cos t, or x = C cos(t - IX), Ct, C2, c, IX

Page 18: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.9 The Integral I = S(X,2 - X2) dt 133

arbitrary constants. Note that by Section 2.7 any minimum for I satisfies the DuBois­Reymond equation 10 - x'fox' = C; hence, X'2 + X2 = c2, and again x = C cos(t - IX). If the end point conditions are x(tJl = XI' x(t2) = X2, tl < t2, t l , Xl> t2, X2 fixed, then we must determine CI, C2 in such a way that CI sin tl + C2 cos tl = XI, CI sin t2 + C2 cos t2 = X2, and the determinant of this linear system is D = sin(tl - t2)' If D #- 0, then

E12 :x(t) = (sin(tl - t2))-I«XI cos t2 - X2 cos tJl sin t + (X2 sin tl - XI sin t2) cos t),

(a) Case 0< t2 - tl < n. Here D #- 0, and then the extremal E 12 :x = xo(t), tl ~ t ~ t2, through 1 = (t1,XI) and 2 = (t2,X2) is uniquely determined, A field ofextremals containing E12 is then defined by Ea:x = x(t,a) = xo(t) + a sin(t - tl + m) for some fixed m > ° chosen so that t2 - tl + m < n. Then Xa = sin(t - tl + m), and as t describes tl ~ t ~ t2, the argument t - tl + m varies between m and t2 - tl + m, with ° < m < t2 - tl + m < n. Thus Xa has a positive minimum Xa = sin(t - tl + m) ~ J1 > ° in [tl , t2]. Thus, for every t, x(t, a) is a strictly increasing function of a, and the extremals Ea:x = x(t,a), tl ~ t ~ t2, -00 < a < +00, cover simply the strip R = [t l ~ t ~ t2, - 00 < X < + 00]. In addition a(t, x) = [x - xo(t)](sin(t - tl + mW I and p(t,x) = x'(t,a) = x~(t) + a(t,x) cos(t - tl + m) are single valued continuous functions in R (possessing continuous derivatives of all orders). Since Ix'x' = 2 > 0, we have E 12(t, x, p,x') > ° for all (t,x) E R, p, x' reals, and by (2.l1.iv) I[C 12] > I[E 12] for every curve C 12 :x = x(t), tl ~ t ~ t2, x(t) AC in [t l,t2], x'(t) essentially bounded, having the same end points as E 12 • Thus, for 0< t2 - tl < n, the problem has one and only one optimal solution E\2' For instance, for 1 = (0,0),2 = (n/3,/3/2), the optimal solution is E12 :

X = xo(t) = sin t, ° ~ t ~ n/3, and I[EnJ = S0'3 (cos2 t - sin2 t)dt = /3/4 = 0.434. For instance, for 1 = (- 15, Xo - m(5), 2 = (15, Xo + m(5), m #- 0, 0< 15 < n/2, the optimal solution is

E 12 :x = m(215/sin 2(5) cos 15 sin t + xo(2 sin 15/sin 2(5) cos t, -15 ~ t ~ 15,

whose derivative is

x'(t) = m(215/sin 2(;) cos (; cos t - xo(2 sin (;/sin 2(;) sin t, -(; ~ t~ (;,

and we see that for m #- ° and Xo in absolute value below a given constant, then Ix'(t) - ml ~ Ki5 for some constant K; in other words, x'(t) can be made as close to m as we want by taking 15 sufficiently small.

(b) Case t2 - tl > n. The system of equations determining CI and C2 mayor may not have solutions, but even ifit has a solution and an extremal E I2 :X = xo(t), tl ~ t ~ t2, joining 1 and 2 exists, I[x] has no minimum, since the Jacobi necessary condition (Section 2.5) is not satisfied. Indeed, Ixx = - 2, Ix'x' = 2, lxx' = 0, and the accessory equation (relative to Ed is r( + '1 = 0. The solution '1 = sin(t - tl)' tl ~ t ~ t2 , is zero at t = tl and at t = tl + n < t2. Thus, the point 3 on E12 of abscissa t3 = tl + n is conjugate to 1 on E2 and is between 1 and 2. For X2 - XI > n, the integral I[x] has no minimum.

For instance, for tl = 0, XI = 0, t2 = a, X2 = 0, n < a < 2n, the only solution of the Euler equation with these boundary data is xo(t) = 0, ° ~ t ~ a, with I[ xo] = 0, This is no minimum, since for x(t) = sin(na-It), ° ~ t ~ a, we have I[ x] = (2a)-I(n2 - a2) < 0.

(c) Case t2 - tl = n, X2 #- -Xl' The integral I has no minimum since the system determining CI and C2 has no solution.

(d) Case t2 - tl = n, X2 = -XI' This case is more difficult to discuss. The system determining CI and C2 reduces to CI sin tl + C2 cos tl = XI and has infinitely many

Page 19: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

134 Chapter 3 Examples and Exercises on Classical Problems

solutions. By a translation we may always assume tl = 0 and t2 = n, and then the system reduces to C2 = Xl' Thus, for tl = 0, t2 = n, Xl = -x2 = k, we have infinitely many extremals joining (0, k) and (n, - k):

E:x = k cos t + c sin t, O:s; t :s; n,

where c is an arbitrary constant. By direct computation we have

I[E] = S:[(c cos t - k sin t)2 - (c sin t + k cos t)2]dt = 0

for every value of the constant c. We shall prove that each one of the extremals E gives an absolute minimum of I[ x], that is, we shall prove that I[ C 12] ~ I[ E] = 0 for every C: X = x(t) joining 1 = (0, k) and 2 = (11:, - k).

Let R > Ikl be any given constant and let use denote by Ao the region Ao = [0,11:] X

[ -R,R]. The integral I = J~(X'2 - x2)dt with the constraints x(O) = k, x(1I:) = -k, (t, x(t» E Ao has an absolute minimum C 12: X = x(t),O :s; t :s; 11:, by the existence theorem (2.20.i) (since now fo = X'2 - x2 ~ X'2 - R2 and (')'1) holds). If Ix(t)1 < R for all t E [0,11:], then by (2.6.iii) X is of class COO and hence an extremal, namely one of the extremals E above. If C12 has two points in common with the straight line x = R, say P' = (t',R), P" = (t",R), then the entire segment P'P" belongs to Cll , since x' ~ 0, Ixl :s; R on Cll

and x' = 0, x = R on the segment, and C 12 is optimal. Let us assume that P'P" is the maximal segment (or single point) that C12 has in common with x = R,O < t' :s; t" < 11:. Let us prove that the two arcs of C12 ending at P' and P" are tangent to x = R at these points. The argument is by contradiction. Indeed, if this is not true for P' and t' < t", then for {) > 0 sufficiently small we could take two points P 1 = (t' - {), x(t' - {)) and P2 = (t' + {), R) in such a way that the arc P 1P' of C12 has only the second end point on the boundary of Ao and hence by (2.6.i) is of class C l with m = X'(t' - 0) > 0, and the segment P1P2 has slope mo as close as we want to m/2 if {) is sufficiently small. But then the minimum of the integral J(X'2 - x 2) dt between P 1 and P' is not the arc P1P'P2 of C12 but the arc indicated at the end of part (a) above whose slope can be made as close to mo as we want by taking {) sufficiently small, and therefore completely contained in Ao, a contradiction. Thus the arc P lP' is tangent to x = R at P' and m = 0, and an analogous argument holds for P". Also an analogous argument holds for the case t' = t", as well as for the contacts of C 12 with x = - R. We see that the minimizing curve C12 is made up of arcs c cos(t - (X) tangent to the straight lines x = R and x = - R and of segments of these lines. The first of such arcs, starting from the point (0, k), must be of the form k cos t + c sin t. But there is only one such curve say E in this family tangent to x = R, and such curve is in Ao. By part (a) C 12 must coincide with E. Since R is arbitrary, and the value of the integral on all curves E is the same, we see that the minimum of the integral in A is given by the infinitely many curves E.

(e) (A remark concerning cases (b) and (c». Let us prove that inf I[ x] = - 00 for t2 - tl > 1I:,aswellasfort2 - tl = 1I:,X2 i= -xl.LetO = tl < 11: < t2 < 211:,Xl = X2 = O. Let {) > 0 be any number, 0 < {) :s; (211: - t2)/2, and let Cll:x = x(t), O:s; t:s; t2, be defined by taking x(t) = sin t for 0 :s; t :s; t3, and x(t) = C sin(t2 - t) for t3 :s; t:s; t2, where t3 = t2 - 11: + {), C = (sin t3)/(sin {). Obviously C joins 1 = (0,0) to 2 = (t2'0), and is made up of two arcs of extremals with a comer point at 3 = (t3,X(t3»' On the other hand, C12 = C13 + C32 , and I[C13] has a finite limit as {) -+ O. By direct com­putations we have

I[C32] = 2- 1C2 sin2(t2 - t3) = -cos {) sin2(t2 + {)/sin {),

and I[C32] -+ - 00 as {) -+ 0 + O. Finally I[C 12] -+ - 00 as {) -+ 0 + 0, and inf I[C] in the class of curves joining 1 and 2 is - 00. An analogous proof holds in the other cases.

Page 20: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.10 The Integral I = JXX,2 dt 135

3.10 The Integral 1= f XX' 2 dt

(a) First, let us consider the class Q of all curves C:x = x(t), t1 ::;; t::;; t2, with x(t) ~ 0, x(t) AC in [t1o t2], joining two given points I = (t1' Xl), 2 = (t2, X2), t1 < t2, Xl ~ 0, X2 ~ O. Then A is the strip [t1 ::;; t::;; t2, 0::;; X < +00], and then fo = XX'2 ~ 0, and fo is continuous with its partial derivatives, fox = X,2, fox' = 2xx', fox'x' = 2x ~ 0, everywhere in A. Let Co denote any optimal solution, if any. Since fox'x' > 0 for X > 0, by (2.6.iii) we conclude that any arc of Co lying above the t-axis, is of class Coo and an extremal arc (corners may occur only on the t-axis). The Euler equation yields 2xx" = _X,2, and since x" = (dx'/dx)(dx/dt), also 2xx'(dx'/dx) = _X'2. Thus, either x' = 0 and x(t) = b,or x' *- o and 2x(dx'/dx) = -x', from whichxx'2 = C, Cconstant,C ~ O. Then, either C = 0 and x(t) = 0, or C> 0 and (2/3)X3/2 = ±at + band x(t) = kit - toI2/3, k, to constants, k ~ O. These are all possible extremals.

If I = (t1o c), 2 = (t2' c), C ~ 0, then E12 :x = xo(t) = C, t1 ::;; t::;; t2, is an extremal arc without corners, and obviously optimal, since I[ x] ~ I[ xo] = O.

If 1= (t1o O), 2 = (t2,X2), X2 > 0, then E12 :x = xo(t) = k(t - td2/3 with k =

X2(t2 - t1)-2/3 is the only possible optimal solution, and E12 is of class Coo with the exception of the point t1 where E12 has a vertical tangent. Analogously, if I = (t1o Xl), Xl> 0, 2 = (t2,0), then E 12 :x = x(t) = k(t2 - t)2/3 with k = X1(t2 - t1)-2/3 is the only possible optimal solution, and E12 possesses a vertical tangent at t2.

Ifl = (t1oxd,2 = (t2,X2),X1 > 0,X2 > 0,X1 *- x2,we consider only the case Xl < X2, the other one being analogous. Here there is an extremal of class Coo joining I and 2, namely E 12 :x = xo(t) = k(t - to)2/3, with k and to determined by k(t1 - to)2/3 = Xl' k(t2 - tof/3 = X2' or (t2 - to)/(t1 - to) = (X2/X1)3/2; hence to < t1 < t2. By computa­tion, we find I[E 12] = (~)k3(t2 - t1) and finally I[E 12] = (~)(X~/2 - X~/2)2(t2 - t1)-1. Thus, I[E 12] is a continuous function ofthe end points I and 2 of E12 •

Actually, in the last considered case 0 < Xl < X2 there are other solutions C12 :

x = x(t), made up of two arcs of extremals with a cusp at some point 3 = (0, t3), t1 < t3 < t2, with vertical tangent there. Indeed for any such t3 we can take x(t) = k1(t3 - t)2/3 for t1 ::;; t::;; t3, and x(t) = k2(t - t3)2/3 for t3 ::;; t::;; t2, provided k1 = X1(t3 - td- 2/3, k2 = X2(t 1 - t3)-2/3.

Let us prove that in all cases E12 is the only optimal solution. We begin with the latter case: Xl> 0, X2 > 0, Xl < X2, where we have to < t1 < t2. The family of extremal arcs Ea:x = x(t, a) = alt - toI2/3, t1 ::;; t::;; t2 (the same to as

for Ed describes the region R = A as a describes [0, + (0). The corresponding functions are a = a(t, x) = x(t - to)-2/3, P = p(t, x) = (2/3)a(t - to)-1/3 = (2x/3)(t - to)-1, and since to is outside [t1,t2], both functions a(t,x), p(t,x) are continuous with all their partial derivatives in A. Here E(t,x,p,x') = 2x(x' - p)2 and therefore we have E ~ 0 for all (t, x) E R, E > 0 for x > 0, x' *- p. By the usual formula (2.11.5) we derive ICC] - I[E 12] > 0 for any curve C12 joining I and 2 in R. In particular this is true for the solutions E12 with a cusp at a point (t3'0).

The same reasoning when Xl, X2 are not both positive runs into trouble, but can be modified. Assume, for instance Xl = 0, X2 > O. Then E12 :X = xo(t) = k(t - t1)2I3, t1 ::;; t ::;; t2, and the family Ea:x = x(t, a) = a(t - t1)2/3 defines a field in every region Rd == [t 1 + ij ::;; t ::;; t 2, 0 ::;; X < + 00], but not in A.

Let C12 :x = x(t), t1 ::;; t::;; t2, be any curve in R joining I = (t1' 0) and 2 = (t2,X2), X2 > 0, with x AC and XX'2 L 1-integrable in [t1o t2]. Let us consider a point 3 =

(t1 + ij, X(t1 + ij)) on C12' with ij > 0 as small as we want. Let E32 denote the optimal

Page 21: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

136 Chapter 3 Examples and Exercises on Classical Problems

solution for the two points 3 and 2 as fixed end points. If C12 is distinct from E12 , so are C32 and E32 for b sufficiently small, and I[C 32] > I[E32]. On the other hand, I[C32] -> I[ C d as b -> 0 by the L 1-integrability of XX'2 in [tt. t2]; and also I[Ed ->

I[E12] as b -> 0 by the remark made above concerning the continuity of I[E12] as a function of the end points 1 and 2. Actually, given e > 0 we can take b so small that both differences I[C 32] - I[C 12] and I[E32] - I[E 12] are in absolute value less than e, while the difference I[C32] - I[E 32] is given by formula (2.11.5) and has a positive limit as b -> O. Thus, I[ C 12] > I[ E 12] for any curve C 12 distinct from E 12' The case 1 = (t1,C), 2 = (t2,c), c;::: 0 does not present difficulties since EI2 is the segment x = C

between 1 and 2, and a field is easily constructed by taking all the segments x = a, Os a < +00. Thus I[C12] > I[E12] for any curve C12 in R between 1 and 2.

(b) The same problem above but in the whole strip t1 s t s t2 , - 00 < x < + 00, has no minimum, since inf ICC] = - 00, as already shown in Section 2.11C, Example 2.

3.11 The Integral I = I x'2(1 + X')2 dt

We shall consider the extrema of this integral in the class Q of all curves C:x = x(t), t1 S t S t2 (x(t) AC in [t1, t2]' x' L 4-integrable) joining two given fixed points 1 = (tt.x1),2 = (t2,X2), t1 < t2. Let m denote the number m = (X2 - X1)/(t2 - td.

(a) Here fo = x'2(l + X')2 = X'4 + 2X'3 + X'2 depends on x' only, fo ;::: 0,

fox' = 4X'3 + 6X'2 + 2x' = 4x'(x' + l)(x' + i), fox'x' = 2(6x'2 + 6x' + 1) = 2(x' - m1)(x' - m2),

with mt.m2=2-1(-1±3-1/2), m1 = -0.7887, m2= -0.2113. Thus fox'x' >O for x' < m1 and x' > m2, fox'x' < 0 for m1 < x' < m2' Since fo does not depend on x, condition (S) of Section 2.7 holds, and hence any optimal solution satisfies the Euler equation foAx') = c. Thus, any extremal arc is a segment x' = a, x(t) = at + b, a, b constants. The Weierstrass function E(t, x, x', X') = fo(X') - fo(x') - (X' - x')foAx') is here E = (X' - X')2[(X' + x' + V + 2x'(x' + 1)].

(b) Assume -1 s m s 0, m = (x2 - xd/(t2 - t1)' Since fo = 0 for x' = 0 and x' = - 1, and fo ;::: 0 otherwise, we see that for any polygonal line E* whose sides have slopes o and -1, as well as for any curve E*:x = x(t), t1 S t s t2, with x(t) AC in [t 1, t2]' and x'(t) = 0 or -1 a.e. in [t 1, t2]' we have I[ E*] = 0, while I[ C] > 0 for any other curve. Since - 1 s m s 0, there is always such a curve, or polygonal line, E 12 joining 1 and 2. Indeed, for m = -lor m = 0, ET2 is necessarily the segment 12. For -1 < m < 0, there are infinitely many polygonal lines and curves ET 2 joining 1 and 2. Indeed, it is enough to take 4>(t) = -ion an arbitrary measurable set H c [t 1, t2] with meas H = (- m)(t2 - t1),

Page 22: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.11 The Integral I = Jx'2(1 + X')2 dt 137

cfJ(t) = 0 on the complementary set, and x(t) = Xl + J:, cfJ("C)d"C. Then x'(t) = -1 or o a.e. in [tl' t2], and x(tl ) = Xl' x(t2) = Xl + (-I)(meas H) = Xl + m(t2 - t l ) = X2' All these curves ET 2 give the absolute minimum to I[ X ] in D, since I [ET J = 0, and I[ C 12] > 0 on every other curve C 12 of D.

Note that along ET2 we have x' = -1 or x' = 0, fo( -1) = fo(O) = 0, fo,..{ -1) = fox-{O) = 0, fox'x'( -1) = 2 > 0, fox'x'(O) = 2 > 0, E = X'2(X' + 1)2 ~ 0, for X' # x', and Euler, DuBois-Reymond, Legendre, and Weierstrass conditions are all satisfied.

(c) Assume that ml < m < m2 where m = (X2 - Xl)!(t2 - tl)' Then E12 :x = Xl + m(t - tl)' tl ::;; t::;; t2, is an extremal arc joining 1 and 2, along whichfox'x' has a constant negative value. The accessory equation reduces here to '1" = 0 and any of its solutions which is zero at t = tl and which is not identically zero has no other zeros. Thus, there are no conjugate points on E12. By (2.11.ii), E12 is a weak local maximum for I[x] in D.

If we assume m < ml' or m > m2' then the same argument proves that E12 :X = Xl + m(t - tl)' tl ::;; t::;; t2, gives a weak local minimum for I[x] in D.

This result can be improved. Again assume first ml < m < m2' and consider the field of extrema Is X = a + m(t - tl)' tl ::;; t::;; t2, a real, covering the strip A = [tl ::;; t::;; t2, - 00 < X < + 00], all of constant slope m. By the usual formula (2.11.5), the difference I[C 12] - I[E 12] is expressed as the integral along C12 of E(t, x, m,x') = 2- 1(x' - m)2 x fox'x,(8), where 8 is some number between m and X'. Since m is in (ml,m2), if we assume ml ::;; X' ::;; m2, then 8 is necessarily in (ml,m2).JOx'x,(8) < 0, and E < 0 for x' # m. We conclude that I[C12] < I[E 12] for every curve C12 in D with slope x'(t) satisfying ml ::;; x'(t) ::;; m2 a.e. in [tl' t2l The same argument shows that if m < ml (or m > m2), then I[C 12] > I[E 12] for every curve C12 of D satisfying x'(t)::;; ml (or x'(t) ~ m2) a.e. in [tl' t2l

This result can be further improved. Again assume ml < m < m2' By using the Weierstrass function E = (x' - m)2[(x' + m + 1)2 + 2m(m + 1)], where the expression in brackets reduces to 6m2 + 6m + 1 = r lfox'x' for x' = m, and hence is negative for x' = m with ml < m < m2. On the other hand, E = (x' - m)2(682 + 68 + 1) for some 8 between m and x'. We see that E < 0 for all x' # m of an interval (PI,P2) containing (m l ,m2)' Thus, I[C 12] < I[E 12] for every curve e12 in D, distinct from E12 , joining 1 and 2 and with slope PI ::;; x'(t) ::;; P2' For instance, for m = -t, we have E = (x' + W«x' + w - t), PI = 2- l ( -1 - 21/2) = -1.2071, P2 = 2- 1( -1 + 21/2) = 0.2071, and (Pt.P2) is an interval larger than (mt.m2)' Note that for ml < m < m2, the extremal E12 is a weak local maximum, but not an absolute maximum, nor a strong local maximum.

Assume now -1 < m < mi' Then E > 0 for all - 00 < x' < ml, x' # m, but the expression for E given above shows that actually E > 0 for all - 00 < x' < PI' x' # m, where PI = -m - 1 - (-2m(m + 1»1/2, and this value of PI must be >ml. Then I[C12] > I[E 12] for all curves C12 ofD distinct from E12 with slopes x'(t)::;; PI' Anal­ogously, we prove that for m2 < m < 0 we have I[C 12] > I[E12] for all curves C12 of D distinct from E12 and slope x'(t) ~ P2' where P2 = -m - 1 + (-2m(m + 1»1/2. In either case E 12 is a weak local minimum, but not an absolute minimum, since I[ E12] > 0, while I[ ET 2] = O.

Assume now m < -lor m > O. Then E > 0 has a constant sign (for x' # m); hence I[C12] > I[E 12] for all curves CI2 in D distinct from E12, and hence E12 is an absolute minimum, and this holds also for m = -1 and m = O.

Summarizing, the extremal arc E12 is an absolute minimum if m ::;; -lor m ~ 0; E12 is a local weak minimum if -1 < m < ml or m2 < m < 0; E12 is a local weak maximum ifml < m < m2'

Page 23: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

138 Chapter 3 Examples and Exercises on Classical Problems

(d) We shall now discuss possible solutions with corner points. At any corner point (t,x) let Sl, S2 be the two slopes, left and right. Then Sl and S2 must satisfy Erdman corner conditions or fOx·(Sl) = fox,(S2)' and fo(sd - sdoAsl) = fO(S2) - Sz/OAS2)' By sUbstitution ofthe expressions above we obtain 4si + 6sI + 2s1 = 4s~ + 6s~ + 2s2, and 3s1 + 4si + si = 3si + 4s~ + s~. By algebraic manipulations we find Sl + S2 = -1, Sl S2 = 0, and hence either Sl = -1, S2 = 0, or Sl = 0, S2 = -1. We obtain again the curves E12 we have discussed in (b), where we proved that each one of these curves gives an absolute minimum for I[ x J in A.

There are other solutions. For instance, assume ml < m < m2, and take any two numbers Sl, S2 such that fOx·(Sl) = foAs2)' and Sl < -1 < ml < m < m2 < S2 < O. This is possible since fox'(s) < 0 for S < -1 and for -t < S < 0, and then c = foAsl) = fox·(S2) < O,fox'x·(Sl) > 0,fox'x,(S2) > O. Take any two complementary subsets El, E2 in [t1> t2J with meas El = (t2 - t l )(S2 - m)(s2 - Sl)-l and then meas E2 = (t2 - t l )(m­Sl)(S2 - Sl) -1. Take cf>(t) = Sl for t EEl' cf>(t) = S2 for t E E2, and let E'12 be the trajectory E~2:X = x(t) = Xl + f~l cf>(r)dr. Then x(td = Xl' x(t2) = x 2, X is AC with x'(t) = Sl for t E El (a.e.), x'(t) = S2 for t E E2 (a.e.), and fox.(t) = c for t E [t l , t2J (a.e.), i.e., the Euler equation is satisfied. If C12 :x = X(t), tl ;<;; t;<;; t2, is any other curve with X(tl) = Xl' X(t2) = x2, 1X'(t) - x'(t) I < c5 and c5 = min[ml - s1> S2 - m2J, then by Taylor's formula

I[C 12J - I[E'12J = tUo(X') - fo(x'»dt

= c i t2 (X'(t) - x'(t»dt + 2- 1 i t2 (X'(t) - x'(t»2fox'x,((})dt, Jtl Jtl

where (} is between X'(t) and x'(t). Hence, either (} < ml or (} > m2. Thus, the first integral is zero and the second one is positive. This proves that any of the curves E~ 2 is a weak local minimum for I[x]. Since fO(SI) > 0, fO(S2) > 0, then I[E12] > 0 = I[E!2]. As a particular case we can take El = [t l , t], E2 = [t, t2J, two consecutive intervals of the indicated measures, and then E'12 is a polygonal line made up of two segments of slopes Sl and S2' In the limiting case with Sl = -1, S2 = 0, then we have the curves E!2 of part (b) which are all absolute minima for I[ x].

Again for ml < m < m2, if we take SI, S2 in such a way that -1 < Sl < ml < m < m2 < 0 < 82 and j~x,(sd = j~x'(S2)' the same construction yields curves E12 which are weak local minima for I[ x].

For -1 < m < ml we take -1 < SI < m < ml < m2 < 0 < S2; for m2 < m < 0 we take -1 < SI < ml < m2 < m < 0 < S2, and each of the corresponding curves E12 is a weak local minimum for I.

Note that these weak local minima do not satisfy the DuBois-Reymond condition in the usual form. This is not contradictory. Indeed we have proved this conditions to be necessary in Section 2.2 Remark 3, as a consequence of the Euler equation for n = 1, fo and X of class C2 ; hence even for weak extrema since the Euler equation holds. We have also proved the necessity of the DuBois-Reymond condition in Section 2.4 and in Section 2.7 for strong extrema for n ~ 1, fo of class Cl and X AC (respectively for x' essentially bounded, and for x' unbounded but under condition (S». In Section 2.4 we have proved the DuBois-Reymond condition even for weak extrema for n ~ 1, fo of class C l , and x' continuous. The weak local minima we have encountered above do not fall in any of these situations.

(e) We shall now discuss the generalized solutions for the problem of the absolute minimum of the same integral I = S~~ j~(x(t»dt. As we have seen in Section 1.14, we are

Page 24: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.12 Brachistochrone, or Path of Quickest Descent 139

concerned with the minimum of the functional

in the class Q* of all AC scalar functions x(t), tl ~ t ~ t2, with X(tl) = Xl> x(t2) = x 2, all measurable functions PI(t) ~ 0, P2(t) ~ 0, PI(t) + P2(t) = 1, and all measurable functions ul(t), u2(t). As mentioned in Section 1.14, the problem reduces to the minimum of the integral

1*[ x ] = r'2 F o(x'(t)) dt, J'l

where F 0 is defined by taking F o(x') = fo(x') for x' ~ -1 and x' ~ 0, and by taking F o(x') = ° for -1 ~ x' ~ 0. Here F 0 is of class CI and the Euler equation Fox' = C is the same as before if C > 0, while for C = ° it leaves x' arbitrary, -1 ~ x' ~ 0, If -1 < m < 0, any curve ET2:X = x(t), tl ~ t ~ t2, X AC with X(tl) = XI' x(t2) = x 2, -1 ~ x'(t) ~ 0, gives the absolute minimum for 1* with common value zero. If m ~ -1 and m ~ ° the absolute minimum of 1* is given by the segment E 12 of slope m joining 1 and 2 as before,

The curves E12 of part (b) are obviously curves ET2' One can also say that any of generalized solutions ET 2, with x' taking arbitrary values between -1 and 0, is equivalent to any ofthe curves E12 of part (b) with x' taking only the values -1 and 0,

3.12 Brachistochrone, or Path of Quickest Descent

(a) Given two points 1 and 2 in a vertical plane n, 2 at a level lower than 1 and not on the same vertical line as 1, we are to determine the curve C joining 1 and 2 such that a material point P, starting at 1 with velocity VI, will glide from 1 to 2 along C, under the force of gravity only and without friction, in a minimum time T, Let us take in n a Cartesian system of reference so that 1 = (Xl> yil, 2 = (X2' Y2), Xl < X2, Yl < Y2' Let Q

be the collection of all (nonparametric) curves C:y = y(x), Xl ~ X ~ x2, with Y(Xl) = Yl, Y(X2) = Y2, and y(x) AC in [Xl> x 2]. Let s denote the arc length on C from 1 to 2, ° ~ s ~ L; let t be time, ° ~ t ~ T; and let V be the instantaneous velocity of P. Then V = ds/dt, and T = S~ ds/v = S;; (1 + y'2)1/2 dx/v, where (') = d/dx, On the other hand, the increase in kinetic energy in the interval [0, t] must be equated to the decrease in potential energy; hence (mv2/2) - (mvi!2) = mg(y - y d, where g is the constant of acceleration. Hence, v2 = 2gy - 2gYl + vi, so v2 = 2g(y - ex), with ex = Yl - vi!2g, Finally

(') = d/dx,

and we have to minimize the integral I[ yJ in the class Q,

(b) Let us assume VI > 0, Then y - ex = y - Yt + vi!2g ~ vi!2g > 0, We may take for A the region A = [Xl ~ X ~ X2, Y ~ Y1J, so that A is a closed set and Y ~ Yl implies (y - ex)-1/2 ~ .J2i/Vl' Thus fo = (y - ex)-1/2(1 + y'2)1/2 does not depend on the independent variable x, but only on y and y' = dy/dx, and fo is continuous in A together with all its partial derivatives,

By the existence theorem (14.3,ii) (see Example 3 at the end of Section 14.4), I[yJ has an absolute minimum in Q, any minimizing curve is of class C2 and satisfies Euler's equation.

Page 25: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

140 Chapter 3 Examples and Exercises on Classical Problems

(c) Here we have foy' = y'(y - 1X)-1/2(l + y'2)-1/2, and the DuBois-Reymond equa­tion (2.2.14), or fo - y'foy' = c, yields

(y _ 1X)-1/2(1 + y'2)1/2 _ y'2(y -1X)-1/2(1 + y'2)-1/2 = c,

and after simplification and taking c = l/ffb,

(3.12.1) (y - 1X)(1 + y'2) = 2b.

The introduction of a parameter T by means of the equation y' = -tan(T/2) simplifies computations. We have y' = -sin T(1 + cos T)-I, and 1 + y'2 = COS- 2(T/2). Then equation (3.12.1) yields

y - IX = 2b(1 + y'2)-1 = 2b COS2(T/2) = b(1 + cos T).

Hence, dy/dT = -b sin T and

dx/dT = (dx/dy)(dy/dT) = (1jy')(dy/dT) = b(1 + cos T) = 2b COS2(T/2).

By immediate integration we have now

(3.12.2) x = a + b(T + sin T), y = IX + b(1 + cos T),

where a and b are constants of integration (b > 0). This shows that any minimizing curve is an arc 12 of a cycloid-precisely, the locus of a point fixed on the circumference of a circle of radius b as the circle rolls on the lower side ofthe line y = IX = YI - vi!2g. The value of the integralI[y] on any arc E12 of these cycloids passing through 1 and 2 can be easily obtained by the use of the variable T, and I[E 12] = (2b)I/2(T2 - TI)' We prove below that there is one and only one of these cycloids passing through 1 and 2. Analyti­cally this boundary value problem reduces to the determination of four real numbers a, b, TI, T2 satisfying the four equations

XI - a = b(TI + sin TI ),

YI -IX = b(1 + cos TI ),

x 2 - a = b(T2 + sin T2)

Y2 - IX = b(1 + cos T2).

Xl

.----r----------------------------~X 0' -,+----------------

" --~ //1

~- -" 2 -------

y

(d) Let us prove that there is one and only one cycloid (3.12.2) passing through 1 and 2. Let us draw an arbitrary one of the cycloids (3.12.2) and intersect it by a line 1'2' parallel to the straight line 12 as shown in the illustration. If we move the line 1'2', keeping it always parallel to 12, from the position L' to the position L", then the ratio of the length of the segment 0'1' to that of the segment 1'2' increases from 0 to 00 and

Page 26: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.12 Brachistochrone, or Path of Quickest Descent 141

r-__ ~O ____ ~ __ TOOr' __ ~ ______ ~(~a,~a)~ ____ ~x

2

y

passes once only through the value of the corresponding ratio 01 to 12. At the position which gives the equality of these ratios, the lengths of 0' l' and 1 '2' are not necessarily equal to 01 and 12 respectively. By changing the value of b, however, the cycloid can be expanded or contracted into another similar to itself, having the same center (a, ex), and the new segments, say 0"1" and 1"2", corresponding to 0'1' and 1'2', will have the same ratio as before. By properly choosing the value of b, the segments 0"1" and 1"2" can be made exactly equal to 01 and 12 respectively. Finally, changing the value of a, we can slide the cycloid along the fixed line Y = Yl' and 1" and 2" can be made to coin­cide with 1 and 2 respectively. This shows that there is one and only one cycloid (3.12.2) through 1 and 2.

(e) Concluding the argument in subsections (b-d), we can finally state that the integral I[y] has an absolute minimum in D, and that this absolute minimum is given by the only cycloid (3.12.2) passing through the two points 1 and 2. This cycloid is thus the curve along which the point P slides without friction from 1 to 2 under the sole action of gravity in a minimum time.

(f) We could reach the same statement by means of sufficient conditions and without the use of existence theorems.

If E12 is the particular arc of the only cycloid (3.12.2) passing through 1 and 2, say for a = ao and b = bo, let us keep ao fixed, and let b vary from 0 to + 00. Then the family of cycloids

(3.12.3) x = ao + b(t + sin t),

y = ex + b(1 + cos t),

Page 27: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

142 Chapter 3 Examples and Exercises on Classical Problems

where - n < "C < n,O < b < + 00, fills once the whole half plane R = [ - 00 < x < + 00, y> IX]. Indeed, for the same "C and different values of b, say b' and b", we have points (x', y'), (x", y") with (x' - ao)/(x" - ao) = b'/b" = (y' - IX)/(Y" - IX). This shows that (x', y'), (x", y") are on the same half straight line issuing from (ao, IX). The cycloids of the family are similar figures, and cover simply R. For each point (x, y) E R, there is a well-determined value b = b(x, y) such that the corresponding cycloid passes through (x, y) for a value "C = "C(x, y) of the parameter "C, and with slope p(x, y) given by p(x, y) = - tan("C(x, y)/2). We have to prove that "C(x, y), p(x, y), b(x, y) are continuous functions of (x, y) in R. Now "C is given by the equation

x - ao "C + sin "C

y-IX l+cos"C'

which has the form g(x, y) = 4>("C), and

4>'("C) = cos- 2("C/2)[1 + ("C/2) tan("C/2)].

For -n < "C < n, 4>'("C) > 0; hence 4>("C) is always increasing, and 4>(-n + 0) = -00, 4>(n - 0) = +00. Therefore s = 4>("C) has an inverse function "C = I/I(s), 1/1(-00) = -n, I/I( + (0) = n, and "C = I/I[(x - ao)/(y - IX)]. This shows that "C is a continuous function of (x, y) in R. On the other hand, 4>'("C) > 0 is continuous in [ - n, n]; hence I/I'(s) exists for all s, - 00 < s < + 00, I/I'(s) = 1/4>'("C), thus 1/1 is continuous with its first derivative, and thus "C(x, y) = I/I[(x - ao)/(y - IX)] has continuous first partial derivatives in R. Finally, b = b(x, y) = (y - 1X)/(1 + cos "C), and thus b(x, y) is also continuous in R with its first partial derivatives.

We have shown that R is a field of extremals. Since h'y' > 0 in R, by (2.1l.iv) we have I[C12] ~ I[E12] for every curve C12 in Q joining 1 and 2, and equality holds if and only ifC12 == E12 .

(g) We shall now discuss the case Vl = O. Then IX = Yl> andjo = (y - Yl)-1/2(1 + y2)1/2 has the singular line y = Yl' The conditions of the existence theorem (14.3.ii) (see Example 3 at the end of Section 14.4) are all satisfied, and therefore, I[y] has an ab­solute minimum. On the other hand, the considerations of subsection (d) above can be repeated with obvious changes, and we conclude as we did there that there is one and only one cycloid E 12 through 1 and 2, but now E 12 has a vertical tangent at 1. We conclude as in subsection (e) that I[y] has an absolute minimum in Q given by a cycloid E12 .

To reach the same result via sufficient conditions in the present case Vl = 0 requires more work.

(h) Again, as in Subsection (d), the equations (3.12.1), where now -n::; "C ::; n, o ::; b < + 00, give a family of cycloid which fills the half plane Ro = [ - 00 < x < 00, y ~ Yl]. Nevertheless, the corresponding functions p(x, y), u(x, y), b(x, y) are continuous only for y > Yl' and hence in each closed region Rl = [ - 00 < x < 00, y ~ Yl + 0], 0> O. The arc E12 does not belong to C1 since its slope is + 00 at the point 1; never­theless, I[E12] has still the finite value I[EnJ = Jib("C2 - "Cl)' Let Ko be the collection of all continuous curves C12 :y = y(x), Xl::; X ::; x 2 , passing through 1 and 2, with y(x) > IX in (Xl' x 2], of class C. when restricted to (Xl' x2], and such that I[y] has a finite value (as a generalized Riemann integral of a nonnegative function having at most one point of infinity at x = Xl)' Then E 12 E Ko. If C 12 is any curve in Ko distinct from E12 , take any point, say 3, on C12 . Through 3 there passes a unique cycloid E3 from 3 to the vertical line through (ao, IX). Let 7 denote the analogous point ofthe cycloid E12 , also continued up to the vertical line through (ao, IX). Then the sum I[ C13] + I[E3]

Page 28: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.13 Surface of Revolution of Minimum Area 143

varies continuously as the point 3 moveS from 1 to 2 along C12 , beginning with the value I[E12] + I[E27] and ending with the value I[C 12] + I[E27]. Hence, if we show that this sum does not decrease, we conclude that 1[ C 12] 2 1[ E 12]. Take a second point 4 on C12 near 3, and let 6 be the corresponding point on the vertical line x = ao. By (2.9.5) we have

and hence (I[C 14] + I[E46]) - (I[C 13] + I[E3s]) = I[C 34] + 1*[Ds6] - 1*[C34].

The vertical line DS6 is orthogonal to the cycloids of the field; hence dx + p dy = 0 along DS6 ' On the other hand

1*[Ds6] = fU dx + (dy - pdx)J;,,]

= f[(y - 1X)-1/2(1 + p2)1/2dx + (dy - pdx)p(y - 1X)-1/2(1 + p2)1/2]

= f(y - 1X)-1/2(1 + p2)-1/2(dx + pdy) = O.

(In other words, D is transversal to the extremals of the family.) Finally, by Section 2.11,

I[C34J - 1*[C34J = rx , E(x, Y,p, Y')dx 2 0, JX 3

and equality holds if and only if the arc C3S is an arc of extremal. If every partial arc C 34 of C 12 is an arc of one of the cycloids of the field, then C 12 is such an arc and C12 = E12. Otherwise, I[C 12] > I[E 12]. Thus E12 E Ko and is the absolute minimum for l[y] in Ko. For the previous analysis, see G. A. Bliss [I].

3.13 Surface of Revolution of Minimum Area

(a) We are to determine the curVe of the xy-plane joining two points 1 = (Xl, Yl), 2 = (X2, Y2), Xl < X2, Yl 20, Y2 20, lying entirely in the half plane y 2 0 and such that the surface of revolution generated by rotation about the x-axis has minimum area. Let Q be the collection of all curves C:y = y(x),x 1 ~ X ~ X2, y(x1) = Yl, y(X2) = Y2, y(x) AC in [Xl> X2], y(x) 2 O. We may as well study the integral

Y E Q, (') = d/dx,

Page 29: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

144 Chapter 3 Examples and Exercises on Classical Problems

since the area above is 211: times the value of this integral. Here A .is the strip A =

[Xl~X~X2'0~y~ +00]. (b) Note that if Yl = Y2 = 0, then l[y] has an absolute minimum given by the trivial

solution Y == 0. Let Yh or Y2' or both be positive. Here we have fo = fo(y, y/) = y(1 + y'2)1/2, foy' = yy'(1 + y/2 )-1/2,foy·,· = y(1 + y'2)-3/2 ~ 0. Here fo does not depend on t; hence, by Section 2.7, if Co:x = x(t), tl ~ t ~ t2, is any optimal solution, then any arc of Co with y(x) > ° satisfies the DuBois-Reymond equation. By Section 2.6 any such arc is also of class C2, and is an extremal.

The DuBois-Reymond equation (2.2.14), fo - Y'!Oy' = C, yields, after simplification, y(1 + y'2)-1/2 = b, and by integration

(3.13.1) y = b cosh(b- 1(x - a»,

where a and b are constants of integration, and as usual cosh z = 2- 1(c + e- Z ),

sinh z = 2- 1(c - e- Z ) for any z. Since there are no comer points with y > 0, we con­clude that there can be a minimum if and only if either Yl = Y2 = 0, and the minimum is given by a segment ofthe x-axis, Ol' both Yh Y2 are positive, and the minimum (if any) is given by an arc of a catenary. Assume Yl, Y2 > 0, and let us discuss whether there is a catenary (3.13.1) through 1 and 2. The passage through 1 gives the equation Yl = b cosh(b- 1(Xl - a»; hence, if we take Yl = b cosh rx, we obtain for a the explicit value a = Xl - Ylrx/cosh rx~ Hence, Xl - a = Ylrx/cosh rx, b- 1(Xl - a) = rx, Yl = b cosh rx = b cosh(b- 1(Xl - a», and we obtain the family of all catenaries (3.13.1) through 1 in the form

(3.13.2)

We shall denote by (') the operation of differentiation with respect to x, and by the subscript rx the operation of differentiation with respect to rx. Also, the subscript 1 denotes the value taken by the corresponding variable when x = Xl' By computation, we have

y' = sinh(rx + (x - Xl)Yl 1 cosh rx), Yl = sinh rx,

(3.13.3) I I ( h )-1 (Y Yl) Ya=YY1COS rx X---Xl+-' Y' Yl

The tangent to any of the catenaries (3.13.2) at the point 1 has the equation Y - Yl =

Y/l(X - Xl), and the tangent at any point P = (x, y) to the same catenary has the equation Y - Y = y'(X - x). Their point of intersection (X, Y) has ordinate Y given by

(3.13.4) Y I I (' ')-1 (Y Yl) = Y Yl Y - Yl X --- Xl + - . Y' y'1

Thus, we obtain from (3.13.3) that Ya = Y(y' - Yl)/cosh rx. We shall see that on each catenary E of the family (3.13.2) there is a point P = (x, Y), or P = (x, y(x, rx», at which Ya = 0, and this point P will generate the envelope r of the family (3.13.2), and will indeed be the point of contact of E with r, and thus the conjugate of 1 on E. From the latter formula we see that Ya and Y vanish together.

(c) This remark justifies the following simple construction of the conjugate point P to 1 on a given extremal E through 1. If M is the intersection with the x-axis of the tangent to E at 1, then P is the point of contact of the tangent to E from M.

If a point (x, y) moves from 1 along a catenary E:y(x, rx), then Ya is at first positive and then changes sign when (x, y) passes the point P conjugate to 1 on E. We shall

Page 30: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.13 Surface of Revolution of Minimum Area 145

need now the second derivatives of y(x, ex) at the conjugate points P. Thus we shall differentiate Y' and y« and use consistently the fact that y« = ° at these points, and hence, by (3.13.3) with y« = 0, X - Xl can be replaced by (y/y') - (YllY't). Also, we know that y'l < 0, and that y' > Oat P. We obtain

y~ = (Y'Yi)-1(y2y~) cosh ex < 0, y" = Y12y cosh2 ex < 0,

Y«« = Y1 2Y'Y'l(CJ - G~Y) > 0.

The last evaluation, Y«« > 0, is in agreement with the fact that at the points of contact P, we have y~ = ° and a minimum for Y as a function of ex.

M w-~--------------------~x

We are now in a position to study further the family of catenaries (3.13.2) and the locus of the conjugate point P to 1 on the catenaries of the family. First of all, we have

Yl (x - Xl ) y(x, ex) = -- cosh ex + -- cosh ex cosh ex Yl

( ) Yt cosh [(_ex + x - Xl) cosh ex]

ex X - Xl cosh ex Yl - --+--- cosh ex Yl (ex X - Xl) h ' --+-- cos ex

cosh ex Yl

where ex/cosh ex- 0, [(ex/cosh ex) + (x - Xl)/Yl] cosh ex - + 00 as ex - 00, and finally the last fraction in the formula above - + 00; hence, y(x, ex) - + 00 as ex - + 00 for every fixed x > Xl> and analogously y(x, ex) - + 00 as ex - - 00. We should nQte that Y« changes from negative to positive whenever it vanishes, since Y«« > 0; hence Y« can vanish only once. Thus, for every x > Xl fixed, y(x, ex) diminishes from - 00 to a mini­mum and then increases to + 00 again, since Y« varies from negative values to zero once and then to positive values. Let us denote by g(x) the minimum of y(x, ex) for any x> Xl' The curve G:y = g(x), X > Xl> is the locus of the conjugate points P. Through every point 2 above G there pass exactly two catenaries of the family, say 132, 124; on one of these (132) there is a conjugate point 3 to 1, and on the other (124) there is none. Every point 3 on the curve G is joined to 1 by one and only one catenary of the family for which 3 is the conjugate point. Every point below G is joined to 1 by no catenary of the family (3.13.2). Thus, if the point 2 is below the curve G, the integral l[y] has no minimum in K, since there is no catenary joining 1 and 2. If the point 2

Page 31: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

146 Chapter 3 Examples and Exercises on Classical Problems

is on the curve G, the integral I[y] has no minimum in K either, since 2 is the conjugate to 1 on E 12 and the reasoning of Section 2.9 applies.

(c) For any given x > Xl the value of IX for which y. = 0 will be denoted by IX(X); hence y.(x,lX(x» = 0, and g(x) = y(x,lX(x». Since y •• > 0 at (X,IX(X», by the implicit function theorem we conclude that IX(X) is differentiable and hence y~ + Y •• IX' = 0, and g'(x) = y' + Y.IX' = y' > O. On the other hand, by using the expressions of y" and y~ above we have

g"(x) = y" + y~IX' = yyd3(cosh 1X)2(yU3 - y3y'?)-1 > O.

The two relations g'(x) > O,g"(x) > o show that Gis concave upward, and thatg(x) ..... + 00

as x ..... + 00. To show that also g'(x) ..... + 00 as x ..... + 00, we should note that the slope g'(x) of Gat P is the same as the slope y'(x, IX (x) ) to the catenary (3.13.2) through P, and this slope is given, because of (3.13.3), by

y'(X,IX(X» = sinh [ IX + Yl1(cosh IX)C,: -7)J = sinh[1X - (sinh IX)-l cosh IX + Yll(yjy') cosh IX J.

This number approaches + 00 as IX ..... + 00. We shall now show that g(x) ..... 0, g'(x) ..... 0 as x ..... Xl + O. The vertex (a, b) of each catenary (3.13.2) is above G, and a = Xl -(YtIX/cosh IX), b = yt/cosh IX. Hence, as IX ..... - 00, we see that a ..... Xl> b ..... 0, and finally (a, b) ..... (Xl' 0). This proves that g(x) ..... 0 as x ..... Xl + O. Also, the slope of the segment joining (Xl' 0) to (a, b) is -l/IX, and -l/IX ..... 0 as IX ..... - 00. This proves that also g'(x) ..... 0 as x"'" Xl +0.

Let 2 be a point above the curve G, and let E 12 be the arc of the unique catenary C:y = bo cosh(bo l(X - ao» joining 1 to 2 and containing no point conjugate to 1. The region above G covered by the curves (3.13.2) is no field (since they all pass through 1). To define a field containing E12 we may take a point 0 on the extension C of the catenary E12 at the left of 1 and so close to 1 that the conjugate point 3 to 0 on C, according to the construction above, is still at the right of 2. The tangents to C at 0 and 3 meet at a point 4 = (x4 ,0) of the real axis. The similarity transformation of center 4, say X - X4 = (bo/b)(X - x4), y = (bo/b)Y, transforms E12 into another catenary

Eb:Y = y(x, b) = b cosh b-l(X - X4 + (b/bo)(X4 - ao»

that is, another catenary of the form y = b cosh(b-l(x - a» with parameters a =

X4 - (b/bo)(X4 - ao) and b. If b is thought of as a parameter, we see that the catenaries Eb fill the V -shaped region V between the radii 40 and 43. These catenaries are all tangent to the two radii 40 and 43. Through each point (x, y) E V there passes one and only one catenary E for b = b(x, y) with slope p(x, y). It is left to the reader to prove

o

Page 32: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

y

3.13 Surface of Revolution of Minimum Area 147

that, in V, b(x, y) and p(x, y) are continuous with their first partial derivatives in V. Since in> 0 in V, by (2.11.iv) we conclude that I[C 12] ~ I[E 12] for every curve C12

ofthe class Q in V, and that equality holds if and only if C12 = E12. Thus E12 is a strong local minimum for I[y] in V.

(d) The reader may ask whether the relation I[C 12] > I[E 12], which is valid for all curves C 12 joining 1 and 2 in the V -shaped region V, is still valid in the larger region R = [Xl ~ X ~ X2, 0 ~ Y < + 00]. The answer, as we shall see below, can be negative.

We shall first consider rectifiable continuous path curves C:x = x(t), y = y(t), tl ~ t ~ t2, in the xy-plane, that is, x(t), y(t) AC in [tl' t2]. Then the integral corre­sponding to I[ C] in this more general class of curves is

I[ C] = Se Y(X'2 + y'2)1/2 dt = Se y ds,

where s is the arc length parameter along C (see Chapter 2 for more details). If we denote by C12 now a segment parallel to the y-axis, say [x = Xl> Y2 ~ Y ~ Yl], and by C13 any given regular parametric curve (X(t), Y(t)) having one end point at 1 and the same length L = Yl - Y2 as s, then

I[C 13] - I[C 12] = S: Y ds - S: yds = S: (Y - y)ds ~ 0

since obviously Y(s) ~ y(s) for every s, and hence equality holds above if and only if Y(s) == y(s), C 13 == C 12, 3 = 2. Now let us consider the usual two points 1 = (Xl> yd, 2 = (X2' Y2), Xl < X2, Yl > 0, Y2 > 0, and denote by P 1342 the polygonal made up of two segments 13, 24 normal to the x-axis from 1 and 2, and of the segment 34 of the x-axis. Then I[P34] = 0, and I[P l342] = (yi + y~)/2. We shall denote by r the region of the xy-plane made up of the three rectangles [Xi - b ~ X ~ Xi + b, 0 ~ y ~ Yi + b], i = 1,2, and [Xl ~ X ~ X2, 0 ~ y ~ b]. The remark above shows that if b > 0 is suffi­ciently small, then for any curve C 12 joining 1 and 2 in r and distinct from P 1342 we have I[C 12] > I[P1342]. In other words, P 1342 gives a strong relative minimum for I[ C], or more precisely, an absolute minimum in the class of all parametric (regular) curves joining 1 to 2 in r.

We can prove that I[C 12] > I[P1342] for all parametric regular curves C12 joining 1 and 2 in the upper half plane y ~ 0 provided 2 = (X2' Y2) is below or on the curve G. Precisely, we prove the relation I[C 12] > I[P1342] for every curve C12 having at least one point, say 5, on G. We may well assume that 5 is the first point of C 12 on G. If 5' is any point on C12 between 1 and 5, then 5' is above G, and we have I[ClS .] ~ I[E151

y

G

Page 33: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

148 Chapter 3 Examples and Exercises on Classical Problems

As 5' approaches 5 on C, we obtain I[C 1S] ~ I[E 1S]. Let 6 be the foot of the perpen­dicular from 5 on the x-axis, and 7 any point on G close to 3. By the string property (Section 2.9), we have

I[C1S] ~ I[E1s] = I[E17 + G7S].

If we take the point 7 sufficiently close to 3, then the length of the path E17 + G7S is larger than Y1 + Y2' and hence I[E17 + G7S] > I[P 1365] and finally I[C1s] > I[P 1365]. On the other hand, if 8 is the point on P 24 having the same ordinate as 5, then the length of CS2 is certainly larger than the length of P 2S ' and hence I[Cs2] > I[Ps2], I[P6s] = I[P64s], and finally

I[Cd = I[C 1S] + I[Cd > I[P136] + I[P64s] + I[Ps2] = I[Pl342].

We conclude that for points 2 = (X2' Y2) on or below the curve G, the polygonal line P 1342 gives the absolute minimum for the integral I in the class of all parametric regular curves joining 1 to 2.

When 2 is above the curve G, then we have already denoted by E12 that one catenary which gives a strong relative minimum for I in a V -shaped region V, and we have now the polygonal line P 1342 which gives a strong relative minimum for I in a region r. We should compare I[E12] and I[P1432]. Their difference LI is given by

I[Ed - I[P l432] = f: 2 yds - r1(YI + y~),

where S2 is the length of E12, that is, the value of the arc length parameter s along E12 at 2. If we denote by E the catenary containing the arc E12, the expression above can be thought of as a function of s (that is, 2 moves along E), and we have dLljds2 =

Y2(1 - dY2/ds2) > O. Thus LI is an increasing function as 2 moves along any catenary E from 1 to 2. Obviously LI < 0 at 1, since I[E 12] is zero there. On the other hand, LI > 0 at the point of contact 5 of E with G. Thus, there is on each E 1s a well-determined point 3 where LI = O. It can be proved that the locus of3 is a curve H:y = h(x), x ~ Xl' with h > 0, h' > 0, h > g, h( + (0) = h'( + (0) = + 00, h(X1 + 0) = h'(x1 + 0) = o.

We conclude that for 2 = (X2' Y2) above or on H, the catenary E12 gives the absolute minimum for I. For 2 = (X2' Y2) between G and H we have I[E 12] > I[P 1432], and thus E12 is only a strong local minimum. For 2 below G the absolute minimum (in the class of all parametric regular curves joining 1 to 2) exists and is given by the polygonal line P1432, while I has no absolute minimum in the class of the nonparametric curves joining 1 and 2.

Various parametric representations of the two curves G and H are known; see W. S. Kimball [I], who also gives numerical tables of g(x) and h(x). (For the previous analysis, see G. A. Bliss [I]). For a discussion of the same problem of the surface of revolution of minimum area, only in terms of parametric problems, see L. Tonelli [I].

H

Page 34: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.14 The Principles of Mechanics 149

3.14 The Principles of Mechanics

A. D'Alembert's Principle and Lagrange's Equations of Motion

As usual in classical mechanics let us consider a finite system of N material points P j of mass mj' variously interconnected, say Pj = (Xj' Yj' z), j = 1, ... , N, with respect to an inertial system Oxyz. The principle of virtual work in statics concerns with real displace­ments dPj and virtual displacements bPj (see below), away from boundaries (reversible displacements), and in contact with boundaries (nonreversible displacements). The principle states that the system is in equilibrium if and only if the virtual work ~L of the active forces F j applied to the points Pj is zero for all possible reversible virtual displacements, and nonpositive for all possible nonreversible virtual displacements. In symbols:

N

bL = L F j ' bPj = 0 j=l

[:::::0].

In dynamics the inertial forces I j = -mjaj = (-mjxj, -mjyj, -mjzj) must be added to the other forces, and the principle then states that the actual motions and changes in the configuration of the system are monitored by the sole requirement that

(3.14.1) N

bL = L (Fj - mja) . bPj = 0 j=l

[:::::OJ,

where aj is the vector acceleration of Pj (D'Alembert's principle). Let us assume that the configurations of the system are representable in terms of n

parameters ql, ... , qn (Lagrangean coordinates) and time t,

i= 1, ... ,N,

in such a way that each Pj' actually each coordinate Xj' Yj' Zj' is a function of class C2 of ql,"" qn' t in some domain Dc Rn+l, and that the representation is 1-1 at least locally. Actually, we shall require that the 3N x nmatrix M = [ox)oq., oYj/oq., oz)oqs], j = 1, ... ,N, S = 1, ... ,n, has maximum rank at every point. Then the velocity Vj of Pj is

n

(3.14.2) Vj = (dx)dt, dYj/dt, dz)dt) = dPj/dt = L (oP)oqs)q~ + oP)ot, s= 1

and the real displacements dPj and the virtual displacements bPj are respectively

n

dPj = L (oP)oqs)dqs + (oP)ot)dt, s= 1

(3.14.3)

bPj = L (oP)oqs)dq., j= 1, .. . ,N, S = 1, ... , n. s= 1

In other words, the virtual displacements bPj are computed disregarding the direct dependence of Pj on t, or equivalently, disregarding instant by instant the possible dependence on time of the constraints. For constraints independent of time, then oP)ot = 0, and real and virtual displacements coincide. Moreover, for q = (qlo ... ,qn), in the interior of D and under the assumptions, the displacements are certainly reversible.

Page 35: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

150 Chapter 3 Examples and Exercises on Classical Problems

We shall now derive the Lagrange equations of motion for systems with reversible displacements (away from boundaries). For the sake of simplicity we shall limit our­selves to systems with constraints independent of time. Then (3.14.3) becomes

(3.14.4) j = 1, ... , n,

and, if we note that the Pj depend only on the qs and that we can interchange the differ­entiations with respect to t and the q" we have

OVj/oq; = oPj/oq" (3.14.5)

(d/dt)(oP)oqs) = (%qs)(dP)dt) = ov)oq" j= 1, ... ,N, s = 1, ... , n.

First we need a few remarks on the kinetic energy T = r 1 Ij mjlvjl2. From (3.14.4) we derive

(3.14.6)

with

IVjl2 = vj ' Vj = (~(OP)Oqs)q;)( POP)Oqk)qi)

= I (oP)oqs)(oP)oqk)q~qi, s.k

T = r 1 I mjlvjl2 = 2- 1 I Tskq~qi,

Tsk = I mioP)oqs)(OP)oqk), j

s,k

Thus, T is a quadratic form in the q; with coefficients Tsk depending only on the qs. We shall also need the derivatives

Ps = j)T/j)q~ = I Tskq~, s = 1, ... , n, k

and we denote by p the n-vector p = (Pl' ... ,P.). By Euler theorem on homogeneous functions we have

• • (3.14.7) 2T = I (oT/oq~)q; = I Psq;·

s= 1 s= 1

Here T = 0 if and only if q' = (q~, ... ,q~) = O. The condition is obviously sufficient. To prove the necessity, note that in the opposite case there would be a system of q; not all zero for which T = 0, that is, Vj = 0 for allj = 1, ... , N, hence dXj/dt, dy)dt, dzj/dt, j = 1, ... , n, all zero, and from relations (3.14.4) we would derive that the matrix M above is not of maximum rank.

Let us prove that det[Tsk] i= O. Indeed, in the opposite case there would be a system of q; not all zero with all Ps equal to zero, hence from (3.14.7), T = 0, a contradiction. Thus, T > 0 for all q' i= 0, in other words, T is a positive definite quadratic form.

Now we have only to use relation (3.14.1) (with equality sign) after some transfor­mations. First we write the second relation (3.14.3) in the form bPj = Is (oP)oqs)bq" with bqs instead of dqs since they represent arbitrary displacements, and then we have

where

I Fj ' bPj = I Fj ' I(oPj/oqs)bqs = I Qsbq" j j

Qs = I Fj . (OPj/oqs)· j

Page 36: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.14 The Principles of Mechanics

These Qs are often called the generalized forces. Also

where

I mpj' bPj = I mjaj . I (oPj/oqs)bqs = IRsbq" j j s

Rs = I mpj . (oP)oqs), j

and now (3.14.1) becomes

151

for all bqs' This implies that Qs = R" s = 1, ... , n. By manipulation and relations (3.14.5) we have now

Rs = I mpj . (oPj/oqs) = I mj(dv)dt) . (oP)oqs) j j

= I ml(d/dt)(Vj . oPj/oqs) - (Vj . (d/dt)(oP)oqs)] j

= I mj[(d/dt)(vj . Ov)oq~) - (Vj . Ovj/oqs)] j

= r 1 I mj(d/dt)(oIVjI2/oq~) - 2- 1 I miolvjI2/oqs) j

= (d/dt)(%q~)2-1 I mjlvjl2 - (%qs)r 1 I mAvjl2. j j

Thus, relations Qs = Rs become

(3.14.8) doT oT ----=Q dt oq~ oqs "

s = 1, ... , n.

These are the Lagrange equations of motion. Whenever the generalized forces Qs derive from a potential V = - U, i.e., Qs = oU/oq" s = 1, ... , n, depending on the qs only, then the Lagrange equations can be written in the form

d o(T + U) iJ(T + U) - - =0 dt oq~ oqs '

s = 1, ... , n.

If L = T + U, then these relations are the Euler equations of the Lagrange, or "action" integral

f, '2 J = Ldt. "

Thus, the Lagrange equations of motion can be reworded by saying that, for the case offorces Fi depending on a potential, then in the actual motion of the system the action integral is stationary (cf. Section 2.2, Remark 1). If we note that

L = T + U = I Tksqkq~ + u, k.s

where all Tks and U depend only on the qs and the quadratic form is definite positive, we see that along the actual motion of the system the Legendre necessary condition for a minimum is certainly satisfied (cf. Section 2.2, part (e) of statement (2.2.i)). This state­ment can be improved. Indeed, if q(t) = (q1,'" ,qn), t1 S t S t2, represents the actual motion, not only TSk(q(t)) are the coefficients of a definite positive quadratic form, but the Tsk(q) are the coefficients of a definite positive quadratic form for all q. Finally, if

Page 37: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

152 Chapter 3 Examples and Exercises on Classical Problems

t1 is fixed and we take t2 > t1 sufficiently close to t1 so that the interval (t1> t2] is free of points conjugate to t1> then by (2.12.v) the actual motion represents a strong local minimum for the action integral J. In other words, along the actual motion the action integral is stationary in the large and a strong local minimum in the small.

The action integral may not have any minimum in the large, as the following example shows. Let P be a material point of mass m free to move along the x-axis under the elastic force -kx. Then T = 2- 1mx,2, U = -2- 1kx2 the Euler equation is mx" + kx = 0, or x" + w2x = ° with w 2 = kim, and J = rIm s:; (X,2 - W 2x 2)dt. The accessory equation is r( + w21] = 0, so that f = t1 + n/w is conjugate to t 1. We know from Section 3.9 that J has no minimum for t2 > f.

The result obtained above that the Lagrange equations of motion are the Euler equations of the action integral is important because so much of Chapter 2 holds in theoretical mechanics as well as in the calculus of variations. In particular the Hamilton­Jacobi partial differential equation holds

as/at + H o(t, q, as/aq) = 0, q = (q1,··· ,qn),

and the same theorems (2.11.x), (2.11.xi) relate this partial differential equation to the Lagrange equations of motion.

B. The Theorems of the Quantities of Motion and of Kinetic Moments

Let us denote here by {Pj} and {P~} the finite collections of all the exterior and all interior forces acting at each instant on the point Pj. Then the equations of motion are

as we know mj d 2Pj /dt2 = Ili + Ih, where L simply denotes the vectorial sum of

all such forces acting on Pj. If Pi = (Xi, Yi, Zi), P~ = (X~, Y~, Z~), then the same equa­tions in component form are

(3.14.9) mjd2x)dt2 = LXi + LX~,

By addition we have L mj d2Pj/dt 2 = LL Fi + LL F~, j j

where the last sum is certainly zero because the interior forces are two by two equal and of opposite signs. Thus, Lj mj d2P)dt2 = LLj Fj, or

L mj d2xj/dt2 = LL X'j, L mj d2y)dt2 = LL Y'j, j j

L mj d 2z)dt2 = LL Z'j. j

These equations can be written in the equivalent form (d/dt) Lj mjdPj/dt = LLj F'j, or

(3.14.10)

(d/dt) I mjdx)dt = II X'j, (d/dt) I mjdy)dt = II Yi, j

(d/dt) I mjdz)dt = II Z'j. j

Since mjdPj/dt is called the quantity of motion of Pj' the last relations express the theorems of quantities of motion:

3.14.i (THEOREM). The derivative with respect to time of the sum of the quantities of motion of the N points of the system is equal to the sum of all exterior forces acting on the system.

Page 38: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.14 The Principles of Mechanics 153

3.14.ii (THEoREM). The derivative with respect to time of the projection of the sum of the quantities of motion with respect to any fixed axis is equal to the projection of the resultant of all the exterior forces on that axis.

Moreover, if M = Ij mj is the total mass of the system, and G = (~, 1'/, C) denote the center of gravity of the system, then

and the relations (3.14.10) yield

M d2~/dt2 = II Xj, j

M d21'//dt2 = II Yj, j

The last relation expresses the theorem:

M d2C/dt2 = II Zj. j

3.14.iii (THEOREM). The motion of the center of mass of a system is the same as if it were a material point of mass M to which are applied forces equal and parallel to all exterior forces applied to the single points Pj of the system.

Finally, from (3.14.9), by multiplications and additions, we have, for instance,

(3.14.11) I mix jd 2Yj/dt2 - yjd2x)dt2) = II (xjYj - YjXj) + II (XjY~ - YjX~), j j j

where the last parentheses represent moments of the interior forces with respect to the Oz axis. Note that, for each force, say (XL Yi) from P2 applied to Ph or (xi, Yi) =

(U(X2 - Xl), U(Y2 - Yl»' there is another force (X~, Y~) = (-U(X2 - Xl), -U(Y2 - Yl» applied to P2, and then (Xly il - ylXil ) + (X2Y~ - Y2X~) = o. Thus, the last sum in (3.14.11) is zero, and (3.14.11) reduces to

(d/dt) I mixjdYj/dt - yjdXj/dt) = II (xjYj - YjXj). j j

This last relation expresses the theorem of the kinetic moments:

3.14.iv (THEoREM). The derivative with respect to time of the sum of the moments of the quantities of motion (kinetic moments) of the points of a system with respect to any fixed axis is equal to the sum of the moments of all external forces with respect to the same axis.

C. Instantaneous Axis of Rotation of a Rigid Body

If Oxyz is a system of coordinates attached to a rigid body B, or a mobile reference system, and OXlYIZl is a fixed system of coordinates, then the usual transformation of coordinates holds at any instant t:

(3.14.12)

Xl = Xo + o(x + O(lY + 0(2Z,

YI = Yo + px + PlY + P2Z,

Zl = Zo + ')IX + ')IIY + ')I2Z,

where (xo,Yo,zo) are the coordinates of 0 with respect to OXlYlZI, and (0(,0(100(2;

p, PI, P2; ')I, ')11' ')12) are the cosines of the nine angles of the axes Oxyz with the axes

Page 39: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

154 Chapter 3 Examples and Exercises on Classical Problems

OX1Y1Z1' If B is in motion, then xo, Yo, zo, a, . .. , Y2 are functions of t, while x, y, z are fixed, since these represent the coordinates of a point P of the rigid body B with respect to Oxyz attached to B. By differentiating, the velocity V = (Vx" Yy" v.J = (dxddt, dyddt, dzddt) of P is given by

(3.14.13) Vx , = x~ + xa' + ya~ + za~,} Yy, = y~ + xfJ' + yfJ~ + zfJ~, v., = z~ + xy' + yy~ + zy~

(') = d/dt).

On the other hand, if Vx, Vy, V. denote the projections of V on the Oxyz axes, we have

Vx = a Vx, + fJYy, + Y v." (3.14.14) Yy = al Vx , + fJl Yy, + Yl v."

V. = a2 Vx , + fJ2 Yy, + Y2 v.,. In substituting (3.14.13) in (3.14.14) we must firsttake note ofthe identity a2 + fJ2 + y2 = 1 and analogous ones, which by differentiation yield aa' + fJfJ' + yy' = 0 etc. On the other hand, from the identities ala2 + fJlfJ2 + Y1Y2 = 0, aa2 + fJfJ2 + YY2 = 0, aa1 + fJfJl + YYl = 0, we obtain by differentiation

p = a2a'1 + fJ2fJ'1 + ydl = -(ala~ + fJlfJ~ + Y1Y~)' (3.14.15) q = aa~ + fJfJ~ + yy~ = - (a2a' + fJ2fJ' + yd),

With these definitions and identities the substitution of (3.14.13) in (3.14.14) yields

(3.14.16)

Vx = V~ + qz - ry,

Yy = V~ + rx - pz,

V. = V~ + py - qx,

where VO = (V~, V~, V~) simply indicates the projections of the velocity V of 0 on the moving axes. The other terms in (3.14.16) indicate that the motion of B with respect to any of its own points 0 can be interpreted as due to an instantaneous rotation vector Ow = (p, q, r) of components p, q, r given by (3.14.15). The vector (qz - ry, rx - pz, py - qx) is said to be the moment of the vector Ow with respect to the point P = (x, y, z) of B.

Concerning the kinetic energy T of the system B, first assume that 0 is fixed, that is, the rigid body B moves around the fixed point O. Then V~ = V~ = V~ = 0, and

where

T = r 1 I mix? + y? + z?) j

= 2- 1 I mJ(qZj - ryi + (rXj - pzi + (PYj - qxiJ j

= rl[Ap2 + Bq2 + Cr2 - 2Dqr - 2Erp - 2Fpq],

A = ImiyJ + zJ), j

D = I mjYjZj, j

B = I mj(zJ + xJ), j

E = I mjZjXj, j

C = I mj(xJ + yJ), j

F = I mjXjYJ. j

Page 40: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.14 The Principles of Mechanics 155

Now let us assume instead that also 0 is moving, but that 0 is the center of mass of P, so that Ij mh = 0, Ij mjYj = 0, Ij mjzj = 0. From the relations (3.14.16) we derive that

T = r 1 Imix/ + Y/ + z/)

= 2- 1 (~mj)(V~2 + V~2 + V~2) + 2 -I [Ap2 + Bq2 + Cr2 - 2Dqr - 2Erp - 2Fpq].

Finally, if 0 is again the center of masses of Band Oxyz are the principal axes, then D = E = F = 0, and if v denotes the velocity of 0 and M = Ij mj' then

T = 2- I Mv2 + rl[Ap2 + Bq2 + Cr2].

D. Euler's Equations for the Motion of a Rigid Body around Its Center of Mass

Let us compute, at any instant t, the resultant moment Oa of the quantities of motion of the points Pj of B relative to O. Let a x, a y' a z denote the components of 00' on the mobile axes Oxyz. Thus, for instance, 0' x is the sum of the moments relative to Ox. The quantities of motion of the points Pj have projections mjvjx, mjvjy, mjvjz. The sum O'x of the moments of motion is, therefore, by (3.14.16),

O'x = I miYjVjz - ZjVjy) j

= I mlp(yJ + zJ) - qXjYj - rXjzJ j

= Ap - Fq - Er.

This expression equals oT/op. Hence, by symmetry we also have

O'x = oT/op, O'y = oT/oq, O'z = oT/or.

Now let L, M, N be the sums of the moments of all exterior forces with respect to the axes Oxyz. We know that the resultant moment of all these forces with respect to 0 is a vector OS whose components are still L, M, N. We also know from (3.14.iv) that the absolute velocity il of the point a is equal to OS. The projections of this velocity il thus are equal to those of OS, that is, are equal to L, M, N.

The point 0' has coordinates O'x, O'y, O'z with respect to Oxyz, and thus, as t varies, the components of il with respect to Oxyz are dO' x/dt, dO' y/dt, dO' z/dt. By increasing these components by the quantities qO'z - ray, rax - pO'z, pO'y - qO'x respectively, we obtain the components of il with respect to the fixed axes. We have then the equations

dax/dt + qO'z - ray = L, dO'y/dt + rO'x - pO'z = M,

dO'z/dt + pO'y - qrx = N. (3.14.17)

By taking for Oxyz a system of axes with origin 0 the center of mass, and for axes the principal axes through 0, then T = 2- I (Ap2 + Bq2 + Cr2), O'x = 8T/8p = Ap, O'y = Bq, O'z = Cr, and the equations (3.14.17) become the Euler equations for the motion of a

Page 41: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

156 Chapter 3 Examples and Exercises on Classical Problems

rigid body around its center of gravity 0 and principal axes Oxyz:

(3.14.18)

Adp/dt + (C - B)qr = L,

B dq/dt + (A - C)rp = M,

C dr/dt + (B - A)pq = N.

We shall see another derivation of these equations in subsection F. We shall use these relations in a number of examples and exercises in Sections 6.1-3.

E. The Euler Angles for the Motion of a Rigid Body about One of Its Points

The Euler angles 8, qJ, t{! can be thought of as independent parameters, or Lagrange coordinates, for the description of the motion of a rigid body B about one of its points.

Let us consider a fixed orthogonal system OX1YIZ1' and another system Oxyz,

similarly oriented, attached to the rigid body and thus in motion with respect to Ox lY 1 z 1.

Let the orientation be chosen in such a way that a rotation of 90° in the positive sense around the z-axis takes the x-axis into the y-axis, and the same occurs for Oxyz.

Let 01 denote the intersection of the xy-plane with the x 1Yt-plane, and on this straight line define the positive direction to be from 0 to I. Let t{! be the angle, measured in the positive direction with respect to the ZI-axis, which takes the axis OXI into the direction 01. Now the axis Oz is orthogonal to the plane lOx; let qJ denote the angle, measured in the positive direction with respect to the z-axis, which takes 01 into Ox.

Finally, the axis 01 is orthogonal to the plane OZZI; let 8 denote the angle, measured in the positive direction with respect to the OI-axis, which takes OZI to Oz. According to all these conventions, the angle lOy is qJ + n/2. The angles 8, qJ, t{! are independent. To each set of values of these three angles, there corresponds one and only one position of Oxyz with respect to OX1YIZI.

We may think of8, qJ, t{! as functions of time, and their derivatives are then 8', qJ', t{!'. Actually, each of these derivatives represents a rotation, and thus a rotation vector, for which we use the same symbol. The instantaneous rotation ill of the system Oxyz with respect to OX1YIZI is the resultant of these rotations t{!', 8', qJ' around OZI' 01, Oz. These three rotation components are represented in the diagram by vectors equal to t{!', 8', qJ'

Page 42: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

3.14 The Principles of Mechanics 157

on the axes OZl' 01, Oz. The resultant rotational W is the geometric sum of these vectors: its projection on any axis is equal to the sum of the projections of the components l/I/, ()/, q>' on the same axis.

First we determine the projection of OJ on the three orthogonal axes 01, OJ, Oz, where OJ is on the xOy plane and forms the angle +n/2 with 01. Let Wi' Wj' w. denote the three projections. Note that the vector tfr/ is in the plane zOJ, and thus tfr/ can be decomposed into its components l/Ii and l/I~ on OJ and z, namely, l/Ii = l/I/ sin (), l/I~ = l/I/ cos (). Thus, the three components Wi' Wj' W. on the axes 01, OJ, Oz are

Wi=(}', Wj = l/I/ sin (), W. = l/I/ cos () + q/.

In order to obtain p and q, it suffices to take the sum of the projections of the components Wi' Wj on Ox and Oy. Since the orthogonal axes 10J are on the plane xOy, and Ox makes an angle qJ with 01, we have

p = Wx = Wi cos qJ + Wj sin qJ,

q = Wy = Wi cOS(qJ + n/2) + Wj cos qJ.

Finally, we have the relations

(3.14.19)

p = l/I/ sin () sin qJ + ()/ cos qJ,

q = l/I/ sin () cos qJ - ()/ sin qJ,

r = l/I/ cos () + qJ/.

F. Derivation of the Euler Equations as Lagrange Equations of Motion

Having determined a system of Lagrange coordinates l/I, (), qJ for the motion of a rigid body B with respect to its center of mass, we have only to apply the Lagrange equations (3.14.8). Choosing for the axes Oxyz the principal axes through the center of mass 0, we have T = 2- 1(Ap2 + Bq2 + Cr2), with p, q, r given by the relations (3.14.19). For l/I, (), qJ as Lagrange coordinates, let 'P, 0, cP denote the corresponding functions Qj of Subsection A.

The Lagrange equation relative to the variable qJ is now

(3.14.20) (d/dt)(8T/8qJ') - 8T/8qJ = CP,

where now, from the expression for T above and the relations (3.14.19), we derive

8T/8qJ' = (8T/8r)(8r/8qJ') = Cr,

8p/8qJ = q, 8q/8qJ = - p,

8T/8qJ = (8T/8p)(8p/8qJ) + (8T/8q)(8q/8qJ)

= Ap(8p/8qJ) + Bq(8q/8qJ) = (A - B)pq.

Thus, (3.14.20) becomes C(dr/dt) + (B - A)pq = N.

It remains to show that cP = N, the moment of all given forces with respect to Oz. To this effect, note that cP fJqJ is the sum of all virtual works of forces in an elementary

Page 43: Optimization—Theory and Applications || Examples and Exercises on Classical Problems

158 Chapter 3 Examples and Exercises on Classical Problems

displacement that leaves 1/1, e constant, i.e. a rotation [)(p around Oz. Thus

tPfHp = I (XjbXj + ljbYj + Zjbzj) = l:{xjlj - yjXj)fHp, j

or

and (3.14.20) becomes C(dr/dt) + (B - A)pq = N.

The quantities 1/1, e, q> do not appear. Since p, q, r play exactly the same roles, by sym­metry the analogous equations must also hold:

A(dp/dt) + (C - B)qr = L,

B(dq/dt) + (A - C)rp = M.

We have again obtained the Euler equations (3.14.18) for the motion of a rigid body.

Bibliographical Notes

For the classical examples in Sections 3.1-4, 3.6-7 the reader may compare R. Weinstock [I], c. Fox [I], I. M. Gelfand and S. V. Fomin [I], L. E. Elsgolc [I]. For the example in Section 3.5 and a great many applications of the calculus of variations to economics we refer to G. Hadley and M. C. Kemp [I], M. D. Intriligator [I], and P. Newman [I]. For the discussion ofthe example in Section 3.9 the reader may compare L. Tonelli [I]. For the examples in Sections 3.12-13 more details are given in G. A. Bliss [I]. In Section 3.14 we have presented a brief introduction to rational mechanics, both as an illustration of a subject whose development has been in the past inextricably related to the calculus of variations, and for the derivation of the equations of motion of rigid bodies, equations which will be used in problems of optimal control in Chapter 6 and which are, after all, the necessary (Euler) conditions for certain problems of the calculus of variations.

The few topics of theoretical mechanics in Section 3.14 are modeled on P. Appell [I, Vol. 2]. For further references on theoretical mechanics we also mention D. Graffi [I], S. Goldstein [I], and C. Lanczos [I].