5.10: Length of a Curve
- Page ID
- 126413
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Lesson 1: Length of a Curve
Learning Objectives
After successful completion of this lesson, you should be able to:
1) estimate the length of a curve given by a continuous function,
2) estimate the length of a curve found by interpolation of discrete data points.
Introduction
The length, \(L\) of a smooth function \(y(x)\) from \((a,f(a))\ to\ (b,f(b))\)is given analytically by Equation \((\PageIndex{1.1})\)
\[L = \int_{a}^{b}{\sqrt{1 + \left( \frac{dy}{dx} \right)^{2}} \ dx}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.1}) \nonumber\]
Over the interval \(\lbrack a,b\rbrack\), the integrand is required to be integrable, \(y(x)\) needs to be differentiable, and \(y'(x)\) has to be continuous. The length can be thought of as the distance one would travel if traversing along the path of the curve. An application would be to calculate the length of the path traversed by a robot that moves from one point to another on a continuous curve.
Plot the function\(\displaystyle \ y = \frac{1}{3}(x^{2} + 2)^{\frac{3}{2}}\) and find the length of the curve from \(x = 0\) to \(x = 3\).
Solution
We use the equation
\[L = \int_{a}^{b}{\sqrt{1 + \left( \frac{dy}{dx} \right)^{2}} \ dx}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E1.1}) \nonumber\]
We have
\[y = \frac{1}{3}(x^{2} + 2)^{\frac{3}{2}} \nonumber\]
So,
\[\begin{split} \frac{dy}{dx} &= \left( \frac{1}{3} \right) \times \left( \frac{3}{2} \right) \times \left( x^{2} + 2 \right)^{\frac{3}{2}\ - 1} \times \left( 2x \right)\\ &= x\sqrt{x^{2} + 2} \end{split} \;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E1.2}) \nonumber\]
Substituting Equation \((\PageIndex{1.E1.2})\) in Equation \((\PageIndex{1.E1.1})\), the length of the curve from \(x = 0\) to \(x = 3\) is
\[\begin{split} L &= \int_{0}^{3}{\sqrt{1 + \left( x\sqrt{x^{2} + 2} \right)^{2}} \ dx}\\ &= \int_{0}^{3}{\sqrt{1 + x^{2}\left( x^{2} + 2 \right)} \ dx}\\ &= \int_{0}^{3}{\sqrt{1 + x^{4} + 2x^{2}} \ dx}\\ &= \int_{0}^{3}{\sqrt{\left( x^{2} + 1 \right)^{2}} \ dx}\\ &= \int_{0}^{3}{(x^{2} + 1) \ dx}\\ &= \left\lbrack \frac{x^{3}}{3} + x \right\rbrack\begin{matrix} 3 \\ \\ 0 \\ \end{matrix}\\ &= 12 \end{split}\]
A function passes consecutively through the three points \({(2,8),~(4,18),~(8,5)}\) and is found by as a \(2^{nd}\) order polynomial interpolant.
a) Find the function.
b) Find the approximate length of the function from \(x = 2\) to \(x = 8\). You can use the Simpson’s \(1/3^{rd}\) rule of integration rule as given by
\[\int_{a}^{b}{f(x)dx} \approx \frac{(b - a)}{6}\left\lbrack f(a) + 4f\left( \frac{a + b}{2} \right) + f(b) \right\rbrack \nonumber\]
c) What is the absolute relative true error for part (b) if the exact length from \((2,8)\) to \((8,5)\) of the second-order polynomial is \(26.042\)?
Solution
a) The second-order interpolant is given by
\[y = a_{0} + a_{1}x + a_{2}x^{2} \nonumber\]
Since it passes through \({(2,8), (4,18), (8,5)}\), we can set up three equations as
\[\begin{split} a_{0} + a_{1}2 + a_{2}2^{2} &= 8\\ a_{0} + a_{1}4 + a_{2}4^{2} &= 18\\ a_{0} + a_{1}8 + a_{2}8^{2} &= 5\end{split}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E2.1} \text{-}1.E2.3) \nonumber\]
Solving Equations \((\PageIndex{1.E2.1})\) through \((\PageIndex{1.E2.3})\) simultaneously gives
\[a_{0} = - 13,\ a_{1} = 13.25,\ a_{2} = - 1.375 \nonumber\]
The resulting interpolant is
\[y = - 13 + 13.25x - 1.375x^{2}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E2.4}) \nonumber\]
b) The length of the curve \(y\left( x \right)\) from \(a\) to \(b\) is given by
\[L = \int_{a}^{b}{\sqrt{1 + \left( \frac{dy}{dx} \right)^{2}} \ dx} \;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E2.5}) \nonumber\]
Since
\[y = - 13 + 13.25x - 1.375x^{2} \nonumber\]
\[\begin{split}\frac{dy}{dx}&= \frac{d}{dx}(- 13 + 13.25x - 1.375x^2)\\ &= 13.25- 1.375(2x) \\ &= 13.25-2.75x\end{split}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E2.6}) \nonumber\]
Substituting Equation \((\PageIndex{1.E2.6})\) in Equation \((\PageIndex{1.E2.5})\), the length of the curve from \(x = 2\) to \(x = 8\) is given by
\[L = \int_{2}^{8}{\sqrt{1 + \left( 13.25 - 2.75x \right)^{2}} \ dx}\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{1.E2.7}) \nonumber\]
Recall from the problem statement of part (b), integral given by Equation \(\PageIndex{1.E2.7})\) is asked to be found approximately by using Simpson’s 1/3rd rule, that is,
\[\int_{a}^{b}{f(x)dx} \approx \frac{(b - a)}{6}\left\lbrack f(a) + 4f\left( \frac{a + b}{2} \right) + f(b) \right\rbrack \nonumber\]
In this case
\[f\left( x \right) = \sqrt{1 + \left( 13.25 - 2.75x \right)^{2}\ },\ a = 2,\ b = 8 \nonumber\]
Hence Equation \((\PageIndex{1.E2.7})\) gives
\[\begin{split}L & = \int_{2}^{8}{f(x)dx} \\ &\approx \frac{(8 - 2)}{6}\left\lbrack f(2) + 4f\left( \frac{2 + 8}{2} \right) + f(8) \right\rbrack\\ &= f(2) + 4f\left( 5 \right) + f(8)\\ &= 7.814 + 4\left( 1.118 \right) + 8.807\\ &= 21.093\end{split}\]
c) \(\text{True value } = 26.042 \nonumber\)
\(\text{Approximate value } = 21.093 \nonumber\)
So the absolute relative true error is
\[\begin{split} \left| \epsilon_{t} \right| &= \left| \frac{\text{True value} - \text{Approximate value}}{\text{True value}} \right| \times 100\\ &= \left| \frac{26.042 - 21.093\ }{26.042} \right| \times 100\\ &= 19.003 \% \end{split}\]
Audiovisual Lecture
Title: Finding an approximate length of a path of an interpolant
Summary: Given discrete data points, you may trace a path through the points by using an interpolant such as polynomial or spline functions. This video goes through an example of finding the length of an interpolant.
Lesson 2: Comparing Lengths of Curves
Learning Objectives
After successful completion of this lesson, you should be able to:
1) find the shortest smooth path through consecutive points, and
2) compare the lengths of different paths.
Introduction
In the previous lesson, you learned how to find the length of a curve numerically. In this lesson, via an example, we apply what we learned to find the shortest but smoothest path to go through consecutive points.
Peter: “Dr. Kaw, I am taking a course in manufacturing. We are solving the following problem. A robot arm with a rapid laser is used to do a quick quality check, such as the radius of a hole, on six holes on a rectangular plate \(15^{\prime\prime} \times 10^{\prime\prime}\) at several points as shown in Table \(\PageIndex{2.1}\) and Figure \(\PageIndex{2.1}\).
| \(x\) | \(y\) |
|---|---|
| \(2.00\) | \(7.2\) |
| \(4.5\) | \(7.1\) |
| \(5.25\) | \(6.0\) |
| \(7.81\) | \(5.0\) |
| \(9.20\) | \(3.5\) |
| \(10.60\) | \(5.0\) |
I am using Excel to fit a fifth-order polynomial through the 6 points. But, when I plot the polynomial, it is taking a long path! (Figure \(\PageIndex{2.2}\))”
Kaw: “Why do you not just join the consecutive points by a straight line, just like the kids do at Pizza Hut™ with those ‘Connect the dots’ activities?”
Peter: “You are making me hungry, and I wish it were that easy. The path of the robot going from one point to another point needs to be smooth so as to avoid sharp jerks in the arm that can otherwise create premature wear and tear of the robot arm.”
Kaw: “As I recall, you took my course in Numerical Methods. What was that …… one year ago?”
Peter: “Yes, your memory is sharp, but my retention from that course – can we not talk about that?!?”
Kaw: “Come into my office. I wrote this program using Maple. See this function, \(f(x)=1/(1+25x^2)\). I am choosing \(7\) points equidistantly (Table \(\PageIndex{2.2}\)) between \(-1\) and \(1\).
Now, look at the sixth order interpolating polynomial and the original function (Figure \(\PageIndex{2.3}\)). See the oscillations in the interpolating polynomial. In 1901, Runge used this example function to show that higher-order interpolation is a bad idea. One of the solutions to your robot path problem is to use quadratic or cubic spline interpolation. That gives a smooth curve with fewer oscillations and a smoother and shorter path.”
| \(x\) | \(\displaystyle y = \frac{1}{1+25x^2}\) |
|---|---|
| \(-1\) | \(0.038462\) |
| \(-0.66667\) | \(0.0826\) |
| \(-0.33333\) | \(0.264706\) |
| \(0\) | \(1\) |
| \(0.333333\) | \(0.264706\) |
| \(0.666667\) | \(0.082569\) |
| \(1\) | \(0.0385\) |
Peter: “Okay. Let’s give that a try.”
Kaw: “Now, let’s try generating a set of cubic splines to go through the data:”
Peter: “Wow! That (Figure \(\PageIndex{2.4}\)) looks much better!”
Kaw: “It may look better, but let’s find out for sure. See if you can combine the two plots (Figure \(\PageIndex{2.5}\)) and compare the lengths of each path.”
Peter: “The length of a path if from to is given by
\[S=\int_{a}^{b}\sqrt{1+\left(\frac{df}{dx}\right)^2} \ dx\;\;\;\;\;\;\;\;\;\;\;\; (\PageIndex{2.1}) \nonumber\]
Right?”
Kaw: “Yes! You solved the problem. See Table \(\PageIndex{2.3}\) for answers.”
Peter: “I guess your class was good for something, after all, Dr. Kaw.”
Kaw: “Are you sure? You could have always fallen back on the connecting-the-dots method. Besides, you don’t want to grow up … you’re a Pizza Hut™ kid, right?”
Peter: “No, that’s a Toys R’ Us™ kid.”
| Type of interpolation | 5th order polynomial | Cubic Spline |
|---|---|---|
| Length of Curve | \(14.919^{\prime\prime}\) | \(11.248^{\prime\prime}\) |
Multiple Choice Test
(1). The arc length of a smooth cartesian curve \(f(x)\) from \(a\) to \(b\) is given by
(A) \(\displaystyle \int_{a}^{b}{f\left( x \right) \ dx}\)
(B) \(\displaystyle {(b - a)}^{2} + {(f\left( b \right) - f\left( a \right))}^{2}\)
(C) \(\displaystyle \int_{a}^{b}\sqrt{\left( 1 + \frac{dy}{dx} \right)} \ dx\)
(D) \(\displaystyle \int_{a}^{b}\sqrt{\left( 1 + \left( \frac{dy}{dx} \right)^{2} \right)} \ dx\)
(2). Which of these definite integrals gives the length of the arc of the function \(f\left( x \right) = 3x + \sin x\) from \(x = 2\) to \(x = 5\)?
(A) \(\displaystyle \int_{2}^{5}{(3x + \sin x) \ dx}\)
(B) \(\displaystyle \int_{2}^{5}{\sqrt{1 + (3 + \cos x)} \ dx}\)
(C) \(\displaystyle \int_{2}^{5}{\sqrt{1 + {(3 + \cos x)}^{2}} \ dx}\)
(D) \(\displaystyle \int_{2}^{5}{\sqrt{1 + {(3x + \sin x)}^{2}} \ dx}\)
(3). The length of the curve \(y = \sqrt{1 - x^{2}}\) from \(x = 0\) to \(x = 1\) is most nearly equal to
(A) \(0.7854\)
(B) \(1.414\)
(C) \(1.518\)
(D) \(1.571\)
(4). A robotic drawing pencil is using linear spline interpolation to trace a path consecutively through three data points \((4,20)\), \((6,10)\), and \((9,25)\). What is the length of the path that the pencil traces if it begins at the first listed data point and ends at the last data point?
(A) \(15.297\)
(B) \(10.198\)
(C) \(7.071\)
(D) \(25.495\)
(5). A path is traversed consecutively through three points \((2,4)\), \((5,11)\), \((8,3)\) using a quadratic polynomial. Estimate the exact length of the path up to at least \(4\) significant digits. Use only MATLAB or some other program to do this problem - do not do it manually as it will take time.
(A) \(16.701\)
(B) \(6.0827\)
(C) \(16.160\)
(D) \(51.387\)
(6). A cable hangs between two poles that are \(50\) feet apart. The shape of the cable is given by
\[f\left( x \right) = a\cosh\left( \frac{x}{a} \right),\ - 25 \leq x \leq 25, \nonumber\]
where \(a\) is dependent on the tension in the cable and the weight per unit length of the cable, and \(x\) is measured between the poles but from the center of the cable. If the length of the cable is 60 feet, the formula that will allow us to find the value of \(a\) is
(A) \(\displaystyle {\ 2\ a^{2}\sinh}\left( \frac{25}{a} \right) = 60\)
(B) \(\displaystyle {\ 2\ a{\ sinh}}\left( \frac{25}{a} \right) = 60\)
(C) \(\displaystyle {\ 2\ a{\ sinh}}\left( \frac{25}{a} \right) = 50\)
(D) \(\displaystyle {\ 2\ a^{2}\sinh}\left( \frac{25}{a} \right) = 50\)
For a complete solution, go to
http://nm.mathforcollege.com/mcquizzes/05inp/quiz_05_10_lengthofcurve_solution.pdf