Envelope (mathematics)

From HandWiki
Short description: Family of curves in geometry
Construction of the envelope of a family of curves.

In geometry, an envelope of a planar family of curves is a curve that is tangent to each member of the family at some point, and these points of tangency together form the whole envelope. Classically, a point on the envelope can be thought of as the intersection of two "infinitesimally adjacent" curves, meaning the limit of intersections of nearby curves. This idea can be generalized to an envelope of surfaces in space, and so on to higher dimensions.

To have an envelope, it is necessary that the individual members of the family of curves are differentiable curves as the concept of tangency does not apply otherwise, and there has to be a smooth transition proceeding through the members. But these conditions are not sufficient – a given family may fail to have an envelope. A simple example of this is given by a family of concentric circles of expanding radius.

Envelope of a family of curves

Let each curve Ct in the family be given as the solution of an equation ft(xy)=0 (see implicit curve), where t is a parameter. Write F(txy)=ft(xy) and assume F is differentiable.

The envelope of the family Ct is then defined as the set [math]\displaystyle{ \mathcal{D} }[/math] of points (x,y) for which, simultaneously,

[math]\displaystyle{ F(t, x, y) = 0~~\mathsf{and}~~{\partial F \over \partial t}(t, x, y) = 0 }[/math]

for some value of t, where [math]\displaystyle{ \partial F/\partial t }[/math] is the partial derivative of F with respect to t.[1]

If t and u, tu are two values of the parameter then the intersection of the curves Ct and Cu is given by

[math]\displaystyle{ F(t, x, y) = F(u, x, y) = 0\, }[/math]

or, equivalently,

[math]\displaystyle{ F(t, x, y) = 0~~\mathsf{and}~~\frac{F(u, x, y)-F(t, x, y)}{u-t} = 0. }[/math]

Letting ut gives the definition above.

An important special case is when F(txy) is a polynomial in t. This includes, by clearing denominators, the case where F(txy) is a rational function in t. In this case, the definition amounts to t being a double root of F(txy), so the equation of the envelope can be found by setting the discriminant of F to 0 (because the definition demands F=0 at some t and first derivative =0 i.e. its value 0 and it is min/max at that t).

For example, let Ct be the line whose x and y intercepts are t and 11−t, this is shown in the animation above. The equation of Ct is

[math]\displaystyle{ \frac{x}{t}+\frac{y}{11-t}=1 }[/math]

or, clearing fractions,

[math]\displaystyle{ x(11-t)+yt-t(11-t)=t^2+(-x+y-11)t+11x=0.\, }[/math]

The equation of the envelope is then

[math]\displaystyle{ (-x+y-11)^2-44x=(x-y)^2-22(x+y)+121=0.\, }[/math]

Often when F is not a rational function of the parameter it may be reduced to this case by an appropriate substitution. For example, if the family is given by Cθ with an equation of the form u(xy)cos θ+v(xy)sin θ=w(xy), then putting t=eiθ, cos θ=(t+1/t)/2, sin θ=(t-1/t)/2i changes the equation of the curve to

[math]\displaystyle{ u{1 \over 2}(t+{1\over t})+v{1 \over 2i}(t-{1\over t})=w }[/math]

or

[math]\displaystyle{ (u-iv)t^2-2wt+(u+iv)=0.\, }[/math]

The equation of the envelope is then given by setting the discriminant to 0:

[math]\displaystyle{ (u-iv)(u+iv)-w^2=0\, }[/math]

or

[math]\displaystyle{ u^2+v^2=w^2.\, }[/math]

Alternative definitions

  1. The envelope E1 is the limit of intersections of nearby curves Ct.
  2. The envelope E2 is a curve tangent to all of the Ct.
  3. The envelope E3 is the boundary of the region filled by the curves Ct.

Then [math]\displaystyle{ E_1 \subseteq \mathcal{D} }[/math], [math]\displaystyle{ E_2 \subseteq \mathcal{D} }[/math] and [math]\displaystyle{ E_3 \subseteq \mathcal{D} }[/math], where [math]\displaystyle{ \mathcal{D} }[/math] is the set of points defined at the beginning of this subsection's parent section.

Examples

Example 1

