Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
0% found this document useful (0 votes)
8 views

Lecture 06 - Vector Functions

,,,,
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
8 views

Lecture 06 - Vector Functions

,,,,
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 13

APSC248 – Engineering Analysis III

Lecture 6 – Vector Functions

𝑃
𝑟⃑(𝑡)

𝑜 𝑦

𝐶
𝑥

Definition of vector function:

Given a vector 𝑟⃑(𝑡) where 𝑡 is a real variable. Suppose that vector 𝑟⃑(𝑡) is defined for each

value of 𝑡 over a given range of 𝑡, i.e., 𝑡 ∈ (𝑎, 𝑏). Then it is said that 𝑟⃑(𝑡) is a “vector

function” of 𝑡 over that interval (𝑎, 𝑏).

Note, 𝑟⃑(𝑡) might represent any physical nature such as position vector a particle or an electric

field vector. Vector 𝑟⃑ can also be a function of arc length 𝑠 along the curve, i.e., 𝑟⃑(𝑠). Or

vector 𝑟⃑ is a function of simply a pure mathematical parameter.

(1) Scalar Function and Vector Function

Scalar functions have forms like 𝑦 = 𝑓(𝑥) or 𝑦 = 𝑓(𝑥! , 𝑥" , … , 𝑥# ). For example,

- 𝑇(𝑥, 𝑦, 𝑧) can be a scalar function for temperature of a location in 3D space;

- 𝑃(𝑥, 𝑦, 𝑧) can be a scalar function for pressure of a location in 3D space;

In 3D space, scalar function tells us how to associate a number (e.g., temperature) with

each point in space.

@Dr. Yang Cao Page 1


APSC248 – Engineering Analysis III

Vector function is a mathematical function of one or more variables whose range is a

set of multidimensional vectors or infinite-dimensional vectors. Often the input of a

vector-valued function is a scalar, but in general the input can be a vector of either

complex or real variables. Vector functions (in 3d space), is a rule which tells us how

to associate a vector with each point (𝑥, 𝑦, 𝑧) in space. For example,

- Velocity of fluid, 𝑣⃑(𝑥, 𝑦, 𝑧). It specifies the speed of fluid and direction of flow

at a location (𝑥, 𝑦, 𝑧).

- Angular velocity of a rotating body, 𝜔


7⃑(𝑥, 𝑦, 𝑧). It represents angular velocity

of a point on the rotating body.

𝑣⃑(𝑥, 𝑦, 𝑧) means it is a vector function depends on three parameters (or variables). It

can be resolved into components as

𝑣⃑(𝑥, 𝑦, 𝑧) = 𝑣$ (𝑥, 𝑦, 𝑧)𝚤̂ + 𝑣% (𝑥, 𝑦, 𝑧)𝚥̂ + 𝑣& (𝑥, 𝑦, 𝑧)𝑘=

In many cases, parameters 𝑥, 𝑦, and 𝑧 are function of a common variable, such as time

𝑡 . This means vector function 𝑣⃑(𝑥, 𝑦, 𝑧) = 𝑣⃑>𝑥(𝑡), 𝑦(𝑡), 𝑧(𝑡)? = 𝑣⃑(𝑡) becomes a

function of parameter 𝑡. And similarly, it can be resolve into

𝑣⃑(𝑡) = 𝑣$ (𝑡)𝚤̂ + 𝑣% (𝑡)𝚥̂ + 𝑣& (𝑡)𝑘=

@Dr. Yang Cao Page 2


APSC248 – Engineering Analysis III

In summary, a scalar function is a function 𝑓(𝑥! , 𝑥" , … , 𝑥# ) of one or more variables

whose range is one-dimensional, as compared to a vector function, whose range is

three-dimensional (or, in general, 𝑛-dimensional).

We can relate scalar function and vector function to system concept. The system can

be represented by the function (e.g., 𝑓(∙), 𝑣⃑(∙) etc.). The system can have single input

or multiple inputs. The number of outputs of the system is referred to as “range” of

the function representing the system. If the system has only one output, the function

representing the system is a scalar function since its “range” is one-dimensional. If the

system has multiple outputs, the function representing the system is a vector function.

Example 1: Vector function 𝑟⃑(𝑡) = 〈2 cos 𝑡 , 4 sin 𝑡 , 𝑡〉

The spiral is the path traced by the


tip of the vector as 𝑡 increases from
𝑟⃑(𝑡) zero through 8π.

The vector 𝑟⃑(𝑡) has its tail at the


origin and its head at the
coordinates evaluated by the
functions
𝑥 = 2 cos 𝑡 , 𝑦 = 4 sin 𝑡 , 𝑧 = 𝑡.

@Dr. Yang Cao Page 3


APSC248 – Engineering Analysis III

(2) Derivative of Vector Function

Given a vector function 𝑟⃑(𝑡), the derivative of vector function 𝑟⃑(𝑡) is defined as

𝑑𝑟⃑(𝑡) 𝑟⃑(𝑡 + Δ𝑡) − 𝑟⃑(𝑡)


= 𝑟⃑ ' (𝑡) ≡ lim
𝑑𝑡 ()→+ Δ𝑡

If the limit on the right side exists for a given value of 𝑡 , we say that 𝑟⃑(𝑡) is

-.⃑())
differentiable at that 𝑡, and that 𝑟⃑ ' (𝑡) or -)
is its derivative.

Let 𝑟⃑(𝑡) = 𝑥(𝑡)𝚤̂ + 𝑦(𝑡)𝚥̂ + 𝑧(𝑡)𝑘=, then its derivative is still a vector.

𝑑𝑟⃑(𝑡) 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝑟⃑ ' (𝑡) = = 𝚤̂ + 𝚥̂ + 𝑘=
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡

Example 2: Position vector function 𝑟⃑(𝑡) = 3 cos 𝑡 𝚤̂ + 3 sin 𝑡 𝚥̂ + 𝑡 ! 𝑘9 (A


spiral)

The velocity vector is

𝑑𝑟⃑(𝑡)
𝑣⃑(𝑡) = = −3 sin 𝑡 𝚤̂ + 3 cos 𝑡 𝚥̂ + 2𝑡𝑘=
𝑑𝑡

The acceleration vector is

𝑑 " 𝑟⃑(𝑡)
𝑎⃑(𝑡) = = −3 cos 𝑡 𝚤̂ − 3 sin 𝑡 𝚥̂ + 2𝑘=
𝑑𝑡 "

The speed is

|𝑣⃑(𝑡)| = P9 cos " 𝑡 + 9 sin" 𝑡 + 4𝑡 " = P9 + 4𝑡 "

To find times when 𝑣⃑ and 𝑎⃑ are orthogonal

𝑣⃑ ∙ 𝑎⃑ = 9 sin 𝑡 cos 𝑡 − 9 cos 𝑡 sin 𝑡 + 4𝑡 = 4𝑡 = 0

Thus, the only time when 𝑣⃑ and 𝑎⃑ are orthogonal is when 𝑡 = 0.

@Dr. Yang Cao Page 4


APSC248 – Engineering Analysis III

(3) Rules of Vector Differentiation

𝑑
7⃑(𝑡)] = 𝑓 ' (𝑡)𝑢
[𝑓(𝑡)𝑢 7⃑' (𝑡)
7⃑(𝑡) + 𝑓(𝑡)𝑢
𝑑𝑡

𝑑 𝑑 𝑑
[𝑢
7⃑(𝑡) ± 𝑣⃑(𝑡)] = 𝑢 7⃑' (𝑡) ± 𝑣⃑ ' (𝑡)
7⃑(𝑡) ± 𝑣⃑(𝑡) = 𝑢
𝑑𝑡 𝑑𝑡 𝑑𝑡

𝑑
[𝑢 7⃑' (𝑡) ∙ 𝑣⃑(𝑡) + 𝑢
7⃑(𝑡) ∙ 𝑣⃑ (𝑡)] = 𝑢 7⃑(𝑡) ∙ 𝑣⃑ ' (𝑡)
Product rule _ 𝑑𝑡
𝑑
[𝑢7⃑(𝑡) × 𝑣⃑(𝑡)] = 𝑢 7⃑' (𝑡) × 𝑣⃑(𝑡) + 𝑢 7⃑(𝑡) × 𝑣⃑ ' (𝑡)
𝑑𝑡

𝑑 𝑑𝑢
7⃑ 𝑑𝑥
Chain rule: e𝑢
7⃑>𝑥(𝑡)?f =
𝑑𝑡 𝑑𝑥 𝑑𝑡

𝑑
[𝑢 7⃑' (𝑡) ∙ 𝑣⃑(𝑡) × 𝑤
77⃑(𝑡)] = 𝑢
7⃑(𝑡) ∙ 𝑣⃑(𝑡) × 𝑤 7⃑(𝑡) ∙ 𝑣⃑ ' (𝑡) × 𝑤
77⃑(𝑡) + 𝑢 77⃑(𝑡) + 𝑢 77⃑ ' (𝑡)
7⃑(𝑡) ∙ 𝑣⃑(𝑡) × 𝑤
𝑑𝑡

;⃑(𝑡) = 𝑡𝚤̂ + 𝑡 ! 𝚥̂ − 3𝑘9 ; 𝑣⃑ (𝑡) = 𝑡𝚤̂ + 𝑡𝚥̂ + 𝑡𝑘9


Example 3: 𝑢

𝑑
[𝑢 7⃑' (𝑡) ∙ 𝑣⃑(𝑡) + 𝑢
7⃑(𝑡) ∙ 𝑣⃑(𝑡)] = 𝑢 7⃑(𝑡) ∙ 𝑣⃑ ' (𝑡)
𝑑𝑡

i ? + e𝑡𝑖h + 𝑡2 =𝑗 − 3𝑘
= >𝑖h + 2𝑡𝑗=? ∙ >𝑡𝑖h + 𝑡𝑗= + 𝑡𝑘 i f ∙ >𝑖h + 𝑗= + 𝑘
i?

= (𝑡 + 2𝑡 " ) + (𝑡 + 𝑡 " − 3)

= 3𝑡 " + 2𝑡 − 3

@Dr. Yang Cao Page 5


APSC248 – Engineering Analysis III

Example 4: vector function of constant length, |𝑟⃑(𝑡)| = 𝑐

Position vector representing a particle moving on the surface of a sphere will always have

a constant length equal to the radius of the sphere.

𝑟⃑(𝑡) ∙ 𝑟⃑(𝑡) = |𝑟⃑(𝑡)|" = 𝑐 "

𝑑
[𝑟⃑(𝑡) ∙ 𝑟⃑(𝑡)] = 0 ⟹ 𝑟⃑ ' (𝑡) ∙ 𝑟⃑(𝑡) + 𝑟⃑(𝑡) ∙ 𝑟⃑ ' (𝑡) = 0
𝑑𝑡

𝑑𝑟⃑
2𝑟⃑ ' (𝑡) ∙ 𝑟⃑(𝑡) = 0 ⇒ 𝑟⃑ ∙ =0
𝑑𝑡

The vectors 𝑟⃑ ' (𝑡) and 𝑟⃑(𝑡) are orthogonal, because their dot product is 0. In summary, if

-.⃑
𝑟⃑(𝑡) is a differentiable vector function of 𝑡 of constant length, then 𝑟⃑ ∙ -) = 0. Similarly, if

