High Dimensional Integration: New Weapons Fighting the Curse of Dimensionality
|
|
- Hilary Robbins
- 6 years ago
- Views:
Transcription
1 High Dimensional Integration: ew Weapons Fighting the Curse of Dimensionality Peter Zinterhof Technical Report 26-4 October 26 Department of Computer Sciences Jakob-Haringer-Straße 2 52 Salzburg Austria Technical Report Series
2 HIGH DIMESIOAL ITEGRATIO: EW WEAPOS FIGHTIG THE CURSE OF DIMESIOALITY P.Zinterhof Department of Computer Sciences, University of Salzburg DEDICATED TO THE MEMORY OF.M. KOROBOV ITRODUCTIO The approximate computation of the definite integral of a function of several variables is one of the basic problems of numerical analysis. The problem is hard because of the so called curse of dimensionality. This curse consists of the following: applying an integration rule with nodes to an univariate function, you will get an integration error, say, ε >. Applying the corresponding cartesian product rule to a s-variate function, you will need s nodes for the same integration error ε >. In mechanics we deal with at least six-dimensional functions, but in contemporary financial mathematics occur 3-variate functions. The probabilistic Monte-Carlo-Methods provide error estimates independent of the dimensionality of the problem. Unfortunately these methods are both slow in convergence end suffer of a lack of effectiveness as well. The Quasi-Monte-Carlo-Methods, based on number theory, are working fast and effectively, at least in the case of finite and smooth integrands. Unfortunately, in reality multivariate functions with singularities do occur. The scope of the present paper is numerical integration of multivariate functions with singularities. In many cases the proposed methods are best possible with respect to the order of convergence. Best possible means an exact order of the error term, essentially not worse than in the univariate case.. THE PROBLEM SETTIG Consider functions f(x, x 2,..., x s ) = y, x ρ, ρ =,..., s. Let I s = (, ) s the open unit cube, and Īs = [, ] s the closed s-dimensional unit cube. We are concerned with the numerical approximation of the integral of the function f by means of finite sums. Given a finite set of points in I s or Īs, (x (), x(2) we consider the integration method R =,..., x(s) ),..., (x() n, x (2) n f(x () n,..., x(s) n )...,..., x (s) n ),..., (x (),..., x(s) ), f(x,..., x s )dx... dx s () One is interested in small values of R, of course. Some known results: If the pointed (x () n,..., x (s) n ), n =,...,, is a set of uniform distributed and independent random variables, one receives the domical estimation of Monte-Carlo-Integration:
3 2 ( ) R = O (2) This convergence rate is rather poor, but independent of the dimensionality of the problem and independent of the smoothness of the function f(x,..., x s ). othing is said about the constants involved. On the other hand we consider the Cartesian Product Rules: Let x, x 2,..., x Īs and y = f(x) be a continuous function on Ī = [, ]. So we have an one-dimensional integration rule R () = f(x n ) f(x)dx (3) The principle of the cartesian product rules consists of a repeated application of the onedimensional rule to a s-variate function: R (s) s = s n = f(x n,..., x ns ) n s= f(x,..., x s )dx... dx s (4) The error term R (s) will not be better than s R(), in general. But the computational complexity is s. This fact is the well known Curse of Dimensionality. There are to remedies: The inequality of Hlawka-Koksma and Korobovs method: Let x, x 2,..., x Īs. Let I( a) = x : x ρ a ρ, ρ =,..., s, a Īs. Definition: D := sup a #{x n I( a)} a a 2... a s (5) is called the *-discrepancy (star discrepancy) of the finite point set x,..., x. The following theorem is essentially due to H. Weyl: Weyl s Criterion: The infinite sequence ( x n ), x n Īs, is uniform distributed if one of the following conditions holds: (a) for all continuous functions f : Īs holds lim f( x n ) = f( x)d x (6) Ī s (b) for all m s holds lim e 2πi m xn = { m m = (7)
4 3 (c) lim D = Weyl s criterion is the guideline for numerical application of number theoretical methods. At first we cite the Hlawka-Koksma-Inequality: Theorem: (Hlawka): Let f( x) be a function with bounded variation in the sense of Hardy-Krause, V (f( x)) <. Then holds the inequality R (f) = f( x) f( x)d x Ī D V (f) (8) s There is a huge number of estimations of the discrepancy of special sequences. We give only two examples. Example : Let x n,...,n s = ( n, n2,..., ) ns, n,..., n s =,...,. Then D s 2s. This also means the Curse of Dimensionality. Example 2: Let (a,..., a s ) s be optimal coefficients modulo in the sense of Korobov. Let x n = ( na,..., ) nas mod n, n =,...,. Then ( ) D = O (ln ) β, β s, (9) holds. Apart of the logarithmic factor this estimation is independent of the dimensionality of the problem. Unfortunately, the Hlawka-Koksma-inequality does not take into account additional smoothness conditions of the function f( x). Korobovs method overcomes this flaw: Let m = max(, m ), m. Consider the Korobov classes } Es {f( x) α (C) = C : C( m) ( m,..., m s ) α, m s, () where C( m) means the Fouriercoefficients of f( x): C( m) = f( x)e 2πi m x d x () Ī s Remark: If f( x) is -periodic in each variable x,..., x s, and if αs f α x,..., α x is continuous and bounded by C, then f E α s (C). This can be shown by αs-fold partial integrations of formular ().
5 4 Theorem: (Korobov): If f( x) Es α(c)and if a = (a,..., a s ) consists of optimal coefficients in the sense of Korobov, then the estimation holds: R (f) = f( n a ) f( x)d x Ī C C(ln αβ ) α, (2) s with an explicit constant C and some β s. This estimation is best possible apart from logarithmic factors: There is always a function f( x) Es α (C), such that f( n a ) f( x)d x Ī C(s)lns () α. (3) s More general, there is no integration rule with R = o ( α ), if f E α s. All this methods are classical and can be found in Korobov [3], Drmota-Tichy [] or iederreiter [2]. The methods described are concerned only with proper integrals of bounded functions. Singularities are not allowed. From the theoretical and also from the practical point of view it is important to develop integration rules for unbounded functions as well: Problem: Let f( x): Ī s or I s. Find classes of unbounded functions f and integration rules g n,f( x n ), such that lim g n, f( x n ) = f( x)d x. (4) I s Furthermore, give estimations for the error term R = g n, f( x n ) f( x)d x. (5) I s 2. SOLUTIOS OF THE PROBLEM We distinguish the two cases: First case: The location of the singularities of f( x) in Īs is unknown. Second case: The location of the singularities is known. We assume, that f( x) is unbounded at most on the bonudary of I s = (, ) s. For the sake of completeness we refer some of own former results [4]. Given a function f : Īs, so we define functions f B, ˆf B, B > such that
6 5 f B ( x) = f( x), if f( x) B =, if f( x) > B ˆf B ( x) =, if f( x) B = f(x), if f( x) > B (6) (7) So we have f( x) = f B ( x) + ˆf B ( x). We gave a suitable class of functions in the following manner: Definition: The class C(β, γ) of s-variate functions f( x), x, consists of all functions which fulfill B > : (a) I( ˆf B ) = O(B β ) for some β > (8) (b) V (f B ) = O(B γ ) for some γ (9) Here V (.) means again the variation of a function in the sense of Hardy and Krause. For dimension s = the definition coincide with the usual total variation of an univariate function. The use of V (.) is natural, because of the functional analytic connection between the spaces of continuous functions and the spaces of Radon measures, i.e. point measures and Lebesgue-measure. We proved the following theorem: Theorem: If f( x) C(β, γ) and if the discrepancy of the set of nodes x, x 2,..., x is D, then for B = (D ) (β+γ) the estimation holds: I(f) = f B ( x n ) + O ( (D) β (β+γ) ) (2) Remark: We also proved, that the order of convergence stated in (2) is best possible even in the case s =, provided f( x) C(β, γ). ow we come to the case two, the new and much more efficient results concerning the case, that the singularities of the integrand are concentrated on the boundary I s of the unit cube. The idea of the method: Consider an univariate function f(x), f : (, ), which has singularities at x = or x =, and which fulfills some smoothness conditions in (, ). We ask
7 6 for an integral-preserving transformation of f(x) which also continuates the differentiability conditions of f(x) to I = [, ]. Let p(t) = x be a function, which is strictly increasing in [, ] and which fulfills differentiability conditions of sufficient high order. Then we have for functions p(t) with p() =, p() = : f(x)dx = f(p(t))p (t)dt = g(t)dt (2) If p(t) does not tend to fast to p() = and p() =, then one will be able to remove singularities at x =, by means of (2). We propose the function p(t) = p γ (t) = p t (τ( τ))γ dτ, p = ( (τ( τ))γ dτ) (22) The connection of p(t) with the incomplete Beta-Integral is clear. We state some important properties of p(t): Lemma: (a) p() =, p() = (b) p () = p () =, p(t) > for t (, ) (c) p (n) () = p (n) () = for n =, 2,..., n < γ (d) p (n) (t) p γ (t( t)) γ+ n for n < γ + and t (e) p γ p i+2j=n n! i!j! (f) p(t)( p(t)) 4+2γ+ (γ+) p (t( t)) γ+ (g) p(t) p γ+ tγ+, p(t) p γ+ ( t)γ+ for t Some proofs of the parts of the Lemma are straight forward, some are not. We now introduce a suitable class of functions, having singularities on I s : Definition: Hs β,α (C) consists of all functions f(x,..., x s ), < x ρ <, ρ =,..., s, such that for all n,..., n s, n ρ α, ρ =,..., s, holds: n+...+ns f(x,..., x s ) C x n xn xns s ( s ρ= (x ρ( x ρ )) β+nρ ) (23)
8 7 whereas all the derivatives are continuous, and < β <. The introduction of the class Hs β,α (C) was motivated by the univariate ëxtremefunction f(x) = (x( x)) β, < β <. We remind (2) for general s =, 2,... : with... f(x, x 2,..., x s )dx dx 2... dx s =... g(t,..., t s )dt dt s (24) g(t,..., t s ) = f(p(t ), p(t 2 ),..., p(t s ))p (t )p (t 2 )... p (t s ) (25) We consider now the reactors of nodes T n = ( 2 + na, 2 + na2,..., 2 + ) nas, mod, where a = (a,..., a s ) are optimal coefficients, a = a(), and n =,...,. We get the integration rule I (f) := f(p(t,n ), p(t 2,n )..., p(t s,n ))p (t,n ), p (t 2,n ),..., p (t s,n ), (26) with t ρ,n = 2 + naρ, ρ =,..., s. ow we are able to state the Theorem: If f Hs β,α (C), and if γ > α+β β, then... f(x,..., x s )dx... dx s I (f) C (α, β, γ, s)c (ln )α,β α, (27) where the constant C (α, β, γ, s) is explicit. The proof makes heavy use of the lemma and makes use of an explicit and complicated estimation of all of the derivatives of g(t,..., t s ). Remark : According to (3), our theorem can not be improved significantly, even in the case of boundedness of f( x). Remark 2: The use of the classical optimal coefficients is only one example of the application of number theoretical methods to improper integrals. We have further methods, using e.g. the Weyl-sequences, (n Θ), especially the sequences n(e r, e r2,..., e rs ), n =, 2,..., r i r k, i k. Estimations of R = fdx I (f) via the Diaphony are available as well.
9 REFERECES 8 LITERATUR Korobovs book is the classical reference, whereas iederreiters book contains most of the resent developments in number theoretical numerics. Drmota-Tichy is the perhaps comprehensive book on uniform distribution of sequences, containing two thousands of references. REFERECES [] M.-Drmota, R.-F.-Tichy, Sequences, Discrepancies and Applications. Lecture otes in Mathematics 65. Springer (997). [2] H.-iederreiter, Random umber Generation and Quasi-Monte Carlo Methods. CBMS- SF REGIOAL COFERECE SERIES I APPLIED MATHEMAICS 63. Capital City Press (992). [3].-M.- Korobov, umber-theoretic Methods in Approximate Analysis, (963). [4] P.-Zinterhof, High Dimension Improper Integration.Proceedings of 7th Intern. Scientific Conference-Section, Applied Mathematics, pg. 9-5, Kosice, Slovakia (22).
c Springer Verlag. The copyright for this contribution is held by Springer Verlag. The original publication is available at
c Springer Verlag. The copyright for this contribution is held by Springer Verlag. The original publication is available at www.springerlink.com. A Pragmatic view on umerical Integration of Unbounded Functions
More informationDyadic diaphony of digital sequences
Dyadic diaphony of digital sequences Friedrich Pillichshammer Abstract The dyadic diaphony is a quantitative measure for the irregularity of distribution of a sequence in the unit cube In this paper we
More informationNumerical Methods in Economics MIT Press, Chapter 9 Notes Quasi-Monte Carlo Methods. Kenneth L. Judd Hoover Institution.
1 Numerical Methods in Economics MIT Press, 1998 Chapter 9 Notes Quasi-Monte Carlo Methods Kenneth L. Judd Hoover Institution October 26, 2002 2 Quasi-Monte Carlo Methods Observation: MC uses random sequences
More informationf(x n ) [0,1[ s f(x) dx.
ACTA ARITHMETICA LXXX.2 (1997) Dyadic diaphony by Peter Hellekalek and Hannes Leeb (Salzburg) 1. Introduction. Diaphony (see Zinterhof [13] and Kuipers and Niederreiter [6, Exercise 5.27, p. 162]) is a
More informationarxiv: v1 [math.co] 8 Feb 2013
ormal numbers and normality measure Christoph Aistleitner arxiv:302.99v [math.co] 8 Feb 203 Abstract The normality measure has been introduced by Mauduit and Sárközy in order to describe the pseudorandomness
More informationNEW SEQUENCES FROM OLD: A SEARCH FOR LOW DISCREPANCY SEQUENCES. Jarvis Haupt
EW SEQUECES FROM OLD: A SEARCH FOR LOW DISCREPACY SEQUECES Jarvis Haupt Electrical and Computer Engineering Department University of Wisconsin, Madison, WI, USA jdhaupt@wisc.edu ABSTRACT Aside from their
More informationQuasi- Monte Carlo Multiple Integration
Chapter 6 Quasi- Monte Carlo Multiple Integration Introduction In some sense, this chapter fits within Chapter 4 on variance reduction; in some sense it is stratification run wild. Quasi-Monte Carlo methods
More informationImproved Discrepancy Bounds for Hybrid Sequences. Harald Niederreiter. RICAM Linz and University of Salzburg
Improved Discrepancy Bounds for Hybrid Sequences Harald Niederreiter RICAM Linz and University of Salzburg MC vs. QMC methods Hybrid sequences The basic sequences Deterministic discrepancy bounds The proof
More informationChapter 8 Indeterminate Forms and Improper Integrals Math Class Notes
Chapter 8 Indeterminate Forms and Improper Integrals Math 1220-004 Class Notes Section 8.1: Indeterminate Forms of Type 0 0 Fact: The it of quotient is equal to the quotient of the its. (book page 68)
More informationuniform distribution theory
Uniform Distribution Theory 2 (2007), no.1, 1 10 uniform distribution theory UNIFORM DISTRIBUTION OF SEQUENCES CONNECTED WITH THE WEIGHTED SUM-OF-DIGITS FUNCTION Friedrich Pillichshammer ABSTRACT. In this
More informationA LeVeque-type lower bound for discrepancy
reprinted from Monte Carlo and Quasi-Monte Carlo Methods 998, H. Niederreiter and J. Spanier, eds., Springer-Verlag, 000, pp. 448-458. A LeVeque-type lower bound for discrepancy Francis Edward Su Department
More informationOn the existence of higher order polynomial. lattices with large figure of merit ϱ α.
On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationDIGITAL EXPANSION OF EXPONENTIAL SEQUENCES
DIGITAL EXPASIO OF EXPOETIAL SEQUECES MICHAEL FUCHS Abstract. We consider the q-ary digital expansion of the first terms of an exponential sequence a n. Using a result due to Kiss und Tichy [8], we prove
More informationQMC methods in quantitative finance. and perspectives
, tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation
More informationarxiv: v5 [math.ca] 23 Jan 2016
ariv: 454756 Calculation of Lebesgue Integrals by Using Uniformly Distributed Sequences ariv:507.02978v5 [math.ca] 23 Jan 206 Gogi Pantsulaia and Tengiz Kiria I.Vekua Institute of Applied Mathematics,
More informationUniform distribution of sequences connected with the weighted sum-of-digits function
Uniform distribution of seuences connected with the weighted sum-of-digits function Friedrich Pillichshammer Abstract In this paper we consider seuences which are connected with the so-called weighted
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationAbstract. 2. We construct several transcendental numbers.
Abstract. We prove Liouville s Theorem for the order of approximation by rationals of real algebraic numbers. 2. We construct several transcendental numbers. 3. We define Poissonian Behaviour, and study
More informationAPPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.
APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,
More informationQuasi-Monte Carlo Algorithms with Applications in Numerical Analysis and Finance
Quasi-Monte Carlo Algorithms with Applications in Numerical Analysis and Finance Reinhold F. Kainhofer Rigorosumsvortrag, Institut für Mathematik Graz, 16. Mai 2003 Inhalt 1. Quasi-Monte Carlo Methoden
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationDigital lattice rules
Digital lattice rules Multivariate integration and discrepancy estimates A thesis presented to The University of New South Wales in fulfilment of the thesis requirement for the degree of Doctor of Philosophy
More informationRandom and Quasi-Random Point Sets
Peter Hellekalek Gerhard Larcher (Editors) Random and Quasi-Random Point Sets With contributions by JözsefBeck Peter Hellekalek Fred J. Hickemell Gerhard Larcher Pierre L'Ecuyer Harald Niederreiter Shu
More informationOn the inverse of the discrepancy for infinite dimensional infinite sequences
On the inverse of the discrepancy for infinite dimensional infinite sequences Christoph Aistleitner Abstract In 2001 Heinrich, ovak, Wasilkowski and Woźniakowski proved the upper bound s,ε) c abs sε 2
More informationepub WU Institutional Repository
epub WU Institutional Repository Wolfgang Hörmann and Josef Leydold Quasi Importance Sampling Working Paper Original Citation: Hörmann, Wolfgang and Leydold, Josef (2005) Quasi Importance Sampling. Preprint
More informationLow Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed?
Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Xiaoqun Wang 1,2 and Ian H. Sloan 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084,
More informationA NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1)-SEQUENCES. Peter Kritzer 1. Abstract
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 5(3 (2005, #A11 A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1-SEQUENCES Peter Kritzer 1 Department of Mathematics, University of Salzburg,
More information1 Which sets have volume 0?
Math 540 Spring 0 Notes #0 More on integration Which sets have volume 0? The theorem at the end of the last section makes this an important question. (Measure theory would supersede it, however.) Theorem
More information1. Introduction. In this paper we consider stochastic optimization problems of the form
O RATES OF COVERGECE FOR STOCHASTIC OPTIMIZATIO PROBLEMS UDER O-I.I.D. SAMPLIG TITO HOMEM-DE-MELLO Abstract. In this paper we discuss the issue of solving stochastic optimization problems by means of sample
More informationPart 2 Continuous functions and their properties
Part 2 Continuous functions and their properties 2.1 Definition Definition A function f is continuous at a R if, and only if, that is lim f (x) = f (a), x a ε > 0, δ > 0, x, x a < δ f (x) f (a) < ε. Notice
More informationIntroduction to Discrepancy Theory
Introduction to Discrepancy Theory Aleksandar (Sasho) Nikolov University of Toronto Sasho Nikolov (U of T) Discrepancy 1 / 23 Outline 1 Monte Carlo 2 Discrepancy and Quasi MC 3 Combinatorial Discrepancy
More informationuniform distribution theory
Uniform Distribution Theory 1 (26), no.1, 45 63 uniform distribution theory Oto Strauch Oľga Blažeková ABSTRACT. In this paper we show that the sequence p n/n mod 1, n = 1, 2,..., where p n is the nth
More informationThe Discrepancy Function and the Small Ball Inequality in Higher Dimensions
The Discrepancy Function and the Small Ball Inequality in Higher Dimensions Dmitriy Georgia Institute of Technology (joint work with M. Lacey and A. Vagharshakyan) 2007 Fall AMS Western Section Meeting
More informationTutorial on quasi-monte Carlo methods
Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo
More informationNORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY
ORMAL UMBERS AD UIFORM DISTRIBUTIO WEEKS -3 OPE PROBLEMS I UMBER THEORY SPRIG 28, TEL AVIV UIVERSITY Contents.. ormal numbers.2. Un-natural examples 2.3. ormality and uniform distribution 2.4. Weyl s criterion
More information11.8 Power Series. Recall the geometric series. (1) x n = 1+x+x 2 + +x n +
11.8 1 11.8 Power Series Recall the geometric series (1) x n 1+x+x 2 + +x n + n As we saw in section 11.2, the series (1) diverges if the common ratio x > 1 and converges if x < 1. In fact, for all x (
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationLECTURE 22: INTEGRATION ON MANIFOLDS. 1. Orientations
LECTURE 22: INTEGRATION ON ANIFOLDS 1. Orientations Let be a smooth manifold of dimension n, and let ω Ω n () be a smooth n-form. We want to define the integral ω. First assume = R n. In calculus we learned
More informationThe Halton sequence and its discrepancy in the Cantor expansion
Period Math Hung 2017 75:128 141 DOI 10.1007/s10998-016-0169-5 The Halton sequence and its discrepancy in the Cantor expansion Alena Haddley 1 Poj Lertchoosakul 2 Radhakrishnan air 1 Published online:
More informationCounting patterns in digital expansions
Counting patterns in digital expansions Manfred G. Madritsch Department for Analysis and Computational Number Theory Graz University of Technology madritsch@math.tugraz.at Séminaire ERNEST Dynamique, Arithmétique,
More informationf(s) e -i n π s/l d s
Pointwise convergence of complex Fourier series Let f(x) be a periodic function with period l defined on the interval [,l]. The complex Fourier coefficients of f( x) are This leads to a Fourier series
More informationTractability of Multivariate Problems
Erich Novak University of Jena Chemnitz, Summer School 2010 1 Plan for the talk Example: approximation of C -functions What is tractability? Tractability by smoothness? Tractability by sparsity, finite
More informationConvergence of greedy approximation I. General systems
STUDIA MATHEMATICA 159 (1) (2003) Convergence of greedy approximation I. General systems by S. V. Konyagin (Moscow) and V. N. Temlyakov (Columbia, SC) Abstract. We consider convergence of thresholding
More informationA NOTE ON ALMOST PERIODIC VARIATIONAL EQUATIONS
A NOTE ON ALMOST PERIODIC VARIATIONAL EQUATIONS PETER GIESL AND MARTIN RASMUSSEN Abstract. The variational equation of a nonautonomous differential equation ẋ = F t, x) along a solution µ is given by ẋ
More informationHyperuniformity on the Sphere
Hyperuniformity on the Sphere Johann S. Brauchart j.brauchart@tugraz.at 14. June 2018 Week 2 Workshop MATRIX Research Program "On the Frontiers of High Dimensional Computation" MATRIX, Creswick [ 1 / 56
More informationNon-linear wave equations. Hans Ringström. Department of Mathematics, KTH, Stockholm, Sweden
Non-linear wave equations Hans Ringström Department of Mathematics, KTH, 144 Stockholm, Sweden Contents Chapter 1. Introduction 5 Chapter 2. Local existence and uniqueness for ODE:s 9 1. Background material
More informationLow-discrepancy constructions in the triangle
Low-discrepancy constructions in the triangle Department of Statistics Stanford University Joint work with Art B. Owen. ICERM 29th October 2014 Overview Introduction 1 Introduction The Problem Motivation
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
More informationProgress in high-dimensional numerical integration and its application to stochastic optimization
Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,
More informationUniformly distributed sequences of partitions
Uniformly distributed sequences of partitions Aljoša Volčič Università della Calabria Uniform distribution and QMC methods A. Volčič (Università della Calabria) U.d. sequences of partitions Linz, October
More informationMathematics for Economists
Mathematics for Economists Victor Filipe Sao Paulo School of Economics FGV Metric Spaces: Basic Definitions Victor Filipe (EESP/FGV) Mathematics for Economists Jan.-Feb. 2017 1 / 34 Definitions and Examples
More informationA List of Problems in Real Analysis
A List of Problems in Real Analysis W.Yessen & T.Ma December 3, 218 This document was first created by Will Yessen, who was a graduate student at UCI. Timmy Ma, who was also a graduate student at UCI,
More informationLS-sequences of points in the unit square
LS-sequences of points in the unit square arxiv:2.294v [math.t] 3 ov 202 Ingrid Carbone, Maria Rita Iacò, Aljoša Volčič Abstract We define a countable family of sequences of points in the unit square:
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More informationConstructing Guaranteed Automatic Numerical Algorithms for U
Constructing Guaranteed Automatic Numerical Algorithms for Univariate Integration Department of Applied Mathematics, Illinois Institute of Technology July 10, 2014 Contents Introduction.. GAIL What do
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationThe Hilbert Transform and Fine Continuity
Irish Math. Soc. Bulletin 58 (2006), 8 9 8 The Hilbert Transform and Fine Continuity J. B. TWOMEY Abstract. It is shown that the Hilbert transform of a function having bounded variation in a finite interval
More informationVerona Course April Lecture 1. Review of probability
Verona Course April 215. Lecture 1. Review of probability Viorel Barbu Al.I. Cuza University of Iaşi and the Romanian Academy A probability space is a triple (Ω, F, P) where Ω is an abstract set, F is
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 informationTHE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES
THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES MICHAEL DRMOTA Abstract. The purpose of this paper is to provide upper bounds for the discrepancy of generalized Van-der-Corput-Halton sequences
More informationStrictly Positive Definite Functions on the Circle
the Circle Princeton University Missouri State REU 2012 Definition A continuous function f : [0, π] R is said to be positive definite on S 1 if, for every N N and every set of N points x 1,..., x N on
More informationIntegral representations in models with long memory
Integral representations in models with long memory Georgiy Shevchenko, Yuliya Mishura, Esko Valkeila, Lauri Viitasaari, Taras Shalaiko Taras Shevchenko National University of Kyiv 29 September 215, Ulm
More informationPhenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 2012
Phenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 202 BOUNDS AND ASYMPTOTICS FOR FISHER INFORMATION IN THE CENTRAL LIMIT THEOREM
More informationPartial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces
Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As
More informationPhysics 411 Lecture 7. Tensors. Lecture 7. Physics 411 Classical Mechanics II
Physics 411 Lecture 7 Tensors Lecture 7 Physics 411 Classical Mechanics II September 12th 2007 In Electrodynamics, the implicit law governing the motion of particles is F α = m ẍ α. This is also true,
More informationScenario generation in stochastic programming
Scenario generation in stochastic programming Werner Römisch Humboldt-University Berlin, Institute of Mathematics, 10099 Berlin, Germany, romisch@math.hu-berlin.de Abstract Stability-based methods for
More informationMonte Carlo Method for Numerical Integration based on Sobol s Sequences
Monte Carlo Method for Numerical Integration based on Sobol s Sequences Ivan Dimov and Rayna Georgieva Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, 1113 Sofia,
More informationNote that in the example in Lecture 1, the state Home is recurrent (and even absorbing), but all other states are transient. f ii (n) f ii = n=1 < +
Random Walks: WEEK 2 Recurrence and transience Consider the event {X n = i for some n > 0} by which we mean {X = i}or{x 2 = i,x i}or{x 3 = i,x 2 i,x i},. Definition.. A state i S is recurrent if P(X n
More informationCSCI-6971 Lecture Notes: Monte Carlo integration
CSCI-6971 Lecture otes: Monte Carlo integration Kristopher R. Beevers Department of Computer Science Rensselaer Polytechnic Institute beevek@cs.rpi.edu February 21, 2006 1 Overview Consider the following
More informationChapter 8. Infinite Series
8.4 Series of Nonnegative Terms Chapter 8. Infinite Series 8.4 Series of Nonnegative Terms Note. Given a series we have two questions:. Does the series converge? 2. If it converges, what is its sum? Corollary
More informationImportance Sampling to Accelerate the Convergence of Quasi-Monte Carlo
Importance Sampling to Accelerate the Convergence of Quasi-Monte Carlo Wolfgang Hörmann, Josef Leydold Department of Statistics and Mathematics Wirtschaftsuniversität Wien Research Report Series Report
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationSINGULAR INTEGRALS ON SIERPINSKI GASKETS
SINGULAR INTEGRALS ON SIERPINSKI GASKETS VASILIS CHOUSIONIS Abstract. We construct a class of singular integral operators associated with homogeneous Calderón-Zygmund standard kernels on d-dimensional,
More informationPseudorandom Sequences II: Exponential Sums and Uniform Distribution
Pseudorandom Sequences II: Exponential Sums and Uniform Distribution Arne Winterhof Austrian Academy of Sciences Johann Radon Institute for Computational and Applied Mathematics Linz Carleton University
More informationL p -convergence of Bernstein Kantorovich-type operators. by Michele Campiti (Bari) and Giorgio Metafune (Lecce)
ANNALES POLONICI MATHEMATICI LXIII.3 (996) L p -convergence of Bernstein Kantorovich-type operators by Michele Campiti (Bari) and Giorgio Metafune (Lecce) Abstract. We study a Kantorovich-type modification
More information2 2 + x =
Lecture 30: Power series A Power Series is a series of the form c n = c 0 + c 1 x + c x + c 3 x 3 +... where x is a variable, the c n s are constants called the coefficients of the series. n = 1 + x +
More informationApproximation numbers of Sobolev embeddings - Sharp constants and tractability
Approximation numbers of Sobolev embeddings - Sharp constants and tractability Thomas Kühn Universität Leipzig, Germany Workshop Uniform Distribution Theory and Applications Oberwolfach, 29 September -
More informationReceived: 2/7/07, Revised: 5/25/07, Accepted: 6/25/07, Published: 7/20/07 Abstract
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 2007, #A34 AN INVERSE OF THE FAÀ DI BRUNO FORMULA Gottlieb Pirsic 1 Johann Radon Institute of Computational and Applied Mathematics RICAM,
More informationSJÄLVSTÄNDIGA ARBETEN I MATEMATIK
SJÄLVSTÄNDIGA ARBETEN I MATEMATIK MATEMATISKA INSTITUTIONEN, STOCKHOLMS UNIVERSITET Gamma function related to Pic functions av Saad Abed 25 - No 4 MATEMATISKA INSTITUTIONEN, STOCKHOLMS UNIVERSITET, 6 9
More informationInfinite Series Summary
Infinite Series Summary () Special series to remember: Geometric series ar n Here a is the first term and r is the common ratio. When r
More informationlim when the limit on the right exists, the improper integral is said to converge to that limit.
hapter 7 Applications of esidues - evaluation of definite and improper integrals occurring in real analysis and applied math - finding inverse Laplace transform by the methods of summing residues. 6. Evaluation
More informationOn the law of the iterated logarithm for the discrepancy of lacunary sequences
On the law of the iterated logarithm for the discrepancy of lacunary sequences Christoph Aistleitner Abstract A classical result of Philipp (1975) states that for any sequence (n k ) k 1 of integers satisfying
More informationIEOR 6711: Stochastic Models I Fall 2013, Professor Whitt Lecture Notes, Thursday, September 5 Modes of Convergence
IEOR 6711: Stochastic Models I Fall 2013, Professor Whitt Lecture Notes, Thursday, September 5 Modes of Convergence 1 Overview We started by stating the two principal laws of large numbers: the strong
More informationBoth these computations follow immediately (and trivially) from the definitions. Finally, observe that if f L (R n ) then we have that.
Lecture : One Parameter Maximal Functions and Covering Lemmas In this first lecture we start studying one of the basic and fundamental operators in harmonic analysis, the Hardy-Littlewood maximal function.
More informationMATH 409 Advanced Calculus I Lecture 7: Monotone sequences. The Bolzano-Weierstrass theorem.
MATH 409 Advanced Calculus I Lecture 7: Monotone sequences. The Bolzano-Weierstrass theorem. Limit of a sequence Definition. Sequence {x n } of real numbers is said to converge to a real number a if for
More informationPenalty and Barrier Methods General classical constrained minimization problem minimize f(x) subject to g(x) 0 h(x) =0 Penalty methods are motivated by the desire to use unconstrained optimization techniques
More informationChapter 1. Introduction to prime number theory. 1.1 The Prime Number Theorem
Chapter 1 Introduction to prime number theory 1.1 The Prime Number Theorem In the first part of this course, we focus on the theory of prime numbers. We use the following notation: we write f( g( as if
More informationA CENTRAL LIMIT THEOREM FOR NESTED OR SLICED LATIN HYPERCUBE DESIGNS
Statistica Sinica 26 (2016), 1117-1128 doi:http://dx.doi.org/10.5705/ss.202015.0240 A CENTRAL LIMIT THEOREM FOR NESTED OR SLICED LATIN HYPERCUBE DESIGNS Xu He and Peter Z. G. Qian Chinese Academy of Sciences
More informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
More informationAnalytic families of multilinear operators
Analytic families of multilinear operators Mieczysław Mastyło Adam Mickiewicz University in Poznań Nonlinar Functional Analysis Valencia 17-20 October 2017 Based on a joint work with Loukas Grafakos M.
More informationDIMENSIONS OF JULIA SETS OF HYPERBOLIC ENTIRE FUNCTIONS
Bull. London Math. Soc. 36 2004 263 270 C 2004 London Mathematical Society DOI: 10.1112/S0024609303002698 DIMENSIONS OF JULIA SETS OF HYPERBOLIC ENTIRE FUNCTIONS GWYNETH M. STALLARD Abstract It is known
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationEE1 and ISE1 Communications I
EE1 and ISE1 Communications I Pier Luigi Dragotti Lecture two Lecture Aims To introduce signals, Classifications of signals, Some particular signals. 1 Signals A signal is a set of information or data.
More informationLower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions
Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationFractional Derivatives of Weierstrass-Type Functions
Fractional Derivatives of Weierstrass-Type Functions M. Zähle H. Ziezold Abstract For a special type of fractional differentiation formulas for the Weierstrass and the Weierstrass-Mandelbrot functions
More information