d ORTHOGONAL POLYNOMIALS AND THE WEYL ALGEBRA

Size: px
Start display at page:

Download "d ORTHOGONAL POLYNOMIALS AND THE WEYL ALGEBRA"

Transcription

1 d ORTHOGONAL POLYNOMIALS AND THE WEYL ALGEBRA Luc VINET Université de Montréal Alexei ZHEDANOV Donetsk for Physics and Technology

2 PLAN 1. Charlier polynomials 2. Weyl algebra 3. Charlier and Weyl 4. d orthogonal polynomials 5. A non linear automorphism of the Weyl algebra 6. d Charlier polynomials

3 7. difference ladder 8. Explicit expression in special case 9. Dual and orthogonality 10. Conclusion

4 1.CHARLIER POLYNOMIALS p n (k) Polynomials of a discrete variable k (used normalized version) 3 term recurrence rela@on p n +1 (k) = (k n a) p n (k) na p n 1 (k) Orthogonality rela@on k= 0 p m (k) p n (k) e a a k k! = a n n!δ m,n

5 2.WEYL ALGEBRA 3 generators: Rela@ons: a, a *, I [ a, a * ] =1, [ a, I] = [ a *, I] = 0 Canonical irreducible representa@on: a n = n n 1 a * n = n +1 n +1 I n = n a * = a + N = a + a N n Vacuum vector a 0 = 0

6 3.CHARLIER POLYNOMIALS AND WEYL ALGEBRA Oscillator a = 1 2 ( d dx + x) a+ = 1 2 ( d dx + x) Shib N = 1 2 x x λ d 2 dx x2 1 2 N(λ) = N λx + λ2 2 a(λ) = a λ 2 a + (λ) = a + λ 2 x = 1 2 (a + a+ ) = 1 2 (a + a+ ) + λ N = N + λ 2 (a + a+ ) + λ2 2

7 Let N n = n n and set k = n k n n N k = n k N + λ 2 (a + a+ ) + λ2 n 2 n k n k = (n + λ2 n k = 2 λ 2 ) n k + λ 2 n +1 n +1 k + λ 2 n n 1 k n 1 n! ν(k) p n(k) p n+1 (k) = (k n λ2 2 ) p n (k) λ2 2 n p n 1 (k) Charlier polynomials

8 obviously amount to linear automorphisms of the Weyl algebra Algebraic context allows to derive various of the Charlier polynomials In its difference and to show that they span a representa@on of the Weyl algebra

9 4. d ORTHOGONAL POLYNOMIALS P n (x) monic, degree n defined by recurrence rela@on of order d+1 Ini@al condi@ons: P n+1 (x) = xp n (x) d s=0 a n, n s P n s (x) P k (x) = 0, k < 0 ; P o (x) =1 d=1 : 3 term P n+1 (x) = xp n (x) a n, n P n (x) a n, n 1 P n 1 (x)

10 Provided u n a n, n 1 0 for d=1, orthogonal P n σ,p n (x)p m (x) = h n δ n, m h n = u 1 u 2 u n func@onal σ defined by moments c n = σ,x n

11 When d > 1 : vector orthogonality, d linear σ k k = 0,1,..., d 1 func@onals such that σ k,p n P m = 0 m > dn + k σ k,p n P dn+ k 0 n 0 E.g. d = 2 σ 0,P n P m = 0 m > 2n σ 1,P n P m = 0 m > 2n +1 σ 0,P n P 2n 0 σ 1,P n P 2n+1 0 n 0

12 Padé problems Lanczos type method for solving linear systems Random matrices

13 5. A NONLINEAR AUTOMORPHISM OF THE WEYL ALGEBRA S = e βa* e Q (a) S 1 = e Q(a) e βa* β: parameter, Q polynomial of degree N 1 SaS 1 = e βa* a e βa* = a β S 1 as = a + β Sa * S 1 = e βa* e Q(a) a * e Q(a) e βa* = e βa* ( a * + [ Q(a),a *]+ 1 [ * Q(a), [ Q(a),a ] 2 ]+ )e βa*

14 = e βa* ( a * + Q (a))e βa* = a * + Q (a β) S 1 a * S = a * Q (a) Note for deg (Q) =1, Q = const S generates transla@ons (linear automorphisms)

15 6. d-charlier POLYNOMIALS Calculate ψ nk = k S n through recurrence relation k a * a S n = k k S n = kψ nk k a * as n = k SS 1 a * as n = k S( a * Q (a))(a + β) n Q (x)(x + β) polynomial of degree N if deg Q=N Write Q (x)(x + β) = N ξ i x i i= 0 ξ i coefficients

