Avram Parter and Szegő limit theorems turned inside out

Size: px
Start display at page:

Download "Avram Parter and Szegő limit theorems turned inside out"

Transcription

1 Avram Parter and Szegő limit theorems turned inside out Egor A. Maximenko joint work with Johan M. Bogoya Ramírez, Albrecht Böttcher, and Sergei M. Grudsky Instituto Politécnico Nacional, ESFM, México Norte-Sur CINVESTAV, México 20 de noviembre de 204 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur 204 / 32

2 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

3 Toeplitz matrices Szegő s limit theorem Quantile function Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

4 Toeplitz matrices Szegő s limit theorem Quantile function Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

5 Toeplitz matrices T 5 (a) = a 0 a a 2 a 3 a 4 a a 0 a a 2 a 3 a 2 a a 0 a a 2 a 3 a 2 a a 0 a a 4 a 3 a 2 a a 0. It is convenient to think that a j are the Fourier coefficients of a certain function a defined on [0, 2π]: a j = 2π a(θ) e jiθ dθ. 2π 0 The function a is called the generating symbol of the matrices T n (a) = [ a j k ] n j,k=. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

6 Hermitian Toeplitz matrices, real bounded symbols We suppose that the generating symbol is bounded and real: a L ([0, 2π], R). In this case the matrices are Hermitian: a k = a k, a 0 R. T 5 (a) = a 0 a a 2 a 3 a 4 a a 0 a a 2 a 3 a 2 a a 0 a a 2 a 3 a 2 a a 0 a a 4 a 3 a 2 a a 0. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

7 stealed from loveiscomix.com Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

8 [ a0 ] a a 2 a 3 a a 0 a a 2 a 2 a a 0 a a 3 a 2 a a 0 To understand an Hermitian matrix is... Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

9 [ a0 ] a a 2 a 3 a a 0 a a 2 a 2 a a 0 a a 3 a 2 a a 0 [ λ ] λ λ λ 4 To understand an Hermitian matrix is to know the behavior of its eigenvalues and eigenvectors. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

10 Eigenvalues of Hermitian Toeplitz matrices For a in L ([0, 2π], R), sp T n (a) [ess inf(a), ess sup(a)]. Widom (994): sp T n (a) in the Hausdorff distance [ess inf(a), ess sup(a)]. λ (n),..., λ(n) n := the eigenvalues of T n (a) written in the ascending order: λ (n)... λ (n) n. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

11 Behavior of the eigenvalues of Hermitian Toeplitz matrices Graph of a Eigenvalues of T 8 (a) 3/4 /4 0 π/2 π 3π/2 2π 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

12 Behavior of the eigenvalues of Hermitian Toeplitz matrices Graph of a Eigenvalues of T 6 (a) 3/4 /4 0 π/2 π 3π/2 2π 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

13 Behavior of the eigenvalues of Hermitian Toeplitz matrices Graph of a Eigenvalues of T 32 (a) 3/4 /4 0 π/2 π 3π/2 2π 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

14 Behavior of the eigenvalues of Hermitian Toeplitz matrices Graph of a Eigenvalues of T 32 (a) 3/4 β=0.4 /4 how many? 0 π/2 π 3π/2 2π α=0. 0 First question: How many eigenvalues are in [α, β]? (Szegő, 920). Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

15 Behavior of the eigenvalues of Hermitian Toeplitz matrices Graph of a Eigenvalues of T 32 (a) 3/4 β=0.4 /4 0 π/2 π 3π/2 2π how many? α=0. 0 λ (32)? 7 First question: How many eigenvalues are in [α, β]? (Szegő, 920). Second question: λ (n) j? Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

16 Toeplitz matrices Szegő s limit theorem Quantile function Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

17 Szegő s limit theorem (920) generating symbol a L ([0, 2π], R) test function ϕ C(R) n n j= ϕ(λ (n) j ) 2π ϕ(a(θ)) dθ 2π 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur 204 / 32

18 Avram Parter limit theorem (988) an analogue of Szegő limit theorem for the singular values generating symbol a L ([0, 2π], C) test function ϕ C(R) n n j= ϕ(s (n) j ) 2π ϕ( a(θ) ) dθ 2π 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

19 Corollary from Szegő s limit theorem distribution of the eigenvalues of Hermitian Toepltiz matrices generating symbol a C([0, 2π], R) segment [α, β] α < β #{j : α λ (n) j β} n µ {θ [0, 2π]: α a(θ) β} 2π Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

20 Example to illustrate the corollary Graph of a Eigenvalues of T 32 (a) π 0. 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

21 Example to illustrate the corollary Graph of a Eigenvalues of T 32 (a) 0.4 eigenvalues 0 2π Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

22 Example to illustrate the corollary Graph of a Eigenvalues of T 32 (a) 0.4 eigenvalues 0 2π Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

23 Example to illustrate the corollary Graph of a Eigenvalues of T 32 (a) 0.4 eigenvalues 0 2π Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

