18.303: Introduction to Green s functions and operator inverses
|
|
- Lauren Nash
- 6 years ago
- Views:
Transcription
1 8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential operator Â, and are led to the concept of a Green s function G(, ) [the ( ), matri element of  ]. Although G is a perfectly well-defined function, in trying to construct it and define its properties we are drawn ineorably to something that appears less well defined: ÂG seems to be some kind of function that is infinitely concentrated at a point, a delta function δ( ) that does not have a proper definition as an ordinary function. In these notes, we show how to solve for G in an eample d problem without resorting to delta functions, by a somewhat awkward process that involves taking a limit after solving the PDE. In principle, we could always apply such contortions, but the inconvenience of this approach will lead us in future lectures to a new definition of function in which delta functions are well defined. Review: Matri inverses The inverse A of a matri A is already familiar to you from linear algebra, where A A = AA = I. et s think about what this means. Each column j of A has a special meaning: AA = I, so A (column j of A ) = (column j of I) = e j, where e j is the unit vector (,,...,,,,..., ) T which has a in the j-th row. That is, (column j of A ) = A e j, the solution to Au = f with e j on the right-hand-side. If we are solving Au = f, the statement u = A f can be written as a superposition of columns of A : u = A f = j (column j of A )f j = A j f j e j. Since j f je j = f is just writing f in the basis of the e j s, we can interpret A f as writing f in the basis of the e j unit vectors and then summing (superimposing) the solutions for each e j. Going one step further, and writing: u i = ( A f ) i = j (A ) i,j f j, we can now have a simple physical interpretation of the matri elements of A : (A ) i,j is the effect (solution) at i of a source at j (a right-hand-side e j ), [and column j of A is the effect (solution) everywhere from the source at j]. et us consider a couple of eamples drawn from the finite-difference matrices for a discretized 2 with Dirichlet () boundary conditions. Recall that one interpretation of the problem 2 u = f is that u is a displacement of a stretched string (d) or membrane(2d) and that f is a force density (pressure). Figure (top) shows several columns of A for 2 2 in d, and it is very much like what you might intuitively epect if you pressed on a string at one point. Figure (bottom) shows two column of A for 2 on a square 2d domain (e.g. a stretched square drum head), and again matches the intuition for what happens if you press at one point. In both cases, the columns are, as epected, the solution from a source e j at one point j (and the minimum points in the plot are indeed the points j).
2 6.5.5 columns of A / 3 3 one A column another A column y y.5 Figure : Columns of A where A is a discretized 2 with Dirichlet boundary conditions. Top: several columns in d (Ω = [, ] = [, ]). Bottom: two columns in 2d (Ω = [, ] [, ]). In both d and 2d, the location of the minimum corresponds to the inde of the column: this is the effect of a unit-vector source or force = at that position (and = elsewhere). 2
3 . Inverses in other bases This is not the only way to write down A, of course, because e j are only one possible choice of basis. For eample, suppose that q,..., q m is any orthonormal basis for our vector space, and that u j solves Au j = q j. Then we can write any f = j q j(q j f) and hence the solution u to Au = f by u = A f = j u j(q j f). In other words, we are writing A = j u jq j = j (A q j )q j in terms of this basis. In matri form, if Q is the (unitary) matri whose columns are the q j, then we are simply writing A = (A Q)Q, since QQ = I. [Similarly, if the columns of B are any basis, not necessarily orthonormal, we can write A = (A B)B : we solve any right-hand-side f by a superposition of the solutions for the columns of B with coefficients B f.] For eample, if A is Hermitian (A = A ), then we can choose the q j to be orthonormal eigenvectors of A (Aq j = λ j q j ), in which case u j = q j /λ j and u = A f = j u j(q j f) = j q j(q j f)/λ j. But this is just writing A = QΛ Q in terms of the (hopefully) familiar diagonalization A = QΛQ. Operator inverses Now we would like to consider the inverse  of an operator Â, if such a thing eists. Suppose N(Â) = {} so that Âu = f has a unique solution u (for all f in some suitable space of functions).  is then whatever operator gives the solution u =  f from any f. This must clearly be a linear operator, since if u solves Âu = f and u 2 solves Âu 2 = f 2 then u = αu + βu 2 solves  = αf + βf 2. In fact, we already know how to write down Â, for self-adjoint Â, in terms of an orthonormal basis of eigenfunctions u n :  f = n u n u n, f /λ n. But now we would like to write  in a different way: analogous to the entries (A ) ij above, we would like to write the solution u() as a superposition of effects at from sources at for all (where just means another point, not a derivative): ˆ u() =  f = d d G(, )f( ), () Ω where G(, ) is the Green s function or the fundamental solution to the PDE. We need to figure out what G is (and whether it eists at all), and how to solve for it. This equation is the direct analogue of u i = j (A ) i,j f j for the matri case above, with sums turned into integrals, so this leads us to interpret G(, ) as the matri element ( ),. So, in analogy to the matri case, we hope that G(, ) is the effect at from a source at. But what does this mean for a PDE?. ÂG and a delta function If G eists, one way of getting at its properties is to consider what PDE G(, ) might satisfy (as opposed to the integral equation above). In particular, we can try to apply the equation Âu = f, noting that  is an operator acting on and not on so that it commutes with the Ω dd integral: ˆ Âu = d d [ÂG(, )]f( ) = f(). Ω This is telling us that ÂG(, ) must be some function that, when integrated against any f( ), gives us f(). What must such an ÂG(, ) look like? It must somehow be sharply peaked around so that all of the integral s contributions come from f(), and in particular the integral must give some sort of average of f( ) in an infinitesimal neighborhood of. How could we get such an average? In one dimension for simplicity, consider the function δ () = { < otherwise (2) 3
4 δ () δ()?? + +? Figure 2: eft: A finite localized function δ () from eq. 2: a rectangle of area on [, ]. Right: We would like to take the + limit to obtain a Dirac delta function δ(), an infinitesimally wide and infinitely high spike of area, but this is not well defined according to our ordinary understanding of function. as depicted in figure 2(left): a bo of width and height /, with area. The integral ˆ d δ ( )f( ) = ˆ f( )d is just the average of f( ) in [, + ]. As, this goes to f() (for continuous f), which is precisely the sort of integral we want. That is, we want ÂG(, ) = lim δ ( ) = δ( ), (3) where δ() is a Dirac delta function that is a bo infinitely high and infinitesimally narrow with area around =, as depicted in figure 2(right). Just one problem: this is not a function, at least as we currently understand function, and the definition of δ() here doesn t actually make much rigorous sense. In fact, it is possible to define δ(), and even lim δ (), in a perfectly rigorous fashion, but it requires a new level of abstraction: we will have to redefine our notion of what a function is, replacing it with something called a distribution. We aren t ready for that, yet, though we should see how far we can get with ordinary functions first. In fact, we will see that it is perfectly possible to obtain G as an ordinary function, even though ÂG will not eist in the ordinary sense (G will have some discontinuity or singularity that prevents it from being classically differentiable). But the process is fairly indirect and painful: once we obtain G the hard way, we will appreciate how much nicer it will be to have δ(), even at the price of an etra layer of abstraction..2 G by a limit of superpositions et s consider the same problem from another angle. We saw for matrices that we could interpret the columns of A as the response to unit vectors e j, and A f as representing f in the basis of unit vectors and then superimposing (summing) the solutions for each e j. et s try to do the same thing here, at least approimately. Consider one dimension for simplicity, with a domain Ω = [, ]. et s approimate f() by a piecewise-constant function as shown in figure 3: a sum of N + rectangles of width = /(N + ). Algebraically, that could be written as a sum in terms of precisely the δ () functions from (2) above: N f() f N () = f(n )δ ( n ), (4) n= since δ ( n ) is just a bo of height on [n, (n + ) ]. In the limit N ( ), f N () clearly goes to the eact f() (for continuous f). But we don t want to 4
5 f() f N () N Figure 3: A function f() (left) can be approimated by a piecewise-constant f N () (right) as in eq. (4), where each segment [n, (n + ) ] takes the value f(n ). take that limit yet. Instead, we will find the solution u N to Âu N = f N, and then take the N limit. For finite N, f N is the superposition of a finite number of bo functions that are localized sources. Our strategy, similar to the matri case, will be to solve [analogous to (3)] Âg n () = δ ( n ) (5) for each of these right-hand sides to obtain solutions g n (), to write u N () = n f(n )g n(), and then to take the limit to obtain an integral where / u() = lim u N() = lim f(n )g n () = N n= ˆ G(, )f( )d, G(, ) = lim g / (). (6) We will find that, even though δ does not have a well-defined limit (at least, not according to the familiar definition of function ), g m does have a perfectly well-defined limit G as an ordinary function. 2 Eample: Green s function of d 2 /d 2 As an eample, consider the familiar  = d2 d 2 on [, ] with Dirichlet boundaries u() = u() =. Following the approach above, we want to solve equation (5) and then take the limit to obtain G(, ). That is, we want to solve g n() = δ ( n ) (7) with the boundary conditions g n () = g n () =, and then take the limit (6) of g n to get G. As depicted in figure 4, we will solve this by breaking the solution up into three regions: (I) [, n ); (II) [n, (n + ) ); and (III) [(n + ), ]. In region II, = /, and the most general g m () with this property is the quadratic function 2 + γ + κ for some unknown constants γ and κ to be determined. In I and III, g n =, so the solutions must be straight lines; combined with the boundary conditions at and, g n 2 this gives a solution of the form: α g n () = γ + κ β( ) < n [n, (n + ) ], (8) > (n + ) for constants α, β, γ, κ to be determined. These unknowns are determined by imposing 5
6 δ ( n ) I II III n (n+) Figure 4: When solving g n() = δ ( n ), we will divide the solution into three regions: regions I and III, where the right-hand side is zero (and g n is a line), and region II where the right-hand side is / (and g n is a quadratic). continuity conditions at the interfaces between the regions, to make sure that the solutions match up with one another. Equation (7) tells us that both g n and g n must be continuous in order to obtain a finite piecewise-continuous g n. Thus, we obtain four equations for our four unknowns: αn = n 2 + γn + κ, 2 α = n + γ, β( [n + ] ) = (n + ) 2 + γ(n + ) + κ, 2 β = (n + ) + γ. After straightforward but tedious algebra, letting = n, one can obtain the solutions: β = + 2 α = β, γ = + α, κ = 2 2. The resulting function g m () is plotted for three values of in figure 5(left). It looks much like what one might epect: linear functions in regions I and III that are smoothly patched together by a quadratic-shaped kink in region II. Physically, this can be interpreted as the shape of a stretched string when you press on it with a localized pressure δ ( ). We are now in a position to take the limit, keeping = n fied. Region II disappears, and we are left with β {( G(, ) =, in region III and α ) < = ( ) ) {( < ( ), in region I: a pleasingly symmetrical function [whose symmetry G(, ) = G(, ) is called reciprocity and, we will eventually show, derives from the self-adjointness of Â]. This function is plotted in figure 5(right), and is continuous but with a discontinuity in its slope at =. [Indeed, it looks just like what we got from the discrete aplacian in figure (top), ecept with a sign flip since we are looking at d 2 /d 2.] Physically, this intuitively corresponds to what happens when you press on a stretched string at one point. Thus, we have obtained the Green s function G(, ), which is indeed a perfectly welldefined, ordinary function. As promised in the previous section, ÂG does not eist according to our ordinary definition of derivative, not even piecewise, since G is discontinuous: the second derivative is infinity at. Correspondingly, δ does not have a well-defined limit 6
7 .25 =.5.25 =.5 g() for = =.25 =.75 eact G(, ).2.5. =.25 = (=) (=) Figure 5: Green s function G(, ) for  = 2 / 2 on [, ] with Dirichlet boundary conditions, for =.25,.5, and.75. eft: approimate Green s function, using our finite-width δ () as the right-hand-side. Right: eact Green s function, from limit. as an ordinary function. So, even though the right-hand side of Âg n = δ ( n ) did not have a well-defined limit, the solution g n did, and indeed we can then write the solution u for any f eactly by our integral () of G. If only there were an easier way to obtain G Deriving the slope discontinuity of G Why is G/ discontinuous at =? It is easy to see this just by integrating both sides of equation (7) for g n (): ˆ (n+) n g n() d = g n(n ) g n(n + ) = ˆ (n+) n δ () d =. Thus, the slope g drops by when going from region I to region II. The limit + of these epressions (the integrals) is perfectly well defined (even though the limit of the integrands is not), and gives G G =, = + = so the slope G/ must drop discontinuously by at =. In fact, we can use this to easily solve for G(, ). For < (region { I) and for > (region III) we must have ÂG = and hence G is a line: G(, α < ) = β( ) >. We have two unknown constants α and β, but we have two equations: continuity of G, α = β( ); and our slope condition α ( β) =. Combining these, one can quickly find that α = and β = as above! 3 To infinity... and beyond! It is perfectly possible to work with Green s functions and  in this way, putting finite δ -like sources on the right-hand sides and only taking the + limit or its equivalent after solving the PDE. Physically, this corresponds to putting in a source in a small but non-infinitesimal region (e.g. a force in a small region for a stretched string, or a small ball of charge for electrostatics), and then taking the limit of the solution to obtain the limiting behavior as the source becomes infinitely concentrated at a point (e.g. to obtain the displacement for a string pressed at one point, or the potential for a point charge). However, this is all terribly cumbersome and complicated. 7
8 It would be so much simpler if we could have an infinitely concentrated source at a point to start with. For eample, when you study freshman mechanics or electromagnetism, you start with point masses and point charges and only later do you consider etended distributions of mass or charge. But to do that in a rigorous way seems problematic because δ() = lim δ () apparently does not eist as a function in any sense that you have learned it does not have a value for every. It turns out that the problem is not with δ(), the problem is with the definition of a function surprisingly, classical functions don t really correspond to what we want in many practical circumstances, and they are often not quite what we want to deal with for PDEs. A clue is provided in section 2., where we saw that even though our integrands might not always have well-defined limits, our integrals did. If we can only define a function in terms of its integrals rather than by its values at points, that might fi our problems. In fact, we can do something very much like this, which leads to the concept of a distribution (or generalized function) or a weak solution. At the cost of an increase in the level of abstraction, we can completely eliminate the tortured contortions that were required above (and only get worse in higher dimensions) to avoid delta functions. 8
Solve Wave Equation from Scratch [2013 HSSP]
1 Solve Wave Equation from Scratch [2013 HSSP] Yuqi Zhu MIT Department of Physics, 77 Massachusetts Ave., Cambridge, MA 02139 (Dated: August 18, 2013) I. COURSE INFO Topics Date 07/07 Comple number, Cauchy-Riemann
More information1 5 π 2. 5 π 3. 5 π π x. 5 π 4. Figure 1: We need calculus to find the area of the shaded region.
. Area In order to quantify the size of a 2-dimensional object, we use area. Since we measure area in square units, we can think of the area of an object as the number of such squares it fills up. Using
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationLinear Algebra, Summer 2011, pt. 2
Linear Algebra, Summer 2, pt. 2 June 8, 2 Contents Inverses. 2 Vector Spaces. 3 2. Examples of vector spaces..................... 3 2.2 The column space......................... 6 2.3 The null space...........................
More informationMatrices, Row Reduction of Matrices
Matrices, Row Reduction of Matrices October 9, 014 1 Row Reduction and Echelon Forms In the previous section, we saw a procedure for solving systems of equations It is simple in that it consists of only
More information2 Generating Functions
2 Generating Functions In this part of the course, we re going to introduce algebraic methods for counting and proving combinatorial identities. This is often greatly advantageous over the method of finding
More informationPage 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03
Page 5 Lecture : Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 008/10/0 Date Given: 008/10/0 Inner Product Spaces: Definitions Section. Mathematical Preliminaries: Inner
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationUNCONSTRAINED OPTIMIZATION PAUL SCHRIMPF OCTOBER 24, 2013
PAUL SCHRIMPF OCTOBER 24, 213 UNIVERSITY OF BRITISH COLUMBIA ECONOMICS 26 Today s lecture is about unconstrained optimization. If you re following along in the syllabus, you ll notice that we ve skipped
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationThe binary entropy function
ECE 7680 Lecture 2 Definitions and Basic Facts Objective: To learn a bunch of definitions about entropy and information measures that will be useful through the quarter, and to present some simple but
More informationGenerating Function Notes , Fall 2005, Prof. Peter Shor
Counting Change Generating Function Notes 80, Fall 00, Prof Peter Shor In this lecture, I m going to talk about generating functions We ve already seen an example of generating functions Recall when we
More informationOne Solution Two Solutions Three Solutions Four Solutions. Since both equations equal y we can set them equal Combine like terms Factor Solve for x
Algebra Notes Quadratic Systems Name: Block: Date: Last class we discussed linear systems. The only possibilities we had we 1 solution, no solution or infinite solutions. With quadratic systems we have
More information1 Introduction. 2 Solving Linear Equations
1 Introduction This essay introduces some new sets of numbers. Up to now, the only sets of numbers (algebraic systems) we know is the set Z of integers with the two operations + and and the system R of
More informationRational Expressions & Equations
Chapter 9 Rational Epressions & Equations Sec. 1 Simplifying Rational Epressions We simply rational epressions the same way we simplified fractions. When we first began to simplify fractions, we factored
More informationCOUNCIL ROCK HIGH SCHOOL MATHEMATICS. A Note Guideline of Algebraic Concepts. Designed to assist students in A Summer Review of Algebra
COUNCIL ROCK HIGH SCHOOL MATHEMATICS A Note Guideline of Algebraic Concepts Designed to assist students in A Summer Review of Algebra [A teacher prepared compilation of the 7 Algebraic concepts deemed
More informationα x x 0 α x x f(x) α x x α x ( 1) f(x) x f(x) x f(x) α x = α x x 2
Quadratic speedup for unstructured search - Grover s Al- CS 94- gorithm /8/07 Spring 007 Lecture 11 01 Unstructured Search Here s the problem: You are given an efficient boolean function f : {1,,} {0,1},
More information1 Rational Exponents and Radicals
Introductory Algebra Page 1 of 11 1 Rational Eponents and Radicals 1.1 Rules of Eponents The rules for eponents are the same as what you saw earlier. Memorize these rules if you haven t already done so.
More informationBasic Linear Algebra in MATLAB
Basic Linear Algebra in MATLAB 9.29 Optional Lecture 2 In the last optional lecture we learned the the basic type in MATLAB is a matrix of double precision floating point numbers. You learned a number
More information22A-2 SUMMER 2014 LECTURE 5
A- SUMMER 0 LECTURE 5 NATHANIEL GALLUP Agenda Elimination to the identity matrix Inverse matrices LU factorization Elimination to the identity matrix Previously, we have used elimination to get a system
More informationLinear Equations in Linear Algebra
1 Linear Equations in Linear Algebra 1.1 SYSTEMS OF LINEAR EQUATIONS LINEAR EQUATION,, 1 n A linear equation in the variables equation that can be written in the form a a a b 1 1 2 2 n n a a is an where
More informationAP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals
AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals 8. Basic Integration Rules In this section we will review various integration strategies. Strategies: I. Separate
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More informationMATH 116, LECTURES 10 & 11: Limits
MATH 6, LECTURES 0 & : Limits Limits In application, we often deal with quantities which are close to other quantities but which cannot be defined eactly. Consider the problem of how a car s speedometer
More informationReview of similarity transformation and Singular Value Decomposition
Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm
More informationIntroduction. So, why did I even bother to write this?
Introduction This review was originally written for my Calculus I class, but it should be accessible to anyone needing a review in some basic algebra and trig topics. The review contains the occasional
More informationCoordinates for Projective Planes
Chapter 8 Coordinates for Projective Planes Math 4520, Fall 2017 8.1 The Affine plane We now have several eamples of fields, the reals, the comple numbers, the quaternions, and the finite fields. Given
More informationSometimes the domains X and Z will be the same, so this might be written:
II. MULTIVARIATE CALCULUS The first lecture covered functions where a single input goes in, and a single output comes out. Most economic applications aren t so simple. In most cases, a number of variables
More informationNIT #7 CORE ALGE COMMON IALS
UN NIT #7 ANSWER KEY POLYNOMIALS Lesson #1 Introduction too Polynomials Lesson # Multiplying Polynomials Lesson # Factoring Polynomials Lesson # Factoring Based on Conjugate Pairs Lesson #5 Factoring Trinomials
More information1.1 Limits & Continuity
1.1 Limits & Continuity What do you see below? We are building the House of Calculus, one side at a time... and we need a solid FOUNDATION. Page 1 of 11 Eample 1: (Calculator) For f ( ) (a) fill in the
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationLecture 10: Powers of Matrices, Difference Equations
Lecture 10: Powers of Matrices, Difference Equations Difference Equations A difference equation, also sometimes called a recurrence equation is an equation that defines a sequence recursively, i.e. each
More information1 Introduction. 2 Solving Linear Equations. Charlotte Teacher Institute, Modular Arithmetic
1 Introduction This essay introduces some new sets of numbers. Up to now, the only sets of numbers (algebraic systems) we know is the set Z of integers with the two operations + and and the system R of
More information3.3 Limits and Infinity
Calculus Maimus. Limits Infinity Infinity is not a concrete number, but an abstract idea. It s not a destination, but a really long, never-ending journey. It s one of those mind-warping ideas that is difficult
More informationCALCULUS I. Integrals. Paul Dawkins
CALCULUS I Integrals Paul Dawkins Table of Contents Preface... ii Integrals... Introduction... Indefinite Integrals... Computing Indefinite Integrals... Substitution Rule for Indefinite Integrals... More
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationC. Finding roots of trinomials: 1st Example: x 2 5x = 14 x 2 5x 14 = 0 (x 7)(x + 2) = 0 Answer: x = 7 or x = -2
AP Calculus Students: Welcome to AP Calculus. Class begins in approimately - months. In this packet, you will find numerous topics that were covered in your Algebra and Pre-Calculus courses. These are
More information3.3.1 Linear functions yet again and dot product In 2D, a homogenous linear scalar function takes the general form:
3.3 Gradient Vector and Jacobian Matri 3 3.3 Gradient Vector and Jacobian Matri Overview: Differentiable functions have a local linear approimation. Near a given point, local changes are determined by
More informationMATH 521, WEEK 2: Rational and Real Numbers, Ordered Sets, Countable Sets
MATH 521, WEEK 2: Rational and Real Numbers, Ordered Sets, Countable Sets 1 Rational and Real Numbers Recall that a number is rational if it can be written in the form a/b where a, b Z and b 0, and a number
More informationk is a product of elementary matrices.
Mathematics, Spring Lecture (Wilson) Final Eam May, ANSWERS Problem (5 points) (a) There are three kinds of elementary row operations and associated elementary matrices. Describe what each kind of operation
More informationQuadratic Equations Part I
Quadratic Equations Part I Before proceeding with this section we should note that the topic of solving quadratic equations will be covered in two sections. This is done for the benefit of those viewing
More informationPractice problems for Exam 3 A =
Practice problems for Exam 3. Let A = 2 (a) Determine whether A is diagonalizable. If so, find a matrix S such that S AS is diagonal. If not, explain why not. (b) What are the eigenvalues of A? Is A diagonalizable?
More informationDefinition 2.3. We define addition and multiplication of matrices as follows.
14 Chapter 2 Matrices In this chapter, we review matrix algebra from Linear Algebra I, consider row and column operations on matrices, and define the rank of a matrix. Along the way prove that the row
More informationRoberto s Notes on Differential Calculus Chapter 1: Limits and continuity Section 7. Discontinuities. is the tool to use,
Roberto s Notes on Differential Calculus Chapter 1: Limits and continuity Section 7 Discontinuities What you need to know already: The concept and definition of continuity. What you can learn here: The
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More information1 Continuity and Limits of Functions
Week 4 Summary This week, we will move on from our discussion of sequences and series to functions. Even though sequences and functions seem to be very different things, they very similar. In fact, we
More informationc 1 v 1 + c 2 v 2 = 0 c 1 λ 1 v 1 + c 2 λ 1 v 2 = 0
LECTURE LECTURE 2 0. Distinct eigenvalues I haven t gotten around to stating the following important theorem: Theorem: A matrix with n distinct eigenvalues is diagonalizable. Proof (Sketch) Suppose n =
More informationLine Integrals and Path Independence
Line Integrals and Path Independence We get to talk about integrals that are the areas under a line in three (or more) dimensional space. These are called, strangely enough, line integrals. Figure 11.1
More informationQuadratic Functions: Global Analysis
Chapter 12 Quadratic Functions: Global Analysis The Essential Question, 191 Concavity-sign, 192 Slope-sign, 194 Etremum, 194 Height-sign, 195 Bounded Graph, 197 0-Concavity Location, 199 0-Slope Location,
More informationCore Connections Algebra 2 Checkpoint Materials
Core Connections Algebra 2 Note to Students (and their Teachers) Students master different skills at different speeds. No two students learn eactly the same way at the same time. At some point you will
More informationFirst Variation of a Functional
First Variation of a Functional The derivative of a function being zero is a necessary condition for the etremum of that function in ordinary calculus. Let us now consider the equivalent of a derivative
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationBasic methods to solve equations
Roberto s Notes on Prerequisites for Calculus Chapter 1: Algebra Section 1 Basic methods to solve equations What you need to know already: How to factor an algebraic epression. What you can learn here:
More informationConceptual Explanations: Radicals
Conceptual Eplanations: Radicals The concept of a radical (or root) is a familiar one, and was reviewed in the conceptual eplanation of logarithms in the previous chapter. In this chapter, we are going
More informationIn the real world, objects don t just move back and forth in 1-D! Projectile
Phys 1110, 3-1 CH. 3: Vectors In the real world, objects don t just move back and forth in 1-D In principle, the world is really 3-dimensional (3-D), but in practice, lots of realistic motion is 2-D (like
More informationDealing with Rotating Coordinate Systems Physics 321. (Eq.1)
Dealing with Rotating Coordinate Systems Physics 321 The treatment of rotating coordinate frames can be very confusing because there are two different sets of aes, and one set of aes is not constant in
More information2 3 x = 6 4. (x 1) 6
Solutions to Math 201 Final Eam from spring 2007 p. 1 of 16 (some of these problem solutions are out of order, in the interest of saving paper) 1. given equation: 1 2 ( 1) 1 3 = 4 both sides 6: 6 1 1 (
More informationSection - 9 GRAPHS. (a) y f x (b) y f x. (c) y f x (d) y f x. (e) y f x (f) y f x k. (g) y f x k (h) y kf x. (i) y f kx. [a] y f x to y f x
44 Section - 9 GRAPHS In this section, we will discuss graphs and graph-plotting in more detail. Detailed graph plotting also requires a knowledge of derivatives. Here, we will be discussing some general
More informationA Quick Algebra Review
1. Simplifying Epressions. Solving Equations 3. Problem Solving 4. Inequalities 5. Absolute Values 6. Linear Equations 7. Systems of Equations 8. Laws of Eponents 9. Quadratics 10. Rationals 11. Radicals
More informationACCUPLACER MATH 0311 OR MATH 0120
The University of Teas at El Paso Tutoring and Learning Center ACCUPLACER MATH 0 OR MATH 00 http://www.academics.utep.edu/tlc MATH 0 OR MATH 00 Page Factoring Factoring Eercises 8 Factoring Answer to Eercises
More informationResonance and response
Chapter 2 Resonance and response Last updated September 20, 2008 In this section of the course we begin with a very simple system a mass hanging from a spring and see how some remarkable ideas emerge.
More informationDifferential Equations
This document was written and copyrighted by Paul Dawkins. Use of this document and its online version is governed by the Terms and Conditions of Use located at. The online version of this document is
More informationIn economics, the amount of a good x demanded is a function of the price of that good. In other words,
I. UNIVARIATE CALCULUS Given two sets X and Y, a function is a rule that associates each member of X with eactly one member of Y. That is, some goes in, and some y comes out. These notations are used to
More informationMITOCW ocw f99-lec23_300k
MITOCW ocw-18.06-f99-lec23_300k -- and lift-off on differential equations. So, this section is about how to solve a system of first order, first derivative, constant coefficient linear equations. And if
More informationLecture I: Vectors, tensors, and forms in flat spacetime
Lecture I: Vectors, tensors, and forms in flat spacetime Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: September 28, 2011) I. OVERVIEW The mathematical description of curved
More informationLESSON 35: EIGENVALUES AND EIGENVECTORS APRIL 21, (1) We might also write v as v. Both notations refer to a vector.
LESSON 5: EIGENVALUES AND EIGENVECTORS APRIL 2, 27 In this contet, a vector is a column matri E Note 2 v 2, v 4 5 6 () We might also write v as v Both notations refer to a vector (2) A vector can be man
More information1 = I I I II 1 1 II 2 = normalization constant III 1 1 III 2 2 III 3 = normalization constant...
Here is a review of some (but not all) of the topics you should know for the midterm. These are things I think are important to know. I haven t seen the test, so there are probably some things on it that
More informationII. The Machinery of Quantum Mechanics
II. The Machinery of Quantum Mechanics Based on the results of the experiments described in the previous section, we recognize that real experiments do not behave quite as we expect. This section presents
More informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationHarmonic Oscillator with raising and lowering operators. We write the Schrödinger equation for the harmonic oscillator in one dimension as follows:
We write the Schrödinger equation for the harmonic oscillator in one dimension as follows: H ˆ! = "!2 d 2! + 1 2µ dx 2 2 kx 2! = E! T ˆ = "! 2 2µ d 2 dx 2 V ˆ = 1 2 kx 2 H ˆ = ˆ T + ˆ V (1) where µ is
More informationSlope Fields: Graphing Solutions Without the Solutions
8 Slope Fields: Graphing Solutions Without the Solutions Up to now, our efforts have been directed mainly towards finding formulas or equations describing solutions to given differential equations. Then,
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationEC5555 Economics Masters Refresher Course in Mathematics September 2014
EC5555 Economics Masters Refresher Course in Mathematics September 4 Lecture Matri Inversion and Linear Equations Ramakanta Patra Learning objectives. Matri inversion Matri inversion and linear equations
More informationChapter 1. Linear Equations
Chapter 1. Linear Equations We ll start our study of linear algebra with linear equations. Lost of parts of mathematics rose out of trying to understand the solutions of different types of equations. Linear
More informationThe following are generally referred to as the laws or rules of exponents. x a x b = x a+b (5.1) 1 x b a (5.2) (x a ) b = x ab (5.
Chapter 5 Exponents 5. Exponent Concepts An exponent means repeated multiplication. For instance, 0 6 means 0 0 0 0 0 0, or,000,000. You ve probably noticed that there is a logical progression of operations.
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationarxiv: v2 [cond-mat.str-el] 20 Apr 2015
Gauging time reversal symmetry in tensor network states ie Chen, 2 and Ashvin Vishwanath 2 Department of Physics and Institute for Quantum Information and Matter, California Institute of echnology, Pasadena,
More informationSolving Equations. Lesson Fifteen. Aims. Context. The aim of this lesson is to enable you to: solve linear equations
Mathematics GCSE Module Four: Basic Algebra Lesson Fifteen Aims The aim of this lesson is to enable you to: solve linear equations solve linear equations from their graph solve simultaneous equations from
More informationCIS 2033 Lecture 5, Fall
CIS 2033 Lecture 5, Fall 2016 1 Instructor: David Dobor September 13, 2016 1 Supplemental reading from Dekking s textbook: Chapter2, 3. We mentioned at the beginning of this class that calculus was a prerequisite
More informationRoberto s Notes on Linear Algebra Chapter 4: Matrix Algebra Section 7. Inverse matrices
Roberto s Notes on Linear Algebra Chapter 4: Matrix Algebra Section 7 Inverse matrices What you need to know already: How to add and multiply matrices. What elementary matrices are. What you can learn
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More informationLecture 4: Constructing the Integers, Rationals and Reals
Math/CS 20: Intro. to Math Professor: Padraic Bartlett Lecture 4: Constructing the Integers, Rationals and Reals Week 5 UCSB 204 The Integers Normally, using the natural numbers, you can easily define
More informationx + 2y + 3z = 8 x + 3y = 7 x + 2z = 3
Chapter 2: Solving Linear Equations 23 Elimination Using Matrices As we saw in the presentation, we can use elimination to make a system of linear equations into an upper triangular system that is easy
More informationMain topics for the First Midterm Exam
Main topics for the First Midterm Exam The final will cover Sections.-.0, 2.-2.5, and 4.. This is roughly the material from first three homeworks and three quizzes, in addition to the lecture on Monday,
More information1 Dirac Notation for Vector Spaces
Theoretical Physics Notes 2: Dirac Notation This installment of the notes covers Dirac notation, which proves to be very useful in many ways. For example, it gives a convenient way of expressing amplitudes
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationThe Particle in a Box
Page 324 Lecture 17: Relation of Particle in a Box Eigenstates to Position and Momentum Eigenstates General Considerations on Bound States and Quantization Continuity Equation for Probability Date Given:
More informationACCUPLACER MATH 0310
The University of Teas at El Paso Tutoring and Learning Center ACCUPLACER MATH 00 http://www.academics.utep.edu/tlc MATH 00 Page Linear Equations Linear Equations Eercises 5 Linear Equations Answer to
More informationUniversity of Colorado at Colorado Springs Math 090 Fundamentals of College Algebra
University of Colorado at Colorado Springs Math 090 Fundamentals of College Algebra Table of Contents Chapter The Algebra of Polynomials Chapter Factoring 7 Chapter 3 Fractions Chapter 4 Eponents and Radicals
More informationRichard Cleve David R. Cheriton School of Computer Science Institute for Quantum Computing University of Waterloo
CS 497 Frontiers of Computer Science Introduction to Quantum Computing Lecture of http://www.cs.uwaterloo.ca/~cleve/cs497-f7 Richard Cleve David R. Cheriton School of Computer Science Institute for Quantum
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationGROUP THEORY PRIMER. New terms: so(2n), so(2n+1), symplectic algebra sp(2n)
GROUP THEORY PRIMER New terms: so(2n), so(2n+1), symplectic algebra sp(2n) 1. Some examples of semi-simple Lie algebras In the previous chapter, we developed the idea of understanding semi-simple Lie algebras
More informationMatrix inversion and linear equations
Learning objectives. Matri inversion and linear equations Know Cramer s rule Understand how linear equations can be represented in matri form Know how to solve linear equations using matrices and Cramer
More informationProbability, Expectation Values, and Uncertainties
Chapter 5 Probability, Epectation Values, and Uncertainties As indicated earlier, one of the remarkable features of the physical world is that randomness is incarnate, irreducible. This is mirrored in
More informationCS168: The Modern Algorithmic Toolbox Lecture #8: How PCA Works
CS68: The Modern Algorithmic Toolbox Lecture #8: How PCA Works Tim Roughgarden & Gregory Valiant April 20, 206 Introduction Last lecture introduced the idea of principal components analysis (PCA). The
More information0.1. Linear transformations
Suggestions for midterm review #3 The repetitoria are usually not complete; I am merely bringing up the points that many people didn t now on the recitations Linear transformations The following mostly
More informationCourse. Print and use this sheet in conjunction with MathinSite s Maclaurin Series applet and worksheet.
Maclaurin Series Learning Outcomes After reading this theory sheet, you should recognise the difference between a function and its polynomial epansion (if it eists!) understand what is meant by a series
More informationLecture 3: Sizes of Infinity
Math/CS 20: Intro. to Math Professor: Padraic Bartlett Lecture 3: Sizes of Infinity Week 2 UCSB 204 Sizes of Infinity On one hand, we know that the real numbers contain more elements than the rational
More informationA BRIEF REVIEW OF ALGEBRA AND TRIGONOMETRY
A BRIEF REVIEW OF ALGEBRA AND TRIGONOMETR Some Key Concepts:. The slope and the equation of a straight line. Functions and functional notation. The average rate of change of a function and the DIFFERENCE-
More information2 Instantaneous Rate of Change: The Derivative
0 Chapter Instantaneous Rate of Change: The Derivative ¾º½ Instantaneous Rate of Change: The Derivative Ì ÐÓÔ Ó ÙÒØ ÓÒ Suppose that y is a function of, say y f() It is often necessary to know how sensitive
More information