MATH 505b Project Random Matrices
|
|
- Adele Greer
- 5 years ago
- Views:
Transcription
1 MATH 505b Project andom Matrices Emre Demirkaya, Mark Duggan, Melike Sirlanci Tuysuzoglu April 30, Introduction Studies of andom Matrix Theory began in J. Wishart studied the distribution of eigenvalues of random matrices with fixed dimesion in his paper: The Generalized Product Moment Distribution in Samples from a Normal Multivariate Population. If the entries of the matrix are Gaussian (nonasymptotic), the joint distribution of the eigenvalues is known. In 1955, E. Wigner introduced asymptotic behavior of eigenvalues of random matrices, Characteristic vectors of bordered matrices with infinite dimensions. His motivation lies in nuclear physics: the statistics of experimentally measured atomic energy levels can be explained by the limiting spectrum of random matrices. That is, he studied the distribution of eigenvalues of random matrices as the matrix dimension tends to infinity. One interesting result, which can be regarded as the central limit theorem of random matrix theory, is that the asymptotic eigenvalue distribution of a random matrix does not depend on the distribution of the entries of the random matrix. In this project, we focus on a theorem about asymptotic behavior of normalized symmetric random matrices, which will be defined in the next chapter. As a further study, one may try to answer the natural question: how fast does the eigenvalue distribution converge as n increases? 1
2 2 Preliminaries In this section, we give the definition of Wigner matrices and other important definitions. Definition 2.1 (Wigner Matrix): A Wigner matrix is a symmetric n n matrix M n with entries Y ] Z i,j i<j M n (j, i) =M n (i, j) = [ Y i i = j. where {Z i,j} 1Æi<j, {Y i } 1Æi are two familes of independent identically distributed real-valued random variables with mean 0 and variance 1. emark 2.2: We call the matrix M n / Ô n a normalized Wigner matrix. We don t divide by Ô n on a whim, notice that E [Mn / Ô n] remains bounded as n æœ, while E [M n ] æœ,as n æœ. Example 2.3: When M i,j is normally distributed (i.e. Gaussian), we call this a Gaussian Wigner matrix. When Y i N(0, 2), and Z i,j is a standard normal, we call this the Gaussian Orthognal Ensemble, which earns its name because it is invariant under orthognal transformation. The Gaussian Orthognal Ensemble is useful, because it is highly symmetric. Definition 2.4 (Empirical Spectral Distribution): For a normalized Wigner matrix X n = M n / Ô n, we consider its n real eigenvalues 1 (X n ) Æ... Æ n (X n ) (the eigenvalues are real, because X n is symmetric). We call µ n the empirical spectral distribution, where µ n := 1 nÿ 1 { n n j Æx}, j=1 where n j = j (X n ). Notice that for each X n (Ê), where Ê œ function., µ n is a cumulative distribution Definition 2.5 (The Empirical Spectral Measure): We define the empirical spectral measure: S T µ n := E U 1 nÿ 1 { n V n j Æx}, defined by de[µ n ]=E dµ n for every compactly supported continuous real function. j=1 Definition 2.6 (Weak Convergence): A sequence n of deterministic probability measures on a space with associated -algebra F is said to converge weakly to a probability measure, if for a continuous, bounded function f, fd n æ fd 2
3 as n æœ.if n is random, then we consider weak convergence to either in probability or almost surely. 3
4 3 Eigenvalue Distribution of Wigner Matrices In this chapter, we aim to prove the following theorem, the main result of this project. Theorem 3.1: Let {M n } Œ n=1 be a sequence of Wigner matrices, and for each n, denote X n = M n / Ô n. Then µ n converges weakly, in probability to the semicircle distribution, (x)dx = 1 Ô 4 x2 1 x Æ2 dx. 2fi Proof: We want to show that c s fdµ n æ s f (x)dx in probability for any continuous bounded function. In other words, we want to show for every >0, P fdµ for any continuous bounded function. 4 f (x)dx - > æ 0 In order to rewrite the integrals in (3) as a linear combination of moments, we can replace f by polynomials using the Weierstrass Approximation Theorem. Let B>0 be a fixed constant, k>0 be any positive integer, and f be a continuous bounded function. By using the Markov s Inequality, we obtain Cl P 1 s x k 1 x >B dµ n > 2 Æ 1 E Ë s x k 1 x >B dµ n È Æ 1 E Ë s È x k x k dµ B k n = 1 E Ë s È B k x2k dµ n = 1 B k s x2k d µ n. By Lemma (4.2), we say that s xk d µ n converges (deterministically) to s xk (x)dx for any positive integer k. Then, we obtain c lim sup næœ P 1 s x k 1 x >B dµ n > 2 Æ 1 s B k x2k (x)dx. On the other hand, by Lemma (4.4), we know c s x2k (x)dx = C k. Also since C k Æ 4 k, combining (3) and (3), we obtain c lim sup næœ P 1 s x2k 1 x >B dµ n > 2 Æ 4k. B k Now, let a k = lim sup næœ P 1 s x2k 1 x >B dµ n > 2, b k = 4k B k, and choose B>4. We claim that a k = 0 for all k œ N. 4
5 Suppose the contrary, that there is k 0 œ N with a k0 > 0. But since a k is an increasing sequence, by definition, a k Ø a k0 for all k Ø k 0. In this case, we have 0 = lim kæœ b k Ø lim kæœ a k Ø a k0 > 0 which is a contradiction. So, a k = 0 for all k œ N. Now, let >0 and f : æ be bounded and continuous. Also, suppose that f is compactly supported on [ B,B]. Then, by the Weierstrass Approximation Theorem, we can choose a polynomial p such that p (x) f(x) Æ /4 for all x œ [ B,B]. Now, remembering f is supported on [ B,B] and using the triangle inequality, Cl s fdµ n s f (x)dx Æ s (f p )dµ n s (f p ) (x)dx + s p dµ n s p (x)dx Æ s f p dµ n + s f p (x)dx + s p dµ n s p (x)dx Æ s f p 1 x ÆB dµ n + s f p 1 x ÆB (x)dx s p 1 x >B dµ n s p dµ n s p (x)dx Since p (x) f(x) Æ /4 for all x œ [ B,B], substitution into (3) yields - fdµ n f (x)dx - Æ p 1 x >B dµ n + - p dµ n p (x)dx -. By using the triangle inequality, we get P fdµ n 4 A --- f (x)dx - > Æ P - Q Æ Pa p 1 x >B dµ n p 1 x >B dµ n + - p dµ n p d µ n p (x)dx - > B 2 p (x)dx p dµ n p d µ n- > b 2 A Æ P - p 1 x >B dµ n- > B 6 A P - p d µ n p (x)dx - > B 6 A P - p dµ n p d µ n- > B. (1) 6 Note that the event s fdµ n s f (x)dx > is contained in the union of the three events - s p x >B dµ n >, 6 s p µ n s p (x)dx >, and 6 s p µ n s p µ n >. 6 Now, from equation (3), the first probabilty in (1) goes to 0 as n æœ. By Lemma (4.2), the second probabilty converges to 0 as n æœ. Finally, the third probabilty also goes to 0 as n æœ, by Lemma (4.3). 5
6 Thus, for any >0, and any bounded function f, we have that c lim næœ P ( s fdµ n int f (x)dx > )=0. Now, by using Lemma 4.5, we can prove a stronger version of the result our main theorem which can be cobsidered as a corollart bof the theorem. Corollary: Let {M n } Œ n=1 be a sequence of Wigner matrices, and for each n, denote X n = M n / Ô n.t henµ n converges weakly, almost surely to the semicircle distribution. Proof: The convergence used in Lemma 4.2 is alsready in deterministic sense. So we just need to check whether the convergence in probabilty in Lemma 4.2 can be updated to the almost sure convergence. Firstly by using the Markov s Inequality and then using the Lemma 4.5, we get the following inequality Cl q Œ n=1 P 1- -s xk dµ n s xk d µ n > Æ 1 q Œn= s E 2 xk dµ n 1 E Ë È s xk dµ n Æ C q Œn=1 C < Œ. n 2 The constant C 1 is used because the inequality in Lemma 4.5 holds when n is su large. Now, by using the Borel-Cantelli Lemma, we immediately get that 3 P lim sup - x k dµ n næœ x k d µ n > =0. Then, this inequality means that s xk dµ n s xk d µ n - -- to0 asn æœalmost surely. ciently Now, again assuming f is compactly supported on [ B,B] and taking the polynomial p which is the approximation of f on [ B,B], we get as in the proof of the main theorem, Cl s fdµ n s f (x)dx Æs f p 1 x ÆB dµ n + s f p 1 x ÆB (x)dx + - s - p -- 1 x >B dµ n + s p d µ n s p (x)dx + s p dµ n s p d µ n. By the choose of p, we can make the first two terms in the above inequalirt arbitrarily small. The third and fourth terms approach 0, deterministically, and the last term converges to 0 almost surely. Therefore, for every bounded function f, s fdµ n s f (x)dx æ0 almost surely which means the emprical spectral distribution converges almost surely, weakly to the semicircle distribution. 6
7 4 Technical Theorems and Lemmas Theorem 4.1 (Weierstrass Approximation Theorem): If f is a continuous, real-valued function on [a, b] and if any >0 is given, there exists a polynomial p on [a, b] such that f(x) p(x) <, for all x œ [a, b]. In other words, any continuous function on a closed and bounded interval can be uniformly approximated on that interval by polynomials to any degree of accuracy. Lemma 4.2: For any positive integer k, s xk d µ n converges (deterministically) to s xk (x)dx as n æœ. Lemma 4.3: Let >0 be given and k be a positive integer. Then, lim P x k dµ n x k d µ n > =0. næœ Lemma 4.4: The moments of the semicircle law are given by Y ] x k C k/2 k is even, (x)dx =. [ 0 k is odd, 2 Note that C n represents the n th Catalan number, namely C n = n n n. Lemma 4.5: Let X n be a Wigner matrix with emprical sprectral distribution µ n. Then for every fixed k, there exists a constant C not depending on n such that for su ciently large n. 4 D 2 - C3 E x k dµ n E x k C dµ n Æ n, 2 Lemma 4.6 (Borel-Cantelli Lemma): Let {E n } Œ n=1 be a sequence of events in some probabilty space. If Œÿ P (E n ) < Œ, then n=1 3 4 P lim sup E n =0. næœ 7
8 1 1 Figure 1: Eigenvalues Gaussian andom Matrix on Complex Plane (n = 100 and n = 1000) Empirical esults 8
9 1 1 Figure 2: Eigenvalues Gaussian andom Matrix on Complex Plane (n = 2000 and n = 4000) 9
10 Figure 3: Eigenvalues of Wigner matrices n= 10, 100 and
11 Figure 4: Eigenvalues of a collection of 1, 10 and 100 Wigner matrices with fixed matrix size n=10 11
12 5 Mathematica Code of Simulations Normalized Symmetric andom Matrix andomsymmetricmatrix[dist_, n_] := Module[{mat = (1/(Sqrt[n*Variance[dist]]))* andomvariate[dist, {n, n}]}, UpperTriangularize[mat, 1] + Transpose[UpperTriangularize[mat]]] dists = {StudentTDistribution[4], NormalDistribution[], LaplaceDistribution[0, 1], ExponentialDistribution[1]}; ev = Eigenvalues[andomSymmetricMatrix[#, 10ˆ2]] & /@ dists; edist = EstimatedDistribution[#, WignerSemicircleDistribution[r]] & /@ ev; h[dist_, data_, i_] := Histogram[data, 20, "PDF", ChartStyle -> (ColorData["Gradients"][[andomInteger[{1, 51}]]]), BaseStyle -> {FontFamily -> "Verdana"}, PlotLabel -> dists[[i]], Plotange -> {0, 1.5 PDF[dist, 0]}, ImageSize -> 280, Epilog -> Inset[Framed[ Style[Grid[{{"Estimated Distribution:"}, {dist}}], 10], oundingadius -> 10, FrameStyle -> GrayLevel@0.3, Background -> LightOrange], {ight, 1.45 PDF[dist, 0]}, {ight, Top}]] distplot[dist_, data_] := Plot[PDF[dist, x], {x, Min[data], Max[data]}, PlotStyle -> {Thick, Black}] Grid[Partition[ Table[Show[h[edist[[i]], ev[[i]], i], distplot[edist[[i]], ev[[i]]]], {i, 4}], 2]] Normalized Gaussian andom Matrix n = {100, 1000, 2000, 4000}; dist = NormalDistribution[]; GaussianandomMatrix = (1/(Sqrt[#*Variance[dist]]))* andomvariate[dist, {#, #}] & /@ n; ev = Eigenvalues[#] & /@ GaussianandomMatrix; Table[Show[ ListPlot[{e[#], Im[#]} & /@ ev[[i]], Aspectatio -> 1], Graphics@Circle[{0, 0}, 1]], {i, 4}] 12
13 Collection of Normalized Symmetric andom Matrices andomsymmetricmatrix[dist_, n_] := Module[{mat = (1/(Sqrt[n*Variance[dist]]))* andomvariate[dist, {n, n}]}, UpperTriangularize[mat, 1] + Transpose[UpperTriangularize[mat]]] dists = NormalDistribution[]; collection = Table[Eigenvalues[andomSymmetricMatrix[dists, 10]], {i, 100}]; ev = Flatten[collection]; edist = EstimatedDistribution[ev, WignerSemicircleDistribution[r]]; h[dist_, data_] := Histogram[data, 20, "PDF", ChartStyle -> (ColorData["Gradients"][[andomInteger[{1, 51}]]]), BaseStyle -> {FontFamily -> "Verdana"}, PlotLabel -> dists, Plotange -> {0, 1.5 PDF[dist, 0]}, ImageSize -> 280, Epilog -> Inset[Framed[ Style[Grid[{{"Estimated Distribution:"}, {dist}}], 10], oundingadius -> 10, FrameStyle -> GrayLevel@0.3, Background -> LightOrange], {ight, 1.45 PDF[dist, 0]}, {ight, Top}]] distplot[dist_, data_] := Plot[PDF[dist, x], {x, Min[data], Max[data]}, PlotStyle -> {Thick, Black}] Show[h[edist, ev], distplot[edist, ev]] 6 eferences Feier, A.. (2012). Methods of Proof in andom Matrix Theory. Unpublished. Weisstein, Eric W. Weierstrass Approximation Theorem. From MathWorld A Wolfram Web esource. Anderson, Guionnet, Zeitouni, An Introduction to andom Matrices, pp
1 Intro to RMT (Gene)
M705 Spring 2013 Summary for Week 2 1 Intro to RMT (Gene) (Also see the Anderson - Guionnet - Zeitouni book, pp.6-11(?) ) We start with two independent families of R.V.s, {Z i,j } 1 i
More informationLectures 6 7 : Marchenko-Pastur Law
Fall 2009 MATH 833 Random Matrices B. Valkó Lectures 6 7 : Marchenko-Pastur Law Notes prepared by: A. Ganguly We will now turn our attention to rectangular matrices. Let X = (X 1, X 2,..., X n ) R p n
More informationThe norm of polynomials in large random matrices
The norm of polynomials in large random matrices Camille Mâle, École Normale Supérieure de Lyon, Ph.D. Student under the direction of Alice Guionnet. with a significant contribution by Dimitri Shlyakhtenko.
More informationFluctuations from the Semicircle Law Lecture 4
Fluctuations from the Semicircle Law Lecture 4 Ioana Dumitriu University of Washington Women and Math, IAS 2014 May 23, 2014 Ioana Dumitriu (UW) Fluctuations from the Semicircle Law Lecture 4 May 23, 2014
More informationRandom Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg
Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws Symeon Chatzinotas February 11, 2013 Luxembourg Outline 1. Random Matrix Theory 1. Definition 2. Applications 3. Asymptotics 2. Ensembles
More informationDistribution of Eigenvalues of Weighted, Structured Matrix Ensembles
Distribution of Eigenvalues of Weighted, Structured Matrix Ensembles Olivia Beckwith 1, Steven J. Miller 2, and Karen Shen 3 1 Harvey Mudd College 2 Williams College 3 Stanford University Joint Meetings
More informationDISTRIBUTION OF EIGENVALUES OF REAL SYMMETRIC PALINDROMIC TOEPLITZ MATRICES AND CIRCULANT MATRICES
DISTRIBUTION OF EIGENVALUES OF REAL SYMMETRIC PALINDROMIC TOEPLITZ MATRICES AND CIRCULANT MATRICES ADAM MASSEY, STEVEN J. MILLER, AND JOHN SINSHEIMER Abstract. Consider the ensemble of real symmetric Toeplitz
More informationHomework # , Spring Due 14 May Convergence of the empirical CDF, uniform samples
Homework #3 36-754, Spring 27 Due 14 May 27 1 Convergence of the empirical CDF, uniform samples In this problem and the next, X i are IID samples on the real line, with cumulative distribution function
More information2.1 Lecture 5: Probability spaces, Interpretation of probabilities, Random variables
Chapter 2 Kinetic Theory 2.1 Lecture 5: Probability spaces, Interpretation of probabilities, Random variables In the previous lectures the theory of thermodynamics was formulated as a purely phenomenological
More informationMath 421, Homework #7 Solutions. We can then us the triangle inequality to find for k N that (x k + y k ) (L + M) = (x k L) + (y k M) x k L + y k M
Math 421, Homework #7 Solutions (1) Let {x k } and {y k } be convergent sequences in R n, and assume that lim k x k = L and that lim k y k = M. Prove directly from definition 9.1 (i.e. don t use Theorem
More informationThe circular law. Lewis Memorial Lecture / DIMACS minicourse March 19, Terence Tao (UCLA)
The circular law Lewis Memorial Lecture / DIMACS minicourse March 19, 2008 Terence Tao (UCLA) 1 Eigenvalue distributions Let M = (a ij ) 1 i n;1 j n be a square matrix. Then one has n (generalised) eigenvalues
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationProblem 1 HW3. Question 2. i) We first show that for a random variable X bounded in [0, 1], since x 2 apple x, wehave
Next, we show that, for any x Ø, (x) Æ 3 x x HW3 Problem i) We first show that for a random variable X bounded in [, ], since x apple x, wehave Var[X] EX [EX] apple EX [EX] EX( EX) Since EX [, ], therefore,
More informationLectures 2 3 : Wigner s semicircle law
Fall 009 MATH 833 Random Matrices B. Való Lectures 3 : Wigner s semicircle law Notes prepared by: M. Koyama As we set up last wee, let M n = [X ij ] n i,j=1 be a symmetric n n matrix with Random entries
More informationLecture I: Asymptotics for large GUE random matrices
Lecture I: Asymptotics for large GUE random matrices Steen Thorbjørnsen, University of Aarhus andom Matrices Definition. Let (Ω, F, P) be a probability space and let n be a positive integer. Then a random
More informationOn the concentration of eigenvalues of random symmetric matrices
On the concentration of eigenvalues of random symmetric matrices Noga Alon Michael Krivelevich Van H. Vu April 23, 2012 Abstract It is shown that for every 1 s n, the probability that the s-th largest
More informationEigenvalue variance bounds for Wigner and covariance random matrices
Eigenvalue variance bounds for Wigner and covariance random matrices S. Dallaporta University of Toulouse, France Abstract. This work is concerned with finite range bounds on the variance of individual
More informationEigenvalues and Singular Values of Random Matrices: A Tutorial Introduction
Random Matrix Theory and its applications to Statistics and Wireless Communications Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction Sergio Verdú Princeton University National
More information3 Integration and Expectation
3 Integration and Expectation 3.1 Construction of the Lebesgue Integral Let (, F, µ) be a measure space (not necessarily a probability space). Our objective will be to define the Lebesgue integral R fdµ
More informationWigner s semicircle law
CHAPTER 2 Wigner s semicircle law 1. Wigner matrices Definition 12. A Wigner matrix is a random matrix X =(X i, j ) i, j n where (1) X i, j, i < j are i.i.d (real or complex valued). (2) X i,i, i n are
More information. Find E(V ) and var(v ).
Math 6382/6383: Probability Models and Mathematical Statistics Sample Preliminary Exam Questions 1. A person tosses a fair coin until she obtains 2 heads in a row. She then tosses a fair die the same number
More informationRandom Matrices: Beyond Wigner and Marchenko-Pastur Laws
Random Matrices: Beyond Wigner and Marchenko-Pastur Laws Nathan Noiry Modal X, Université Paris Nanterre May 3, 2018 Wigner s Matrices Wishart s Matrices Generalizations Wigner s Matrices ij, (n, i, j
More informationSequences. Chapter 3. n + 1 3n + 2 sin n n. 3. lim (ln(n + 1) ln n) 1. lim. 2. lim. 4. lim (1 + n)1/n. Answers: 1. 1/3; 2. 0; 3. 0; 4. 1.
Chapter 3 Sequences Both the main elements of calculus (differentiation and integration) require the notion of a limit. Sequences will play a central role when we work with limits. Definition 3.. A Sequence
More informationLecture 6 Basic Probability
Lecture 6: Basic Probability 1 of 17 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 6 Basic Probability Probability spaces A mathematical setup behind a probabilistic
More informationA Note on the Central Limit Theorem for the Eigenvalue Counting Function of Wigner and Covariance Matrices
A Note on the Central Limit Theorem for the Eigenvalue Counting Function of Wigner and Covariance Matrices S. Dallaporta University of Toulouse, France Abstract. This note presents some central limit theorems
More informationProblem set 1, Real Analysis I, Spring, 2015.
Problem set 1, Real Analysis I, Spring, 015. (1) Let f n : D R be a sequence of functions with domain D R n. Recall that f n f uniformly if and only if for all ɛ > 0, there is an N = N(ɛ) so that if n
More informationRectangular Young tableaux and the Jacobi ensemble
Rectangular Young tableaux and the Jacobi ensemble Philippe Marchal October 20, 2015 Abstract It has been shown by Pittel and Romik that the random surface associated with a large rectangular Young tableau
More informationA Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices
A Remark on Hypercontractivity and Tail Inequalities for the Largest Eigenvalues of Random Matrices Michel Ledoux Institut de Mathématiques, Université Paul Sabatier, 31062 Toulouse, France E-mail: ledoux@math.ups-tlse.fr
More informationTHE N-VALUE GAME OVER Z AND R
THE N-VALUE GAME OVER Z AND R YIDA GAO, MATT REDMOND, ZACH STEWARD Abstract. The n-value game is an easily described mathematical diversion with deep underpinnings in dynamical systems analysis. We examine
More informationarxiv: v2 [math.pr] 2 Nov 2009
arxiv:0904.2958v2 [math.pr] 2 Nov 2009 Limit Distribution of Eigenvalues for Random Hankel and Toeplitz Band Matrices Dang-Zheng Liu and Zheng-Dong Wang School of Mathematical Sciences Peking University
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More informationConvergence of spectral measures and eigenvalue rigidity
Convergence of spectral measures and eigenvalue rigidity Elizabeth Meckes Case Western Reserve University ICERM, March 1, 2018 Macroscopic scale: the empirical spectral measure Macroscopic scale: the empirical
More informationLectures 2 3 : Wigner s semicircle law
Fall 009 MATH 833 Random Matrices B. Való Lectures 3 : Wigner s semicircle law Notes prepared by: M. Koyama As we set up last wee, let M n = [X ij ] n i,j= be a symmetric n n matrix with Random entries
More informationLecture 2: Convergence of Random Variables
Lecture 2: Convergence of Random Variables Hyang-Won Lee Dept. of Internet & Multimedia Eng. Konkuk University Lecture 2 Introduction to Stochastic Processes, Fall 2013 1 / 9 Convergence of Random Variables
More informationPCA with random noise. Van Ha Vu. Department of Mathematics Yale University
PCA with random noise Van Ha Vu Department of Mathematics Yale University An important problem that appears in various areas of applied mathematics (in particular statistics, computer science and numerical
More informationConcentration inequalities for non-lipschitz functions
Concentration inequalities for non-lipschitz functions University of Warsaw Berkeley, October 1, 2013 joint work with Radosław Adamczak (University of Warsaw) Gaussian concentration (Sudakov-Tsirelson,
More information5 Birkhoff s Ergodic Theorem
5 Birkhoff s Ergodic Theorem Birkhoff s Ergodic Theorem extends the validity of Kolmogorov s strong law to the class of stationary sequences of random variables. Stationary sequences occur naturally even
More information9 Brownian Motion: Construction
9 Brownian Motion: Construction 9.1 Definition and Heuristics The central limit theorem states that the standard Gaussian distribution arises as the weak limit of the rescaled partial sums S n / p n of
More informationAnalysis III. Exam 1
Analysis III Math 414 Spring 27 Professor Ben Richert Exam 1 Solutions Problem 1 Let X be the set of all continuous real valued functions on [, 1], and let ρ : X X R be the function ρ(f, g) = sup f g (1)
More informationConcentration Inequalities for Random Matrices
Concentration Inequalities for Random Matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify the asymptotic
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 21, 2015 Francesco Franco Empirical Macroeconomics 1/33 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationSemicircle law on short scales and delocalization for Wigner random matrices
Semicircle law on short scales and delocalization for Wigner random matrices László Erdős University of Munich Weizmann Institute, December 2007 Joint work with H.T. Yau (Harvard), B. Schlein (Munich)
More informationOperator-Valued Free Probability Theory and Block Random Matrices. Roland Speicher Queen s University Kingston
Operator-Valued Free Probability Theory and Block Random Matrices Roland Speicher Queen s University Kingston I. Operator-valued semicircular elements and block random matrices II. General operator-valued
More informationConvergence of Random Variables
1 / 15 Convergence of Random Variables Saravanan Vijayakumaran sarva@ee.iitb.ac.in Department of Electrical Engineering Indian Institute of Technology Bombay March 19, 2014 2 / 15 Motivation Theorem (Weak
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationPart II Probability and Measure
Part II Probability and Measure Theorems Based on lectures by J. Miller Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationRandom Matrix Theory Lecture 3 Free Probability Theory. Symeon Chatzinotas March 4, 2013 Luxembourg
Random Matrix Theory Lecture 3 Free Probability Theory Symeon Chatzinotas March 4, 2013 Luxembourg Outline 1. Free Probability Theory 1. Definitions 2. Asymptotically free matrices 3. R-transform 4. Additive
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 5, 2016 Francesco Franco Empirical Macroeconomics 1/39 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
More informationLecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1
Random Walks and Brownian Motion Tel Aviv University Spring 011 Lecture date: May 0, 011 Lecture 9 Instructor: Ron Peled Scribe: Jonathan Hermon In today s lecture we present the Brownian motion (BM).
More information1 Introduction. 2 Measure theoretic definitions
1 Introduction These notes aim to recall some basic definitions needed for dealing with random variables. Sections to 5 follow mostly the presentation given in chapter two of [1]. Measure theoretic definitions
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationSTAT C206A / MATH C223A : Stein s method and applications 1. Lecture 31
STAT C26A / MATH C223A : Stein s method and applications Lecture 3 Lecture date: Nov. 7, 27 Scribe: Anand Sarwate Gaussian concentration recap If W, T ) is a pair of random variables such that for all
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 informationSTA205 Probability: Week 8 R. Wolpert
INFINITE COIN-TOSS AND THE LAWS OF LARGE NUMBERS The traditional interpretation of the probability of an event E is its asymptotic frequency: the limit as n of the fraction of n repeated, similar, and
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationDISTRIBUTION OF EIGENVALUES FOR THE ENSEMBLE OF REAL SYMMETRIC TOEPLITZ MATRICES
DISTRIBUTIO OF EIGEVALUES FOR THE ESEMBLE OF REAL SYMMETRIC TOEPLITZ MATRICES CHRISTOPHER HAMMOD AD STEVE J. MILLER Abstract. Consider the ensemble of real symmetric Toeplitz matrices, each independent
More informationTHE DYNAMICS OF SUCCESSIVE DIFFERENCES OVER Z AND R
THE DYNAMICS OF SUCCESSIVE DIFFERENCES OVER Z AND R YIDA GAO, MATT REDMOND, ZACH STEWARD Abstract. The n-value game is a dynamical system defined by a method of iterated differences. In this paper, we
More informationEdexcel GCE A Level Maths Further Maths 3 Matrices.
Edexcel GCE A Level Maths Further Maths 3 Matrices. Edited by: K V Kumaran kumarmathsweebly.com kumarmathsweebly.com 2 kumarmathsweebly.com 3 kumarmathsweebly.com 4 kumarmathsweebly.com 5 kumarmathsweebly.com
More informationConstructing Approximations to Functions
Constructing Approximations to Functions Given a function, f, if is often useful to it is often useful to approximate it by nicer functions. For example give a continuous function, f, it can be useful
More informationMeasurable Choice Functions
(January 19, 2013) Measurable Choice Functions Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/fun/choice functions.pdf] This note
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationEigenvalue Statistics for Toeplitz and Circulant Ensembles
Eigenvalue Statistics for Toeplitz and Circulant Ensembles Murat Koloğlu 1, Gene Kopp 2, Steven J. Miller 1, and Karen Shen 3 1 Williams College 2 University of Michigan 3 Stanford University http://www.williams.edu/mathematics/sjmiller/
More informationFrom random matrices to free groups, through non-crossing partitions. Michael Anshelevich
From random matrices to free groups, through non-crossing partitions Michael Anshelevich March 4, 22 RANDOM MATRICES For each N, A (N), B (N) = independent N N symmetric Gaussian random matrices, i.e.
More informationChapter 6. Convergence. Probability Theory. Four different convergence concepts. Four different convergence concepts. Convergence in probability
Probability Theory Chapter 6 Convergence Four different convergence concepts Let X 1, X 2, be a sequence of (usually dependent) random variables Definition 1.1. X n converges almost surely (a.s.), or with
More informationDefinitions & Theorems
Definitions & Theorems Math 147, Fall 2009 December 19, 2010 Contents 1 Logic 2 1.1 Sets.................................................. 2 1.2 The Peano axioms..........................................
More informationRandom matrices: Distribution of the least singular value (via Property Testing)
Random matrices: Distribution of the least singular value (via Property Testing) Van H. Vu Department of Mathematics Rutgers vanvu@math.rutgers.edu (joint work with T. Tao, UCLA) 1 Let ξ be a real or complex-valued
More informationThe Free Central Limit Theorem: A Combinatorial Approach
The Free Central Limit Theorem: A Combinatorial Approach by Dennis Stauffer A project submitted to the Department of Mathematical Sciences in conformity with the requirements for Math 4301 (Honour s Seminar)
More informationContinuity. Chapter 4
Chapter 4 Continuity Throughout this chapter D is a nonempty subset of the real numbers. We recall the definition of a function. Definition 4.1. A function from D into R, denoted f : D R, is a subset of
More informationRandom matrices: A Survey. Van H. Vu. Department of Mathematics Rutgers University
Random matrices: A Survey Van H. Vu Department of Mathematics Rutgers University Basic models of random matrices Let ξ be a real or complex-valued random variable with mean 0 and variance 1. Examples.
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More informationContinuity. Chapter 4
Chapter 4 Continuity Throughout this chapter D is a nonempty subset of the real numbers. We recall the definition of a function. Definition 4.1. A function from D into R, denoted f : D R, is a subset of
More informationACO Comprehensive Exam October 18 and 19, Analysis of Algorithms
Consider the following two graph problems: 1. Analysis of Algorithms Graph coloring: Given a graph G = (V,E) and an integer c 0, a c-coloring is a function f : V {1,,...,c} such that f(u) f(v) for all
More informationOn the distinguishability of random quantum states
1 1 Department of Computer Science University of Bristol Bristol, UK quant-ph/0607011 Distinguishing quantum states Distinguishing quantum states This talk Question Consider a known ensemble E of n quantum
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationRANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS
RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationExercises in Extreme value theory
Exercises in Extreme value theory 2016 spring semester 1. Show that L(t) = logt is a slowly varying function but t ǫ is not if ǫ 0. 2. If the random variable X has distribution F with finite variance,
More informationMath 104: Homework 7 solutions
Math 04: Homework 7 solutions. (a) The derivative of f () = is f () = 2 which is unbounded as 0. Since f () is continuous on [0, ], it is uniformly continous on this interval by Theorem 9.2. Hence for
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
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 informationFluctuations from the Semicircle Law Lecture 1
Fluctuations from the Semicircle Law Lecture 1 Ioana Dumitriu University of Washington Women and Math, IAS 2014 May 20, 2014 Ioana Dumitriu (UW) Fluctuations from the Semicircle Law Lecture 1 May 20, 2014
More informationSolutions Final Exam May. 14, 2014
Solutions Final Exam May. 14, 2014 1. Determine whether the following statements are true or false. Justify your answer (i.e., prove the claim, derive a contradiction or give a counter-example). (a) (10
More informationLocal semicircle law, Wegner estimate and level repulsion for Wigner random matrices
Local semicircle law, Wegner estimate and level repulsion for Wigner random matrices László Erdős University of Munich Oberwolfach, 2008 Dec Joint work with H.T. Yau (Harvard), B. Schlein (Cambrigde) Goal:
More informationBALANCING GAUSSIAN VECTORS. 1. Introduction
BALANCING GAUSSIAN VECTORS KEVIN P. COSTELLO Abstract. Let x 1,... x n be independent normally distributed vectors on R d. We determine the distribution function of the minimum norm of the 2 n vectors
More informationLUCK S THEOREM ALEX WRIGHT
LUCK S THEOREM ALEX WRIGHT Warning: These are the authors personal notes for a talk in a learning seminar (October 2015). There may be incorrect or misleading statements. Corrections welcome. 1. Convergence
More informationarxiv: v1 [math-ph] 19 Oct 2018
COMMENT ON FINITE SIZE EFFECTS IN THE AVERAGED EIGENVALUE DENSITY OF WIGNER RANDOM-SIGN REAL SYMMETRIC MATRICES BY G.S. DHESI AND M. AUSLOOS PETER J. FORRESTER AND ALLAN K. TRINH arxiv:1810.08703v1 [math-ph]
More informationLocal Kesten McKay law for random regular graphs
Local Kesten McKay law for random regular graphs Roland Bauerschmidt (with Jiaoyang Huang and Horng-Tzer Yau) University of Cambridge Weizmann Institute, January 2017 Random regular graph G N,d is the
More informationA Generalization of Wigner s Law
A Generalization of Wigner s Law Inna Zakharevich June 2, 2005 Abstract We present a generalization of Wigner s semicircle law: we consider a sequence of probability distributions (p, p 2,... ), with mean
More informationMonotonic ɛ-equilibria in strongly symmetric games
Monotonic ɛ-equilibria in strongly symmetric games Shiran Rachmilevitch April 22, 2016 Abstract ɛ-equilibrium allows for worse actions to be played with higher probability than better actions. I introduce
More informationAn almost sure invariance principle for additive functionals of Markov chains
Statistics and Probability Letters 78 2008 854 860 www.elsevier.com/locate/stapro An almost sure invariance principle for additive functionals of Markov chains F. Rassoul-Agha a, T. Seppäläinen b, a Department
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationhere, this space is in fact infinite-dimensional, so t σ ess. Exercise Let T B(H) be a self-adjoint operator on an infinitedimensional
15. Perturbations by compact operators In this chapter, we study the stability (or lack thereof) of various spectral properties under small perturbations. Here s the type of situation we have in mind:
More informationMA 123 (Calculus I) Lecture 3: September 12, 2017 Section A2. Professor Jennifer Balakrishnan,
What is on today Professor Jennifer Balakrishnan, jbala@bu.edu 1 Techniques for computing limits 1 1.1 Limit laws..................................... 1 1.2 One-sided limits..................................
More informationRandom Bernstein-Markov factors
Random Bernstein-Markov factors Igor Pritsker and Koushik Ramachandran October 20, 208 Abstract For a polynomial P n of degree n, Bernstein s inequality states that P n n P n for all L p norms on the unit
More information1.1 Limits and Continuity. Precise definition of a limit and limit laws. Squeeze Theorem. Intermediate Value Theorem. Extreme Value Theorem.
STATE EXAM MATHEMATICS Variant A ANSWERS AND SOLUTIONS 1 1.1 Limits and Continuity. Precise definition of a limit and limit laws. Squeeze Theorem. Intermediate Value Theorem. Extreme Value Theorem. Definition
More informationSTAT 7032 Probability Spring Wlodek Bryc
STAT 7032 Probability Spring 2018 Wlodek Bryc Created: Friday, Jan 2, 2014 Revised for Spring 2018 Printed: January 9, 2018 File: Grad-Prob-2018.TEX Department of Mathematical Sciences, University of Cincinnati,
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
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 informationLimit Laws for Random Matrices from Traffic Probability
Limit Laws for Random Matrices from Traffic Probability arxiv:1601.02188 Slides available at math.berkeley.edu/ bensonau Benson Au UC Berkeley May 9th, 2016 Benson Au (UC Berkeley) Random Matrices from
More information