Spectral analysis of two doubly infinite Jacobi operators
|
|
- Aron Walsh
- 5 years ago
- Views:
Transcription
1 Spectral analysis of two doubly infinite Jacobi operators František Štampach jointly with Mourad E. H. Ismail Stockholm University Spectral Theory and Applications conference in memory of Boris Pavlov March 13, 2016 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
2 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
3 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
4 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
5 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
6 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
7 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. B with n 0 [Ismail, Mulla 87]: q-chebyshev polynomials. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
8 Two doubly-infinite Jacobi matrices We analyze the spectral properties of Jacobi operators A and B acting on vectors of the standard basis of l 2 Z) as: and Ae n = q n+1 e n 1 + q n e n+1, n Z, where q 0, 1) and α R. Be n = e n 1 + αq n e n + e n+1, n Z, The semi-infinite parts of A and B has been studied before mainly within the theory of Orthogonal Polynomials: A with n < 0 [Al-Salam, Ismail 83]: Rogers Ramanujan s functions and continued fractions. A with n > 0 [Chen, Ismail 98, FS 16]: Indeterminate moment problem, Nevanlinna functions. B with n 0 [Ismail, Mulla 87]: q-chebyshev polynomials. The spectrum of any associated semi-infinite Jacobi operator is never known explicitly α 0) but is expressible in terms of zeros of certain special functions.functions František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
9 Special functions 1/2 - basic hypergeometric series Let 0 < q < 1, r, s Z +. Recall the basic hypergeometric function [ ] a1, a r φ 2,... a r s ; q, z b 1, b 2,... b s is defined by the power series a 1 ; q) na 2 ; q) n... a r ; q) n 1) s r+1)n q s r+1)nn 1)/2 z n b n=0 1 ; q) nb 2 ; q) n... b s; q) n q; q) n where z, a 1, a 2,..., a r C, b 1, b 2,..., b s C \ q Z and a; q) 0 = 1, a; q) n = 1 a)1 aq)... 1 aq n 1 ) is the q-pochhammer symbol. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
10 Special functions 1/2 - basic hypergeometric series Let 0 < q < 1, r, s Z +. Recall the basic hypergeometric function [ ] a1, a r φ 2,... a r s ; q, z b 1, b 2,... b s is defined by the power series a 1 ; q) na 2 ; q) n... a r ; q) n 1) s r+1)n q s r+1)nn 1)/2 z n b n=0 1 ; q) nb 2 ; q) n... b s; q) n q; q) n where z, a 1, a 2,..., a r C, b 1, b 2,..., b s C \ q Z and a; q) 0 = 1, a; q) n = 1 a)1 aq)... 1 aq n 1 ) is the q-pochhammer symbol. Here we will need only 0 φ 1 and 1 φ 1. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
11 Special functions 2/2 - theta functions The theta function: 1 θ qz) := z; q) q/z; q) = q nn 1)/2 z) n q; q) n= František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
12 Special functions 2/2 - theta functions The theta function: 1 θ qz) := z; q) q/z; q) = q nn 1)/2 z) n q; q) n= Jacobi s theta functions: ϑ 1 z q) = iq 1/4 e iz q 2 ; q 2 ) θ q 2 e 2iz) ϑ 2 z q) = q 1/4 e iz q 2 ; q 2 ) θ q 2 e 2iz) ϑ 3 z q) = q 2 ; q 2 ) θ q 2 qe 2iz) ϑ 4 z q) = q 2 ; q 2 ) θ q 2 qe 2iz) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
13 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
14 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as has deficiency indices 1, 1). Ae n = q n+1 e n 1 + q n e n+1 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
15 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
16 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
17 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: Proposition self-adjoint extensions) For t R }, operators A t, acting as A t ψ = Aψ, with domains Dom A t = ψ D if t R, or are all self-adjoint extensions of A. lim n q n ψ 2n+1 + tψ 2n ) = 0 Dom A = ψ D } lim n q n qψ 2n 1 tψ 2n ) = 0, } lim n q n ψ 2n = 0, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
18 Operators associated with the Jacobi matrix A Operator A with Dom A = spane n n Z} acting as Ae n = q n+1 e n 1 + q n e n+1 has deficiency indices 1, 1). Thus, there is a one-parameter family of self-adjoint extensions A t, t R }. Let D := ψ l 2 Aψ l 2 }. By using the theory of self-adjoint extensions and simple structure of matrix A one gets: Proposition self-adjoint extensions) For t R }, operators A t, acting as A t ψ = Aψ, with domains Dom A t = ψ D if t R, or are all self-adjoint extensions of A. In addition, lim n q n ψ 2n+1 + tψ 2n ) = 0 Dom A = ψ D } lim n q n qψ 2n 1 tψ 2n ) = 0, } lim n q n ψ 2n = 0, σ ca t ) = σessa t ) = 0}, t R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
19 Solutions of the eigenvalue equation The q-exponential function: E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
20 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
21 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. By inspection of the asymptotic behavior of ψ ± n, as n ±, one gets: ψ ± l 2 + ), however, ψ ± / l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
22 Solutions of the eigenvalue equation The q-exponential function: Sequences ψ ±, where E qz) = n=0 q n2 /4 ) z n = 1 φ 1 0; q 1/2 ; q 1/2, q 1/4 z q; q) n ψ n ± := ±i) n q n/2 E q 2 ±ixq n ), are two linearly independent solutions of the difference equation Aψ = xψ. By inspection of the asymptotic behavior of ψ ± n, as n ±, one gets: ψ ± l 2 + ), however, ψ ± / l 2 ). Hence, one expects there are non-trivial coefficients a = ax) and b = bx) such that aψ + + bψ l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
23 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
24 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. Moreover, ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2) and František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
25 The l 2 -solution Proposition For all x C \ 0}, the sequence ) ) ϕx) := θ q iq 1/2 x ψ ) x) + θ q iq 1/2 x ψ +) x), is the non-trivial solution of Aφ = xφ which belongs to l 2 Z). In addition, within the space l 2 Z), this solution is given uniquely up to a multiplicative constant. Moreover, and ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2) q 2 ϕx) 2 ; q 2) 2 l 2 = 4 ) θ q; q 2 2 q 2 z 2) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
26 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
27 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
28 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. How to solve the secular equation in general? František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
29 The secular equation Theorem secular equation) For t R }, one has spec c A t ) = 0} and spec p A t ) coincides with the set of roots of the secular equation: xθ q 4 q 2 x 2) + tθ q 4 x 2) = 0, for t R, and θ q 4 x 2) = 0, for t =. In addition, all eigenvalues of A t are simple. Corollary: } spec p A 0 ) = ±q 2n+1 n Z } and spec p A ) = ±q 2n n Z. How to solve the secular equation in general? Reparametrize t = Φs) and use nice properties of Jacobi s theta functions... František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
30 Reparametrization t = Φs) The function Φs) := iq 1/2 ϑ 4 is q 2 ) ϑ ) 1 is q 2 is real-valued, strictly decreasing on 0, 2 ln q), and maps [0, 2 ln q) onto R }. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
31 Reparametrization t = Φs) The function Φs) := iq 1/2 ϑ 4 is q 2 ) ϑ ) 1 is q 2 is real-valued, strictly decreasing on 0, 2 ln q), and maps [0, 2 ln q) onto R }. For the inverse function, one has Φ 1 dx t) = Cq) Dq) ) t + x 2 q + x 2) where and Cq) = q 1/2 ϑ 2 0 q 2 ) ϑ 3 0 q 2 ) 0 q 2 ) Dq) = qϑ2 3 ϑ 2 ). 2 0 q 2 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
32 Spectrum fully explicitly Using the reparametrization t = Φs), the secular equation simplifies to θ q 2 e s x ) θ q 2 e s x ) = 0. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
33 Spectrum fully explicitly Using the reparametrization t = Φs), the secular equation simplifies to θ q 2 e s x ) θ q 2 e s x ) = 0. Theorem Let t R }, then where spec p A t ) = e s q 2Z e s q 2Z dx s = Cq) Dq) ) t + x 2 q + x 2). In addition, the family of corresponding eigenvectors ϕ±e s q 2N ) N Z}, where forms an orthogonal basis of l 2 Z). ϕ nx) = 1; q) x n q nn 1)/2 0φ 1 ; 0; q 2, q 2n+4 x 2), František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
34 Contents 1 Introduction 2 Spectral resolution of A 3 Spectral resolution of B František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
35 The discrete Schrödinger operator B The second Jacobi matrix B determines the unique operator B = U + U + αv where U is the forward shift operator and V is the self-adjoint diagonal operator: Ue n = e n+1 and Ve n = q n e n, n Z. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
36 The discrete Schrödinger operator B The second Jacobi matrix B determines the unique operator B = U + U + αv where U is the forward shift operator and V is the self-adjoint diagonal operator: Ue n = e n+1 and Ve n = q n e n, n Z. Proposition essential spectrum) The operator B is self-adjoint and one has σessb) = [ 2, 2]. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
37 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Bψ = z + z 1 )ψ. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
38 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Note that z z + z 1 : Bψ = z + z 1 )ψ. z 0 < z < 1} C \ [ 2, 2], outside σessa) ) e iθ θ [0, π]} [ 2, 2], inside σessa) ) František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
39 Solutions of the eigenvalue equation Using properties of the Hahn-Exton q-bessel functions one verifies: f nz) := 1) n α n q 2 1 nn+1) z 1 α 1 q n+1 ; q ) 1 φ 1 0; z 1 α 1 q n+1 ; q, zα 1 q n+1) and g nz) := z n zαq 1 n ; q ) 1 φ 1 0; zαq 1 n ; q, qz 2). are two solutions of the equation for all α, z 0. Note that z z + z 1 : Bψ = z + z 1 )ψ. z 0 < z < 1} C \ [ 2, 2], outside σessa) ) e iθ θ [0, π]} [ 2, 2], inside σessa) ) The solutions f z) and gz) are linearly independent iff z / α 1 q Z 0} since W f, g) = z 1 θ q αz). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
40 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
41 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z).. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
42 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B.. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
43 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B. Further, eigenvectors v m corresponding to eigenvalues α 1 q m + αq m, can be chosen as v m = v m,j } j=, with v m,j = f j α 1 q m) = 1) j α j q 1 2 jj+1) q m+j+1 ; q) 1 φ 1 0; q m+j+1 ; q, α 2 q m+j+1).. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
44 The point spectrum Detailed asymptotic analysis of solutions f and g yields: If 0 < z < 1 and z / α 1 q Z 0}, then l 2 + ) / l 2 + ) f z) / l 2 gz) ) l 2 ) If z = 1 the asymptotic behavior of solutions is very different and, in the end, it implies that: For α R and x [ 2, 2], there is no non-trivial solution of Bψ = xψ belonging to l 2 Z). Theorem point spectrum) If α 0, then } σb) \ [ 2, 2] = σ pb) = α 1 q m + αq m m > log q α and all points from this set are simple eigenvalues of B. Further, eigenvectors v m corresponding to eigenvalues α 1 q m + αq m, can be chosen as v m = v m,j } j=, with v m,j = f j α 1 q m) = 1) j α j q 1 2 jj+1) q m+j+1 ; q) 1 φ 1 0; q m+j+1 ; q, α 2 q m+j+1).. In addition, v m l 2 Z) = α m q mm+1)/2 1 α 2 q 2m q; q), m > log q α. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
45 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
46 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. To determine the spectral measure in the essential spectrum we use the formula E k,l a, b)) = lim lim 1 b δ Gk,l x + iɛ) G k,l x iɛ) ) dx, δ 0+ ɛ 0+ 2πi a+δ where G k,l z) := e k, B z) 1 e l = 1 g k z)f l z), k l, W f, g) g l z)f k z), k l, František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
47 The absolutely continuous part of the spectral measure Let us denote E k,l ) := e k, E B )e l, k, l Z, where E B stands for the spectral measure of B. To determine the spectral measure in the essential spectrum we use the formula E k,l a, b)) = lim lim 1 b δ Gk,l x + iɛ) G k,l x iɛ) ) dx, δ 0+ ɛ 0+ 2πi a+δ where G k,l z) := e k, B z) 1 e l = 1 g k z)f l z), k l, W f, g) g l z)f k z), k l, Proposition Let α 0 and 2 a < b 2. Then for any k, l Z, it holds E k,l [a, b]) = 1 2π φa f l φ b e iφ) f k e iφ) e 2iφ ; q ) 2 αe iφ, qα 1 e iφ ; q ) dφ where φ a = arccos a/2) and φ b = arccos b/2). Consequently, σ acb) = [ 2, 2]. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
48 Summary Theorem If α 0, then and σ essb) = σ acb) = [ 2, 2], } σ pb) = σ d B) = α 1 q m + αq m m > log q α. František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
49 Summary Theorem If α 0, then and σ essb) = σ acb) = [ 2, 2], } σ pb) = σ d B) = α 1 q m + αq m m > log q α. In addition, for M R a Borel set, we have E k,l M) = 1 f l e iφ) f k e iφ) e 2iφ ; q ) 2 2π 2 cos φ [ 2,2] A αe iφ, qα 1 e iφ ; q ) dφ α 2 q 2m) α 2m q mm+1) f l α 1 q m) f k α 1 q m). q; q) 2 m> log α α 1 q m +αq m M František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
50 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
51 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). Elements of the eigenvectors v m are expressible in terms of Hanhn-Exton q-bessel function and the formula for the norm of the eigenvectors yields n Z Jn 2 z; q) = 1, z < 1. 1 z2 František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
52 One consequence for special functions Recall the Hanhn-Exton or third Jackson s) q-bessel function is defined as J νz; q) = z ν qν+1 ; q) q; q) 1φ 1 0; q ν+1 ; q, qz 2). Elements of the eigenvectors v m are expressible in terms of Hanhn-Exton q-bessel function and the formula for the norm of the eigenvectors yields n Z Jn 2 z; q) = 1, z < 1. 1 z2 This formula seems to be new Really?) and it generalizes the well-known summation formula for the Bessel functions of the first kind: Jn 2 z) = 1, z < 1. n Z František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
53 The end Thank you! František Štampach Stockholm University) Jacobi operators Spectr. Theor. Appl.) March 13, / 20
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationOn rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro
On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational
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 informationThe roots of The Third Jackson q-bessel Function
The roots of The Third Jacson q-bessel Function L. D. Abreu Departamento de Matemática Universidade de Coimbra Coimbra, Portugal J. Bustoz Department of Mathematics, Arizona State University Tempe, Arizona
More informationKaczmarz algorithm in Hilbert space
STUDIA MATHEMATICA 169 (2) (2005) Kaczmarz algorithm in Hilbert space by Rainis Haller (Tartu) and Ryszard Szwarc (Wrocław) Abstract The aim of the Kaczmarz algorithm is to reconstruct an element in a
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationBIRTH-DEATH PROCESSES AND q-continued FRACTIONS
BIRTH-DEATH PROCESSES AND q-continued FRACTIONS TONY FENG, RACHEL KIRSCH, ELISE MCCALL, AND MATT WAGE Abstract. In the 997 paper of Parthasarathy, Lenin, Schoutens, and Van Assche [6], the authors study
More informationOrthogonal Polynomial Ensembles
Chater 11 Orthogonal Polynomial Ensembles 11.1 Orthogonal Polynomials of Scalar rgument Let wx) be a weight function on a real interval, or the unit circle, or generally on some curve in the comlex lane.
More informationGaussian Random Fields
Gaussian Random Fields Mini-Course by Prof. Voijkan Jaksic Vincent Larochelle, Alexandre Tomberg May 9, 009 Review Defnition.. Let, F, P ) be a probability space. Random variables {X,..., X n } are called
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 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 informationSpectral Analysis of Matrices - An Introduction for Engineers
Spectral Analysis of Matrices - An Introduction for Engineers Timo Weidl Department of Mathematics Royal Institute of Technology S-144 Stockholm weidl@mathkthse January 13, 26 Plan of the Talk 1 A practical
More informationMATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003
MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003 1. True or False (28 points, 2 each) T or F If V is a vector space
More informationOperator algebras related to bounded positive operator. operator with simple spectrum
Operator algebras related to bounded positive operator with simple spectrum Institute of Mathematics, University of Bia lystok joint work with A. Odzijewicz, A. Tereszkiewicz 13 June 2007 Varna Coherent
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationPseudospectra 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 informationSPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT
SPECTRAL THEOREM FOR SYMMETRIC OPERATORS WITH COMPACT RESOLVENT Abstract. These are the letcure notes prepared for the workshop on Functional Analysis and Operator Algebras to be held at NIT-Karnataka,
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationMatrix orthogonal polynomials and group representations
Matrix orthogonal polynomials and group representations Inés Pacharoni Universidad Nacional de Córdoba, Argentina (Joint work with Juan Tirao) Discrete Systems and Special Functions Cambridge, UK. July
More informationJUST THE MATHS UNIT NUMBER 9.9. MATRICES 9 (Modal & spectral matrices) A.J.Hobson
JUST THE MATHS UNIT NUMBER 9.9 MATRICES 9 (Modal & spectral matrices) by A.J.Hobson 9.9. Assumptions and definitions 9.9.2 Diagonalisation of a matrix 9.9.3 Exercises 9.9.4 Answers to exercises UNIT 9.9
More informationMatrix Lie groups. and their Lie algebras. Mahmood Alaghmandan. A project in fulfillment of the requirement for the Lie algebra course
Matrix Lie groups and their Lie algebras Mahmood Alaghmandan A project in fulfillment of the requirement for the Lie algebra course Department of Mathematics and Statistics University of Saskatchewan March
More information13 Maximum Modulus Principle
3 Maximum Modulus Principle Theorem 3. (maximum modulus principle). If f is non-constant and analytic on an open connected set Ω, then there is no point z 0 Ω such that f(z) f(z 0 ) for all z Ω. Remark
More information1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i
Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationAsymptotic distribution of eigenvalues of Laplace operator
Asymptotic distribution of eigenvalues of Laplace operator 23.8.2013 Topics We will talk about: the number of eigenvalues of Laplace operator smaller than some λ as a function of λ asymptotic behaviour
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
More informationEigenvalues and Eigenvectors
Sec. 6.1 Eigenvalues and Eigenvectors Linear transformations L : V V that go from a vector space to itself are often called linear operators. Many linear operators can be understood geometrically by identifying
More informationSymmetric and self-adjoint matrices
Symmetric and self-adjoint matrices A matrix A in M n (F) is called symmetric if A T = A, ie A ij = A ji for each i, j; and self-adjoint if A = A, ie A ij = A ji or each i, j Note for A in M n (R) that
More informationAnalysis IV : Assignment 3 Solutions John Toth, Winter ,...). In particular for every fixed m N the sequence (u (n)
Analysis IV : Assignment 3 Solutions John Toth, Winter 203 Exercise (l 2 (Z), 2 ) is a complete and separable Hilbert space. Proof Let {u (n) } n N be a Cauchy sequence. Say u (n) = (..., u 2, (n) u (n),
More informationDETERMINANT IDENTITIES FOR THETA FUNCTIONS
DETERMINANT IDENTITIES FOR THETA FUNCTIONS SHAUN COOPER AND PEE CHOON TOH Abstract. Two proofs of a theta function identity of R. W. Gosper and R. Schroeppel are given. A cubic analogue is presented, and
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationQuadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.
Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric
More informationKrein s formula and perturbation theory
ISSN: 40-567 Krein s formula and perturbation theory P.Kurasov S.T.Kuroda Research Reports in athematics Number 6, 2000 Department of athematics Stockholm University Electronic versions of this document
More informationA Right Inverse Of The Askey-Wilson Operator arxiv:math/ v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail
A Right Inverse Of The Askey-Wilson Operator arxiv:math/9310219v1 [math.ca] 5 Oct 1993 B. Malcolm Brown and Mourad E. H. Ismail Abstract We establish an integral representation of a right inverse of the
More informationACM 104. Homework Set 4 Solutions February 14, 2001
ACM 04 Homework Set 4 Solutions February 4, 00 Franklin Chapter, Problem 4, page 55 Suppose that we feel that some observations are more important or reliable than others Redefine the function to be minimized
More informationFredholm Theory. April 25, 2018
Fredholm Theory April 25, 208 Roughly speaking, Fredholm theory consists of the study of operators of the form I + A where A is compact. From this point on, we will also refer to I + A as Fredholm operators.
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationTrace Class Operators and Lidskii s Theorem
Trace Class Operators and Lidskii s Theorem Tom Phelan Semester 2 2009 1 Introduction The purpose of this paper is to provide the reader with a self-contained derivation of the celebrated Lidskii Trace
More informationBasic Calculus Review
Basic Calculus Review Lorenzo Rosasco ISML Mod. 2 - Machine Learning Vector Spaces Functionals and Operators (Matrices) Vector Space A vector space is a set V with binary operations +: V V V and : R V
More informationProblem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.
Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)
More informationRIESZ BASES AND UNCONDITIONAL BASES
In this paper we give a brief introduction to adjoint operators on Hilbert spaces and a characterization of the dual space of a Hilbert space. We then introduce the notion of a Riesz basis and give some
More informationReal Variables # 10 : Hilbert Spaces II
randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,
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 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 informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More informationARTICLE IN PRESS. JID:YAIMA AID:6155 /FLA [m1l; v1.231; Prn:21/02/2018; 9:28] P.1 (1-56) Advances in Mathematics ( )
JID:YAIMA AID:655 /FLA [ml; v.3; Prn:/0/08; 9:8] P. -56 Advances in Mathematics Contents lists available at ScienceDirect Advances in Mathematics 5 www.elsevier.com/locate/aim 5 Nevanlinna theory of the
More informationAn identity of Andrews and the Askey-Wilson integral
Ramanujan J DOI 0.007/s39-008-922-4 An identity of Andrews and the Askey-Wilson integral Zhi-Guo Liu Received: 6 July 2007 / Accepted: 7 January 2008 Springer Science+Business Media, LLC 2008 Abstract
More informationHeisenberg s inequality
Heisenberg s inequality Michael Walter 1 Introduction 1 Definition. Let P be a probability measure on R. Then its mean or expectation is defined as E(P) := x dp(x) Its variance is the expected square deviation
More informationCentral Limit Theorems for linear statistics for Biorthogonal Ensembles
Central Limit Theorems for linear statistics for Biorthogonal Ensembles Maurice Duits, Stockholm University Based on joint work with Jonathan Breuer (HUJI) Princeton, April 2, 2014 M. Duits (SU) CLT s
More informationMORE NOTES FOR MATH 823, FALL 2007
MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C
More informationMath Homework 8 (selected problems)
Math 102 - Homework 8 (selected problems) David Lipshutz Problem 1. (Strang, 5.5: #14) In the list below, which classes of matrices contain A and which contain B? 1 1 1 1 A 0 0 1 0 0 0 0 1 and B 1 1 1
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More informationFormalism of quantum mechanics
Formalism of quantum mechanics August 22, 204 Contents Introduction 2 Classical physics 3 Quantum physics 3 4 Harmonic oscillator 0 Introduction Quite a bit of the serious mathematical theory of self-adjoint
More informationInteger Partitions With Even Parts Below Odd Parts and the Mock Theta Functions
Integer Partitions With Even Parts Below Odd Parts and the Mock Theta Functions by George E. Andrews Key Words: Partitions, mock theta functions, crank AMS Classification Numbers: P84, P83, P8, 33D5 Abstract
More information5 Compact linear operators
5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.
More informationBorn-Oppenheimer Approximation in a singular system Prague. O.T.Turgut Bogazici Univ.
Born-Oppenheimer Approximation in a singular system Prague O.T.Turgut Bogazici Univ. June 5, 206 How can we use Born-Oppenheimer in a singular system (if we can)? What is the usual approach: Let us consider
More informationReal Analysis Problems
Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More informationFurther Mathematical Methods (Linear Algebra)
Further Mathematical Methods (Linear Algebra) Solutions For The Examination Question (a) To be an inner product on the real vector space V a function x y which maps vectors x y V to R must be such that:
More informationPhysics 505 Homework No. 1 Solutions S1-1
Physics 505 Homework No s S- Some Preliminaries Assume A and B are Hermitian operators (a) Show that (AB) B A dx φ ABψ dx (A φ) Bψ dx (B (A φ)) ψ dx (B A φ) ψ End (b) Show that AB [A, B]/2+{A, B}/2 where
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 informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationSome Combinatorial and Analytical Identities
Some Combinatorial and Analytical Identities Mourad EH Ismail Dennis Stanton November 6, 200 Abstract We give new proofs and explain the origin of several combinatorial identities of Fu and Lascoux, Dilcher,
More informationComputationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:
Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication
More informationChapter 2 Spectra of Finite Graphs
Chapter 2 Spectra of Finite Graphs 2.1 Characteristic Polynomials Let G = (V, E) be a finite graph on n = V vertices. Numbering the vertices, we write down its adjacency matrix in an explicit form of n
More informationEigenvalues and Eigenvectors A =
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
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 informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationBergman 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 informationThe L p -dissipativity of first order partial differential operators
The L p -dissipativity of first order partial differential operators A. Cialdea V. Maz ya n Memory of Vladimir. Smirnov Abstract. We find necessary and sufficient conditions for the L p -dissipativity
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationCalculus C (ordinary differential equations)
Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear
More informationLINEAR ALGEBRA MICHAEL PENKAVA
LINEAR ALGEBRA MICHAEL PENKAVA 1. Linear Maps Definition 1.1. If V and W are vector spaces over the same field K, then a map λ : V W is called a linear map if it satisfies the two conditions below: (1)
More informationScattering matrices and Weyl functions
Scattering matrices and Weyl functions Jussi Behrndt a, Mark M. Malamud b, Hagen Neidhardt c April 6, 2006 a)technische Universität Berlin, Institut für Mathematik, Straße des 7. Juni 36, D 0623 Berlin,
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationIn this section again we shall assume that the matrix A is m m, real and symmetric.
84 3. The QR algorithm without shifts See Chapter 28 of the textbook In this section again we shall assume that the matrix A is m m, real and symmetric. 3.1. Simultaneous Iterations algorithm Suppose we
More informationRemarks on Some Basic Hypergeometric Series
Remarks on Some Basic Hypergeometric Series Changgui Zhang Many results in Mathematical Analysis seem to come from some obvious computations. For a few years, we have been interested in the analytic theory
More informationCommutative 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 informationMATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)
MATH 20F: LINEAR ALGEBRA LECTURE B00 (T KEMP) Definition 01 If T (x) = Ax is a linear transformation from R n to R m then Nul (T ) = {x R n : T (x) = 0} = Nul (A) Ran (T ) = {Ax R m : x R n } = {b R m
More informationChapter 3. Vector spaces
Chapter 3. Vector spaces Lecture notes for MA1111 P. Karageorgis pete@maths.tcd.ie 1/22 Linear combinations Suppose that v 1,v 2,...,v n and v are vectors in R m. Definition 3.1 Linear combination We say
More informationEigenvalues, Eigenvectors, and an Intro to PCA
Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Changing Basis We ve talked so far about re-writing our data using a new set of variables, or a new basis.
More informationFredholm 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 informationOWN ZEROS LUIS DANIEL ABREU
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 4 3 FUNCTIONS q-orthogonal WITH RESPECT TO THEIR OWN ZEROS LUIS DANIEL ABREU Abstract: In [4], G. H. Hardy proved that,
More information. The following is a 3 3 orthogonal matrix: 2/3 1/3 2/3 2/3 2/3 1/3 1/3 2/3 2/3
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. An n n matrix
More informationExplicit inverse of a tridiagonal (p, r) Toeplitz matrix
Explicit inverse of a tridiagonal (p, r Toeplitz matrix A.M. Encinas, M.J. Jiménez Dpt. Matemàtiques, Universitat Politècnica de Catalunya - BarcelonaTech Abstract Tridiagonal matrices appears in many
More informationProblems 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 informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationLecture 10: Solving the Time-Independent Schrödinger Equation. 1 Stationary States 1. 2 Solving for Energy Eigenstates 3
Contents Lecture 1: Solving the Time-Independent Schrödinger Equation B. Zwiebach March 14, 16 1 Stationary States 1 Solving for Energy Eigenstates 3 3 Free particle on a circle. 6 1 Stationary States
More informationDRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS
DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary
More informationOn Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable
Communications in Mathematics and Applications Volume (0), Numbers -3, pp. 97 09 RGN Publications http://www.rgnpublications.com On Tricomi and Hermite-Tricomi Matrix Functions of Complex Variable A. Shehata
More informationApplications of Diagonalization
Applications of Diagonalization Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw Apr 2, 200 The notes cover applications of matrix diagonalization Boas 3.2. Quadratic curves Consider the quadratic curve, 5x 2 4xy
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More information