24 Example to illustrate the corollary Graph of a Eigenvalues of T 32 (a) 0.4 eigenvalues 0 2π 0. 0 µ {θ : 0. a(θ) 0.4} 2π = Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

25 Szegő found an approximate answer to the first question (the number of the eigenvalues on any segment [α, β]). The second question is open: λ (n) j? Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

26 Corollary by Trench (202): the reordering idea a C([0, 2π], R) n n λ (n) j j= v (n) j 0 where v (n),..., v (n) ( ) ( ) ( ) n are the values a 2π n, a 4π n,..., a 2πn n written in the ascending order: v (n) v n (n) Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

27 Toeplitz matrices Szegő s limit theorem Quantile function Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

28 Quantile function of a list of numbers The same numbers in the ascending order: QuantileFunction(/3) = 8 because 8 is the minimal number v such that at least /3 of the elements are less or equal to v. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

29 Quantile function associated to a L ([0, 2π], R) F a := the cumulative distribution function of a : F a (v) := µ {θ [0, 2π]: a(θ) v}, v R. 2π Q a := the corresponding quantile function : Q a (p) := inf{v R: F a (v) p}, p (0, ]. Q a increases and has the same distribution as a. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

30 Construction of the quantile function associated to a piecewise-linear generating symbol Graph of a Graph of Q a π 2 π 3π 2 2π Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

31 Construction of the quantile function associated to a piecewise-linear generating symbol Graph of a Graph of Q a 3 4 v v 4 0 π 2 π 3π 2 2π a and Q a are identically distributed: 2π µ{θ [0, 2π]: a(θ) v} = µ{p (0, ]: Q a(p) v} Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

32 Construction of the quantile function associated to a piecewise-linear generating symbol Graph of a Graph of Q a π 2 π 3π 2 2π a reordering in the Lebesgue-style Q a Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

33 Q a may be continuous even when a is not 0 π 2π Graph of a Graph of Q a Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

34 Criterion of continuity of the quantile function Let a : [0, 2π] R be a measurable function. R(a) := the essential range of a. Then the following conditions are equivalent. R(a) is connected F a is strictly increasing Q a is continuous Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

35 Toeplitz matrices Szegő s limit theorem Quantile function Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

36 Main result: uniform convergence of the eigenvalues a L ([0, 2π], R) R(a) is connected max j n λ (n) j Q a ( j/n) 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

37 Idea of the proof Szegő s limit theorem convergence in distribution portmanteau lemma convergence of the quantile functions Uniform continuity of Q a Monotonicity: Q a λ (n)... λ (n) n Uniform convergence Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

38 Uniform convergence of the singular values (our corollary from Avram Parter theorem) a L ([0, 2π], C) R( a ) = [0, a ] max j n s (n) j Q a ( j/n) 0 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

39 First example continuous piecewise-linear generating symbol Graph of a Graph of Q a Eigenvalues of T 64 (a) 3/4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 9/48 /6 64 Every eigenvalue λ (n) j is shown as a point The third picture mimics the second one. ( ) j n, λ(n) j. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

40 First example continuous piecewise-linear generating symbol Graph of a Graph of Q a 3/4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 9/48 /6 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

41 First example continuous piecewise-linear generating symbol Graph of a Graph of Q a 3/4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 9/48 /6 and the points ( j/n, λ (n) j ) Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

42 First example continuous piecewise-linear generating symbol Graph of a 3/4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 9/48 /6 ( ) 3 64, λ(64) 3 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

43 First example continuous piecewise-linear generating symbol Graph of a ( 3 64, Q a ( )) /4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 9/48 /6 ( ) 3 64, λ(64) 3 Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

44 Second example trigonometric polynomial 8/9 Graph of a a(θ) = cos(θ) 2 9 cos(2θ) 8/9 Graph of Q a Eigenvalues of T 64 (a) 0 2π 3 π 4π 3 2π 0 /2 64 If the generating symbol is a trigonometric polynomial, then the maximal error is of the order O(/n): C > 0 max λ (n) j n j Q a ( j/n) C n Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

45 Third example a is not continuous, but R(a) is connected Graph of a Graph of Q a Eigenvalues of T 64 (a) 3/4 3/4 /4 /4 0 π/2 π 3π/2 2π 0 /6 5/6 64 In this example, λ (n) j is also uniformly approximated by Q a ( j/n) as n. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

46 Fourth example If R(a) is not connected, then the uniform convergence fails Graph of a Graph of Q a Eigenvalues of T 64 (a) 3/4 3/4 /4 /4 0 π 2π 0 /2 64 In this example, λ (n) n/2 can not be approximated by values of Q a. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

47 Summary The Szegő limit theorem combined with the notion of quantile function yields the main term of the individual asymptotics of the eigenvalues: λ (n) j Q a ( j n ) assuming that a L ([0, 2π], R) and R(a) is connected. Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

