LECTURE 5: THE METHOD OF STATIONARY PHASE
|
|
- Gerard Sherman
- 5 years ago
- Views:
Transcription
1 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 α = x α x αn n. x n. n. α = α x αn D α = D α D αn For any ϕ C ( ), Gradient vector ϕ = ( ϕ,, n ϕ). Hessian matrix d 2 ϕ = [ i j ϕ] n n. Schwartz functions. Definition.. A function ϕ C ( ) is called a Schwartz function (or a rapidly decreasing function) if () sup x α β ϕ < for all multi-indices α and β. We will denote the set of all Schwartz functions by S ( ), or by S for simplicity. Obviously S is an infinitely dimensional vector space. It is easy to see C 0 ( ) S ( ) L p ( ) for any p. For any λ C with Re(λ) > 0, ϕ(x) = e λ x 2 S ( ). If ϕ is Schwartz, so are x α D β ϕ, D α x β ϕ, where α, β are any multi-indices. Recall that a semi-norm on a vector space V is function : V R so that for all λ C and all u, v V, () (Absolute homogeneity) λv = λ v. (2) (Triangle inequality) u + v u + v. On S one can define, for each pair of multi-indices α, β a semi-norm (2) ϕ α,β = sup x α β ϕ.
2 2 LECTURE 5: THE METHOD OF STATIONARY PHASE We say that ϕ j ϕ in S if ϕ j ϕ α,β 0 as j for all α, β. The topology induced by these semi-norms turns S into a Fréchet space. Equivalently, one can define a metric d on S via (3) d(ϕ, ψ) = k=0 2 k ϕ ψ k + ϕ ψ k, where for each k 0, ϕ k = max sup x α D β ϕ α + β k x is a norm on S. Then we have a metric topology on S which coincides with the one described above. Exercise. Under the topology described above, C 0 ( ) is dense in S. Fourier transform on S. Let ϕ S be a Schwartz function. By definition its Fourier transform Fϕ is (4) Fϕ(ξ) = ˆϕ(ξ) = e ix ξ ϕ(x)dx. For example, one can prove (5) F(e x 2 2 ) = (2π) n 2 e ξ 2 2. By defintion we have the obvious inequality Fϕ L ϕ L. It is easy to check that for any ϕ S and any multi-index α, (6) F(x α ϕ) = ( ) α D α ξ (F(ϕ)) and (7) F(D α x ϕ) = ξ α F(ϕ). As a consequence, we see ϕ S = Fϕ S. From definition it is also easy to check (8) Fϕ(ξ)ψ(ξ)dξ = ϕ(x)fψ(x)dx. Lemma.2. Suppose ϕ S. () For any λ R +, if we let ϕ λ = ϕ(λx), then ϕ λ S and (9) Fϕ λ (ξ) = λ n Fϕ( ξ λ ).
3 LECTURE 5: THE METHOD OF STATIONARY PHASE 3 (2) For any a, if we let T a ϕ(x) = ϕ(x + a), then T a ϕ S and (0) (FT a ϕ)(ξ) = e ia ξ Fϕ(ξ). Now we can prove Theorem.3. F : S S is an isomorphism with inverse () F ψ(x) = e ix ξ ψ(ξ)dξ. (2π) n Proof. We have Fϕ(ξ)ψ(λξ) = λ n ϕ(x)fψ( x λ )(x)dx = ϕ(λx)fψ(x)dx. Letting λ 0, we get ψ(0) Fϕ(ξ)dξ = ϕ(0) Fψ(x)dx. In particular, if we take ψ(ξ) = e ξ 2 2, then we get (2) ϕ(0) = Fϕ(ξ)dξ. (2π) n Finally we replace ϕ by T a ϕ in the above formula, then ϕ(a) = T a ϕ(0) = e ia ξ Fϕ(ξ)dξ. (2π) n As a consequence, we get ϕ 2 L 2 = (2π) n ˆϕ 2 L 2. The Fourier transform of the Gaussian e 2 xt Qx. Theorem.4. Let Q be any real, symmetric, positive definite n n matrix, then (3) F(e 2 xt Qx )(ξ) = (2π)n/2 (det Q) e /2 2 ξt Q ξ. Proof. Let s first assume n =, so that Q = q is a positive number. By definition, F(e 2 xt qx )(ξ) = e 2 qx2 ixξ dx. R
4 4 LECTURE 5: THE METHOD OF STATIONARY PHASE Taking ξ derivative, we get d dξ F(e 2 xt qx )(ξ) = ix e 2 qx2 ixξ dx = i d R q R dx (e 2 qx2 )e ixξ dx = i e d 2 qx2 q dx (e ixξ )dx = ξ q F(e 2 xt qx )(ξ). It follows where R F(e 2 xt qx )(ξ) = Ce ξ2 2q, C = F(e 2 xt qx )(0) = R e 2 qx2 dx = 2π q. Next suppose Q is a diagonal matrix. Then the left hand side of (3) is a product of n one-dimensional integrals, and the result follows from the one-dimensional case. For the general case, one only need to choose an orthogonal matrix O so that O T QO = D = diag(λ,, λ n ), where λ i s are the eigenvalues of Q, then if we set y = O x, F(e 2 xt Qx )(ξ) = 2 xt Qx ix ξ dx = 2 yt Dy iy O T ξ dy e = (2π)n/2 (det D) e /2 e 2 ξt OD O T ξ = (2π)n/2 (det Q) e 2 ξt Q ξ. /2 Now suppose Q is an n n complex matrix such that Q = Q T and Re(Q) is positive definite. We observe that the integral on the left hand side of (3) is an analytic function of the entries Q ij = Q ji of Q in the region ReQ > 0. The same is true for the right hand side of (3), as long as we choose det /2 Q to be the branch such that det /2 Q > 0 for real positive definite Q. So we end up with Theorem.5. Let Q be any n n complex matrix such that Q = Q T and Re(Q) is positive definite, then (4) F(e 2 xt Qx )(ξ) = (2π)n/2 (det Q) e /2 2 ξt Q ξ. Another way to describe this choice of det /2 : Since Q = Q T, for any complex vector w C n we have Re( w T Qw) = w T (ReQ)w. It follows that all eigenvalues λ,, λ n of Q has positive real part. Then where λ /2 j (det Q) /2 = λ /2 λ /2 n, is the square root of λ j that has positive real part.
5 LECTURE 5: THE METHOD OF STATIONARY PHASE 5 The Fourier transform on S. Definition.6. Any linear continuous map u : S C is called a tempered distribution. The space of tempered distributions is denoted by S. So S is the dual of S. The topology of the space S is defined to be th weak-* topology, so that u j u in S if For example, u j (ϕ) u(ϕ), ϕ S. For any a, the Dirac distribution δ a : S C defined by δ a (ϕ) = ϕ(a) is a tempered distribution. More generally, for any a C n and any multi-index α, the map u : S C defined by u(ϕ) = D α ϕ(a) is a tempered distribution. For any bounded continuous function ψ on defines a tempered distribution by u(ϕ) = ϕ(x)ψ(x)dx. Note that if ψ S, then the distribution defined by the above formula is non-zero unless ψ = 0. It follows that S S. Definition.7. Let u S. We define () x α u S via (x α u)(ϕ) = u(x α ϕ). (2) D α u S via (D α u)(ϕ) = ( ) α u(d α ϕ). (3) Fu S via (Fu)(ϕ) = u(fϕ). Now let Q be any real, symmetric, non-singular n n matrix. Then the function e i 2 xt Qx is not a Schwartz function, but we can think of it as a distribution. We would like to calculate its Fourier transform. Recall that the signature of Q is (5) sgn(q) = N + (Q) N (Q), where N ± (Q) = the number of positive/negative eigenvalues of Q. Theorem.8. For any real, symmetric, non-singular n n matrix Q, (6) F(e i 2 xt Qx ) = (2π)n/2 e i π 4 sgn(q) e i det Q 2 ξt Q ξ. 2
6 6 LECTURE 5: THE METHOD OF STATIONARY PHASE Proof. Let s first assume n =, so that Q = q 0 is a non-zero real number. For any ε > 0 we let q ε = q + iε. Applying theorem.5 we get F(e i 2 qεx2 ) = F(e 2 (ε iq)x2 ) = (2π)/2 (ε iq) e 2 (ε iq)ξ2, /2 the square root is described after theorem.5. Note that when ε 0+, we have { qe (ε iq) /2 iπ 4 q > 0 iπ qe 4 q < 0 = q e iπ 4 sgn(q), so the conclusion follows. For general n, one just proceed via the diagonalization method as before. 2. The method of stationary phase Oscillatory integrals. The method of stationary phase is a method for estimating the small behavior of the oscillatory integrals (7) I = e i ϕ(x) a(x)dx, where ϕ C ( ) is called the amplitude, and a Cc ( ) is called the phase. For simplicity we will assume ϕ is real-valued. To illustrate, let s start with two extremal cases: Suppose ϕ(x) = c is a constant. Then I = e ic/ A, where A = a(x)dx is a constant independent of. Suppose n = and ϕ(x) = x. Then by definition I = F(a)( ). Since a is compactly supported, F(a) is Schwartz. It follows I = O( N ), Intuitively, in the second case the exponential exp(i x ) oscillates fast in any interval of x for small, so that many cancellations take place and thus we get a function rapidly decreasing in. This is in fact the case at any point x which is not a critical point of ϕ. On the other hand, near a critical point of x the phase function ϕ doesn t change much, i.e. looks like a constant, so that we are in case. According to these intuitive observation, we expect that the main contributions to I should arise from the critical points of the phase function ϕ. We will use the following notions: Definition 2.. Let f = f(). N.
7 LECTURE 5: THE METHOD OF STATIONARY PHASE 7 () We say f k=0 a k k if for each nonnegative integer N, (8) f N a k k = O( N+ ), 0. k=0 (2) We say f = O( ) if f 0, i.e. (9) f = O( N ), 0 for all N. So in the second extremal case we discussed above, we can write I = O( ). Remark. We don t require the series k=0 a k k to be convergent! Non-stationary phase. Proposition 2.2. Suppose the phase function ϕ has no critical point in a neighborhood of the support of a. Then (20) I = O( ). Proof. Let χ be a smooth cut-off function which is identically one on the support of a so that ϕ has no critical point on the support of χ. Then a(x) = χ(x)a(x), and is a smooth function. Define an operator L by χ(x) ϕ(x) 2 Then it is easy to check It follows Lf(x) = L(e i ϕ(x) )(x) = χ I = i L(e i χ(x) ϕ(x) 2 j ϕ(x) j f(x). ϕ 2 j j j ϕ i jϕ ϕ(x) ei = i ϕ(x) χei. ϕ(x) )a(x)dx = e i ϕ(x) (L a)(x)dx. i Repeating this N times, we get ( ) N I = e i ϕ(x) (L ) N a(x)dx = O( N ). i
8 8 LECTURE 5: THE METHOD OF STATIONARY PHASE Morse lemma. Recall that a point p is called a critical point of a smooth function ϕ if (2) ϕ(p) = 0. A critical point p of ϕ is called non-degenerate if the matrix d 2 ϕ(p) is non-degenerate: [ ] (22) det d 2 2 ϕ ϕ(p) = det (p) 0. x i x j A smooth function all of whose critical points are non-degenerate is called a Morse function. Theorem 2.3 (Morse lemma, form ). Let ϕ C ( ). Suppose p is a nondegenerate critical point of ϕ. Then there exists a neighborhood U of 0, a neighborhood V of p and a diffeomorphism f : U V so that f(0) = p and (23) f ϕ(x) = ϕ(p) + 2 (x2 + + x 2 r x 2 r+ x 2 n), where r is the number of positive eigenvalues of the Hessian matrix d 2 ϕ(p). We shall prove the following form of the Morse lemma which is obviously equivalent to the above statement: Theorem 2.4 (Morse lemma, form 2). Let ϕ 0 and ϕ be smooth functions on such that ϕ i (0) = 0, ϕ 0 (0) = ϕ (0) = 0. d 2 ϕ 0 (0) = d 2 ϕ (0) is non-degenerate. Then there exist neighborhoods U 0 and U of 0 and a diffeomorphism f : U 0 U such that f(0) = 0 and f ϕ = ϕ 0. Proof. Let ϕ t = ( t)ϕ 0 + tϕ. Applying Moser s trick, it is enough to find a time-dependent vector field X t such that X t (0) = 0 and (24) L Xt ϕ t = ϕ 0 ϕ. Let X t = A j (x, t) x j, where A j (0, t) = 0. According to the second and the third conditions, [ 2 ϕ t x i x j (0)] is non-degenerate. It follows that the system of functions { ϕ t x i i =, 2,, n} form a system of coordinates near 0 with ϕt x i (0) = 0. According to the next lemma, one can find functions A ij (x, t) with A ij (0, t) = 0 so that A j (x, t) = A ij (x, t) ϕ t x i.
9 LECTURE 5: THE METHOD OF STATIONARY PHASE 9 Inserting this into the equation (24), we get (25) Aij (x, t) ϕ t ϕ t = ϕ 0 ϕ. x i x j On the other hand, the three conditions implies that ϕ 0 ϕ vanishes to 2nd order at 0, so according the next lemma, one can find functions B ijk (x, t) so that ϕ 0 ϕ = B ijk (x, t) ϕ t x i ϕ t x j ϕ t x k. In other words, the equation (25) is solvable with A ij (0, t) = 0. This proves the existence of the diffeomorphism f. Lemma 2.5. Let ϕ be a smooth function with ϕ(0) = 0 and α ϕ(0) = 0 for all α < k. Then there exists smooth functions ϕ α so that (26) ϕ(x) = x α ϕ α (x). α =k Proof. For k =, one just take ϕ ϕ i (x) = (tx)dt. 0 x i For larger k, apply the above formula and induction. By exactly the same method, one can prove the following Morse lemma with parameters: Theorem 2.6 (Morse lemma with parameter). Let ϕ s i C ( ) be two family of smooth functions, depending smoothly on the parameter s R k. Suppose ϕ s i (0) = 0, ϕ s 0(0) = ϕ s (0) = 0, d 2 ϕ s 0(0) = d 2 ϕ s (0) are non-degenerate. Then there exist an ε > 0, a neighborhood U of 0 and for all s < ε an open embedding f s : U, depending smoothly on s, so that f s (0) = 0 and f s ϕ s = ϕ s 0. The stationary phase formula for quadratic phase. Associated to any polynomial p(ξ) = p α ξ α we have a constant coefficient differential operator p(d) = p α D α. According to the Fourier transform formula that for any a S, we immediately get F(D α a) = ξ α Fa
10 0 LECTURE 5: THE METHOD OF STATIONARY PHASE Lemma 2.7. For any a S, (27) F(p(D)a)(ξ) = p(ξ)fa(ξ). In particular, associated to any n n symmetric matrix A = (a ij ) one has a quadratic function p A (ξ) = ξ T Aξ = a ij ξ i ξ j and thus a second order constant coefficient differential operator (28) p A (D) := a kl D k D l. It follows that F(p A (D)a)(ξ) = (ξ T Aξ)Fa(ξ). More generally, for any nonnegative integer k, (29) F(p A (D) k a)(ξ) = (ξ T Aξ) k Fa(ξ). Suppose the phase function is a quadratic function, i.e. ϕ(x) = 2 xt Qx for some symmetric real n n matrix Q. Then the only critical point of ϕ is the origin x = 0. More over, the non-degeneracy condition implies that Q is non-singular. Theorem 2.8. Let Q be a real, symmetric, non-singular n n matrix. For any a Cc ( ), one has (30) e i 2 xt Qx a(x)dx (2π) n/2 ei π 4 sgn(q) k ( i ) k R det Q n /2 k! 2 p Q (D) a(0). Proof. We start with the formula Fϕ(ξ)ψ(ξ)dξ = Recall from theorem.5 that if ϕ(x) = e i 2 xt Qx, then k ϕ(x)fψ(x)dx. Fϕ(ξ) = (2π)n/2 e i π 4 sgn(q) e i det(q) /2 2 ξt Q ξ. It follows e i 2 xt Qx a(x)dx = n/2 e i π 4 sgn(q) e i R (2π) n n/2 det Q /2 2 ξt Q ξ Fa(ξ)dξ. Using the Taylor s expansion formula for the exponential function, we see that for any non-negative integer N, the difference is bounded by e i 2 ξt Q ξ Fa(ξ)dξ N k=0 k! ( i 2 )k (ξ T Q ξ) k Fa(ξ)dξ N+ 2 N+ (N + )! (ξ T Q ξ) k Fa(ξ) dξ.
11 LECTURE 5: THE METHOD OF STATIONARY PHASE Note that so by the Fourier inversion formula, The conclusion follows. (ξ T Q ξ) k Fa(ξ) = F(p Q (D) k a)(ξ), p Q (D) k a(0) = (2π) n (ξ T Q ξ) k Fa(ξ)dξ. The stationary phase formula for general phase. As we have seen, only the critical points of the phase function ϕ give an essential contribution to the oscillatory integral I = e i ϕ(x) a(x)dx. In what follows we will assume that ϕ C ( ) admits only non-degenerate critical points in a neighborhood of the support of a. Since non-degenerate critical points must be discrete, ϕ has only finitely many critical points in the support of a. Thus one can find a partition of unity {U i i N + } of the support of a so that each U i, i N, contains exactly one critical point p i of ϕ, and U N+ contains no critical points of ϕ. A partition of unity arguments converts the asymptotic behavior of I to the sum of the asymptotic behavior of I p i = e i ϕ(x) U i χi (x)a(x)dx, where χ i is a function compactly supported in U i and equals identically near p i. According to the Morse lemma, for each U i (one can always shrink U i in the above integral if necessary) one can find a diffeomorphism f : V i U i so that f(0) = p i and f ϕ(x) = ϕ(p i ) + 2 (x2 + + x 2 r i x 2 r i + x 2 n) =: ψ pi (x), where r i is the number of positive eigenvalues of the matrix d 2 ϕ(p i ). It follows from the change of variable formula that = a(f(x))χ i (f(x)) det df(x) dx. I p i e i ψp i (x) V i Moreover, since ψ pi (x) has a unique non-degenerate critical point at 0, modulo O( ) one has (3) I p i = a(f(x)) det df(x) dx + O( ). e i ψp i (x) The asymptotic of this integral is basically given by theorem (2.8). In particular, to get the leading term in the asymptotic expansion above, one only need to evaluate det df(0), which follows from
12 2 LECTURE 5: THE METHOD OF STATIONARY PHASE Lemma 2.9. Let ψ, ϕ be smooth functions defined on V and U respectively, and suppose 0 is a non-degenerate critical point of both ψ and p a non-degenerate critical point of ϕ. If f : V U a diffeomorphism so that f(0) = p and f ϕ = ψ. Then (32) df T (0)d 2 ϕ(p)df(0) = d 2 ψ(0). Proof. We start with the equation ϕ f = ψ. Taking derivatives of both sides one gets ϕ (f(x)) f k = ψ. y k x j x j Taking the second order derivatives at 0 of both sides, and using the conditions f(0) = p and ϕ(p) = 0, one gets in other words It follows that in our case, 2 ϕ y k y l (p) f l x i (0) f k x j (0) = 2 ψ x i x j (0), df T (0)d 2 ϕ(p)df(0) = d 2 ψ(0). (33) det df(0) = det d 2 ϕ(p i ) /2. In conclusion, we have Theorem 2.0. As 0, the oscillatory integral has the asymptotic expansion (34) I I p i, where dϕ(p i )=0 (35) I p i (2π)n/2 e i ϕ(p i) e iπ 4 sgn(d2 ϕ(p i )) a(p i ) j L det d 2 ϕ(p i ) /2 j (a)(p i ), where L j = L j (x, D) is a differential operator of order 2j (which depends on the phase function ϕ) with L 0 =. In particular, the leading term is (36) I = (2π) n/2 e i ϕ(p i ) e iπ 4 sgn(d2 ϕ(p i )) a(p i ) det d 2 ϕ(p i ) + O( n /2 2 + ). dϕ(p i )=0 Remark. There are more complicated version of the stationary phase method where the phase function is allowed to have degenerate critical points. j 0
MATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More 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 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 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 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 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 informationIntroduction to Pseudodifferential Operators
Introduction to Pseudodifferential Operators Mengxuan Yang Directed by Prof. Dean Baskin May, 206. Introductions and Motivations Classical Mechanics & Quantum Mechanics In classical mechanics, the status
More informationIntroduction to Microlocal Analysis
Introduction to Microlocal Analysis First lecture: Basics Dorothea Bahns (Göttingen) Third Summer School on Dynamical Approaches in Spectral Geometry Microlocal Methods in Global Analysis University of
More informationMATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS. Beforehand we constructed distributions by taking the set Cc
MATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS Beforehand we constructed distributions by taking the set Cc ( ) as the set of very nice functions, and defined distributions as continuous linear
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 informationChapter 4. Inverse Function Theorem. 4.1 The Inverse Function Theorem
Chapter 4 Inverse Function Theorem d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d dd d d d d This chapter
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 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 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 informationINTRODUCTION TO DISTRIBUTIONS
INTRODUCTION TO DISTRIBUTIONS TSOGTGEREL GANTUMUR Abstract. We introduce locally convex spaces by the seminorms approach, and present the fundamentals of distributions. Contents 1. Introduction 1 2. Locally
More informationLECTURE 10: THE ATIYAH-GUILLEMIN-STERNBERG CONVEXITY THEOREM
LECTURE 10: THE ATIYAH-GUILLEMIN-STERNBERG CONVEXITY THEOREM Contents 1. The Atiyah-Guillemin-Sternberg Convexity Theorem 1 2. Proof of the Atiyah-Guillemin-Sternberg Convexity theorem 3 3. Morse theory
More information2. Function spaces and approximation
2.1 2. Function spaces and approximation 2.1. The space of test functions. Notation and prerequisites are collected in Appendix A. Let Ω be an open subset of R n. The space C0 (Ω), consisting of the C
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 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 informationDEVELOPMENT OF MORSE THEORY
DEVELOPMENT OF MORSE THEORY MATTHEW STEED Abstract. In this paper, we develop Morse theory, which allows us to determine topological information about manifolds using certain real-valued functions defined
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 information13. Fourier transforms
(December 16, 2017) 13. Fourier transforms Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2017-18/13 Fourier transforms.pdf]
More informationCoordinate Systems and Canonical Forms
Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
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 informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationNormal form for the non linear Schrödinger equation
Normal form for the non linear Schrödinger equation joint work with Claudio Procesi and Nguyen Bich Van Universita di Roma La Sapienza S. Etienne de Tinee 4-9 Feb. 2013 Nonlinear Schrödinger equation Consider
More informationNotes for Elliptic operators
Notes for 18.117 Elliptic operators 1 Differential operators on R n Let U be an open subset of R n and let D k be the differential operator, 1 1 x k. For every multi-index, α = α 1,...,α n, we define A
More informationNotes on Distributions
Notes on Distributions Functional Analysis 1 Locally Convex Spaces Definition 1. A vector space (over R or C) is said to be a topological vector space (TVS) if it is a Hausdorff topological space and the
More informationFOURIER TRANSFORMS OF SURFACE MEASURE ON THE SPHERE MATH 565, FALL 2017
FOURIER TRANSFORMS OF SURFACE MEASURE ON THE SPHERE MATH 565, FALL 17 1. Fourier transform of surface measure on the sphere Recall that the distribution u on R n defined as u, ψ = ψ(x) dσ(x) S n 1 is compactly
More informationChapter III Beyond L 2 : Fourier transform of distributions
Chapter III Beyond L 2 : Fourier transform of distributions 113 1 Basic definitions and first examples In this section we generalize the theory of the Fourier transform developed in Section 1 to distributions.
More informationTOPICS IN FOURIER ANALYSIS-IV. Contents
TOPICS IN FOURIER ANALYSIS-IV M.T. NAIR Contents 1. Test functions and distributions 2 2. Convolution revisited 8 3. Proof of uniqueness theorem 13 4. A characterization of distributions 14 5. Distributions
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationLocally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem
56 Chapter 7 Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem Recall that C(X) is not a normed linear space when X is not compact. On the other hand we could use semi
More information1 Fourier Integrals of finite measures.
18.103 Fall 2013 1 Fourier Integrals of finite measures. Denote the space of finite, positive, measures on by M + () = {µ : µ is a positive measure on ; µ() < } Proposition 1 For µ M + (), we define the
More information2 Infinite products and existence of compactly supported φ
415 Wavelets 1 Infinite products and existence of compactly supported φ Infinite products.1 Infinite products a n need to be defined via limits. But we do not simply say that a n = lim a n N whenever the
More informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More informationMath 6455 Oct 10, Differential Geometry I Fall 2006, Georgia Tech. Integration on Manifolds, Volume, and Partitions of Unity
Math 6455 Oct 10, 2006 1 Differential Geometry I Fall 2006, Georgia Tech Lecture Notes 13 Integration on Manifolds, Volume, and Partitions of Unity Suppose that we have an orientable Riemannian manifold
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 informationHomework I, Solutions
Homework I, Solutions I: (15 points) Exercise on lower semi-continuity: Let X be a normed space and f : X R be a function. We say that f is lower semi - continuous at x 0 if for every ε > 0 there exists
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 information10. Smooth Varieties. 82 Andreas Gathmann
82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationMath 172 Problem Set 8 Solutions
Math 72 Problem Set 8 Solutions Problem. (i We have (Fχ [ a,a] (ξ = χ [ a,a] e ixξ dx = a a e ixξ dx = iξ (e iax e iax = 2 sin aξ. ξ (ii We have (Fχ [, e ax (ξ = e ax e ixξ dx = e x(a+iξ dx = a + iξ where
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationFourier transforms, I
(November 28, 2016) Fourier transforms, I Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/fun/notes 2016-17/Fourier transforms I.pdf]
More information1. Bounded linear maps. A linear map T : E F of real Banach
DIFFERENTIABLE MAPS 1. Bounded linear maps. A linear map T : E F of real Banach spaces E, F is bounded if M > 0 so that for all v E: T v M v. If v r T v C for some positive constants r, C, then T is bounded:
More informationGENERALIZED FUNCTIONS LECTURES. Contents
GENERALIZED FUNCTIONS LECTURES Contents 1. The space of generalized functions on R n 2 1.1. Motivation 2 1.2. Basic definitions 3 1.3. Derivatives of generalized functions 5 1.4. The support of generalized
More informationLECTURE 5: COMPLEX AND KÄHLER MANIFOLDS
LECTURE 5: COMPLEX AND KÄHLER MANIFOLDS Contents 1. Almost complex manifolds 1. Complex manifolds 5 3. Kähler manifolds 9 4. Dolbeault cohomology 11 1. Almost complex manifolds Almost complex structures.
More informationTransformations from R m to R n.
Transformations from R m to R n 1 Differentiablity First of all because of an unfortunate combination of traditions (the fact that we read from left to right and the way we define matrix multiplication
More informationfy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))
1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical
More informationHermite functions. Jordan Bell Department of Mathematics, University of Toronto. September 9, 2015
Hermite functions Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto September 9, 25 Locally convex spaces If V is a vector space and {p α : α A} is a separating family
More informationLECTURE 5: SMOOTH MAPS. 1. Smooth Maps
LECTURE 5: SMOOTH MAPS 1. Smooth Maps Recall that a smooth function on a smooth manifold M is a function f : M R so that for any chart 1 {ϕ α, U α, V α } of M, the function f ϕ 1 α is a smooth function
More informationINTEGRATION OF FORMS. 3.1 Introduction. This is page 105 Printer: Opaque this CHAPTER 3
This is page 105 Printer: Opaque this CHAPTER 3 INTEGRATION OF FORMS 3.1 Introduction The change of variables formula asserts that if U and V are open subsets of R n and f : U V a C 1 diffeomorphism then,
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 informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationLECTURE 3 Functional spaces on manifolds
LECTURE 3 Functional spaces on manifolds The aim of this section is to introduce Sobolev spaces on manifolds (or on vector bundles over manifolds). These will be the Banach spaces of sections we were after
More informationINVERSE FUNCTION THEOREM and SURFACES IN R n
INVERSE FUNCTION THEOREM and SURFACES IN R n Let f C k (U; R n ), with U R n open. Assume df(a) GL(R n ), where a U. The Inverse Function Theorem says there is an open neighborhood V U of a in R n so that
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationMA651 Topology. Lecture 9. Compactness 2.
MA651 Topology. Lecture 9. Compactness 2. This text is based on the following books: Topology by James Dugundgji Fundamental concepts of topology by Peter O Neil Elements of Mathematics: General Topology
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationHandlebody Decomposition of a Manifold
Handlebody Decomposition of a Manifold Mahuya Datta Statistics and Mathematics Unit Indian Statistical Institute, Kolkata mahuya@isical.ac.in January 12, 2012 contents Introduction What is a handlebody
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationVector fields Lecture 2
Vector fields Lecture 2 Let U be an open subset of R n and v a vector field on U. We ll say that v is complete if, for every p U, there exists an integral curve, γ : R U with γ(0) = p, i.e., for every
More informationConvergence Rates on Root Finding
Convergence Rates on Root Finding Com S 477/577 Oct 5, 004 A sequence x i R converges to ξ if for each ǫ > 0, there exists an integer Nǫ) such that x l ξ > ǫ for all l Nǫ). The Cauchy convergence criterion
More informationINTRODUCTION TO NONLINEAR WAVE EQUATIONS
INTRODUCTION TO NONLINEAR WAVE EQUATIONS JONATHAN LUK Contents 1. Introduction 1. Linear wave equation via Fourier transform 4.1. Review of Fourier theory 4.. Returning to the linear wave equation 4 3.
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More information7.1 Tangent Planes; Differentials of Maps Between
Chapter 7 Tangent Planes Reading: Do Carmo sections 2.4 and 3.2 Today I am discussing 1. Differentials of maps between surfaces 2. Geometry of Gauss map 7.1 Tangent Planes; Differentials of Maps Between
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
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 informationLecture 5 - Hausdorff and Gromov-Hausdorff Distance
Lecture 5 - Hausdorff and Gromov-Hausdorff Distance August 1, 2011 1 Definition and Basic Properties Given a metric space X, the set of closed sets of X supports a metric, the Hausdorff metric. If A is
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 informationSolutions to Problem Set 5 for , Fall 2007
Solutions to Problem Set 5 for 18.101, Fall 2007 1 Exercise 1 Solution For the counterexample, let us consider M = (0, + ) and let us take V = on M. x Let W be the vector field on M that is identically
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationThe wave front set of a distribution
The wave front set of a distribution The Fourier transform of a smooth compactly supported function u(x) decays faster than any negative power of the dual variable ξ; that is for every number N there exists
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 informationThe theory of manifolds Lecture 2
The theory of manifolds Lecture 2 Let X be a subset of R N, Y a subset of R n and f : X Y a continuous map. We recall Definition 1. f is a C map if for every p X, there exists a neighborhood, U p, of p
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
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 informationTips and Tricks in Real Analysis
Tips and Tricks in Real Analysis Nate Eldredge August 3, 2008 This is a list of tricks and standard approaches that are often helpful when solving qual-type problems in real analysis. Approximate. There
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 informationLagrangian submanifolds and generating functions
Chapter 4 Lagrangian submanifolds and generating functions Motivated by theorem 3.9 we will now study properties of the manifold Λ φ X (R n \{0}) for a clean phase function φ. As shown in section 3.3 Λ
More informationWAVELET EXPANSIONS OF DISTRIBUTIONS
WAVELET EXPANSIONS OF DISTRIBUTIONS JASSON VINDAS Abstract. These are lecture notes of a talk at the School of Mathematics of the National University of Costa Rica. The aim is to present a wavelet expansion
More informationSmooth Structure. lies on the boundary, then it is determined up to the identifications it 1 2
132 3. Smooth Structure lies on the boundary, then it is determined up to the identifications 1 2 + it 1 2 + it on the vertical boundary and z 1/z on the circular part. Notice that since z z + 1 and z
More informationHW1 solutions. 1. α Ef(x) β, where Ef(x) is the expected value of f(x), i.e., Ef(x) = n. i=1 p if(a i ). (The function f : R R is given.
HW1 solutions Exercise 1 (Some sets of probability distributions.) Let x be a real-valued random variable with Prob(x = a i ) = p i, i = 1,..., n, where a 1 < a 2 < < a n. Of course p R n lies in the standard
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 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 informationThe Poisson summation formula, the sampling theorem, and Dirac combs
The Poisson summation ormula, the sampling theorem, and Dirac combs Jordan Bell jordanbell@gmailcom Department o Mathematics, University o Toronto April 3, 24 Poisson summation ormula et S be the set o
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
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 informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
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 informationProblem: A class of dynamical systems characterized by a fast divergence of the orbits. A paradigmatic example: the Arnold cat.
À È Ê ÇÄÁ Ë ËÌ ÅË Problem: A class of dynamical systems characterized by a fast divergence of the orbits A paradigmatic example: the Arnold cat. The closure of a homoclinic orbit. The shadowing lemma.
More informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationDifferentiable Functions
Differentiable Functions Let S R n be open and let f : R n R. We recall that, for x o = (x o 1, x o,, x o n S the partial derivative of f at the point x o with respect to the component x j is defined as
More information(d). Why does this imply that there is no bounded extension operator E : W 1,1 (U) W 1,1 (R n )? Proof. 2 k 1. a k 1 a k
Exercise For k 0,,... let k be the rectangle in the plane (2 k, 0 + ((0, (0, and for k, 2,... let [3, 2 k ] (0, ε k. Thus is a passage connecting the room k to the room k. Let ( k0 k (. We assume ε k
More informationChapter 8. P-adic numbers. 8.1 Absolute values
Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.
More information