-.⃑ -.⃑ -.⃑


𝑟⃑ is a function of the arc length, we have 𝑟⃑(𝑠) ∙ -2 = 0. -2
and -)
point to the same

#⃑
!" #⃑
!"
-.⃑ -.⃑
direction, but have different length. Note that -2 is a unit vector since -2 = !%
!& = !%
#⃑
!"
.
!%
3 3
!%

(4) Integrals of Vector Function

Let 𝑟⃑(𝑡) = 𝑥(𝑡)𝚤̂ + 𝑦(𝑡)𝚥̂ + 𝑧(𝑡)𝑘=, then

m 𝑟⃑(𝑡)𝑑𝑡 = m 𝑥(𝑡)𝑑𝑡 𝚤̂ + m 𝑦(𝑡)𝑑𝑡 𝚥̂ + m 𝑧(𝑡)𝑑𝑡 𝑘=

Or if 𝑅7⃑ is any antiderivative of 𝑟⃑(𝑡), then

m 𝑟⃑(𝑡)𝑑𝑡 = 𝑅7⃑ (𝑡) + 𝐶⃑

@Dr. Yang Cao Page 6


APSC248 – Engineering Analysis III

Example 5: 𝑟⃑(𝑡) = 2𝑡𝚤̂ − 𝚥̂ + 3𝑡 ! 𝑘9 , ∫ 𝑟⃑(𝑡) 𝑑𝑡 =?

