Arc length

Determining the length of an irregular arc segment is also called rectification of a curve. Historically, many methods were used for specific curves. The advent of infinitesimal calculus led to a general formula that provides closed-form solutions in some cases.

General approach

Approximation by multiple linear segments

A curve in the plane can be approximated by connecting a finite number of points on the curve using line segments to create a polygonal path. Since it is straightforward to calculate the length of each linear segment (using the Pythagorean theorem in Euclidean space, for example), the total length of the approximation can be found by summing the lengths of each linear segment; that approximation is known as the (cumulative) chordal distance.[1]

If the curve is not already a polygonal path, using a progressively larger number of segments of smaller lengths will result in better approximations. The lengths of the successive approximations will not decrease and may keep increasing indefinitely, but for smooth curves they will tend to a limit as the lengths of the segments get arbitrarily small.

For some curves there is a smallest number L that is an upper bound on the length of any polygonal approximation. These curves are called rectifiable and the number L is defined as the arc length.

Definition

Let C be a curve in a metric space X. This means C is the image of a continuous function f:[a,b] \rightarrow X mapping the interval [a,b] into X. From a partition     a=t_0<t_1<\dots < t_{n-1}<t_n=b of the interval [a,b] we obtain a finite collection of points f(t_0), f(t_1),\dots, f(t_n) on the curve C. Denote the distance from f(t_i) to f(t_{i+1}) by d(f(t_i), f(t_{i+1})), which is the length of the line segment connecting the two points. The arc length L of C is then defined to be

L(C) = \sup_{a=t_0 < t_1 < \cdots < t_n = b} \sum_{i = 0}^{n - 1} d(f(t_i), f(t_{i+1}))

where the supremum is taken over all possible partitions of [a,b] and n is any positive integer. If L<\infty then C is rectifiable and non-rectifiable otherwise. This definition of arc length does not require f to be differentiable. The notion of differentiability is not necessarily even defined on a metric space.

A curve may be parametrized in infinitely many ways. Let \varphi be any continuous bijection mapping [a,b] onto [c,d] and define g:[c,d]\rightarrow X by g=f\circ \phi^{-1}. Then g is another parametrization of C. Let u_i = \varphi(t_i) (and hence     t_i = \varphi^{-1}(u_i)) for     i=0,1,\dots,n-1. Then f(t_i) = (g\circ \varphi)(t_i) = g(u_i) and

    \sum_{i = 0}^{n - 1} d(f(t_i), f(t_{i+1})) = \sum_{i = 0}^{n - 1} d(g(u_i), g(u_{i+1})).

This equation means that for any partition     a=t_0<t_1<\dots < t_{n-1}<t_n=b of     [a,b] there is a partition     c=u_0<u_1<\dots <u_{n-1}<u_n = d of     [c,d] that produces exactly the same approximate arc length (and vice versa). So the value of     L(C) will be the same when parametrizing     C by     g as it is when using     f. In other words,     L(C) has the same value regardless of the parametrization of     C.

Finding arc lengths by integrating

Consider a real function f(x) such that f(x) and f'(x)=\frac{dy}{dx} (its derivative with respect to x) are continuous on [a, b]. The length s of the part of the graph of f between x = a and x = b can be found as follows:

Consider an infinitesimal part of the curve ds (or consider this as a limit in which the change in s approaches ds). According to Pythagoras' theorem (ds)^2=(dx)^2+(dy)^2, from which:

