Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich
|
|
- Blanche Ryan
- 5 years ago
- Views:
Transcription
1 Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008
2 MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. A linear functional µ[p ] = P (x) dµ(x), µ : R[x] R. Positive: µ[p (x) 2 ] 0. Inner product P, Q = µ[p Q]. Gram-Schmidt {, x, x 2, x 3,... } monic orthogonal polynomials {P 0 =, P, P 2, P 3,...}.
3 Theorem. (Favard, Stone, etc.) For some β i R, γ i 0 xp n = P n+ + β n P n + γ n P n. 2nd order recursion relation. Two independent solutions {P n, Q n }. Initial conditions P = 0, P 0 =, P = x β 0, Q 0 = 0, Q =, Q 2 = x β. Exercise. [ Pn (x) P n (y) Q n (x) = (I µ) x y ]. 2
4 µ { (β0, β, β 2,...) (γ, γ 2, γ 3,...) } equivalent. More explicit relation? If know { (β0, β, β 2,...) (γ, γ 2, γ 3,...) }, how to recover µ? Without going through {P n }. Cauchy transform G µ (z)= z x dµ(x) = µ [ z x ] = µ[] z + µ[x] z 2 + µ[x2 ] z 3 + µ[x3 ] z
5 G µ (z) = µ[] z + µ[x] z 2 + µ[x2 ] z 3 + µ[x3 ] z Theorem. Also G µ (z) = z β 0 z β Same coefficients as in the recursion. γ γ 2 z β 2 γ 3 z... Note: µ 0 all γ 0, no determinants. Proof II. Flajolet (980): lattice paths. 4
6 Proof I. Look at G(z) = z β 0 z β γ z β 2 γ 2 γ 3... z β n γ n H G = polynomial polynomial. Claim. G = Q n γ n Q n H P n γ n P n H, where {P, Q} with recursion {β i, γ i }. 5
7 Proof. By induction. Assume and say Then G = G= z β 0 γ H 0 = Q γ Q 0 H 0 P γ P 0 H 0. G = Q n γ n Q n H P n γ n P n H H = z β n γ n+ K. Q n γ n Q n z β n γ n+ K P n γ n P n z β n γ n+ K = zq n β n Q n γ n+ Q n K γ n Q n zp n β n P n γ n+ P n K γ n P n = Q n+ γ n+ Q n K P n+ γ n+ P n K. 6
8 FINITE CONTINUED FRACTIONS. Let G n = G cut off at level n. G n (z) = z β 0 z β γ z β 2 γ 2 γ 3... H = 0. z β n 0 G n = Q n γ n Q n H P n γ n P n H = Q n P n. 7
9 G n = Q n P n. P n monic, n real roots G n (z)= P n (x) = n i= (x x i ). z x dµ n(x) = Q n(z) (z xi ) = (partial fractions) = a i z x i. G n = Cauchy transform of µ n = a i δ xi, x i = roots of P n, a i = Q n(x i ) P n(x i ). 8
10 G n (z) = z β 0 z β γ z β 2 γ 2 γ 3... G n approximate G. z β n 0 G n = G µn, µ n = n i= a i δ xi. Do µ n approximate µ? Yes. In fact, µ n [P (x)] = µ[p (x)] for deg P 2n. 9
11 GAUSSIAN QUADRATURE. Want to evaluate f(x) dµ(x) Riemann sums. n i= a i f(x i ). Want P (x) dµ(x) = n i= a i P (x i ) for P of low degree. How to choose a i, x i? Answer: take x i = roots of P n. Choose a i so that x k dµ(x) = a i x k i, k = 0,,..., n (n equations, n unknowns). Our a i = Q n(x i ) P n(x i ) work. 0
12 Proof. Lagrange interpolation: for any P with deg P < n, Note so µ [ Pn (x) µ[p (x)] = x x i P (x) = n ] i= = µ n i= P (x i )P n (x) P n (x i)(x x i ). [ Pn (x) P n (x i ) x x i P (x i ) P n (x i) Q n(x i ) = n ] i= = Q n (x i ) Q n (x i ) P n (x i) δ x i [P ] µ[p (x)] = n i= a i δ xi [P (x)] µ[p (x)] = µ n [P (x)] for deg P < n. In fact, same x i, a i work for k = n, n +,..., 2n.
13 For P k, n k 2n, P k (x) = A(x)P n (x) + B(x), deg A, B n. µ[p k ] =, P k = 0. To show: µ n [P k ] = 0. µ[ap n ] = A, P n = 0 deg A < n µ[b] = 0 µ n [B] = 0. Finally, so µ n [P k ] = 0. µ n [AP n ] = a i A(x i )P n (x i ) = 0, So µ n µ, G n G, and therefore G µ = G. 2
14 If know {β i, γ i }, can find µ? Usually hard: G µ (z) = z β 0 an infinite expression. z β γ γ 2 z β 2 γ 3 z... Class of explicit examples. Semicircle law: 2πt 4t x 2 [ 2 t,2 t] dx y x 2 3 3
15 Semicircle law: 2πt 4t x 2 [ 2 t,2 t] dx y x 2 3 Marchenko-Pastur distributions: 2π 4t (x t) 2 x [+t 2 t,+t+2 t] (x) dx + max( t, 0)δ y 0.5 y x x 3 4 4
16 Semicircular, Marchenko-Pastur orthogonal polynomials satisfy n=0 P n (x)u n = A(u) B(u)x. In general: free Meixner distributions 2πt 4(t + c) (x b) 2 + (b/t)x + (c/t 2 )x 2 [b 2 t+c,b+2 t+c] dx +0,, 2 atoms. AC support an interval polynomial polynomial limited atoms outside of the AC support Not SC part 5
17 PERIODIC CONTINUED FRACTIONS. β i+n = β i, γ i+n = γ i. H = G in G(z) = z β 0 z β γ z β 2 γ 2 γ 3... z β n γ n H G = Q n γ n Q n G P n γ n P n G. γ n P n G 2 (γ n Q n + P n )G + Q n = 0. Quadratic equation! D = (γ n Q n + P n ) 2 4γ n Q n P n. 6
18 G = (γ nq n + P n ) D 2γ n P n = Stieltjes inversion formula: z x dµ(x). dµ(x) > 0 if dµ(x) = π lim y 0 Im G(x + iy). D(x) ir, i.e. D < 0. D degree 2n, n intervals for D 0. No SC. atoms: roots of P n, at most (n ). Recall D = (γ n Q n + P n ) 2 4γ n Q n P n. So if P n (a) = 0, then D(a) 0. Atoms outside of the AC support. 7
19 Eventually constant continued fractions (n = ) β i = β, γ i = γ for i N. polynomial on one interval. polynomial If β = 0, γ = (Peherstorfer?) Bernstein-Szegő class 4 x 2 polynomial on [ 2, 2]. Weyl s Theorem. If β i 0, γ i, then σ ess (µ) = [ 2, 2]. Denisov-Rakhmanov Theorem. If σ ess (µ) = AC support of µ = [ 2, 2], then β i 0, γ i. 8
20 Eventually periodic polynomial polynomial on n intervals. polynomial polynomial on n intervals eventually periodic. Weyl: if {β i, γ i } asymptotically periodic, same essential spectrum as for actually periodic. Converse false. Last, Simon: if {β i, γ i } approaches the isospectral torus of a periodic sequence, same essential spectrum. Damanik, Killip, Simon: converse true. 9
21 QUESTIONS. Connection between random matrices and (eventually) periodic continued fractions (Pastur). Multivariate (non-commutative) orthogonal polynomials, states, continued fractions. All exist. Continued fractions matricial. Formulas for states with periodic continued fractions. Free Meixner states known; constant after step 2. If not formulas for states, description of their operator algebras. Connection between multi-matrix models and states with multivariate (eventually) periodic continued fractions. 20
Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich
Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. 2 MEASURES
More informationCharacterizations of free Meixner distributions
Characterizations of free Meixner distributions Texas A&M University March 26, 2010 Jacobi parameters. Matrix. β 0 γ 0 0 0... 1 β 1 γ 1 0.. m n J =. 0 1 β 2 γ.. 2 ; J n =. 0 0 1 β.. 3............... A
More informationIntroduction to orthogonal polynomials. Michael Anshelevich
Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)
More informationFree Meixner distributions and random matrices
Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2
More informationSpectral Theory of Orthogonal Polynomials
Spectral Theory of Orthogonal Polynomials Barry Simon IBM Professor of Mathematics and Theoretical Physics California Institute of Technology Pasadena, CA, U.S.A. Lecture 1: Introduction and Overview Spectral
More informationRatio Asymptotics for General Orthogonal Polynomials
Ratio Asymptotics for General Orthogonal Polynomials Brian Simanek 1 (Caltech, USA) Arizona Spring School of Analysis and Mathematical Physics Tucson, AZ March 13, 2012 1 This material is based upon work
More informationLinearization coefficients for orthogonal polynomials. Michael Anshelevich
Linearization coefficients for orthogonal polynomials Michael Anshelevich February 26, 2003 P n = monic polynomials of degree n = 0, 1,.... {P n } = basis for the polynomials in 1 variable. Linearization
More informationContraction Principles some applications
Contraction Principles some applications Fabrice Gamboa (Institut de Mathématiques de Toulouse) 7th of June 2017 Christian and Patrick 59th Birthday Overview Appetizer : Christian and Patrick secret lives
More informationMath 4310 Solutions to homework 7 Due 10/27/16
Math 4310 Solutions to homework 7 Due 10/27/16 1. Find the gcd of x 3 + x 2 + x + 1 and x 5 + 2x 3 + x 2 + x + 1 in Rx. Use the Euclidean algorithm: x 5 + 2x 3 + x 2 + x + 1 = (x 3 + x 2 + x + 1)(x 2 x
More informationApproximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationCharacterizations of free Meixner distributions
Characterizations of free Meixner distributions Texas A&M University August 18, 2009 Definition via Jacobi parameters. β, γ, b, c R, 1 + γ, 1 + c 0. Tridiagonal matrix {(β, b, b,...), (1 + γ, 1 + c, 1
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationThe Gram matrix in inner product modules over C -algebras
The Gram matrix in inner product modules over C -algebras Ljiljana Arambašić (joint work with D. Bakić and M.S. Moslehian) Department of Mathematics University of Zagreb Applied Linear Algebra May 24 28,
More informationSection 33 Finite fields
Section 33 Finite fields Instructor: Yifan Yang Spring 2007 Review Corollary (23.6) Let G be a finite subgroup of the multiplicative group of nonzero elements in a field F, then G is cyclic. Theorem (27.19)
More informationORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS
ORTHOGONAL POLYNOMIALS WITH EXPONENTIALLY DECAYING RECURSION COEFFICIENTS BARRY SIMON* Dedicated to S. Molchanov on his 65th birthday Abstract. We review recent results on necessary and sufficient conditions
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
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 informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More information13. Examples of measure-preserving tranformations: rotations of a torus, the doubling map
3. Examples of measure-preserving tranformations: rotations of a torus, the doubling map 3. Rotations of a torus, the doubling map In this lecture we give two methods by which one can show that a given
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationModern Computer Algebra
Modern Computer Algebra Exercises to Chapter 25: Fundamental concepts 11 May 1999 JOACHIM VON ZUR GATHEN and JÜRGEN GERHARD Universität Paderborn 25.1 Show that any subgroup of a group G contains the neutral
More informationPolynomial Review Problems
Polynomial Review Problems 1. Find polynomial function formulas that could fit each of these graphs. Remember that you will need to determine the value of the leading coefficient. The point (0,-3) is on
More informationElliptic Curves and Public Key Cryptography
Elliptic Curves and Public Key Cryptography Jeff Achter January 7, 2011 1 Introduction to Elliptic Curves 1.1 Diophantine equations Many classical problems in number theory have the following form: Let
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week: December 4 8 Deadline to hand in the homework: your exercise class on week January 5. Exercises with solutions ) Let H, K be Hilbert spaces, and A : H K be a linear
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 informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationChapter 8. P-adic numbers. 8.1 Absolute values
Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.
More informationAlgebra Review 2. 1 Fields. A field is an extension of the concept of a group.
Algebra Review 2 1 Fields A field is an extension of the concept of a group. Definition 1. A field (F, +,, 0 F, 1 F ) is a set F together with two binary operations (+, ) on F such that the following conditions
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationLEAST SQUARES APPROXIMATION
LEAST SQUARES APPROXIMATION One more approach to approximating a function f (x) on an interval a x b is to seek an approximation p(x) with a small average error over the interval of approximation. A convenient
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationYORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #2 Solutions
YORK UNIVERSITY Faculty of Science Department of Mathematics and Statistics MATH 3. M Test # Solutions. (8 pts) For each statement indicate whether it is always TRUE or sometimes FALSE. Note: For this
More informationSolutions for Math 225 Assignment #5 1
Solutions for Math 225 Assignment #5 1 (1) Find a polynomial f(x) of degree at most 3 satisfying that f(0) = 2, f( 1) = 1, f(1) = 3 and f(3) = 1. Solution. By Lagrange Interpolation, ( ) (x + 1)(x 1)(x
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 informationKernel families of probability measures. Saskatoon, October 21, 2011
Kernel families of probability measures Saskatoon, October 21, 2011 Abstract The talk will compare two families of probability measures: exponential, and Cauchy-Stjelties families. The exponential families
More information0.1 Rational Canonical Forms
We have already seen that it is useful and simpler to study linear systems using matrices. But matrices are themselves cumbersome, as they are stuffed with many entries, and it turns out that it s best
More informationHermite Interpolation and Sobolev Orthogonality
Acta Applicandae Mathematicae 61: 87 99, 2000 2000 Kluwer Academic Publishers Printed in the Netherlands 87 Hermite Interpolation and Sobolev Orthogonality ESTHER M GARCÍA-CABALLERO 1,, TERESA E PÉREZ
More information10/22/16. 1 Math HL - Santowski SKILLS REVIEW. Lesson 15 Graphs of Rational Functions. Lesson Objectives. (A) Rational Functions
Lesson 15 Graphs of Rational Functions SKILLS REVIEW! Use function composition to prove that the following two funtions are inverses of each other. 2x 3 f(x) = g(x) = 5 2 x 1 1 2 Lesson Objectives! The
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationi x i y i
Department of Mathematics MTL107: Numerical Methods and Computations Exercise Set 8: Approximation-Linear Least Squares Polynomial approximation, Chebyshev Polynomial approximation. 1. Compute the linear
More informationOrthogonal Polynomials on the Unit Circle
American Mathematical Society Colloquium Publications Volume 54, Part 2 Orthogonal Polynomials on the Unit Circle Part 2: Spectral Theory Barry Simon American Mathematical Society Providence, Rhode Island
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 informationZeros and ratio asymptotics for matrix orthogonal polynomials
Zeros and ratio asymptotics for matrix orthogonal polynomials Steven Delvaux, Holger Dette August 25, 2011 Abstract Ratio asymptotics for matrix orthogonal polynomials with recurrence coefficients A n
More informationHankel determinants, continued fractions, orthgonal polynomials, and hypergeometric series
Hankel determinants, continued fractions, orthgonal polynomials, and hypergeometric series Ira M. Gessel with Jiang Zeng and Guoce Xin LaBRI June 8, 2007 Continued fractions and Hankel determinants There
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
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 informationFinal A. Problem Points Score Total 100. Math115A Nadja Hempel 03/23/2017
Final A Math115A Nadja Hempel 03/23/2017 nadja@math.ucla.edu Name: UID: Problem Points Score 1 10 2 20 3 5 4 5 5 9 6 5 7 7 8 13 9 16 10 10 Total 100 1 2 Exercise 1. (10pt) Let T : V V be a linear transformation.
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More information2a 2 4ac), provided there is an element r in our
MTH 310002 Test II Review Spring 2012 Absractions versus examples The purpose of abstraction is to reduce ideas to their essentials, uncluttered by the details of a specific situation Our lectures built
More informationChapter 1. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
More informationMathematical Methods for Engineers and Scientists 1
K.T. Tang Mathematical Methods for Engineers and Scientists 1 Complex Analysis, Determinants and Matrices With 49 Figures and 2 Tables fyj Springer Part I Complex Analysis 1 Complex Numbers 3 1.1 Our Number
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 informationMTH310 EXAM 2 REVIEW
MTH310 EXAM 2 REVIEW SA LI 4.1 Polynomial Arithmetic and the Division Algorithm A. Polynomial Arithmetic *Polynomial Rings If R is a ring, then there exists a ring T containing an element x that is not
More informationLecture 7.5: Euclidean domains and algebraic integers
Lecture 7.5: Euclidean domains and algebraic integers Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4120, Modern Algebra M. Macauley
More information1. Group Theory Permutations.
1.1. Permutations. 1. Group Theory Problem 1.1. Let G be a subgroup of S n of index 2. Show that G = A n. Problem 1.2. Find two elements of S 7 that have the same order but are not conjugate. Let π S 7
More information6.3 Partial Fractions
6.3 Partial Fractions Mark Woodard Furman U Fall 2009 Mark Woodard (Furman U) 6.3 Partial Fractions Fall 2009 1 / 11 Outline 1 The method illustrated 2 Terminology 3 Factoring Polynomials 4 Partial fraction
More informationChapter 2 Orthogonal Polynomials and Weighted Polynomial Approximation
Chapter 2 Orthogonal Polynomials and Weighted Polynomial Approximation 2.1 Orthogonal Systems and Polynomials 2.1.1 Inner Product Space and Orthogonal Systems Suppose that X is a complex linear space of
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationFree Probability Theory and Random Matrices. Roland Speicher Queen s University Kingston, Canada
Free Probability Theory and Random Matrices Roland Speicher Queen s University Kingston, Canada We are interested in the limiting eigenvalue distribution of N N random matrices for N. Usually, large N
More informationMath 2: Algebra 2, Geometry and Statistics Ms. Sheppard-Brick Chapter 4 Test Review
Chapter 4 Test Review Students will be able to (SWBAT): Write an explicit and a recursive function rule for a linear table of values. Write an explicit function rule for a quadratic table of values. Determine
More informationPolynomials. Chapter 4
Chapter 4 Polynomials In this Chapter we shall see that everything we did with integers in the last Chapter we can also do with polynomials. Fix a field F (e.g. F = Q, R, C or Z/(p) for a prime p). Notation
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 informationφ(xy) = (xy) n = x n y n = φ(x)φ(y)
Groups 1. (Algebra Comp S03) Let A, B and C be normal subgroups of a group G with A B. If A C = B C and AC = BC then prove that A = B. Let b B. Since b = b1 BC = AC, there are a A and c C such that b =
More informationMeasurable functions are approximately nice, even if look terrible.
Tel Aviv University, 2015 Functions of real variables 74 7 Approximation 7a A terrible integrable function........... 74 7b Approximation of sets................ 76 7c Approximation of functions............
More information1 Linear Algebra Problems
Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and
More informationBMT 2016 Orthogonal Polynomials 12 March Welcome to the power round! This year s topic is the theory of orthogonal polynomials.
Power Round Welcome to the power round! This year s topic is the theory of orthogonal polynomials. I. You should order your papers with the answer sheet on top, and you should number papers addressing
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationLinear DifferentiaL Equation
Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace
More information(U) =, if 0 U, 1 U, (U) = X, if 0 U, and 1 U. (U) = E, if 0 U, but 1 U. (U) = X \ E if 0 U, but 1 U. n=1 A n, then A M.
1. Abstract Integration The main reference for this section is Rudin s Real and Complex Analysis. The purpose of developing an abstract theory of integration is to emphasize the difference between the
More informationDEPARTMENT OF MATHEMATICS
DEPARTMENT OF MATHEMATICS Ma322 - Final Exam Spring 2011 May 3,4, 2011 DO NOT TURN THIS PAGE UNTIL YOU ARE INSTRUCTED TO DO SO. Be sure to show all work and justify your answers. There are 8 problems and
More informationPERTURBATIONS OF ORTHOGONAL POLYNOMIALS WITH PERIODIC RECURSION COEFFICIENTS arxiv:math/ v2 [math.sp] 6 Dec 2008
PERTURBATIONS OF ORTHOGONAL POLYNOMIALS WITH PERIODIC RECURSION COEFFICIENTS arxiv:math/0702388v2 [math.sp] 6 Dec 2008 DAVID DAMANIK 1, ROWAN KILLIP 2, AND BARRY SIMON 3 Abstract. We extend the results
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS
ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS arxiv:083666v [mathpr] 8 Jan 009 Abstract We characterize, up to a conjecture, probability distributions of finite all order moments
More informationMoreover this binary operation satisfies the following properties
Contents 1 Algebraic structures 1 1.1 Group........................................... 1 1.1.1 Definitions and examples............................. 1 1.1.2 Subgroup.....................................
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More information(a + b)c = ac + bc and a(b + c) = ab + ac.
2. R I N G S A N D P O LY N O M I A L S The study of vector spaces and linear maps between them naturally leads us to the study of rings, in particular the ring of polynomials F[x] and the ring of (n n)-matrices
More informationProblem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that
Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume
More informationJacobi-Angelesco multiple orthogonal polynomials on an r-star
M. Leurs Jacobi-Angelesco m.o.p. 1/19 Jacobi-Angelesco multiple orthogonal polynomials on an r-star Marjolein Leurs, (joint work with Walter Van Assche) Conference on Orthogonal Polynomials and Holomorphic
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationLECTURE 7. k=1 (, v k)u k. Moreover r
LECTURE 7 Finite rank operators Definition. T is said to be of rank r (r < ) if dim T(H) = r. The class of operators of rank r is denoted by K r and K := r K r. Theorem 1. T K r iff T K r. Proof. Let T
More informationOrthogonal Polynomials and Gaussian Quadrature
Orthogonal Polynomials and Gaussian Quadrature 1. Orthogonal polynomials Given a bounded, nonnegative, nondecreasing function w(x) on an interval, I of the real line, we consider the Hilbert space L 2
More informationSolving Quadratic Equations
Solving Quadratic Equations MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: solve quadratic equations by factoring, solve quadratic
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More informationGQE ALGEBRA PROBLEMS
GQE ALGEBRA PROBLEMS JAKOB STREIPEL Contents. Eigenthings 2. Norms, Inner Products, Orthogonality, and Such 6 3. Determinants, Inverses, and Linear (In)dependence 4. (Invariant) Subspaces 3 Throughout
More informationP AC COMMUTATORS AND THE R TRANSFORM
Communications on Stochastic Analysis Vol. 3, No. 1 (2009) 15-31 Serials Publications www.serialspublications.com P AC COMMUTATORS AND THE R TRANSFORM AUREL I. STAN Abstract. We develop an algorithmic
More information1. General Vector Spaces
1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule
More informationDefinition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition
6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More information18.S34 (FALL 2007) PROBLEMS ON ROOTS OF POLYNOMIALS
18.S34 (FALL 2007) PROBLEMS ON ROOTS OF POLYNOMIALS Note. The terms root and zero of a polynomial are synonyms. Those problems which appeared on the Putnam Exam are stated as they appeared verbatim (except
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationTropical Polynomials
1 Tropical Arithmetic Tropical Polynomials Los Angeles Math Circle, May 15, 2016 Bryant Mathews, Azusa Pacific University In tropical arithmetic, we define new addition and multiplication operations on
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More information1 Fourier transform as unitary equivalence
Tel Aviv University, 009 Intro to functional analysis 1 1 Fourier transform as unitary equivalence 1a Introduction..................... 1 1b Exponential map................... 1c Exponential map as an
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More information