m 𝑟⃑(𝑡)𝑑𝑡 = m>2𝑡𝚤̂ − 𝚥̂ + 3𝑡 " 𝑘= ?𝑑𝑡

= 𝑡 " 𝚤̂ − 𝑡𝚥̂ + 𝑡 4 𝑘= + 𝐶⃑

where 𝐶⃑ = 𝑐! 𝚤̂ + 𝑐" 𝚥̂ + 𝑐4 𝑘= is a constant vector.

(5) Tangent Vector, Unit Tangent Vector, and Unit Normal Vector

A tangent vector is a vector that follows the direction of a curve at a given point.

𝑧
7⃑(𝑡)
𝑇

𝑃
𝑟⃑(𝑡)

𝑜 𝑦

𝑥 𝐶

Suppose that curve 𝐶 is represented by vector function 𝑟⃑(𝑡) = 𝑥(𝑡)𝚤̂ + 𝑦(𝑡)𝚥̂ + 𝑧(𝑡)𝑘=.

Then the derivative of 𝑟⃑(𝑡) is the tangent vector of curve 𝐶 at point 𝑃 (in the direction of

increasing 𝑡). Note that if you think of 𝑟⃑(𝑡) as a position vector, derivative of 𝑟⃑(𝑡) gives the

velocity vector and velocity vector is always tangent to the trajectory.

