arxiv: v3 [math.fa] 23 Sep 2016


 Arleen Robinson
 2 years ago
 Views:
Transcription
1 AN INVERSE PROBLEM FOR A CLASS OF CANONICAL SYSTEMS AND ITS APPLICATIONS TO SELFRECIPROCAL POLYNOMIALS MASATOSHI SUZUKI arxiv: v3 [mathfa] 23 Sep 2016 Abstract A canonical system is a kind of firstorder system of ordinary differential equations on an interval of the real line parametrized by complex numbers It is known that any solution of a canonical system generates an entire function of the HermiteBiehler class In this paper, we deal with the inverse problem to recover a canonical system from a given entire function of the HermiteBiehler class satisfying appropriate conditions This inverse problem was solved by de Branges in 1960s However his results are often not enough to investigate a Hamiltonian of recovered canonical system In this paper, we present an explicit way to recover a Hamiltonian from a given exponential polynomial belonging to the HermiteBiehler class After that, we apply it to study distributions of roots of selfreciprocal polynomials 1 Introduction Let H(a) be a 2 2 symmetric matrixvalued function defined almost everywhere on a finite interval I = [a 1,a 0 ) (0 < a 1 < a 0 < ) with respect to the Lebesgue measure da We refer to a firstorder system of differential equations a d [ A(a,z) 0 1 A(a,z) A(a,z) 1 = z H(a), lim = (11) da B(a, z) 1 0 B(a, z) aրa 0 B(a, z) 0] on I parametrized by z C as a quasicanonical system (on I) A column vectorvalued function (A(,z),B(,z)) : I C 2 1 is called a solution if it consists of absolutely continuous functions and satisfies (11) almost everywhere on I for every fixed z C A quasicanonical system (11) is called a canonical system if (H1) H(a) is a real positive semidefinite symmetric matrix for almost every a I, (H2) H(a) 0 on any subset of I with positive Lebesgue measure, (H3) H(a) is locally integrable on I with respect to da/a, and H(a) is called its Hamiltonian Even if a quasicanonical system (11) is not a canonical system, a matrix valued function H(a) is often called a Hamiltonian in this paper by abuse of language if there is no confusion A number of different secondorder differential equations, like Schrödinger equations and Sturm Liouville equations of appropriate form, and systems of firstorder differential equations like Dirac type systems of appropriate form are reduced to a canonical system Fundamental results on the spectral theory of canonical systems was established in works of Gohberg Kreĭn [5], de Branges [3] and many other authors (see the survey articles Winkler [27], Woracek [28] and references there in for historical details on canonical systems) Note that the variable a in (11) is the variable on the multiplicative group R >0 By the change of variable a = e x, (11) is transformed into the equation for the variable x on the additive group R, and the right endpoint of a for the initial value is transformed into the left endpoint of x for the initial value The transformed equation is the one treated by de Branges, Kreĭn, Kac and others (See the final paragraph of the introduction for the reason why we use a multiplicative variable) The subject of the present paper is an inverse spectral problem for canonical systems In order to state it, we review the theory of the Hermite Biehler class 2000 Mathematics Subject Classification 30C15, 34A55, 34L40 1
2 2 M SUZUKI We use the notation F (z) = F( z) for functions of the complex variable z and denote by C + the (open) upper halfplane {z = x +iy C y > 0} An entire function E(z) satisfying E (z) < E(z) for every z C + (12) and having no real zeros is said to be a function of the Hermite Biehler class, or the class HB for short (This definition of the class HB is equivalent to the definition of Levin [11, 1 of Chap VII] if we replace the word the upper halfplane by the lower halfplane, because (12) implies that E(z) has no zeros in C + We adopt the above definition for the convenience to use the theory of canonical systems via the theory of de Branges spaces) Suppose that (11) is a canonical system endowed with a solution (A(a,z),B(a,z)) Then E(a,z) := A(a,z) ib(a,z) is a function of the class HB for every fixed regular point a I Inparticular, lim aրa0 E(a,z) = 1 and E(a 1,z) is a function of theclass HB Therefore, an inverse problem for canonical systems is to recover their Hamiltonians from given entire functions of the class HB satisfying appropriate conditions Usually such inverse problem is difficult to solve in general, because it is an inverse spectral problem ([8, Section 2], [17, Section 7]) However it is already solved by the theory of de Branges in 1960s ([3]) For example, for fixed I and an entire function E(z) of exponential type belonging to the class HB and satisfying (1+x2 ) 1 E(x) 2 dx < and E(0) = 1, there exists a canonical system (that is, there exists a Hamiltonian H(a) on I) such that A(a 1,z) = 1 2 (E(z)+E (z)), B(a 1,z) = i 2 (E(z) E (z)) for its solution (A(a,z),B(a,z)) (see [17, Theorem 73] with [4, Lemma 33]) Moreover, such canonical system is uniquely determined by E(z) and I under appropriate normalizations More general situation is treated in Kac [7] (see also [27] and [28]) As above, de Bragnes s theory ensures the existence of canonical systems or Hamiltonians for given functions of the class HB, but it does not provide explicit or useful expressions of Hamiltonians In fact, an explicit form of H(a) is not known except for a few examples of E(z) as in [3, Chapter 3], [4, Section 8] and some additional examples constructed from such known examples using transformation rules for Hamiltonians and Weyl functions ([26]) In this paper, we deal with the above inverse problem for a special class of exponential polynomials together with the problem of explicit constructions of H(a), and apply the results to the studying of the distribution of roots of selfreciprocal polynomials Let R := R\{0}, g Z >0 and q > 1 We denote by C a vector of length 2g +1 of the form C = (C g,c g 1,,C g ) R R 2g 1 R and consider exponential polynomials along with two associated functions E(z) := E q (z;c) := C m q imz (13) A(z) := A q (z;c) := 1 2 (E q(z;c)+e q(z;c)), (14) B(z) := B q (z;c) := i 2 (E q(z;c) Eq (z;c)) (15) AbasicfactisthatanexponentialpolynomialE(z) of(13)belongstotheclasshbifand only if it has no zeros on the closed upper halfplane C + R = {z = x+iy C y 0} (see [11, Chapter VII, Theorem 6], for example) Therefore, there exists a Hamiltonian H(a) of a canonical system corresponding to E(z) if E(z) has no zeros on C + R
3 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 3 Beyond such results standing on a general theory, we show that there exists a symmetric matrix valued function H(a) of a quasicanonical system corresponding to an exponential polynomial E(z) if E(z) has no zeros on the real line at least It is constructed explicitly as follows We define two lower triangular matrices E + and E of size 2g +1 by E ± = E ± (C) := C g C (g 1) C g C 0 C 1 C g and define square matrices J n of size 2g +1 by J n = J n (2g+1) := C ±(g 1) C ±(g 2) C g C ±g C ±(g 1) C 0 C (g 1) C g J (n) 0 0 0, J(n) := for 1 n 2g, where J (n) is the antidiagonal matrix of size n with 1 s on the antidiagonal line By using the above matrices, we define the numbers n by n = n (C) := det(e+ (C)+E (C)J n ) det(e + (C) E (C)J n ) 1 1 (16) for 1 n 2g and 0 = 0 (C) := 1 We write n = if det(e + E J n ) = 0 In addition, we define the realvalued locally constant function γ : [1,q g ) R by γ(a) = γ(a;c) := (C) n (C) for q ()/2 a < q n/2 (17) Then, we obtain the following results for the inverse problem associated with the exponential polynomial (13) Theorem 11 Let q > 1 and C R R 2g 1 R Let E(z) be the exponential polynomial defined by (13) Define the square matrix D n (C) := C g C g+1 C g C g 1 C g C g+ C g+n 2 C g C g n+1 C g n+2 C g C g C g 1 C g n+1 C g C g+1 C g+ C g C g n+2 C g C g+n 2 C g C g C g of size 2n Suppose that E(z) has no zeros on the real line and that detd n (C) 0 for every 1 n 2g Then, (1) det(e + ±E J n ) 0 for every 1 n 2g, (2) the pair of functions (A(a,z),B(a,z)) defined in (213) below satisfies a d A(a,z) 0 1 A(a,z) = z H(a) (a [1,q da B(a, z) 1 0 B(a, z) g ), z C) (18)
4 4 M SUZUKI together with the boundary condition A(1,z) A(z) = B(1, z) B(z), lim aրq g A(a,z) = B(a, z) E(0), (19) 0 where A(z) and B(z) are functions of (14) and (15), respectively, and γ(a;c) 1 0 H(a) = H(a;C) := (110) 0 γ(a; C) We mention another way of constructing (γ(a),a(a,z),b(a,z)) in Section 6 The positive definiteness of the Hamiltonian H(a) in (110) is characterized by the following tobeexpected way Theorem 12 Let q > 1 and C R R 2g 1 R Let E(z) be the exponential polynomial defined by (13) and let H(a) be the matrix valued function defined by (110) (1) Suppose that E(z) belongs to the class HB Then H(a) is welldefined and is positive definite for every 1 a < q g Hence the quasicanonical system attached to (18) and (19) is a canonical system and (A(a,z),B(a,z))/E(0) is its solution (2) Suppose that H(a) is welldefined and is positive definite for every 1 a < q g Then E(z) belongs to the class HB As mentioned before, E(z) of (13) belongs to the class HB if and only if it has no zeros on C + R On the other hand, the following conditions are equivalent to each other by 0 = 1 and definitions (16), (17) and (110): (i) H(a) is positive definite for every 1 a < q g, (ii) γ(a) > 0 for every 1 a < q g, (iii) 0 < n < for every 1 n 2g Therefore, we obtain the following corollary: Corollary 13 An exponential polynomial E(z) of (13) has no zeros on C + R if and only if 0 < n < for every 1 n 2g The converse of Theorem 11 is the direct problem for quasicanonical systems (11) with the Hamiltonians of the form (110) It is easier than the inverse problem, because the Hamiltonians have a simple form: H(a) = diag(γ(a) 1,γ(a)), and γ(a) is a (special) locally constant function Theorem 14 Let q > 0, g Z >0, and let γ(a) be a locally constant function on [1,q g ) such that γ(a) = γ n R for q ()/2 a < q n/2 (1 n 2g) Then the quasicanonical system (11) with H(a) = diag(γ(a) 1,γ(a)) on [1,q g ) has the unique solution (A(a, z), B(a, z)) whose components have the forms [ (q m) iz q m) ] iz A(a,z) = α m (n) +(, a a +n [ (q ) m iz ) ] (111) q m iz ib(a,z) = β m (n) ( a a +n for q ()/2 a < q n/2 and 1 n 2g, where α ± m(n) and β m(n) ± are real constants depending only on the set {γ n } of values of γ(a) Therefore, for q ()/2 a < q n/2, E(a,z) := A(a,z) ib(a,z) is the exponential polynomial [ ( ) q m iz ) ] q m iz E(a,z) = (α m (n)+β m (n)) +(α m (n) β m (n))( a a +n Moreover, E(a,0) = 1 and E(a,z) has no real zeros for any fixed 1 a q g In particular, each Hamiltonian of the form (110) with (17) yields an exponential polynomial E(1,z) = A(1,z) ib(1,z) having no real zeros
5 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 5 Theorem 14 does not guarantee that E(1,z) has the form (13) In fact, H(a) on [0,q) with γ(a) = 1 for 1 a < q 1/2 and γ(a) = 1 for q 1/2 a < q yields the constant function E(1,z) = 1, and H(a) on [0,q) with γ(a) = 1 for 1 a < q yields the function E(1,z) = q iz In these examples, E(1,z) does not have the form (13) A sufficient condition for which E(1,z) has the form (13) is the following Theorem 15 Under the notation of Theorem 14, E(1, z) is an exponential polynomial of the form (13) if γ n > 0 for every 1 n 2g and γ 1 1 If the assumption of Theorem 15 is satisfied, it is proved that E(1,z) belongs to the class HB in a way similar to the proof of Theorem 12 (2) in Section 52 On the other hand, γ n > 0 for every 1 n 2g and γ 1 1 if we start from the exponential polynomial of the form (13) by Theorem 12 and (58) below Hence, as a consequence of the above theorems, the exponential polynomials (13) belonging to the class HB are characterized in terms of positive definiteness of Hamiltonians Now we turn to an application of Corollary 13 A nonzero polynomial P(x) = c 0 x n +c 1 x + +c x+c n with real coefficients is called a selfreciprocal polynomial of degree n if c 0 0 and it satisfies the selfreciprocal condition P(x) = x n P(1/x), or equivalently, c 0 0 and c k = c n k for every 0 k n The roots of a selfreciprocal polynomial either lie on the unit circle T = {z C z = 1} or are distributed symmetrically with respect to T Therefore, a basic problem is to finda nice condition on coefficients of a selfreciprocal polynomial for which all its roots lie on T There are quite many results for this problem in literatures For example, see books of Marden[14], Milovanović Mitrinović Rassias [15], Takagi [24, Section 10] and the survey paper of Milovanović Rassias [16] for several systematic treatments on roots of polynomials As an application of Corollary 13, we study roots selfreciprocal polynomials of even degree The restriction on the degree is not essential, because if P(x) is a selfreciprocal polynomial of odd degree, then there exists a selfreciprocal polynomial P(x) of even degree and an integer r 1 such that P(x) = (x+1) r P(x) In contrast, the realness of coefficients is essential We denote by P g (x) a selfreciprocal polynomial of degree 2g of the form g 1 P g (x) = c k (x (2g k) +x k )+c g x g, c 0 0 (112) k=0 and identify the polynomial with the vector c = (c 0,c 1,,c g ) R R g consisting of its coefficients For a vector c R R g and a real number q > 1, we define the numbers δ n (c) (1 n 2g) by δ n (c) := det(e+ (C)+E (C)J n ) det(e + (C) E (C)J n ) where C is the vector defined by C m T m := c g m (1 logq m )T m + m=0 { 1 if n is even, glogq if n is odd, (113) c g m (1+logq m )T m (114) Theorem 16 Let g Z >0, q > 1 Let c R R g be coefficients of a selfreciprocal polynomial P g (x) of the form (112) Define C R R 2g 1 R by (114) Then, for every 1 n 2g, it is independent of q whether det(e + ±E J n ) is zero, and the numbers δ n (c) of (113) are independent of q if det(e + ±E J n ) is not zero Moreover, a necessary and sufficient condition for all roots of P g (x) to be simple roots on T is that 0 < δ n (c) < for every 1 n 2g Remark As a function of indeterminate elements (c 0,,c g ), δ n (c) is a rational function of (c 0,,c g ) over Q m=1
6 6 M SUZUKI Thecriterion of Theorem 16 itself may not benew, becauseitseems that thequantity δ n (c)intheorem16probablyessentially coincidewithquantities inaclassical theoryon the studying of the roots of polynomials in Section 75 The author has not yet thought rigorously whether δ n (c) and the classical quantities indeed have the same meaning However, at least, the theory of the presented paper provides a new bridge between the studying on the roots of polynomials and the theory of canonical systems The author believe that this new connection deserve to be recorded in a literature In order to deal with the case that all roots of P g (x) lie on T but P g (x) may have a multiple zero, we modify the above definition of δ n (c) as follow For a vector c R R g and real numbers q > 1, ω > 0, we define the numbers δ n (c;q ω ) (1 n 2g) by δ n (c;q ω ) = det(e+ (C)+E (C)J n ) 1 if n is even, det(e + (C) E (C)J n ) q gω q gω (115) q gω +q gω if n is odd, where C is the vector defined by C m T m := c g m q mω T m + m=0 c g m q mω T m (116) Theorem 17 Let g Z >0, q > 1 Let c R R g be coefficients of a selfreciprocal polynomial P g (x) of the form (112) Then a necessary and sufficient condition for all roots of P g (x) to lie on T is that 0 < δ n (c;q ω ) < for every 1 n 2g and ω > 0 m=1 Furthermore, the quantities δ n (c) and δ n (c;q ω ) are related as follows: Theorem 18 Let δ n (c) and δ n (c;q ω ) be as above Then lim δ n(c;q ω ) = δ n (c) q ω ց1 as a rational function of c = (c 0,,c g ) over Q Suppose that all roots of a selfreciprocal polynomial (112) lie on T and are simple Then we have where implied constants depend only on c δ n (c;q ω ) = δ n (c)+o(logq ω ) as q ω ց 1, Finally, we comment on the reason why we use a multiplicative variable on R >0 and the right endpoint for the initial value in (11) It comes from the author s personal motivation for the work of this paper: it is the inverse problem for entire functions of the class HB obtained by Mellin transforms of functions on [1, ) This problem was stimulated by Burnol [1], in particular Sections 5 8, and was partially treated in [23] under the number theoretic setting In [1], Burnol often use a multiplicative variable on R >0 and the left endpoint (= + ) corresponds to the initial value In order to apply the method of [1] to the above inverse problem, the author noted on the formulas and X 1 1 f(x)x izdx x = lim X X 1 f(x)x izdx x f(x)x izdx logx/logq x = lim 1 f(q n )q inz q 1+ logq Here exponential polynomials appear in the left hand side of the second formula If we fix X > 1, q > 1 and put 2g = logx/logq, we obtain an exponential polynomial of degree 2g Recall that the results of the paper corresponds the exponential polynomial to the canonical system on [1,q g ) Then, if we imagine the limiting situation q 1 and n=0
7 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 7 X +, it is expected (but not rigorous) that a canonical system corresponds to the Mellin transform 1 f(x)x iz dx x should be a system on [1, ) To realize the above heuristic discussion to the original motivation, the right endpoint may be more useful and convenient than the left endpoint for the initial value, although the left endpoint is useful for the initial value in the usual theory of canonical systems Anyway, the author believe that the above number theoretic aspect of de Branges theory is an important one The paper is organized as follows In Section 2, we describe the outline of the proof of Theorem 11 In Section 3, we prepare several lemmas and notations in order to prove statements of Section 2 and Theorem 11 in Section 4 In Section 5, we prove Theorems 12, 14 and 15 In Section 6, we mention a way of constructing (γ(a),a(a,z),b(a,z)) which is different from the way of Section 2 and 4 In Section 7, we prove Theorems 16, 17 and 18 and compare these with classical results on the roots of selfreciprocal polynomials Acknowledgments The author thanks Shigeki Akiyama for suggesting the book of Takagi ([24, Section 10]) The author also thanks the referee for a number of helpful suggestions for improvement in the article and for careful reading In particular, Theorems 14 and 15 were added by the suggestion of the referee This work was supported by KAKENHI (GrantinAid for Young Scientists (B)) No and No Outline of the proof of Theorem Hilbert spaces and operators Let q > 1 and let L 2 (T q ) be the completion of the space of (2π/logq)periodic continuous functions on R with respect to the L 2 norm f 2 L 2 (T := f,f q) L 2 (T q),where f,g L 2 (T q) := logq 2π I q f(z)g(z)dz andi q = [0,2π/(logq)) Every f L 2 (T q ) has the Fourier expansion f(z) = k Z u(k)qikz with {u(k)} k Z l 2 (Z) and f 2 L 2 (T q) = k Z u(k) 2, where l 2 (Z) is the Hilbert space of sequences {u(k) C k Z} satisfying k Z u(k) 2 < ([18, Chapter 4]) For 0 < a q Z/2 = {q n/2 n Z}, we define the vector space V a := { φ(z) = a iz f(z)+a iz g(z) f, g L 2 (T q ) } of functions of z R As a vector space, V a is isomorphic to the direct sum L 2 (T q ) L 2 (T q ), since a iz f(z) + a iz g(z) = 0 if and only if (f,g) = (0,0) In fact, if one of f and g is zero and a iz f(z)+a iz g(z) = 0, the another one should be zero On the other hand, if f 0, g 0 and a iz f(z)+a iz g(z) = 0, we have a 2iz = f(z)/g(z), and hence a q Z/2 The maps p 1 : (a iz f(z)+a iz g(z)) a iz f(z) and p 2 : (a iz f(z)+a iz g(z)) a iz g(z) are projections from V a to the firstand the second components of the direct sum, respectively In addition, we define the inner product in V a by φ 1,φ 2 = f 1,g 1 L 2 (T q) + g 1,g 2 L 2 (T q), where φ j (z) = a iz f j (z) + a iz g j (z) (j = 1,2) Then, V a forms a Hilbert space with respect to the inner product and is isomorphic to the (orthogonal) direct sum L 2 (T q ) L 2 (T q ) of Hilbert spaces ([9, Chapter 5]) We put X(k) := (q k /a) iz, Y(l) := X(l) 1 = (q l /a) iz (21) for k,l Z and a > 0 We regard X(k) and Y(l) as functions of z, functions of (a,z) or symbols depending on the situation For a fixed 0 < a q Z/2, the countable set consisting of all X(k) and Y(l) is linearly independent over C as a set of functions of z, since the linear dependence of {X(k), Y(l)} k,l Z implies the existence of a nontrivial
8 8 M SUZUKI pair of f,g L 2 (T q ) satisfying a iz f(z)+a iz g(z) = 0 Using these vectors, we have { V a = φ = u(k)x(k)+ } v(l)y(l) {u(k)} k Z, {v(l)} l Z l 2 (Z) k Z l Z and L 2 (T q ) L 2 (T q ), { } p 1 V a = u(k)x(k) {u(k)} k Z l 2 (Z) L 2 (T q ), k Z { } p 2 V a = v(l)y(l) {v(l)} l Z l 2 (Z) L 2 (T q ) l Z φ,φ = u(k)u (k)+ v(l)v (l) (φ,φ V a ), k Z l Z φ 2 = φ,φ = u(k) 2 + v(l) 2 (φ V a ), (22) k Z l Z where φ = k Z u(k)x(k)+ l Z v(l)y(l) and φ = k Z u (k)x(k)+ l Z v (l)y(l) On the other hand, we have φ 2 = logq p 1 φ(z)p 1 φ(z)dz + logq p 2 φ(z)p 2 φ(z)dz (23) 2π I q 2π I q for φ V a, since logq (a ±iz f 1 (z))(a 2π ±iz f 2 (z))dz = u 1 (k)u 2 (k) = f 1,f 2 L 2 (T q) I q k Z for f j (z) = k Z u j(k)q ikz L 2 (T q ) (j = 1,2) Note that, for φ V a, p 1 φ and p 2 φ are not periodic functions of z, but the integrals I p jφ(z)p j φ (z)dz (j = 1,2, φ,φ V a ) are independent of the intervals I = [α,α+2π/logq] (α R) We write φ V a as φ(z) (resp φ(a,z)) to emphasize that φ is a function of z (resp (a,z)) If we regard X(k) and Y(l) as symbols, V a is abstract Hilbert space endowed with the norm defined by (22) and is isomorphic to l 2 (Z) l 2 (Z) For each nonnegative integer n, we define the closed subspace V a,n of V a by { } V a,n = φ n = u n (k)x(k) + v n (l)y(l) {u n(k)} k=0, {v n(l)} l= l2 (Z) k=0 k= and denote by P n the projection operator V a V a,n Then, for the conjugate operator P n := JP nj : V a V a of P n by the involution J : φ(z) φ( z), we have P n φ = JP njφ = u(k)x(k) + v(l)y(l) (φ V a ), k= because J acts on V a as exchanging of coefficients between X(k) and Y(l): Jφ = v(k)x(k) + u(l)y(l) (φ V a ) k= l= Therefore, P n maps V a,n into V a,n for every nonnegative integer n In addition, JP n also maps V a,n into V a,n for every nonnegative n, because l=0 l=0 JP n φ n = v n (l)x(l)+ u n (k)y(k) (24) for φ n V a,n Note that P 0 Va,0 = JP 0 Va,0, since P 0 J V a,0 = 0 by definition k=0
9 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 9 For an exponential polynomial E(z) of (13), we have E (z) = C m q imz = E q (z;c ) for C := (C g,c (g 1),,C g ) R R 2g 1 R by a simple calculation By using E(z) and E (z), we define two multiplication operators E : φ(z) E(z)φ(z), E : φ(z) E (z)φ(z) (25) These operators map V a into V a, because E and E are expressed as E = C m T m, E = C m T m by using shift operators T m : V a V a (m Z) defined by T m v = u(k)x(k +m)+ v(l)y(l m) k= l= Both E and E are bounded on V a, since E g C m T m M and E g C m T m M for M = max{c m g m g} Suppose that E(z) has no zeros on the real line Then the operator E is invertible on V a (Lemma 41(1)) Thus the operator is welldefined We have Θ := E 1 E (26) (Θφ)(z) = Θ(z)φ(z) for Θ(z) := E (z) E(z), φ V a (27) and Θ(V a,n ) V a,n for each nonnegative integer n (Lemma 41(2)) 22 Quasicanonical systems associated with exponential polynomials Under the above settings, a quasicanonical system associated with an exponential polynomial E(z) of (13) is constructed starting from solutions of the set of linear equations (I±ΘJP n )φ ± n = ΘX(0) (φ ± n V a,n, 0 n 2g), (28) where I is the identity operator The set of 4g+2 equations (28) is a discrete analogue of the (right) Mellin transform of differential equations (117a) and (117b) in [1, Section 6] (see also [23, Section 4]) Under the assumption that both I±ΘJP n are invertible on V a,n for every 0 n 2g, that is, (I±ΘJP n ) 1 exists as a bounded operator on V a,n, we define A n (a,z) := 1 2 (I+J)Eφ+ n (a,z), B n (a,z) := 1 2i (I J)Eφ n (a,z) (29) by using unique solutions of (28) Then they are entire functions of z and are extended to functions of a on (0, ) (by formula (45)) Here we define functions A (a,z) and B (a,z) of (a,z) [1,q g ) C by A (a,z) : = A n (a,z) B (a,z) : = B n(a,z) for q ()/2 a < q n/2 (210) The function A (a,z) is discontinuous at a [1,q g ) q Z/2 in general, because A n(q n/2,z) = A n+1(q n/2,z)
10 10 M SUZUKI may not be hold It is similar about B (a,z) However, we will see that α n = α n (C) := A (q()/2,z) A n(q ()/2,z), β n = β n (C) := B (q()/2,z) B n(q ()/2,z) (211) are independent of z for every 1 n 2g (Proposition 48) Therefore, we obtain functions A(a,z) and B(a,z) of (a,z) [1,q g ) C which are continuous for a and entire for z by the modification for 1 n 2g and A n (a,z) : = α 1 α n A n(a,z), B n (a,z) : = β 1 β n B n(a,z) A(a,z) : = A n (a,z) B(a,z) : = B n (a,z) (212) for q ()/2 a < q n/2 (213) Moreover, (A(a, z), B(a, z)) satisfies the system (18) endowed with the boundary condition (19) for the locally constant function γ(a) defined by and γ n = γ n (C) := α 1(C) α n (C) β 1 (C) β n (C) (214) γ(a) := γ(a;c) := γ n if q ()/2 a < q n/2 (215) (Proposition 45 and Theorem 418) This γ(a) is equal to the function of (17) (Proposition 48) Therefore, (A(a, z), B(a, z))/e(0) is a solution of a quasicanonical system on [1,q g ) if E(0) 0 and α n,β n 0 for every 1 n 2g As a summary of the above argument, we obtain Theorem 11 See Section 4 for details In order to prove Theorems 12, 14 and 15, we use the theory of de Branges spaces which is a kind of reproducing kernel Hilbert spaces consisting of entire functions Roughly, the positivity of H(a) corresponds to the positivity of the reproducing kernels of de Branges spaces See Section 5 for details 3 Preliminaries 31 Identities for matrices The following basic facts of linear algebra are used often in later sections Lemma 31 Let A and B be square matrices of the same size Then, A B det = det(a+b)det(a B) B A Lemma 32 Let A, B, C, D be square matrices of the same size If deta 0, detd 0 and det(d CA 1 B) 0, 1 [ A B (A BD = 1 C) 1 A 1 B(D CA 1 B) 1 ] C D D 1 C(A BD 1 C) 1 (D CA 1 B) 1 Lemma 33 (Laplace s expansion formula) Let A be a square matrix of size n Let i = (i 1,i 2,,i k ) (resp j = (j 1,j 2,,j k )) be a list of indices of k rows (resp columns), where 1 k < n and 0 i 1 < i 2 < i k < n (resp 1 j 1 < j 2 < < j k < n) Denote by A(i,j) the submatrix of A obtained by keeping the entries in the intersection of any row and column that are in the lists Denote by A c (i,j) the submatrix of A obtained
11 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 11 by removing the entries in the rows and columns that are in the lists By Laplace s formula for determinants, we get deta = j ( 1) i + j deta(i,j)deta c (i,j), where i = i 1 +i 2 + +i k, j = j 1 +j 2 + +j k, and the summation is taken over all ktuples j = (j 1,j 2,,j k ) for which 1 j 1 < j 2 < < j k < n Proof See [6, Chap IV], for example 32 Definition of special matrices, I We define several special matrices for a convenience of later arguments The matrix E 0 ± is the square matrix of size 8g defined by E 0 ± = E± 0 (C) := e ± 0 (C) 0 0 e± 0 (C), where e ± 0 (C) are lower triangular matrices of size 4g defined by C g C (g 1) C g e ± 0 = C e± 0 (C) := ±(g 1) C (g 1) C g C ±g C (g 1) C g 0 C±(g 1) C±g C ±(g 1) C (g 1) C g 0 0 C ±g C ±(g 1) C (g 1) C g Replacing the column at the left edge of e 0 by the zero column vector, C g C (g 1) Cg e 1 = e 1 (C) := C (g 1) C(g 1) C g C g C(g 1) C g 0 C (g 1) 0 C g C (g 1) C(g 1) C g C g C (g 1) C (g 1) C g We define and E nj = E nj(c) := E 0 J = E 0 J(C) := 0, E 0,1 J = E 0,1J(C) := 0 0 e 1,n (C) e 2,n (C) 0, E n,1 J = E n,1j(c) := 0 e 3,n (C) e 2,n (C) 0
12 12 M SUZUKI by setting e 1,n = e 1,n (C) := e 2,n = e 2,n (C) := C g C (g 1) C g 0 C (g 1) C (g 1) 0 C g C (g 1) C g C g C (g 1) C g C (g 1) C (g 1) 0 C g C (g 1) C g = = 2g n n 2g n 4g n 2g n 2g +n 2g +n, 2g n and e 3,n = e 3,n (C) := C g 0 0 C (g 1) 0 C g 0 C (g 1) C (g 1) 0 0 C g 0 C (g 1) C g = 2g n n 2g 2g +n 2g n for 1 n 2g, where the righthand sides mean the size of each block of matrices in middle terms In addition, we define square matrices e 0 of size 4g by e 0 = e 0 (C) := C g C g 1 C (g 1) C g C g 1 C g C g C (g 1) C g C g 1 C (g 1) C g C g 1 C g C (g 1) C (g 1) C g C g
13 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 13 so that J (8g) E + 0 = 0 e 0 e 0 0 Replacing n columns from the left edge of e 0 by zero column vectors, C g C g 1 C (g 1) C g C g 1 C g C g C (g 1) e 0,n = e 0,n (C) : = 0 Cg 1 C (g 1) C g C g C (g 1) C g = n 4g n 4g We denote by I (m) the identity matrix of size m and by J n (m) of size m: J n (m) = We also use the notation J (n) 0 0 0, J (n) = 1 the following square matrix 1 χ n = t[ ] = the unit column vector of length n For a square matrix M of size N, we denote by [M] տn (resp [M] րn ) the square matrix of size n( N) in the top left (resp the top right) corner, and denote by [M] ցn (resp [M] ւn ) the square matrix of size n( N) in the lower right (resp the lower left) corner: M = [ [M]տn = [M] ցn n [M] ւn n [M] րn 33 Definition of special matrices, II We define the square matrix P k (m k ) of size 2k+2 endowed with the parameter m k and the (2k+2) (2k+4) matrix Q k for every nonnegative integer k as follows For k = 0,1, P 1 (m 1 ) := P 0 := [ m 1 ], Q 0 :=, Q 1 := [ ], ] For k 2, we define P k (m k ) and Q k blockwisely as follows V + k 0 W + P k (m k ) := 0 V k 0 k, Q k := 0 W k 0I k m k 0I k 0 k,k+2 0 k,k+2,
14 14 M SUZUKI where 0I k := 0 k,1 I k, mk 0I k = 0 k,1 m k I k, 0k,l is the k l zero matrix, I k is the identity matrix of size k, 1 0 V + k := 1 1 ( ) k+3 (k+1) if k is odd, 2 1 V k := := := ( k +2 2 ( k +1 ) (k +1) if k is even, 2 ) (k +1) if k is odd, ( ) k+2 (k+1) if k is even, 2 and matrices W ± k are defined by adding column vectors t (1 0 0) to the rightside end of matrices V ± k : W + k := 1 1 ( ) k+3 (k+2) if k is odd, ( ) := k +2 (k +2) if k is even, W k := := Lemma 34 For k 1, we have where ε 2j+1 = detp k (m k ) = { { ( k +1 ( k+2 ε 2j+1 2 j m j+1 2j+1 if k = 2j +1, ε 2j 2 j m j 2j if k = 2j, 2 2 ) (k +2) if k is odd, ) (k+2) if k is even { +1 j 2,3 mod4 1 j 0,1 mod4, ε +1 j 0,1 mod4 2j = 1 j 2,3 mod4 In particular, P k (m k ) is invertible if and only if m k 0
15 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 15 Proof This is trivial for k = 1 Suppose that k = 2j and write P k (m k ) as (v 1 v 2k+2 ) by its column vectors v l At first, we make the identity matrix I k+2 at the leftupper corner by exchanging the columns v (k+5)/2,,v k+1 and v k+2,,v (3k+3)/2 so that detp k (m k ) = det(v 1 v (k+3)/2 v k+2 v (3k+3)/2 v (k+5)/2 v k+1 v (3k+5)/2 v 2k+2 ) Then, by eliminating every 1 and m k under I k+2 of the leftupper corner, we have Ik+2 detp k (m k ) = det 0 k,k+2 Z k Here Z k is the k k matrix Z k = Zk,1 Z k,2 0 j+1,j Z k,3 for which Z k,1 is the j j antidiagonal matrix with 1 on the antidiagonal line and m k j (j +1) Zk,2 m k = Z k,3 m k 2m k (j +1) (j +1) 2mk The above formula of detp k (m k ) implies the desired result The case of even k is proved in a way similar to the case of odd k Lemma 35 We have P k (m k ) 1 Q k = Mk,1 M k,2 M k,3 M k,4 (2k +2) (2k +4) (31) for every 0 k 2g, where M k,1, M k,2, M k,3, M k,4 are (k+1) (k+2) matrices defined by M0,1 M 0, M1,1 M =, 1,2 = M 0,3 M 0, M 1,3 M 1, /m and M k,1 = if k 3 is odd, = if k 2 is even,
16 16 M SUZUKI M k,4 = if k 3 is odd, = if k 2 is even, M k,2 = m k if k 3 is odd, = m k if k 2 is even, M k,3 = m k 2 if k 3 is odd, = 1 2m k 1 1 if k 2 is even Proof According to definitions of P k (m k ), Q k and M k,l (1 l 4), the identity Mk,1 M P k (m k ) k,2 = Q M k,3 M k k,4 is checked by an elementary way
17 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 17 4 Proof of Theorem 11 In this section, we complete the proof of Theorem 11 by showing each statement in Section 2 We fix g Z >0, q > 1 and C R R 2g 1 R throughout this section Lemma 41 Let E be the multiplication operator defined by (25) for E(z) of (13) Suppose that E(z) has no zeros on the real line Then (1) E is invertible on V a Moreover, for the operator Θ of (26), we have (2) Θ(V a,n ) V a,n for each nonnegative integer n Proof (1) The multiplication by 1/E(z) maps a sum of X(k) (resp Y(l)) to another sum of X(k) (resp Y(l)) as well as multiplication by E(z), since 1 E(z) = q igz 2g m=0 C m gq = C m q imz imz for some sequence { C m } of real numbers Therefore, it is sufficient to prove that E is invertible on p 1 V a and p 2 V a We have 1/E(z) L (T q ) by the assumption Therefore, multiplication by1/e(z) definesaboundedoperatore 1 onp i V a withthenorm E 1 = 1/E L (T q) for i = 1,2 (2) We have Θ(z) = 1 on the real line by definition (27) Thus, multiplication by Θ(z) defines a bounded operator on V a with the norm Θ = 1 We can write 2g m=0 Θ(z) = C g mq imz 2g m=0 C m gq = R m q imz (R m R) imz m=0 Therefore, Θ(z)φ n (z) = R m u n (k)x(k +m)+ R m v n (l)y(l m) = m=0 k=0 ( m=0 k=m This shows Θ(V a,n ) V a,n Let us consider the equations ) R m u n (k m) X(m)+ m=g m=0l= 0 m= (+m l= ) R m v n (l m) Y(m) (E±E JP n )φ ± n = E X(0) (φ ± n V a,n, 1 n 2g) (41) instead of (28) They are equivalent to (41) if E is invertible Lemma 42 Let 0 < a q Z/2 Let D n (C) be matrices define in Theorem 11 Suppose that E(z) of (13) has no zeros on the real line and that detd n (C) 0 for every 1 n 2g Then ΘJP n defines a compact operator on V a,n for each 0 n 2g and the resolvent set of ΘJP n Va,n contains both ±1 In particular, I ±ΘJP n are invertible on V a,n and (41) has unique solutions in V a,n for each 0 n 2g Remark If a q Z/2 (0, ), I±ΘJP n may not be invertible Proof It is trivial for n = 0, since P 0 Va,0 = 0 On V a,n for n 1, the HilbertSchmidt norm of ΘJP n is finite In fact, ΘJP n Va,n 2 2 = ΘJP n X(k) 2 + ΘJP n Y(l) 2 k=0 JP n X(k) 2 + k=0 l= l= JP n Y(l) 2 k=0 1 < 1+ l=0
18 18 M SUZUKI by(22)and(24) Therefore, ΘJP n Va,n isahilbertschmidtoperatoronv a,n, andhence it is compact Thus, either λ C\{0} is an eigenvalue of ΘJP n Va,n or an element of the resolvent set of ΘJP n Va,n Assume that ΘJP n φ n = ±φ n Then ΘJP n φ n = φ n Because Θ and p i (i = 1,2) commute (see the proof of Lemma 41(2)), ΘJP n φ n 2 = logq 2π by (23) While, φ n 2 = logq 2π = logq 2π k=0 by (23) Thus, it must be In this case, and Θ(z)p 1 JP n φ n (z) 2 dz + logq Θ(z)p 2 JP n φ n (z) 2 dz I q 2π I q p 1 JP n φ n (z) 2 dz + logq p 2 JP n φ n (z) 2 dz I q 2π I q = u(k) 2 + v(l) 2 E JP n φ n = ±Eφ n = ± l=0 p 1 φ n (z) 2 dz + logq I q 2π k=0 k=0 I q p 2 φ n (z) 2 dz = φ n = u(k)x(k) + v(l)y(l) l=0 C m v(k)x(k +m)+ k=0 C m u(k)x(k +m)± l=0 l=0 u(k) 2 + k=0 l= C m u(l)y(l m) C m v(l)y(l m) v(l) 2 Comparing the coefficients of {X(k)} g k g+ and {X(k)} g k g+ in the equality E JP n φ n ±Eφ n = 0, yields where A ± = A ± t[ u n (0) u n () v n (0) v n () ] = 0, (42) ±C g C g ±C g+1 ±C g C g 1 C g ±C g+ ±C g+n 2 ±C g C g n+1 C g n+2 C g ±C g ±C g 1 ±C g n+1 C g C g+1 C g+ ±C g ±C g n+2 C g C g+n 2 ±C g C g Here deta ± 0 by detd n (C) 0 Therefore (42) has no nontrivial solutions, thus φ n = 0 Consequently, both ±1 are not eigenvalue, and hence they belong to the resolvent set ofθjp n Va,n whichmeansthati±θjp n Va,n existsandboundedonv a,n In the remaining part of the section, we assume that E(z) 0 for z R and detd n (C) 0 for 1 n 2g so that both E±E JP n are invertible on V a,n for every 0 n 2g (cf Lemma 41 and Lemma 42)
19 Let φ ± n = k=0 u± n (k)x(k) + Using coefficient of φ ± n and putting INVERSE PROBLEMS FOR CANONICAL SYSTEMS 19 l= v± n (l)y(l) be solutions of (41) for 0 n 2g v ± n(n) = v ± n(n+1) = v ± n(2g 1) = 0 if 0 n 2g 1, we define the column vectors Φ ± n of length 8g by Φ ± n = t[ u ± n (0) u± n (1) u± n (4g 1) v± n (2g 1) v± n (2g 2) v± n ( 2g)] Substituting φ ± n in (41) and then comparing coefficients of X(k) and Y(l), we obtain and (E + 0 ±E n J)Φ± n = E 0 χ 8g, (Φ ± n V a,n, 0 n 2g) (43) 2 j=0 C g j u ± n(k +j) = 0, for every K 4g and L 2g 1 2 j=0 C g j v ± n(l j) = 0 (44) Lemma 43 Let 0 n 2g We have det(e + 0 ±E nj) 0 if E±E JP n is invertible Proof For fixed n, equations (41) have unique solutions φ ± n, if E ± E JP n are invertible On the other hand, the solutions φ ± n correspond to solutions of the pair of linear equations (43) and (44) Therefore, det(e 0 + ±E nj) 0 On the other hand, Eφ ± n = E X(0) E JP n φ ± n ( = C m (X(m) v ± n (k)x(k +m)+u ± n (k)y(k m))) by (41) Therefore, we can write k=0 g+ Eφ ± n = ) (p ± n (k)x(k)+q± n (k)y(k) k= g (45) for some real numbers p ± n(k) and q n(k) ± Hence Eφ ± n(a,z) are extended to smooth functions of a on (0, ) by the righthand side of (45) We use the same notation for such extended functions We put p ± n(k) = q n(k) ± = 0 for every g+n k 3g 1 if 0 n 2g 1 and define the column vectors Ψ ± n of length 8g by Ψ ± n = t[ p ± n ( g) p± n ( g +1) p± n (3g 1) q± n (3g 1) q± n (3g 2) q± n ( g)] Then, we have by g+ Eφ ± n = ) (p ± n (k)x(k)+q± n (k)y(k) = k= g k=0 C m u n (k)x(k +m)+ Ψ ± n = E+ 0 Φ± n (46) l= C m v n (l)y(l m) Lemma 44 Let 0 n 2g We have p ± n(k)±q n(k) ± = 0 if g +1 k g + or g k g+, and p ± n(g)±q n(g) ± = C g Therefore, (I±J)Eφ ± n = (p ± n(k)±q n(k))(x(k) ± ±Y(k)) k= g+n
20 20 M SUZUKI Proof By definition of matrices, the equality (E 0 + ±E nj)φ ± n = E0 χ 8g is written as e + 0 ±e 1,n ±e 2,n e + Φ± n = E 0 χ 8g (47) 0 On the other hand, e + 0 ±e 0 ±e 0 e + 0 Φ± n = p ± n( g)±q n( g) ± p ± n( g+1)±q n( g ± +1) p ± n(3g 1)±q n(3g ± 1) ±(p ± n(3g 1)±q n(3g ± 1)), ±(p ± n(3g 2)±q n(3g ± 2)) ±(p ± n( g)±q ± n( g)) since E + 0 ±J(8g) E + 0 = e + 0 ±e 0 ±e 0 e + 0 by definition of e 0 and p ± n ( g)±q± n ( g) p ± n ( g +1)±q± n ( g +1) (I (8g) ±J (8g) )E 0 + Φ± n = (I (8g) ±J (8g) )Ψ ± n = p ± n (3g 1)±q± n (3g 1) ±(p ± n (3g 1)±q± n (3g 1)) ±(p ± n (3g 2)±q± n (3g 2)) ±(p ± n ( g)±q± n ( g)) by (46) In addition, e + 0 ±e 0 ±e 0 e + 0 Φ± n = e + 0 ±e 0,n ±e 0 e + 0 Φ± n, since v n(n) ± = v n(n ± +1) = v n(2g ± 1) = 0 for 1 n 2g 1 by definition of Φ ± n Therefore, p ± n ( g)±q± n ( g) p ± n ( g +1)±q± n ( g +1) e + 0 ±e 0,n ±e 0 e + Φ± n = p ± n(3g 1)±q n(3g ± 1) ±(p ± n(3g 1)±q n(3g ± 1)) (48) 0 ±(p ± n (3g 2)±q± n (3g 2)) ±(p ± n( g)±q n( g)) ±
21 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 21 Here we find that 2g rows of both E 0 + e + ±E nj and 0 e 0,n e 0 e + with indices (2g+1,2g+ 0 2,,4g) and n rows of both E 0 + e + ±E nj and 0 e 0,n e 0 e + with indices (8g n+1,8g 0 n+2,,8g) have the same entries Therefore, by comparing 2g rows of (47) and (48) with indices (2g +1,2g +2,,4g), we have p ± n(g)±q ± n(g) = C g and p ± n(k)±q n(k) ± = 0 (g +1 k 3g 1) Similarly, by comparing n rows of (47) and (48) with indices (8g n + 1,8g n + 2,,8g), we have Hence we complete the proof p ± n(k)±q ± n(k) = 0 ( g k g+) We define the column vectors A n and B n of length 8g by A n = A n(c) := (I +J)Ψ + n = (I +J)E 0 + Φ+ n, B n = B n(c) := (I J)Ψ n = (I J)E + 0 Φ n, where I = I (8g) and J = J (8g), and define the row vector F(a,z) of length 8g by F(a,z) := [ X( g) X( g +1) X(g) 0 0 Y(g) Y(g 1) Y( g) ] Then, we obtain by (210) and Lemma 44 Proposition 45 We have a d (49) A n(a,z) = 1 2 F(a,z) A n, B n(a,z) = 1 2i F(a,z) B n (410) da A n (a,z) = zb n for every 0 n 2g (a,z), a d da B n (a,z) = za n (a,z) We provide two lemmas in order to prove Proposition 45 Lemma 46 We have for every 0 n 2g u + n(k) = u n(k) (0 k 4g 1), v n + (k) = v n (k) ( 2g k 2g 1) Proof We have Φ ± n = (E+ 0 ±E nj) 1 E0 χ 8g with E 0 + ±E n J = e + 0 ±e 1,n ±e 2,n e + 0 by (43) and definition of EnJ Applying Lemma 32 to A = D = P := e + 0, B = ±Q := ±e 1,n and C = ±R := ±e 2,n, we have [ Φ ± n = (P QP 1 R) 1 P 1 Q(P RP 1 Q) 1 ] P 1 R(P QP 1 R) 1 (P RP 1 Q) 1 E0 χ 8g This shows Lemma 46, since all 4g entries of E0 χ 8g with indices (4g+1,4g+2,,8g) are zero
22 22 M SUZUKI Lemma 47 We have p + n (k)+q+ n (k) = p n (k) q n (k) ( g k g +) for every 0 n 2g, where p ± n (k)±q± n (k) = 0 if g k g + or g +1 k g + by Lemma 44 Proof By (41), Eφ ± n = E X(0) E JP n φ ± n = C m X(m) Therefore, (I+J)Eφ + n = + On the other hand, (I J)Eφ n = = C m X(m) C m Y(m) C m X(m) C m Y(m)+ C m X(m) C m Y(m)+ C m v n ± (k)x(k +m) k=0 ( ( ( ( C m C m u ± n (l)y(l m) l=0 ) u + n(k)x(k +m)+c m v n(k)x(k + +m) k=0 C m k=0 ) u + n (l)y(l+m)+c m v n + (l)y(l+m) l=0 C m ( l=0 ) u n(k)x(k +m) C m vn(k)x(k +m) k=0 C m C m ( k=0 ) u n(l)y(l+m) C m vn(l)y(l+m) l=0 l=0 ) u + n(k)x(k +m)+c m v n(k)x(k + +m) k=0 C m k=0 ) u + n(l)y(l+m)+c m v n(l)y(l+m) + by Lemma 46 Comparing the righthand sides of the above formulas of (I ± J)Eφ ± n with formulas of (I±J)Eφ ± n in Lemma 44, we obtain Lemma 47 Proof of Proposition 45 By definition of X(k) and Y(l), (45), and Lemma 44, (I±J)Eφ ± n (a,z) = (p ± n (k)±q± n (k))((qk /a) iz ±(q k /a) iz ) k= g+n Therefore, the differentiability of A n(a,z) and Bn(a,z) with respect to a is trivial, and a d da (I+J)Eφ+ n(a,z) = iz (p + n(k)+q n(k))((q + k /a) iz (q k /a) iz ), a d 1 da i (I J)Eφ n(a,z) = z k= g+n k= g+n l=0 l=0 (p n(k) q n(k))((q k /a) iz +(q k /a) iz ) Applying Lemma 47 to the righthand sides, we obtain Proposition 45 by definition (29) Asmentionedbefore, A n (a,z)andb n (a,z)of(29)aresmoothfunctionsofaon(0, ) for every 0 n 2g However A (a,z) and B (a,z) of (210) may be discontinuous at a (0, ) q Z/2 Therefore, as mentioned in Section 2, the next aim is to make modifications so that A (a,z) and B (a,z) become continuous functions of a on [1,q g ) The essential part is the following proposition
23 INVERSE PROBLEMS FOR CANONICAL SYSTEMS 23 Proposition 48 Let 1 n 2g Suppose that E±E JP is invertible on V a, for every q (n 2)/2 < a < q n/2 and E±E JP n is invertible on V a,n for every q ()/2 < a < q n/2 Define α n and β n by (211) Then, we have det(e + +E J 1 ) if n = 1, α n = det(e + ) det(e + +E J 2 ) det(e + ) det(e + 0 ) det(e + 0 +E 1 J) if n = 2, det(e + +E J n ) det(e 0 + +E n 2 J) det(e + +E J n 2 ) det(e 0 + +E J) if 3 n 2g, det(e + E J 1 ) det(e + if n = 1, ) det(e + E J 2 ) det(e 0 + ) β n = det(e + ) det(e 0 + +E 1 J) if n = 2, det(e + E J n ) det(e 0 + E n 2 J) det(e + E J n 2 ) det(e 0 + E J) if 3 n 2g In particular, α n and β n depend only on C We prove Proposition 48 after preparing several lemmas Lemma 49 For every 0 n 2g, we have (411) (412) det(e + 0 +E nj) = det(e + 0 E nj) (413) Proof Suppose that n 1, since EnJ = 0 by definition We have E 0 + ±E n J = e + 0 ±e 1,n ±e 2,n e + 0 A 0 ±B C D 0 = ±B E A = 2g +n 4g 2n 2g +n 2g +n 4g 2n 2g +n by definition of e 1,n and e 2,n, where [ A = [e + 0 ] տ2g+n, B = [e [e + 2,n ] ւ2g+n, D = 0 ] ց2g n [e + 0 ] տ2g n and the righthand side means the size of each block of the matrix of the second line Applying Lemma 33 to the (4g 2g n) columns of E 0 + ± E nj with indices (2g + n+1,,4g), and then applying Lemma 33 to the (4g 2g n) rows of the resulting matrix with indices (2g +n+1,,4g), we obtain det(e 0 + A ±B ±E nj) = det(d)det ±B A By Lemma 31, det This equality shows (413) A ±B = det(d)det(a+b)det(a B) ±B A ]
24 24 M SUZUKI Lemma 410 We have adj(e + 0 ±e 0 J(4g) n )e 0 χ 4g = adj(e + 0 ±e 1 J(4g) n )e 0 χ 4g for every 1 n 2g, where adj(a) = (deta)a 1 is the adjugate matrix of A Proof Let w = t[ C g C g 0 0 ] be the column vector of length 4g Denote by u n,j and v n,j the column vectors of e + 0 ±e 0 J(4g) n e + 0 ±e 0 J(4g) n = [ u n,1 u n,4g ], e + 0 ±e 1 J(4g) n = [ v n,1 v n,4g ] By Cramer s formula, it is sufficient to prove that det[u n,1 u n,j 1 j w un,j+1 u n,4g ] = det[v n,1 v n,j 1 j w vn,j+1 v n,4g ] and e + 0 ±e 1 J(4g) n, respectively: (414) for every 1 j 2g By definition of e 0 and e 1, u n,j = v n,j for every j n, 1 j 4g and u n,n = v n,n +w Therefore, (414) is trivial if j = n If 1 j < n, det[u n,1 u n,j 1 j w un,j+1 u n,4g ] = det[v n,1 v n,j 1 j w vn,j+1 v n,4g ]+det[v n,1 j w n w vn,4g ] = det[v n,1 v n,j 1 j w vn,j+1 v n,4g ]+0 Hence we obtain (414) Lemma 411 We have and for 3 n 2g det(e + 0 ±e 1 J(4g) 1 ) = det(e + 0 ±e 1 J(4g) 2 ) = det(e + 0 ) det(e + 0 ±e 1 J(4g) n ) = det(e + 0 ±e 0 J(4g) n 2 ) Proof We have det(e + 0 ± e 1 J(4g) 1 ) = det(e + 0 ), since e 1 J(4g) 1 = 0 On the other hand, det(e + 0 ±e 1 J(4g) 2 ) = C 4g g = det(e+ 0 ), since e+ 0 is a lower triangular matrix of size 4g and For 3 n 2g, we have e 1 J(4g) 2 = e + 0 ±e 1 J(4g) n = C g 0 0 [e + 0 ±e 0 J(4g) n 2 ] տ4g 1 The 4gth row and the 4gth column of e + 0 ± e 0 J(4g) n 2 are [0 0 C g C g ] and t [0 0 C g ], respectively Therefore, det([e + 0 ±e 0 J(4g) n 2 ] տ4g 1) = C 1 g det(e+ 0 ±e 0 J(4g) n 2 ) Hence, we obtain the equality of Lemma 411
25 Lemma 412 For 1 n 2g, we have det(e + 0 ±E n INVERSE PROBLEMS FOR CANONICAL SYSTEMS 25 J) = C8g 2n g det([e 0 + +E 0 J(8g) n ] տn )det([e 0 + E 0 J(8g) n ] տn ) (415) Proof By the proof of Lemma 49, det(e + 0 ±E n Here J) = C4g 2n g det([e + 0 ] տ2g+n +[e 2,n ] ւ2g+n)det([e + 0 ] տ2g+n [e 2,n ] ւ2g+n) [e + 0 ] տ2g+n ±[e 2,n ] ւ2g+n = [E + 0 ±E 0 J(8g) n ] տ2g+n, since [e + 0 ] տ2g+n = [E + 0 ] տ2g+n and [e 2,n ] ւ2g+n = [E 0 J(8g) n ] տ2g+n by definition of matrices Applying Lemma 33 to the 2g columns of [E + 0 ± E 0 J(8g) n ] տ2g+n with indices (n+1,n+2,,2g +n), we have det[e + 0 ±E 0 J(8g) n ] տ2g+n = C 2g g det([e+ 0 ±E 0 J(8g) n ] տn ), since all entries above [E 0 + ±E 0 J(8g) n ] ց2g in [E 0 + ±E 0 J(8g) n ] տ2g+n are zero and [E 0 + ± E0 J(8g) n ] ց2g is lower triangular with diagonal (C g,,c g ) Therefore, we obtain (415) Lemma 413 We have det(e + 0 ±E 0 J(8g) n ) ( 1 (u ± n ()+v± n (0))) = det(e + 0 ±E 0 J(8g) n 2 ) Proof Let w = t[ C g C g 0 0 ] be the column vector of length 8g Denote by u ± n,j the columns of E+ 0 ±E nj By Cramer s formula, det(e + 0 ±E nj)(1 (u ± n()+v ± n(0))) = det(e + 0 ±E nj) det[u ± n,1 n w u ± n,8g ] det[u ± n,1 6g w u ± n,8g ] = det(e + 0 ±E n J) det[u± n,1 n w u ± n,8g ] det[u ± n,1 6g ±w u ± n,8g ] = det(e + 0 ±E nj) det[u ± n,1 n w u ± n,8g ] det[u ± n,1 6g u ± n,6g u± n,8g ]+det[u± n,1 6g u ± n,6g w u± n,8g ] = det[u ± n,1 n w u ± n,8g ]+det[u ± n,1 6g u ± n,6g w u± n,8g ] n 6g 6g n w = det[u ± n,1 w ±w u ± n,8g ] det[u± n,1 u ± n,6g w u± n,8g ] +det[u ± n,1 6g u ± n,6g w u± n,8g ] n 6g = det[u ± n,1 u ± n,6g w Therefore, it is sufficient to prove det(e + 0 ±E 0 J(8g) n 2 )det(e+ 0 ±E n J) = det(e + 0 ±E 0 J(8g) n )det[u ± n,1 u ± n,6g w u± n,8g ] u ± n,6g n w 6g Applying Lemma 33 to the columns (n+1,n+2,,8g), (416) u ± n,6g w u± n,8g ] det(e + 0 ±E 0 J(8g) n 2 ) = det([e+ 0 ±E 0 J(8g) n 2 ] տn 2)det([E + 0 ±E 0 J(8g) n 2 ] ց8g n+2) = C 8g n+2 g det([e 0 + ±E 0 J(8g) n 2 ] տn 2),
26 26 M SUZUKI since all entries above [E 0 + ±E 0 J(8g) n 2 ] ց8g n+2 are zero and [E 0 + ±E 0 J(8g) n 2 ] ց8g n+2 is lower triangular with diagonal (C g,,c g ) By the proof of Lemma 49, det(e + 0 ±E n J) = C4g n g det([e 0 + ] տ2g+n +[e 1,n ] ր2g+n)det([e 0 + ] տ2g+n [e 1,n ] ր2g+n), since [e 1,n ] ր2g+n = [e 2,n ] ւ2g+n We observe that By Lemma 412, [E + 0 ] տ2g+n ±[e 1,n ] ր2g+n = [E + 0 ±E 0 J(8g) n ] տ2g+n (LHS of (416)) = C 16g 2n+2 g det([e 0 + ±E 0 J(8g) n 2 ] տn 2) On the right hand side of (416), we have det([e + 0 +E 0 J(8g) n ] տn )det([e + 0 E 0 J(8g) n ] տn ) det(e 0 + ±E 0 J(8g) n ) = C 8g n g det([e + 0 ±E 0 J(8g) n ] տn ) by applying Lemma 33 to the columns (n+1,n+2,,8g) of E + 0 ±E 0 J(8g) n Therefore, it is sufficient to prove D ± n := det[u ± n,1 u ± n,6g n w 6g u ± n,6g w u± n,8g ] = C 8g n+2 g det([e 0 + E 0 J(8g) n ] տn )det([e 0 + ±E 0 J(8g) n 2 ] տn 2) Subtracting the (6g j)th column from the (n j)th column for 1 j, D ± n = det[v± n,1 v± n, u ± n,6g n w 6g u ± n,6g w u± n,8g ], (417) where v ± n,n j = u± n,n j u n,6g j for 1 j Successively, adding the (n j+1)th column to the (6g j 1)th column for 1 j n 2, D ± n = det[v ± n,1 v± n, n u ± n,6g w v n,6g n+1 v n,6g 2 6g u ± n,6g w u± n,8g ], where v n,6g j 1 ± = u± n,6g j 1 v± n,n j+1 for 1 j n 2 Note that the columns (n+1,,6g n) are not changed Put M ± n := [ v ± n,1 v± n, n u ± n,6g w v n,6g n+1 v n,6g 2 By the deformation of the matrix, we observe that [M ± n ] տ4g = [E + 0 E 0 J(8g) n ] տ4g Therefore, by applying Lemma 33 to the columns (n+1,n+2,,4g), D n ± = C4g n g det([e 0 + E 0 J(8g) n ] տn )det([m n ± ] ց4g) 6g u ± n,6g w u± n,8g] Here we note that [M n ± ] ց2g+n = [E 0 + E 0 J(8g) n 2 ] տ2g+n and that the columns (4g+1,,6g n) are not changed by the deformation of matrix Then, by applying Lemma 33 to the rows (1,,2g n) of [M n ± ] ց2g+n, D n ± = C 6g 2n g det([e 0 + E 0 J(8g) n ] տn )det([e 0 + E 0 J(8g) n 2 ] տ2g+n) Now (417) is proved = C 8g 2n+2 g det([e 0 + E 0 J(8g) n ] տn )det([e 0 + E 0 J(8g) n 2 ] տn 2)
27 Lemma 414 We have for every 1 n 2g Proof If 1 n 2g, we have Therefore, INVERSE PROBLEMS FOR CANONICAL SYSTEMS 27 det(e 0 + ±E 0 J(8g) n ) = C 4g g det(e+ 0 ±e 0 J(4g) n ) = C 6g 1 g E 0 + ±E 0 J(8g) n = det(e + ±E J (2g+1) n ) e + 0 ±e 0 J(4g) n 0 0 e + 0 det(e 0 + ±E 0 J(8g) n ) = C 4g g det(e+ 0 ±e 0 J(4g) n ), since det(e + 0 ) = C4g g Further, if 1 n 2g, we have e + 0 ±e 0 J(4g) n = E + ±E J (2g+1) n 0 [e + 0 ] ց(2g 1) Hence det(e 0 + ±E 0 J(8g) n ) = C 6g 1 g det(e + ±E J n (2g+1) ), and we complete the proof Lemma 415 We have det(e 0 + ±E 1,1 J) = det(e+ 0 ) and det(e 0 + ±E n,1 J) = det(e+ 0 ±E J) for 2 n 2g Proof By the proof of Lemma 49, for det(e 0 + ±E J) = C4g 2n+2 g det(a +B )det(a B ) (418) A = [e + 0 ] տ2g+, B = [e 2, ] ւ2g+ Ontheotherhand,applyingLemma33tothecolumns(2g+n+1,,4g) ofe 0 + ±E n,1 J, and then applying Lemma 33 to the rows (2g +n+1,4g) of the resulting matrix, [ det(e 0 + An ±B ±E n,1j) = C4g 2n g det ], ±B n A n [ where B = [e 2, ] An ±B ւ2g+n Applying Lemma 33 to the first row of ] ±B n A n and then applying Lemma 33 to the (4g +2)th column of the resulting matrix, An ±B det[ ] = C ±B n A g 2 det A ±B n ±B A Hence = C 2 g det(a +B )det(a B ) det(e 0 + ±E n,1j) = C4g 2n+2 g det(a +B )det(a B ) (419) We complete the proof by comparing (418) and (419)
Linear Systems and Matrices
Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) a 11 x 1 + a 12 x 2 + + a 1n x n = b 1 a 21 x 1 + a 22 x 2 + + a 2n x n = b 2.......
More informationDeterminants Chapter 3 of Lay
Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
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 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 informationMATH2210 Notebook 2 Spring 2018
MATH2210 Notebook 2 Spring 2018 prepared by Professor Jenny Baglivo c Copyright 2009 2018 by Jenny A. Baglivo. All Rights Reserved. 2 MATH2210 Notebook 2 3 2.1 Matrices and Their Operations................................
More informationa 11 a 12 a 11 a 12 a 13 a 21 a 22 a 23 . a 31 a 32 a 33 a 12 a 21 a 23 a 31 a = = = = 12
24 8 Matrices Determinant of 2 2 matrix Given a 2 2 matrix [ ] a a A = 2 a 2 a 22 the real number a a 22 a 2 a 2 is determinant and denoted by det(a) = a a 2 a 2 a 22 Example 8 Find determinant of 2 2
More informationDeterminants. Beifang Chen
Determinants Beifang Chen 1 Motivation Determinant is a function that each square real matrix A is assigned a real number, denoted det A, satisfying certain properties If A is a 3 3 matrix, writing A [u,
More informationMATH 2030: EIGENVALUES AND EIGENVECTORS
MATH 2030: EIGENVALUES AND EIGENVECTORS Determinants Although we are introducing determinants in the context of matrices, the theory of determinants predates matrices by at least two hundred years Their
More informationChapter 2:Determinants. Section 2.1: Determinants by cofactor expansion
Chapter 2:Determinants Section 2.1: Determinants by cofactor expansion [ ] a b Recall: The 2 2 matrix is invertible if ad bc 0. The c d ([ ]) a b function f = ad bc is called the determinant and it associates
More informationChapter SSM: Linear Algebra Section Fails to be invertible; since det = 6 6 = Invertible; since det = = 2.
SSM: Linear Algebra Section 61 61 Chapter 6 1 2 1 Fails to be invertible; since det = 6 6 = 0 3 6 3 5 3 Invertible; since det = 33 35 = 2 7 11 5 Invertible; since det 2 5 7 0 11 7 = 2 11 5 + 0 + 0 0 0
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationLinear Algebra M1  FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8.
Linear Algebra M1  FIB Contents: 5 Matrices, systems of linear equations and determinants 6 Vector space 7 Linear maps 8 Diagonalization Anna de Mier Montserrat Maureso Dept Matemàtica Aplicada II Translation:
More informationEquality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.
Introduction Matrix Operations Matrix: An m n matrix A is an mbyn array of scalars from a field (for example real numbers) of the form a a a n a a a n A a m a m a mn The order (or size) of A is m n (read
More informationLinear Algebra: Lecture notes from Kolman and Hill 9th edition.
Linear Algebra: Lecture notes from Kolman and Hill 9th edition Taylan Şengül March 20, 2019 Please let me know of any mistakes in these notes Contents Week 1 1 11 Systems of Linear Equations 1 12 Matrices
More information1 Matrices and Systems of Linear Equations. a 1n a 2n
March 31, 2013 161 16. Systems of Linear Equations 1 Matrices and Systems of Linear Equations An m n matrix is an array A = (a ij ) of the form a 11 a 21 a m1 a 1n a 2n... a mn where each a ij is a real
More information1 Matrices and Systems of Linear Equations
March 3, 203 66. Systems of Linear Equations Matrices and Systems of Linear Equations An m n matrix is an array A = a ij of the form a a n a 2 a 2n... a m a mn where each a ij is a real or complex number.
More informationELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS. 1. Linear Equations and Matrices
ELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS KOLMAN & HILL NOTES BY OTTO MUTZBAUER 11 Systems of Linear Equations 1 Linear Equations and Matrices Numbers in our context are either real numbers or complex
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
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 informationRepresentations of Sp(6,R) and SU(3) carried by homogeneous polynomials
Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Govindan Rangarajan a) Department of Mathematics and Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560 012,
More information0.1 Rational Canonical Forms
We have already seen that it is useful and simpler to study linear systems using matrices. But matrices are themselves cumbersome, as they are stuffed with many entries, and it turns out that it s best
More informationBasic Concepts of Group Theory
Chapter 1 Basic Concepts of Group Theory The theory of groups and vector spaces has many important applications in a number of branches of modern theoretical physics. These include the formal theory of
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationLINEAR ALGEBRA WITH APPLICATIONS
SEVENTH EDITION LINEAR ALGEBRA WITH APPLICATIONS Instructor s Solutions Manual Steven J. Leon PREFACE This solutions manual is designed to accompany the seventh edition of Linear Algebra with Applications
More informationDeterminant lines and determinant line bundles
CHAPTER Determinant lines and determinant line bundles This appendix is an exposition of G. Segal s work sketched in [?] on determinant line bundles over the moduli spaces of Riemann surfaces with parametrized
More informationThe Determinant: a Means to Calculate Volume
The Determinant: a Means to Calculate Volume Bo Peng August 16, 2007 Abstract This paper gives a definition of the determinant and lists many of its wellknown properties Volumes of parallelepipeds are
More informationProblem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that
Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume
More information5.3 Determinants and Cramer s Rule
304 53 Determinants and Cramer s Rule Unique Solution of a 2 2 System The 2 2 system (1) ax + by = e, cx + dy = f, has a unique solution provided = ad bc is nonzero, in which case the solution is given
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
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 information1. General Vector Spaces
1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule
More informationLecture 7. Econ August 18
Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar
More informationFinitedimensional spaces. C n is the space of ntuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a prehilbert space, or a unitary space) if there is a mapping (, )
More informationChapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations
Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2
More informationVector Spaces. Chapter 1
Chapter 1 Vector Spaces Linear algebra is the study of linear maps on finitedimensional vector spaces. Eventually we will learn what all these terms mean. In this chapter we will define vector spaces
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 informationNOTES FOR LINEAR ALGEBRA 133
NOTES FOR LINEAR ALGEBRA 33 William J Anderson McGill University These are not official notes for Math 33 identical to the notes projected in class They are intended for Anderson s section 4, and are 2
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 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 informationDeterminants  Uniqueness and Properties
Determinants  Uniqueness and Properties 222008 In order to show that there s only one determinant function on M(n, R), I m going to derive another formula for the determinant It involves permutations
More informationMath 215B: Solutions 1
Math 15B: Solutions 1 Due Thursday, January 18, 018 (1) Let π : X X be a covering space. Let Φ be a smooth structure on X. Prove that there is a smooth structure Φ on X so that π : ( X, Φ) (X, Φ) is an
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationSUMS OF ENTIRE FUNCTIONS HAVING ONLY REAL ZEROS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 00029939(XX)00000 SUMS OF ENTIRE FUNCTIONS HAVING ONLY REAL ZEROS STEVEN R. ADAMS AND DAVID A. CARDON (Communicated
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationUpper triangular matrices and Billiard Arrays
Linear Algebra and its Applications 493 (2016) 508 536 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Upper triangular matrices and Billiard Arrays
More informationNOTES ON LINEAR ALGEBRA. 1. Determinants
NOTES ON LINEAR ALGEBRA 1 Determinants In this section, we study determinants of matrices We study their properties, methods of computation and some applications We are familiar with the following formulas
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationChapter 4  MATRIX ALGEBRA. ... a 2j... a 2n. a i1 a i2... a ij... a in
Chapter 4  MATRIX ALGEBRA 4.1. Matrix Operations A a 11 a 12... a 1j... a 1n a 21. a 22.... a 2j... a 2n. a i1 a i2... a ij... a in... a m1 a m2... a mj... a mn The entry in the ith row and the jth column
More informationThe Riemann Hypothesis Project summary
The Riemann Hypothesis Project summary The spectral theory of the vibrating string is applied to a proof of the Riemann hypothesis for the Hecke zeta functions in the theory of modular forms. A proof of
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationLECTURE13 : GENERALIZED CAUCHY S THEOREM
LECTURE3 : GENERALIZED CAUCHY S THEOREM VED V. DATAR The aim of this lecture to prove a general form of Cauchy s theorem applicable to multiply connected domains. We end with computations of some real
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 information0 Sets and Induction. Sets
0 Sets and Induction Sets A set is an unordered collection of objects, called elements or members of the set. A set is said to contain its elements. We write a A to denote that a is an element of the set
More informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
More informationIntroduction to Groups
Introduction to Groups S F Ellermeyer November 2, 2006 A group, G, is a set, A, endowed with a single binary operation,, such that: The operation is associative, meaning that a (b c) = (a b) c for all
More informationLecture Notes in Linear Algebra
Lecture Notes in Linear Algebra Dr. Abdullah AlAzemi Mathematics Department Kuwait University February 4, 2017 Contents 1 Linear Equations and Matrices 1 1.2 Matrices............................................
More informationMAT Linear Algebra Collection of sample exams
MAT 342  Linear Algebra Collection of sample exams Ax. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationRemark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationProduct distance matrix of a tree with matrix weights
Product distance matrix of a tree with matrix weights R B Bapat StatMath Unit Indian Statistical Institute, Delhi 7SJSS Marg, New Delhi 110 016, India email: rbb@isidacin Sivaramakrishnan Sivasubramanian
More informationChapter 5 Eigenvalues and Eigenvectors
Chapter 5 Eigenvalues and Eigenvectors Outline 5.1 Eigenvalues and Eigenvectors 5.2 Diagonalization 5.3 Complex Vector Spaces 2 5.1 Eigenvalues and Eigenvectors Eigenvalue and Eigenvector If A is a n n
More information12 CHAPTER 1. PRELIMINARIES Lemma 1.3 (CauchySchwarz inequality) Let (; ) be an inner product in < n. Then for all x; y 2 < n we have j(x; y)j (x; x)
1.4. INNER PRODUCTS,VECTOR NORMS, AND MATRIX NORMS 11 The estimate ^ is unbiased, but E(^ 2 ) = n?1 n 2 and is thus biased. An unbiased estimate is ^ 2 = 1 (x i? ^) 2 : n? 1 In x?? we show that the linear
More informationWeek 910: Recurrence Relations and Generating Functions
Week 910: Recurrence Relations and Generating Functions April 3, 2017 1 Some number sequences An infinite sequence (or just a sequence for short is an ordered array a 0, a 1, a 2,..., a n,... of countably
More informatione j = Ad(f i ) 1 2a ij/a ii
A characterization of generalized KacMoody algebras. J. Algebra 174, 10731079 (1995). Richard E. Borcherds, D.P.M.M.S., 16 Mill Lane, Cambridge CB2 1SB, England. Generalized KacMoody algebras can be
More informationEXERCISES IN MODULAR FORMS I (MATH 726) (2) Prove that a lattice L is integral if and only if its Gram matrix has integer coefficients.
EXERCISES IN MODULAR FORMS I (MATH 726) EYAL GOREN, MCGILL UNIVERSITY, FALL 2007 (1) We define a (full) lattice L in R n to be a discrete subgroup of R n that contains a basis for R n. Prove that L is
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
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 informationELEMENTARY SUBALGEBRAS OF RESTRICTED LIE ALGEBRAS
ELEMENTARY SUBALGEBRAS OF RESTRICTED LIE ALGEBRAS J. WARNER SUMMARY OF A PAPER BY J. CARLSON, E. FRIEDLANDER, AND J. PEVTSOVA, AND FURTHER OBSERVATIONS 1. The Nullcone and Restricted Nullcone We will need
More informationMATH 240 Spring, Chapter 1: Linear Equations and Matrices
MATH 240 Spring, 2006 Chapter Summaries for Kolman / Hill, Elementary Linear Algebra, 8th Ed. Sections 1.1 1.6, 2.1 2.2, 3.2 3.8, 4.3 4.5, 5.1 5.3, 5.5, 6.1 6.5, 7.1 7.2, 7.4 DEFINITIONS Chapter 1: Linear
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
More informationQ N id β. 2. Let I and J be ideals in a commutative ring A. Give a simple description of
Additional Problems 1. Let A be a commutative ring and let 0 M α N β P 0 be a short exact sequence of Amodules. Let Q be an Amodule. i) Show that the naturally induced sequence is exact, but that 0 Hom(P,
More informationTheorem A.1. If A is any nonzero m x n matrix, then A is equivalent to a partitioned matrix of the form. k k nk. mk k mk nk
I. REVIEW OF LINEAR ALGEBRA A. Equivalence Definition A1. If A and B are two m x n matrices, then A is equivalent to B if we can obtain B from A by a finite sequence of elementary row or elementary column
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationLINEAR ALGEBRA REVIEW
LINEAR ALGEBRA REVIEW JC Stuff you should know for the exam. 1. Basics on vector spaces (1) F n is the set of all ntuples (a 1,... a n ) with a i F. It forms a VS with the operations of + and scalar multiplication
More informationRecursion 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 informationII. Determinant Functions
Supplemental Materials for EE203001 Students II Determinant Functions ChungChin Lu Department of Electrical Engineering National Tsing Hua University May 22, 2003 1 Three Axioms for a Determinant Function
More informationQuestion: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of AcI?
Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of AcI? Property The eigenvalues
More informationIntroduction to Linear Algebra. Tyrone L. Vincent
Introduction to Linear Algebra Tyrone L. Vincent Engineering Division, Colorado School of Mines, Golden, CO Email address: tvincent@mines.edu URL: http://egweb.mines.edu/~tvincent Contents Chapter. Revew
More informationDETERMINANTS. , x 2 = a 11b 2 a 21 b 1
DETERMINANTS 1 Solving linear equations The simplest type of equations are linear The equation (1) ax = b is a linear equation, in the sense that the function f(x) = ax is linear 1 and it is equated to
More informationimplies that if we fix a basis v of V and let M and M be the associated invertible symmetric matrices computing, and, then M = (L L)M and the
Math 395. Geometric approach to signature For the amusement of the reader who knows a tiny bit about groups (enough to know the meaning of a transitive group action on a set), we now provide an alternative
More informationLecture Summaries for Linear Algebra M51A
These lecture summaries may also be viewed online by clicking the L icon at the top right of any lecture screen. Lecture Summaries for Linear Algebra M51A refers to the section in the textbook. Lecture
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationand let s calculate the image of some vectors under the transformation T.
Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationLinear Algebra Primer
Introduction Linear Algebra Primer Daniel S. Stutts, Ph.D. Original Edition: 2/99 Current Edition: 4//4 This primer was written to provide a brief overview of the main concepts and methods in elementary
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More information2 TwoPoint Boundary Value Problems
2 TwoPoint Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationContinued fractions for complex numbers and values of binary quadratic forms
arxiv:110.3754v1 [math.nt] 18 Feb 011 Continued fractions for complex numbers and values of binary quadratic forms S.G. Dani and Arnaldo Nogueira February 1, 011 Abstract We describe various properties
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationA matrix A is invertible i det(a) 6= 0.
Chapter 4 Determinants 4.1 Definition Using Expansion by Minors Every square matrix A has a number associated to it and called its determinant, denotedbydet(a). One of the most important properties of
More informationVector Spaces ปร ภ ม เวกเตอร
Vector Spaces ปร ภ ม เวกเตอร 1 5.1 Real Vector Spaces ปร ภ ม เวกเตอร ของจ านวนจร ง Vector Space Axioms (1/2) Let V be an arbitrary nonempty set of objects on which two operations are defined, addition
More informationHASSEMINKOWSKI THEOREM
HASSEMINKOWSKI THEOREM KIM, SUNGJIN 1. Introduction In rough terms, a localglobal principle is a statement that asserts that a certain property is true globally if and only if it is true everywhere locally.
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS
More informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationREPRESENTATION THEORY OF S n
REPRESENTATION THEORY OF S n EVAN JENKINS Abstract. These are notes from three lectures given in MATH 26700, Introduction to Representation Theory of Finite Groups, at the University of Chicago in November
More information