These definitions E1, E2, and E3 of the envelope may be different sets. Consider for instance the curve y = x3 parametrised by γ : RR2 where γ(t) = (t,t3). The one-parameter family of curves will be given by the tangent lines to γ.

First we calculate the discriminant [math]\displaystyle{ \mathcal D }[/math]. The generating function is

[math]\displaystyle{ F(t,(x,y)) = 3t^2x - y - 2t^3. }[/math]

Calculating the partial derivative Ft = 6t(xt). It follows that either x = t or t = 0. First assume that x = t and t ≠ 0. Substituting into F: [math]\displaystyle{ F(t,(t,y)) = t^3 - y \, }[/math] and so, assuming that t ≠ 0, it follows that F = Ft = 0 if and only if (x,y) = (t,t3). Next, assuming that t = 0 and substituting into F gives F(0,(x,y)) = −y. So, assuming t = 0, it follows that F = Ft = 0 if and only if y = 0. Thus the discriminant is the original curve and its tangent line at γ(0):

[math]\displaystyle{ \mathcal{D} = \{(x,y) \in \R^2 : y = x^3\} \cup \{(x,y) \in \R^2 : y = 0 \} \ . }[/math]

Next we calculate E1. One curve is given by F(t,(x,y)) = 0 and a nearby curve is given by F(t + ε,(x,y)) where ε is some very small number. The intersection point comes from looking at the limit of F(t,(x,y)) = F(t + ε,(x,y)) as ε tends to zero. Notice that F(t,(x,y)) = F(t + ε,(x,y)) if and only if

[math]\displaystyle{ L := F(t,(x,y)) - F(t+\varepsilon,(x,y)) = 2\varepsilon^3+6\varepsilon t^2+6\varepsilon^2t-(3\varepsilon^2+6\varepsilon t)x = 0. }[/math]

If t ≠ 0 then L has only a single factor of ε. Assuming that t ≠ 0 then the intersection is given by

[math]\displaystyle{ \lim_{\varepsilon \to 0} \frac{1}{\varepsilon} L = 6t(t-x) \ . }[/math]

Since t ≠ 0 it follows that x = t. The y value is calculated by knowing that this point must lie on a tangent line to the original curve γ: that F(t,(x,y)) = 0. Substituting and solving gives y = t3. When t = 0, L is divisible by ε2. Assuming that t = 0 then the intersection is given by

[math]\displaystyle{ \lim_{\varepsilon \to 0} \frac{1}{\varepsilon^2} L = 3x \ . }[/math]

It follows that x = 0, and knowing that F(t,(x,y)) = 0 gives y = 0. It follows that

[math]\displaystyle{ E_1 = \{(x,y) \in \R^2 : y = x^3 \} \ . }[/math]

Next we calculate E2. The curve itself is the curve that is tangent to all of its own tangent lines. It follows that

[math]\displaystyle{ E_2 = \{(x,y) \in \R^2 : y = x^3 \} \ . }[/math]

Finally we calculate E3. Every point in the plane has at least one tangent line to γ passing through it, and so region filled by the tangent lines is the whole plane. The boundary E3 is therefore the empty set. Indeed, consider a point in the plane, say (x0,y0). This point lies on a tangent line if and only if there exists a t such that

[math]\displaystyle{ F(t,(x_0,y_0)) = 3t^2x_0 - y_0 - 2t^3 = 0 \ . }[/math]

This is a cubic in t and as such has at least one real solution. It follows that at least one tangent line to γ must pass through any given point in the plane. If y > x3 and y > 0 then each point (x,y) has exactly one tangent line to γ passing through it. The same is true if y < x3 y < 0. If y < x3 and y > 0 then each point (x,y) has exactly three distinct tangent lines to γ passing through it. The same is true if y > x3 and y < 0. If y = x3 and y ≠ 0 then each point (x,y) has exactly two tangent lines to γ passing through it (this corresponds to the cubic having one ordinary root and one repeated root). The same is true if yx3 and y = 0. If y = x3 and x = 0, i.e., x = y = 0, then this point has a single tangent line to γ passing through it (this corresponds to the cubic having one real root of multiplicity 3). It follows that

[math]\displaystyle{ E_3 = \varnothing. }[/math]

