One Method for Construction of Inverse Orthogonal Matrices
Total Page:16
File Type:pdf, Size:1020Kb
ONE METHOD FOR CONSTRUCTION OF INVERSE ORTHOGONAL MATRICES, P. DIÞÃ National Institute of Physics and Nuclear Engineering, P.O. Box MG6, Bucharest, Romania E-mail: [email protected] Received September 16, 2008 An analytical procedure to obtain parametrisation of complex inverse orthogonal matrices starting from complex inverse orthogonal conference matrices is provided. These matrices, which depend on nonzero complex parameters, generalize the complex Hadamard matrices and have applications in spin models and digital signal processing. When the free complex parameters take values on the unit circle they transform into complex Hadamard matrices. 1. INTRODUCTION In a seminal paper [1] Sylvester defined a general class of orthogonal matrices named inverse orthogonal matrices and provided the first examples of what are nowadays called (complex) Hadamard matrices. A particular class of inverse orthogonal matrices Aa()ij are those matrices whose inverse is given 1 t by Aa(1ij ) (1 a ji ) , where t means transpose, and their entries aij satisfy the relation 1 AA nIn (1) where In is the n-dimensional unit matrix. When the entries aij take values on the unit circle A–1 coincides with the Hermitian conjugate A* of A, and in this case (1) is the definition of complex Hadamard matrices. Complex Hadamard matrices have applications in quantum information theory, several branches of combinatorics, digital signal processing, etc. Complex orthogonal matrices unexpectedly appeared in the description of topological invariants of knots and links, see e.g. Ref. [2]. They were called two- weight spin models being related to symmetric statistical Potts models. These matrices have been generalized to two-weight spin models, also called , Paper presented at the National Conference of Physics, September 10–13, 2008, Bucharest-Mãgurele, Romania. Rom. Journ. Phys., Vol. 54, Nos. 5–6, P. 433–440, Bucharest, 2009 434 P. Diþã 2 generalized spin models, by removing the symmetry condition [3]. After that these models have been extended to four-weight spin models, see [4]. These matrices are also known as type II matrices, W, with non-zero entries which satisfy the relations n W B ik nij…n 12 (2) W ij k1 jk see Ref. [5]. It is easily seen that the relations (1) and (2) are equivalent. In the case of spin models the above relation is supplemented with two other relations which take into account the star-triangle relations, which impose further constraints on pairs of W matrices. Particular cases of type II matrices also appeared as a generalized Hadamard transform for processing multiphase or multilevel signals, see [6] and [7], which includes the classical Fourier, Walsh- Hadamard and Reverse Jacket transforms. They are defined for 2n 2n-dimen- sional matrices and depend on a pth root of unity and/or one complex non-zero parameter. The aim of the paper is to provide an analytic method for the construction of inverse orthogonal matrices. It is well known, see [5], that a complete solution was given only for dimensions, n " 5, and here we give a few examples for dimensions n 8. 2. INVERSE ORTHOGONAL MATRICES Our method for construction of inverse orthogonal matrices uses as starting point an analog of the complex n n conference matrices Cn, which are matrices with aii = 0, i = 1, …, n and aij 1 ij that satisfy $ CCnn(1) n I n (3) $ where Cn is the Hermitian conjugate of Cn. Conference matrices Cn are important because by construction the matrix CSCICI $ H DTnnnn (4) 2n DT $ EUCICInn nn is a complex Hadamard matrix, as one can easily verify. To extend the above construction to complex orthogonal conference matrices we have to prescribe a recipe for a proper treating of the zero entries on the main diagonal. Thus the complex inverse orthogonal conference matrices are defined by a similar relation to the relation (3), with aij, ij, complex non-zero numbers, and aii = 0. Our formula for the inverse is given by 3 Method for construction of inverse orthogonal matrices 435 1 t CCII(1 (nn ) n ) (5) where 1/ A is the matrix whose elements are 1aij . We remark that the above relation makes sense since all the entries of Cn + In are complex non-null numbers. In fact one can generalize the above inverse for an arbitrary complex matrix, M, with arbitrary located zeros, the entries of the inverse matrix M–1 being given by I1for0MM 1 J ji ij Mij (6) K 0for0Mij In the case of complex orthogonal matrices the formula (4) takes the form CSCICI 1 O DTnnnn (7) 2n DT 1 EUCICInn nn To see how the method works we start with the simplest example, i.e. we consider the case n = 2, when the orthogonal conference matrix has the form CS0 a C DT (8) n EUa 0 with a+C$ an arbitrary non-zero number. According to (5), or (6), the inverse is given by CS01a C 1 DT (9) n EU10a 1 and one simply verifies that CC22 I 2. By using the relation (4) one gets a complex inverse orthogonal matrix which depends on a single complex parameter, a. CS111aa DT DTaa11 1 D (10) 4 DT 111aa DT EUaa 11 1 Multiplying at right and at left by the diagonal matrices da 11,1 and 1 a d 11,11 generated by the inverse of the first row, respectively the first 2 aa column, one finds, modulo a transposition, the jacket form of D4 as 436 P. Diþã 4 CS11 1 1 DT DT11bb d (11) 4 DT11 bb DT EU1111 where ba1.2 Jacket matrices have entries on the first row and column equal to 1, and on the last row and column 1. When be it with i2 1 and t a real number one finds a matrix which depends on a continuous parameter, t +R, matrix that was found by Hadamard, [8]. Since there are no inverse complex orthogonal conference matrices for n = 3, we consider the case n = 4 and start with CS0111 DT DT10 aa C (12) 4 DT10aa DT EU10 aa where a 0 is an arbitrary complex number. It was obtained from the complex conference matrix given in [9], by changing the phase into the complex number a. To see what is the general procedure for obtaining inverse orthogonal matrices we multiply (12) at right and respectively at left by the diagonal matrices d1 = (a1, a2, a3, a4) and d2 = (1, b1, b2, b3), where ai and bi are non-zero complex numbers to obtain a more general form of matrix C4. By using formula (7) one gets CS11111aaa21134 ababab 1213 DT DTab1111111 aab 31 aab41 a 2 aab 2 2 aab 23 DTab aab11111 aab a aab ab DT1222 42 3 31 33 DTab aab aab11 a 1 aab 1 aab 1 O 1323 33 4 41 42 (13) 8 DT 11111aaa ababab DT21134 1213 DTab111111 aab 31 aab41 a 2 aab22 1 aab23 DT DTab12 aab2211111 aab 42 a 3 aab 31 ab 33 DT EUa13 b aa23 b aa 33 b11 a 4 1 aa 41 b 1 aa 42 b 1 Multiplying O8 at right and, respectively, at left by diagonal matrices generated by the inverse of the first row, respectively, of the first column, one finds after the substitution awaxbaybazb121324311 1 (14) 5 Method for construction of inverse orthogonal matrices 437 the final form of the matrix O8, up to a few permutations necessary to put it in a more symmetric form CS11111111 DT DT11wx aw x a aw wx x a DT11aw y a aw wy wy y a DT DT11 aw z a wz z a aw wz D (15) 8 DT11 aw z a wz z a aw wz DT DT11 wx aw x a aw wx x a DT DT11aw y a aw wy wy y a DT EU11 1 1 1 1 1 1 matrix which depends on five complex arbitrary non-zero parameters. When a, w, x, y, z take values on the 5-dimensional torus T5, i.e. values on the unit circle, D8 gets a complex Hadamard matrix. If in (15) one makes the identification x = y = a one finds, in the terminology used in [7], a jacket matrix which depends on three arbitrary non-zero parameters CS11111111 DT DT111aw aw aw aw 1 DT11aw aw 1 aw aw 1 DT DT11 aw z a wz z a aw wz d (16) 8 DT11 aw z a wz z a aw wz DT DT11 aw aw 1 aw aw 1 DT DT11aw aw 1 aw aw 1 DT EU11111111 The matrix d8 is the generalization of the jacket matrix K4(i) from [7] which depends only on 1 and i. It is easily seen that the above construction can be generalised such that the following proposition holds: Proposition 1. Let us suppose that the n n complex inverse orthogonal conference matrix Cn has the standard form, i.e. all the entries on the first row and the first column are equal to unity, excepting that on the main diagonal which is zero. If such a Cn depends on p complex parameters, then the complex inverse orthogonal matrix D2n obtained by using formula (7) depends on n + p complex parameters. Similarly, by starting with the n = 5 conference matrix 438 P. Diþã 6 CS01111 DT10 babab2 DT2 C5 DT10babab (17) 10abab2 b DT EU10aabb2 where ae 23+i , and b+C$ , one finds by using the above procedure the matrix CS1111111111 DT 22 DT11yy12 y1122 ab aby 1 yy 12 aby y b by 1 y ab DT11by y a2 b a2 by y b y y aby y y y ab DT11131333 1 3 2 DT11aby114141 y y y4144 ab y y by y ab a by y b DT11a2 by y b by y a2 b y y y y y ab aby D DT1115155 155 (18) 10 2 2 DT1 aby1115 y55 b by y ab yy yy 15y5 ab aby1 1 DT 22 DT11y12 y w y 2 ab aby 1 y 12 y a by12 y b by 1 y 2 a b DT11by yy a2 b a2 by y b y y aby y y y ab DT111313331 3 2 2 DT11aby114141 y y y4144 ab y y by y a b a by y b DT EU11 1 1 1 1 1 1 1 1 which depends on six complex parameters: b, y1, y2, y3, y4, y5.