Probability Density Under Transformation
|
|
- Madison Daniels
- 5 years ago
- Views:
Transcription
1 Probability Density Under Transformation Pramook Khungurn September 5, 5 Introduction In creating an algorithm that samples points from some domain, a problem that always comes up is the following: Let A and B be sets, p A ( ) be a probability density on A, and f be a function from A to B. If one samples x from A according to p A, then what is the probability density of f(x)? This document discusses the solution to the above problem and its application to construction of sampling algorithm. One-Dimensional Case. The Main Theorem We first start with the simplest case where A and B are both subsets of the real line R. Let x A. The number p A (x) means that, in the infinitesimal interval [x, x + δx), there exists p A (x)δx amount of probability mass. Here, δx is a differential quantity such that (δx) =. Assume that f is continuous and infinitely differentiable. The function f sends the interval [x, x + δx) to the interval [f(x), f(x + δx)). By Taylor expansion, f(x + δx) = f(x) + f (x)δx + O((δx) ) = f(x) + f (x)δx. So, the resuling interval is [f(x), f(x) + f (x)δx), which as width f (x) δx. This means that the mass p A (x)δx gets distributed to an interval of width f (x)δx. As a result: Density at point f(x) = p A(x)δx f (x) δx = p A(x) f (x). This density is defined only when f (x), which means that f is one-to-one in a neighborhood of x. As such, we have the following theorem. Theorem. Let A and B be subsets of R, p A be a probability density on A, f : A B be continuous and differentiable and f (x) for all x A. The induced probability density p B ( ) arisen from the process of sampling x according to p A and then computing f(x) is given by: p B (f(x)) = p A(x) f (x).
2 . The Inversion Method The above theorem can be used to create sampling algorithm for any integrable density function on the real line from a uniformly random sample from the interval [, ). In this situation, A = [, ) and p A (x) = for all x A. The density p B ( ) is given to us. We want to find function f : A B such that, for any x A: Multiply both sides by f (x), we have: Let P B be the CDF of p B : We have that: p B (f(x)) = p A(x) f (x) = f (x). p B (f(x)) f (x) =. P B (y) = y p B (t) dt. {P B (f(x))} = p B (f(x))f (x) = p B (f(x)) f (x) given that f is an increasing function. Let us assume that f is increasing for now. We have that Integrating both sides from t = to t = x, we have: x {P B (f(x))} = {P B (f(t))} dt = x P B (f(x)) P B (f()) = x. With the assumption that f() should correspond to the lowest number in the set B, we can safely set P B (f()) =. So, P B (f(x)) = x f(x) = P B (x). The CDF is an increasing function, so is its inverse. Moreover, P B () maps to the lowest number in the set B. So, it is a valid choice for f. In other words, to generate a point on the real line with probability distribution p B, simply apply the inverse of the CDF to a point x picked uniformly randomly from the interval [, )..3 Sampling from the Exponential Distribution dt We present a simple application of the inversion method. The exponential distribution with parameter λ is defined on [, ) with The CDF is given by: So, p(x) = λe λx. P (x) = e λx. P (y) = ln( y). Hence, to sample x acoording to the exponential distribution, we simply set: x := ln( ξ) where ξ is a randomly and uniformly sampled from the interval [, ).
3 3 Multi-Dimensional Case 3. The Main Theorem Let A, B R n, and p A ( ) be a probability density on A. Let f be given by: f (x, x,..., x n ) f (x, x,..., x n ) f(x, x,..., x n ) =. f n (x, x,..., x n ) be a function from A to B. The induced probability distribution p B arisen from the process of sampling a point x = (x, x,..., x n ) according to p A can then computing f(x) can again be computed by finding the volume of the image of the interval This volume is given by: [x, x + δx ) [x, x + δx ) [x n, x n + δx n ). Df(x) δx δx... δx n where f / x f / x f / x n f / x f / x f / x n Df(x) = f n / x f n / x f n / x n where all the partial derivatives are evaluated at x. Thus, p B (f(x)) = p A(x) Df(x). Notice that Df(x) is the factor that shows up when we perform change of variables during an integration. In two-dimensional space, we may write: [ ] x(u, v) f(u, v) =. y(u, v) In this case: Thus, Df(u, v) = x/ y/ x/ v y/ v. p A (u, v) p B (x, y) = x/ x/ v. y/ y/ v 3. The Polar Coordinate Transform The polar coordinate transforms two numbers (r, φ) to a point (x, y) on the plane as follows: x = r cos φ y = r sin φ, 3
4 which gives: x r = cos φ x = r sin φ φ y r = sin θ y φ = r cos φ So, if we sample a polar coordinate (r, φ) with probability distribution p A, then the distribution p B of the point (x, y) is given by: p A (r, φ) p p B (x, y) = cos φ A (r, φ) r sin φ = r cos φ + r sin φ = p A(r, φ). r sin θ r cos φ 3.3 Sampling Uniformly from the Unit Disk The unit disk is given by the polar coordinates in the set [, ] [, π). How should we be sampling the polar coordinates so that the resulting point distribution is uniform on the disk? In our case, we have that p B (x, y) = /π. So, we want p A such that: π = p A(r, φ) r p A (r, φ) = r π. A common strategy is to sample r and φ independently so that p A (r, φ) = p r (r)p φ (φ). Moreover, we shall sample φ uniformly from the interval [, π) so that p φ (φ) = /(π). Thus, p r (r) = r. The above distribution can be sampled with the inversion method. The CDF is given by: The inverse CDF is then: P r (r) = r r dr = [r ] r = r. P r (t) = t. So, we can sample points uniformly from the unit disk by setting: r := ξ φ = πξ where ξ and ξ are two independent random samples chosen uniformly from the interval [, ). 3.4 Sampling Uniformly from a Triangle Suppose we have a triangle in a plane with point A = (x A, y A ), B = (x B, y B ), C = (x C, y C ). Let us assume further that (B A) (C A) is pointing in the positive z-direction so that: area(abc) = (B A) (C A) = x B x A x C x A y B y A y C y A. 4
5 We wish to find a transformation f that takes a point (u, v) uniformly and randomly picked from the rectangle [, ) so that the distribution of (x, y) = f(u, v) is uniform on the triangle ABC. In this setting, we have that p A (u, v) =, and p B (x, y) = /area(abc). In other words, area(abc) = Df(u, v) Df(u, v) = x B x A y B y A x C x A y C y A. One way to generate a point on a triangle is to generate barycentric coordinates (α, β, γ) such that α, β, γ and α + β + γ =. Then, we can get a point on the triangle by computing In other words, (x, y) = αa + βb + γc = ( β γ)a + βb + γc = A + (B A)β + (C A)γ. x = x A + (x B x A )β + (x C x A )γ y = y A + (y B y A )β + (y C y A )γ. Our task is to figure out what β and γ are as functions of u and v. We have that x = (x B x A ) β + (x C x A ) γ x v = (x B x A ) β v + (x C x A ) γ v y = (y B y A ) β + (y C y A ) γ y v = (y B y A ) β v + (y C y A ) γ v. So, the matrix Df(u, v) is given by: [ (xb x Df(u, v) = A ) β + (x C x A ) γ (x B x A ) β v + (x C x A ) γ ] v (y B y A ) β + (y C y A ) γ (y B y A ) β v + (y C y A ) γ v [ ] [ ] xb x = A x C x A β/ β/ v. y B y A y C y A γ/ γ/ v Thus, Df(u, v) = x B x A x C x A y B y A y C y A β/ β/ v γ/ γ/ v x B x A x C x A y B y A y C y A = x B x A x C x A y B y A y C y A β/ β/ v γ/ γ/ v = β/ β/ v γ/ γ/ v β γ v β γ v =. What should β and γ be as functions of u and v? We have the constraint that β + γ. This condition is satisfied if we let β = g(u)( v) γ = g(u)v 5
6 where g(u) is a function such that g(u). With this choice of β and γ, we have that = β γ v β γ v = [g (u)( v)]g(u) [ g(u)][g (u)v] = g(u)g (u). It remains to find the function g with makes the above equation holds: g dg du = g dg = du g dg = du g = u g = u. Hence, a uniform distribution of points on triangle ABC can be generated by computing: ( u( v) uv)a + u( v)b + uvc where (u, v) is randomly and uniformly sampled from the rectangle [, ). 4 Dealing with 3D Manifolds 4. The Main Theorem Suppose that we have a differentiable function f that maps a set A R to a surface B R 3. We shall write: f x (u, v) f(u, v) = f y (u, v). f z (u, v) Again, let p A be a probability distribution on A. Given point (u, v) A, consider the rectangle [u + δu) [v + δv), which has area δuδv. This rectangle has probability mass p A (u, v)δuδv in it. We have that: where (u, v) f(u, v) (u + δu, v) f(u + δu, v) = f(u, v) + f u (u, v)δu (u, v + δv) f(u, v + δv) = f(u, v) + f v (u, v)δv (u + δu, v + δv) f(u + δu, v + δv) = f(u, v) + f u (u, v)δu + f v (u, v)δv f u (u, v) = f v (u, v) = f x f y f z f x v f y v f z v (u, v) (u, v), and (u, v) (u, v) (u, v). (u, v) In other words, the rectangle gets mapped to a parallelogram with sides defined by the vector f u (u, v)δu and f v (u, v)δv. The area of this parallelogram is given by: f u (u, v)δu f u (u, v)δv = f u (u, v) f u (u, v) δuδv 6
7 (Since the notation is getting a little unwieldy, let us drop the (u, v) arguments from the function from now on.) To compute the cross product, we make use of the following identity: So, Define We have that: a b + (a b) = (a a)(b b). f u f v = (f u f u )(f v f v ) (f u f v ) f u f v = (f u f u )(f v f v ) (f u f v ). E(u, v) = f u (u, v) f u (u, v) F (u, v) = f u (u, v) f v (u, v) G(u, v) = f v (u, v) f v (u, v). area of parallelogram = f u f v = EG F. In differential geometry, E, F, and G are called the coefficients of the first fundamental form. As a result, we have that the induced probability distribution is given by: p B (f(u, v)) = p A(u, v)δuδv f u f v δuδv = p A(u, v) EG F. 4. The Spherical Coordinate Transform The spherical coordinate is the transformation from (θ, φ) (, π) [, π) to a point ω on a 3D sphere S with: sin θ cos φ ω = sin θ sin φ. cos θ We then have that: So, cos θ cos φ ω θ = cos θ sin φ, sin θ sin θ sin φ ω φ = sin θ cos φ. E = cos θ cos φ + cos θ sin φ + sin θ = cos θ + sin θ = F = cos θ cos φ sin θ sin φ + cos θ sin φ sin θ cos φ + = G = sin θ sin φ + sin θ cos φ = sin θ EG F = sin θ = sin θ. 7
8 The inducted probability distribution is given by: p B (ω(θ, φ)) = p A(θ, φ) sin θ. However, since θ (, θ), we have that sin θ >. So, we can write: 4.3 Uniformly Sampling a Sphere p B (ω(θ, φ)) = p A(θ, φ) sin θ. We will use the identity to construct a sampling algorithm to sample a point on the unit sphere uniformly. The idea is to pick a probability distribution p A on (θ, φ) (, π) [, π) such that the induced probability distribution p B is the constant distribution /(4π). In other words: In other words: 4π = p A(θ, φ) sin θ. p A (θ, φ) = sin θ 4π. A common strategy is to sample φ independenty from θ so that p A (θ, φ) = p θ (θ)p φ (φ). Moreover, let us sample φ uniformly from [, π) so that p φ (φ) = /(π). In other words, p θ (θ) π = sin θ 4π p θ (θ) = sin θ. We can sample p θ (θ) using the inversion method. The CDF of p Θ is given by: P θ (θ) = θ So, the inverse function is given by: sin θ dθ = [ cos θ ] θ = cos cos θ = cos θ. P θ (u) = cos ( u). In conclusion, we compute θ and φ as: θ := cos ( ξ ) φ := πξ where ξ, ξ are two independent random numbers sampled uniformly from the interval [, ). Notice, however, that if the end goal is to get a point ω, there is no need to compute θ because θ never appears directly in the expression for ω. More specifically, sin θ cos φ ( ξ ) cos φ ω = sin θ sin φ = ( ξ ) sin φ. cos θ ξ 8
9 4.4 Sampling a Cosine-Weighted Hemisphere In this section, we want to sample the z-positive unit hemisphere such that the probability density being proportional to cos θ at each point. In this case: cos θ C = p A(θ, φ) sin θ C cos θ sin θ = p A(θ, φ), where C is the constant such that cos θ C is a probability distribution on the sphere. Again, we sample θ and φ independently with φ being uniform in [, π). So, The CDF then is given by: P θ (θ) = π C θ To determine C, note that P θ (π/) =, so π C cos θ sin θ = p θ(θ). cos θ sin θ dθ = π C [ cos θ = π C ( cos π ) = π C. In other words, C = π, and P θ (θ) = cos θ. Hence, we can sample the cosine-weighted hemisphere by setting: 4.5 From Area to Solid Angle cos θ := ξ φ := πξ. ] θ = π C ( cos θ). When shading from an area light source, a way to sample the incoming light direction is to sample a point on the light source s surface with some probability density p A and then convert the vector from the shaded point to the sampled point to a unit vector ω. In this section, we find the relation between p A and the induced probability density. For simplicity, let us say that the shaded point is at the origin and lying on the xy-plane so that the normal is the z-axis. Let r = (r x, r y, r z ) denote the point on the light source. Let n be the normal at r, and let s and t be the basis of the tangent plane at r in such a way that (s, t, n) is an orthonormal basis. The tangent plane is the set {r + us + vt (u, v) R }. The function f that maps the tangent plane to the direction is given by: ω = f(u, v) = Hence, using Lemma (proven in the appendix), we have: r + us + vt r + us + vt s f u (u, v) = r + us + vt r + us + vt (r s + u) r + us + vt 3 t f v (u, v) = r + us + vt r + us + vt (r t + v) r + us + vt 3 9
10 At (u, v) = (, ), we have that So, Next, f u (, ) = f v (, ) = s r r r 3 (r s) = s r r(r s) r 3 t r r r 3 (r t) = t r r(r t) r 3 E = r 4 r (r s) + r (r s) r 6 = r 4 r (r s) r 6 = r (r s) r 4 F = r (r s)(r t) (r s)(r t) r 6 = r 4 G = r (r t) r 4 EG F = r 4 r (r s) r (r t) + (r s) (r t) r 8 = r 4 r (r s) r (r t) r 8 = [ ( ) ( ) ] r r r 4 r s r t = r 4 [ (ˆr s) (ˆr t) ] (r s) (r t) r 8 where ˆr is the unit vector in the direction of r. Because s, t, n forms an orthonormal basis and ˆr =, we have that = ˆr = (ˆr s) + (ˆr t) + (ˆr n). So, Thus, EG F = r 4 [ (ˆr s) (ˆr t) ] = (ˆr n) r 4 EG F = (ˆr n) ˆr n r 4 = r. In conclusion, p B (f(r)) = r ˆr n p A(r) = r cos θ p A(r) 4.6 The Hair Coordinate System Transform The hair coordinate system maps (θ, φ) (π/, π/) [, π) to a sphere with the following transformation function: sin θ ω = cos θ cos φ. cos θ sin φ
11 So, cos θ ω θ = sin θ cos φ sin θ sin φ ω φ = cos θ sin φ cos θ cos φ E = cos θ + sin θ cos φ + sin θ sin φ = F = sin θ cos θ cos φ sin φ sin θ cos θ cos φ sin φ = G = cos θ sin φ + cos θ cos φ = cos θ EG F = cos θ = cos θ. However, since θ ( π/, π/), we have that cos θ >, so EG F = cos θ. So, the probability density transformation formula is: 4.7 Sampling for Diffuse Hair p B (ω(θ, φ)) = p A(θ, φ) cos θ. In this section, we want to sample the sphere so that p B (ω) cos θ. Applying the main theorem in this section, we have: cos θ C = p A(θ, φ) cos θ p A (θ, φ) = cos θ C. Again, we sample φ uniformly from [, π), and then sample θ independently from φ. So, To find C, we note that P θ (π/) =, so p θ (θ) = π C cos θ θ P θ (θ) = π cos θ dθ C π/ = π [ θ + sin θ cos θ ] θ C = π [ ] θ θ + sin(θ ) C π/ = π ( θ + sin(θ) + π ). C = π C ( π + + π ) = π C π/
12 So, C = π, and P θ (θ) = π ( θ + sin(θ) + π ). The above function cannot be inverted symbolically. So, in Mitsuba s implementation, they solve for it using Brent s method. 5 Appendix Lemma. Proof. a a = a a a ( a 3 a a ) a a = ( a a a ) a a = ( a a a ) a a a = ( a a a ( a a a )) a a = ( a a a a ( a a )) a = a a a ( a 3 a a )
SOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More information3. On the grid below, sketch and label graphs of the following functions: y = sin x, y = cos x, and y = sin(x π/2). π/2 π 3π/2 2π 5π/2
AP Physics C Calculus C.1 Name Trigonometric Functions 1. Consider the right triangle to the right. In terms of a, b, and c, write the expressions for the following: c a sin θ = cos θ = tan θ =. Using
More information1 Integration in many variables.
MA2 athaye Notes on Integration. Integration in many variables.. Basic efinition. The integration in one variable was developed along these lines:. I f(x) dx, where I is any interval on the real line was
More informationFigure 25:Differentials of surface.
2.5. Change of variables and Jacobians In the previous example we saw that, once we have identified the type of coordinates which is best to use for solving a particular problem, the next step is to do
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
More informationUNIVERSITY OF HOUSTON HIGH SCHOOL MATHEMATICS CONTEST Spring 2018 Calculus Test
UNIVERSITY OF HOUSTON HIGH SCHOOL MATHEMATICS CONTEST Spring 2018 Calculus Test NAME: SCHOOL: 1. Let f be some function for which you know only that if 0 < x < 1, then f(x) 5 < 0.1. Which of the following
More informationCampus Academic Resource Program Chain Rule
This handout will: Provide a strategy to identify composite functions Provide a strategy to find chain rule by using a substitution method. Identifying Composite Functions This section will provide a strategy
More informationMATH1120 Calculus II Solution to Supplementary Exercises on Improper Integrals Alex Fok November 2, 2013
() Solution : MATH Calculus II Solution to Supplementary Eercises on Improper Integrals Ale Fok November, 3 b ( + )( + tan ) ( + )( + tan ) +tan b du u ln + tan b ( = ln + π ) (Let u = + tan. Then du =
More informationAPPENDIX 2.1 LINE AND SURFACE INTEGRALS
2 APPENDIX 2. LINE AND URFACE INTEGRAL Consider a path connecting points (a) and (b) as shown in Fig. A.2.. Assume that a vector field A(r) exists in the space in which the path is situated. Then the line
More information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More information1 The Differential Geometry of Surfaces
1 The Differential Geometry of Surfaces Three-dimensional objects are bounded by surfaces. This section reviews some of the basic definitions and concepts relating to the geometry of smooth surfaces. 1.1
More informationJim Lambers MAT 169 Fall Semester Practice Final Exam
Jim Lambers MAT 169 Fall Semester 2010-11 Practice Final Exam 1. A ship is moving northwest at a speed of 50 mi/h. A passenger is walking due southeast on the deck at 4 mi/h. Find the speed of the passenger
More informationMath 10C - Fall Final Exam
Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 8.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.02 Lecture 8. hange of variables.
More informationC4 "International A-level" (150 minute) papers: June 2014 and Specimen 1. C4 INTERNATIONAL A LEVEL PAPER JUNE 2014
C4 "International A-level" (150 minute) papers: June 2014 and Specimen 1. C4 INTERNATIONAL A LEVEL PAPER JUNE 2014 1. f(x) = 2x 3 + x 10 (a) Show that the equation f(x) = 0 has a root in the interval [1.5,
More informationNote: Each problem is worth 14 points except numbers 5 and 6 which are 15 points. = 3 2
Math Prelim II Solutions Spring Note: Each problem is worth points except numbers 5 and 6 which are 5 points. x. Compute x da where is the region in the second quadrant between the + y circles x + y and
More informationMath 147, Homework 1 Solutions Due: April 10, 2012
1. For what values of a is the set: Math 147, Homework 1 Solutions Due: April 10, 2012 M a = { (x, y, z) : x 2 + y 2 z 2 = a } a smooth manifold? Give explicit parametrizations for open sets covering M
More informationChapter 3. Differentiable Mappings. 1. Differentiable Mappings
Chapter 3 Differentiable Mappings 1 Differentiable Mappings Let V and W be two linear spaces over IR A mapping L from V to W is called a linear mapping if L(u + v) = Lu + Lv for all u, v V and L(λv) =
More informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
More informationSolutions to old Exam 3 problems
Solutions to old Exam 3 problems Hi students! I am putting this version of my review for the Final exam review here on the web site, place and time to be announced. Enjoy!! Best, Bill Meeks PS. There are
More informationCS Tutorial 5 - Differential Geometry I - Surfaces
236861 Numerical Geometry of Images Tutorial 5 Differential Geometry II Surfaces c 2009 Parameterized surfaces A parameterized surface X : U R 2 R 3 a differentiable map 1 X from an open set U R 2 to R
More informationf(p i )Area(T i ) F ( r(u, w) ) (r u r w ) da
MAH 55 Flux integrals Fall 16 1. Review 1.1. Surface integrals. Let be a surface in R. Let f : R be a function defined on. efine f ds = f(p i Area( i lim mesh(p as a limit of Riemann sums over sampled-partitions.
More informationMATH 12 CLASS 23 NOTES, NOV Contents 1. Change of variables: the Jacobian 1
MATH 12 CLASS 23 NOTES, NOV 11 211 Contents 1. Change of variables: the Jacobian 1 1. Change of variables: the Jacobian So far, we have seen three examples of situations where we change variables to help
More informationFigure 21:The polar and Cartesian coordinate systems.
Figure 21:The polar and Cartesian coordinate systems. Coordinate systems in R There are three standard coordinate systems which are used to describe points in -dimensional space. These coordinate systems
More information1. (a) (5 points) Find the unit tangent and unit normal vectors T and N to the curve. r (t) = 3 cos t, 0, 3 sin t, r ( 3π
1. a) 5 points) Find the unit tangent and unit normal vectors T and N to the curve at the point P 3, 3π, r t) 3 cos t, 4t, 3 sin t 3 ). b) 5 points) Find curvature of the curve at the point P. olution:
More informationCourse 2BA1: Hilary Term 2006 Section 7: Trigonometric and Exponential Functions
Course 2BA1: Hilary Term 2006 Section 7: Trigonometric and Exponential Functions David R. Wilkins Copyright c David R. Wilkins 2005 Contents 7 Trigonometric and Exponential Functions 1 7.1 Basic Trigonometric
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More information4. CONTINUOUS RANDOM VARIABLES
IA Probability Lent Term 4 CONTINUOUS RANDOM VARIABLES 4 Introduction Up to now we have restricted consideration to sample spaces Ω which are finite, or countable; we will now relax that assumption We
More informationChange of Variables, Parametrizations, Surface Integrals
Chapter 8 Change of Variables, Parametrizations, Surface Integrals 8. he transformation formula In evaluating any integral, if the integral depends on an auxiliary function of the variables involved, it
More informationMcGill University April Calculus 3. Tuesday April 29, 2014 Solutions
McGill University April 4 Faculty of Science Final Examination Calculus 3 Math Tuesday April 9, 4 Solutions Problem (6 points) Let r(t) = (t, cos t, sin t). i. Find the velocity r (t) and the acceleration
More informationHandout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration
1 Handout 6: Rotational motion and moment of inertia Angular velocity and angular acceleration In Figure 1, a particle b is rotating about an axis along a circular path with radius r. The radius sweeps
More informationHOMEWORK 2 - RIEMANNIAN GEOMETRY. 1. Problems In what follows (M, g) will always denote a Riemannian manifold with a Levi-Civita connection.
HOMEWORK 2 - RIEMANNIAN GEOMETRY ANDRÉ NEVES 1. Problems In what follows (M, g will always denote a Riemannian manifold with a Levi-Civita connection. 1 Let X, Y, Z be vector fields on M so that X(p Z(p
More informationMATH 32A: MIDTERM 1 REVIEW. 1. Vectors. v v = 1 22
MATH 3A: MIDTERM 1 REVIEW JOE HUGHES 1. Let v = 3,, 3. a. Find e v. Solution: v = 9 + 4 + 9 =, so 1. Vectors e v = 1 v v = 1 3,, 3 b. Find the vectors parallel to v which lie on the sphere of radius two
More informationSurfaces JWR. February 13, 2014
Surfaces JWR February 13, 214 These notes summarize the key points in the second chapter of Differential Geometry of Curves and Surfaces by Manfredo P. do Carmo. I wrote them to assure that the terminology
More informationWeek 7: Integration: Special Coordinates
Week 7: Integration: Special Coordinates Introduction Many problems naturally involve symmetry. One should exploit it where possible and this often means using coordinate systems other than Cartesian coordinates.
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. M. A. and M.Sc. Scholarship Test. September 25, Time Allowed: 150 Minutes Maximum Marks: 30
NATIONAL BOARD FOR HIGHER MATHEMATICS M. A. and M.Sc. Scholarship Test September 25, 2010 Time Allowed: 150 Minutes Maximum Marks: 30 Please read, carefully, the instructions on the following page 1 INSTRUCTIONS
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationDouble Integrals. Advanced Calculus. Lecture 2 Dr. Lahcen Laayouni. Department of Mathematics and Statistics McGill University.
Lecture Department of Mathematics and Statistics McGill University January 9, 7 Polar coordinates Change of variables formula Polar coordinates In polar coordinates, we have x = r cosθ, r = x + y y = r
More informationMath 425 Notes 9. We state a rough version of the fundamental theorem of calculus. Almost all calculations involving integrals lead back to this.
Multiple Integrals: efintion Math 425 Notes 9 We state a rough version of the fundamental theorem of calculus. Almost all calculations involving integrals lead back to this. efinition 1. Let f : R R and
More informationMultiple Choice. Compute the Jacobian, (u, v), of the coordinate transformation x = u2 v 4, y = uv. (a) 2u 2 + 4v 4 (b) xu yv (c) 3u 2 + 7v 6
.(5pts) y = uv. ompute the Jacobian, Multiple hoice (x, y) (u, v), of the coordinate transformation x = u v 4, (a) u + 4v 4 (b) xu yv (c) u + 7v 6 (d) u (e) u v uv 4 Solution. u v 4v u = u + 4v 4..(5pts)
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationFinal exam (practice 1) UCLA: Math 32B, Spring 2018
Instructor: Noah White Date: Final exam (practice 1) UCLA: Math 32B, Spring 218 This exam has 7 questions, for a total of 8 points. Please print your working and answers neatly. Write your solutions in
More informationMATH Final Review
MATH 1592 - Final Review 1 Chapter 7 1.1 Main Topics 1. Integration techniques: Fitting integrands to basic rules on page 485. Integration by parts, Theorem 7.1 on page 488. Guidelines for trigonometric
More informationRule ST1 (Symmetry). α β = β α for 1-forms α and β. Like the exterior product, the symmetric tensor product is also linear in each slot :
2. Metrics as Symmetric Tensors So far we have studied exterior products of 1-forms, which obey the rule called skew symmetry: α β = β α. There is another operation for forming something called the symmetric
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
More informationGood Problems. Math 641
Math 641 Good Problems Questions get two ratings: A number which is relevance to the course material, a measure of how much I expect you to be prepared to do such a problem on the exam. 3 means of course
More informationDr. Allen Back. Nov. 5, 2014
Dr. Allen Back Nov. 5, 2014 12 lectures, 4 recitations left including today. a Most of what remains is vector integration and the integral theorems. b We ll start 7.1, 7.2,4.2 on Friday. c If you are not
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationPart IA Vector Calculus
Part IA ector alculus Based on lectures by B. Allanach Notes taken by Dexter hua Lent 215 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures. They
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES 3. The Product and Quotient Rules In this section, we will learn about: Formulas that enable us to differentiate new functions formed from old functions by
More informationStudy guide for Exam 1. by William H. Meeks III October 26, 2012
Study guide for Exam 1. by William H. Meeks III October 2, 2012 1 Basics. First we cover the basic definitions and then we go over related problems. Note that the material for the actual midterm may include
More informationMath 396. Map making
Math 396. Map making 1. Motivation When a cartographer makes a map of part of the surface of the earth, this amounts to choosing a C isomorphism from a 2-submanifold with corners S in the sphere S 2 (r)
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationX. Linearization and Newton s Method
163 X. Linearization and Newton s Method ** linearization ** X, Y nls s, f : G X Y. Given y Y, find z G s.t. fz = y. Since there is no assumption about f being linear, we might as well assume that y =.
More informationPart IB Geometry. Theorems. Based on lectures by A. G. Kovalev Notes taken by Dexter Chua. Lent 2016
Part IB Geometry Theorems Based on lectures by A. G. Kovalev Notes taken by Dexter Chua Lent 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More informationInstructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.
Exam 3 Math 850-007 Fall 04 Odenthal Name: Instructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.. Evaluate the iterated integral
More informationHOMEWORK 3 - GEOMETRY OF CURVES AND SURFACES. where ν is the unit normal consistent with the orientation of α (right hand rule).
HOMEWORK 3 - GEOMETRY OF CURVES AND SURFACES ANDRÉ NEVES ) If α : I R 2 is a curve on the plane parametrized by arc length and θ is the angle that α makes with the x-axis, show that α t) = dθ dt ν, where
More informationIntroduction to Vector Calculus
Introduction to Vector Calculus Phillip Zenor Department of Mathematics Auburn University Edward E. Slaminka Department of Mathematics Auburn University Donald Thaxton Department of Physics Auburn University
More information1. For each function, find all of its critical points and then classify each point as a local extremum or saddle point.
Solutions Review for Exam # Math 6. For each function, find all of its critical points and then classify each point as a local extremum or saddle point. a fx, y x + 6xy + y Solution.The gradient of f is
More informationLAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More informationMA2AA1 (ODE s): The inverse and implicit function theorem
MA2AA1 (ODE s): The inverse and implicit function theorem Sebastian van Strien (Imperial College) February 3, 2013 Differential Equations MA2AA1 Sebastian van Strien (Imperial College) 0 Some of you did
More informationNote: Final Exam is at 10:45 on Tuesday, 5/3/11 (This is the Final Exam time reserved for our labs). From Practice Test I
MA Practice Final Answers in Red 4/8/ and 4/9/ Name Note: Final Exam is at :45 on Tuesday, 5// (This is the Final Exam time reserved for our labs). From Practice Test I Consider the integral 5 x dx. Sketch
More informationn=0 ( 1)n /(n + 1) converges, but not
Math 07H Topics for the third exam (and beyond) (Technically, everything covered on the first two exams plus...) Absolute convergence and alternating series A series a n converges absolutely if a n converges.
More informationCore Mathematics C12
Write your name here Surname Other names Core Mathematics C12 SWANASH A Practice Paper Time: 2 hours 30 minutes Paper - E Year: 2017-2018 The formulae that you may need to answer some questions are found
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 4 (Products and quotients) & (Logarithmic differentiation) A.J.Hobson
JUST THE MATHS UNIT NUMBER 104 DIFFERENTIATION 4 (Products and quotients) & (Logarithmic differentiation) by AJHobson 1041 Products 1042 Quotients 1043 Logarithmic differentiation 1044 Exercises 1045 Answers
More informationJim Lambers MAT 280 Fall Semester Practice Final Exam Solution
Jim Lambers MAT 8 Fall emester 6-7 Practice Final Exam olution. Use Lagrange multipliers to find the point on the circle x + 4 closest to the point (, 5). olution We have f(x, ) (x ) + ( 5), the square
More informationIntroduction to the Fourier transform. Computer Vision & Digital Image Processing. The Fourier transform (continued) The Fourier transform (continued)
Introduction to the Fourier transform Computer Vision & Digital Image Processing Fourier Transform Let f(x) be a continuous function of a real variable x The Fourier transform of f(x), denoted by I {f(x)}
More informationPartial Fractions. Combining fractions over a common denominator is a familiar operation from algebra: 2 x 3 + 3
Partial Fractions Combining fractions over a common denominator is a familiar operation from algebra: x 3 + 3 x + x + 3x 7 () x 3 3x + x 3 From the standpoint of integration, the left side of Equation
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
More informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationGradient, Divergence and Curl in Curvilinear Coordinates
Gradient, Divergence and Curl in Curvilinear Coordinates Although cartesian orthogonal coordinates are very intuitive and easy to use, it is often found more convenient to work with other coordinate systems.
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationProblem Set 4. f(a + h) = P k (h) + o( h k ). (3)
Analysis 2 Antti Knowles Problem Set 4 1. Let f C k+1 in a neighborhood of a R n. In class we saw that f can be expressed using its Taylor series as f(a + h) = P k (h) + R k (h) (1) where P k (h).= k r=0
More informationTopic 5.6: Surfaces and Surface Elements
Math 275 Notes Topic 5.6: Surfaces and Surface Elements Textbook Section: 16.6 From the Toolbox (what you need from previous classes): Using vector valued functions to parametrize curves. Derivatives of
More informationHaus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN:
MIT OpenCourseWare http://ocw.mit.edu Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, 989. ISBN: 978032490207. Please use the following citation
More informationMaterial for review. By Lei. May, 2011
Material for review. By Lei. May, 20 You shouldn t only use this to do the review. Read your book and do the example problems. Do the problems in Midterms and homework once again to have a review. Some
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More informationSolutions to Homework 11
Solutions to Homework 11 Read the statement of Proposition 5.4 of Chapter 3, Section 5. Write a summary of the proof. Comment on the following details: Does the proof work if g is piecewise C 1? Or did
More informationIntroduction to Differentials
Introduction to Differentials David G Radcliffe 13 March 2007 1 Increments Let y be a function of x, say y = f(x). The symbol x denotes a change or increment in the value of x. Note that a change in the
More informationNotes on Cartan s Method of Moving Frames
Math 553 σιι June 4, 996 Notes on Cartan s Method of Moving Frames Andrejs Treibergs The method of moving frames is a very efficient way to carry out computations on surfaces Chern s Notes give an elementary
More information7. Let X be a (general, abstract) metric space which is sequentially compact. Prove X must be complete.
Math 411 problems The following are some practice problems for Math 411. Many are meant to challenge rather that be solved right away. Some could be discussed in class, and some are similar to hard exam
More informationChapter 4. The First Fundamental Form (Induced Metric)
Chapter 4. The First Fundamental Form (Induced Metric) We begin with some definitions from linear algebra. Def. Let V be a vector space (over IR). A bilinear form on V is a map of the form B : V V IR which
More informationCLASS X FORMULAE MATHS
Real numbers: Euclid s division lemma Given positive integers a and b, there exist whole numbers q and r satisfying a = bq + r, 0 r < b. Euclid s division algorithm: This is based on Euclid s division
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationMonte Carlo Techniques Basic Concepts
Monte Carlo Techniques Basic Concepts Image Synthesis Torsten Möller Reading Chapter 12, 13, 14.1 14.3 of Physically Based Rendering by Pharr&Humphreys Chapter 7 in Principles of Digital Image Synthesis,
More informationPractice Exam 1 Solutions
Practice Exam 1 Solutions 1a. Let S be the region bounded by y = x 3, y = 1, and x. Find the area of S. What is the volume of the solid obtained by rotating S about the line y = 1? Area A = Volume 1 1
More information(b) g(x) = 4 + 6(x 3) (x 3) 2 (= x x 2 ) M1A1 Note: Accept any alternative form that is correct. Award M1A0 for a substitution of (x + 3).
Paper. Answers. (a) METHOD f (x) q x f () q 6 q 6 f() p + 8 9 5 p METHOD f(x) (x ) + 5 x + 6x q 6, p (b) g(x) + 6(x ) (x ) ( + x x ) Note: Accept any alternative form that is correct. Award A for a substitution
More informationChapter 8: Continuous Probability Distributions
Chapter 8: Continuous Probability Distributions 8.1 Introduction This chapter continued our discussion of probability distributions. It began by describing continuous probability distributions in general,
More informationCore Mathematics C34
Write your name here Surname Other names Pearson Edexcel International Advanced Level Centre Number Candidate Number Core Mathematics C34 Advanced Monday 16 June 2014 Morning Time: 2 hours 30 minutes You
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationCalculus III. Math 233 Spring Final exam May 3rd. Suggested solutions
alculus III Math 33 pring 7 Final exam May 3rd. uggested solutions This exam contains twenty problems numbered 1 through. All problems are multiple choice problems, and each counts 5% of your total score.
More informationTHEODORE VORONOV DIFFERENTIAL GEOMETRY. Spring 2009
[under construction] 8 Parallel transport 8.1 Equation of parallel transport Consider a vector bundle E B. We would like to compare vectors belonging to fibers over different points. Recall that this was
More informationVectors, dot product, and cross product
MTH 201 Multivariable calculus and differential equations Practice problems Vectors, dot product, and cross product 1. Find the component form and length of vector P Q with the following initial point
More informationNotes on multivariable calculus
Notes on multivariable calculus Jonathan Wise February 2, 2010 1 Review of trigonometry Trigonometry is essentially the study of the relationship between polar coordinates and Cartesian coordinates in
More informationFurther Mathematical Methods (Linear Algebra) 2002
Further Mathematical Methods (Linear Algebra) 22 Solutions For Problem Sheet 3 In this Problem Sheet, we looked at some problems on real inner product spaces. In particular, we saw that many different
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More information