Example 2

This plot gives the envelope of the family of lines connecting points (t,0), (0,k - t), in which k takes the value 1.

In string art it is common to cross-connect two lines of equally spaced pins. What curve is formed?

For simplicity, set the pins on the x- and y-axes; a non-orthogonal layout is a rotation and scaling away. A general straight-line thread connects the two points (0, kt) and (t, 0), where k is an arbitrary scaling constant, and the family of lines is generated by varying the parameter t. From simple geometry, the equation of this straight line is y = −(k − t)x/t + k − t. Rearranging and casting in the form F(x,y,t) = 0 gives:

[math]\displaystyle{ F(x,y,t)=-\frac{kx}{t} - t + x + k -y = 0\, }[/math]

 

 

 

 

(1)

Now differentiate F(x,y,t) with respect to t and set the result equal to zero, to get

[math]\displaystyle{ \frac{\partial F(x,y,t)}{\partial t}= \frac{kx}{t^2} - 1 = 0\, }[/math]

 

 

 

 

(2)

These two equations jointly define the equation of the envelope. From (2) we have:

[math]\displaystyle{ t = \sqrt{kx} \, }[/math]

Substituting this value of t into (1) and simplifying gives an equation for the envelope:

[math]\displaystyle{ y=(\sqrt{x}-\sqrt{k})^2\, }[/math]

 

 

 

 

(3)

Or, rearranging into a more elegant form that shows the symmetry between x and y:

[math]\displaystyle{ \sqrt{x}+\sqrt{y}=\sqrt{k} }[/math]

 

 

 

 

(4)

We can take a rotation of the axes where the b axis is the line y=x oriented northeast and the a axis is the line y=−x oriented southeast. These new axes are related to the original x-y axes by x=(b+a)/2 and y=(ba)/2 . We obtain, after substitution into (4) and expansion and simplification,

[math]\displaystyle{ b = \frac{1}{k\sqrt{2}} a^2 + \frac{k}{2\sqrt{2}}, }[/math]

 

 

 

 

(5)

which is apparently the equation for a parabola with axis along a=0, or y=x.

Example 3

Let IR be an open interval and let γ : IR2 be a smooth plane curve parametrised by arc length. Consider the one-parameter family of normal lines to γ(I). A line is normal to γ at γ(t) if it passes through γ(t) and is perpendicular to the tangent vector to γ at γ(t). Let T denote the unit tangent vector to γ and let N denote the unit normal vector. Using a dot to denote the dot product, the generating family for the one-parameter family of normal lines is given by F : I × R2R where

[math]\displaystyle{ F(t,{\mathbf x}) = ({\mathbf x} - \gamma(t)) \cdot {\mathbf T}(t) \ . }[/math]

Clearly (x − γ)·T = 0 if and only if x − γ is perpendicular to T, or equivalently, if and only if x − γ is parallel to N, or equivalently, if and only if x = γ + λN for some λ ∈ R. It follows that

[math]\displaystyle{ L_{t_0} := \{ {\mathbf x} \in \R^2 : F(t_0,{\mathbf x}) = 0 \} }[/math]

is exactly the normal line to γ at γ(t0). To find the discriminant of F we need to compute its partial derivative with respect to t:

[math]\displaystyle{ \frac{\partial F}{\partial t}(t,{\mathbf x}) = \kappa (t) ({\mathbf x}-\gamma(t))\cdot {\mathbf N}(t) - 1 \ , }[/math]

where κ is the plane curve curvature of γ. It has been seen that F = 0 if and only if x - γ = λN for some λ ∈ R. Assuming that F = 0 gives

[math]\displaystyle{ \frac{\partial F}{\partial t} = \lambda \kappa(t) - 1 \ . }[/math]

Assuming that κ ≠ 0 it follows that λ = 1/κ and so

[math]\displaystyle{ \mathcal{D} = \gamma(t) + \frac{1}{\kappa(t)}{\mathbf N}(t) \ . }[/math]

This is exactly the evolute of the curve γ.

Example 4

An astroid as the envelope of the family of lines connecting points (s,0), (0,t) with s2 + t2 = 1

