球坐标系下的向量运算(英文)
Spherical Coordinates
Transforms
The forward and reverse coordinate transformations are
r =
x =r sin !cos "y =r sin !sin " z =r cos !
z $ !=arctan %
&=arctan (y , x )
where we formally take advantage of the two argument arctan
function to eliminate quadrant confusion.
Unit Vectors
The unit vectors in the spherical coordinate system are functions of position. It is convenient to express them in terms of the spherical coordinates and the unit vectors of the rectangular coordinate system which are not themselves functions of position.
!
ˆ ˆ ˆ r ˆ ==x x +y y +z z =x ˆ sin !cos "+y ˆ s in !sin "+z ˆ cos !r r r ˆ #r ˆ z ˆ ˆ s in "+y ˆ cos " " ==$x sin !ˆ ˆ ! =" #r ˆ =x ˆ cos !cos "+y ˆ cos !sin "$z ˆ sin !
Variations of unit vectors with the coordinates
Using the expressions obtained above it is easy to derive the following handy relationships: ˆ !r =0!r ˆ !r ˆ ˆ cos "cos #+y ˆ cos "sin #$ˆ =x z sin "=" !"ˆ !r ˆ ˆ sin "sin #+ˆ ˆ sin #+ˆ =$x y s in "cos #=($x y c os #)sin "=# sin "!#
ˆ !" =0!r ˆ !"
=0
!#ˆ !" ˆ ˆ cos "$y ˆ sin "=$r ˆ sin #+#=$x cos #!"
()
ˆ !" =0!r ˆ !"
ˆ s in "cos $#y ˆ sin "sin $#z ˆ cos "=#r ˆ
=#x !"ˆ !" ˆ ˆ cos "sin $+y ˆ cos "cos $=$=#x c os "
Path increment
!
d r expressed in spherical coordinates: We will have many uses for the path increment
$ !r ˆ ˆ ˆ ' ! !r !r
d r =d (r r ˆ )=r ˆ d r +rd r ˆ =r ˆ d r +r & dr +d "+d #)
!"!#( % !r
ˆ ˆ ˆ =r d r +" r d "+#r sin "d #
Time derivatives of the unit vectors
We will also have many uses for the time derivatives of the unit vectors expressed in spherical coordinates:
ˆ ˆ ˙ !r ˆ ˙ ˆ ˙ ˆ ˙ !r ˙ ˆ =!r r r ˙ +" +# =" " +## sin "!r !"!#ˆ ˆ ˆ ˙ !" !" ˙ !" ˙ ˆ ˙ ˆ ˙ ˙ ˆ " =r +" +# =$r " +# # cos " !r !"!#ˆ ˆ ˆ ˙ !# !# ˙ !# ˙ ˆ ˆ ˙ ˙ ˆ sin "+"# =r +" +# =$r cos "# !r !"!#
()
Velocity and Acceleration
The velocity and acceleration of a particle may be expressed in spherical coordinates by taking into account the associated
rates of change in the unit vectors:
! ! ˙ =r ˙ ˆ ˆ ˙ v =r r +r r ˙ ˆ ˙ +!v r r ! +"r " sin ! =r
! ! ˙ ˙ ˙ ˙ ˆ ˙ ˆ ˆ ˙ ˙ ˆ ˙ ˆ ˆ ˙ ˙ ˆ ˙ ˙ ˆ ˆ r a =v =r r ˙ +r ˙ ˙ +!r ! +!r ˙! ˙ +!r ! +"r " sin !+"r ˙" ˙ sin !+" r " sin !+"r " ! cos !
ˆ ˆ ˙ ˆ ˙ ˆ ˆ ˙ ˙ ˙ ˙ +r ˆ ˙˙ +#r ˆ ! ˙ =!! ˙ +" " sin !r r +"" ˙ cos !r ! +!r ˙! +!r !
ˆ ˙ ˙ ˆ ˆ ˙ ˙ ˆ ˙ +#r ˆ sin !+! cos !" r " sin !+"r ˙" ˙ sin !+"r " sin !+"r "! ˙ cos !
()()
˙ ˙ ˙ ˙ ˙ ˆ r ˙ ˙ !r "˙ ˙ a =r 2!r # 2sin "+" r " +2r " !r # 2sin "cos "+# r # sin "+2r "# ˙ cos "+2r # sin "[()]
The del operator from the definition of the gradient
Any (static) scalar field u may be considered to be a function of the spherical coordinates ! r , θ, and φ. The value of u
d r . That change may be determined changes by an infinitesimal amount du when the point of observation is changed by
from the partial derivatives as
du =
!u !u !u dr +d "+d #. !r !"!#
But we also define the gradient in such a way as to obtain the result ! ! du =!u "d r Therefore,
! ! !u !u !u dr +d "+d #=$u %d r !r !"!#
or, in spherical coordinates,
! ! ! !u !u !u
dr +d "+d #=$u dr +$u rd "+$ u r sin "d #
r "#!r !"!#
()()
=
()
and we demand that this hold for any choice of dr , d θ, and d φ. Thus,
from which we find
()
! !u
r
"u =, "r
()
! !u
#
1"u =, r "#
()
! !u
$
1"u
,
r sin #"$
Divergence
! !
! "A is carried out taking into account, once again, that the unit vectors themselves are functions of the The divergence
coordinates. Thus, we have
! ! & #$ˆ ˆ #% #) ˆ ˆ +ˆ +A $$+ +A %% ! "A =( + "A r r
#r r #$r sin $#%' *
where the derivatives must be taken before the dot product so that
()
! ! & #$ˆ ˆ #% #) !
! "A =( +++ "A #r r #$r sin $#%' * ! ! !
ˆ ˆ #A $ #A % #A
ˆ "=r +"+"#r r #$r sin $#%
ˆ ˆ & #A ˆ #A ˆ #A %ˆ #r #$ #% )
ˆ "( r ˆ +$=r +% +A r +A $+A %+
#r #r #r #r #r * ' #r
ˆ ˆ ˆ ˆ $ & #A ˆ #A ˆ #A %ˆ #r #$ #% ) +"( r r +$ +% +A r +A $+A % r ' #$#$#$#$#$#$* ˆ ˆ ˆ & #A r ˆ % #A ˆ #A %ˆ #r #$ #% )
ˆ ++"( r $ +% +A r +A $+A %
r sin $' #%#%#%#%#%#%*
With the help of the partial derivatives previously obtained, we find
! ! & #A ) #A ˆ #A %ˆ ! "A =r ˆ "( r ˆ +$ +% +0+0+0+
#r #r ' #r *
ˆ ) $ & #A ˆ #A ˆ #A %ˆ ˆ ˆ )+0+ +"( r r +$ +% +A r $ +A $(,r r ' #$#$#$* ˆ & #A ) % #A ˆ #A %ˆ ˆ ˆ ˆ ˆ +$ˆ sin $+$+"( r +% +A r sin $% +A $cos $% +A %,r cos $+ r sin $' #%#%#%*
[()]
& #A ) & 1#A A r ) & 1#A %A r A cos $) =( r +( + +( ++ ' #r * ' r #$r * ' r sin $#%r r sin $*
2A ) & 1#A A cos $) 1#A %& #A
=( ++ +( ++ +' #r r * ' r #$r sin $* r sin $#%
Curl
! !
! "A is also carried out taking into account that the unit vectors themselves are functions of the coordinates. The curl
Thus, we have
! ! & #$ˆ ˆ #% #) ˆ ˆ ˆ +A $$+ +A %% ! "A =( r ++ "A r r
#r r #$r sin $#%' *
where the derivatives must be taken before the cross product so that
()
! ! & #$ˆ ˆ #% #) !
+! "A =( ++ "A #r r #$r sin $#%' * ! ! !
ˆ ˆ
ˆ "#A +$ "#A +% "#A =r
#r r #$r sin $#%
& #A ˆ ˆ +=r "( r
#r ' ˆ $ & #A +"( r r ' #$+
ˆ ˆ ˆ #A ˆ #A %ˆ #r #$ #% )
$ +% +A r +A $+A %+ #r #r #r #r #r *
ˆ ˆ ˆ #A ˆ #A %ˆ #r #$ #% ) ˆ r +$ +% +A r +A $+A %
#$#$#$#$#$*
With the help of the partial derivatives previously obtained, we find
ˆ ˆ ˆ & #A ˆ % #A ˆ #A %ˆ #r #$ #% )
ˆ +$"( r r +% +A r +A $+A %
r sin $' #%#%#%#%#%#%*
! ! & #A ) #A ˆ #A %ˆ ˆ "( r ˆ +$! "A =r +% +0+0+0+
#r #r ' #r *
ˆ ) $ & #A #A ˆ #A %ˆ ˆ +"( r ˆ +$ +% +A r $ +A $(,r ˆ )+0+ r ' #$#$#$* ˆ & #A #A ˆ #A %ˆ ) % ˆ ˆ ˆ ˆ +$ˆ sin $+$+"( r r +% +A r sin $% +A $cos $% +A %,r cos $+ r sin $' #%#%#%*
[()]
& #A ˆ #A %=( % ,
#r #r ' ) & 1#A ˆ 1#A %) A ˆ ˆ $ + +( ,% +r ˆ +% +r #$r #$r * ' *
)
ˆ r + *
A %ˆ A %cos $& 1#A ˆ 1#A ˆ ,+( $ ,r $ +
r sin $#%r r sin $' r sin $#%
& 1#A %1#A A %cos $)
ˆ ( =r ,++
r sin $* ' r #$r sin $#%& #A 1#A A %) ˆ +$ ( ,%+, #r r sin $#%r ' *
1#A r A ) & #A ˆ +% ( ,++
' #r r #$r *
Laplacian
The Laplacian is a scalar operator that can be determined from its definition as
ˆ ˆ ! ! ˆ ˆ & #$) & #u $ #% # #u % #u ) +ˆ !u =! "!u =( r +"++ (
' #r r #$r sin $#%* ' #r r #$r sin $#%*
2
()
ˆ "=r
##r
ˆ ˆ & #u $ #u % #u )
ˆ r ++( ' #r r #$r sin $#%*
With the help of the partial derivatives previously obtained, we find
22* ˆ ˆ ˆ ˆ ' #2u % #u % #u & #u & #u
ˆ ") r $!u =r +$+,
r #%r #%#r r sin %#&r sin %#&#r + ( #r 2
ˆ ˆ ˆ $ #& #u $ #u % #u )
+"+( ++ r #$' #r r #$r sin $#%* ˆ ˆ ˆ % #& #u $ #u % #u ) +"r ++( + r sin $#%' #r r #$r sin $#%*
ˆ ˆ ˆ ˆ ˆ #u %% ' #u #2u r #2u & cos %#u & #2u *
+") % +$+$+,
r ( #r #r #%r #%r #%r sin %#&r sin %#&#%+
ˆ ˆ ˆ ˆ ˆ ' ˆ sin %+%& #u #2u & cos %#u % #2u r cos %#u & #2u * ˆ +") & sin %+r ++$+, r sin %( #r #r #&r #%r #%#&r sin %#&r sin %#&+ ' #2u * ' 1#u 1#2u * ' 1#u cos %#u 1#2u * =) +) +++, +) , ( #r + ( r #r r #%+ ( r #r r sin %#%r sin %#&+ ' #2u 2#u * ' 1#2u cos %#u * ' 1#2u * =) ++ +) , +) ,
r #r + ( r #%r sin %#%+ ( r sin %#&+ ( #r
1#' 2#u * 1#' #u * 1#u
=) r , +) sin %, +r #r ( #r + r sin %#%( #%+ r sin %#&Thus, the Laplacian operator can be written as
2