MATH 205C: STATIONARY PHASE LEMMA
|
|
- Gyles McCormick
- 5 years ago
- Views:
Transcription
1 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(ω) C(K) sup u ; we are interested in the asymptotic behavior as ω. If f has no critical points, then as ω, the exponential becomes highly oscillatory, and one expects that I(ω) 0 rapidly. Indeed, I(ω) C ω for all in this case. To see this, we note that if f never vanishes then, with f (x) 2 = j f 2, e iωf(x) = iω Leiωf(x), L = f (x) 2 j f, and so, with L t the transpose of L, i.e. L t = f j f (x) 2, with the relevant factors acting as multiplication operators, I(ω) = (Le iωf(.) )(x) u(x) dx = e iωf(x) (L t u)(x) dx, iω iω so by induction I(ω) = (iω) e iωf(x) ((L t ) u)(x) dx, leading to the conclusion that () I(ω) C(K, f, )ω sup D α u. α Of course, here we only needed u C with support in K for this estimate. For a moment, also consider non-real valued f with Im f 0 so e iωf(x) = e ω Im f(x) for ω. Then, as long as f 0, the above calculation goes through if we replace f by its complex conjugate in the definition of L, so () also holds. Moreover, if f does have some critical points, but at these Im f > 0, then a partition of unity argument (noting that the set of critical points is closed, as is the set of points where f is real, while supp u is compact) allows one to reduce consideration of the integral the two cases where either Im f 0 and f 0, which we just analyzed, or instead Im f > 0 (hence bounded below by a positive constant). In the latter case, one actually gets I(ω) e ω inf Im f u L, i.e. one has exponential decay. Thus, if Im f 0, and in addition f = 0 implies Im f > 0, then () still holds. Returning to real valued f, the interesting case is if f has some critical points, and the simplest setting is if these are non-degenerate, i.e. if f (x 0 ) = 0 implies that the Hessian f is invertible at x 0. Thus, one assumes that if f (x 0 ) = 0 then (2) f(x) = f(x 0 ) + 2 Q(x x 0, x x 0 ) + R(x, x 0 ),
2 2 where R(x, x 0 ) C x x 0 3, with Q a non-degenerate symmetric bilinear quadratic form; writing it as Q(x x 0, x x 0 ) = A(x x 0 ), x x 0 for the standard inner product on R n, A is a symmetric invertible linear operator on R n. The signature of Q, and of A, can be defined as the pair (, n ) where is the maximal dimension of a subspace on which Q is positive definite, or the number of positive eigenvalues of A. It is also useful to have a single number sign Q = sign A = (n ), the difference of the number of positive and negative eigenvalues of A. Note that non-degenerate critical points are isolated, for f (x) A(x x 0 ) vanishes quadratically at x 0, and A is invertible. Thus, by using a partition unity, when the critical points are non-degenerate, one can easily arrange that there is a single critical point x 0 over supp u, which one may arrange to be at 0 by a translation of the coordinates. In this case, one has f (x) C x x 0 for a suitable C > 0, so f (x) C x x 0. It is useful to observe that even if f has a non-degenerate critical point at 0, thus the above integration by parts argument breas down in general, it continues to wor if u vanishes to sufficiently high order at 0. It is convenient to be somewhat more general (due to the singularities the integration by parts argument induces), so suppose that u C (R n \ {0}), with support in K, with D α u C x 2 α for α. Then, by the divergence theorem in the last step, I(ω) = lim e x ρ iωf(x) u(x) dx = iω lim = ( iω lim e iωf(x) (L t u)(x) dx + x ρ x ρ x =ρ (Le iωf(.) )(x) u(x) dx f 2 ( j f ) ν j )e iωf(x) u(x) ds(x), where ν is the outward unit normal of R n \ {x : x ρ}. Now, by the assumptions on u, the estimate on f and the vanishing of f at 0, the integrand in the surface integral is bounded by Cρ, and thus the surface integral has limit 0 as ρ 0. Correspondingly, I(ω) = iω lim x ρ e iωf(x) (L t u)(x) dx = iω as L t u is bounded by the assumptions. Thus, I(ω) C(K, f)ω sup (D α u) x 2+ α. α e iωf(x) (L t u)(x) dx, Further, if u C (R n \ {0}), with support in K, with D α u C x 2 α for α, then L t u C (R n \ {0}), with support in K, with D α L t u C x 2( ) α for α. Thus, the argument can be applied iteratively to conclude that (3) I(ω) C(K, f)ω sup (D α u) x 2+ α. α We now show how to use this to reduce the case of general f to that of a quadratic form, 2 Q(x x 0, x x 0 ). With the notation of (2), let f s (x) = f(x 0 ) + 2 Q(x x 0, x x 0 ) + sr(x, x 0 ), s [0, ], so f 0 is quadratic (plus a constant), and f = f. Let I(ω, s) = e iωfs(x) u(x) dx.
3 3 Then differentiating in s 2 times yields I (2) (ω, s) = (iω) 2 e iωfs(x) R(x, x 0 ) 2 u(x) dx. Since R is C and vanishes cubically at x = x 0, v x0 (x) = R(x, x 0 ) 2 u(x) satisfies (D α v x0 )(x) C x x 0 6 α when u C 3 (R n ) with compact support in K and α 3. Thus, applying (3) with 3 in place of, we deduce that By Taylor s theorem, (4) I(ω, ) j<2 I (2) (ω, s) Cω. j! I(j) (ω, 0) sup s [0,] (2)! I(2) (ω, s), so modulo ω decay, I(ω) can be calculated merely by calculating I (j) (ω, 0) = (iω) j e iωf0(x) R(x, x 0 ) j u(x) dx (5) = (iω) j e iωf(x0) e iωq(x x0,x x0)/2 R(x, x 0 ) j u(x) dx for j < 2, i.e. we are reduced to the case of purely quadratic phase. We now want to compute integrals of the form I 0 (ω) = e iω 2 Ax,x u(x) dx, with u being C. Note that we could replace u by any function that has the same Taylor series to order 2 in view of our preceding estimates, and one way to proceed would be to compute explicitly such integrals. Instead, we rewrite using Parseval s formula I 0 (ω) = (F e iω 2 A.,. )(ξ)(fu)(ξ) dξ We mae a general definition: if P C (R n ) with D α P (ξ) C α ξ Mα, i.e. if P is polynomially bounded with all derivatives, let P (D)u = F (P (.)(Fu)(.)), so P (D) maps S to itself and S to itself. Thus, we have (6) I 0 (ω) = (2π) n ((F e iω 2 A.,. )(D)u)(0). Our first tas is to compute the inverse Fourier transform of the imaginary Gaussian. First, recall that and (F(t e t2 /2 ))(τ) = 2πe τ 2 /2, (F (t e t2 /2 ))(τ) = 2π e τ 2 /2. A change of variables gives that for α > 0 (F (t e αt2 /2 ))(τ) = 2πα e τ 2 /(2α). Now, both sides are analytic functions of α with values in S(R) in Re α > 0, and both sides are continuous functions of α with values in S (R) in Re α 0, α 0, so the formula continues to hold there (with the square root being the standard one
4 4 for positive reals, with a branch cut along the negative reals). In particular, when α = βi, β > 0, we deduce that (F (t e ±iβt2 /2 ))(τ) = 2πβ e ±iπ/4 e iτ 2 /(2β). Now, if A is a non-degenerate symmetric linear map on R n, then there is an orthogonal transformation O diagonalizing A as Λ = OAO t = diag(λ j ), with rows of O given by orthonormal eigenvectors e j of A, and with λ j 0 the eigenvalues of A, so Ax, x = j λ j e j, x 2, i.e. in coordinates y j = e j, x, y = Ox, eeping in mind that O has determinant of absolute value, (2π) n e ix ξ e i Ax,x /2 dx = (2π) n e iy Oξ e i j λjy2 j /2 dy Thus, (7) (F e iω 2 A.,. )(ξ) = = Now, by Taylor s formula, f(y) = (2π) n/2 λ /2 ei(π/4)... λ n /2 j sign λj e i j (Oξ)2 j /(2λj) = (2π) n/2 det A /2 ei(π/4) sign A e i (OAO t ) Oξ,Oξ /2 = (2π) n/2 det A /2 ei(π/4) sign A e i A ξ,ξ /2. (2πω) n/2 det A /2 ei(π/4) sign A e i ω A ξ,ξ /2 f (j) (x) (y x) j (y x) + j! ( )! we have for w R (or indeed for w with Im w 0) e iw ( iw) j j! ( )! w 0 ( t) f () (ty + ( t)x) dt, Thus, with P (.) = B.,., using Plancherel s theorem, 0 ( t) dt = w! e ip (D) u ( ip (D)) j u j! L 2! P (D) u L 2.. Now, we would lie to replace the L 2 norm on the left by just the absolute value of the value of the argument at 0 in order to obtain an expansion for (6) (while strengthening the norm on the right, of course); this goal is reached if we can replace the L 2 norm by the L norm, for the argument is Schwartz when u is in C c (R n ). But by the Sobolev embedding, for s > n/2 integer, there is C > 0 such that v L C( v L 2 + D α v L 2), α =s
5 5 so as P (D) and e ip (D) commute with D α (for they are all Fourier multipliers) e ip (D) u ( ip (D)) j u j! L ( C e ip (D) u ( ip (D)) j u j! L 2 + D α( e ip (D) u ( ip (D)) j u ) L 2 j!) α =s ( C e ip (D) u ( ip (D)) j u j! L 2 + e ip (D) (D α u) ( ip (D)) j Dα u j! α =s C ( P (D) u L 2 + ) P (D) D α u L 2! α =s C! B ( u L 2 + ) D α u L 2 C B α =s β 2+s D β u L 2, where in the penultimate step we used that P (ξ) B ξ 2 for ξ R n so P (ξ) Fu L 2 B ξ 2 Fu L 2 and thus P (D) u L 2 B u L 2 by Plancherel, and then in the last step we expanded = n j= D2 j. Here C depends on the Sobolev inequality parameters C, s and on only. Applying this with B = ω A /2, (6) and (7) give (8) (2π) n/2 I 0 (ω) ω n/2 det A /2 ei(π/4) sign A ( iω A D, D ) j u j! (0) π n/2 det A /2 (2ω) n/2 C A D β u L 2, β 2+s Combining this with (4) and (5), we deduce that for L j differential operators of order 2j, and L 0 the identity, (9) I(ω) e iωf(x0) (2π) n/2 sign f (x ω n/2 det f ei(π/4) 0) ω j L (x 0 ) /2 j u(x 0 ) ω n/2 C f (x 0 ) D β u L 2, β 2+s which is the stationary phase lemma. Notice that if (D α u)(x 0 ) vanishes for α 2l +, l <, then the first l terms in the sum on the left hand side vanish, so the expansion starts with ω n/2 l. Assuming f(x 0 ) = 0, since ω l ( ωi)(ω) l = ω l e iωf(x) f(x) l u(x) dx, and f vanishes quadratically at x 0, derivatives of I in ω possess a similar expansion, starting with a multiple if ω n/2. Factoring out e iωf(x0) in general, (0) ω l ( ω(e l if(x0)ω (2π) n/2 I))(ω) sign f (x ω n/2 det f ei(π/4) 0) ω j L (x 0 ) /2 l,j u(x 0 ) ω n/2 C l f (x 0 ) D β u L 2, β 2+2l+s L 2 )
6 6 where L l,j are differential operators of the form L l+j f(x) l. One can integrate this expression to obtain the already-nown expansion for I, which shows that the coefficients in the expansion of ω(e l if(x0)ω I) are the term-by-term differentiated coefficients of e if(x0)ω I. Finally, ( ωi)(ω) l = ( ω(e l if(x0)ω (e if(x0)ω I))(ω) can be computed by the product rule. One can also allow parametric dependence on another variable y R m, i.e. for f C (R n+m ) with Im f 0, consider the integral I(ω, y) = e iωf(x,y) u(x, y) dx, u Cc (R n+m ). R n As before I decays rapidly in ω, uniformly in y unless f has critical points in x at which it is real, i.e. unless there exists (x 0, y 0 ) such that f x(x 0, y 0 ) = 0 and Im f(x 0, y 0 ) = 0. In addition, differentiation under the integral sign shows that in this case I is C in (ω, y), with Dy α DωI l still rapidly decreasing. If f is real, and has some critical point (x 0, y 0 ) in x which is non-degenerate, i.e. f xx is invertible, so d x f x,... d x f x n are linearly independent at (x 0, y 0 ), then the implicit function theorem guarantees that the joint zero set of f x,..., f (x n ) is a C graph over a neighborhood of y 0 in R m, i.e. there is a C function X defined on a neighborhood of y 0 such that in a neighborhood of (x 0, y 0 ), the critical points of f are given by (X(y), y), and these are non-degenerate. Thus, the previous arguments are applicable, uniformly in y, so the stationary phase expansion (0) holds with the y dependence added in. Writing I(ω, y) = e iωf(x(y),y) R n e iωf (x,y) u(x, y) dx, F (x, y) = f(x, y) f(x(y), y), notice that, with yj hand side, standing for the jth component of the second slot on the right ( yj F )(x, y) = ( yj f)(x, y) ( yj f)(x(y), y) (( x f)(x(y), y)) X (y), y j and thus vanishes at x = X(y) since that is a critical point of f in x. Since x f, and thus x F, have a non-degenerate zero at x = X(y), we actually obtain ( yj F )(x, y) = G j (x, y)( x F )(x, y). Differentiation under the integral sign gives, with Ĩ(ω, y) = e iωf (x,y) u(x, y) dx, R n yj Ĩ(ω, y) = (iω yj F (x, y))e iωf (x,y) u(x, y) dx + e iωf (x,y) yj u(x, y) dx. R n R n Now, the second term has a stationary phase expansion as before. On the other hand, the first term is G j (x, y)(iω x F (x, y))e iωf (x,y) u(x, y) dx R n = G j (x, y)( x e iωf (x,y) ) u(x, y) dx = e iωf (x,y) x (G j (x, y)u(x, y)) dx, R n R n so the stationary phase expansion is again applicable. Thus, yj I possesses a stationary phase expansion. Since the expansion can be integrated in y j, this proves
7 7 that the coefficients of the expansion of the derivative are the differentiated coefficients of the expansion of I (including the statement that the latter are differentiable). Iterating the argument shows that ω l y α ωi l still possesses an expansion as in (0), with the expansion given by term-by-term differentiation. We finally discuss an interpretation of the stationary phase result by partially compactifying (or bordifying) the (ω, y)-space. Thus, consider the new variable h = ω (0, ]. Then the stationary phase lemma, together with ω ω = h h and that ω l ω l can be rewritten as a linear combination of (ω ω ) j for j l, and similarly for h l h l, shows that I(h, y) = e if(x,y)/h u(x, y) dx = e if(x(y),y)/h h n/2 J(h, y), with h l l h α y J j h j Lj,l,α u, where means that if one sums over j <, then one obtains an error terms bounded by C l,α, h, and the series on the right is the term-by-term derivative of h j Lj,0,0 u. But notice that h l j<l hj Lj,0,0 u = 0, so in particular, taing = l, h l h l α y J C l,α,l h l, i.e. h l α y J C l,α,l for all l and α, i.e. all partial derivatives of the C function J on (0, ] h R m y are bounded. This implies that J has a unique C extension to [0, ] R m : uniqueness is automatic as (0, ] R m is dense in [0, ] R m. On the other hand, h l α y J extends continuously to h = 0 since l h α y J(h 0, y) = l h α y J(, y) l+ h h 0 α y J(h, y) dh, and l+ h y α J(h, y) is bounded, and then repeated integration shows that for l l, α α, h l α y J extends as well (which of course would follow from the same argument with (l, α) replaced by (l, α )), with extension given by the appropriate integral of l+ h y α J(h, y), so in particular this extension is differentiable (l l, α α ) times. As l and α are arbitrary, we deduce that the extension of J we defined is C ; in view of the uniqueness, one usually simply writes J for the extension as well. A slight generalization is then obtained by allowing u to depend on h as well, i.e. consider u Cc (R n R m [0, ]), and consider J(h, y) = e if(x,y)/h u(x, y, h) dx. Regarding the h in u as one of the parameters y, we have then J(h, y) = Ĵ(h, y, h), where Ĵ(h, y, h) = e if(x,y)/h u(x, y, h) dx = e if(x(y),y)/h h n/2 Ĵ (h, y, h) with Ĵ C ([0, ] R m [0, ]), so restricting to h = h yields with J C ([0, ] R m ). J(h, y) = e if(x(y),y)/h h n/2 J(h, y),
LECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationWe denote the derivative at x by DF (x) = L. With respect to the standard bases of R n and R m, DF (x) is simply the matrix of partial derivatives,
The derivative Let O be an open subset of R n, and F : O R m a continuous function We say F is differentiable at a point x O, with derivative L, if L : R n R m is a linear transformation such that, for
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationMathematics Department Stanford University Math 61CM/DM Inner products
Mathematics Department Stanford University Math 61CM/DM Inner products Recall the definition of an inner product space; see Appendix A.8 of the textbook. Definition 1 An inner product space V is a vector
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
More informationx λ ϕ(x)dx x λ ϕ(x)dx = xλ+1 λ + 1 ϕ(x) u = xλ+1 λ + 1 dv = ϕ (x)dx (x))dx = ϕ(x)
1. Distributions A distribution is a linear functional L : D R (or C) where D = C c (R n ) which is linear and continuous. 1.1. Topology on D. Let K R n be a compact set. Define D K := {f D : f outside
More informationOscillatory integrals
Chapter Oscillatory integrals. Fourier transform on S The Fourier transform is a fundamental tool in microlocal analysis and its application to the theory of PDEs and inverse problems. In this first section
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More information(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y
() Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open
More informationTOOLS FROM HARMONIC ANALYSIS
TOOLS FROM HARMONIC ANALYSIS BRADLY STADIE Abstract. The Fourier transform can be thought of as a map that decomposes a function into oscillatory functions. In this paper, we will apply this decomposition
More informationThe method of stationary phase
Chapter The method of stationary phase In this chapter we introduce a very useful analytical tool which will allow us to find asymptotic expansions for integrals that cannot, in many cases, be calculated
More informationPremiliminary Examination, Part I & II
Premiliminary Examination, Part I & II Monday, August 9, 16 Dept. of Mathematics University of Pennsylvania Problem 1. Let V be the real vector space of continuous real-valued functions on the closed interval
More informationTEST CODE: PMB SYLLABUS
TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
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 informationQuadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.
Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationx R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).
Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x,
More informationAlgebra II. Paulius Drungilas and Jonas Jankauskas
Algebra II Paulius Drungilas and Jonas Jankauskas Contents 1. Quadratic forms 3 What is quadratic form? 3 Change of variables. 3 Equivalence of quadratic forms. 4 Canonical form. 4 Normal form. 7 Positive
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More information1 Inner Product and Orthogonality
CSCI 4/Fall 6/Vora/GWU/Orthogonality and Norms Inner Product and Orthogonality Definition : The inner product of two vectors x and y, x x x =.., y =. x n y y... y n is denoted x, y : Note that n x, y =
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationLinear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions
Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Main problem of linear algebra 2: Given
More informationCalculus 2502A - Advanced Calculus I Fall : Local minima and maxima
Calculus 50A - Advanced Calculus I Fall 014 14.7: Local minima and maxima Martin Frankland November 17, 014 In these notes, we discuss the problem of finding the local minima and maxima of a function.
More information17 The functional equation
18.785 Number theory I Fall 16 Lecture #17 11/8/16 17 The functional equation In the previous lecture we proved that the iemann zeta function ζ(s) has an Euler product and an analytic continuation to the
More informationSingular Integrals. 1 Calderon-Zygmund decomposition
Singular Integrals Analysis III Calderon-Zygmund decomposition Let f be an integrable function f dx 0, f = g + b with g Cα almost everywhere, with b
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, February 2, Time Allowed: Two Hours Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Scholarships Screening Test Saturday, February 2, 2008 Time Allowed: Two Hours Maximum Marks: 40 Please read, carefully, the instructions on the following
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationReal Variables # 10 : Hilbert Spaces II
randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,
More informationPOINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO
POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO HART F. SMITH AND MACIEJ ZWORSKI Abstract. We prove optimal pointwise bounds on quasimodes of semiclassical Schrödinger
More informationLecture 4: Fourier Transforms.
1 Definition. Lecture 4: Fourier Transforms. We now come to Fourier transforms, which we give in the form of a definition. First we define the spaces L 1 () and L 2 (). Definition 1.1 The space L 1 ()
More informationCALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =
CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.
More informationLecture 8 : Eigenvalues and Eigenvectors
CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationj=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p.
LINEAR ALGEBRA Fall 203 The final exam Almost all of the problems solved Exercise Let (V, ) be a normed vector space. Prove x y x y for all x, y V. Everybody knows how to do this! Exercise 2 If V is a
More informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationDot Products. K. Behrend. April 3, Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem.
Dot Products K. Behrend April 3, 008 Abstract A short review of some basic facts on the dot product. Projections. The spectral theorem. Contents The dot product 3. Length of a vector........................
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationPutzer s Algorithm. Norman Lebovitz. September 8, 2016
Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationWhitney s Extension Problem for C m
Whitney s Extension Problem for C m by Charles Fefferman Department of Mathematics Princeton University Fine Hall Washington Road Princeton, New Jersey 08544 Email: cf@math.princeton.edu Supported by Grant
More informationDetailed Proof of The PerronFrobenius Theorem
Detailed Proof of The PerronFrobenius Theorem Arseny M Shur Ural Federal University October 30, 2016 1 Introduction This famous theorem has numerous applications, but to apply it you should understand
More informationNotes on basis changes and matrix diagonalization
Notes on basis changes and matrix diagonalization Howard E Haber Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, CA 95064 April 17, 2017 1 Coordinates of vectors and matrix
More informationMATH 235. Final ANSWERS May 5, 2015
MATH 235 Final ANSWERS May 5, 25. ( points) Fix positive integers m, n and consider the vector space V of all m n matrices with entries in the real numbers R. (a) Find the dimension of V and prove your
More informationMath 24 Spring 2012 Sample Homework Solutions Week 8
Math 4 Spring Sample Homework Solutions Week 8 Section 5. (.) Test A M (R) for diagonalizability, and if possible find an invertible matrix Q and a diagonal matrix D such that Q AQ = D. ( ) 4 (c) A =.
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationWeek Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,
Math 051 W008 Margo Kondratieva Week 10-11 Quadratic forms Principal axes theorem Text reference: this material corresponds to parts of sections 55, 8, 83 89 Section 41 Motivation and introduction Consider
More informationHW3 - Due 02/06. Each answer must be mathematically justified. Don t forget your name. 1 2, A = 2 2
HW3 - Due 02/06 Each answer must be mathematically justified Don t forget your name Problem 1 Find a 2 2 matrix B such that B 3 = A, where A = 2 2 If A was diagonal, it would be easy: we would just take
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationARE211 FINAL EXAM DECEMBER 12, 2003
ARE2 FINAL EXAM DECEMBER 2, 23 This is the final exam for ARE2. As announced earlier, this is an open-book exam. Try to allocate your 8 minutes in this exam wisely, and keep in mind that leaving any questions
More informationMORE NOTES FOR MATH 823, FALL 2007
MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3-space, with a point a S where S looks smooth, i.e., without any fold or cusp or self-crossing, we can intuitively define the tangent
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationHere we used the multiindex notation:
Mathematics Department Stanford University Math 51H Distributions Distributions first arose in solving partial differential equations by duality arguments; a later related benefit was that one could always
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More information1. General Vector Spaces
1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationOn rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro
On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More information