16 Comparing ( a n = n n 1, a * n = n +1 n +1 ) k k S n = k S a * a + βa * Gives n=0: N i= 0 ξ i a i n = (n ξ 0 ) k S n + β n +1 k S n +1 ξ i n(n 1) (n i +1) k S n i kψ n,k = (n ξ 0 )ψ n,k + β n +1ψ n +1,k ξ i n(n 1) (n i +1)ψ n i,k i=1 recursively from N N i=1 ψ n,k ψ 0,k βψ 1,k = (k + ξ 0 )ψ o,k

17 N kψ n, k = (n ξ 0 )ψ n, k + β n +1ψ n+1, k ξ i i =1 n(n 1) (n i +1)ψ n i, k n=1 β 2ψ 2, k = (k = 1+ ξ 0 )ψ 1, k + ξ 1 ψ 0, k = [(k 1+ ξ 0 )(k + ξ 0 )+ ξ 1 ] ψ 0, k Hence P n (k) ψ n, k = P n (k) ψ 0, k polynomial of degree n in the discrete variable k sa@sfying the recurrence rela@ons

18 β n +1 P n+1 (k) = (k n + ξ 0 )P n (k)+ ξ i N i =1 n(n 1) (n i +1) P n i (k) P 0 (k) = 1, P n (k) = 0 n < 0 Recurrence rela@ons of order N+1 ψ 0,k = k S 0 = k e βa* e Q(a) 0 e Q(a) 0 = e Q(0) 0 k e βa* 0 = k β k (a * ) k k! 0 = β k k!

19 ψ 0, k = e Q (0) β k k! For N=1 recover ordinary case of 3 term recurrence rela@on β n +1 P n+1 (k) = (k n + ξ 0 )+ ξ 1 P n 1 (k) To make polynomials monic, set ˆ P n (k) = β n n! P n (k)

20 β n +1 β n 1 (n +1)! ˆ P n+1 = (k n + ξ 0 ) β n n! ˆ P n + N ξ i n(n 1) (n i +1) i=1 n +i β (n i)! ˆ P n i P ˆ n+1 (K ) = (k n + ξ 0 ) P ˆ n (k)+ N i =1 β i ξ i d orthogonal Charlier polynomials n(n 1) (n i +1) ˆ P n i (k) For N=1, ordinary orthogonal Charlier polynomials P ˆ n+1 (k) = (k n + ξ 0 ) P ˆ n (k)+ βξ 1 n P ˆ n 1 (k) up to shift k k + βξ 1 ξ 0

21 7. PROPERTIES Difference Start again from ψ n, k = k S n = k e βa* e Q (a) n = β n now need a recurrence relation in variable k nψ n, k = k S a * a n k S a * a n = k S a * a S 1 S n n! ˆ P n (k)ψ 0, k = k ( a * + Q (a β) )(a β)s n using the dual rela@ons

22 (x β) Q (x β) = N x i i =0 η i nψ n, k = k a * a β a * + N i =0 η i a i S n = kψ n, k β k ψ n, k 1 + η i N i =0 (k +1) (k + i)ψ n, k+ i using ψ n, k = β n n! ˆ P n (k)ψ 0, k ψ 0, k = e Q (0) β k k! we find (n k) ˆ P n (k) + k ˆ P n (k 1) = N η i β i i= 0 ˆ P n (k + i)

23 Use the standard difference operators Δf (k) = f (k +1) f (k) f(k) = f (k) f (k 1) and note that F(x + m) = m i =0 m Δ i F(x) i (n k) P ˆ n (k)+ k P ˆ n (k 1) = N β i i =0 η i ˆ P n (k + i) becomes n P ˆ n (k) = k P ˆ n (k)+ N i i η i β i Δ l P ˆ n (k) l i =0 l =0

24 n P ˆ n (k) = k P ˆ n (k)+ N i i η i β i Δ l P ˆ l n (k) i =0 l =0 Recall (x β) Q (x β) = N x i i =0 η i y Q (y) = N η i (y + β) i = i= 0 N i η i i= 0 l= 0 i N y l β i l = η l i β i i= 0 i l= 0 i l y β l n P ˆ n (k) = k P ˆ n (k)+ N µ i Δ i P ˆ n (k) i =0 N j =0 µ j y β j Eigenvalue equation: H ˆ P n = n ˆ P n

25 Ladder Operators Relation n ψ n 1, k = k S a n = k (a β)s n = k +1ψ n, k+1 βψ n, k Δ P ˆ n (k) P ˆ n (k +1) P ˆ n (k) = n P ˆ n 1 (k) Difference Appell property Similarly from n +1 ψ n+1, k = k S a * n = k ( a * + Q (a β) ) S n