48 Summary The Szegő limit theorem combined with the notion of quantile function yields the main term of the individual asymptotics of the eigenvalues: λ (n) j Q a ( j n ) assuming that a L ([0, 2π], R) and R(a) is connected. Thanks for attention! Egor Maximenko (IPN-ESFM, Mexico) Avram Parter and Szegő theorems Norte-Sur / 32

Uniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices

Uniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices Uniform individual asymptotics for the eigenvalues and eigenvectors of large Toeplitz matrices Sergei Grudsky CINVESTAV, Mexico City, Mexico The International Workshop WIENER-HOPF METHOD, TOEPLITZ OPERATORS,

More information

Toeplitz matrices. Niranjan U N. May 12, NITK, Surathkal. Definition Toeplitz theory Computational aspects References

Toeplitz matrices. Niranjan U N. May 12, NITK, Surathkal. Definition Toeplitz theory Computational aspects References Toeplitz matrices Niranjan U N NITK, Surathkal May 12, 2010 Niranjan U N (NITK, Surathkal) Linear Algebra May 12, 2010 1 / 15 1 Definition Toeplitz matrix Circulant matrix 2 Toeplitz theory Boundedness

More information

Cayley-Hamilton Theorem

Cayley-Hamilton Theorem Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A

More information

arxiv: v1 [math.na] 21 Jan 2019

arxiv: v1 [math.na] 21 Jan 2019 Approximating the Perfect Sampling Grids for Computing the Eigenvalues of Toeplitz-like Matrices Using the Spectral Symbol Sven-Erik Ekström see@2pi.se Athens University of Economics and Business January

More information

Pseudospectra and Nonnormal Dynamical Systems

Pseudospectra and Nonnormal Dynamical Systems Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures

More information

RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS

RANDOM 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 information

Weak convergence. Amsterdam, 13 November Leiden University. Limit theorems. Shota Gugushvili. Generalities. Criteria

Weak convergence. Amsterdam, 13 November Leiden University. Limit theorems. Shota Gugushvili. Generalities. Criteria Weak Leiden University Amsterdam, 13 November 2013 Outline 1 2 3 4 5 6 7 Definition Definition Let µ, µ 1, µ 2,... be probability measures on (R, B). It is said that µ n converges weakly to µ, and we then

More information

Notes: 1. Regard as the maximal output error and as the corresponding maximal input error

Notes: 1. Regard as the maximal output error and as the corresponding maximal input error Limits and Continuity One of the major tasks in analysis is to classify a function by how nice it is Of course, nice often depends upon what you wish to do or model with the function Below is a list of

More information

The Fisher-Hartwig Conjecture and Toeplitz Eigenvalues

The Fisher-Hartwig Conjecture and Toeplitz Eigenvalues The Fisher-Hartwig Conjecture and Toeplitz Eigenvalues Estelle L. Basor Kent E. Morrison Department of Mathematics California Polytechnic State University San Luis Obispo, CA 93407 6 September 99 Abstract

More information

Problems Session. Nikos Stylianopoulos University of Cyprus

Problems Session. Nikos Stylianopoulos University of Cyprus [ 1 ] University of Cyprus Problems Session Nikos Stylianopoulos University of Cyprus Hausdorff Geometry Of Polynomials And Polynomial Sequences Institut Mittag-Leffler Djursholm, Sweden May-June 2018

