A complex exponential Fourier transform approach to gradient density estimation
|
|
- George Byrd
- 5 years ago
- Views:
Transcription
1 A complex onential Fourier transform approach to gradient density estimation Karthik S. Gurumoorthy a,, Anand Rangarajan a, Arunava Banerjee a a Department of Computer and Information Science and Engineering, University of Florida, Gainesville, Florida, USA Abstract We prove a novel result wherein the density function of the gradients of a thrice differentiable function S obtained via a random variable transformation of a uniformly distributed random variable defined on a closed, bounded interval Ω R is accurately approximated by the normalized power spectrum of ϕ = is as the free parameter 0. The result is shown using the well known stationary phase approximation and standard integration techniques and requires proper ordering of limits. Experimental results provide anecdotal visual evidence corroborating the result.. Introduction The literature is replete with techniques which attempt to estimate a non-observable probability density function using observed data believed to be sampled from an underlying distribution. Density estimation techniques have a long history and run the gamut of histogramming, Parzen windows [see Parzen 962], vector quantization, wavelets etc. In the present work, we prove a novel mathematical result relating the density function of the gradients of a function S [viewed as a random variable transformation Y = S X where X is uniformly distributed] with the normalized power spectrum of is as the quantity 0. Corresponding Author Address: E30, CSE Building, University of Florida, P.O. Box 620, Gainesville, FL , USA. Ph: , Fax: addresses: ksg@cise.ufl.edu Karthik S. Gurumoorthy, anand@cise.ufl.edu Anand Rangarajan, arunava@cise.ufl.edu Arunava Banerjee Preprint submitted to Statistics and Probability Letters March, 20
2 We provide a simple Fourier transform-based technique to determine the density of Y. Our method computes the gradient density directly from the function S, circumventing the need to compute its derivative S. The approach is based on ressing the function S as the phase of a wave function ϕ, specifically ϕx = for isx small values of and then considering the normalized power spectrum magnitude squared of the Fourier transform of ϕ [see Bracewell 999]. Using the stationary phase approximation a well known technique in asymptotic analysis [please see Olver 974a] we show that in the limiting case as 0, the power spectrum of ϕ converges to the density of Y and hence can serve as its density estimator at small, non-zero values of. Our new mathematical relationship is motivated by the classical-quantum relation, wherein classical physics is ressed as a limiting case of quantum mechanics [please refer to Griffiths 2004; Feynman and Hibbs 965]. When S is treated as the Hamilton-Jacobi scalar field, the gradient of S corresponds to the classical momentum of a particle [see Goldstein et al. 200]. In the parlance of quantum mechanics, the magnitude square of the wave function ressed either in its position or momentum basis corresponds to its position or momentum density respectively. Since these representations either in the position or momentum basis are simply the suitably scaled Fourier transforms of each other, the magnitude square of the Fourier transform of the wave function ressed in its position basis, is its quantum momentum density [see Griffiths 2004]. The principal theorem proved in the article [Theorem 2.5] states that the classical momentum density denoted by P can be ressed as a limiting case as 0 of its corresponding quantum momentum density denoted by P, in complete agreement with the correspondence principle [see Griffiths 2004]. 2. Equivalence of the gradient density and the power spectrum As stated above, the density function for the gradients of S denoted by Y can be obtained via a random variable transformation of a uniformly distributed random variable X using the derivative S as the transformation function, namely, Y = S X. We assume that S is thrice differentiable on a closed, bounded interval Ω = [b, b 2 ] with length L = b 2 b and has a non-vanishing second derivative almost everywhere 2
3 on Ω, i.e. µ {x : S x = 0} = 0, 2. where µ denotes the Lebesgue measure. The assumption in 2. is made in order to ensure that the density function of Y exists almost everywhere. This is further clarified in Lemma 2.2 below. Define the following sets: B {x : S x = 0}, C {S x : x B} {S b, S b 2 }, and A u {x : S x = u}. 2.2 Here, S b = lim x b + S x and S b 2 = lim x b S x. The higher derivatives 2 of S at the end points b, b 2 are also defined along similar lines using one-sided limits. The main purpose of defining these one-sided limits is to exactly determine the set C where the density of Y is not defined. Since µb = 0, we also have µc = 0. Lemma 2.. [Finiteness Lemma] A u is finite for every u / C. Lemma 2.2. [Density Lemma] The probability density of Y on R C exists and is given by P = L Nu 0 k= S x k, 2.3 where the summation is over A u0 which is the finite set of locations x k Ω where S x k = as per Lemma 2., with A u0 = N. Since the density is based on the transformation Y = S X, the probability density function in 2.3 assumes the existence of the inverse of the transformation function S [see Billingsley 995]. This is made licit in the following lemma which is required by the main theorem. Lemma 2.3. [Interval Lemma] For every u / C, η > 0 and a closed interval J η = [u η, u + η] such that J η C is empty. The proofs of the lemmas are available in AppendixA. We now prove the main result which relates the normalized power spectrum of in the limit as 0 with the probability density of the random variable isx Y = S X denoted by P. 3
4 Define a function F : R R + C as ˆ b2 isx F u, 2πL For a fixed value of, define a function F : R C as b iux dx. 2.4 F u F u,. 2.5 Observe that F is closely related to the Fourier transform of isx. The scale factor 2πL is the normalizing term such that the L 2 norm of F is one as seen in the following lemma whose proof is straightforward and omitted for lack of space. Lemma 2.4. With F defined as above, F L 2 R and F =. Define a function P : R C R + as By definition, P P u F u 2 = F uf u Since µc = 0, from Lemma 2.4, P udu =. Hence, treating P u as a density function, we have the following theorem statement. Theorem 2.5. If P and P are defined as above, then the ˆ u0 +α lim α 0 α lim P udu = P, / C Before embarking on the proof, we would like to emphasize that the ordering of the limits and the integral as given in the theorem statement is crucial and cannot be arbitrarily interchanged. To press this point home, we show below that after solving for P, the lim 0 P does not exist. Hence, the order of the integral followed by the limit 0 cannot be interchanged. Furthermore, when we swap the limits between α and, we get lim lim 0 α 0 α ˆ u0 +α P udu = lim 0 P 2.8 which does not exist. Hence, the theorem statement can be valid only for the specified sequence of limits and the integral. 4
5 2.. Brief osition of the result To understand the result in simpler terms, let us reconsider the definition of the scaled Fourier transform given in 2.4. The first onential is is a varying complex sinusoid, whereas the second onential iux is a fixed complex sinusoid at frequency u. When we multiply these two complex onentials, at low values of, the two sinusoids are usually not in sync and tend to cancel each other out. However, around the locations where S x = u, the two sinusoids are in perfect sync as the combined onent is stationary with the approximate duration of this resonance depending on S x. The value of the integral in 2.4 can be approximated via the stationary phase approximation [please see Olver 974a] as F u ±i π L 4 Nu k= { } i Sx k ux k S x k 2.9 where Nu = A u. The approximation is increasingly tight as 0. The squared Fourier transform P gives us the required result L Nu k= S x k except for the cross phase factors Sx k Sx l ux k x l obtained as a byproduct of two or more remote locations x k and x l indexing into the same frequency bin u, i.e, x k x l, but S x k = S x l = u. Integrating the squared Fourier transform over a small frequency range [u, u + α] removes these cross phase factors and we obtain the desired result Formal proof We shall now provide the proof by considering different cases. case i: Let us consider the case in which no stationary points exist for the given, i.e, there is no x Ω for which S x =. Let tx = Sx x. Then, t x is of constant sign in [b, b 2 ] and hence tx is strictly monotonic. Defining v = tx, we have from 2.4, ˆ tb2 iv F = gvdv πL tb 5
6 Here, gv = t x where x = t v. Integrating by parts, we get F 2πL = i i [ ˆ tb2 tb itb2 g tb 2 Then F 2πL S b 2 + S b + itb ] g tb iv g vdv. 2. ˆ tb2 tb g v dv. 2.2 Hence, F γ, where γ > 0 is a continuous function of. Then P γ 2. Since S x is continuous and / C, we can find a ρ > 0 such that for every u [ ρ, + ρ], no stationary points exist. The value ρ can also be chosen appropriately such that [ ρ, + ρ] C =. If α < ρ, then ˆ u0 +α lim P udu = Furthermore, from 2.3 we have P = 0. The result immediately follows. case ii: We now consider the case where stationary points exist. Since we are interested only in the situation as 0, the stationary phase method in Olver 974a,b can be used to obtain a good approximation for F defined in 2.4 and 2.5. The phase term in this function, Sx x, is stationary only when S x =. Consider the set A u0 defined in 2.2. Since it is finite by Lemma 2., let A u0 = {x, x 2,..., x Nu0} with x k < x k+, k. We break Ω into disjoint intervals such that each interval has utmost one stationary point. To this end, consider numbers {c, c 2,..., c Nu0 +} such that b < c < x, x k tx = Sx x. Then, < c k+ < x k+ and x Nu0 < c Nu0+ < b 2. Let F Nu 0 2πL = G + G 2 + K k + K k 2.4 k= 6
7 where G = G 2 = K k = K k = ˆ c b ˆ b2 itx dx, itx itx itx c Nu0 + ˆ xk c ˆ k ck+ x k dx, dx, and dx. 2.5 Note that the integrals G and G 2 do not have any stationary points. From case i above, we get G + G 2 = ϵ, = O 2.6 as 0. Furthermore, ϵ, γ 2, where γ 2 > 0 is a continuous function of. In order to evaluate K k and K k, observe that when we and the phase term up to second order, tx tx k Qx k x x k 2 as x x k, where Qx k = S x k 2. Furthermore, in the open intervals c k, x k and x k, c k+, t x = S x is continuous and is of constant sign. From Theorem 3. in Olver 974a, we get K k = K k = 2 ± iπ 4 Γ 2 itxk 2 S x k + ϵ 2,. 2.7 From Lemma 2.3 in Olver 974a, it can be verified that ϵ 2, = o as 0 and can also be uniformly bounded by a function of independent of for small values of. In 2.7, Γ is the Gamma function and the sign in the phase term depends on whether S x k > 0 or S x k < 0. Plugging the values of these integrals in 2.4 and noting that Γ 2 = π, we get F 2πL = N k= i [Sx 2π k x k ] S x k ± iπ 4 +ϵ, + ϵ 2,. 2.8 Hence, F = N L k= i [Sx k x k ] ± iπ S x k ϵ 3,,
8 where ϵ 3, = ϵ,+ϵ 2, 2πL = o as 0. From the definition of P u in 2.6, we have P = L + L N k= Nu 0 k= S x k Nu 0 l=;l k cos [Sx k Sx l x k x l ] + θx k, x l S x k 2 S x l 2 +ϵ 4,, 2.2 where ϵ 4, includes both the magnitude square of ϵ 3, and the cross terms involving the main first term in 2.9 and ϵ 3,. Notice that the main term in 2.9 can be bounded by a function of independent of as i [Sx k x k ] =, 2.22 and S x k 0, k. Since ϵ 3, = o, we get ϵ 4, = o as 0. Additionally, θx k, x l = 0, π 2 or π 2 and θx l, x k = θx k, x l. The first term in 2.2 exactly matches the ression for P as seen from Lemma 2.2. But, since lim 0 cos is not defined, lim 0 P does not exist and hence the cross cosine terms do not vanish when we take the limit. We now show that integrating P u over a small non-zero interval [, + α] and then taking the limit with respect to followed by the limit with respect to α does yield the density of Y. From Lemmas 2. and 2.3, we see that for a given a A u0, when u is varied over the interval J η = [ η, + η], the inverse function S u is well defined with S u N a, where N a is some small neighborhood around a. For each a A u0, define the inverse function S a u : J η N a as S a u = x iff u = S x and x N a Unfortunately, when we move from a fixed value to a variable u defined in the interval J η, the locations x k and x l which previously were fixed in 2.2 now also vary over the interval J η. This makes the notation somewhat cumbersome and unwieldy. Using the inverse functions defined in 2.23 and defining a k x k [and consequently 8
9 a l = x l ], we define x k u S a k u [and therefore xl u = S a l u] for u J η. Finally, define the functions p u and q u over J η as Observe that p u S x k u S x l u u x k u x l u, and 2.24 q u. S x k u 2 S x l u p u = S x k ux ku S x l ux lu x k u x l u u x ku x lu = x l u x k u 2.26 as u = S x k u = S x l u. In particular, if x l > x k, then x l u > x k u and vice versa. Hence, p u never vanishes and is also of constant sign over J η. Then, it follows that p u is strictly monotonic and specifically bijective on J η. We will use this result in the subsequent steps. Now, let α < η. Since the additional error term ϵ 4, in 2.2 converges to zero as 0 and can also be uniformly bounded by a function of for small values of, we have Then, we get where I L I 2 L Here, I 3 k, l is given by I 3 k, l ˆ u0+α ˆ u0+α lim 0 ϵ 4, du = ˆ u0 +α lim 0 P udu = I + I N k= N k= ˆ u0 +α S du, and 2.29 x k u N l=;l k lim I 3k, l [ ] p u q u cos + θ x k u, x l u du. 2.3 When α < η, the sign of S x k u around x k and the sign of S x l u around x l do not change over the interval [, + α]. Since θx k u, x l u depends on the sign of S, θ is constant on [, +α] and equals θ = θx k, x l. 9
10 Now, anding the cosine term in 2.3, we get ˆ u0+α p u I 3 k, l = cosθ q u cos ˆ u0 +α p u sinθ q u sin Since p u is bijective, we get via a standard change of variables: du du I 3 k, l = cosθ I 4 k, l sinθ I 5 k, l 2.33 where I 4 k, l = I 5 k, l = ˆ β2 β ˆ β2 β v cos g vdv, and 2.34 v sin g vdv Here, β = p, β 2 Integrating I 4 k, l by parts, we get β 2 I 4 k, l = sin g Then, I 4 k, l g = p + α and g v = q p v p p v. ˆ β2 β β 2 β 2 sin β g β v sin g vdv g β ˆ β2 + g v dv β It is worth mentioning that q and hence g are differentiable over their respective intervals as the sign of S x k us x l u does not change over the interval [, + α]. We then have I 4 k, l M where M is some constant independent of. Hence, lim 0 I 4 k, l = 0, k, l. By a similar argument, lim 0 I 5 k, l = 0, k, l. From 2.30 and 2.33, we get I 2 = 0. Since lim α 0 α I = Nu 0 L S x k = P, 2.38 k= the main result ressed in Theorem 2.5 follows. 0
11 Below, we show comparisons between our Fourier transform approach and the standard histogramming technique on some trigonometric and onential functions sampled on a regular grid between [ 0.25, 0.25] at a grid spacing of 2 5. For the sake of convenience, we normalized the functions such that its maximum gradient value is. Using the sampled values Ŝ, we computed the Fast Fourier transform of iŝ at = , took its magnitude squared and then normalized it to compute the gradient density. We also computed the discrete derivative of S at the grid locations and then determined its gradient density by histogramming using around 220 bins. The plots shown in Figure provide anecdotal empirical evidence supporting the mathematical result stated in Theorem 2.5. Notice the near-perfect match between the gradient densities computed via standard histogramming and the gradient densities determined using our Fourier transform method. Figure : Comparison results. i Left: Gradient densities obtained from histogramming, ii Right : Gradient densities obtained from squared Fourier transform of the wave function 3. Discussion Observe that the integrals I = ˆ u0+α ˆ u0+α P udu, I = P udu 3.
12 give the interval measures of the density functions P and P respectively. Theorem 2.5 states that at small values of, both the interval measures are approximately equal, with the difference between them being o. Recall that by definition, P ϕx = isx is the normalized power spectrum [see Bracewell 999] of the wave function. Hence, we conclude that the power spectrum of ϕx can potentially serve as a density estimator for the gradients of S at small values of. The erimental results shown above serve as a demonstration, anecdotally attesting to the verity of the result. Our result is directly inspired by the three-way relationships between the classical momentum S, the quantum momentum operator i x and its spatial frequency spectrum. Since these relationships hold in higher dimensions as well, we are likewise interested in extending our density estimation result to higher dimensions. This is a fruitful avenue for future work. References Billingsley, P., 995. Probability and measure, 3rd Edition. Wiley-Interscience. Bracewell, R., 999. The Fourier transform and its applications, 3rd Edition. McGraw- Hill Science and Engineering. Feynman, R., Hibbs, A., 965. Quantum mechanics and path integrals. McGraw-Hill. Goldstein, H., Poole, C., Safko, J., 200. Classical mechanics, 3rd Edition. Addison- Wesley. Griffiths, D., Introduction to quantum mechanics, 2nd Edition. Benjamin- Cummings. Olver, F., 974a. Asymptotics and special functions. Academic Press, New York and London. Olver, F., 974b. Error bounds for stationary phase approximations. SIAM Journal of Mathematical Analysis 5,
13 Parzen, E., 962. On the estimation of a probability density function and the mode. Annals of Mathematical Statistics 33, Rudin, W., 976. Principles of mathematical analysis, 3rd Edition. McGraw-Hill. AppendixA. Proof of Lemmas. Proof of Finiteness Lemma We prove the result by contradiction. Observe that A u is a subset of the compact set Ω. If A u is not finite, then by Theorem 2.37 in Rudin 976, A u has a limit point x 0 Ω. Consider a sequence {x n } n=, with each x n A u, converging to x 0. Since S x n = u, n, from the continuity of S we get S x 0 = u and hence x 0 A u. Then S x 0 S x n lim = 0 = S x 0. A. n x 0 x n Based on the definitions given in 2.2, we have x 0 B and hence u C resulting in a contradiction. 2. Proof of Interval Lemma Observe that B is closed because if x 0 is a limit point of B, from the continuity of S we have S x 0 = 0 and hence x 0 B. B is also compact as it is a closed subset of Ω. Since S is continuous, C = S B {S b, S b 2 } is also compact and hence R C is open. Then for u / C, there exists an open neighborhood N r u for some r > 0 around u such that N r u C =. By defining η = r 2, the proof is complete. 3. Proof of Density Lemma Since the random variable X is assumed to have a uniform distribution on Ω, its density is given by f X x = L for every x Ω. Recall that the random variable Y is obtained via a random variable transformation from X, using the function S. Hence, its density function exists on R C where we have banished the image under S of the measure zero set of points where S vanishes and is given by 2.3. The reader may refer to Billingsley 995 for a detailed lanation. 3
A Schrödinger formalism for simultaneously computing the Euclidean distance transform and its gradient density
A Schrödinger formalism for simultaneously computing the Euclidean distance transform and its gradient density ABSTRACT Karthik S. Gurumoorthy International Center for Theoretical Sciences, Tata Institute
More informationDISTANCE TRANSFORM GRADIENT DENSITY ESTIMATION USING THE STATIONARY PHASE APPROXIMATION
DISTANCE TRANSFORM GRADIENT DENSITY ESTIMATION USING THE STATIONARY PHASE APPROXIMATION KARTHIK S. GURUMOORTHY AND ANAND RANGARAJAN Abstract. The complex wave representation (CWR converts unsigned 2D distance
More informationFact: Every matrix transformation is a linear transformation, and vice versa.
Linear Transformations Definition: A transformation (or mapping) T is linear if: (i) T (u + v) = T (u) + T (v) for all u, v in the domain of T ; (ii) T (cu) = ct (u) for all scalars c and all u in the
More informationFrame Wavelet Sets in R d
Frame Wavelet Sets in R d X. DAI, Y. DIAO Department of Mathematics University of North Carolina at Charlotte Charlotte, NC 28223 xdai@uncc.edu Q. GU Department of Mathematics Each China Normal University
More informationProbability and Measure
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability
More informationDifferential of the Exponential Map
Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements
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 informationSolution. 1 Solution of Homework 7. Sangchul Lee. March 22, Problem 1.1
Solution Sangchul Lee March, 018 1 Solution of Homework 7 Problem 1.1 For a given k N, Consider two sequences (a n ) and (b n,k ) in R. Suppose that a n b n,k for all n,k N Show that limsup a n B k :=
More informationENERGY DECAY ESTIMATES FOR LIENARD S EQUATION WITH QUADRATIC VISCOUS FEEDBACK
Electronic Journal of Differential Equations, Vol. 00(00, No. 70, pp. 1 1. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp ENERGY DECAY ESTIMATES
More informationLecture Characterization of Infinitely Divisible Distributions
Lecture 10 1 Characterization of Infinitely Divisible Distributions We have shown that a distribution µ is infinitely divisible if and only if it is the weak limit of S n := X n,1 + + X n,n for a uniformly
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 informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationAnother Riesz Representation Theorem
Another Riesz Representation Theorem In these notes we prove (one version of) a theorem known as the Riesz Representation Theorem. Some people also call it the Riesz Markov Theorem. It expresses positive
More informationIntroduction to Real Analysis Alternative Chapter 1
Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces
More informationEnergy and Equations of Motion
Energy and Equations of Motion V. Tanrıverdi tanriverdivedat@googlemail.com Physics Department, Middle East Technical University, Ankara / TURKEY Abstract. From the total time derivative of energy an equation
More informationIMA Preprint Series # 2143
A BASIC INEQUALITY FOR THE STOKES OPERATOR RELATED TO THE NAVIER BOUNDARY CONDITION By Luan Thach Hoang IMA Preprint Series # 2143 ( November 2006 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationOnly Intervals Preserve the Invertibility of Arithmetic Operations
Only Intervals Preserve the Invertibility of Arithmetic Operations Olga Kosheleva 1 and Vladik Kreinovich 2 1 Department of Electrical and Computer Engineering 2 Department of Computer Science University
More informationFixed Points & Fatou Components
Definitions 1-3 are from [3]. Definition 1 - A sequence of functions {f n } n, f n : A B is said to diverge locally uniformly from B if for every compact K A A and K B B, there is an n 0 such that f n
More informationFilters in Analysis and Topology
Filters in Analysis and Topology David MacIver July 1, 2004 Abstract The study of filters is a very natural way to talk about convergence in an arbitrary topological space, and carries over nicely into
More information= F (b) F (a) F (x i ) F (x i+1 ). a x 0 x 1 x n b i
Real Analysis Problem 1. If F : R R is a monotone function, show that F T V ([a,b]) = F (b) F (a) for any interval [a, b], and that F has bounded variation on R if and only if it is bounded. Here F T V
More informationRearrangement on Conditionally Convergent Integrals in Analogy to Series
Electronic Proceedings of Undergraduate Mathematics Day, Vol. 3 (8), No. 6 Rearrangement on Conditionally Convergent Integrals in Analogy to Series Edward J Timko University of Dayton Dayton OH 45469-36
More informationEstimates for probabilities of independent events and infinite series
Estimates for probabilities of independent events and infinite series Jürgen Grahl and Shahar evo September 9, 06 arxiv:609.0894v [math.pr] 8 Sep 06 Abstract This paper deals with finite or infinite sequences
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationVectors. The standard geometric definition of vector is as something which has direction and magnitude but not position.
Vectors The standard geometric definition of vector is as something which has direction and magnitude but not position. Since vectors have no position we may place them wherever is convenient. Vectors
More information8 Singular Integral Operators and L p -Regularity Theory
8 Singular Integral Operators and L p -Regularity Theory 8. Motivation See hand-written notes! 8.2 Mikhlin Multiplier Theorem Recall that the Fourier transformation F and the inverse Fourier transformation
More informationColoring Vertices and Edges of a Path by Nonempty Subsets of a Set
Coloring Vertices and Edges of a Path by Nonempty Subsets of a Set P.N. Balister E. Győri R.H. Schelp April 28, 28 Abstract A graph G is strongly set colorable if V (G) E(G) can be assigned distinct nonempty
More informationMeasure Theory and Lebesgue Integration. Joshua H. Lifton
Measure Theory and Lebesgue Integration Joshua H. Lifton Originally published 31 March 1999 Revised 5 September 2004 bstract This paper originally came out of my 1999 Swarthmore College Mathematics Senior
More informationSMSTC (2017/18) Geometry and Topology 2.
SMSTC (2017/18) Geometry and Topology 2 Lecture 1: Differentiable Functions and Manifolds in R n Lecturer: Diletta Martinelli (Notes by Bernd Schroers) a wwwsmstcacuk 11 General remarks In this lecture
More informationVANDERBILT UNIVERSITY. MATH 3120 INTRO DO PDES The Schrödinger equation
VANDERBILT UNIVERSITY MATH 31 INTRO DO PDES The Schrödinger equation 1. Introduction Our goal is to investigate solutions to the Schrödinger equation, i Ψ t = Ψ + V Ψ, 1.1 µ where i is the imaginary number
More informationDerivation of the Nonlinear Schrödinger Equation. from First Principles. Theodore Bodurov
Annales de la Fondation Louis de Broglie, Volume 30, no 3-4, 2005 343 Derivation of the Nonlinear Schrödinger Equation from First Principles Theodore Bodurov Eugene, Oregon, USA, email: bodt@efn.org ABSTRACT.
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationTHE CANTOR GAME: WINNING STRATEGIES AND DETERMINACY. by arxiv: v1 [math.ca] 29 Jan 2017 MAGNUS D. LADUE
THE CANTOR GAME: WINNING STRATEGIES AND DETERMINACY by arxiv:170109087v1 [mathca] 9 Jan 017 MAGNUS D LADUE 0 Abstract In [1] Grossman Turett define the Cantor game In [] Matt Baker proves several results
More informationOptimality Conditions for Constrained Optimization
72 CHAPTER 7 Optimality Conditions for Constrained Optimization 1. First Order Conditions In this section we consider first order optimality conditions for the constrained problem P : minimize f 0 (x)
More informationStationary mean-field games Diogo A. Gomes
Stationary mean-field games Diogo A. Gomes We consider is the periodic stationary MFG, { ɛ u + Du 2 2 + V (x) = g(m) + H ɛ m div(mdu) = 0, (1) where the unknowns are u : T d R, m : T d R, with m 0 and
More informationEigenvalues and Eigenfunctions of the Laplacian
The Waterloo Mathematics Review 23 Eigenvalues and Eigenfunctions of the Laplacian Mihai Nica University of Waterloo mcnica@uwaterloo.ca Abstract: The problem of determining the eigenvalues and eigenvectors
More informationA note on some approximation theorems in measure theory
A note on some approximation theorems in measure theory S. Kesavan and M. T. Nair Department of Mathematics, Indian Institute of Technology, Madras, Chennai - 600 06 email: kesh@iitm.ac.in and mtnair@iitm.ac.in
More informationMATH 202B - Problem Set 5
MATH 202B - Problem Set 5 Walid Krichene (23265217) March 6, 2013 (5.1) Show that there exists a continuous function F : [0, 1] R which is monotonic on no interval of positive length. proof We know there
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationAn asymptotic ratio characterization of input-to-state stability
1 An asymptotic ratio characterization of input-to-state stability Daniel Liberzon and Hyungbo Shim Abstract For continuous-time nonlinear systems with inputs, we introduce the notion of an asymptotic
More informationGaussian vectors and central limit theorem
Gaussian vectors and central limit theorem Samy Tindel Purdue University Probability Theory 2 - MA 539 Samy T. Gaussian vectors & CLT Probability Theory 1 / 86 Outline 1 Real Gaussian random variables
More informationMATH31011/MATH41011/MATH61011: FOURIER ANALYSIS AND LEBESGUE INTEGRATION. Chapter 2: Countability and Cantor Sets
MATH31011/MATH41011/MATH61011: FOURIER ANALYSIS AND LEBESGUE INTEGRATION Chapter 2: Countability and Cantor Sets Countable and Uncountable Sets The concept of countability will be important in this course
More informationPart III. 10 Topological Space Basics. Topological Spaces
Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.
More informationASYMPTOTICALLY NONEXPANSIVE MAPPINGS IN MODULAR FUNCTION SPACES ABSTRACT
ASYMPTOTICALLY NONEXPANSIVE MAPPINGS IN MODULAR FUNCTION SPACES T. DOMINGUEZ-BENAVIDES, M.A. KHAMSI AND S. SAMADI ABSTRACT In this paper, we prove that if ρ is a convex, σ-finite modular function satisfying
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationNotes on Complex Analysis
Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................
More informationVolume in n Dimensions
Volume in n Dimensions MA 305 Kurt Bryan Introduction You ve seen that if we have two vectors v and w in two dimensions then the area spanned by these vectors can be computed as v w = v 1 w 2 v 2 w 1 (where
More informationLEBESGUE MEASURE AND L2 SPACE. Contents 1. Measure Spaces 1 2. Lebesgue Integration 2 3. L 2 Space 4 Acknowledgments 9 References 9
LBSGU MASUR AND L2 SPAC. ANNI WANG Abstract. This paper begins with an introduction to measure spaces and the Lebesgue theory of measure and integration. Several important theorems regarding the Lebesgue
More informationFIRST YEAR CALCULUS W W L CHEN
FIRST YER CLCULUS W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
More informationLECTURE 3 RANDOM VARIABLES, CUMULATIVE DISTRIBUTION FUNCTIONS (CDFs)
OCTOBER 6, 2014 LECTURE 3 RANDOM VARIABLES, CUMULATIVE DISTRIBUTION FUNCTIONS (CDFs) 1 Random Variables Random experiments typically require verbal descriptions, and arguments involving events are often
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationMA1131 Lecture 15 (2 & 3/12/2010) 77. dx dx v + udv dx. (uv) = v du dx dx + dx dx dx
MA3 Lecture 5 ( & 3//00) 77 0.3. Integration by parts If we integrate both sides of the proct rule we get d (uv) dx = dx or uv = d (uv) = dx dx v + udv dx v dx dx + v dx dx + u dv dx dx u dv dx dx This
More information2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is
More informationMATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5
MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have
More informationOn the Set of Limit Points of Normed Sums of Geometrically Weighted I.I.D. Bounded Random Variables
On the Set of Limit Points of Normed Sums of Geometrically Weighted I.I.D. Bounded Random Variables Deli Li 1, Yongcheng Qi, and Andrew Rosalsky 3 1 Department of Mathematical Sciences, Lakehead University,
More informationAN EFFECTIVE METRIC ON C(H, K) WITH NORMAL STRUCTURE. Mona Nabiei (Received 23 June, 2015)
NEW ZEALAND JOURNAL OF MATHEMATICS Volume 46 (2016), 53-64 AN EFFECTIVE METRIC ON C(H, K) WITH NORMAL STRUCTURE Mona Nabiei (Received 23 June, 2015) Abstract. This study first defines a new metric with
More informationReal and Complex Analysis, 3rd Edition, W.Rudin Elementary Hilbert Space Theory
Real and Complex Analysis, 3rd Edition, W.Rudin Chapter 4 Elementary Hilbert Space Theory Yung-Hsiang Huang. It s easy to see M (M ) and the latter is a closed subspace of H. Given F be a closed subspace
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationFeynman Path Integrals in Quantum Mechanics
Feynman Path Integrals in Quantum Mechanics Christian Egli October, 2004 Abstract This text is written as a report to the seminar course in theoretical physics at KTH, Stockholm. The idea of this work
More informationTensor Analysis in Euclidean Space
Tensor Analysis in Euclidean Space James Emery Edited: 8/5/2016 Contents 1 Classical Tensor Notation 2 2 Multilinear Functionals 4 3 Operations With Tensors 5 4 The Directional Derivative 5 5 Curvilinear
More informationLecture 3f Polar Form (pages )
Lecture 3f Polar Form (pages 399-402) In the previous lecture, we saw that we can visualize a complex number as a point in the complex plane. This turns out to be remarkable useful, but we need to think
More informationHausdorff Measure. Jimmy Briggs and Tim Tyree. December 3, 2016
Hausdorff Measure Jimmy Briggs and Tim Tyree December 3, 2016 1 1 Introduction In this report, we explore the the measurement of arbitrary subsets of the metric space (X, ρ), a topological space X along
More informationBernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of
More informationSPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS
SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS BRYAN P. RYNNE Abstract. We consider the m-point boundary value problem consisting of the equation u = f(u), on
More informationA Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets
A Minimal Uncertainty Product for One Dimensional Semiclassical Wave Packets George A. Hagedorn Happy 60 th birthday, Mr. Fritz! Abstract. Although real, normalized Gaussian wave packets minimize the product
More informationFourier Series. 1. Review of Linear Algebra
Fourier Series In this section we give a short introduction to Fourier Analysis. If you are interested in Fourier analysis and would like to know more detail, I highly recommend the following book: Fourier
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationIn English, this means that if we travel on a straight line between any two points in C, then we never leave C.
Convex sets In this section, we will be introduced to some of the mathematical fundamentals of convex sets. In order to motivate some of the definitions, we will look at the closest point problem from
More informationLECTURE 3: Quantization and QFT
LECTURE 3: Quantization and QFT Robert Oeckl IQG-FAU & CCM-UNAM IQG FAU Erlangen-Nürnberg 14 November 2013 Outline 1 Classical field theory 2 Schrödinger-Feynman quantization 3 Klein-Gordon Theory Classical
More informationMath212a1413 The Lebesgue integral.
Math212a1413 The Lebesgue integral. October 28, 2014 Simple functions. In what follows, (X, F, m) is a space with a σ-field of sets, and m a measure on F. The purpose of today s lecture is to develop the
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 225 Estimates of the second-order derivatives for solutions to the two-phase parabolic problem
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationAPPROXIMATING BOCHNER INTEGRALS BY RIEMANN SUMS
APPROXIMATING BOCHNER INTEGRALS BY RIEMANN SUMS J.M.A.M. VAN NEERVEN Abstract. Let µ be a tight Borel measure on a metric space, let X be a Banach space, and let f : (, µ) X be Bochner integrable. We show
More information01 Harmonic Oscillations
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-2014 01 Harmonic Oscillations Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More informationDIFFERENTIAL CRITERIA FOR POSITIVE DEFINITENESS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX)- DIFFERENTIAL CRITERIA FOR POSITIVE DEFINITENESS J. A. PALMER Abstract. We show how the Mellin transform can be used to
More information1 Orthonormal sets in Hilbert space
Math 857 Fall 15 1 Orthonormal sets in Hilbert space Let S H. We denote by [S] the span of S, i.e., the set of all linear combinations of elements from S. A set {u α : α A} is called orthonormal, if u
More information2 Sequences, Continuity, and Limits
2 Sequences, Continuity, and Limits In this chapter, we introduce the fundamental notions of continuity and limit of a real-valued function of two variables. As in ACICARA, the definitions as well as proofs
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 08: Stochastic Convergence
Introduction to Empirical Processes and Semiparametric Inference Lecture 08: Stochastic Convergence Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationClosest Moment Estimation under General Conditions
Closest Moment Estimation under General Conditions Chirok Han and Robert de Jong January 28, 2002 Abstract This paper considers Closest Moment (CM) estimation with a general distance function, and avoids
More informationEmpirical Processes: General Weak Convergence Theory
Empirical Processes: General Weak Convergence Theory Moulinath Banerjee May 18, 2010 1 Extended Weak Convergence The lack of measurability of the empirical process with respect to the sigma-field generated
More informationMATH 6605: SUMMARY LECTURE NOTES
MATH 6605: SUMMARY LECTURE NOTES These notes summarize the lectures on weak convergence of stochastic processes. If you see any typos, please let me know. 1. Construction of Stochastic rocesses A stochastic
More informationThe infinite square well in a reformulation of quantum mechanics without potential function
The infinite square well in a reformulation of quantum mechanics without potential function A.D. Alhaidari (a), T.J. Taiwo (b) (a) Saudi Center for Theoretical Physics, P.O Box 32741, Jeddah 21438, Saudi
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationColoring Vertices and Edges of a Path by Nonempty Subsets of a Set
Coloring Vertices and Edges of a Path by Nonempty Subsets of a Set P.N. Balister E. Győri R.H. Schelp November 8, 28 Abstract A graph G is strongly set colorable if V (G) E(G) can be assigned distinct
More informationComplex Analysis Qualifying Exam Solutions
Complex Analysis Qualifying Exam Solutions May, 04 Part.. Let log z be the principal branch of the logarithm defined on G = {z C z (, 0]}. Show that if t > 0, then the equation log z = t has exactly one
More informationOn the logarithm of the minimizing integrand for certain variational problems in two dimensions
On the logarithm of the minimizing integrand for certain variational problems in two dimensions University of Kentucky John L. Lewis University of Kentucky Andrew Vogel Syracuse University 2012 Spring
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationLinear Programming Methods
Chapter 11 Linear Programming Methods 1 In this chapter we consider the linear programming approach to dynamic programming. First, Bellman s equation can be reformulated as a linear program whose solution
More informationPattern generation, topology, and non-holonomic systems
Systems & Control Letters ( www.elsevier.com/locate/sysconle Pattern generation, topology, and non-holonomic systems Abdol-Reza Mansouri Division of Engineering and Applied Sciences, Harvard University,
More informationarxiv: v1 [cs.dm] 12 Jun 2016
A Simple Extension of Dirac s Theorem on Hamiltonicity Yasemin Büyükçolak a,, Didem Gözüpek b, Sibel Özkana, Mordechai Shalom c,d,1 a Department of Mathematics, Gebze Technical University, Kocaeli, Turkey
More informationConvexity in R n. The following lemma will be needed in a while. Lemma 1 Let x E, u R n. If τ I(x, u), τ 0, define. f(x + τu) f(x). τ.
Convexity in R n Let E be a convex subset of R n. A function f : E (, ] is convex iff f(tx + (1 t)y) (1 t)f(x) + tf(y) x, y E, t [0, 1]. A similar definition holds in any vector space. A topology is needed
More informationThe concentration of the chromatic number of random graphs
The concentration of the chromatic number of random graphs Noga Alon Michael Krivelevich Abstract We prove that for every constant δ > 0 the chromatic number of the random graph G(n, p) with p = n 1/2
More informationPart V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory
Part V 7 Introduction: What are measures and why measurable sets Lebesgue Integration Theory Definition 7. (Preliminary). A measure on a set is a function :2 [ ] such that. () = 2. If { } = is a finite
More informationSOME MEASURABILITY AND CONTINUITY PROPERTIES OF ARBITRARY REAL FUNCTIONS
LE MATEMATICHE Vol. LVII (2002) Fasc. I, pp. 6382 SOME MEASURABILITY AND CONTINUITY PROPERTIES OF ARBITRARY REAL FUNCTIONS VITTORINO PATA - ALFONSO VILLANI Given an arbitrary real function f, the set D
More informationFourier Theory on the Complex Plane II Weak Convergence, Classification and Factorization of Singularities
Fourier Theory on the Complex Plane II Weak Convergence, Classification and Factorization of Singularities Jorge L. delyra Department of Mathematical Physics Physics Institute University of São Paulo April
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 09: Stochastic Convergence, Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and
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 informationAsymptotic Behavior of Infinity Harmonic Functions Near an Isolated Singularity
Savin, O., and C. Wang. (2008) Asymptotic Behavior of Infinity Harmonic Functions, International Mathematics Research Notices, Vol. 2008, Article ID rnm163, 23 pages. doi:10.1093/imrn/rnm163 Asymptotic
More information