San José State University |
---|
applet-magic.com Thayer Watkins Silicon Valley & Tornado Alley U.S.A. |
---|
The Change of Variables Theorem for Integrals |
---|
The case for double integrals will be considered first; i.e., those over two dimensional regions. However there is one preliminary that must be established before they are considered.
Let u_{x} and u_{y} be the unit vectors in the x and y directions, respectively. Let two vectors be given by
The area P of the parallelogram for which the vectors A and B are two adjacent sides is given by
where |A| and |B| are the lengths of the two vectors and |θ| is the positive angle between them. Suppose A is the base of their parallelogram and B is the slanted side.
The area P can also be expressed as the value of the determinant created from the components of the vectors; i.e.,
| a_{x} | b_{x} | | ||
P = | det | | | | |
| a_{y} | b_{y} | |
For more on the topic of determinants and area see determinants and area/volumes.
Let ∫∫_{R}f(x, y)dxdy be an integral over Euclidean 2-space. Let T be a transformation from (x, y) space to (u, v) space such that
and such that the inverse transformation
exists. This requirement limits the transformations to continuous ones. Let the inverse transformation be denoted as W.
Then
∫∫_{R}f(x, y)dxdy = ∫∫_{S}|det(J)|f(x(u,v), y(u, v))dudv
where |det(J)| is the absolute value of the determinant of the derivative matrix J for the inverse transformation W. The derivative matrix for W is defined as
|(∂x/∂u) | (∂x/∂v)| | ||
J | = | | | | |
|(∂y/∂u) | (∂y/∂v)| |
This is known as the Jacobian matrix of the transformation. The Jacobian K for the transformation T is
|(∂u/∂x) | (∂u/∂y)| | ||
K | = | | | | |
|(∂v/∂x) | (∂v/∂y)| |
It is the inverse transformation that is relevant for the change of variables for integrals.
The relationship between J and K is that the composition of T and W or W and T is the identity transformation and hence the products KJ and JK are equal to the 2×2 identity matrix I. Thus
Furthermore, of course,
The region S is the image of the transformation T of the region R.
Now consider the four corners of an area element ΔxΔy in (x, y) space
Under the transformation T and for small Δx and Δy the four corners are transformed to
The two vectors emanating from (u_{0}, v_{0}) are:
When a determinant is formed from the components of these vectors its value can be reduced to
|(∂u/∂x) | (∂u/∂y)| | ||||
ΔuΔv | = | det | | | | | ΔxΔy |
|(∂v/∂x) | (∂v/∂y)| |
This is none other than
This means that
So, in changing the variables of integration in ∫∫_{R}f(x, y)dxdy from (x, y) to (u, v) through the transformation T what is needed is the Jacobian of the inverse transformation W.
Let (x, y) be the Cartesian coordinates of the plane and (r, θ) the polar coordinates. The transformation of (x, y) to (r, θ) would be given by
The inverse transformation, that is to say the transformation of (r, θ) to (x, y), is
It is the Jacobian of this inverse transformation that is relevant in changing the variables in a double integral from (x, y) to (r, θ). That Jacobian is r so
Consider an integral of the form
If the variables of integration are to be changed from (x, y, z) to (u, v, w) then the integral in terms of (u, v, w) that is equal to the above triple integral is
where J is the Jacobian derivative matrix of the transformation of (u,v,w) to (x,y,z); i.e.,
That Jacobian derivative matrix is
| (∂x/∂u) | (∂x/∂v) | (∂x/∂w)| | |
J = | |(∂y/∂u) | (∂y/∂v) | (∂y/∂w)| |
|(∂z/∂u) | (∂z/∂v) | (∂z/∂w)| |
The transformation from (x,y,z) to (u,v,w) is also relevant because the region of integration S in the (u,v,w) space is the image of R under that transformation.
Now consider the transformation of the volume element ΔuΔvΔw into the (x,y,z) space. There are eight corner points to the volume element, but it is only the following four that are relevant:
These can be viewed as the point (u_{0}, v_{0}, w_{0}) with three vectors (Δu, 0, 0), (0, Δv, 0), (0, 0, Δw).
The point (u_{0}, v_{0}, w_{0}) goes to the point (x_{0}, y_{0}, z_{0}). The vector (Δu, 0, 0) induces the vector (((∂x/∂u)Δu, ((∂y/∂u)Δu, ((∂z/∂u)Δu) of deviations from (x_{0}, y_{0}, z_{0}). Likewise the vector (0, Δv, 0) induces the vector ((∂x/∂v)Δv, ((∂y/∂v)Δv, ((∂z/∂v)Δv) of deviations. And finally the vector (0, 0, Δw) induces the vector (((∂x/∂w)Δw, ((∂y/∂w)Δw, ((∂z/∂u)Δw). These three vectors of deviations from the point (x_{0}, y_{0}, z_{0}) determine a parallelpiped volume element in the (x, y, z) space. The volume of the parallelpiped is equal to the determinant of the matrix formed using the three vectors of deviations as columns. That matrix is
| (∂x/∂u) | (∂x/∂v) | (∂x/∂w)| |
|(∂y/∂u) | (∂y/∂v) | (∂y/∂w)| |
|(∂z/∂u) | (∂z/∂v) | (∂z/∂w)| |
This matrix is simply the Jacobian of the transformation from (u,v,w) to (x,y,z), above denoted as J. Thus
To rigorously establish the equality it is necessary to go back to the limiting process which defines the integrals and show that the partitioning of region R by the images of partitioning elements of S is valid.
After going through the two and three dimensional cases it is clear how the general case of transforming the variables of integration from (x_{1}, …, x_{n}) to (u_{1}, …, u_{n}) should be carried out; i.e.,
where J is the Jacobian matrix of the inverse transformation of (u_{1}, …, u_{n}) to (x_{1}, …, x_{n})
The inverse transformation is represented as
When u_{k} changes by an amount Δu_{k} the x variables change by the amounts given by
The transformation of the parallelotopes (parallelogram, parallelopiped, etc.) of an area/volume element Δu_{1}*…*Δu_{n} to the (x_{1}, …, x_{n}) space creates an area/volume element whose generalized volume is
where the j,k element of the matrix M is
But this means that the matrix M is none other than the Jacobian matrix for the transformation of the u variables to the x variables. Therefore the area/volume element of dx_{l}*…*dx_{n} is transformed to
(To be continued.)
HOME PAGE OF applet-magic |
---|