More information

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x BASIC ALGORITHMS IN LINEAR ALGEBRA STEVEN DALE CUTKOSKY Matrices and Applications of Gaussian Elimination Systems of Equations Suppose that A is an n n matrix with coefficents in a field F, and x = (x,,

More information

Orthogonal Symmetric Toeplitz Matrices

Orthogonal Symmetric Toeplitz Matrices Orthogonal Symmetric Toeplitz Matrices Albrecht Böttcher In Memory of Georgii Litvinchuk (1931-2006 Abstract We show that the number of orthogonal and symmetric Toeplitz matrices of a given order is finite

More information

the configuration space of Branched polymers R. Kenyon, P. Winkler

the configuration space of Branched polymers R. Kenyon, P. Winkler the configuration space of Branched polymers R. Kenyon, P. Winkler artificial blood catalyst recovery Branched polymers (dendrimers) in modern science artificial photosynthesis A branched polymer is a

More information

Zeros of Polynomials: Beware of Predictions from Plots

Zeros of Polynomials: Beware of Predictions from Plots [ 1 / 27 ] University of Cyprus Zeros of Polynomials: Beware of Predictions from Plots Nikos Stylianopoulos a report of joint work with Ed Saff Vanderbilt University May 2006 Five Plots Fundamental Results

More information

Sobolev Duals of Random Frames

Sobolev Duals of Random Frames Sobolev Duals of Random Frames C. Sinan Güntürk, Mark Lammers, 2 Alex Powell, 3 Rayan Saab, 4 Özgür Yılmaz 4 Courant Institute of Mathematical Sciences, New York University, NY, USA. 2 University of North

More information

COFACTORS AND EIGENVECTORS OF BANDED TOEPLITZ MATRICES: TRENCH FORMULAS VIA SKEW SCHUR POLYNOMIALS. 1. Introduction

COFACTORS AND EIGENVECTORS OF BANDED TOEPLITZ MATRICES: TRENCH FORMULAS VIA SKEW SCHUR POLYNOMIALS. 1. Introduction Operators and Matrices Volume, Number 4 (207), 49 69 doi:0.753/oam-207--79 COFACTORS AND EIGENVECTORS OF BANDED TOEPLITZ MATRICES: TRENCH FORMULAS VIA SKEW SCHUR POLYNOMIALS EGOR A. MAXIMENKO AND MARIO

More information

Spectral properties of Toeplitz+Hankel Operators

Spectral properties of Toeplitz+Hankel Operators Spectral properties of Toeplitz+Hankel Operators Torsten Ehrhardt University of California, Santa Cruz ICM Satellite Conference on Operator Algebras and Applications Cheongpung, Aug. 8-12, 2014 Overview

More information

v = 1(1 t) + 1(1 + t). We may consider these components as vectors in R n and R m : w 1. R n, w m

v = 1(1 t) + 1(1 + t). We may consider these components as vectors in R n and R m : w 1. R n, w m 20 Diagonalization Let V and W be vector spaces, with bases S = {e 1,, e n } and T = {f 1,, f m } respectively Since these are bases, there exist constants v i and w such that any vectors v V and w W can

More information

Computation of Rational Szegő-Lobatto Quadrature Formulas

Computation of Rational Szegő-Lobatto Quadrature Formulas Computation of Rational Szegő-Lobatto Quadrature Formulas Joint work with: A. Bultheel (Belgium) E. Hendriksen (The Netherlands) O. Njastad (Norway) P. GONZÁLEZ-VERA Departament of Mathematical Analysis.

More information

2 Eigenvectors and Eigenvalues in abstract spaces.

2 Eigenvectors and Eigenvalues in abstract spaces. MA322 Sathaye Notes on Eigenvalues Spring 27 Introduction In these notes, we start with the definition of eigenvectors in abstract vector spaces and follow with the more common definition of eigenvectors

More information

2. The Concept of Convergence: Ultrafilters and Nets

2. The Concept of Convergence: Ultrafilters and Nets 2. The Concept of Convergence: Ultrafilters and Nets NOTE: AS OF 2008, SOME OF THIS STUFF IS A BIT OUT- DATED AND HAS A FEW TYPOS. I WILL REVISE THIS MATE- RIAL SOMETIME. In this lecture we discuss two

More information

List of Symbols, Notations and Data

List of Symbols, Notations and Data List of Symbols, Notations and Data, : Binomial distribution with trials and success probability ; 1,2, and 0, 1, : Uniform distribution on the interval,,, : Normal distribution with mean and variance,,,

More information

On prediction and density estimation Peter McCullagh University of Chicago December 2004

On prediction and density estimation Peter McCullagh University of Chicago December 2004 On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating

More information

Preconditioners for ill conditioned (block) Toeplitz systems: facts a

Preconditioners for ill conditioned (block) Toeplitz systems: facts a Preconditioners for ill conditioned (block) Toeplitz systems: facts and ideas Department of Informatics, Athens University of Economics and Business, Athens, Greece. Email:pvassal@aueb.gr, pvassal@uoi.gr

More information

A Generalized Eigenmode Algorithm for Reducible Regular Matrices over the Max-Plus Algebra

A Generalized Eigenmode Algorithm for Reducible Regular Matrices over the Max-Plus Algebra International Mathematical Forum, 4, 2009, no. 24, 1157-1171 A Generalized Eigenmode Algorithm for Reducible Regular Matrices over the Max-Plus Algebra Zvi Retchkiman Königsberg Instituto Politécnico Nacional,

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 43 9) 57 69 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa On the best constants in inequalities

More information

One Szegő Widom Limit Theorem

One Szegő Widom Limit Theorem One Szegő Widom Limit Theorem Vladimir A. Vasil ev, Egor A. Maksimenko, Igor B. Simonenko Doklady Mathematics 68:3 (2003), 361 362. Translated from Doklady Akademii Nauk 393:3 (2003), 307 308. English

More information

A Lower Bound for the Size of Syntactically Multilinear Arithmetic Circuits

A Lower Bound for the Size of Syntactically Multilinear Arithmetic Circuits A Lower Bound for the Size of Syntactically Multilinear Arithmetic Circuits Ran Raz Amir Shpilka Amir Yehudayoff Abstract We construct an explicit polynomial f(x 1,..., x n ), with coefficients in {0,

More information

Lecture 8 : Eigenvalues and Eigenvectors

Lecture 8 : Eigenvalues and Eigenvectors CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with

More information

First-order trace formulas for the iterates of the Fox Li operator

First-order trace formulas for the iterates of the Fox Li operator First-order trace formulas for the iterates of the Fox Li operator Albrecht Böttcher, Sergei Grudsky, Daan Huybrechs and Arieh Iserles To the memory of Israel Gohberg, one of the pioneers of present day

More information

k is a product of elementary matrices.

k is a product of elementary matrices. Mathematics, Spring Lecture (Wilson) Final Eam May, ANSWERS Problem (5 points) (a) There are three kinds of elementary row operations and associated elementary matrices. Describe what each kind of operation

More information

INTERLACING PROPERTIES FOR HERMITIAN MATRICES WHOSE GRAPH IS A GIVEN TREE

INTERLACING PROPERTIES FOR HERMITIAN MATRICES WHOSE GRAPH IS A GIVEN TREE INTERLACING PROPERTIES FOR HERMITIAN MATRICES WHOSE GRAPH IS A GIVEN TREE C M DA FONSECA Abstract We extend some interlacing properties of the eigenvalues of tridiagonal matrices to Hermitian matrices

More information

Fredholm determinant with the confluent hypergeometric kernel

Fredholm determinant with the confluent hypergeometric kernel Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2

More information

Applied Linear Algebra in Geoscience Using MATLAB

Applied 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 information

The Singular Values of the Exponentiated Adjacency Matrices of Broom-Tree Graphs

The Singular Values of the Exponentiated Adjacency Matrices of Broom-Tree Graphs The Singular Values of the Exponentiated Adjacency Matrices of Broom-Tree Graphs Tracy Powell Lesley A. Ward, Advisor Estelle Basor, Reader May, 2006 Department of Mathematics Copyright c 2006 Tracy Powell.

More information

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1)

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1) Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1) Travis Schedler Thurs, Nov 17, 2011 (version: Thurs, Nov 17, 1:00 PM) Goals (2) Polar decomposition

