Calculus/Vector Functions

From testwiki
Jump to navigation Jump to search

Introduction

If we have a function f:n , we say that f's image (the set {f(t)|t} - or some subset of ) is a curve in n and f is its parametrization.

Parameterizations are not necessarily unique - for example, f(t)=(cos(t),sin(t)) such that t[0,2π) is one parametrization of the unit circle, and g(t)=(cos(at),sin(at)) such that t[0,2πa) is a whole family of parameterizations of that circle.

Collision and intersection points

Say we have two different curves. It may be important to consider

  • points the two curves share - where they intersect
  • intersections which occur for the same value of t - where they collide.

Intersection points

Firstly, we have two parameterizations f(t) and g(t) , and we want to find out when they intersect, this means that we want to know when the function values of each parametrization are the same. This means that we need to solve

f(t)=g(s)

because we're seeking the function values independent of the times they intersect.

For example, if we have f(t)=(t,3t) and g(t)=(t,t2) , and we want to find intersection points:

f(t)=g(s)
(t,3t)=(s,s2)
t=s , 3t=s2

with solutions (t,s)=(0,0) , (3,9)

So, the two curves intersect at the points (0,0) , (3,9) .

Collision points

However, if we want to know when the points "collide", with f(t) and g(t) , we need to know when both the function values and the times are the same, so we need to solve instead

f(t)=g(t)

For example, using the same functions as before, f(t)=(t,3t) and g(t)=(t,t2) , and we want to find collision points:

f(t)=g(t)
(t,3t)=(t,t2)
t=t , 3t=t2

which gives solutions t=0,3 So the collision points are (0,0) , (3,9) .

We may want to do this to actually model physical problems, such as in ballistics.

Intersection vector functions

We can also use vector functions to represent the curve of intersection of two surfaces. For example, we want to know the curve of intersection of the cylinder

x2+y2=1

and the plane

y+z=2

.

Vector functions rely on parameterizations, so we can rewrite the equation of the cylinder into:

cos2t+sin2t=1

, where

x=cost,y=sint

.


From the equation of the plane, we know that z=2y=2sint. Thus the corresponding vector equation is:


𝐫(t)=cost,sint,2sint

Limits and Continuity

The limit of a vector function 𝐫 is defined by taking the limits of its component functions. Template:Calculus/Def

And the requirement for continuity is also simple:

A vector function

𝐫

is continuous at

a

if

limta𝐫(t)=𝐫(a)

.

Derivatives and Integrals

Differentiation

Recall that the first derivative of a scalar function

f(x)

is defined as:

f(x)=limΔx0f(x+Δx)f(x)Δx

The first derivative of a vector function

𝐫(t)

is defined in much the same way:

𝐫(t)=limΔt0𝐫(t+Δt)𝐫(t)Δt

We can use this definition to prove that the derivative of a vector function can be presented as the derivative of its component functions.

𝐫(t)=limΔt01Δt[𝐫(t+Δt)𝐫(t)]definition=limΔt01Δt[f(t+Δt),g(t+Δt),h(t+Δt)f(t),g(t),h(t)]component form=limΔt0f(t+Δt)f(t)Δt,g(t+Δt)g(t)Δt,h(t+Δt)h(t)Δtvector addition and multiplication=limΔt0f(t+Δt)f(t)Δt,limΔt0g(t+Δt)g(t)Δt,limΔt0h(t+Δt)h(t)Δtdefinition=f(t),g(t),h(t)definition

Thus, using the same method, we can derive the second derivative and so on. Template:Calculus/DefSo the

n

th-order derivative should look like this:

𝐫(n)(t)=f(n)(t),g(n)(t),h(n)(t)

Differentiation rules

Just like real-valued functions, there are some differentiation rules in the world of vector functions. The factor that makes vector differentiation rules slight more complicated is the product rule because there are two kinds of multiplication in vectors: dot product and cross product. Template:Calculus/Def Naturally, we will prove that those rules are correct. Let us assume that 𝐮(t)=f1(t),f2(t),f3(t) and 𝐯(t)=g1(t),g2(t),g3(t).


Rule 1: the addition rule