7⃑(𝑡), then
Let tangent vector be denoted as 𝑇

𝑑𝑟⃑(𝑡) 𝑑𝑥 𝑑𝑦 𝑑𝑧
7⃑ (𝑡) =
𝑇 = 𝚤̂ + 𝚥̂ + 𝑘=
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡

@Dr. Yang Cao Page 7


APSC248 – Engineering Analysis III

Unit tangent vector is then

7⃑(𝑡)
𝑇 𝑟⃑ ' (𝑡)
𝑇=(𝑡) = = '
7⃑(𝑡)o |𝑟⃑ (𝑡)|
o𝑇

-.⃑
We can also show that 𝑇=(𝑡) = -2
. Proof is given below

𝑑𝑟⃑ 𝑑𝑟⃑ 𝑑𝑡 1
= = 𝑣⃑ = 𝑇=
𝑑𝑠 𝑑𝑡 𝑑𝑠 |𝑣⃑|

Principal unit normal vector

There is one vector of particular significant among the vectors orthogonal to the unit

tangent vector (as shown in the diagram below). This vector points in the direction in

which the curve is turning. Since 𝑇= has constant length, from Example 4, we have

𝑑𝑇=(𝑡)
⋅ 𝑇=(𝑡) = 0
𝑑𝑡
-56
7⃑(𝑡) =
We obtain a vector 𝑁 which orthogonal to 𝑇=(𝑡). Its unit vector is
-)

7⃑
𝑁 𝑑𝑇=⁄𝑑𝑡
i=
𝑁 =
7⃑ o o𝑑𝑇=⁄𝑑𝑡o
o𝑁

𝑁 is called principal unit normal vector for the smooth curve in plane. The following

diagram illustrates the relationship between the curve, the unit tangent vector and the unit

normal vector.

𝑑𝑇=⁄𝑑𝑡 𝑟⃑(𝑡)
i=
𝑁
o𝑑𝑇=⁄𝑑𝑡o

𝑇=
𝑇=

𝑑𝑇=⁄𝑑𝑡
i=
𝑁
o𝑑𝑇=⁄𝑑𝑡o

@Dr. Yang Cao Page 8


APSC248 – Engineering Analysis III

Binormal vector

The binormal vector is defined as

𝐵= = 𝑇= × 𝑁
i

𝐵= is a unit vector.

Normal Plane of curve 𝐶

The plane through a given point on curve 𝐶 and normal to vector 𝑇=, is called the normal

plane of curve 𝐶.