26 we get P ˆ n+1 (k) = k P ˆ n (k 1)+ Q (x β) = With N 1 ν i x i i =0 N 1 ν i β i+1 P ˆ n (k + i) i =0 let N 1 T ± F(k) = F(k ±1) R = k T + ν i β i+1 i T + i= 0 Rodrigues formula : Factorization : Realization : H = RΔ [ Δ,R] = 1 ˆ P n (k) = R n 1 { }

27 Generating function Consider the generating function F(z;k) = e Q (0) β k k! n=0 ψ n, k z n n! recalling that ψ n, k = e Q (0) β k n n!k! ˆ P n (k) it yields F(z,k) = n=0 z β n! n ˆ P n (k)

28 Coherent states : z = e za* 0 n z = zn n! a z = z z hence F(a) z = F(z) z Coming back n=0 F(z,k) = e Q (0) β k k! = e Q (0) β k k! k S n n z = e Q (0) β k k! k S z n=0 ψ n, k z n n! k S z = k e βa* e Q(a) z

29 is calculated using e Q(a) z = e Q(z) z e βa* z = e βa*. e za* 0 = e (β + z) a* 0 = z + β to find We thus have (z + β) k k S z = e Q(z) k z + β = e Q(z) k! F(z;k) = e Q(z) Q (0) 1+ z β k n=0 = z β n ˆ P n (k) n! as genera@ng func@on for the d Charlier polynomials.

30 Characterization Theorem: N+1 order recurrence relation Matrix elements of S (non linear automorphisms) Appell difference property eigenvectors of H All equivalent characterization of d-orthogonal Charlier polynomials

31 8. EXPLICIT EXPRESSION IN A SPECIAL CASE Take Expanding F and some transformations yield a m = (q k + m) N m = 0,,N 1 b m = (q + m +1) N m = 0,,N 1 excluding (q m+1) N ( k) q = ( k)( k +1) ( k + q 1)

32 For N=1, we retrieve the standard expression for the Charlier polynomials:

33 9. DUAL FUNCTIONS AND ORTHOGONALITY Much more can be obtained from framework, in particular orthogonality properties Assume S invertible and define ψ n, k = k S n φ n,k = n S 1 k k =0 m S 1 k k S n = δ m, n implies the biorthogonality property φ m, k k =0 ψ n, k = δ m, n

34 φ m, k - also expressible in terms of d-charlier polynomials - Generating function for can be obtained Noteworthy: Can relate orthogonality of and to vector orthogonality of d-orthogonal polynomials φ φ m, k ψ In fact, from the biorthogonality relations and the knowledge of φ m, k can construct explicit realization of linear functionals expressing vector orthogonality.

35 10. CONCLUSIONS Thorough of d orthogonal Charlier polynomials from matrix elements of nonlinear automorphism of Weyl algebra Can establish rela@on between Meixner polynomials and linear automorphism of SU(1,1) much like the one between Charlier and Weyl. (Floreanini, LeTourneux, Vinet) One uses dila@on instead of transla@ons Can expect other d orthogonal polynomials to appear as matrix elements of non unitary automorphisms of algebra. We have worked out the case of Krawtchouk d orthogonal polynomials

arxiv: v3 [math-ph] 9 Feb 2012

arxiv: v3 [math-ph] 9 Feb 2012 d-orthogonal polynomials and su(2) Vincent X. Genest 1, Luc Vinet 1, and Alexei Zhedanov 2 arxiv:1110.1902v3 [math-ph] 9 Feb 2012 1 Centre de recherches mathématiques, Université de Montréal, C.P. 6128,

More information

PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS

PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS ELLIPTIC INTEGRABLE SYSTEMS PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS A.S. ZHEDANOV Abstract. We study recurrence relations and biorthogonality properties for polynomials and rational

More information

Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013

Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv: v2 [math.ca] 8 Nov 2013 Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv:1310.0982v2 [math.ca] 8 Nov 2013 F Ndayiragie 1 and W Van Assche 2 1 Département de Mathématiques, Université du Burundi, Campus

More information

Spectral Transformations of the Laurent Biorthogonal Polynomials

Spectral Transformations of the Laurent Biorthogonal Polynomials Spectral Transformations of the Laurent Biorthogonal Polynomials Luc Vinet Alexei Zhedanov CRM-2616 June 1999 Centre de recherches mathématiques, Université de Montréal, C.P. 6128, succ. Centre-Ville,

More information

Spectral Transformations of the Laurent Biorthogonal Polynomials. II. Pastro Polynomials

Spectral Transformations of the Laurent Biorthogonal Polynomials. II. Pastro Polynomials Spectral Transformations of the Laurent Biorthogonal Polynomials. II. Pastro Polynomials Luc Vinet Alexei Zhedanov CRM-2617 June 1999 Centre de recherches mathématiques, Université de Montréal, C.P. 6128,

