Vector Functions - Study Guide

13 Vector Functions 13.1 Spa e Curves We have already seen that a convenient way to describe a line in three dimensions is to provide a vector that “points to” every point on the line as a parameter t varies, like h1, 2, 3i + th1, −2, 2i = h1 + t, 2 − 2t, 3 + 2ti. Except that this gives a particularly simple geometric object, there is nothing special about the individual functions of t that make up the coordinates of this vector—any vector with a parameter, like hf (t), g(t), h(t)i, will describe some curve in three dimensions as t varies through all possible values. EXAMPLE 13.1.1 Describe the curves hcos t, sin t, 0i, hcos t, sin t, ti, and hcos t, sin t, 2ti. As t varies, the first two coordinates in all three functions trace out the points on the unit circle, starting with (1, 0) when t = 0 and proceeding counter-clockwise around the circle as t increases. In the first case, the z coordinate is always 0, so this describes precisely the unit circle in the x-y plane. In the second case, the x and y coordinates still describe a circle, but now the z coordinate varies, so that the height of the curve matches the value of t. When t = π, for example, the resulting vector is h−1, 0, πi. A bit of thought should convince you that the result is a helix. In the third vector, the z coordinate varies twice as fast as the parameter t, so we get a stretched out helix. Both are shown in figure 13.1.1. On the left is the first helix, shown for t between 0 and 4π; on the right is the second helix, shown for t between 0 and 2π. Both start and end at the same point, but the first helix takes two full “turns” to get there, because its z coordinate grows more slowly. 329 330 Chapter 13 Vector Functions Figure 13.1.1 Two helixes. (AP) A vector expression of the form hf (t), g(t), h(t)i is called a vector function; it is a function from the real numbers R to the set of all three-dimensional vectors. We can alternately think of it as three separate functions, x = f (t), y = g(t), and z = h(t), that describe points in space. In this case we usually refer to the set of equations as parametric equations for the curve, just as for a line. While the parameter t in a vector function might represent any one of a number of physical quantities, or be simply a “pure number”, it is often convenient and useful to think of t as representing time. The vector function then tells you where in space a particular object is at any time. Vector functions can be difficult to understand, that is, difficult to picture. When available, computer software can be very helpful. When working by hand, one useful approach is to consider the “projections” of the curve onto the three standard coordinate planes. We have already done this in part: in example 13.1.1 we noted that all three curves project to a circle in the x-y plane, since hcos t, sin ti is a two dimensional vector function for the unit circle. EXAMPLE 13.1.2 Graph the projections of hcos t, sin t, 2ti onto the x-z plane and the y-z plane. The two dimensional vector function for the projection onto the x-z plane is hcos t, 2ti, or in parametric form, x = cos t, z = 2t. By eliminating t we get the equation x = cos(z/2), the familiar curve shown on the left in figure 13.1.2. For the projection onto the y-z plane, we start with the vector function hsin t, 2ti, which is the same as y = sin t, z = 2t. Eliminating t gives y = sin(z/2), as shown on the right in figure 13.1.2. Exercises 13.1. 1. Describe the curve r = hsin t, cos t, cos 8ti. 2. Describe the curve r = ht cos t, t sin t, ti. 13.2 Calculus with vector functions z... z... ...... ....... 4π . ... ... .... . . . . . ..... ..... ...... ....... . . . . . . ....... ...... ...... ....... . . . . . . ...... ..... ..... ..... . . ... ... ... .... .... ..... ..... ...... ...... ....... ....... ....... ....... ....... ....... ...... ...... ..... ..... .... ... ... . .. ........ .. 4π 2π −1 Figure 13.1.2 0 331 1 ...... ....... ...... ...... . . . . . ..... ..... .... ... . . .. ... ... .... ..... ..... ..... ...... ...... ....... ....... ....... ....... ....... ...... ...... ..... ..... ..... .... ... ... ... .. . . . . . . ..... ..... ..... ...... . . . . . ... ...... .... ....... ....... 2π x −1 0 1 y The projections of hcos t, sin t, 2ti onto the x-z and y-z planes. 3. Describe the curve r = ht, t2 , cos ti. √ √ 4. Describe the curve r = hcos(20t) 1 − t2 , sin(20t) 1 − t2 , ti 5. Find a vector function for the curve of intersection of x2 + y2 = 9 and y + z = 2. ⇒ 6. A bug is crawling outward along the spoke of a wheel that lies along a radius of the wheel. The bug is crawling at 1 unit per second and the wheel is rotating at 1 radian per second. Suppose the wheel lies in the y-z plane with center at the origin, and at time t = 0 the spoke lies along the positive y axis and the bug is at the origin. Find a vector function r(t) for the position of the bug at time t. ⇒ 7. What is the difference between the parametric curves f (t) = ht, t, t2 i, g(t) = ht2 , t2 , t4 i, and h(t) = hsin(t), sin(t), sin2 (t)ias t runs over all real numbers? 8. Plot each of the curves below in 2 dimensions, projected onto each of the three standard planes (the x-y, x-z, and y-z planes). a. f (t) = ht, t3 , t2 i, t ranges over all real numbers b. f (t) = ht2 , t − 1, t2 + 5i for 0 ≤ t ≤ 3 9. Given points A = (a1 , a2 , a3 ) and B = (b1 , b2 , b3 ), give parametric equations for the line segment connecting A and B. Be sure to give appropriate t values. 10. With a parametric plot and a set of t values, we can associate a ‘direction’. For example, the curve hcos t, sin ti is the unit circle traced counterclockwise. How can we amend a set of given parametric equations and t values to get the same curve, only traced backwards? 13.2 Cal ulus with ve tor fun tions A vector function r(t) = hf (t), g(t), h(t)i is a function of one variable—that is, there is only one “input” value. What makes vector functions more complicated than the functions y = f (x) that we studied in the first part of this book is of course that the “output” values are now three-dimensional vectors instead of simply numbers. It is natural to wonder if there is a corresponding notion of derivative for vector functions. In the simpler case of 332 Chapter 13 Vector Functions a function y = s(t), in which t represents time and s(t) is position on a line, we have seen that the derivative s′ (t) represents velocity; we might hope that in a similar way the derivative of a vector function would tell us something about the velocity of an object moving in three dimensions. One way to approach the question of the derivative for vector functions is to write down an expression that is analogous to the derivative we already understand, and see if we can make sense of it. This gives us r(t + ∆t) − r(t) ∆t→0 ∆t hf (t + ∆t) − f (t), g(t + ∆t) − g(t), h(t + ∆t) − h(t)i = lim ∆t→0 ∆t f (t + ∆t) − f (t) g(t + ∆t) − g(t) h(t + ∆t) − h(t) = lim h , , i ∆t→0 ∆t ∆t ∆t r′ (t) = lim = hf ′ (t), g ′(t), h′ (t)i, if we say that what we mean by the limit of a vector is the vector of the individual coordinate limits. So starting with a familiar expression for what appears to be a derivative, we find that we can make good computational sense out of it—but what does it actually mean? We know how to interpret r(t + ∆t) and r(t)—they are vectors that point to locations in space; if t is time, we can think of these points as positions of a moving object at times that are ∆t apart. We also know what ∆r = r(t + ∆t) − r(t) means—it is a vector that points from the head of r(t) to the head of r(t + ∆t), assuming both have their tails at the origin. So when ∆t is small, ∆r is a tiny vector pointing from one point on the path of the object to a nearby point. As ∆t gets close to 0, this vector points in a direction that is closer and closer to the direction in which the object is moving; geometrically, it approaches a vector tangent to the path of the object at a particular point. ..... ...... ...... ........ . . . . . . . . . ..... ..................................................................... ...................... .............................................. .. ........... ....... ........ ............................... . . . . . . . . . . . . . .. . ....................... ....... ...... ......... ...... ......... ...... ........... . . . . . . . . . . . . . ....... ..... .... ...... ..... ... ....... .... ... ...... . . .... ...... . . . . . .. . ....... ... .... ...... ... ..... ............ .. . . . . . . . . . . . ............ ......... ∆r r(t) Figure 13.2.1 r(t + ∆t) Approximating the derivative. Unfortunately, the vector ∆r approaches 0 in length; the vector h0, 0, 0i is not very informative. By dividing by ∆t, when it is small, we effectively keep magnifying the length 13.2 Calculus with vector functions 333 of ∆r so that in the limit it doesn’t disappear. Thus the limiting vector hf ′ (t), g ′ (t), h′ (t)i will (usually) be a good, non-zero vector that is tangent to the curve. What about the length of this vector? It’s nice that we’ve kept it away from zero, but what does it measure, if anything? Consider the length of one of the vectors that approaches the tangent vector: |r(t + ∆t) − r(t)| r(t + ∆t) − r(t) = ∆t |∆t| The numerator is the length of the vector that points from one position of the object to a “nearby” position; this length is approximately the distance traveled by the object between times t and t + ∆t. Dividing this distance by the length of time it takes to travel that distance gives the average speed. As ∆t approaches zero, this average speed approaches the actual, instantaneous speed of the object at time t. So by performing an “obvious” calculation to get something that looks like the derivative of r(t), we get precisely what we would want from such a derivative: the vector r′ (t) points in the direction of travel of the object and its length tells us the speed of travel. In the case that t is time, then, we call v(t) = r′ (t) the velocity vector. Even if t is not time, r′ (t) is useful—it is a vector tangent to the curve. We have seen that r = hcos t, sin t, ti is a helix. We compute p √ sin2 t + cos2 t + 1 = 2. So thinking of this as a r = h− sin t, cos t, 1i, and |r′ | = √ description of a moving object, its speed is always 2; see figure 13.2.2. EXAMPLE 13.2.1 ′ Figure 13.2.2 A tangent vector on the helix. (AP) EXAMPLE 13.2.2 The velocity vector for hcos t, sin t, cos ti is h− sin t, cos t, − sin ti. As before, the first two coordinates mean that from above this curve looks like a circle. The 334 Chapter 13 Vector Functions z coordinate is now also periodic, so that as the object moves around the curve its height oscillates up and down. In fact it turns out that the curve is a tilted ellipse, as shown in figure 13.2.3. Figure 13.2.3 The ellipse r = hcos t, sin t, cos ti. (AP) EXAMPLE 13.2.3 The velocity vector for hcos t, sin t, cos 2ti is h− sin t, cos t, −2 sin 2ti. The z coordinate is now oscillating twice as fast as in the previous example, so the graph is not surprising; see figure 13.2.4. Figure 13.2.4 hcos t, sin t, cos 2ti. (AP) EXAMPLE 13.2.4 Find the angle between the curves ht, 1−t, 3+t2 i and h3−t, t−2, t2 i where they meet. 13.2 Calculus with vector functions 335 The angle between two curves at a point is the angle between their tangent vectors— any tangent vectors will do, so we can use the derivatives. We need to find the point of intersection, evaluate the two derivatives there, and finally find the angle between them. To find the point of intersection, we need to solve the equations t=3−u 1−t=u−2 3 + t 2 = u2 Solving either of the first two equations for u and substituting in the third gives 3 + t2 = (3 − t)2 , which means t = 1. This together with u = 2 satisfies all three equations. Thus the two curves meet at (1, 0, 4), the first when t = 1 and the second when t = 2. The derivatives are h1, −1, 2ti and h−1, 1, 2ti; at the intersection point these are h1, −1, 2i and h−1, 1, 4i. The cosine of the angle between them is then cos θ = −1 − 1 + 8 1 √ √ = √ , 6 18 3 √ so θ = arccos(1/ 3) ≈ 0.96. The derivatives of vector functions obey some familiar looking rules, which we will occasionally need. THEOREM 13.2.5 Suppose r(t) and s(t) are differentiable functions, f (t) is a differentiable function, and a is a real number. d a. ar(t) = ar′ (t) dt d b. (r(t) + s(t)) = r′ (t) + s′ (t) dt d f (t)r(t) = f (t)r′ (t) + f ′ (t)r(t) c. dt d (r(t) · s(t)) = r′ (t) · s(t) + r(t) · s′ (t) d. dt d e. (r(t) × s(t)) = r′ (t) × s(t) + r(t) × s′ (t) dt d r(f (t)) = r′ (f (t))f ′ (t) f. dt Note that because the cross product is not commutative you must remember to do the three cross products in formula (e) in the correct order. 336 Chapter 13 Vector Functions When the derivative of a function f (t) is zero, we know that the function has a horizontal tangent line, and may have a local maximum or minimum point. If r′ (t) = 0, the geometric interpretation is quite different, though the interpretation in terms of motion is similar. Certainly we know that the object has speed zero at such a point, and it may thus be abruptly changing direction. In three dimensions there are many ways to change direction; geometrically this often means the curve has a cusp or a point, as in the path of a ball that bounces off the floor or a wall. EXAMPLE 13.2.6 Suppose that r(t) = h1 + t3 , t2 , 1i, so r′ (t) = h3t2 , 2t, 0i. This is 0 at t = 0, and there is indeed a cusp at the point (1, 0, 1), as shown in figure 13.2.5. 2.0 1.5 1.0 0.5 0.0 0.0 0.0 0.5 0.25 1.0 0.5 1.5 0.75 1.0 Figure 13.2.5 2.0 h1 + t3 , t2 , 1i has a cusp at h1, 0, 1i. (AP) Sometimes we will be interested in the direction of r′ but not its length. In some cases, we can still work with r′ , as when we find the angle between two curves. On other occasions it will be useful to work with a unit vector in the same direction as r′ ; of course, we can compute such a vector by dividing r′ by its own length. This standard unit tangent vector is usually denoted by T: r′ T= ′ . |r | In a sense, when we computed the angle between two tangent vectors we have already made use of the unit tangent, since cos θ = r′ · s′ r′ s′ = · |r′ ||s′ | |r′ | |s′ | 13.2 Calculus with vector functions 337 Now that we know how to make sense of r′ , we immediately know what an antiderivative must be, namely Z r(t) dt = h Z f (t) dt, Z g(t) dt, Z h(t) dti, if r = hf (t), g(t), h(t)i. What about definite integrals? Suppose that v(t) gives the velocity of an object at time t. Then v(t)∆t is a vector that approximates the displacement of the object over the time ∆t: v(t)∆t points in the direction of travel, and |v(t)∆t| = |v(t)||∆t| is the speed of the object times ∆t, which is approximately the distance traveled. Thus, if we sum many such tiny vectors: n−1 X v(ti )∆t i=0 we get an approximation to the displacement vector over the time interval [t0 , tn ]. If we take the limit we get the exact value of the displacement vector: lim n−1 X v(ti )∆t = i=0 Z tn t0 v(t) dt = r(tn ) − r(t0 ). Denote r(t0 ) by r0 . Then given the velocity vector we can compute the vector function r giving the location of the object: r(t) = r0 + Z t v(u) du. t0 EXAMPLE 13.2.7 An object moves with velocity vector hcos t, sin t, cos ti, starting at (1, 1, 1) at time 0. Find the function r giving its location. r(t) = h1, 1, 1i + Z 0 t hcos u, sin u, cos ui du t = h1, 1, 1i + hsin u, − cos u, sin ui|0 = h1, 1, 1i + hsin t, − cos t, sin ti − h0, −1, 0i = h1 + sin t, 2 − cos t, 1 + sin ti See figure 13.2.6. 338 Chapter 13 Vector Functions Figure 13.2.6 Path of the object with its initial velocity vector. (AP) Exercises 13.2. 1. Find r′ and T for r = ht2 , 1, ti. ⇒ 2. Find r′ and T for r = hcos t, sin 2t, t2 i. ⇒ 3. Find r′ and T for r = hcos(et ), sin(et ), sin ti. ⇒ 4. Find a vector function for the line tangent to the helix hcos t, sin t, ti when t = π/4. ⇒ 5. Find a vector function for the line tangent to hcos t, sin t, cos 4ti when t = π/3. ⇒ 6. Find the cosine of the angle between the curves h0, t2 , ti and hcos(πt/2), sin(πt/2), ti where they intersect. ⇒ 7. Find the cosine of the angle between the curves hcos t, − sin(t)/4, sin ti and hcos t, sin t, sin(2t)i where they intersect. ⇒ 8. Suppose that |r(t)| = k, for some constant k. This means that r describes some path on the sphere of radius k with center at the origin. Show that r is perpendicular to r′ at every point. Hint: Use Theorem 13.2.5, part (d). 9. A bug is crawling along the spoke of a wheel that lies along a radius of the wheel. The bug is crawling at 1 unit per second and the wheel is rotating at 1 radian per second. Suppose the wheel lies in the y-z plane with center at the origin, and at time t = 0 the spoke lies along the positive y axis and the bug is at the origin. Find a vector function r(t) for the position of the bug at time t, the velocity vector r′ (t), the unit tangent T(t), and the speed of the bug |r′ (t)|. ⇒ 10. An object moves with velocity vector hcos t, sin t, ti, starting at h0, 0, 0i when t = 0. Find the function r giving its location. ⇒ 11. An object moves with velocity vector ht, t2 , −ti, starting at h1, 2, 3i when t = 0. Find the function r giving its location. ⇒ 12. An object moves with velocity vector ht2 , sin t, cos ti, starting at h−1, 1, 2i when t = 1. Find the function r giving its location. ⇒ 13. The position function of a particle is given by r(t) = ht2 , 5t, t2 − 16ti, t ≥ 0. When is the speed of the particle a minimum? ⇒ 13.3 Arc length and curvature 339 14. A particle moves so that its position is given by hcos t, sin t, cos(6t)i. Find the maximum and minimum speeds of the particle. ⇒ 15. An object moves with velocity vector ht, t2 , cos ti, starting at h0, 0, 0i when t = 0. Find the function r giving its location. ⇒ 16. What is the physical interpretation of the dot product of two vector valued functions? What is the physical interpretation of the cross product of two vector valued functions? 17. Show, using the rules of cross products and differentiation, that d (r(t) × r′ (t)) = r(t) × r′′ (t). dt 18. Determine the point at which f (t) = ht, t2 , t3 i and g(t) = hcos(t), cos(2t), t + 1i intersect, and find the angle between the curves at that point. (Hint: You’ll need to set this one up like a line intersection problem, writing one in s and one in t.) If these two functions were the trajectories of two airplanes on the same scale of time, would the planes collide at their point of intersection? Explain. ⇒ 19. Find the equation of the plane perpendicular to the curve r(t) = h2 sin(3t), t, 2 cos(3t)i at the point (0, π, −2). ⇒ 20. Find the equation of the plane perpendicular to hcos t, sin t, cos(6t)i when t = π/4. ⇒ 21. At what point on the curve r(t) = ht3 , 3t, t4 i is the plane perpendicular to the curve also parallel to the plane 6x + 6y − 8z = 1? ⇒ 22. Find the equation of the line tangent to hcos t, sin t, cos(6t)i when t = π/4. ⇒ 13.3 Ar length and urvature Sometimes it is useful to compute the length of a curve in space; for example, if the curve represents the path of a moving object, the length of the curve between two points may be the distance traveled by the object between two times. Recall that if the curve is given by the vector function r then the vector ∆r = r(t + ∆t) − r(t) points from one position on the curve to another, as depicted in figure 13.2.1. If the points are close together, the length of ∆r is close to the length of the curve between the two points. If we add up the lengths of many such tiny vectors, placed head to tail along a segment of the curve, we get an approximation to the length of the curve over that segment. In the limit, as usual, this sum turns into an integral that computes precisely the length of the curve. First, note that |∆r| = |∆r| ∆t ≈ |r′ (t)| ∆t, ∆t when ∆t is small. Then the length of the curve between r(a) and r(b) is lim n→∞ n−1 X i=0 |∆r| = lim n→∞ n−1 X i=0 Z b n−1 X |∆r| ′ |r′ (t)| dt. |r (t)| ∆t = ∆t = lim n→∞ ∆t a i=0 340 Chapter 13 Vector Functions (Well, sometimes. This works if between a and b the segment of curve is traced out exactly once.) Let’s find the length of one turn of the helix r = hcos t, sin t, ti (see p √ figure 13.1.1). We compute r′ = h− sin t, cos t, 1i and |r′ | = sin2 t + cos2 t + 1 = 2, so the length is Z 2π √ √ 2 dt = 2 2π. EXAMPLE 13.3.1 0 EXAMPLE 13.3.2 Suppose y = ln x; what is the length of this curve between x = 1 √ and x = 3? Although this problem does not appear to involve vectors or three dimensions, we can interpret it in those terms: let r(t) = ht, ln t, 0i. This vector function traces out precisely p y = ln x in the x-y plane. Then r′ (t) = h1, 1/t, 0i and |r′ (t)| = 1 + 1/t2 and the desired length is Z √3 r √ √ 1 1 1 + 2 dt = 2 − 2 + ln( 2 + 1) − ln 3. t 2 1 (This integral is a bit tricky, but requires only methods we have learned.) Notice that there is nothing special about y = ln x, except that the resulting integral can be computed. In general, given any y = f (x), we can think of this as the vector p function r(t) = ht, f (t), 0i. Then r′ (t) = h1, f ′ (t), 0i and |r′ (t)| = 1 + (f ′ )2 . The length of the curve y = f (x) between a and b is thus Z a b p 1 + (f ′ (x))2 dx. Unfortunately, such integrals are often impossible to do exactly and must be approximated. One useful application of arc length is the arc length parameterization. A vector function r(t) gives the position of a point in terms of the parameter t, which is often time, but need not be. Suppose s is the distance along the curve from some fixed starting point; if we use s for the variable, we get r(s), the position in space in terms of distance along the curve. We might still imagine that the curve represents the position of a moving object; now we get the position of the object as a function of how far the object has traveled. EXAMPLE 13.3.3 Suppose r(t) = hcos t, sin t, 0i. We know that this curve is a circle of radius 1. While t might represent time, it can also in this case represent the usual angle between the positive x-axis and r(t). The distance along the circle from (1, 0, 0) to 13.3 Arc length and curvature 341 (cos t, sin t, 0) is also t—this is the definition of radian measure. Thus, in this case s = t and r(s) = hcos s, sin s, 0i. EXAMPLE 13.3.4 Suppose r(t) = hcos t, sin t, ti. We know that this curve is a helix. The distance along the helix from (1, 0, 0) to (cos t, sin t, t) is s= Z t ′ 0 |r (u)| du = Z tp 2 cos2 u + sin u + 1 du = 0 Z t√ 2 du = √ 2t. 0 √ Thus, the value of t that gets us distance s along the helix is t = s/ 2, and so the same √ √ √ curve is given by r̂(s) = hcos(s/ 2), sin(s/ 2), s/ 2i. In general, if we have a vector function r(t), to convert it to a vector function in terms of arc length we compute Z t |r′ (u)| du = f (t), s= a solve s = f (t) for t, getting t = g(s), and substitute this back into r(t) to get r̂(s) = r(g(s)). Suppose that t is time. By the Fundamental Theorem of Calculus, if we start with arc length Z t s(t) = |r′ (u)| du a and take the derivative, we get s′ (t) = |r′ (t)|. Here s′ (t) is the rate at which the arc length is changing, and we have seen that |r′ (t)| is the speed of a moving object; these are of course the same. Suppose that r(s) is given in terms of arc length; what is |r′ (s)|? It is the rate at which arc length is changing relative to arc length; it must be 1! In the case of the helix, for √ √ √ example, the arc length parameterization is hcos(s/ 2), sin(s/ 2), s/ 2i, the derivative √ √ √ √ √ is h− sin(s/ 2)/ 2, cos(s/ 2)/ 2, 1/ 2i, and the length of this is s r √ √ sin2 (s/ 2) cos2 (s/ 2) 1 1 1 + + = + = 1. 2 2 2 2 2 So in general, r′ is a unit tangent vector. Given a curve r(t), we would like to be able to measure, at various points, how sharply curved it is. Clearly this is related to how “fast” a tangent vector is changing direction, so a first guess might be that we can measure curvature with |r′′ (t)|. A little thought shows that this is flawed; if we think of t as time, for example, we could be tracing out the curve 342 Chapter 13 Vector Functions more or less quickly as time passes. The second derivative |r′′ (t)| incorporates this notion of time, so it depends not simply on the geometric properties of the curve but on how quickly we move along the curve. EXAMPLE 13.3.5 Consider r(t) = hcos t, sin t, 0i and s(t) = hcos 2t, sin 2t, 0i. Both of these vector functions represent the unit circle in the x-y plane, but if t is interpreted as time, the second describes an object moving twice as fast as the first. Computing the second derivatives, we find |r′′ (t)| = 1, |s′′ (t)| = 4. To remove the dependence on time, we use the arc length parameterization. If a curve is given by r(s), then the first derivative r′ (s) is a unit vector, that is, r′ (s) = T(s). We now compute the second derivative r′′ (s) = T′ (s) and use |T′ (s)| as the “official” measure of curvature, usually denoted κ. EXAMPLE 13.3.6 We have seen that the arc length parameterization of a particu√ √ √ lar helix is r(s) = hcos(s/ 2), sin(s/ 2), s/ 2i. Computing the second derivative gives √ √ r′′ (s) = h− cos(s/ 2)/2, − sin(s/ 2)/2, 0i with length 1/2. What if we are given a curve as a vector function r(t), where t is not arc length? We have seen that arc length can be difficult to compute; fortunately, we do not need to convert to the arc length parameterization to compute curvature. Instead, let us imagine that we have done this, so we have found t = g(s) and then formed r̂(s) = r(g(s)). The first derivative r̂′ (s) is a unit tangent vector, so it is the same as the unit tangent vector T(t) = T(g(s)). Taking the derivative of this we get dt d T(g(s)) = T′ (g(s))g ′(s) = T′ (t) . ds ds The curvature is the length of this vector: κ = |T′ (t)|| |T′ (t)| |T′ (t)| dt |= = ′ . ds |ds/dt| |r (t)| (Recall that we have seen that ds/dt = |r′ (t)|.) Thus we can compute the curvature by computing only derivatives with respect to t; we do not need to do the conversion to arc length. EXAMPLE 13.3.7 Returning to the helix, suppose we start with the parameteri√ zation r(t) = hcos t, sin t, ti. Then r′ (t) = h− sin t, cos t, 1i, |r′ (t)| = 2, and T(t) = √ √ √ h− sin t, cos t, 1i/ 2. Then T′ (t) = h− cos t, − sin t, 0i/ 2 and |T′ (t)| = 1/ 2. Finally, √ √ κ = 1/ 2/ 2 = 1/2, as before. 13.3 Arc length and curvature 343 EXAMPLE 13.3.8 Consider this circle of radius a: r(t) = ha cos t, a sin t, 1i. Then r′ (t) = h−a sin t, a cos t, 0i, |r′ (t)| = a, and T(t) = h−a sin t, a cos t, 0i/a. Now T′ (t) = h−a cos t, −a sin t, 0i/a and |T′ (t)| = 1. Finally, κ = 1/a: the curvature of a circle is everywhere the reciprocal of the radius. It is sometimes useful to think of curvature as describing what circle a curve most resembles at a point. The curvature of the helix in the previous example is 1/2; this means that a small piece of the helix looks very much like a circle of radius 2, as shown in figure 13.3.1. Figure 13.3.1 A circle with the same curvature as the helix. (AP) EXAMPLE 13.3.9 Consider r(t) = hcos t, sin t, cos 2ti, as shown in figure 13.2.4. r′ (t) = q ′ h− sin t, cos t, −2 sin(2t)i and |r (t)| = 1 + 4 sin2 (2t), so T(t) = * − sin t cos t −2 sin 2t q ,q ,q 1 + 4 sin2 (2t) 1 + 4 sin2 (2t) 1 + 4 sin2 (2t) + . Computing the derivative of this and then the length of the resulting vector is possible but unpleasant. Fortunately, there is an alternate formula for the curvature that is often simpler than the one we have: |r′ (t) × r′′ (t)| . κ= |r′ (t)|3 EXAMPLE 13.3.10 Returning to the previous example, we compute the second derivative r′′ (t) = h− cos t, − sin t, −4 cos(2t)i. Then the cross product r′ (t) × r′′ (t) is h−4 cos t cos 2t − 2 sin t sin 2t, 2 cos t sin 2t − 4 sin t cos 2t, 1i. 344 Chapter 13 Vector Functions Computing the length of this vector and dividing by |r′ (t)|3 is still a bit tedious. With the aid of a computer we get √ 48 cos4 t − 48 cos2 t + 17 . κ= (−16 cos4 t + 16 cos2 t + 1)3/2 Graphing this we get 4 2 . . . . ... .. .... .... .... ... .. ..... ..... ..... .. .. .. . .. . .. . ... ... ... .... ... .... ... .... ... . .. .. .. .. .. ... .. ... .. ... .. ... .. .. ... .. ... .. ... ... ... ... .... ... .... ... .... . . . . . . . .. ... ... ... ... ... ... ... ... .. .. .... .. .... .. .... ... . . . . . . . . . . . .. .. .. .. ... .. ... .. ... ... .. .. ... .. ... .. ... .. .. .. .. .. .. .. .. ... .. .. ... .. ... .. ... . . . . ... ... ... ... .. .. .. .. ... ... ... ... .. .. .. .. ... ... ... ... .. .. .. .. ... ... ... ... ... ... ... ... ... ... ... ... .. .. .. ... .. .. .. .. .. .. .. .. ... ... .. .. .. .. .. .. . . . . . .. ... . . . .. .. . .. .. ... ... .. .. .. .. ... ... ... ... .. ... ... ... ... ... ... ... .... ... ... ... .... ...... . . . . . . . . . . . . . . . . . . . ......... ......... ......... ......... π 2 π 3π 2 2π Compare this to figure 13.2.4—you may want to load the Java applet there so that you can see it from different angles. The highest curvature occurs where the curve has its highest and lowest points, and indeed in the picture these appear to be the most sharply curved portions of the curve, while the curve is almost a straight line midway between those points. Let’s see why this alternate formula is correct. Starting with the definition of T, r = |r′ |T so by the product rule r′′ = |r′ |′ T + |r′ |T′ . Then by Theorem 12.4.2 the cross product is r′ × r′′ = |r′ |T × (|r′ |′ T + |r′ |T′ ) ′ = |r′ |T × |r′ |′ T + |r′ |T × |r′ |T′ = |r′ ||r′ |′ (T × T) + |r′ |2 (T × T′ ) = |r′ |2 (T × T′ ) because T × T = 0, since T is parallel to itself. Then |r′ × r′′ | = |r′ |2 |T × T′ | = |r′ |2 |T||T′ | sin θ = |r′ |2 |T′ | using exercise 8 in section 13.2 to see that θ = π/2. Dividing both sides by |r′ |3 then gives the desired formula. We used the fact here that T′ is perpendicular to T; the vector N = T′ /|T′ | is thus a unit vector perpendicular to T, called the unit normal to the curve. Occasionally of use is the unit binormal B = T × N, a unit vector perpendicular to both T and N. 13.4 Motion along a curve 345 Exercises 13.3. 1. Find the length of h3 cos t, 2t, 3 sin ti, t ∈ [0, 2π]. ⇒ 2. Find the length of ht2 , 2, t3 i, t ∈ [0, 1]. ⇒ 3. Find the length of ht2 , sin t, cos ti, t ∈ [0, 1]. ⇒ 4. Find the length of the curve y = x3/2 , x ∈ [1, 9]. ⇒ 5. Set up an integral to compute the length of hcos t, sin t, et i, t ∈ [0, 5]. (It is tedious but not too difficult to compute this integral.) ⇒ 6. Find the curvature of ht, t2 , ti. ⇒ 7. Find the curvature of ht, t2 , t2 i. ⇒ 8. Find the curvature of ht, t2 , t3 i. ⇒ 9. Find the curvature of y = x4 at (1, 1). ⇒ 13.4 Motion along a urve We have already seen that if t is time and an object’s location is given by r(t), then the derivative r′ (t) is the velocity vector v(t). Just as v(t) is a vector describing how r(t) changes, so is v′ (t) a vector describing how v(t) changes, namely, a(t) = v′ (t) = r′′ (t) is the acceleration vector. EXAMPLE 13.4.1 Suppose r(t) = hcos t, sin t, 1i. Then v(t) = h− sin t, cos t, 0i and a(t) = h− cos t, − sin t, 0i. This describes the motion of an object traveling on a circle of radius 1, with constant z coordinate 1. The velocity vector is of course tangent to the curve; note that a · v = 0, so v and a are perpendicular. In fact, it is not hard to see that a points from the location of the object to the center of the circular path at (0, 0, 1). Recall that the unit tangent vector is given by T(t) = v(t)/|v(t)|, so v = |v|T. If we take the derivative of both sides of this equation we get a = |v|′ T + |v|T′ . (13.4.1) Also recall the definition of the curvature, κ = |T′ |/|v|, or |T′ | = κ|v|. Finally, recall that we defined the unit normal vector as N = T′ /|T′ |, so T′ = |T′ |N = κ|v|N. Substituting into equation 13.4.1 we get a = |v|′ T + κ|v|2 N. (13.4.2) The quantity |v(t)| is the speed of the object, often written as v(t); |v(t)|′ is the rate at which the speed is changing, or the scalar acceleration of the object, a(t). Rewriting equation 13.4.2 with these gives us a = aT + κv 2 N = aT T + aN N; aT is the tangential component of acceleration and aN is the normal component of acceleration. We have already seen that aT measures how the speed is changing; if 346 Chapter 13 Vector Functions you are riding in a vehicle with large aT you will feel a force pulling you into your seat. The other component, aN , measures how sharply your direction is changing with respect to time. So it naturally is related to how sharply the path is curved, measured by κ, and also to how fast you are going. Because aN includes v 2 , note that the effect of speed is magnified; doubling your speed around a curve quadruples the value of aN . You feel the effect of this as a force pushing you toward the outside of the curve, the “centrifugal force.” In practice, if want aN we would use the formula for κ: aN = κ|v|2 = |r′ × r′′ | ′ 2 |r′ × r′′ | |r | = . |r′ |3 |r′ | To compute aT we can project a onto v: aT = v·a r′ · r′′ = . |v| |r′ | EXAMPLE 13.4.2 Suppose r = ht, t2 , t3 i. Compute v, a, aT , and aN . Taking derivatives we get v = h1, 2t, 3t2 i and a = h0, 2, 6ti. Then 4t + 18t3 aT = √ 1 + 4t2 + 9t4 and aN √ 4 + 36t2 + 36t4 = √ . 1 + 4t2 + 9t4 Exercises 13.4. 1. Let r = hcos t, sin t, ti. Compute v, a, aT , and aN . ⇒ 2. Let r = hcos t, sin t, t2 i. Compute v, a, aT , and aN . ⇒ 3. Let r = hcos t, sin t, et i. Compute v, a, aT , and aN . ⇒ 4. Let r = ht2 , 2t − 3, 3t2 − 3ti. Compute v, a, aT , and aN . ⇒ 5. Let r = het , sin t, et i. Compute v, a, aT , and aN . ⇒ 6. Suppose an object moves so that its acceleration is given by a = h−3 cos t, −2 sin t, 0i. At time t = 0 the object is at (3, 0, 0) and its velocity vector is h0, 2, 0i. Find v(t) and r(t) for the object. ⇒ 7. Suppose an object moves so that its acceleration is given by a = h−3 cos t, −2 sin t, 0i. At time t = 0 the object is at (3, 0, 0) and its velocity vector is h0, 2.1, 0i. Find v(t) and r(t) for the object. ⇒ 8. Suppose an object moves so that its acceleration is given by a = h−3 cos t, −2 sin t, 0i. At time t = 0 the object is at (3, 0, 0) and its velocity vector is h0, 2, 1i. Find v(t) and r(t) for the object. ⇒ 9. Suppose an object moves so that its acceleration is given by a = h−3 cos t, −2 sin t, 0i. At time t = 0 the object is at (3, 0, 0) and its velocity vector is h0, 2.1, 1i. Find v(t) and r(t) for the object. ⇒ 13.4 Motion along a curve 347 10. Describe a situation in which the normal component of acceleration is 0 and the tangential component of acceleration is non-zero. Is it possible for the tangential component of acceleration to be 0 while the normal component of acceleration is non-zero? Explain. Finally, is it possible for an object to move (not be stationary) so that both the tangential and normal components of acceleration are 0? Explain.