ddt[𝐮(t)+𝐯(t)]=ddt[f1(t),f2(t),f3(t)+g1(t),g2(t),g3(t)]component form=ddtf1(t)+g1(t),f2(t)+g2(t),f3(t)+g3(t)vector addition=ddt(f1(t)+g1(t)),ddt(f2(t)+g2(t)),ddt(f3(t)+g3(t))distribution=f1(t)+g1(t),f2(t)+g2(t),f3(t)+g3(t)addition rule for real-valued functions=f1(t),f2(t),f3(t)+g1(t),g2(t),g3(t)vector addition=𝐮(t)+𝐯(t)vector form

Rule 2: the scalar multiplication rule

ddt[f(t)𝐮(t)]=ddt[f(t)(f1(t)f2(t)f3(t))]component form=ddt(f(t)f1(t)f(t)f2(t)f(t)f3(t))scalar multiplication=(ddt(f(t)f1(t))ddt(f(t)f2(t))ddt(f(t)f3(t)))distribution=(f(t)f1(t)+f(t)f1(t)f(t)f2(t)+f(t)f2(t)f(t)f3(t)+f(t)f3(t))multiplication rule for real-valued functions=(f(t)f1(t)f(t)f2(t)f(t)f3(t))+(f(t)f1(t)f(t)f2(t)f(t)f3(t))vector addition=f(t)(f1(t)f2(t)f3(t))+f(t)(f1(t)f2(t)f3(t))scalar multiplication=f(t)𝐮(t)+f(t)𝐮(t)vector form

Rule 3: the dot product rule

ddt[𝐮(t)𝐯(t)]=ddt[f1(t),f2(t),f3(t)g1(t),g2(t),g3(t)]component form=ddt(f1(t)g1(t)+f2(t)g2(t)+f3(t)g3(t))dot product=ddti=13fi(t)gi(t)simplification=i=13ddtfi(t)gi(t)distribution=i=13[fi(t)gi(t)+fi(t)gi(t)]multiplication rule for real-valued functions=i=13fi(t)gi(t)+i=13fi(t)gi(t)vector addition=𝐮(t)𝐯(t)+𝐮(t)𝐯(t)vector form

Rule 4: the cross product rule

ddt[𝐮(t)×𝐯(t)]=ddt[f1(t),f2(t),f3(t)×g1(t),g2(t),g3(t)]component form=ddt|𝐢^𝐣^𝐤^f1(t)f2(t)f3(t)g1(t)g2(t)g3(t)|cross product=ddt(f2(t)g3(t)f3(t)g2(t)f3(t)g1(t)f1(t)g3(t)f1(t)g2(t)f2(t)g1(t))=(ddt(f2(t)g3(t)f3(t)g2(t))ddt(f3(t)g1(t)f1(t)g3(t))ddt(f1(t)g2(t)f2(t)g1(t)))distribution=(f2(t)g3(t)+f2(t)g3(t)f3(t)g2(t)f3(t)g2(t)f3(t)g1(t)+f3(t)g1(t)f1(t)g3(t)f1(t)g3(t)f2(t)g3(t)+f2(t)g3(t)f3(t)g2(t)f3(t)g2(t))multiplication rule for real-valued functions=(f2(t)g3(t)f3(t)g2(t)+f2(t)g3(t)f3(t)g2(t)f3(t)g1(t)f1(t)g3(t)+f3(t)g1(t)f1(t)g3(t)f2(t)g3(t)f3(t)g2(t)+f2(t)g3(t)f3(t)g2(t))rearrangement=(f2(t)g3(t)f3(t)g2(t)f3(t)g1(t)f1(t)g3(t)f2(t)g3(t)f3(t)g2(t))+(f2(t)g3(t)f3(t)g2(t)f3(t)g1(t)f1(t)g3(t)f2(t)g3(t)f3(t)g2(t))vector addtion=|𝐢^𝐣^𝐤^f1(t)f2(t)f3(t)g1(t)g2(t)g3(t)|+|𝐢^𝐣^𝐤^f1(t)f2(t)f3(t)g1(t)g2(t)g3(t)|cross product=(f1(t)f2(t)f3(t))×(g1(t)g2(t)g3(t))+(f1(t)f2(t)f3(t))×(g1(t)g2(t)g3(t))=𝐮(t)×𝐯(t)+𝐮(t)×𝐯(t)vector form

Rule 5: the chain rule