More information

Markov operators, classical orthogonal polynomial ensembles, and random matrices

Markov operators, classical orthogonal polynomial ensembles, and random matrices Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random

More information

MATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,

MATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3, MATH 205 HOMEWORK #3 OFFICIAL SOLUTION Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. a F = R, V = R 3, b F = R or C, V = F 2, T = T = 9 4 4 8 3 4 16 8 7 0 1

More information

Problem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that

Problem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume

More information

Linearization coefficients for orthogonal polynomials. Michael Anshelevich

Linearization coefficients for orthogonal polynomials. Michael Anshelevich Linearization coefficients for orthogonal polynomials Michael Anshelevich February 26, 2003 P n = monic polynomials of degree n = 0, 1,.... {P n } = basis for the polynomials in 1 variable. Linearization

More information

arxiv: v2 [math-ph] 7 Nov 2013

arxiv: v2 [math-ph] 7 Nov 2013 Interbasis expansions for the isotropic 3D harmonic oscillator and bivariate Krawtchouk polynomials arxiv:1307.0692v2 [math-ph] 7 Nov 2013 Vincent X. Genest E-mail: genestvi@crm.umontreal.ca Centre de

More information

Alphonse Magnus, Université Catholique de Louvain.

Alphonse Magnus, Université Catholique de Louvain. Rational interpolation to solutions of Riccati difference equations on elliptic lattices. Alphonse Magnus, Université Catholique de Louvain. http://www.math.ucl.ac.be/membres/magnus/ Elliptic Riccati,

More information

Elliptic Biorthogonal Polynomials Connected with Hermite s Continued Fraction

Elliptic Biorthogonal Polynomials Connected with Hermite s Continued Fraction Symmetry, Integrability and Geometry: Methods and Applications Elliptic Biorthogonal Polynomials Connected with Hermite s Continued Fraction Luc VINET and Alexei ZHEDANOV SIGMA 3 (27), 3, 18 pages Université

More information

Jordan Normal Form. Chapter Minimal Polynomials

Jordan Normal Form. Chapter Minimal Polynomials Chapter 8 Jordan Normal Form 81 Minimal Polynomials Recall p A (x) =det(xi A) is called the characteristic polynomial of the matrix A Theorem 811 Let A M n Then there exists a unique monic polynomial q

More information

Generalized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.

Generalized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems. Generalized Burgers equations and Miura Map in nonabelian rings as integrable systems. Sergey Leble Gdansk University of Technology 05.07.2015 Table of contents 1 Introduction: general remarks 2 Remainders

More information

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

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

Lecture 4.6: Some special orthogonal functions

Lecture 4.6: Some special orthogonal functions Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics

More information

MIT Final Exam Solutions, Spring 2017

MIT Final Exam Solutions, Spring 2017 MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of

More information

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v ) Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O

More information

Introduction to orthogonal polynomials. Michael Anshelevich

Introduction to orthogonal polynomials. Michael Anshelevich Introduction to orthogonal polynomials Michael Anshelevich November 6, 2003 µ = probability measure on R with finite moments m n (µ) = R xn dµ(x)

More information

Math 123 Homework Assignment #2 Due Monday, April 21, 2008

Math 123 Homework Assignment #2 Due Monday, April 21, 2008 Math 123 Homework Assignment #2 Due Monday, April 21, 2008 Part I: 1. Suppose that A is a C -algebra. (a) Suppose that e A satisfies xe = x for all x A. Show that e = e and that e = 1. Conclude that e

More information

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. A linear functional µ[p ] = P (x) dµ(x), µ

More information

Legendre s Equation. PHYS Southern Illinois University. October 18, 2016

Legendre s Equation. PHYS Southern Illinois University. October 18, 2016 Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying

More information

A classification of sharp tridiagonal pairs. Tatsuro Ito, Kazumasa Nomura, Paul Terwilliger

A classification of sharp tridiagonal pairs. Tatsuro Ito, Kazumasa Nomura, Paul Terwilliger Tatsuro Ito Kazumasa Nomura Paul Terwilliger Overview This talk concerns a linear algebraic object called a tridiagonal pair. We will describe its features such as the eigenvalues, dual eigenvalues, shape,

More information

Power Series Solutions to the Legendre Equation

Power Series Solutions to the Legendre Equation Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial

More information

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich

Measures, orthogonal polynomials, and continued fractions. Michael Anshelevich Measures, orthogonal polynomials, and continued fractions Michael Anshelevich November 7, 2008 MEASURES AND ORTHOGONAL POLYNOMIALS. MEASURES AND ORTHOGONAL POLYNOMIALS. µ a positive measure on R. 2 MEASURES