(ds)^2=(dx)^2+(dy)^2
\frac{(ds)^2}{(dx)^2}=1+\frac{(dy)^2}{(dx)^2}
\sqrt{\frac{(ds)^2}{(dx)^2}}=\sqrt{1+\frac{(dy)^2}{(dx)^2}}
\frac{ds}{dx}=\sqrt{1+\frac{(dy)^2}{(dx)^2}}
ds=\sqrt{1+\left(\frac{dy}{dx}\right)^2}dx
s = \int_{a}^{b} \sqrt { 1 + [f'(x)]^2 }\, dx.

If a curve is defined parametrically by x = X(t) and y = Y(t), then its arc length between t = a and t = b is

s = \int_{a}^{b} \sqrt { [X'(t)]^2 + [Y'(t)]^2 }\, dt.

This is more clearly a consequence of the distance formula where instead of a \Delta x and \Delta y, we take the limit. An equivalent expression is

s = \lim \sum_a^b \sqrt { (\Delta x)^2 + (\Delta y)^2 } = \int_{a}^{b} \sqrt { (dx)^2 + (dy)^2 } = \int_{a}^{b} \sqrt { \left(\frac{dx}{dt}\right)^2 + \left(\frac{dy}{dt}\right)^2 }\,dt.

If a function is defined as a function of x by  y=f(x) then it is simply a special case of a parametric equation where x = t and  y = f(t) , and the arc length is given by:

s = \int_{a}^{b} \sqrt{ 1 + \left(\frac{dy}{dx}\right)^2 } \, dx.

In most cases, including even simple curves, there are no closed-form solutions of arc length and numerical integration is necessary.

Curves with closed-form solutions for arc length include the catenary, circle, cycloid, logarithmic spiral, parabola, semicubical parabola and (mathematically, a curve) straight line. The lack of a closed form solution for the arc length of an elliptic arc led to the development of the elliptic integrals.

Derivation

For a small piece of curve, ∆s can be approximated with the Pythagorean theorem
A representative linear element of the function y=t 5, x = t 3

In order to approximate the arc length of the curve, it is split into many linear segments. To make the value exact, and not an approximation, infinitely many linear elements are needed. This means that each element is infinitely small. This fact manifests itself later on when an integral is used.

Begin by looking at a representative linear segment (see image) and observe that its length (element of the arc length) will be the differential ds. We will call the horizontal element of this distance dx, and the vertical element dy.

The Pythagorean theorem tells us that

ds = \sqrt{(dx)^2 + (dy)^2}.\,

Since the function is defined in time, segments (ds) are added up across infinitesimally small intervals of time (dt) yielding the integral

\int_a^b \sqrt{\bigg(\frac{dx}{dt}\bigg)^2+\bigg(\frac{dy}{dt}\bigg)^2}\,dt,

If y is a function of x, so that we could take t = x, then we have:

\int_a^b \sqrt{1+\bigg(\frac{dy}{dx}\bigg)^2}\,dx,

which is the arc length from x = a to x = b of the graph of the function ƒ.

For example, the curve in this figure is defined by

\begin{cases} y = t^5, \\ x = t^3. \end{cases}

Subsequently, the arc length integral for values of t from -1 to 1 is

\int_{-1}^1 \sqrt{(3t^2)^2 + (5t^4)^2}\,dt = \int_{-1}^1 \sqrt{9t^4 + 25t^8}\,dt.

Using computational approximations, we can obtain a very accurate arc length of 2.905.

Another way to obtain the integral formula

Approximation by multiple hypotenuses.

Suppose that there exists a rectifiable curve given by a function f(x). To approximate the arc length S along f between two points a and b in that curve, construct a series of right triangles whose concatenated hypotenuses "cover" the arc of curve chosen as shown in the figure. For convenience, the bases of all those triangles can be set equal to  \Delta x , so that for each one an associated  \Delta y exists. The length of any given hypotenuse is given by the Pythagorean Theorem:

\sqrt {(\Delta x)^2 + (\Delta y)^2}

The summation of the lengths of the n hypotenuses approximates S:

S \sim \sum_{i=1}^n \sqrt { (\Delta x_i)^2 + (\Delta y_i)^2 }

Multiplying the radicand by \frac{(\Delta x)^2}{(\Delta x)^2} produces:

\sqrt { (\Delta x)^2 + (\Delta y)^2 }=\sqrt{ {(\Delta x)^2 + ((\Delta y)^2})\,\frac{(\Delta x)^2}{(\Delta x)^2}}=\sqrt { 1 + \frac{(\Delta y)^2}{(\Delta x)^2}}\,\Delta x=\sqrt { 1 + \left(\frac{\Delta y} {\Delta x} \right)^2 }\,\Delta x

Then, our previous result becomes:

S \sim \sum_{i=1}^n \sqrt { 1 + \left(\frac{\Delta y_i} {\Delta x_i} \right)^2 }\,\Delta x_i

As the length  \Delta x of these segments decreases, the approximation improves. The limit of the approximation, as  \Delta x goes to zero, is equal to S:

S = \lim_{\Delta x_i \to 0} \sum_{i=1}^\infty \sqrt { 1 + \left(\frac{\Delta y_i}{\Delta x_i} \right)^2 }\,\Delta x_i = \int_{a}^{b} \sqrt { 1 + \left(\frac{dy}{dx}\right)^2 } \,dx = \int_{a}^{b} \sqrt{1 + \left [ f' \left ( x \right ) \right ] ^2} \, dx.

Curve on a surface

Let \mathbf{x}(u,v) be a surface mapping and let \mathbf{C}(t) = (u(t), v(t)) be a curve on this surface. The integrand of the arc length integral is |(\mathbf{x}\circ\mathbf{C})'(t)|. Evaluating the derivative requires the chain rule for vector fields:

D(\mathbf{x} \circ \mathbf{C}) = (\mathbf{x}_u \ \mathbf{x}_v)\binom{u'}{v'} = \mathbf{x}_u u' + \mathbf{x}_v v'.

The squared norm of this vector is (\mathbf{x}_u u' + \mathbf{x}_v v') \cdot (\mathbf{x}_u u' + \mathbf{x}_v v') = g_{11}(u')^2 + 2g_{12}u'v' + g_{22}(v')^2 (where g_{ij} is the first fundamental form coefficient), so the integrand of the arc length integral can be written as \sqrt{g_{ab} (u^a)' (u^b)'} (where u^1 = u and u^2 = v ).

Other coordinate systems

Let \mathbf{C}(t) = (r(t), \theta(t)) be a curve expressed in polar coordinates. The mapping that transforms from polar coordinates to rectangular coordinates is

\mathbf{x}(r,\theta) = (r\cos\theta, r\sin\theta ).

The integrand of the arc length integral is |(\mathbf{x}\circ\mathbf{C})'(t)|. The chain rule for vector fields shows that D(\mathbf{x}\circ \mathbf{C}) = \mathbf{x}_r r' + \mathbf{x}_{\theta} \theta'. So the squared integrand of the arc length integral is

(\mathbf{x_r}\cdot\mathbf{x}_r)(r')^2 + 2(\mathbf{x}_r\cdot\mathbf{x}_{\theta})r'\theta' + (\mathbf{x}_{\theta}\cdot\mathbf{x}_{\theta})(\theta')^2 = (r')^2 + r^2(\theta')^2.

So for a curve expressed in polar coordinates, the arc length is

\int_{t_1}^{t_2} \sqrt{\left({\operatorname{d}\!r\over\operatorname{d}\!t}\right)^2 + r^2\left({\operatorname{d}\!\theta\over\operatorname{d}\!t}\right)^2 } \operatorname{d}\!t = \int_{\theta(t_1)}^{\theta(t_2)} \sqrt{\left({\operatorname{d}\!r\over\operatorname{d}\!\theta}\right)^2 + r^2 } \operatorname{d}\!\theta.

Now let \mathbf{C}(t) = (r(t), \theta(t), \phi(t)) be a curve expressed in spherical coordinates where \theta is the polar angle measured from the positive z-axis and \phi is the azimuthal angle. The mapping that transforms from spherical coordinates to rectangular coordinates is

\mathbf{x}(r,\theta,\phi) = (r\sin\theta\cos\phi, r\sin\theta\sin\phi, r\cos\theta).

Using the chain rule again shows that D(\mathbf{x}\circ\mathbf{C}) = \mathbf{x}_r r' + \mathbf{x}_{\theta}\theta' + \mathbf{x}_{\phi}\phi'. All dot products \mathbf{x}_i \cdot \mathbf{x}_j where i and j differ are zero, so the squared norm of this vector is

(\mathbf{x}_r\cdot \mathbf{x}_r )(r'^2) + (\mathbf{x}_{\theta} \cdot \mathbf{x}_{\theta})(\theta')^2 + (\mathbf{x}_{\phi}\cdot \mathbf{x}_{\phi})(\phi')^2 = (r')^2 + r^2(\theta')^2 + r^2 \sin^2\theta (\phi')^2.

So for a curve expressed in spherical coordinates, the arc length is

\int_{t_1}^{t_2} \sqrt{\left({\operatorname{d}\!r\over\operatorname{d}\!t}\right)^2 + r^2\left({\operatorname{d}\!\theta\over\operatorname{d}\!t}\right)^2 + r^2\sin^2\theta \left({\operatorname{d}\!\phi\over\operatorname{d}\!t}\right)^2 } \operatorname{d}\!t.

A very similar calculation shows that the arc length of a curve expressed in cylindrical coordinates is

\int_{t_1}^{t_2} \sqrt{\left({\operatorname{d}\!r\over\operatorname{d}\!t}\right)^2 + r^2\left({\operatorname{d}\!\theta\over\operatorname{d}\!t}\right)^2 + \left({\operatorname{d}\!z\over\operatorname{d}\!t}\right)^2 } \operatorname{d}\!t.

Simple cases

Arcs of circles

Arc lengths are denoted by s, since the Latin word for length (or size) is spatium.

In the following lines, r represents the radius of a circle, d is its diameter, C is its circumference, s is the length of an arc of the circle, and \theta is the angle which the arc subtends at the centre of the circle. The distances r, d, C, and s are expressed in the same units.

Arcs of great circles on the Earth

Two units of length, the nautical mile and the metre (or kilometre), were originally defined so the lengths of arcs of great circles on the Earth's surface would be simply numerically related to the angles they subtend at its centre. The simple equation s=\theta applies in the following circumstances:

  • if s is in nautical miles, and \theta is in arcminutes (160 degree), or
  • if s is in kilometres, and \theta is in centigrades (1100 grad).

The lengths of the distance units were chosen to make the circumference of the Earth equal 40,000 kilometres, or 21,600 nautical miles. These are the numbers of the corresponding angle units in one complete turn.

These definitions of the metre and nautical mile have been superseded by more precise ones, but the original definitions are still accurate enough for conceptual purposes, and for some calculations. For example, they imply that one kilometre is exactly 0.54 nautical miles. Using official modern definitions, one nautical mile is exactly 1.852 kilometres,[2] which implies that 1 kilometre ≈ 0.53995680 nautical miles.[3] This modern ratio differs from the one calculated from the original definitions by less than one part in ten thousand.

Length of an arc of a parabola

For a calculation of the length of a parabolic arc, see Parabola#Length of an arc of a parabola.

Historical methods

Antiquity

For much of the history of mathematics, even the greatest thinkers considered it impossible to compute the length of an irregular arc. Although Archimedes had pioneered a way of finding the area beneath a curve with his "method of exhaustion", few believed it was even possible for curves to have definite lengths, as do straight lines. The first ground was broken in this field, as it often has been in calculus, by approximation. People began to inscribe polygons within the curves and compute the length of the sides for a somewhat accurate measurement of the length. By using more segments, and by decreasing the length of each segment, they were able to obtain a more and more accurate approximation. In particular, by inscribing a polygon of many sides in a circle, they were able to find approximate values of π.

17th century

In the 17th century, the method of exhaustion led to the rectification by geometrical methods of several transcendental curves: the logarithmic spiral by Evangelista Torricelli in 1645 (some sources say John Wallis in the 1650s), the cycloid by Christopher Wren in 1658, and the catenary by Gottfried Leibniz in 1691.

In 1659, Wallis credited William Neile's discovery of the first rectification of a nontrivial algebraic curve, the semicubical parabola.[4]

Integral form

Before the full formal development of calculus, the basis for the modern integral form for arc length was independently discovered by Hendrik van Heuraet and Pierre de Fermat.

In 1659 van Heuraet published a construction showing that the problem of determining arc length could be transformed into the problem of determining the area under a curve (i.e., an integral). As an example of his method, he determined the arc length of a semicubical parabola, which required finding the area under a parabola.[5] In 1660, Fermat published a more general theory containing the same result in his De linearum curvarum cum lineis rectis comparatione dissertatio geometrica (Geometric dissertation on curved lines in comparison with straight lines).[6]

Fermat's method of determining arc length

Building on his previous work with tangents, Fermat used the curve

 y = x^{3/2} \,

whose tangent at x = a had a slope of

 \textstyle {3 \over 2} a^{1/2}

so the tangent line would have the equation

 y = \textstyle {3 \over 2} {a^{1/2}}(x - a) + f(a).

Next, he increased a by a small amount to a + ε, making segment AC a relatively good approximation for the length of the curve from A to D. To find the length of the segment AC, he used the Pythagorean theorem:

\begin{align}
AC^2 &{}= AB^2 + BC^2 \\
&{} = \textstyle \varepsilon^2 + {9 \over 4} a \varepsilon^2 \\
&{}= \textstyle \varepsilon^2 \left (1 + {9 \over 4} a \right )
\end{align}

which, when solved, yields

AC = \textstyle \varepsilon \sqrt { 1 + {9 \over 4} a\ }.

In order to approximate the length, Fermat would sum up a sequence of short segments.

Curves with infinite length

The Koch curve.
The graph of xsin(1/x).

As mentioned above, some curves are non-rectifiable. That is, there is no upper bound on the lengths of polygonal approximations; the length can be made arbitrarily large. Informally, such curves are said to have infinite length. There are continuous curves on which every arc (other than a single-point arc) has infinite length. An example of such a curve is the Koch curve. Another example of a curve with infinite length is the graph of the function defined by f(x) = x sin(1/x) for any open set with 0 as one of its delimiters and f(0) = 0. Sometimes the Hausdorff dimension and Hausdorff measure are used to quantify the size of such curves.

Generalization to (pseudo-)Riemannian manifolds

Let M be a (pseudo-)Riemannian manifold, γ : [0, 1] → M a curve in M and g the (pseudo-) metric tensor.

The length of γ is defined to be

\ell(\gamma)=\int_{0}^{1} \sqrt{ \pm g(\gamma'(t),\gamma '(t)) } \, dt,

where γ'(t)Tγ(t)M is the tangent vector of γ at t. The sign in the square root is chosen once for a given curve, to ensure that the square root is a real number. The positive sign is chosen for spacelike curves; in a pseudo-Riemannian manifold, the negative sign may be chosen for timelike curves . Thus the length of a curve in a non-negative real number. Usually no curves are considered which are partly spacelike and partly timelike.

In theory of relativity, arc length of timelike curves (world lines) is the proper time elapsed along the world line, and arc length of a spacelike curve the proper distance along the curve.

See also

References and notes

  1. p.51 in Ahlberg & Nilson (1967) The theory of splines and their applications, Academic Press, 1967
  2. physics.nist.gov/Pubs/SP811/appenB8.html
  3. CRC Handbook of Chemistry and Physics, page F-254
  4. John Wallis, Tractatus Duo. Prior, De Cycloide et de Corporibus inde Genitis. … (Oxford, England: University Press, 1659), pages 91-96; the accompanying figures appear on page 145. On page 91, William Neile is mentioned as "Gulielmus Nelius".
  5. Henricus van Heuraet, "Epistola de transmutatione curvarum linearum in rectas" (Letter on the transformation of curved lines into right ones [i.e., Letter on the rectification of curves]), Renati Des-Cartes Geometria, 2nd ed. (Amsterdam ["Amstelædami"], (Netherlands): Louis & Daniel Elzevir, 1659), pages 517-520.
  6. "M.P.E.A.S." (pseudonym of Fermat) De Linearum Curvarum cum Lineis Rectis Comparatione Dissertatio Geometrica (Toulouse [Tolosæ], France: Arnaud Colomer, 1660).

External links

This article is issued from Wikipedia - version of the Monday, February 15, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.