ddt[𝐮(f(t))]=ddt(f1(f(t))f2(f(t))f3(f(t)))component form=(ddtf1(f(t))ddtf2(f(t))ddtf3(f(t)))distribution=(f(t)f1(f(t))f(t)f2(f(t))f(t)f3(f(t)))chain rule for real-valued functions=f(t)(f1(f(t))f2(f(t))f3(f(t)))scalar multiplication=f(t)𝐮(f(t))vector form

Angle between curves

We can then formulate the concept of the angle between two curves by considering the angle between the two tangent vectors. If two curves, parametrized by f1 and f2 intersect at some point, which means that

f1(s)=f2(t)=c

the angle between these two curves at c is the angle between the tangent vectors f1(s) and f2(t), using the dot product, is given by

arccos(f1(s)f2(t)f1(s)f2(t))

Integration

Similar to real-valued functions, the definite integral of a vector function

𝐫(t)

is defined as:

ab𝐫(t)dt=limni=1n𝐫(ti*)Δt=limn[(i=1nf(ti*)Δt)𝐢+(i=1ng(ti*)Δt)𝐣+(i=1nh(ti*)Δt)𝐤]=(abf(t)dt)𝐢+(abg(t)dt)𝐣+(abh(t)dt)𝐤

We can extend the Fundamental Theorem of Calculus to continuous vector functions as follows:

Template:Calculus/Def

For indefinite integrals, the definition is:

𝐫(t)dt=(f(t)dt)𝐢+(g(t)dt)𝐣+(h(t)dt)𝐤+𝐂

Arc length

Recall that in Chapter Template:Calculus/map page, we deduced the length of a curve with parametric equations

{x=f(t)y=g(t)

,

atb

should be:

L=ab(dxdt)2+(dydt)2dt

Since vector functions are fundamentally parametric equations with directions, we can utilize the formula above into the length of a space curve.

Template:Calculus/Def

For those who prefer simplicity, the formula can be rewritten into:

L=ab|𝐫(t)|dt

or

dLdt=|𝐫(t)|

Reparametrization

Suppose that there is a vector function which describes the displacement of a particle with respect to time and has the equation:

𝐫(t)=cost 𝐢+sint 𝐣+t 𝐤

However, for some reason, we do not want to know the displacement of this particle with respect to time. Instead, we want to know its displacement with respect to its traveled distance (

s

) from

(1,0,0)

in the direction of increasing

t

. In order to do so, we need to find a way to describe time as a function of distance. In other words, we need to find

t=t(s)

. We can use the formula for the arc length to establish the relationship between time and distance because arc length in this case describes the distance traveled by the particle.


Before we start calculating, we need to introduce the arc length function. Template:Calculus/Def

According to the definition, the arc length function for our curve

𝐫(t)

should be

s(t)=0t|𝐫(u)|du=0t(sinu)2+(cosu)2+12du=2u]0t=2t

Note that a=0 because the initial point (1,0,0) corresponds to the parameter value t=0. Since it is in the direction of increasing t, the integration direction should be from 0 to t.

t(s)=s2

Then we substitute the original function and get the answer:

𝐫(s)=coss2 𝐢+sins2 𝐣+s2 𝐤

Reparametrization has important applications in real life because sometimes we want to know a value with respect to different variables. In this case, instead of describing the path of a particle with respect to time, we described its path with respect to its distance, which will be very useful in certain situations.

Curvature

Terminology

Before we start discussing curvature, there are some important vectors and concepts we need to be at least aware of.

The unit tangent vector

In the differentiation section of this chapter, we discussed the derivatives of a vector function. We know that

𝐯(t)=𝐫(t)

at

t=t0

is tangent to the curve

𝐫(t)

at

t=t0

.

𝐫(t)

is called the tangent vector. The unit tangent vector, however, eliminates the aspect of magnitude because it is defined as:

𝐓(t)=𝐫(t)|𝐫(t)|

As we can see, the magnitude of the unit tangent vector is always

1

. We can imagine that

𝐫(t)

as the displacement of a particle with respect to time. So, the unit tangent vector can be perceived as the direction of the velocity of the particle with respect to time. It can also be perceived as the direction of the tangential acceleration of the particle with respect to time. We will discuss motion in space in the next section, but this is a useful method to intuitively understand some vectors.

The unit normal vector

The unit normal vector is defined as

𝐍(t)=𝐓(t)|𝐓(t)|