More information

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1)

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1) Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1) Travis Schedler Thurs, Nov 17, 2011 (version: Thurs, Nov 17, 1:00 PM) Goals (2) Polar decomposition

More information

Given the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w.

Given the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w. Selected problems from the tetbook J. Neustupa, S. Kračmar: Sbírka příkladů z Matematiky I Problems in Mathematics I I. LINEAR ALGEBRA I.. Vectors, vector spaces Given the vectors u, v, w and real numbers

More information

Matrix structures and applications Luminy, November, Dario A. Bini

Matrix structures and applications Luminy, November, Dario A. Bini Matrix structures and applications Luminy, November, 3 7 2014 Dario A Bini September 30, 2014 Introduction Matrix structures are encountered in various guises in a wide variety of mathematical problems

More information

Iterative solvers for linear equations

Iterative solvers for linear equations Spectral Graph Theory Lecture 23 Iterative solvers for linear equations Daniel A. Spielman November 26, 2018 23.1 Overview In this and the next lecture, I will discuss iterative algorithms for solving

More information

Week 5: Functions and graphs

Week 5: Functions and graphs Calculus and Linear Algebra for Biomedical Engineering Week 5: Functions and graphs H. Führ, Lehrstuhl A für Mathematik, RWTH Aachen, WS 07 Motivation: Measurements at fixed intervals 1 Consider a sequence

More information

Sparse Legendre expansions via l 1 minimization

Sparse Legendre expansions via l 1 minimization Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery

More information

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient

More information

Toeplitz matrices: spectral properties and preconditioning in the CG method

Toeplitz matrices: spectral properties and preconditioning in the CG method Toeplitz matrices: spectral properties and preconditioning in the CG method Stefano Serra-Capizzano December 5, 2014 Abstract We consider multilevel Toeplitz matrices T n (f) generated by Lebesgue integrable

More information

linearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice

linearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice 3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is

More information

Homework 2 Foundations of Computational Math 2 Spring 2019

Homework 2 Foundations of Computational Math 2 Spring 2019 Homework 2 Foundations of Computational Math 2 Spring 2019 Problem 2.1 (2.1.a) Suppose (v 1,λ 1 )and(v 2,λ 2 ) are eigenpairs for a matrix A C n n. Show that if λ 1 λ 2 then v 1 and v 2 are linearly independent.

More information

INTRODUCTION TO LIE ALGEBRAS. LECTURE 10.

INTRODUCTION TO LIE ALGEBRAS. LECTURE 10. INTRODUCTION TO LIE ALGEBRAS. LECTURE 10. 10. Jordan decomposition: theme with variations 10.1. Recall that f End(V ) is semisimple if f is diagonalizable (over the algebraic closure of the base field).

More information

New York Journal of Mathematics. Spectral Approximation of Multiplication Operators

New York Journal of Mathematics. Spectral Approximation of Multiplication Operators New York Journal of Mathematics New York J. Math. 1 (1995) 75 96. Spectral Approximation of Multiplication Operators Kent E. Morrison Abstract. A linear operator on a Hilbert space may be approximated

More information