More information

Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35

Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35 Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35 1. Let R be a commutative ring with 1 0. (a) Prove that the nilradical of R is equal to the intersection of the prime

More information

Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015

Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic

More information

Decomposition of a recursive family of polynomials

Decomposition of a recursive family of polynomials Decomposition of a recursive family of polynomials Andrej Dujella and Ivica Gusić Abstract We describe decomposition of polynomials f n := f n,b,a defined by f 0 := B, f 1 (x := x, f n+1 (x = xf n (x af

More information

Math 554 Qualifying Exam. You may use any theorems from the textbook. Any other claims must be proved in details.

Math 554 Qualifying Exam. You may use any theorems from the textbook. Any other claims must be proved in details. Math 554 Qualifying Exam January, 2019 You may use any theorems from the textbook. Any other claims must be proved in details. 1. Let F be a field and m and n be positive integers. Prove the following.

More information

Simple Harmonic Oscillator

Simple Harmonic Oscillator Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in

More information

Linear Algebra Review. Vectors

Linear Algebra Review. Vectors Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors

More information

The Simple Harmonic Oscillator

The Simple Harmonic Oscillator The Simple Harmonic Oscillator Asaf Pe er 1 November 4, 215 This part of the course is based on Refs [1] [3] 1 Introduction We return now to the study of a 1-d stationary problem: that of the simple harmonic

More information

January 2016 Qualifying Examination

January 2016 Qualifying Examination January 2016 Qualifying Examination If you have any difficulty with the wording of the following problems please contact the supervisor immediately. All persons responsible for these problems, in principle,

More information

The Lorentz and Poincaré groups. By Joel Oredsson

The Lorentz and Poincaré groups. By Joel Oredsson The Lorentz and Poincaré groups By Joel Oredsson The Principle of Special Rela=vity: The laws of nature should be covariant with respect to the transforma=ons between iner=al reference frames. x µ x' µ

More information

Numerical Linear Algebra Homework Assignment - Week 2

Numerical Linear Algebra Homework Assignment - Week 2 Numerical Linear Algebra Homework Assignment - Week 2 Đoàn Trần Nguyên Tùng Student ID: 1411352 8th October 2016 Exercise 2.1: Show that if a matrix A is both triangular and unitary, then it is diagonal.

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

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors 5 Eigenvalues and Eigenvectors 5.2 THE CHARACTERISTIC EQUATION DETERMINANATS nn Let A be an matrix, let U be any echelon form obtained from A by row replacements and row interchanges (without scaling),

More information

Frame Diagonalization of Matrices

Frame Diagonalization of Matrices Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)

More information

Unit 5. Matrix diagonaliza1on

Unit 5. Matrix diagonaliza1on Unit 5. Matrix diagonaliza1on Linear Algebra and Op1miza1on Msc Bioinforma1cs for Health Sciences Eduardo Eyras Pompeu Fabra University 218-219 hlp://comprna.upf.edu/courses/master_mat/ We have seen before

More information

1. Foundations of Numerics from Advanced Mathematics. Linear Algebra

1. Foundations of Numerics from Advanced Mathematics. Linear Algebra Foundations of Numerics from Advanced Mathematics Linear Algebra Linear Algebra, October 23, 22 Linear Algebra Mathematical Structures a mathematical structure consists of one or several sets and one or

More information

The oscillator model for the Lie superalgebra sh(2 2) and Charlier polynomials

The oscillator model for the Lie superalgebra sh(2 2) and Charlier polynomials The oscillator model for the Lie superalgebra sh( ) and Charlier polynomials E.I. Jafarov 1 and J. Van der Jeugt Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 81-S9,

More information

Jordan normal form notes (version date: 11/21/07)

Jordan normal form notes (version date: 11/21/07) Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let

More information

Tutorial 5 Clifford Algebra and so(n)

Tutorial 5 Clifford Algebra and so(n) Tutorial 5 Clifford Algebra and so(n) 1 Definition of Clifford Algebra A set of N Hermitian matrices γ 1, γ,..., γ N obeying the anti-commutator γ i, γ j } = δ ij I (1) is the basis for an algebra called

More information

Chapter 4 Euclid Space

Chapter 4 Euclid Space Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,

More information

X n = c n + c n,k Y k, (4.2)

X n = c n + c n,k Y k, (4.2) 4. Linear filtering. Wiener filter Assume (X n, Y n is a pair of random sequences in which the first component X n is referred a signal and the second an observation. The paths of the signal are unobservable

More information

ECE Unit 4. Realizable system used to approximate the ideal system is shown below: Figure 4.47 (b) Digital Processing of Analog Signals

ECE Unit 4. Realizable system used to approximate the ideal system is shown below: Figure 4.47 (b) Digital Processing of Analog Signals ECE 8440 - Unit 4 Digital Processing of Analog Signals- - Non- Ideal Case (See sec8on 4.8) Before considering the non- ideal case, recall the ideal case: 1 Assump8ons involved in ideal case: - no aliasing

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

Eigenvalues and Eigenvectors A =

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

arxiv: v1 [math-ph] 26 May 2017

arxiv: v1 [math-ph] 26 May 2017 arxiv:1705.09737v1 [math-ph] 26 May 2017 An embedding of the Bannai Ito algebra in U (osp(1, 2)) and 1 polynomials Pascal Baseilhac Laboratoire de Mathématiques et Physique Théorique CNRS/UMR 7350, Fédération

More information

Linear Algebra. Workbook

Linear Algebra. Workbook Linear Algebra Workbook Paul Yiu Department of Mathematics Florida Atlantic University Last Update: November 21 Student: Fall 2011 Checklist Name: A B C D E F F G H I J 1 2 3 4 5 6 7 8 9 10 xxx xxx xxx

More information

Matrix products in integrable probability

Matrix products in integrable probability Matrix products in integrable probability Atsuo Kuniba (Univ. Tokyo) Mathema?cal Society of Japan Spring Mee?ng Tokyo Metropolitan University 27 March 2017 Non-equilibrium sta?s?cal mechanics Stochas?c

More information

Free Meixner distributions and random matrices

Free Meixner distributions and random matrices Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2

More information

GG612 Lecture 3. Outline

GG612 Lecture 3. Outline GG612 Lecture 3 Strain 11/3/15 GG611 1 Outline Mathema8cal Opera8ons Strain General concepts Homogeneous strain E (strain matri) ε (infinitesimal strain) Principal values and principal direc8ons 11/3/15

More information

From Bernstein approximation to Zauner s conjecture

From Bernstein approximation to Zauner s conjecture From Bernstein approximation to Zauner s conjecture Shayne Waldron Mathematics Department, University of Auckland December 5, 2017 Shayne Waldron (University of Auckland) Workshop on Spline Approximation

More information

arxiv:math/ v1 [math.ca] 9 Jul 1993

arxiv:math/ v1 [math.ca] 9 Jul 1993 The q-harmonic Oscillator and the Al-Salam and Carlitz polynomials Dedicated to the Memory of Professor Ya. A. Smorodinskiĭ R. Askey and S. K. Suslov arxiv:math/9307207v1 [math.ca] 9 Jul 1993 Abstract.

More information

NONCOMMUTATIVE POLYNOMIAL EQUATIONS. Edward S. Letzter. Introduction

NONCOMMUTATIVE POLYNOMIAL EQUATIONS. Edward S. Letzter. Introduction NONCOMMUTATIVE POLYNOMIAL EQUATIONS Edward S Letzter Introduction My aim in these notes is twofold: First, to briefly review some linear algebra Second, to provide you with some new tools and techniques

More information

CHAPTER VIII HILBERT SPACES

CHAPTER VIII HILBERT SPACES CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)

More information

Quantum Computing Lecture 2. Review of Linear Algebra

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

Chapter 12. Algebraic numbers and algebraic integers Algebraic numbers

Chapter 12. Algebraic numbers and algebraic integers Algebraic numbers Chapter 12 Algebraic numbers and algebraic integers 12.1 Algebraic numbers Definition 12.1. A number α C is said to be algebraic if it satisfies a polynomial equation with rational coefficients a i Q.

More information

Field Theory Qual Review

Field Theory Qual Review Field Theory Qual Review Robert Won Prof. Rogalski 1 (Some) qual problems ˆ (Fall 2007, 5) Let F be a field of characteristic p and f F [x] a polynomial f(x) = i f ix i. Give necessary and sufficient conditions

More information

Economics 472. Lecture 10. where we will refer to y t as a m-vector of endogenous variables, x t as a q-vector of exogenous variables,

Economics 472. Lecture 10. where we will refer to y t as a m-vector of endogenous variables, x t as a q-vector of exogenous variables, University of Illinois Fall 998 Department of Economics Roger Koenker Economics 472 Lecture Introduction to Dynamic Simultaneous Equation Models In this lecture we will introduce some simple dynamic simultaneous

More information

Model reduction of large-scale dynamical systems

Model reduction of large-scale dynamical systems Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/

More information

1 Unitary representations of the Virasoro algebra

1 Unitary representations of the Virasoro algebra Week 5 Reading material from the books Polchinski, Chapter 2, 15 Becker, Becker, Schwartz, Chapter 3 Ginspargs lectures, Chapters 3, 4 1 Unitary representations of the Virasoro algebra Now that we have

More information

Linear Algebra Practice Problems

Linear Algebra Practice Problems Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,

More information

arxiv: v2 [math-ph] 10 Dec 2015

arxiv: v2 [math-ph] 10 Dec 2015 A superintegrable discrete harmonic oscillator based on bivariate Charlier polynomials Vincent X. Genest, 1, Hiroshi Miki,, Luc Vinet, 3,4, and Guofu Yu 4, 1 Department of Mathematics, Massachusetts Institute

More information

Group Theory and Vibrational Spectroscopy

Group Theory and Vibrational Spectroscopy Group Theory and Vibrational Spectroscopy Pamela Schleissner Physics 251 Spring 2017 Outline Molecular Symmetry Representations of Molecular Point Groups Group Theory and Quantum Mechanics Vibrational

More information

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs Ma/CS 6b Class 3: Eigenvalues in Regular Graphs By Adam Sheffer Recall: The Spectrum of a Graph Consider a graph G = V, E and let A be the adjacency matrix of G. The eigenvalues of G are the eigenvalues

More information

ACM 104. Homework Set 5 Solutions. February 21, 2001

ACM 104. Homework Set 5 Solutions. February 21, 2001 ACM 04 Homework Set 5 Solutions February, 00 Franklin Chapter 4, Problem 4, page 0 Let A be an n n non-hermitian matrix Suppose that A has distinct eigenvalues λ,, λ n Show that A has the eigenvalues λ,,

More information

SQUARE ROOTS OF 2x2 MATRICES 1. Sam Northshield SUNY-Plattsburgh

SQUARE ROOTS OF 2x2 MATRICES 1. Sam Northshield SUNY-Plattsburgh SQUARE ROOTS OF x MATRICES Sam Northshield SUNY-Plattsburgh INTRODUCTION A B What is the square root of a matrix such as? It is not, in general, A B C D C D This is easy to see since the upper left entry

More information

Diagonal Representation of Density Matrix Using q-coherent States

Diagonal Representation of Density Matrix Using q-coherent States Proceedings of Institute of Mathematics of NAS of Ukraine 24, Vol. 5, Part 2, 99 94 Diagonal Representation of Density Matrix Using -Coherent States R. PARTHASARATHY and R. SRIDHAR The Institute of Mathematical

More information

1 2 2 Circulant Matrices

1 2 2 Circulant Matrices Circulant Matrices General matrix a c d Ax x ax + cx x x + dx General circulant matrix a x ax + x a x x + ax. Evaluating the Eigenvalues Find eigenvalues and eigenvectors of general circulant matrix: a

More information

THE JOHNS HOPKINS UNIVERSITY Faculty of Arts and Sciences FINAL EXAM - SPRING SESSION ADVANCED ALGEBRA II.

THE JOHNS HOPKINS UNIVERSITY Faculty of Arts and Sciences FINAL EXAM - SPRING SESSION ADVANCED ALGEBRA II. THE JOHNS HOPKINS UNIVERSITY Faculty of Arts and Sciences FINAL EXAM - SPRING SESSION 2006 110.402 - ADVANCED ALGEBRA II. Examiner: Professor C. Consani Duration: 3 HOURS (9am-12:00pm), May 15, 2006. No

More information

Introduc)on to linear algebra