Remarks:
6
-5 6
-5⁄-)
• 𝑇= is an unit vector, but one can’t conclude that -) is a unit vector. Of course, |-56⁄-)|

guarantees to be a unit vector. You can use 𝑇= = cos 2𝑡 𝚤̂ + sin 2𝑡 𝚥̂ to verify.

• To relate to a physical meaning, we can also say tangent vector gives us the velocity

7⃑ = -.⃑()) and acceleration would be 𝑎⃑ = -9:⃑ = -5:⃑.


𝑣⃑ = 𝑇 -) -) -)

:⃑
-5 -56
7⃑ (𝑡) ≠
Note that normal vector 𝑁 7⃑(𝑡) =
, instead, 𝑁 .
-) -)

As pointed out in Textbook on page 735, we can decompose acceleration vector as

𝑎⃑ = 𝑎 5 𝑇= + 𝑎; 𝑁
i

@Dr. Yang Cao Page 9


APSC248 – Engineering Analysis III

Example 6: A roller coaster based on the parametric equations of a 3D


Lissajous curve
𝑟⃑(𝑡) = sinE𝑠(𝑡)G 𝚤̂ + sinE2𝑠(𝑡)G 𝚥̂ + cosE𝑠(𝑡)G 𝑘9
𝑠(𝑡) represents the distance along the track. Obtain the velocity of the
roller coaster.

Solution:

𝑥(𝑡) = sin>𝑠(𝑡)? ; 𝑦(𝑡) = sin>2𝑠(𝑡)? ; 𝑧(𝑡) = cos>𝑠(𝑡)?

𝑑𝑥(𝑡) 𝑑𝑠
= cos>𝑠(𝑡)? = cos>𝑠(𝑡)? 𝑠̇
𝑑𝑡 𝑑𝑡

𝑑𝑦(𝑡) 𝑑𝑠
= 2 cos>2𝑠(𝑡)? = 2 cos>2𝑠(𝑡)? 𝑠̇
𝑑𝑡 𝑑𝑡

𝑑𝑧(𝑡) 𝑑𝑠
= − sin>𝑠(𝑡)? = − sin>𝑠(𝑡)? 𝑠̇
𝑑𝑡 𝑑𝑡

The velocity is then

𝑑𝑟⃑(𝑡) i ?𝑠̇
𝑣⃑(𝑡) = = >cos>𝑠(𝑡)? 𝑖h + 2 cos>2𝑠(𝑡)? 𝑗= − sin>𝑠(𝑡)? 𝑘
𝑑𝑡

@Dr. Yang Cao Page 10


APSC248 – Engineering Analysis III

Example 7: 𝑟⃑(𝑡 ) = 𝑡 ! 𝚤̂ + 3𝑡𝚥̂ + 5𝑘9 , Obtain the unit tangent vector and
unit normal vector at 𝑡 = 0.

Solution:
𝑑𝑟⃑(𝑡)
7⃑(𝑡) =
𝑇 = 2𝑡𝚤̂ + 3𝚥̂
𝑑𝑡
7⃑(𝑡)
𝑇 1
𝑇=(𝑡) = = (2𝑡𝚤̂ + 3𝚥̂)
7⃑(𝑡)o √4𝑡 " + 9
o𝑇

𝑑𝑇=(𝑡) 2𝑡 '
3 '
=y z 𝚤̂ + y z 𝚥̂
𝑑𝑡 √4𝑡 " + 9 √4𝑡 " + 9

< ' <' 9=9 ' <


