Changing Variables and Jacobians

Download as pdf or txt
Download as pdf or txt
You are on page 1of 17

Lecture 3

Changing variables & Jacobians

3.1 Introduction
The slope of the f = f (x, y ) surface depends on the direction in which one moves from
the point (x, y ), and in Lecture 1 we learnt that the partial derivative must be defined
in terms of the change along a particular direction or axis.
For a function f (x, y ) we decided that the obvious directions to choose are along the
x and y axes, keeping y and x fixed, respectively.
The question now is are these really the obvious directions?
For example, consider the function
2
+y 2 )
f (x, y ) = e −(x cos(4(x 2 + y 2 ))
shown in Figure 3.1. Does it really make sense to impose a square “x-constant, y -
constant” cake-rack mesh onto this function? Probably not! It would be in more
sympathy with the function to use a mesh with radial symmetry as in Figure 3.2.

0.5
f(x,y)

−0.5
5
4
0 2
0
−2
y −5 −4
x

Figure 3.1: The function exp −(x 2 + y 2 ) cos(4(x 2 + y 2 )).

1
3/2 LECTURE 3. CHANGING VARIABLES & JACOBIANS

φ
lines of
constant y

circles of radial lines of


lines of
constant r constant φ
constant x
Figure 3.2: A radial “dartboard” mesh reflects the symmetry rather better than the x, y constant cake
rack.

To achieve this we must make an appropriate transformation to a new set of variables.


This raises the questions of
• How to choose the new variables and thus describe the transformation.
• How to describe the function in the new variables.
• How to find the partial derivatives with respect to these new variables.

3.1.1 The general problem


The more general problem is illustrated in Figure 3.3. Subfigure (a) shows a surface
f (x, y ) with lines of constant x and y underneath. The partial derivatives wrt x and y
are found by slicing or moving along these directions.
Figure 3.3(b) shows the exactly the same surface with lines of constant u and v un-
derneath. The partial derivatives wrt u and v are found by slicing along these new
directions.

f (x , y ) F (u , v )

(a) (b)

Figure 3.3: (a) Surface f (x, y ) above a set of x-constant and y -constant lines. (b) Surface F (u, v )
above a set of u-constant and v -constant curves but plotted in the x, y frame. The surface has exactly
the same shape.
3.2. REWRITING THE FUNCTION AND FINDING THE DERIVATIVES 3/3

3.1.2 Choosing the new variables


There is no routine method for doing this, but usually the problem symmetry drops
very large hints as to the best transformation.
Later we will consider some of the most standard transformations, viz between Cartesian
and

• Plane Polar Coordinates (2D): for radial symmetry


• Spherical Polar Coordinates (3D): for spherical symmetry, and
• Cylindrical Polars Coordinates (3D): for cylindrical symmetry.

We also consider what would make a bad transformation.

3.2 Rewriting the function and finding the derivatives


First though we should consider the general problems of rewriting the function and
finding the partial derivatives with respect to arbitrary new variables (u, v ).
There are two main cases to consider:

1. The straightforward case. OLD variables in terms of NEW.


x = x(u, v ) and y = y (u, v ) . (3.1)

2. The more complicated case. NEW variables in terms of OLD.


u = u(x, y ) and v = v (x, y ) . (3.2)

3.2.1 CASES 1A and 1B: OLD variables in terms of NEW


We are transforming from (x, y ) to (u, v ) coordinates, and the transformation is given
as “old in terms of new” variables. That is:
x = x(u, v ) and y = y (u, v ) . (3.3)

Case 1A:
If we know what the function f (x, y ) is explicitly, we can replace x and y in the function
to get the new function as

f = f (x(u, v ), y (u, v )) = F (u, v ) . (3.4)

Then one can work out the partial derivatives with respect to u and v directly from the
new function F (u, v ).
3/4 LECTURE 3. CHANGING VARIABLES & JACOBIANS

♣ Example of Case 1A.


Let f (x, y ) = x/y , and let the transformation be x = u and y = u/v .
u ∂f ∂f
⇒f = F (u, v ) = =v and ⇒ = 0, =1. (3.5)
(u/v ) ∂u ∂v

Case 1B: In this case we don’t have an explicit form for f (x, y ), and, therefore, cannot
find an explicit form for function F (u, v ).
But we can still write expressions for the partial derivatives using the Chain Rule for
partials from Lecture 2:
Case 1B: Use the Chain Rule for partials to obtain the derivatives
       
∂f ∂f ∂x ∂f ∂y
= + (3.6)
∂u ∂x ∂u ∂y ∂u
       
∂f ∂f ∂x ∂f ∂y
= + (3.7)
∂v ∂x ∂v ∂y ∂v
♣ Example of Case 1B.
We are told f = f (x, y ), and x = u and y = u/v .
Given this transformation, the best we can do is to write for any function f (x, y )
          
∂f ∂f ∂f 1 ∂f ∂f 1
= (1) + = + (3.8)
∂u ∂x ∂y v ∂x ∂y v
       
∂f ∂f ∂f u ∂f u
= (0) + − 2 =− . (3.9)
∂v ∂x ∂y v ∂y v 2

To check this result, let’s use (as in the example of Case 1A) f (x, y ) = x/y . Then
           2    
∂f 1 x 1 v u v v v
= (1) + − 2 = − = − = 0 (3.10)
∂u y y v u v u2 u u
x  u  u 2v
   
∂f 1
= (0) + − 2 − 2 = 2 =1 (3.11)
∂v y y v u v

This is exactly as found earlier at Equation 3.5.

3.2.2 A source of confusion introduced, and cleared up


Before continuing, it is worth checking that you are happy that the functions f () and
F () have different names. They work differently on their parameters, so should be
named differently.
3.2. REWRITING THE FUNCTION AND FINDING THE DERIVATIVES 3/5

But notice that the chain rule for partials seems to deliver an expression for ∂f /∂u and
∂f /∂v .
But you might complain that f is a function of (x, y ) NOT (u, v ), and so we should
have written ∂F/∂u and ∂F/∂v . You are correct. So how are we getting away with it?
We are using f as a “go-between” value f = f (x, y ) and f = F (u, v ). This is correct
because at corresponding values of (x, y ) and (u, v ) the functions do return exactly
the same value. (Look back at Figure 3.3.)
Some books introduce another symbol (z, say) as the “go-between”. They write z =
f (x, y ) and z = F (u, v ) and then write
∂z ∂z ∂x ∂z ∂y ∂z ∂z ∂x ∂z ∂y
= + , = + , (3.12)
∂u ∂x ∂u ∂y ∂u ∂v ∂x ∂v ∂y ∂v
and so on — but this seems over-fussy.

3.2.3 CASES 2A and 2B: NEW variables in terms of OLD


We still wish to transform the function f = f (x, y ) to (u, v ) coordinates. But now the
transformation is given as “new in terms of old” variables, that is:

u = u(x, y ) and v = v (x, y ) . (3.13)

One’s first thought is that this should be the better way of defining a transformation,
but it turns out to be really inconvenient!
New-ito-Old prevents us from finding either F (u, v ) or using the Chain Rule for Partials.
We cannot immediately work out ∂x/∂u, etc.
We proceed in one of two ways.
Case 2A: Try to invert the transformation.
Use u = u(x, y ) and v = v (x, y ) as simultaneous equations from which to find x =
x(u, v ) and y = y (u, v ).
If you can invert, the problem becomes a Case 1 problem.
♣ Example of Case 2A
Suppose you are told that

u=x and v = x/y . (3.14)

These are new in terms of old. But it is obvious that

x = u, and ⇒v = u/y , ⇒ y = u/v . (3.15)

We now have Old-in-terms-of-New, and can move to Case 1 above.


3/6 LECTURE 3. CHANGING VARIABLES & JACOBIANS

Case 2B: When the transformation is not invertible Write down the Chain Rule for
Partials the “wrong” way round! That is, we want ∂f /∂u, so start with ∂f /∂x.
   
∂f ∂f ∂u ∂f ∂v
= + (3.16)
∂x ∂u ∂x ∂v ∂x
   
∂f ∂f ∂u ∂f ∂v
= + (3.17)
∂y ∂u ∂y ∂v ∂y
Treat these as simultaneous equations in the unknown ∂f /∂u and ∂f /∂v and re-
arrange. Ie:
   
∂f ∂f ∂v ∂f ∂v ∂u ∂v ∂u ∂v
= − / − (3.18)
∂u ∂x ∂y ∂y ∂x ∂x ∂y ∂y ∂x
   
∂f ∂f ∂u ∂f ∂u ∂u ∂v ∂u ∂v
= − / − (3.19)
∂v ∂y ∂x ∂x ∂y ∂x ∂y ∂y ∂x
So we can recover the partial derivatives with respect to the new variables without an
explicit expression for the function in the new variables.

3.2.4 SUMMARY
To summarize:
1: If you have the transformation as OLD variables in terms of NEW
either
1A: find the function explicity and find the partials from it directly, or
1B: find just the partials using the Chain-rule-for-partials directly.
2: If you have the transformation as NEW variables in terms of OLD
either
2A: invert the given transformation to get OLD variables in terms of NEW, and
then go to CASE 1, or
2B: use the Chain-rule-for-partials for the OLD variables to give simultaneous
equations for the partial derivatives with respect to to the NEW variables.
3.2. REWRITING THE FUNCTION AND FINDING THE DERIVATIVES 3/7

♣ Example #1
Q: The elevation of a hill is described by
z = z(x, y ) = exp −(x 2 + y 2 ) .

(3.20)
By deriving the function in plane polars, find ∂z/∂r and ∂z/∂φ and comment on their
values. Determine the value of r where the hill is steepest.
A: Plane polars are (r, φ) where x = r cos φ, y = r sin φ.
This is old in terms of new, and we know the function, ⇒CASE 1A.
The hill is therefore
z = exp −r 2 .

(3.21)
So
∂z ∂z
= −2r exp −r 2

and =0. (3.22)
∂r ∂φ
The change of height is all in the radial direction. ∂z/∂φ = 0 means that if you move
at constant r (that is, “round” the hill) you will not change height at all.
The gradient is a function of r alone, so we can find the total derivative
d d
Slope = (−2r exp −r 2 ) = (−2 + 4r 2 ) exp −r 2
 
(3.23)
dr dr
√ √
This is zero when r = 1/ 2 and Slope = − 2 exp {−1/2}.
♣ Example #2.
Here is one that can be solved using all approaches.
Q: Find the partial derivatives ∂f /∂u and ∂f /∂v when
f = f (x, y ) = x 2 − y 2 and u = (x + y ) v = (x − y ) . (3.24)

A: Transformation is new in terms of old. ⇒CASE 2.


Try CASE 2A: Can we invert the transformation? Yes! Adding then substracting
we find
1 1
x = (u + v ) y = (u − v ) . (3.25)
2 2

Now goto CASE 1. We know the function, ⇒CASE 1A.


1
f = F (u, v ) = x 2 − y 2 = (u + v )2 − (u − v )2 = uv

(3.26)
4
3/8 LECTURE 3. CHANGING VARIABLES & JACOBIANS

So the derivatives are


∂f ∂f
= v and =u (3.27)
∂u ∂v
That’s the problem solved. But now notice that if CASE *A works, we can always find
just the derivatives using CASE *B.
For example, having inverted the transformation to get
1 1
x = (u + v ) y = (u − v ) .
2 2
we could use CASE 1B to find just the derivatives:
∂f ∂f ∂x ∂f ∂y 1 1
= + = (2x). + (−2y ). = x − y = v (3.28)
∂u ∂x ∂u ∂y ∂u 2 2
∂f ∂f ∂x ∂f ∂y 1 −1
= + = (2x). + (−2y ). =x +y =u (3.29)
∂v ∂x ∂v ∂y ∂v 2 2

Now pretend that we did not (or could not) invert the transformation. We could still
find the derivative using Case 2B.
Write the chain rule for partials the other way around:
   
∂f ∂f ∂u ∂f ∂v
= + (3.30)
∂x ∂u ∂x ∂v ∂x
   
∂f ∂f ∂u ∂f ∂v
= + (3.31)
∂y ∂u ∂y ∂v ∂y

and for simplicity fill in the values you know


   
∂f ∂f
2x = (1) + (1) (3.32)
∂u ∂v
   
∂f ∂f
−2y = (1) + (−1) (3.33)
∂u ∂v

Now use as simultaneous equations for ∂f /∂u and ∂f /∂v :


   
∂f ∂f
2x − 2y = 2 ⇒ = (x − y ) = v , (3.34)
∂u ∂u
and subtracting gives, once again, the same results:
   
∂f ∂f
2x + 2y = 2 ⇒ = (x + y ) = u . (3.35)
∂v ∂v
3.2. REWRITING THE FUNCTION AND FINDING THE DERIVATIVES 3/9

♣ Example #3. Harder!


Q: A quantity V = V (x, y ) satisfies Laplace’s equation

∂ 2V ∂ 2V
+ =0. (3.36)
∂x 2 ∂y 2

Express this equation in plane polar coordinates (r, φ) where x = r cos φ, y = r sin φ.
A: We want first to write down expressions for ∂V /∂x and ∂V /∂y which will involve
∂r /∂x, ∂r /∂y , ∂φ/∂x, and ∂φ/∂y . This is the wrong way round for the transformation.
So

• Either invert the transformation (Case 2A), or

• Write the CRfP the “other way round” (Case 2B) and use sim eqs.

Both will work, but let’s use 2B ...

       
∂V ∂V ∂x ∂V ∂y ∂V ∂V
= + = cos φ + sin φ (3.37)
∂r ∂x ∂r ∂y ∂r ∂x ∂y
       
∂V ∂V ∂x ∂V ∂y ∂V ∂V
= + = (−r sin φ) + (r cos φ) (3.38)
∂φ ∂x ∂φ ∂y ∂φ ∂x ∂y

Sorting out the simultaneous equations (do check!) gives

∂V ∂V sin φ ∂V
= cos φ − (3.39)
∂x ∂r r ∂φ
∂V ∂V cos φ ∂V
= sin φ + (3.40)
∂y ∂r r ∂φ


We can drop the V and consider ∂x as an operator

∂ ∂ sin φ ∂
= cos φ − (3.41)
∂x ∂r r ∂φ
∂ ∂ cos φ ∂
= sin φ + (3.42)
∂y ∂r r ∂φ
3/10 LECTURE 3. CHANGING VARIABLES & JACOBIANS

So the operator

∂2
  
∂ sin φ ∂ ∂ sin φ ∂
= cos φ − cos φ − (3.43)
∂x 2 ∂r r ∂φ ∂r r ∂φ
2  2 
∂ 1 ∂ 1 ∂
= cos2 φ 2 − cos φsin φ − 2 + (3.44)
∂r r ∂φ r ∂φ∂r
∂2 sin φ ∂ 2
   
1 ∂ −sin φ cos φ ∂
− sin φ −sin φ + cos φ + − −
r ∂r ∂φ∂r r r ∂φ r ∂φ2

which can be tidied up a little. Note how the operators move through to the right,
operating as they go.
Let’s follow through a couple of terms. The following is straightforward:
∂2
  
∂ ∂ 2
cos φ cos φ = cos φ 2 (3.45)
∂r ∂r ∂r
but this requires the product rule
       2
∂ sin φ ∂ 1 ∂ 1 ∂
cos φ − = −cos φ sin φ − 2 − cos φ sin φ (3.46)
∂r r ∂φ r ∂φ r ∂r ∂φ

The tidied expression for y is (you should check!)

∂2 2
1 ∂2
 
2 ∂ 2 ∂
= sin φ − cos φsin φ 2 − (3.47)
∂y 2 ∂r 2 r ∂φ r ∂φ∂r
∂2 2
 
1 ∂ 1 2 ∂
+ cos φ cos φ + sin φ + 2 cos φ (3.48)
r ∂r ∂φ∂r r ∂φ2

Adding them up:


∂2 ∂2
 2
1 ∂2

∂ 1∂
+ = + + 2 2 . (3.49)
∂x 2 ∂y 2 ∂r 2 r ∂r r ∂φ
So Laplace’s equation in plane polar coordinates is
∂ 2V 1 ∂V 1 ∂ 2V
+ + 2 2 = 0. (3.50)
∂r 2 r ∂r r ∂φ
3.3. JACOBIANS (JACOBIAN DETERMINANTS) 3/11

3.3 Jacobians (Jacobian determinants)


Recall the equations that appeared from solving the simultaneous equations in case 2B:
   
∂f ∂f ∂v ∂f ∂v ∂u ∂v ∂u ∂v
= − / − (3.51)
∂u ∂x ∂y ∂y ∂x ∂x ∂y ∂y ∂x
   
∂f ∂f ∂u ∂f ∂u ∂u ∂v ∂u ∂v
= − / − (3.52)
∂v ∂y ∂x ∂x ∂y ∂x ∂y ∂y ∂x
Notice that the denominator is the same in both expressions. It has special significance,
and is called the Jacobian.
It is common to write the Jacobian as a determinant, but there is also another useful
notation.
  ∂u ∂v

∂u ∂v ∂v ∂u ∂x ∂x
∂(u, v )
J= − = ∂u = (3.53)

∂x ∂y ∂x ∂y ∂v ∂(x, y )
∂y ∂y

As we shall see later, the Jacobian is especially useful when changing variables from
(x, y ) to (u, v ) in multiple integrals.

3.3.1 Multiple definitions of the Jacobian determinant?


In some books one reads that
∂u ∂v

∂(u, v ) ∂x ∂x
= (3.54)
∂(x, y ) ∂u ∂v

∂y ∂y
while in others
∂u ∂u
∂(u, v ) ∂x ∂y
= . (3.55)
∂(x, y ) ∂v ∂v
∂x ∂y
This is puzzling, until you realize that they are identical! However, perhaps the first
way is to be preferred for reasons we see now.

3.3.2 Jacobians and their inverses


In examples given later you will see that

∂(x, y ) ∂(u, v )
=1 (3.56)
∂(u, v ) ∂(x, y )

We shall gain an intuitive understanding of why this must be so when we come to use
Jacobians in multiple integration. However it is quite easy to prove using the Jacobian
matrix. This is not on the syllabus.
3/12 LECTURE 3. CHANGING VARIABLES & JACOBIANS

The Chain rule for partials says


∂f ∂f ∂x ∂f ∂y ∂f ∂f ∂x ∂f ∂y
= + = + (3.57)
∂u ∂x ∂u ∂y ∂u ∂v ∂x ∂v ∂y ∂v
or, switching the order of the products,
∂f ∂x ∂f ∂y ∂f ∂f ∂x ∂f ∂y ∂f
= + = + (3.58)
∂u ∂u ∂x ∂u ∂y ∂v ∂v ∂x ∂v ∂y
We can treat these as operator equations: in steps ...
∂ ∂x ∂ ∂y ∂ ∂ ∂x ∂ ∂y ∂
f = f + f f = f + f (3.59)
∂u ∂u ∂x ∂u ∂y ∂v ∂v ∂x ∂v ∂y
∂ ∂x ∂ ∂y ∂ ∂ ∂x ∂ ∂y ∂
⇒ f ≡ f + f f ≡ f + f (3.60)
∂u ∂u ∂x ∂u ∂y ∂v ∂v ∂x ∂v ∂y
In matrix notation, this operator is
      
∂/∂u ∂x/∂u ∂y /∂u ∂/∂x ∂/∂x
= = J u,v
x,y (3.61)
∂/∂v ∂x/∂v ∂y /∂v ∂/∂y ∂/∂y
The matrix J u,v
x,y is a Jacobian matrix.

Now repeat with the Chain Rule around the other way. You will arrive at the following
      
∂/∂x ∂u/∂x ∂v /∂x ∂/∂u ∂/∂u
= = J u,v (3.62)
∂/∂y ∂u/∂y ∂v /∂y ∂/∂v x,y ∂/∂v
But it follows that
     
∂/∂u ∂/∂x ∂/∂u
= J u,v
x,y x,y J u,v
= J u,v (3.63)
∂/∂v ∂/∂y x,y ∂/∂v
so that
h i−1
J u,v
x,y J u,v = I
x,y
⇒ J u,v
x,y = J u,v
x,y
. (3.64)

But a result from linear algebra tells us that, if A and B are square, |AB| = |A||B|,
and hence |A||A−1 | = 1. The result for Jacobian determinants follows immediately.

3.4 Three standard transformations


There are certain transformations which occur very frequently, viz:
(1) Cartesian to plane polar coordinates;
(2) Cartesian to spherical polar coordinates; and
(3) Cartesian to cylindrical polar coordinates.
3.4. THREE STANDARD TRANSFORMATIONS 3/13

3.4.1 Cartesian to plane polars


This is useful for problems with radial symmetry in the plane.

y (x,y) (r, φ )
r
φ
x
Figure 3.4: Cartesian to plane polars. Many texts use θ rather than φ.

The transformation is
x = r cos φ, y = r sin φ (3.65)
∂x ∂x ∂y ∂y
⇒ = cos φ; = −r sin φ; = sin φ; = r cos φ (3.66)
∂r ∂φ ∂r ∂φ
Writing the CRfP in operator form gives
∂ ∂x ∂ ∂y ∂ ∂ ∂
= + = cos φ + sin φ (3.67)
∂r ∂r ∂x ∂r ∂y ∂x ∂y
∂ ∂x ∂ ∂y ∂ ∂ ∂
= + = −r sin φ + r cos φ . (3.68)
∂φ ∂φ ∂x ∂φ ∂y ∂x ∂y
Using these are simutaneous equations gives
∂ ∂ sin φ ∂ ∂ ∂ cos φ ∂
= cos φ − and = sin φ + . (3.69)
∂x ∂r r ∂φ ∂y ∂r r ∂φ
(You could also find these by inverting the transformation.)
The Jacobian
∂x ∂x
∂(x, y ) ∂r ∂φ
= ∂y ∂y = r cos2 φ + r sin2 φ = r . (3.70)
∂(r, φ)
∂r ∂φ
and
∂r ∂r
∂(r, φ) ∂x ∂y cos φ sin φ 1
= ∂φ ∂φ = (cos φ)( ) − (sin φ)(− )= . (3.71)
∂(x, y ) r r r
∂x ∂y
3/14 LECTURE 3. CHANGING VARIABLES & JACOBIANS

3.4.2 3D Cartesian to cylindrical polars


Cylindrical polars add a z-axis, perpendicular to the plane, which is identical to the
Cartesian z-axis.

z r (x,y,z) (r, φ ,z)

y
φ
x
Figure 3.5: Cartesian x, y , z to cylindrical polars r, φ, z

The transformation is
x = r cos φ; y = r sin φ; z = z. (3.72)
So
∂ ∂ ∂
= cos φ + sin φ (3.73)
∂r ∂x ∂y
∂ ∂ ∂ ∂ ∂
= −r sin φ + r cos φ and = (3.74)
∂φ ∂x ∂y ∂z ∂z
and using these as simultaneous equations (or inverting the transformation):
∂ ∂ sin φ ∂
= cos φ − (3.75)
∂x ∂r r ∂φ
∂ ∂ cos φ ∂ ∂ ∂
= sin φ + and = (3.76)
∂y ∂r r ∂φ ∂z ∂z

The Jacobian is now a 3 × 3 determinant



∂x ∂y ∂z cos φ sin φ 0


∂r ∂r ∂r
∂(x, y , z) ∂x
∂y ∂z

= ∂φ ∂φ = −r sin φ r cos φ 0 = r, (3.77)

∂φ
∂(r, φ, z) ∂x ∂y ∂z

∂z ∂z ∂z 0 0 1
and (do check!)

∂r ∂φ ∂z
∂(r, φ, z) ∂x ∂x ∂x 1
∂r ∂φ ∂z
= = . (3.78)

∂y ∂y ∂y
∂(x, y , z) ∂r ∂φ ∂z r
∂z ∂z ∂z

Most entries can be copied over from the plane polar case.
3.4. THREE STANDARD TRANSFORMATIONS 3/15

3.4.3 Cartesian to spherical polars

Polar angle (x,y,z) (r, θ, φ)


z
θ
r
y

x φ

Azimuthal angle

Figure 3.6: Cartesian to spherical polars

The transformation involves the radius r of a sphere r 2 = x 2 + y 2 + z 2 , the polar angle


0 ≤ θ ≤ π and the azimuthal angle 0 ≤ φ ≤ 2π
x = r sin θ cos φ; y = r sin θ sin φ; z = r cos θ. (3.79)
Hence
∂ ∂ ∂ ∂
= sin θ cos φ + sin θ sin φ + cos θ (3.80)
∂r ∂x ∂y ∂z
∂ ∂ ∂ ∂
= r cos θ cos φ + r cos θ sin φ − r sin θ (3.81)
∂θ ∂x ∂y ∂z
∂ ∂ ∂
= −r sin θ sin φ + r sin θ cos φ (3.82)
∂φ ∂x ∂y
and
∂ ∂ 1 ∂ sin φ ∂
= sin θ cos φ + cos θ cos φ − (3.83)
∂x ∂r r ∂θ r sin θ ∂φ
∂ ∂ 1 ∂ cos φ ∂
= sin θ sin φ + cos θ sin φ + (3.84)
∂y ∂r r ∂θ r sin θ ∂φ
∂ ∂ sin θ ∂
= cos θ − (3.85)
∂z ∂r r ∂θ
The Jacobian

∂x ∂x ∂x sin θ cos φ r cos θ cos φ −r sin θ sin φ
∂r ∂θ ∂φ

∂(x, y , z) ∂y ∂y ∂y

= r 2 sin θ,
= = sin θ cos φ r cos θ sin φ r sin θ cos φ

∂r ∂θ ∂φ
∂(r, θ, φ) ∂z ∂z ∂z

cos θ −r sin θ 0

∂r ∂θ ∂φ

3/16 LECTURE 3. CHANGING VARIABLES & JACOBIANS

(3.86)

and as an exercise you should show that



∂r ∂r ∂r
∂(r, θ, φ) ∂x ∂y ∂z 1
∂θ ∂θ ∂θ
= ∂x ∂y ∂z = 2 . (3.87)
∂(x, y , z) ∂φ ∂φ ∂φ r sin θ
∂x ∂y ∂z

3.5 A “good” transformation


Cylindrical and spherical polars are particularly useful when solving problems involving
cylindrical and spherical regions. But what in general makes a “good” transformation.
If you were asked what is the nicest region shape to fit into Cartesian coordinates Oxy
you would answer a rectangle.
Similarly, if you drew your Ouv coordinate system at right angles and asked the same
question, you would again answer a square.
So given some arbitrary region shape in the xy plane, the best sort of transformation
is one which maps it onto a rectangle in the uv plane, as illustrated in Figure 3.7(a).
A specific example is given in Figure 3.7(b) for plane polars. The quadrant’s extent is
awkard to describe in Cartesians, but easy in describe in plane polars.

y v

x u

y φ

π/2

x r
a a

Figure 3.7: Useful transformations tend to take an awkward region shape in Oxy into a rectangle in
Ouv .
3.5. A “GOOD” TRANSFORMATION 3/17

3.5.1 Will any transformation work? When is one bad?


Another question is “will any transformation work?” The answer is no. If we start with
two (or n) independent variables (x, y ), we must ensure that the new two (or n) are
also independent.
The test for independence is rather simple.
If u(x, y ) and v (x, y ) are functionally INDEPENDENT then
∂(u, v )
6= 0. (3.88)
∂(x, y )
Recall that the Jacobian is the denominator in a relationship relating partial derivatives
with respect to different sets of variables. Trouble is expected if this is zero.
Proof: For reading only

The proof is straightforward. If u(x, y ) and v (x, y ) are functionally dependent, then
u = u(v ), which in turn means that a function z = z(u, v ) = 0 exists. Differentiating
z with respect to x and y gives:
∂z ∂u ∂z ∂v
+ = 0 (3.89)
∂u ∂x ∂v ∂x
∂z ∂u ∂z ∂v
+ = 0 (3.90)
∂u ∂y ∂v ∂y
For consistency between these two equations we must have ux = βuy and vx = βvy ,
where β is some constant. Thus

∂(u, v ) ux vx
= = (ux vy − vx uy ) = (βuy vy − βvy uy ) = 0 . (3.91)
∂(x, y ) uy vy
♣ Example
Q: Test whether there is a functional dependence between the variables (u, v ) if
u = x 2 + y + 1 and v = x 4 + 2x 2 y + y 2 + x 2 − y . (3.92)
A:
∂(u, v )
= (ux vy − vx uy ) (3.93)
∂(x, y )
= (2x)(2x 2 + 2y − 1) − (4x 3 + 4xy + 2x)(1)
= (4x 3 + 4xy − 2x) − (4x 3 + 4xy + 2x)
= 0
So, yes the variables are not independent. With a little work you’ll spot that the
dependence is v = u 2 − 3u + 2.

You might also like