Introduc)on to linear algebra Introduc)on to linear algebra Vector A vector, v, of dimension n is an n 1 rectangular array of elements v 1 v v = 2 " v n % vectors will be column vectors. They may also be row vectors, when transposed

More information

On rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro

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

Math Ordinary Differential Equations

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

Ket space as a vector space over the complex numbers

Ket space as a vector space over the complex numbers Ket space as a vector space over the complex numbers kets ϕ> and complex numbers α with two operations Addition of two kets ϕ 1 >+ ϕ 2 > is also a ket ϕ 3 > Multiplication with complex numbers α ϕ 1 >

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

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

The Sommerfeld Polynomial Method: Harmonic Oscillator Example

The Sommerfeld Polynomial Method: Harmonic Oscillator Example Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic

More information

GG612 Lecture 3. Outline

GG612 Lecture 3. Outline GG61 Lecture 3 Strain and Stress Should complete infinitesimal strain by adding rota>on. Outline Matrix Opera+ons Strain 1 General concepts Homogeneous strain 3 Matrix representa>ons 4 Squares of line

More information

Some Explicit Biorthogonal Polynomials

Some Explicit Biorthogonal Polynomials Some Explicit Biorthogonal Polynomials D. S. Lubinsky and H. Stahl Abstract. Let > and ψ x) x. Let S n, be a polynomial of degree n determined by the biorthogonality conditions Z S n,ψ,,,..., n. We explicitly