Recall the quotient rule, e9 f = 9(

1 !
2𝑡 ' 2√4𝑡 " + 9 − y2 (4𝑡 " + 9)=" 8𝑡z (2𝑡) 2(4𝑡 " + 9) − 8𝑡 " 18
y z = = 4 = 4
√4𝑡 " + 9 4𝑡 " +9 (4𝑡 " + (4𝑡 " + 9)"
9)"
'
3 𝑑 ! 1 4 −12𝑡
y z =3 y(4𝑡 " + 9 )=" z = 3 y− z (4𝑡 " + 9 )=" (8𝑡) = 4
√4𝑡 " + 9 𝑑𝑡 2
(4𝑡 " + 9)"

𝑑𝑇=(𝑡) 1 6
= 4
(18𝚤̂ − 12𝑡𝚥̂) = 4
(3𝚤̂ − 2𝑡𝚥̂)
𝑑𝑡 (4𝑡 + 9)
" " (4𝑡 + 9)
" "

𝑑𝑇=(𝑡) 6
} }= P
4 9 + 4𝑡
"
𝑑𝑡 (4𝑡 " + 9)"

𝑑𝑇=⁄𝑑𝑡 1
i=
𝑁 = (3𝚤̂ − 2𝑡𝚥̂)
o𝑑𝑇=⁄𝑑𝑡 o √9 + 4𝑡 "

At 𝑡 = 0, we have

1 1
𝑇=(𝑡) = (2𝑡𝚤̂ + 3𝚥̂)~ = (0𝚤̂ + 3𝚥̂) = 𝚥̂
√4𝑡 " + 9 )>+
3

1
i=
𝑁 (18𝚤̂ − 12𝑡𝚥̂)~ = 𝚤̂
√324 + 144𝑡 " )>+

@Dr. Yang Cao Page 11


APSC248 – Engineering Analysis III

(6) Curvature of a Curve (optional)

Consider the unit tangent vector 𝑇=(𝑡) of a curve at point 𝑃 which moves on a smooth

curve. Since 𝑇=(𝑡) is unit vector, its length remains constant and only its direction changes

as the point 𝑃 is moving along the curve. The rate at which 𝑇=(𝑡) turns per unit of length

along the curve is call the curvature.

Definition:

If 𝑇=(𝑡) is the unit vector of a smooth curve, the curvature function of the curve is

𝑑𝑇= 𝑑𝑇= 𝑑𝑡 1 𝑑𝑇= 1 𝑑𝑇=


𝜅=} }=} }= } }= } }
𝑑𝑠 𝑑𝑡 𝑑𝑠 |𝑑𝑠/𝑑𝑡| 𝑑𝑡 |𝑣⃑| 𝑑𝑡

𝜅 is a Greek letter pronounced “kappa”.

Comments:
6
-5
• If • -2 • is large, 𝑇= turns sharply around point 𝑃 location, and the curvature at 𝑃 is

large.

-56
• If • -2 • is close to zero, 𝑇= turns slowly and the curvature at 𝑃 is small.

Example: Curvature of a circle, 𝑟⃑(𝑡 ) = (𝑎 cos 𝑡 )𝚤̂ + (𝑎 sin 𝑡)𝚥̂

𝑑𝑟⃑
𝑣⃑ = = −(𝑎 sin 𝑡)𝚤̂ + (𝑎 cos 𝑡)𝚥̂
𝑑𝑡

𝑣⃑
𝑇= = = −(sin 𝑡)𝚤̂ + (cos 𝑡)𝚥̂
|𝑣⃑|

𝑑𝑇=
= −(cos 𝑡)𝚤̂ − (sin 𝑡)𝚥̂
𝑑𝑡

1 𝑑𝑇= 1 1
𝜅= } }= =
|𝑣⃑| 𝑑𝑡 𝑎 radius

@Dr. Yang Cao Page 12


APSC248 – Engineering Analysis III

i in terms of radius of curvature 𝜅.


We can also obtain the principal unit normal vector 𝑁

Consider the unit tangent vector 𝑇= as a function of arc length 𝑠. Since 𝑇=(𝑠) has constant

-56
length, then the derivative -2
is orthogonal to 𝑇= . The principal unit normal vector 𝑁
i is

-56
obtained by normalizing -2
, i.e.,

𝑑𝑇=⁄𝑑𝑠 1 𝑑𝑇=
i=
𝑁 =
o𝑑𝑇=⁄𝑑𝑠o 𝜅 𝑑𝑠

The vector 𝑑𝑇= ⁄𝑑𝑠 points in the direction in which 𝑇= turns as the curve bends. Therefore, if

we face in the direction of increasing arc length, the vector 𝑑𝑇=⁄𝑑𝑠 points toward the right if

𝑇= turns clockwise and toward the left if 𝑇= turns counterclockwise. In other words, the

i will point toward the concave side of the curve.


principal normal vector 𝑁

@Dr. Yang Cao Page 13

You might also like