The following example shows that in some cases the envelope of a family of curves may be seen as the topologic boundary of a union of sets, whose boundaries are the curves of the envelope. For [math]\displaystyle{ s\gt 0 }[/math] and [math]\displaystyle{ t\gt 0 }[/math] consider the (open) right triangle in a Cartesian plane with vertices [math]\displaystyle{ (0,0) }[/math], [math]\displaystyle{ (s,0) }[/math] and [math]\displaystyle{ (0,t) }[/math]

[math]\displaystyle{ T_{s,t}:=\left\{(x,y)\in\R_+^2:\ \frac{x}{s}+\frac{y}{t}\lt 1\right\}. }[/math]

Fix an exponent [math]\displaystyle{ \alpha\gt 0 }[/math], and consider the union of all the triangles [math]\displaystyle{ T_{s,t} }[/math] subjected to the constraint [math]\displaystyle{ \textstyle s^\alpha+t^\alpha=1 }[/math], that is the open set

[math]\displaystyle{ \Delta_\alpha:=\bigcup_ {s^\alpha+t^\alpha=1} T_{s,t}. }[/math]

To write a Cartesian representation for [math]\displaystyle{ \textstyle\Delta_\alpha }[/math], start with any [math]\displaystyle{ \textstyle s\gt 0 }[/math], [math]\displaystyle{ \textstyle t\gt 0 }[/math] satisfying [math]\displaystyle{ \textstyle s^\alpha+t^\alpha=1 }[/math] and any [math]\displaystyle{ \textstyle(x,y)\in\R_+^2 }[/math]. The Hölder inequality in [math]\displaystyle{ \textstyle\R^2 }[/math] with respect to the conjugated exponents [math]\displaystyle{ p:=1+\frac{1}{\alpha} }[/math] and [math]\displaystyle{ \textstyle q:={1+\alpha} }[/math] gives:

[math]\displaystyle{ x^\frac{\alpha}{\alpha+1}+y^\frac{\alpha}{\alpha+1}\leq \left(\frac{x}{s}+\frac{y}{t}\right)^\frac{\alpha}{\alpha+1}\Big(s^\alpha+t^\alpha\Big)^\frac{1}{\alpha+1}=\left(\frac{x}{s}+\frac{y}{t}\right)^\frac{\alpha}{\alpha+1} }[/math],

with equality if and only if [math]\displaystyle{ \textstyle s:\,t=x^\frac{1}{1+\alpha}:\,y^\frac{1}{1+\alpha} }[/math]. In terms of a union of sets the latter inequality reads: the point [math]\displaystyle{ (x,y)\in\R_+^2 }[/math] belongs to the set [math]\displaystyle{ \textstyle\Delta_\alpha }[/math], that is, it belongs to some [math]\displaystyle{ \textstyle T_{s,t} }[/math] with [math]\displaystyle{ \textstyle s^\alpha+t^\alpha=1 }[/math], if and only if it satisfies

[math]\displaystyle{ x^\frac{\alpha}{\alpha+1}+y^\frac{\alpha}{\alpha+1}\lt 1. }[/math]

Moreover, the boundary in [math]\displaystyle{ \R_+^2 }[/math] of the set [math]\displaystyle{ \textstyle \Delta_\alpha }[/math] is the envelope of the corresponding family of line segments

[math]\displaystyle{ \left\{(x,y)\in\R_+^2:\ \frac{x}{s}+\frac{y}{t}=1\right\}\ ,\qquad s^\alpha+t^\alpha=1 }[/math]

(that is, the hypotenuses of the triangles), and has Cartesian equation

[math]\displaystyle{ x^\frac{\alpha}{\alpha+1}+y^\frac{\alpha}{\alpha+1}=1. }[/math]

Notice that, in particular, the value [math]\displaystyle{ \alpha=1 }[/math] gives the arc of parabola of the Example 2, and the value [math]\displaystyle{ \alpha=2 }[/math] (meaning that all hypotenuses are unit length segments) gives the astroid.

Example 5

The orbits' envelope of the projectiles (with constant initial speed) is a concave parabola. The initial speed is 10 m/s. We take g = 10 m/s2.