More information

The character table of a split extension of the Heisenberg group H 1 (q) by Sp(2, q), q odd

The character table of a split extension of the Heisenberg group H 1 (q) by Sp(2, q), q odd Università di Milano Bicocca Quaderni di Matematica The character table of a split extension of the Heisenberg group H (q by Sp(, q, q odd Marco Antonio Pellegrini Quaderno n. 7/8 (arxiv:math/85.48 Stampato

More information

MATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial.

MATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial. MATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial. Geometric properties of determinants 2 2 determinants and plane geometry

More information

Calculus C (ordinary differential equations)

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

Notes on Linear Algebra and Matrix Theory

Notes on Linear Algebra and Matrix Theory Massimo Franceschet featuring Enrico Bozzo Scalar product The scalar product (a.k.a. dot product or inner product) of two real vectors x = (x 1,..., x n ) and y = (y 1,..., y n ) is not a vector but a

More information

Section 33 Finite fields

Section 33 Finite fields Section 33 Finite fields Instructor: Yifan Yang Spring 2007 Review Corollary (23.6) Let G be a finite subgroup of the multiplicative group of nonzero elements in a field F, then G is cyclic. Theorem (27.19)

More information

Problem 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. 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 information

Lecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University

Lecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1

More information

Spectra of Digraph Transformations

Spectra of Digraph Transformations Spectra of Digraph Transformations Aiping Deng a,, Alexander Kelmans b,c a Department of Applied Mathematics, Donghua University, 201620 Shanghai, China arxiv:1707.00401v1 [math.co] 3 Jul 2017 b Department

More information

INTRODUCTION TO LIE ALGEBRAS. LECTURE 2.

INTRODUCTION TO LIE ALGEBRAS. LECTURE 2. INTRODUCTION TO LIE ALGEBRAS. LECTURE 2. 2. More examples. Ideals. Direct products. 2.1. More examples. 2.1.1. Let k = R, L = R 3. Define [x, y] = x y the cross-product. Recall that the latter is defined

More information

Numerical Methods in Physics

Numerical Methods in Physics Numerical Methods in Physics Numerische Methoden in der Physik, 515.421. Instructor: Ass. Prof. Dr. Lilia Boeri Room: PH 03 090 Tel: +43-316- 873 8191 Email Address: l.boeri@tugraz.at Room: TDK Seminarraum

More information

Triangular matrices and biorthogonal ensembles

Triangular matrices and biorthogonal ensembles /26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n

More information

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions

Recursion Systems and Recursion Operators for the Soliton Equations Related to Rational Linear Problem with Reductions GMV The s Systems and for the Soliton Equations Related to Rational Linear Problem with Reductions Department of Mathematics & Applied Mathematics University of Cape Town XIV th International Conference

More information

CS 246 Review of Linear Algebra 01/17/19

CS 246 Review of Linear Algebra 01/17/19 1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector

More information

Solutions Preliminary Examination in Numerical Analysis January, 2017

Solutions Preliminary Examination in Numerical Analysis January, 2017 Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which

More information

Linear algebra II Homework #1 solutions A = This means that every eigenvector with eigenvalue λ = 1 must have the form

Linear algebra II Homework #1 solutions A = This means that every eigenvector with eigenvalue λ = 1 must have the form Linear algebra II Homework # solutions. Find the eigenvalues and the eigenvectors of the matrix 4 6 A =. 5 Since tra = 9 and deta = = 8, the characteristic polynomial is f(λ) = λ (tra)λ+deta = λ 9λ+8 =

More information

Difference Equations for Multiple Charlier and Meixner Polynomials 1

Difference Equations for Multiple Charlier and Meixner Polynomials 1 Difference Equations for Multiple Charlier and Meixner Polynomials 1 WALTER VAN ASSCHE Department of Mathematics Katholieke Universiteit Leuven B-3001 Leuven, Belgium E-mail: walter@wis.kuleuven.ac.be

More information

A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit

A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship

More information

Counting Perfect Polynomials

Counting Perfect Polynomials Enrique Treviño joint work with U. Caner Cengiz and Paul Pollack 49th West Coast Number Theory December 18, 2017 49th West Coast Number Theory 2017 1 Caner (a) Caner Cengiz (b) Paul Pollack 49th West Coast

More information

Representations of quantum SU(2) operators on a local chart

Representations of quantum SU(2) operators on a local chart Journal of Physics: Conference Series PAPER OPEN ACCESS Representations of quantum SU() operators on a local chart To cite this article: Elmar Wagner 016 J. Phys.: Conf. Ser. 670 01053 View the article

More information