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 finite 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 for a smooth curve

Let f:[a,b]\to\R^n be a continuously differentiable function. The length of the curve defined by f can be defined as the limit of the sum of line segment lengths for a regular partition of [a,b] as the number of segments approaches infinity. This means

L(f)=\lim_{N\to\infty}\sum_{i=1}^N \bigg|f(t_i)-f(t_{i-1})\bigg|

where t_i=a+i(b-a)/N=a+i\Delta t for i=0,1,\dots,N. This definition is equivalent to the standard definition of arc length as an integral:

\lim_{N\to\infty}\sum_{i=1}^N \bigg|f(t_i)-f(t_{i-1})\bigg|=\lim_{N\to\infty}\sum_{i=1}^N \left|\frac{f(t_i)-f(t_{i-1})}{\Delta t}\right|\Delta t=\int_a^b \Big|f'(t)\Big|\ dt.

The last equality above is true because the definition of the derivative as a limit implies that there is a positive real function \delta(\epsilon) of positive real \epsilon such that \Delta t<\delta(\epsilon) implies \left|\bigg|\frac{f(t_i)-f(t_{i-1})}{\Delta t}\bigg|-\Big|f'(t_i)\Big|\right|<\epsilon. This means


\sum_{i=1}^N \left|\frac{f(t_i)-f(t_{i-1})}{\Delta t}\right|\Delta t-\sum_{i=1}^N \Big|f'(t_i)\Big|\Delta t

has absolute value less than \epsilon (b-a) for N>(b-a)/\delta(\epsilon). This means that in the limit N\rightarrow\infty, the left term above equals the right term, which is just the Riemann integral of |f'(t)| on [a,b]. This definition of arc length shows that the length of a curve f:[a,b]\rightarrow\mathbb{R}^n continuously differentiable on [a,b] is always finite. In other words, the curve is always rectifiable.

The definition of arc length of a smooth curve as the integral of the norm of the derivative is equivalent to the definition

L(f)=\sup\sum_{i=1}^N \bigg|f(t_i)-f(t_{i-1})\bigg|

where the supremum is taken over all possible partitions a=t_0<t_1<\dots <t_{N-1}<t_N=b of [a,b].[2] This definition is also valid if f is merely continuous, not differentiable.

A curve can be parameterized in infinitely many ways. Let \varphi:[a,b]\to [c,d] be any continuously differentiable bijection. Then g=f\circ\varphi^{-1}:[c,d]\to\R^n is another continuously differentiable parameterization of the curve originally defined by f. The arc length of the curve is the same regardless of the parameterization used to define the curve:

\begin{align}
L(f) &= \int_a^b \Big|f'(t)\Big|\ dt=\int_a^b \Big|g'(\varphi(t))\varphi'(t)\Big|\ dt \\
&= \int_a^b \Big|g'(\varphi(t))\Big|\varphi'(t)\ dt \quad \textrm{since} \ \varphi \ \textrm{must} \ \textrm{be} \ \textrm{non-decreasing} \\
&= \int_c^d \Big|g'(u)\Big|\ du \quad \textrm{using} \ \textrm{integration} \ \textrm{by} \ \textrm{substitution}\\
&= L(g).
\end{align}

Finding arc lengths by integrating

Quarter circle

If a planar curve in \mathbb{R}^2 is defined by the equation  y=f(x), where f is continuously differentiable, 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.

Curves with closed-form solutions for arc length include the catenary, circle, cycloid, logarithmic spiral, parabola, semicubical parabola and 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.

Numerical integration

In most cases, including even simple curves, there are no closed-form solutions for arc length and numerical integration is necessary. Numerical integration of the arc length integral is usually very efficient. For example, consider the problem of finding the length of a quarter of the unit circle by numerically integrating the arc length integral. The upper half of the unit circle can be parameterized as y=\sqrt{1-x^2}. The interval x\in [-\sqrt{2}/2, \sqrt{2}/2] corresponds to a quarter of the circle. Since dy/dx=-x/\sqrt{1-x^2} and 1+(dy/dx)^2 = 1/(1-x^2), the length of a quarter of the unit circle is

\int_{-\sqrt{2}/2}^{\sqrt{2}/2}\frac{1}{\sqrt{1-x^2}} dx.

The 15-point Gauss-Kronrod rule estimate for this integral of 1.570796326808177 differs from the true length of \pi/2 by 1.3e-11 and the 16-point Gaussian quadrature rule estimate of 1.570796326794727 differs from the true length by only 1.7e-13. This means it is possible to evaluate this integral to almost machine precision with only 16 integrand evaluations.

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(\frac{dr}{dt}\right)^2 + r^2\left(\frac{d\theta}{dt}\right)^2 }dt = \int_{\theta(t_1)}^{\theta(t_2)} \sqrt{\left(\frac{dr}{d\theta}\right)^2 + r^2}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(\frac{dr}{dt}\right)^2 + r^2\left(\frac{d\theta}{dt}\right)^2 + r^2\sin^2\theta \left(\frac{d\phi}{dt}\right)^2}dt.

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

\int_{t_1}^{t_2} \sqrt{\left(\frac{dr}{dt}\right)^2 + r^2\left(\frac{d\theta}{dt}\right)^2 + \left(\frac{dz}{dt}\right)^2}dt.

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

Main article: Great-circle distance
Further information: Geodesics on an ellipsoid

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,[3] which implies that 1 kilometre ≈ 0.53995680 nautical miles.[4] 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.[5]

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.[6] 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).[7]

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. Rudin, Walter (1976). Principles of Mathematical Analysis. McGraw-Hill, Inc. p. 137. ISBN 0-07-054235-X.
  3. physics.nist.gov/Pubs/SP811/appenB8.html
  4. CRC Handbook of Chemistry and Physics, page F-254
  5. 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".
  6. 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.
  7. "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 Tuesday, May 03, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.