The unit normal is orthogonal to the unit tangent because since

|𝐓(t)|=1

, we can get that:

ddt|𝐓(t)|2=0=ddt[𝐓(t)𝐓(t)]=2𝐓(t)𝐓(t)


𝐓(t)𝐓(t)=0

This means that

𝐓(t)

is orthogonal to

𝐓(t)

. Therefore,

𝐍(t)

is orthogonal to

𝐓(t)

. We can imagine that the unit normal vector is the direction of the normal acceleration of the particle with respect to time.

The binormal vector

The binormal vector is defined as

𝐁(t)=𝐓(t)×𝐍(t)

The binormal vector is perpendicular to both the unit tangent and the unit normal because of the properties of the cross product. The magnitude of the binormal is always 1 because

|𝐁(t)|=|𝐓(t)||𝐍(t)|sinθ=1

The normal plane, osculating plane, and the osculating circle

  • The normal plane is the plane determined by the normal and binormal vectors 𝐍 and 𝐁. The normal plane consists of all lines that are orthogonal to the tangent vector 𝐓.
  • The osculating plane is the plane determined by the unit tangent and unit normal 𝐓 and 𝐍. It is the plane that comes closest to containing the part of the curve near a point where t=t0.
  • The osculating circle is the circle that lies in the osculating plane towards the direction of 𝐍 with a radius r=1κ (the inverse of the curvature, which we will immediately discuss after this). It best describes how the curve behaves near the point where t=t0 because it shares the same tangent, normal, and curvature at that point.

These concepts are very important in the branch of differential geometry and in its applications to the motion of spacecraft.

Curvature

The curvature of a curve at a given point is a measure of how quickly the curve changes direction at that point. We define it to be the magnitude of the rate of change of the unit tangent with respect to arc length. We use arc length so that the curvature will be independent of parametrization.

Suppose that a space curve has vector function

𝐫

, unit tangent vector

𝐓

, and arc length

s

. The curvature of this curve is:

κ=|d𝐓(s)ds|

.

There are two other ways to express the curvature. We can express curvature in terms of

t

instead of

s

by utilizing the chain rule (recall that

dLdt=|𝐫(t)|

):

κ=|d𝐓ds|=|d𝐓dsdsdtdsdt|=|d𝐓(t)dtdsdt|=|𝐓(t)||𝐫(t)|

The third way is more complicated to deduce, but it is often more convenient to apply because it only requires

𝐫(t)

and its derivatives.

κ(t)=|𝐫(t)×𝐫(t)||𝐫(t)|3

And now for the proof for this theorem:

According to the definition of the unit tangent vector, we know that

𝐫=|𝐫|𝐓

. So the second derivative of

𝐫

should be:

𝐫=(|𝐫| 𝐓)=|𝐫| 𝐓+|𝐫| 𝐓the product rule

Now we calculate

𝐫×𝐫

.

𝐫×𝐫=(|𝐫|𝐓)×(|𝐫| 𝐓+|𝐫| 𝐓)substitution=(|𝐫|𝐓×|𝐫|𝐓)+(|𝐫|𝐓×|𝐫|𝐓)distribution=|𝐫||𝐫|(𝐓×𝐓)+|𝐫|2(𝐓×𝐓)rearrangement=|𝐫|2(𝐓×𝐓)realizing that 𝐓×𝐓=0

And then we calculate

|𝐫×𝐫|

.

|𝐫×𝐫|=|𝐫|2|𝐓×𝐓|substitution=|𝐫|2|𝐓||𝐓|sinθthe magnitude for the cross product=|𝐫|2|𝐓|when you realize |𝐓|=1 and 𝐓𝐓

We rearrange the equation into:

|𝐓|=|𝐫×𝐫||𝐫|2

Since

κ=|𝐓||𝐫|

, we can substitute

|𝐓|

with

|𝐫×𝐫||𝐫|2

and get:

κ=|𝐫×𝐫||𝐫|3

Here is a little summary on ways to calculate the curvature.

Definition With parametrization with respect to 𝐭 In terms of 𝐫(t) and its derivatives
κ=|d𝐓ds| κ=|𝐓||𝐫| κ=|𝐫×𝐫||𝐫|3

Motion in space

Velocity and acceleration

Remember in 2-dimensional calculus, we mentioned that a particle with displacement function