We consider the following example of envelope in motion. Suppose at initial height 0, one casts a projectile into the air with constant initial velocity v but different elevation angles θ. Let x be the horizontal axis in the motion surface, and let y denote the vertical axis. Then the motion gives the following differential dynamical system:

[math]\displaystyle{ \frac{d^2 y}{dt^2} = -g,\; \frac{d^2 x}{dt^2} = 0, }[/math]

which satisfies four initial conditions:

[math]\displaystyle{ \frac{dx}{dt}\bigg|_{t=0} = v \cos \theta,\; \frac{dy}{dt}\bigg|_{t=0} = v \sin \theta,\; x\bigg|_{t=0} = y\bigg|_{t=0} = 0. }[/math]

Here t denotes motion time, θ is elevation angle, g denotes gravitational acceleration, and v is the constant initial speed (not velocity). The solution of the above system can take an implicit form:

[math]\displaystyle{ F(x,y,\theta) = x\tan \theta - \frac{gx^2}{2v^2 \cos^2 \theta} - y = 0. }[/math]

To find its envelope equation, one may compute the desired derivative:

[math]\displaystyle{ \frac{\partial F}{\partial \theta} = \frac{x}{\cos^2 \theta} - \frac{gx^2 \tan \theta}{v^2 \cos^2 \theta} = 0. }[/math]

By eliminating θ, one may reach the following envelope equation:

[math]\displaystyle{ y = \frac{v^2}{2g} - \frac{g}{2v^2}x^2. }[/math]

Clearly the resulted envelope is also a concave parabola.

Envelope of a family of surfaces

A one-parameter family of surfaces in three-dimensional Euclidean space is given by a set of equations

[math]\displaystyle{ F(x,y,z,a)=0 }[/math]

depending on a real parameter a.[2] For example, the tangent planes to a surface along a curve in the surface form such a family.

Two surfaces corresponding to different values a and a' intersect in a common curve defined by

[math]\displaystyle{ F(x,y,z,a)=0,\,\,{F(x,y,z,a^\prime)-F(x,y,z,a)\over a^\prime -a}=0. }[/math]

In the limit as a' approaches a, this curve tends to a curve contained in the surface at a

[math]\displaystyle{ F(x,y,z,a)=0,\,\,{\partial F\over \partial a}(x,y,z,a)=0. }[/math]

This curve is called the characteristic of the family at a. As a varies the locus of these characteristic curves defines a surface called the envelope of the family of surfaces.

The envelope of a family of surfaces is tangent to each surface in the family along the characteristic curve in that surface.

Generalisations

The idea of an envelope of a family of smooth submanifolds follows naturally. In general, if we have a family of submanifolds with codimension c then we need to have at least a c-parameter family of such submanifolds. For example: a one-parameter family of curves in three-space (c = 2) does not, generically, have an envelope.

Applications

Ordinary differential equations

Envelopes are connected to the study of ordinary differential equations (ODEs), and in particular singular solutions of ODEs.[3] Consider, for example, the one-parameter family of tangent lines to the parabola y = x2. These are given by the generating family F(t,(x,y)) = t2 – 2tx + y. The zero level set F(t0,(x,y)) = 0 gives the equation of the tangent line to the parabola at the point (t0,t02). The equation t2 – 2tx + y = 0 can always be solved for y as a function of x and so, consider

[math]\displaystyle{ t^2 - 2tx + y(x) = 0. \ }[/math]

Substituting

[math]\displaystyle{ t = \left(\frac{dy}{dx}\right)/2 }[/math]

gives the ODE

[math]\displaystyle{ \left(\frac{dy}{dx}\right)^2 \!\! - 4x\frac{dy}{dx} + 4y = 0. }[/math]

Not surprisingly y = 2tx − t2 are all solutions to this ODE. However, the envelope of this one-parameter family of lines, which is the parabola y = x2, is also a solution to this ODE. Another famous example is Clairaut's equation.

Partial differential equations

Envelopes can be used to construct more complicated solutions of first order partial differential equations (PDEs) from simpler ones.[4] Let F(x,u,Du) = 0 be a first order PDE, where x is a variable with values in an open set Ω ⊂ Rn, u is an unknown real-valued function, Du is the gradient of u, and F is a continuously differentiable function that is regular in Du. Suppose that u(x;a) is an m-parameter family of solutions: that is, for each fixed a ∈ A ⊂ Rm, u(x;a) is a solution of the differential equation. A new solution of the differential equation can be constructed by first solving (if possible)