Bergman shift operators for Jordan domains

Bergman shift operators for Jordan domains [ 1 ] University of Cyprus Bergman shift operators for Jordan domains Nikos Stylianopoulos University of Cyprus Num An 2012 Fifth Conference in Numerical Analysis University of Ioannina September 2012

More information

Random Toeplitz Matrices

Random Toeplitz Matrices Arnab Sen University of Minnesota Conference on Limits Theorems in Probability, IISc January 11, 2013 Joint work with Bálint Virág What are Toeplitz matrices? a0 a 1 a 2... a1 a0 a 1... a2 a1 a0... a (n

More information

Large Deviations for Random Matrices and a Conjecture of Lukic

Large Deviations for Random Matrices and a Conjecture of Lukic Large Deviations for Random Matrices and a Conjecture of Lukic Jonathan Breuer Hebrew University of Jerusalem Joint work with B. Simon (Caltech) and O. Zeitouni (The Weizmann Institute) Western States

More information

Economics 205, Fall 2002: Final Examination, Possible Answers

Economics 205, Fall 2002: Final Examination, Possible Answers Economics 05, Fall 00: Final Examination, Possible Answers Comments on the Exam Grades: 43 possible; high: 413; median: 34; low: 36 I was generally happy with the answers to questions 3-8, satisfied with

More information

Review of some mathematical tools

Review of some mathematical tools MATHEMATICAL FOUNDATIONS OF SIGNAL PROCESSING Fall 2016 Benjamín Béjar Haro, Mihailo Kolundžija, Reza Parhizkar, Adam Scholefield Teaching assistants: Golnoosh Elhami, Hanjie Pan Review of some mathematical

More information

Engineering Jump Start - Summer 2012 Mathematics Worksheet #4. 1. What is a function? Is it just a fancy name for an algebraic expression, such as

Engineering Jump Start - Summer 2012 Mathematics Worksheet #4. 1. What is a function? Is it just a fancy name for an algebraic expression, such as Function Topics Concept questions: 1. What is a function? Is it just a fancy name for an algebraic expression, such as x 2 + sin x + 5? 2. Why is it technically wrong to refer to the function f(x), for

More information

Journées Nationales de Calcul Formel

Journées Nationales de Calcul Formel Journées Nationales de Calcul Formel Rencontre organisée par : Paola Boito, Alin Bostan, Adrien Poteaux et Mohab Safey el Din 2014 Matrix structures and applications Vol 4, n o 1 (2014), Course n o I,

More information

sources of program support: MRSEC Toeplitz at Brandeis April

sources of program support: MRSEC Toeplitz at Brandeis April Landau Quasiparticle Methods for Calculating Eigenvalues and Eigenvectors of Large Matrices Hui Dai, Zachary Geary, and Leo Kadanoff # # e-mail LeoP@UChicago.edu Matrices are used to describe the modes

More information

DIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix

DIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix DIAGONALIZATION Definition We say that a matrix A of size n n is diagonalizable if there is a basis of R n consisting of eigenvectors of A ie if there are n linearly independent vectors v v n such that

More information

AN ALGORITHM FOR COMPUTING FUNDAMENTAL SOLUTIONS OF DIFFERENCE OPERATORS

AN ALGORITHM FOR COMPUTING FUNDAMENTAL SOLUTIONS OF DIFFERENCE OPERATORS AN ALGORITHM FOR COMPUTING FUNDAMENTAL SOLUTIONS OF DIFFERENCE OPERATORS HENRIK BRANDÉN, AND PER SUNDQVIST Abstract We propose an FFT-based algorithm for computing fundamental solutions of difference operators

More information

5.4 Continuity: Preliminary Notions

5.4 Continuity: Preliminary Notions 5.4. CONTINUITY: PRELIMINARY NOTIONS 181 5.4 Continuity: Preliminary Notions 5.4.1 Definitions The American Heritage Dictionary of the English Language defines continuity as an uninterrupted succession,

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Journal of Inequalities in Pure and Applied Mathematics NEWER APPLICATIONS OF GENERALIZED MONOTONE SEQUENCES L. LEINDLER Bolyai Institute, University of Szeged, Aradi vértanúk tere 1 H-6720 Szeged, Hungary

More information

TRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES

TRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES TRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES JOSEPH A. CIMA, WILLIAM T. ROSS, AND WARREN R. WOGEN Abstract. In this paper, we study the matrix representations of compressions of Toeplitz operators

More information

SPECTRAHEDRA. Bernd Sturmfels UC Berkeley

SPECTRAHEDRA. Bernd Sturmfels UC Berkeley SPECTRAHEDRA Bernd Sturmfels UC Berkeley Mathematics Colloquium, North Carolina State University February 5, 2010 Positive Semidefinite Matrices For a real symmetric n n-matrix A the following are equivalent:

More information

The spectral decomposition of near-toeplitz tridiagonal matrices

The spectral decomposition of near-toeplitz tridiagonal matrices Issue 4, Volume 7, 2013 115 The spectral decomposition of near-toeplitz tridiagonal matrices Nuo Shen, Zhaolin Jiang and Juan Li Abstract Some properties of near-toeplitz tridiagonal matrices with specific

More information

Commutative Banach algebras 79

Commutative Banach algebras 79 8. Commutative Banach algebras In this chapter, we analyze commutative Banach algebras in greater detail. So we always assume that xy = yx for all x, y A here. Definition 8.1. Let A be a (commutative)

More information

Applications and fundamental results on random Vandermon

Applications and fundamental results on random Vandermon Applications and fundamental results on random Vandermonde matrices May 2008 Some important concepts from classical probability Random variables are functions (i.e. they commute w.r.t. multiplication)

More information

Mathematical Methods for Physics and Engineering

Mathematical 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 information

A Faber-Krahn type inequality for regular trees

A Faber-Krahn type inequality for regular trees A Faber-Krahn type inequality for regular trees Josef Leydold leydold@statistikwu-wienacat Institut für Statistik, WU Wien Leydold 2003/02/13 A Faber-Krahn type inequality for regular trees p0/36 The Faber-Krahn

More information

Reconstruction of sparse Legendre and Gegenbauer expansions

Reconstruction of sparse Legendre and Gegenbauer expansions Reconstruction of sparse Legendre and Gegenbauer expansions Daniel Potts Manfred Tasche We present a new deterministic algorithm for the reconstruction of sparse Legendre expansions from a small number

More information

EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants

EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For

More information

Iyad T. Abu-Jeib and Thomas S. Shores

Iyad T. Abu-Jeib and Thomas S. Shores NEW ZEALAND JOURNAL OF MATHEMATICS Volume 3 (003, 9 ON PROPERTIES OF MATRIX I ( OF SINC METHODS Iyad T. Abu-Jeib and Thomas S. Shores (Received February 00 Abstract. In this paper, we study the determinant

More information

Spring 2019 Exam 2 3/27/19 Time Limit: / Problem Points Score. Total: 280

Spring 2019 Exam 2 3/27/19 Time Limit: / Problem Points Score. Total: 280 Math 307 Spring 2019 Exam 2 3/27/19 Time Limit: / Name (Print): Problem Points Score 1 15 2 20 3 35 4 30 5 10 6 20 7 20 8 20 9 20 10 20 11 10 12 10 13 10 14 10 15 10 16 10 17 10 Total: 280 Math 307 Exam

More information

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons

Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Magda Khalile (joint work with Konstantin Pankrashkin) Université Paris Sud /25 Robin Laplacians on infinite sectors 1 /

More information

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th

More information

Toeplitz band matrices with exponentially growing condition numbers

Toeplitz band matrices with exponentially growing condition numbers Electronic Journal of Linear Algebra Volume 5 ELA Volume 5 (1999) Article 8 1999 Toeplitz band matrices with exponentially growing condition numbers Albrecht Boettcher aboettch.mathe@mathematik.tu-chemnitz.de

More information

MA 1: PROBLEM SET NO. 7 SOLUTIONS

MA 1: PROBLEM SET NO. 7 SOLUTIONS MA 1: PROBLEM SET NO. 7 SOLUTIONS 1. (1) Obtain the number r = 15 3 as an approimation to the nonzero root of the equation 2 = sin() by using the cubic Taylor polynomial approimation to sin(). Proof. The

More information

The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In

The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed

More information

Linear Algebra. Session 12

Linear Algebra. Session 12 Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)

More information

arxiv: v2 [math.fa] 8 Jan 2014

arxiv: v2 [math.fa] 8 Jan 2014 A CLASS OF TOEPLITZ OPERATORS WITH HYPERCYCLIC SUBSPACES ANDREI LISHANSKII arxiv:1309.7627v2 [math.fa] 8 Jan 2014 Abstract. We use a theorem by Gonzalez, Leon-Saavedra and Montes-Rodriguez to construct

More information

Rolle s Theorem for Polynomials of Degree Four in a Hilbert Space 1

Rolle s Theorem for Polynomials of Degree Four in a Hilbert Space 1 Journal of Mathematical Analysis and Applications 265, 322 33 (2002) doi:0.006/jmaa.200.7708, available online at http://www.idealibrary.com on Rolle s Theorem for Polynomials of Degree Four in a Hilbert

More information

g 2 (x) (1/3)M 1 = (1/3)(2/3)M.

g 2 (x) (1/3)M 1 = (1/3)(2/3)M. COMPACTNESS If C R n is closed and bounded, then by B-W it is sequentially compact: any sequence of points in C has a subsequence converging to a point in C Conversely, any sequentially compact C R n is

More information

Eigenvalue PDFs. Peter Forrester, M&S, University of Melbourne

Eigenvalue PDFs. Peter Forrester, M&S, University of Melbourne Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices

More information

Sign changes of Fourier coefficients of cusp forms supported on prime power indices

Sign changes of Fourier coefficients of cusp forms supported on prime power indices Sign changes of Fourier coefficients of cusp forms supported on prime power indices Winfried Kohnen Mathematisches Institut Universität Heidelberg D-69120 Heidelberg, Germany E-mail: winfried@mathi.uni-heidelberg.de

More information

Szegö Limit Theorems With Varying Coefficients. Optimal trace ideal properties of the restriction operator and applications

Szegö Limit Theorems With Varying Coefficients. Optimal trace ideal properties of the restriction operator and applications Szegö Limit Theorems With Varying Coefficients Alain Bourget California State University, Fullerton, abourget@fullerton.edu Szegö s first limit theorem is a remarkable result concerning the asymptotic

More information

Eigenvalue and Eigenvector Homework

Eigenvalue and Eigenvector Homework Eigenvalue and Eigenvector Homework Olena Bormashenko November 4, 2 For each of the matrices A below, do the following:. Find the characteristic polynomial of A, and use it to find all the eigenvalues

More information

REPRESENTATION THEORY WEEK 7

REPRESENTATION THEORY WEEK 7 REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable

More information

Convex Functions and Optimization

Convex Functions and Optimization Chapter 5 Convex Functions and Optimization 5.1 Convex Functions Our next topic is that of convex functions. Again, we will concentrate on the context of a map f : R n R although the situation can be generalized

More information

Singular Value Decomposition

Singular Value Decomposition The Bulletin of Society for Mathematical Services and Standards Online: 2014-09-01 ISSN: 2277-8020, Vol. 11, pp 13-20 doi:10.18052/www.scipress.com/bsmass.11.13 2014 SciPress Ltd., Switzerland Singular

More information

Sample Quantum Chemistry Exam 2 Solutions

Sample Quantum Chemistry Exam 2 Solutions Chemistry 46 Fall 7 Dr. Jean M. Standard Name SAMPE EXAM Sample Quantum Chemistry Exam Solutions.) ( points) Answer the following questions by selecting the correct answer from the choices provided. a.)

More information

City Suburbs. : population distribution after m years

City Suburbs. : population distribution after m years Section 5.3 Diagonalization of Matrices Definition Example: stochastic matrix To City Suburbs From City Suburbs.85.03 = A.15.97 City.15.85 Suburbs.97.03 probability matrix of a sample person s residence

More information

Vector fields on the plane

Vector fields on the plane Vector fields on the plane Sasha Patotski Cornell University ap744@cornell.edu December 15, 2014 Sasha Patotski (Cornell University) Vector fields on the plane December 15, 2014 1 / 11 Vector fields Definition

More information

9 Radon-Nikodym theorem and conditioning

9 Radon-Nikodym theorem and conditioning Tel Aviv University, 2015 Functions of real variables 93 9 Radon-Nikodym theorem and conditioning 9a Borel-Kolmogorov paradox............. 93 9b Radon-Nikodym theorem.............. 94 9c Conditioning.....................

More information

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues

More information

T.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003

T.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003 T.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003 A vector x R n is called positive, symbolically x > 0,

More information

Determinants of Perturbations of Finite Toeplitz Matrices. Estelle Basor American Institute of Mathematics

Determinants of Perturbations of Finite Toeplitz Matrices. Estelle Basor American Institute of Mathematics Determinants of Perturbations of Finite Toeplitz Matrices Estelle Basor American Institute of Mathematics October 2010 One of the purposes of this talk is to describe a Helton observation that fundamentally

More information

Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks

Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks Course Description - Master in of Mathematics Comprehensive exam& Thesis Tracks 1309701 Theory of ordinary differential equations Review of ODEs, existence and uniqueness of solutions for ODEs, existence

More information

The Geometry of Semidefinite Programming. Bernd Sturmfels UC Berkeley

The Geometry of Semidefinite Programming. Bernd Sturmfels UC Berkeley The Geometry of Semidefinite Programming Bernd Sturmfels UC Berkeley Positive Semidefinite Matrices For a real symmetric n n-matrix A the following are equivalent: All n eigenvalues of A are positive real

More information

Modeling and Stability Analysis of a Communication Network System

Modeling and Stability Analysis of a Communication Network System Modeling and Stability Analysis of a Communication Network System Zvi Retchkiman Königsberg Instituto Politecnico Nacional e-mail: mzvi@cic.ipn.mx Abstract In this work, the modeling and stability problem

More information

BOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS

BOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS BOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS Alessandra Luati Dip. Scienze Statistiche University of Bologna Tommaso Proietti S.E.F. e ME. Q. University

More information

Linear Algebra, part 3 QR and SVD

Linear Algebra, part 3 QR and SVD Linear Algebra, part 3 QR and SVD Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2012 Going back to least squares (Section 1.4 from Strang, now also see section 5.2). We

More information

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero. Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called

More information