Calculus of Variations and Computer Vision
|
|
- Laurence Riley
- 6 years ago
- Views:
Transcription
1 Calculus of Variations and Computer Vision Sharat Chandran Page 1 of 23 Computer Science & Engineering Department, Indian Institute of Technology, Bombay. sharat January 8, 2002
2 Page 2 of 23 The Top Ten Algorithms of the Century Metropolis Algorithm for Monte Carlo Simplex Method for Linear Programming Krylov Subspace Iteration Methods The Decompositional Approach to Matrix Computations QR Algorithm for Computing Eigenvalues The Fortran Optimizing Compiler Quicksort Algorithm for Sorting Fast Fourier Transform Integer Relation Detection Fast Multipole Method
3 Overview Page 3 of Sample problems from computer vision (fade out computer vision gurus) The problem of shape The problem of obtaining 3D information The problem of motion 2. The beauty of mathematics (fade out math gurus) The method of the calculus of variations Allied mathematical methods: Singular Value Decomposition, and the Method of Lagrange Multipliers 3. Sample solution to problems posed (or how computer scientists teach/cheat) 4. Concluding remarks
4 Shape Page 4 of 23 When placed out of context, images can be quite involved What you get is not what you see (WYGINWYS) What you want to compute: A contour (useful for obtaining properties such as size) Key mathematical concept: parametric curve x(s), y(s) Idealized image
5 Motion Page 5 of 23 Optical flow is an intensitybased approximation to image motion from sequential timeordered images Key mathematical concept: Two functions: u(x, y) and v(x, y) What does it look like? How do we compute optical flow?
6 A short video segment Specifications to help. Application of Optical Flow Page 6 of 23
7 Two Problems Page 7 of 23 Determine the equation of the curve joining two points (x 1, y 1 ) and (x 2, y 2 ) on the plane such that the length of the curve joining them is minimum Goal: Prove that your answer is correct (x 1, y 1 ) (x 2, y 2 ) Similar problem: Determine the equation of the support so that a ball placed at (x 1, y 1 ) reaches (due to gravity) (x 2, y 2 ) in the least possible time
8 Page 8 of 23 Functions and Functionals Differential segment length is ( dx 2 + dy 2 ) 1/2 where dx and dy are infinitesimal lengths in the x and y directions along the curve. Therefore the length of the curve joining the two end points (x 1, y 1 ) and (x 2, y 2 ) is J = = x2 x 1 x2 x 1 (1 + y ) dx F (x, y, y ) dx For the first problem, what should the function y be so that J is minimized? Compare: Find x 0 such that y = f(x) is minimum. J is a functional, a quantity that depends on functions rather than dependent variables.
9 Solution to the Shortest Length Problem The Euler-Lagrange equation is a necessary condition for minimising J F y d dx F y = 0 (1) Page 9 of 23 In our problem, F = 1 + y. Applying Equation 1 we get Solution: y is a constant. d dx ( y 1 + y ) = 0 The shortest distance curve is a straight line.
10 Detour: Why the Euler equation Page 10 of 23 J = x 2 x 1 F (x, f, f ) dx with boundary conditions f(x 1 ) = f 1 and f(x 2 ) = f 2. Let η(x) be a test function. Deform f by ɛη(x). Then, dj/ dɛ = 0 at the right place. This is true for all test functions η(x). The boundary conditions assert η(x 1 ) = η(x 2 ) = 0 If f(x) is replaced by f(x) + ɛη(x) then f (x) will be replaced by f (x) + ɛη (x). The integral then becomes J = x2 x 1 F (x, f + ɛη, f + ɛη ) dx
11 Page 11 of 23 Detour: Why the Euler equation If F is differentiable, expand the integrand in a Taylor series Differentiate with respect to ɛ and set to zero Apply integration by parts to get x2 x 1 η (x)f f dx = [η(x)f f] x 2 x 1 x2 The first term is zero due to the boundary conditions Hence x2 η(x)(f f df f x 1 dx )dx = 0 x 1 η(x) d(f f) dx dx,
12 Bottom Line Euler Equations Various generalizations are possible Higher order derivatives: J = x 2 x 1 F (x, f, f, f,...)dx Integrand may depend on several functions instead of only one Finding functions that have two independent variables J = F (x, y, f, f x, f y )dxdy Page 12 of 23 The bottom line Function to optimize F (x, ux )dx F (x, ux, u xx )dx F (x, ux, v x )dx F (x, y, ux, u y )dxdy The Euler-Lagrange equations F u d F dx u x = 0 F u d F dx u x d2 F dx 2 uxx = 0 F u d F dx u x = 0 F v d F dx v x = 0 F u d F dx u x d F dy u y = 0
13 The Snake Formulation Page 13 of 23 Start with any old curve (a snake ). Mathematically deform it to be the desired shape (based on image data). Curve is v(s) = (x(s), y(s)). Define energy of the curve to be E(v(s)) = 1 where E snake = E int + E ext 0 E snake (v(s)) ds (2) Find the curve that minimizes Equation 2 Remember, actual image is not so clean!
14 Page 14 of 23 More on Formulation The energy functional is given by E = 1 0 (E int(v(s)) + E img (v(s)) + E con (v(s))) ds E int = (α(s) v s (s) 2 + β(s) v ss (s) 2 )/2 E con = k(v(s) x) 2 E img = w line E line (v(s)) + w edge E edge (v(s)) + w term E term (v(s)) Line energy is E line (v(s)) = I(x, y) Edge energy is E edge (v(s)) = I(x, y) 2
15 Applying Euler-Lagrange to Snakes We minimize J 1 = (α(s)x 2 s + β(s)x 2 ss)/2 + E x (s)ds = X(s, x, x, x )ds Page 15 of 23 and J 2 = (α(s)y 2 s + β(s)y 2 ss)/2 + E y (s)ds = Y (s, y, y, y )ds For X(s, x, x, x ) the Euler equation is X x d ds X x + d2 ds 2X x = 0
16 More on Application Page 16 of 23 X x = E s, X x = α(s)x (s), and X x = β(s)x (s) d ds X x = α x + αx and d ds (X (s)) = β x + βx and d2 (X (s)) = β x + 2β x + βx ds 2 For purpose of illustration use only two terms: E s α x αx = 0 Numerically approximating we have E x E i+1 E i at location i α (s) α i+1 α i x (s) x i+1 x i α x = (α i+1 α i )(x i+1 x i ) = (α i+1 α i )x i+1 (α i+1 α i )x i x (s) x i+1 2x i + x i 1
17 Page 17 of 23 Reducing Snakes To A Matrix Equation E E E n (α 2 α 1 ) (α 2 α 1 ) (α 3 α 2 ) (α 3 α 2 ) (α n+1 α n ) T α x 1 α x = 0... α n x n This can be written as Ax = B, and solved using SVD. x 1 x 2. x n
18 Page 18 of 23 The Method of Lagrange Multipliers Minimize a function f(x, y) subject to g(x, y) = 0. For example, find a point on the line x sin θ y cos θ + p = 0 closest to origin That is, minimize x 2 + y 2 subject to the given constraint. In MOLM, we generate a new function E = f + λg and set the partial derivatives with respect to x, y and λ to zero. In the example The solution is the obvious one 2x + λ sin θ = 0 2y + λ cos θ = 0 x = p sin θ y = +p cos θ x sin θ y cos θ + p = 0 The method can be generalized to larger number of constraints.
19 Page 19 of 23 Formulation for Motion Let E(x, y, t) be the intensity value at point (x, y) in the image Due to motion, let this point correspond to a point (x + δx, y + δy) at time instance t + δt Key assumption: E(x, y, t) = E(x + δx, y + δy, t + δt) Two parts to determine u(x, y) and v(x, y) Use Taylor s theorem to get E x u + E y v + E t = 0 Introduce a smoothness constraint: f(u, v, u x, v x, u y, v y ) = u x2 + u y2 + v x2 + v y 2 Using MOLM we have a functional that needs to be minimized f(u, v, u x, v x, u y, v y ) + λg(u, v, t) 2 dxdy u x 2 + u y 2 + v x 2 + v y 2 + λ(e x u + E y v + E t ) 2 dxdy
20 Page 20 of 23 Applying Euler-Lagrange for Optical Flow We have F u x u x F y u y = 0 F v x v x F y v y = 0 Applying these we get F u = 2λE x (E x u + E y v + E t ) F ux = 2u x, x F u x = 2u xx F uy = 2u y, And F v = 2λE x (E x v + E y v + E t ) F vx = 2v x, x F v x = 2v xx F vy = 2v y, This implies a coupled system of equations 2 u = λ(e x u + E y v + E t )E x 2 v = λ(e x u + E y v + E t )E y y F u y y F v y = 2u yy = 2v yy
21 Page 21 of 23 Discretizing we have Solution using Euler Lagrange u k,l, u k+1,l u k,l+1 u k 1,l u k,l 1 : u component of the velocity at points (k, l), (k + 1, l), (k, l + 1), (k 1, l) and (k, l 1) respectively. v k,l, v k+1,l v k,l+1 v k 1,l v k,l 1 : v component of the velocity at points (k, l), (k + 1, l), (k, l + 1), (k 1, l) and (k, l 1) respectively. u k,l = u k,l λe x E t 1 + λe x 2 v k,le x E t 1 + λe x 2 v k,l = v k,l λe x E t 1 + λe x 2 u k,le x E t 1 + λe x 2
22 Uses of Euler Lagrange Page 22 of 23 P roblem Regularization principle Contours Esnake (v(s))ds Area based Optical flow 2 [(ux + u 2 y + v 2 x + v 2 y ) + λ(e xu + E yv + i t) 2 )]dxdy Edge detection [(Sf i) 2 + λ (f xx) 2 ]dx Contour based Optical flow [(V N V N ) 2 + λ( δv δx )2 ] Surface reconstruction [(S f d 2 + λ(f xx + 2f 2 xy + f 2 yy )]dxdy Spatiotemporal approximation [(S f i) 2 + λ( f V + ft) 2 ]dxdydt Colour I y Ax 2 + λ P z 2 Shape from shading [(E R(f, g)) 2 + λ(f 2 x + f 2 y + g 2 x + g 2 y )]dxdy Stereo {[ 2 G (L(x, y) R(x + d(x, y), y))] 2 + λ( d) 2 }dxdy
23 Page 23 of 23 Concluding Remarks Two interesting problems have been described Defined and derived the Euler Lagrange equations Used mathematics to solve the computer vision problem
Introduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 45 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics October 11, 25 Chapter
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
More informationTaylor Series and stationary points
Chapter 5 Taylor Series and stationary points 5.1 Taylor Series The surface z = f(x, y) and its derivatives can give a series approximation for f(x, y) about some point (x 0, y 0 ) as illustrated in Figure
More informationSuggested Solution to Assignment 7
MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationErrata for Robot Vision
Errata for Robot Vision This is a list of known nontrivial bugs in Robot Vision 1986) by B.K.P. Horn, MIT Press, Cambridge, MA ISBN 0-6-08159-8 and McGraw-Hill, New York, NY ISBN 0-07-030349-5. Thanks
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationErrata for Robot Vision
Errata for Robot Vision This is a list of known nontrivial bugs in Robot Vision 1986) by B.K.P. Horn, MIT Press, Cambridge, MA ISBN 0-6-08159-8 and McGraw-Hill, New York, NY ISBN 0-07-030349-5. If you
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 informationLecture Note 18: Duality
MATH 5330: Computational Methods of Linear Algebra 1 The Dual Problems Lecture Note 18: Duality Xianyi Zeng Department of Mathematical Sciences, UTEP The concept duality, just like accuracy and stability,
More informationDerivatives in 2D. Outline. James K. Peterson. November 9, Derivatives in 2D! Chain Rule
Derivatives in 2D James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 9, 2016 Outline Derivatives in 2D! Chain Rule Let s go back to
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationMATHEMATICS (2) Before you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IA Wednesday, 13 June, 2018 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: The paper has two sections, A and B. Section A contains
More informationCalculus - II Multivariable Calculus. M.Thamban Nair. Department of Mathematics Indian Institute of Technology Madras
Calculus - II Multivariable Calculus M.Thamban Nair epartment of Mathematics Indian Institute of Technology Madras February 27 Revised: January 29 Contents Preface v 1 Functions of everal Variables 1 1.1
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationMath 210, Final Exam, Spring 2012 Problem 1 Solution. (a) Find an equation of the plane passing through the tips of u, v, and w.
Math, Final Exam, Spring Problem Solution. Consider three position vectors (tails are the origin): u,, v 4,, w,, (a) Find an equation of the plane passing through the tips of u, v, and w. (b) Find an equation
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationMath Refresher Course
Math Refresher Course Columbia University Department of Political Science Fall 2007 Day 2 Prepared by Jessamyn Blau 6 Calculus CONT D 6.9 Antiderivatives and Integration Integration is the reverse of differentiation.
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationMath 127C, Spring 2006 Final Exam Solutions. x 2 ), g(y 1, y 2 ) = ( y 1 y 2, y1 2 + y2) 2. (g f) (0) = g (f(0))f (0).
Math 27C, Spring 26 Final Exam Solutions. Define f : R 2 R 2 and g : R 2 R 2 by f(x, x 2 (sin x 2 x, e x x 2, g(y, y 2 ( y y 2, y 2 + y2 2. Use the chain rule to compute the matrix of (g f (,. By the chain
More informationMath 361: Homework 1 Solutions
January 3, 4 Math 36: Homework Solutions. We say that two norms and on a vector space V are equivalent or comparable if the topology they define on V are the same, i.e., for any sequence of vectors {x
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationApplications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations
Applications of Symmetries and Conservation Laws to the Study of Nonlinear Elasticity Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationRobotics. Islam S. M. Khalil. November 15, German University in Cairo
Robotics German University in Cairo November 15, 2016 Fundamental concepts In optimal control problems the objective is to determine a function that minimizes a specified functional, i.e., the performance
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
More informationPhysics 129a Calculus of Variations Frank Porter Revision
Physics 129a Calculus of Variations 71113 Frank Porter Revision 171116 1 Introduction Many problems in physics have to do with extrema. When the problem involves finding a function that satisfies some
More informationMath 265H: Calculus III Practice Midterm II: Fall 2014
Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
More informationContents. 2 Partial Derivatives. 2.1 Limits and Continuity. Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v. 2.
Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v 260) Contents 2 Partial Derivatives 1 21 Limits and Continuity 1 22 Partial Derivatives 5 23 Directional Derivatives and the Gradient
More informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
More informationMathematical Analysis Outline. William G. Faris
Mathematical Analysis Outline William G. Faris January 8, 2007 2 Chapter 1 Metric spaces and continuous maps 1.1 Metric spaces A metric space is a set X together with a real distance function (x, x ) d(x,
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationCalculus Overview. f(x) f (x) is slope. I. Single Variable. A. First Order Derivative : Concept : measures slope of curve at a point.
Calculus Overview I. Single Variable A. First Order Derivative : Concept : measures slope of curve at a point. Notation : Let y = f (x). First derivative denoted f ʹ (x), df dx, dy dx, f, etc. Example
More informationISE I Brief Lecture Notes
ISE I Brief Lecture Notes 1 Partial Differentiation 1.1 Definitions Let f(x, y) be a function of two variables. The partial derivative f/ x is the function obtained by differentiating f with respect to
More informationMath Exam IV - Fall 2011
Math 233 - Exam IV - Fall 2011 December 15, 2011 - Renato Feres NAME: STUDENT ID NUMBER: General instructions: This exam has 16 questions, each worth the same amount. Check that no pages are missing and
More information(A) Opening Problem Newton s Law of Cooling
Lesson 55 Numerical Solutions to Differential Equations Euler's Method IBHL - 1 (A) Opening Problem Newton s Law of Cooling! Newton s Law of Cooling states that the temperature of a body changes at a rate
More informationBasic Theory of Linear Differential Equations
Basic Theory of Linear Differential Equations Picard-Lindelöf Existence-Uniqueness Vector nth Order Theorem Second Order Linear Theorem Higher Order Linear Theorem Homogeneous Structure Recipe for Constant-Coefficient
More informationVariational Principles
Part IB Variational Principles Year 218 217 216 215 214 213 212 211 21 218 Paper 1, Section I 4B 46 Variational Principles Find, using a Lagrange multiplier, the four stationary points in R 3 of the function
More informationMath 234 Final Exam (with answers) Spring 2017
Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial ifferential Equations «Viktor Grigoryan 3 Green s first identity Having studied Laplace s equation in regions with simple geometry, we now start developing some tools, which will lead
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs 2e: Section 3 - Exercise Page of 6 Exercise Carefully derive the equation of a string in a medium in which the resistance is proportional to the velocity Solution There are two ways (among
More informationIndian Institute of Technology Bombay
Indian Institute of Technology Bombay MA 5 Calculus Autumn SRG/AA/VDS/KSK Solutions and Marking Scheme for the Mid-Sem Exam Date: September, Weightage: 3 % Time: 8.3 AM.3 AM Max. Marks:. (i) Consider the
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationFaculty of Engineering, Mathematics and Science School of Mathematics
Faculty of Engineering, Mathematics and Science School of Mathematics GROUPS Trinity Term 06 MA3: Advanced Calculus SAMPLE EXAM, Solutions DAY PLACE TIME Prof. Larry Rolen Instructions to Candidates: Attempt
More informationA Tutorial on Active Contours
Elham Sakhaee March 4, 14 1 Introduction I prepared this technical report as part of my preparation for Computer Vision PhD qualifying exam. Here we discuss derivation of curve evolution function for some
More informationMotion Estimation (I)
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationLinear DifferentiaL Equation
Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace
More informationHouston Journal of Mathematics. c 2015 University of Houston Volume 41, No. 4, 2015
Houston Journal of Mathematics c 25 University of Houston Volume 4, No. 4, 25 AN INCREASING FUNCTION WITH INFINITELY CHANGING CONVEXITY TONG TANG, YIFEI PAN, AND MEI WANG Communicated by Min Ru Abstract.
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 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 informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationMath 115 ( ) Yum-Tong Siu 1. General Variation Formula and Weierstrass-Erdmann Corner Condition. J = x 1. F (x,y,y )dx.
Math 115 2006-2007 Yum-Tong Siu 1 General Variation Formula and Weierstrass-Erdmann Corner Condition General Variation Formula We take the variation of the functional F x,y,y dx with the two end-points,y
More informationMathematical Analysis II, 2018/19 First semester
Mathematical Analysis II, 208/9 First semester Yoh Tanimoto Dipartimento di Matematica, Università di Roma Tor Vergata Via della Ricerca Scientifica, I-0033 Roma, Italy email: hoyt@mat.uniroma2.it We basically
More informationDerivatives and Integrals
Derivatives and Integrals Definition 1: Derivative Formulas d dx (c) = 0 d dx (f ± g) = f ± g d dx (kx) = k d dx (xn ) = nx n 1 (f g) = f g + fg ( ) f = f g fg g g 2 (f(g(x))) = f (g(x)) g (x) d dx (ax
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
More informationDirectional Derivative and the Gradient Operator
Chapter 4 Directional Derivative and the Gradient Operator The equation z = f(x, y) defines a surface in 3 dimensions. We can write this as z f(x, y) = 0, or g(x, y, z) = 0, where g(x, y, z) = z f(x, y).
More informationMath 240 Calculus III
Calculus III Summer 2015, Session II Monday, August 3, 2015 Agenda 1. 2. Introduction The reduction of technique, which applies to second- linear differential equations, allows us to go beyond equations
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 0
CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 0 GENE H GOLUB 1 What is Numerical Analysis? In the 1973 edition of the Webster s New Collegiate Dictionary, numerical analysis is defined to be the
More informationChapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010
Section 6.1: Rational Expressions and Functions Definition of a rational expression Let u and v be polynomials. The algebraic expression u v is a rational expression. The domain of this rational expression
More informationSolutions to Sample Questions for Final Exam
olutions to ample Questions for Final Exam Find the points on the surface xy z 3 that are closest to the origin. We use the method of Lagrange Multipliers, with f(x, y, z) x + y + z for the square of the
More informationa k 0, then k + 1 = 2 lim 1 + 1
Math 7 - Midterm - Form A - Page From the desk of C. Davis Buenger. https://people.math.osu.edu/buenger.8/ Problem a) [3 pts] If lim a k = then a k converges. False: The divergence test states that if
More information2.8 Linear Approximations and Differentials
Arkansas Tech University MATH 294: Calculus I Dr. Marcel B. Finan 2.8 Linear Approximations and Differentials In this section we approximate graphs by tangent lines which we refer to as tangent line approximations.
More informationJim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt
Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain
More informationReview for the Final Exam
Math 171 Review for the Final Exam 1 Find the limits (4 points each) (a) lim 4x 2 3; x x (b) lim ( x 2 x x 1 )x ; (c) lim( 1 1 ); x 1 ln x x 1 sin (x 2) (d) lim x 2 x 2 4 Solutions (a) The limit lim 4x
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationFind the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x
Assignment 5 Name Find the indicated derivative. ) Find y(4) if y = sin x. ) A) y(4) = cos x B) y(4) = sin x y(4) = - cos x y(4) = - sin x ) y = (csc x + cot x)(csc x - cot x) ) A) y = 0 B) y = y = - csc
More informationJUST THE MATHS UNIT NUMBER INTEGRATION APPLICATIONS 10 (Second moments of an arc) A.J.Hobson
JUST THE MATHS UNIT NUMBER 13.1 INTEGRATION APPLICATIONS 1 (Second moments of an arc) by A.J.Hobson 13.1.1 Introduction 13.1. The second moment of an arc about the y-axis 13.1.3 The second moment of an
More informationMath512 PDE Homework 2
Math51 PDE Homework October 11, 009 Exercise 1.3. Solve u = xu x +yu y +(u x+y y/ = 0 with initial conditon u(x, 0 = 1 x. Proof. In this case, we have F = xp + yq + (p + q / z = 0 and Γ parameterized as
More informationIntegration - Past Edexcel Exam Questions
Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point
More informationMATH529 Fundamentals of Optimization Constrained Optimization I
MATH529 Fundamentals of Optimization Constrained Optimization I Marco A. Montes de Oca Mathematical Sciences, University of Delaware, USA 1 / 26 Motivating Example 2 / 26 Motivating Example min cost(b)
More informationMATH 312 Section 8.3: Non-homogeneous Systems
MATH 32 Section 8.3: Non-homogeneous Systems Prof. Jonathan Duncan Walla Walla College Spring Quarter, 2007 Outline Undetermined Coefficients 2 Variation of Parameter 3 Conclusions Undetermined Coefficients
More informationA Note on Integral Transforms and Partial Differential Equations
Applied Mathematical Sciences, Vol 4, 200, no 3, 09-8 A Note on Integral Transforms and Partial Differential Equations Adem Kılıçman and Hassan Eltayeb Department of Mathematics and Institute for Mathematical
More informationMA 441 Advanced Engineering Mathematics I Assignments - Spring 2014
MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff
More informationReview for the Final Exam
Calculus 3 Lia Vas Review for the Final Exam. Sequences. Determine whether the following sequences are convergent or divergent. If they are convergent, find their limits. (a) a n = ( 2 ) n (b) a n = n+
More informationLecture 1 Numerical methods: principles, algorithms and applications: an introduction
Lecture 1 Numerical methods: principles, algorithms and applications: an introduction Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationCITY UNIVERSITY. London
611.51 CITY UNIVERSITY London BSc Honours Degrees in Mathematical Science BSc Honours Degree in Mathematical Science with Finance and Economics BSc Honours Degree in Actuarial Science BSc Honours Degree
More informationExamples: Solving nth Order Equations
Atoms L. Euler s Theorem The Atom List First Order. Solve 2y + 5y = 0. Examples: Solving nth Order Equations Second Order. Solve y + 2y + y = 0, y + 3y + 2y = 0 and y + 2y + 5y = 0. Third Order. Solve
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 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 informationMATH20602 Numerical Analysis 1
M\cr NA Manchester Numerical Analysis MATH20602 Numerical Analysis 1 Martin Lotz School of Mathematics The University of Manchester Manchester, January 27, 2014 Outline General Course Information Introduction
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More informationLösning: Tenta Numerical Analysis för D, L. FMN011,
Lösning: Tenta Numerical Analysis för D, L. FMN011, 090527 This exam starts at 8:00 and ends at 12:00. To get a passing grade for the course you need 35 points in this exam and an accumulated total (this
More informationHybrid Analog-Digital Solution of Nonlinear Partial Differential Equations
Hybrid Analog-Digital Solution of Nonlinear Partial Differential Equations Yipeng Huang, Ning Guo, Kyle Mandli, Mingoo Seok, Yannis Tsividis, Simha Sethumadhavan Columbia University Hybrid Analog-Digital
More information