[math]\displaystyle{ D_a u(x;a) = 0\, }[/math]

for a = φ(x) as a function of x. The envelope of the family of functions {u(·,a)}aA is defined by

[math]\displaystyle{ v(x) = u(x;\varphi(x)),\quad x\in\Omega, }[/math]

and also solves the differential equation (provided that it exists as a continuously differentiable function).

Geometrically, the graph of v(x) is everywhere tangent to the graph of some member of the family u(x;a). Since the differential equation is first order, it only puts a condition on the tangent plane to the graph, so that any function everywhere tangent to a solution must also be a solution. The same idea underlies the solution of a first order equation as an integral of the Monge cone.[5] The Monge cone is a cone field in the Rn+1 of the (x,u) variables cut out by the envelope of the tangent spaces to the first order PDE at each point. A solution of the PDE is then an envelope of the cone field.

In Riemannian geometry, if a smooth family of geodesics through a point P in a Riemannian manifold has an envelope, then P has a conjugate point where any geodesic of the family intersects the envelope. The same is true more generally in the calculus of variations: if a family of extremals to a functional through a given point P has an envelope, then a point where an extremal intersects the envelope is a conjugate point to P.

Caustics

Reflective caustic generated from a circle and parallel rays

In geometrical optics, a caustic is the envelope of a family of light rays. In this picture there is an arc of a circle. The light rays (shown in blue) are coming from a source at infinity, and so arrive parallel. When they hit the circular arc the light rays are scattered in different directions according to the law of reflection. When a light ray hits the arc at a point the light will be reflected as though it had been reflected by the arc's tangent line at that point. The reflected light rays give a one-parameter family of lines in the plane. The envelope of these lines is the reflective caustic. A reflective caustic will generically consist of smooth points and ordinary cusp points.

From the point of view of the calculus of variations, Fermat's principle (in its modern form) implies that light rays are the extremals for the length functional

[math]\displaystyle{ L[\gamma] = \int_a^b |\gamma'(t)|\,dt }[/math]

among smooth curves γ on [a,b] with fixed endpoints γ(a) and γ(b). The caustic determined by a given point P (in the image the point is at infinity) is the set of conjugate points to P.[6]

Huygens's principle

Light may pass through anisotropic inhomogeneous media at different rates depending on the direction and starting position of a light ray. The boundary of the set of points to which light can travel from a given point q after a time t is known as the wave front after time t, denoted here by Φq(t). It consists of precisely the points that can be reached from q in time t by travelling at the speed of light. Huygens's principle asserts that the wave front set Φq0(s + t) is the envelope of the family of wave fronts Φq(s) for q ∈ Φq0(t). More generally, the point q0 could be replaced by any curve, surface or closed set in space.[7]

See also

References

  1. Bruce, J. W.; Giblin, P. J. (1984), Curves and Singularities, Cambridge University Press, ISBN 0-521-42999-4 
  2. Eisenhart, Luther P. (2008), A Treatise on the Differential Geometry of Curves and Surfaces, Schwarz Press, ISBN 1-4437-3160-9 
  3. Forsyth, Andrew Russell (1959), Theory of differential equations, Six volumes bound as three, New York: Dover Publications , §§100-106.
  4. Evans, Lawrence C. (1998), Partial differential equations, Providence, R.I.: American Mathematical Society, ISBN 978-0-8218-0772-9 .
  5. John, Fritz (1991), Partial differential equations (4th ed.), Springer, ISBN 978-0-387-90609-6, https://archive.org/details/partialdifferent00john_0 .
  6. Born, Max (October 1999), Principle of Optics, Cambridge University Press, ISBN 978-0-521-64222-4 , Appendix I: The calculus of variations.
  7. Arnold, V. I. (1997), Mathematical Methods of Classical Mechanics, 2nd ed., Berlin, New York: Springer-Verlag, ISBN 978-0-387-96890-2, https://archive.org/details/mathematicalmeth0000arno , §46.

External links