f(x)

has velocity

f(x)

and acceleration

f(x)

. In vector functions, the definition is basically the same. Suppose a particle moves through space so that its position vector at time

t

is

𝐫(t)

, its velocity function and acceleration function are:

𝐯(t)=limh0𝐫(t+h)𝐫(t)h=𝐫(t)

and

𝐚(t)=limh0𝐯(t+h)𝐯(t)h=𝐯(t)=𝐫(t)

Simply put:

𝐫(t)=𝐯(t)=𝐚(t)

.


The speed of the particle ignores the direction. It is the magnitude of the velocity vector: |𝐯(t)|. The distance traveled by the particle from t=a to b is ab|𝐯(t)|dt, which is also the formula for the arc length.


With the help of the Fundamental Theorem of Calculus, we can deduce the velocity function and position function when we know the particle's acceleration.

𝐯(t)=𝐯(t0)+t0t𝐚(u)du

and

𝐫(t)=𝐫(t0)+t0t𝐯(u)du

Tangential and normal acceleration

We can split the acceleration vector into two components: the tangential acceleration

aT

and the normal acceleration

aN

. The tangential acceleration faces the same direction as the unit tangent vector (

𝐓

), and the normal acceleration faces the same direction as the unit normal vector (

𝐍

). Since both

𝐓

and

𝐍

are unit vectors, the acceleration vector can be written as the sum of two vectors:

𝐚=aT𝐓+aN𝐍

Our goal is to figure out how to describe the two components.

Recall that

𝐓=𝐯|𝐯|

, thus

𝐯=|𝐯|𝐓

Now we differentiate both sides of the equation,

𝐯=(|𝐯|𝐓)=|𝐯|𝐓+|𝐯|𝐓

So we get

𝐚=|𝐯|𝐓+|𝐯|𝐓

.

Recall that

κ=|𝐓||𝐫|=|𝐓||𝐯|

, thus

|𝐓|=κ|𝐯|

.

Recall that

𝐍=𝐓|𝐓|

, thus

𝐓=|𝐓|𝐍=κ|𝐯|𝐍

.

We substitute

𝐓

for

κ|𝐯|𝐍

to yield:

𝐚=|𝐯|𝐓+κ|𝐯|2𝐍

That leaves us with:

aT=|𝐯|

and

aN=κ|𝐯|2

Of course, it will be more convenient if those components can be written in terms if

𝐫(t)

and its derivatives. Suppose that

θ

is the angle between

𝐓

and

𝐚

, then we can write

aT,aN

like this:

aT=|𝐚(t)|cosθ=|𝐯(t)||𝐚(t)|cosθ|𝐯(t)|algebraic manipulation=𝐯(t)𝐚(t)|𝐯(t)|dot product=𝐫(t)𝐫(t)|𝐫(t)|in terms of 𝐫(t) and its derivativesaN=|𝐚(t)|sinθ=|𝐯(t)||𝐚(t)|sinθ|𝐯(t)|algebraic manipulation=|𝐯(t)×𝐚(t)||𝐯(t)|cross product=|𝐫(t)×𝐫(t)||𝐫(t)|in terms of 𝐫(t) and its derivatives

To sum it up:

aT=|𝐯|=𝐫(t)𝐫(t)|𝐫(t)|

and

aN=κ|𝐯|2=|𝐫(t)×𝐫(t)||𝐫(t)|

Expansion

We only discussed vector functions with three variables

𝐫(t)=f(t),g(t),h(t)

. How about expanding our understanding of vector functions into

n

variables? Suppose we have a curve with vector function:

𝐫(t)=(f1(t)f2(t)fn(t))

Limits

The limit of a vector function is defined as:

limtc𝐫(t)=(limtcf1(t)limtcf2(t)limtcfn(t))

Differentiation and integration

The derivative of a vector function is defined as:

𝐫(t)=limΔt0𝐫(t+Δt)𝐫(t)Δt

, thus

𝐫(t)=(f1(t)f2(t)fn(t))

All differentiation rules apply.


We can expand integration into:

𝐫(t)dt=(f1(t)dtf2(t)dtfn(t)dt)+𝐂

Then, the arc length will become:

L=ab|𝐫(t)|dt

or

L=ab[f1(t)]2+[f2(t)]2++[fn(t)]2dt

Template:BookCat