Pascal Eigenspaces and Invariant Sequences of the First or Second Kind
|
|
- Cassandra Booth
- 6 years ago
- Views:
Transcription
1 Pascal Eigenspaces and Invariant Sequences of the First or Second Kind I-Pyo Kim a,, Michael J Tsatsomeros b a Department of Mathematics Education, Daegu University, Gyeongbu, 38453, Republic of Korea b Department of Mathematics and Statistics, Washington State University, Pullman, WA 99164, USA Abstract An infinite real sequence a n } is called an invariant sequence of the first (resp, second ind if a n = n ( n = ( 1 a (resp, a n = ( =n n ( 1 a We review and investigate invariant sequences of the first and second ind, and study their relationships using similarities of Pascal-type matrices and their eigenspaces Keywords: Invariant sequence, Pascal matrix, Eigenvalue, Eigenvector 1 MSC: 15B18, 11B39, 11B65 1 Introduction Inverse relations play an important role in combinatorics [11] The binomial inversion formula, which states that for sequences a n } and b n } (n =, 1,,, a n = n = ( n ( 1 b if and only if b n = n = ( n ( 1 a, (11 is a typical inverse relation considered in [6, 8, 1, 14, 15, 16] Specifically, (11 motivated Sun [14] to investigate the following sequences Definition 11 Let a n } (n =, 1,, be a sequence such that ( 1 s 1 a n = n = ( n ( 1 a ; s = 1 or s = (1 5 We refer to a n } as an invariant sequence (when s = 1 or an inverse invariant sequence (when s = of the first ind Research supported by Daegu University Research Grant 13 Corresponding author addresses: imipyo@daeguacr (I-Pyo Kim, tsat@mathwsuedu (Michael J Tsatsomeros Preprint submitted to Journal of LATEX Templates September 1, 17
2 Several examples of invariant sequences of the first ind can be found in [14], including 1 n }, nf n 1}, L n }, ( 1 n B n } (n, where F 1 = and F n }, L n }, and B n } are the Fibonacci sequence, Lucas sequence, and Bernoulli numbers [7], respectively In this paper, we will establish (see Lemma 1 the modified binomial inversion formula, that is, a n = =n ( ( 1 b if and only if b n = n Motivated by (13, we will introduce and consider the following sequences =n Definition 1 Let a n } (n =, 1,, be a sequence such that ( 1 s 1 a n = =n ( ( 1 a (13 n ( ( 1 a ; s = 1 or s = (14 n We refer to a n } as an invariant sequence (when s = 1 or an inverse invariant sequence (when s = of the second ind 1 15 Naturally arising are the questions of existence, identification, and construction of (inverse invariant sequences of the second ind, as well as the problem of characterizing such sequences and examining their relationship to their counterparts of the first ind Invariant sequences, which are also called self-inverse sequences in [15], have indeed been studied by several authors [6, 8, 14, 15] They are naturally connected to involutory (also nown as involution or self-invertible matrices [9] and to Riordan involutions [5] Involutory matrices find use in numerical methods for differential equations [, 9] They are also useful in cryptography, information theory, and computer security by providing convenient encryption and decryption methods [1] Motivated by Shapiro s open questions [1], Riordan involutions have been intensely investigated as a combinatorial concept [4, 5] In this paper, we investigate invariant sequences by means of the eigenspaces of P D and P T D, where P is the Pascal matrix and D an infinite diagonal matrix with alternating diagonal entries in 1, 1} (see Sections, 3 In fact, P D and P T D are involutory matrices and P D is a Riordan involution Our investigation follows the ideas and connections of invariant sequences to the eigenspaces of P D and P T D developed in Choi et al [6] This will allow us to associate (inverse invariant sequences of the first and second ind, as well as identify and construct such sequences (Section 4 Notation and preliminaries 5 The following notation and conventions are used throughout the manuscript The infinite matrices in this paper have infinite numbers of rows i and columns, with i,, 1,, }
3 E λ (A denotes the eigenspace of a (finite or infinite matrix A corresponding to its eigenvalue λ For a matrix A with columns A ( =, 1,, [ and ] with denoting the vector of zeros in R, A denotes the matrix whose th column is A, where is by convention vacuous For a matrix A, its (possibly infinite row and column index sets are J and K, respectively For J J, K K, let A(J K denote the matrix obtained from A by deleting rows in J and columns in K, and let A[J K ] denote the matrix A(J K, where J = J \J, K = K \K For brevity, write A( K and A(J in place of A( K and A(J, respectively Further, for m, n, 1,, }, we let A m,n = A[, 1,,, m}, 1,,, n}]; A m,m is abbreviated by A m The binomial coefficient ( i choose is denoted by ( i with the convention that it equals when i < or < [ (i P = ] (i, =, 1,, denotes the (infinite Pascal matrix D = diag(1, 1, 1, 1, Infinite real sequences x n } are identified with the infinite dimensional real vector space R consisting of column vectors x = [x, x 1, x, ] T Notice that as a consequence of the binomial inversion formula (11, we have that [ P 1 = DP D = ( 1 i ( i ] (i, =, 1,, Thus (P D 1 = P D and (11 can be converted [6] into a vector equation for x = [a, a 1, a, ] T and y = [b, b 1, b, ] T R, as follows: P Dx = y if and only if P Dy = x (1 Lemma 1 Let P and D be the Pascal matrix and the diagonal matrix defined above, and let x, y R Then P T Dx = y if and only if P T Dy = x, ( and the modified binomial inversion formula (13 holds 45 Proof As P 1 = DP D, we have (P T D 1 = D(P T 1 = P T D As a consequence, ( holds Letting x = [a, a 1, a, ] T and y = [b, b 1, b, ] T implies (13 Let F = [F, F 1, F, ] T and L = [L, L 1, L, ] T denote the vectors in R whose entries are the members of the Fibonacci and Lucas sequences, respectively; that is F =, F 1 = 1, F n = F n 1 + F n (n, L =, L 1 = 1, L n = L n 1 + L n (n 3
4 x The generating functions of F and L are h 1 (x = 1 x x and h (x = x 1 x x, respectively [3, 7] The following fact is nown, however, we include a proof for completeness Lemma P D F = F and P D L = L 5 Proof Let g(x = 1 x 1 x and f(x = 1 x For =, 1,,, the generating function of the th column of P D is g(xf(x [13] Thus the generating function of P DF is [g(x, g(xf(x, g(xf(x, ][F, F 1, F, ] T = g(x(f + F 1 f(x + F f(x + = g(xh 1 (f(x = h 1 (x, which implies that P D F = F The proof of P D L = L is similar That is, 1 and 1 are eigenvalues of P D and consequently of P T D In fact, these are the only eigenvalues of P D and P T D; see [6] The corresponding eigenspaces are infinite dimensional Indeed, if we consider the Pascal-type matrices P and Q constructed via the Pascal matrix P and the matrix [ ] 1 T Q = P +, P then, as shown in [6], the columns of [ T P ] and form bases for E 1 (P D and E 1 (P D, respectively The following observation follows directly from the definitions and properties mentioned above: Observation 3 The entries of x R form an invariant sequence of the first ind if and only if x E 1 (P D; Q 55 an inverse invariant sequence of the first ind if and only if x E 1 (P D; an invariant sequence of the second ind if and only if x E 1 (P T D; an inverse invariant sequence of the second ind if and only if x E 1 (P T D Based on Observation 3, our goal is to study the eigenspaces E λ (P D and E λ (P T D (λ 1, 1} and discover their relationships Our approach entails showing the existence of an infinite invertible matrix N such that N(P T DN 1 = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 (3 6 and D(N 1 T D (P D DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, (4 which are infinite direct sums of copies of P1 T D 1 and P 1 D 1, respectively This result will be applied to characterize E λ (P D and E λ (P T D Extending the wor in [6], we will also show that the columns of P T and Q T ( form bases for E λ (P T D This will indeed allow us to investigate the relationships between invariant sequences of the first and second ind 4
5 3 The Eigenspaces of P T D and P D Let B = [b i ] (i, =, 1,, be the matrix defined by ( 1 i, if i, b i =, if i >, (31 and let J(a denote the infinite Jordan bloc of the form a 1 a 1 O a 1 a 1 O a (3 It readily follows that B 1 = J(1 65 In the next two lemmas, we will construct an infinite matrix N and its inverse M = N 1, which will give rise to similarity transformations of P T D and P D into direct sums as in (3 and (4 Lemma 31 Let m be a positive integer and let N (m = H (m H (m 1 H (1 with H (l = I l B (l = 1,,, m, where B is the matrix defined in (31 Then lim N (m = N = [n i ] m is the infinite matrix defined by n = 1, n = and n i = for i, = 1,,, and n ( 1 i( i 1 + i i =, if 1 i, i 1, if i > 1 Proof Let m be a positive integer and let N (m = H (m H (m 1 H (1 = [n m i ] with nm = 1, n m = and n m i = for i, = 1,, We will prove that n m ( 1 i( + i i =, if i,, if i > by induction on m, where = i 1, if i = 1,,, m; = 1,,, m 1, if i = m + 1, m +, ; = 1,, The claim is clear for m = 1 For m =, by the construction of H (1 and H (, we have n ( 1 i( + i i =, if i,, if i >, 5
6 where = i 1, if i = 1,, 3, 4; = 1,,, 1, if i = 5, 6, ; = 1,,, since for each pair of indices i and with i = 3, 4, and = 1,,, n i = ( l ( l ( 1 l i ( 1 l = ( 1 i l=i l=i ( 1 + i = ( 1 i 1 Let now m 3 Then by the construction of H (m, we have n m i = nm 1 i for all i = 1,,, m and all = 1,, By the induction hypothesis, n m i = ( 1 i( + i, if i,, if i >, where because = n m i = i 1, if i = 1,,, m; = 1,,, m 1, if i = m + 1, m +, ; = 1,,, ( m + l ( 1 l i ( 1 l m l=i ( m + l = ( 1 i m l=i ( m 1 + i = ( 1 i m 1 (33 for each pair of indices i and with i = m 1, m, and = 1,, Thus, by (33, we obtain lim m N (m = N = [n i ] given by n ( 1 i( i 1 + i i =, if i, i 1, if i > for each i, = 1,, Clearly, we have n = 1, n = and n i construction of H (l (l = 1,,, and the proof is complete = for i, = 1,, by the The difference sequence a = [ a, a 1, a, ] T of a sequence a = [a, a 1, a, ] T is defined by a i = a i+1 a i for each i =, 1,, Let a = [ a, a 1, a, ] T ( =, 1,, be the th difference sequence defined inductively by a = ( 1 a, where a = a The infinite matrix a a 1 a a a 1 a a a 1 a 6
7 is called the difference matrix of a It is well nown [3] that for each n =, 1,,, ( ( ( ( n n n n a n = a + a + a + + n a, (34 1 n which is used in the proof of the following lemma Lemma 3 Let M = [m i ] be the matrix with m = 1, m = and m i and ( m i = i, if i, (i, = 1,,, if i > = for i, = 1,,, 7 Then M = N 1, where N is the matrix defined in Lemma 31 Proof Let n i denote the (i, entry of N (i, =, 1,, We would lie to show that n il m l = δ i, l= the Kronecer delta For i = or =, there is nothing to show, so let i 1 and 1 If i = resp i >, then clearly n il m l = l= i ( i 1 n il m li = ( i i 1 = 1 resp l=i n il m l = So it is enough to show that if = i + r with r >, then l= n il m i 1 l = Let = [ and (+r ( consider the sequence z =, +r 1 (, +r (,, ( 1 r 1 +1 ( ] T, ( 1 r,, We can construct the difference matrix A = [a i ] having z as its first column as follows: A = ( +r ( +r 1 ( 1 r ( +r 1 ( +r 1 ( r 1 ( +r ( +r 3 ( 1 r Y ( 1 r ( +1 1 ( 1 r ( ( 1 r 1( 1 ( 1 r 1( 1 ( 1 r( ( 1 r( 1 1 ( 1 r( ( 1 r ( + ( 1 r 1( +1 l=, 7
8 where Y is the (r + 1 (r + 1 matrix given by Y = Thus a = for = + 1,, + r Since + 1 i+r i 1+r+r = + r, by using (34 we obtain where = i 1 n il m l = l= = i+r ( ( + l i n il m l = ( 1 l i t i + r l l=i t ( ( + r s t ( 1 r+s = ( 1 r a t =, s l=i s= i+r and t =, completing the proof Let M (m = U (1 U ( U (m, where U (l = I l J(1 (l = 1,,, m where J(1 is the matrix defined in (3 Then, by Lemmas 31 and 3, we have lim m M (m = M, which is the matrix defined in Lemma 3 In fact we have N =
9 and M = N 1 = We can now state and prove the similarity transformations of P D and P T D claimed in (3 and (4 [ (i Theorem 33 Let P = ] (i, =, 1,, and D = diag(1, 1, 1, 1, Then, 75 (a N(P T DM = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1, (b (DM T D (P D (DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, where P 1 resp D 1 are the leading submatrices of P resp D, and N resp M are the matrices defined in Lemmas 31 resp 3 Proof (a For each = 1,,, let H ( = I B and U( = I J(1, where B is the matrix defined in (31 Let m, n be arbitrary positive integers with n = m + 1 First, we will show by induction on m that ( H(m H (m 1 H (1 P T D ( U (1 U ( U (m = Z (P T D, where Z is an n n matrix such that Z = m }} (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 When m = 1, the first row of H (1 P T DU (1 is clearly [1, 1,,, ] Since for a pair of indices i and with i, = 1,,, (( ( ( ( 1 (H (1 P T D i = ( ( 1 = ( 1, i i + 1 i 1 the second row of H (1 P T DU (1 is [, 1, 1, 1, 1, ]U (1 = [, 1,,,, ] For a pair of indices i and with i,, we have ( ( ( 1 (H (1 P T DU (1 i = ( ( 1 = ( 1, i 1 i 1 i and so H (1 P T DU (1 = (P T 1 D 1 (P T D By induction on m, it follows N (m P T DU (m = m }} (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 (P T D, 9
10 where N (m = H (m H (m 1 H (1 and M (m = U (1 U ( U (m Since m was an arbitrary positive integer, NP T DM = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1, where N = lim m N (m and M = lim m M (m (b It follows directly from (a that (DM T DP D(DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, completing the proof 8 Let e i denote the ith column of the identity matrix I (i =, 1, Then E λ (P T D and E λ (P D (λ 1, 1} can be characterized via Theorem 33 as follows: Theorem 34 Let N resp M = N 1 be the matrices in Lemma 31 resp Lemma 3 Then the following hold: (a Me, Me, Me 4, } is a basis for E 1 (P T D 85 (b M(e + e 1, M(e + e 3, } is a basis for E 1 (P T D (c DN T D(e + e 1, DN T D(e + e 3 } is a basis for E 1 (P D (d DN T De 1, DN T De 3, DN T De 5, } is a basis for E 1 (P D 9 Proof (a Let B (1 = e, e, e 4, } and x = [x, x 1, x, ] T E 1 (NP T DM Then, by Theorem 33, we have ( [ ] (NP T 1 DM Ix = x =, i=1 which implies that x i = t i, x i+1 = for each i =,, 4, and t i R So B (1 spans E 1 (NP T DM and since e, e, e 4, are linearly independent, B (1 is a basis for E 1 (NP T DM Therefore Me, Me, Me 4, } is a basis for E 1 (P T D (b Let B ( 1 = e + e 1, e + e 3, } and y = [y, y 1, y, ] T E 1 (NP T DM Then, by Theorem 33, we have ( [ ] (NP T 1 DM + Iy = y =, i=1 [ ] [ ] y i 1 which implies that = s i for each i =,, 4, and s i R So B ( 1 spans y i+1 E 1 (NP T DM and since e + e 1, e + e 3, are linearly independent, B ( 1 is a basis for E 1 (NP T DM Therefore M(e + e 1, M(e + e 3, } is a basis for E 1 (P T D Clauses (c and (d can be proven similarly 1
11 Consider now P T = The following is a matrix expression of Theorem 34; the last two clauses appeared in [6] [ ] 1 T [ (i Corollary 35 Let Q = P + and D = diag(( 1, ( 1 1,, where P = P ] (i, =, 1,, Then the following hold: 95 (a The columns of P T form a basis for E 1 (P T D (b The columns of Q T ( form a basis for E 1 (P T D [ ] (c The columns of T P form a basis for E 1 (P D (d The columns of Q form a basis for E 1 (P D Proof (a Let (i, be a pair of indices i and with i, The ith component of Me equals ( }} i when i, and equals, otherwise; Me = [,,, ( (, ( 1,,,,, ] T is the th column of P T, which implies (a (b The ith component of M(e +e +1 equals ( ( ( i + i 1 = +1 ( i + i 1, when i +1, and equals otherwise; }} ( ( ( M(e + e +1 = [,,,,,,,,, ] T }} ( ( ( + [,,,,,,,,, ] T 1 [ ] T is the th column of P T 1 T ( + ( = Q T (, which implies (b Using the fact P T that for a pair of indices i and with i 1, ( 1 (DN T D i = ( 1 i ( 1 i + i 1 ( 1, clauses (c and (d can be proven similarly 11
12 1 4 Invariant sequences of two inds: Relations and examples We begin by some basic examples of (inverse invariant sequences Example 41 It follows from Corollary 35 (c and (d that the Fibonacci sequence F is an inverse invariant sequence of the first ind and the Lucas sequence L is an invariant sequence of the first ind The ith row of P T (i =, 1, is i }} ( i [,,, ( ( i i 1 i,,, 1, ( i,,, ], (41 15 from which we can get that J(F is an invariant sequence of the second ind and J(L is an inverse invariant sequence of the second ind Recall that J( is the infinite Jordan bloc with in the main diagonal, which is the matrix defined in (3 In the following theorem, we provide a general mechanism for transforming invariant sequences into inverse invariant sequences, and vice versa 11 Theorem 4 Let J(λ denote the matrix defined in (3 Then the following hold: (a If x is an invariant sequence of the second ind, then (J(1 + J( T x is an inverse invariant sequence of the second ind (b If x is an inverse invariant sequence of the second ind, then J( 1 J(x is an invariant sequence of the second ind (c If x is an invariant sequence of the first ind, then ( J(J( 1 T x is an inverse invariant 115 sequence of the first ind (d If x is an inverse invariant sequence of the first ind, then (J( 1 + J(x is an invariant sequence of the first ind Proof (a If x is an invariant sequence of the second ind, then by Corollary 35 (a, there exists b R such that P T b = x Let L = [l i ] (i, =, 1,, be the infinite lower triangular matrix defined by ( 1 i+, if i, l i =, if i < Then LP T ( = P T because for a pair of indices i and with i,, the (i, entry of LP T ( equals ( ( ( ( ( 1 i+ + ( 1 i ( 1 i++i = 1 i i 1
13 1 when i, and equals otherwise, which coincides with the (i, entry of P T Thus, [ ] LQ T ( = P T + L T P T From the fact that L 1 = J(1 T and by Corollary 35 (b, it follows that [ ] Q T ( b = J(1 T P T b + T P T (4 b = J(1 T x + [, x T ] T = (J(1 + J( T x is an inverse invariant sequence of the second ind (b Let x be an inverse invariant sequence of the second ind From (4 and J(J(1 T = J(1, it readily follows that J( 1 J(Q T ( = P T and so by Corollary 35 (b, we conclude that J( 1 J(x is an invariant sequence of the second ind 1 T T (c Since Q = P + T, we obtain J( T Q = J( T P 1 T + P T For nonnegative P integers i and with i, =, 1,,, since ( ( i i + i ( i i ( i+ = i+1 i++1, it implies that ([ ] [ ] [ ] T Q T P T P Ω =, where Ω is the infinite (, 1-matrix with 1 s everywhere on and [ ] [ ] below its main diagonal Since J( 1 T = Ω 1, we get (I J( 1 T 1 T T =, which implies that ( J(J( 1 T 15 x is an inverse invariant sequence of the first ind Clause (d easily follows from (c similarly Q P Let τ 1 = 1+ 5 and τ = 1 5 It is well nown that F n = 1 5 τ n τ n and L n = τ n 1 + τ n where F n and L n are the nth terms of F and L, respectively (n =, 1,, [3] Since J(F is an invariant sequence of the second ind by (41, it follows from Theorem 4 (a that J(1 T J(F + J( T J(F = J(1 T J(F + F is an inverse invariant sequence of the second ind In fact, the nth term (n =, 1,, of J(1 T J(F + F is 1 τ1 n+ 1 τ n+ + 1 τ1 n 1 τ n = L n+1, which implies that J(1 T J(F + F = J(L On the other hand, since for a pair of nonnegative integers i and with i, =, 1,,, (J( 1 ( 1 i ( 1 i = i+1, if i,, if i >, 13
14 and the n-th term of J( 1 J( L is ( 1 (1/ +1 (τ1 ++n + τ ++n = 1 τ1 n+1 1 τ n+1 = F n+1, 5 5 = we have J( 1 J( L = J(F, which is an invariant sequence of the second ind by Theorem 4 (b, as stated in Example 41 It directly follows from Theorem 4 (d that (J(+J( 1F, namely L, is an invariant sequence of the first ind Since for a pair of nonnegative integers i and with i, =, 1,,, ( J( 1 ( 1 i = i+1, if i,, if i >, we obtain ( J( 1 T [, L T ] T = F, which is an inverse invariant sequence of the first ind by Theorem 4 (c, because for n =, 1,,, n 1 n 1 (1/ n L = (1/ n ((τ 1 + (τ = 1 τ n 1 τ n = F n 5 5 = = 13 The sequence B = (B, B 1, B, T defined by B = 1 and n ( n+1 = B = (n 1 comprises the Bernoulli numbers and DB is an invariant sequence of the first ind [14], which also follows directly from the fact that P DDB = DB A new inverse invariant sequence of the first ind from the Bernoulli numbers B is provided next See Table 1 for explicit members of these sequences Corollary 43 Let B = (B, B 1, B, T be the Bernoulli numbers Then the sequence K = (K, K 1, K, T defined by n 1 K =, K n = (1/ n ( 1 B (n = 1,, is an inverse invariant sequence of the first ind Proof It follows from Theorem 4 (c that ( J( 1 T J( T DB = ( J( 1 T [, DB T ] T = 135 is an inverse invariant sequence of the first ind Notice now that the sequence K in the statement is indeed equal to ( J( 1 T [, DB T ] T By Theorem 4 (d, we get DB = (J(+J( 1K, since the first component of (J(+J( 1K is clearly ( 1 B, and for each positive integer i with i = 1,,, ith component of (J(+J( 1K is (1/ i ( 1 B + i 1 = i (1/ i+1 ( 1 B = ( 1 i B i = By Corollary 35 and Theorem 4, we can directly get more (inverse invariant sequences of the first and second ind as follows: 14
15 Table 1: (Inverse invariant sequences of the first ind associated with the Bernoulli numbers n B DB 1 1 K Corollary 44 For a positive integer n, let Φ n = S ln S ln 1 S l3 S l S l1, where P T, if l i = 1, S li = J(1 T + J( T, if l i = 1, and let Φ n = S ln Sln 1 S l3 Sl Sl1 where J( S 1 J(, if l i = 1, li = Q T (, if l i = 1 for i = 1,,, n Then for x R, we have the following: 14 (a If n is odd (even and l i = ( 1 i+1 for i = 1,,, n, then Φ n x is an invariant (inverse invariant sequence of the second ind (b If n is odd (even and l i = ( 1 i for i = 1,,, n, then Φ n x is an inverse invariant (invariant sequence of the second ind Example 45 It follows from (41 and Corollary 44 (a that for x = [x, x 1, x, ] T R, Φ x = y is an inverse invariant sequence of the second ind, where y = [y, y 1, y, ] T with y i = i t= i (( t + i 1 t ( t + 1 x t (i =, 1,, i t For example, let x = [,,,,,,, 1,,, ] T This results to the inverse invariant sequence of the second ind y = [y, y 1, y, ] T, where ( ( 7 8 y i = + (i =, 1,, i 1 7 i 7 That is, y = = y 6 =, y 16 = y 17 = =, and by direct calculation, one can compute the nonzero components y 7, y 8,, y 15 ; eg, y 1 = ( ( = 77 Corollary 46 For a positive integer n, let Ψ n = T ln T ln 1 T l3 T l T l1, where T li = Q, if l i = 1, ( J( 1 T J( T, if l i = 1 15
16 and let Ψ n = T ln Tln 1 T l3 Tl Tl1, where J( + J( 1, if l i = 1, [ ] T li = T, if l i = 1 P 145 for i = 1,,, n Then for x R, we have the following: (a If n is odd (even and l i = ( 1 i+1 for i = 1,,, n, then Ψ n x is an invariant (inverse invariant sequence of the first ind (b If n is odd (even and l i = ( 1 i for i = 1,,, n, then Ψ n x is an inverse invariant (invariant sequence of the first ind [ ] T Example 47 For i = 1,,, the ith row of is [ ( ( i 1, i ( i 1 1,,,,, ] It follows i 1 P that for x = [x, x 1, x, ] T R, Ψ x = y is an invariant sequence of the first ind, where by Corollary 46 (b, y = (y, y 1, y, T satisfies i (( i 1 t y i = + t 1 t= ( i t t x t For example, let x = [,,,,,,, 1,,, ] T This results to the invariant sequence of the first ind y = [y, y 1, y, ] T, where ( ( i 1 7 i 7 y i = + (i = 1,, Thus y = y 1 = = y 13 = and one can find by direct calculation, the nonzero components y 14, y 15, ; eg, y 14 = ( ( = In the next two theorems, we obtain direct relationships among (inverse invariant sequences of the first ind and (inverse invariant sequences of the second ind 155 Theorem 48 Let x and y be, respectively, either (i an invariant sequence of the first ind and an inverse invariant sequence of the second ind, or (ii an inverse invariant sequence of the first ind and an invariant sequence of the second ind Then x T Dy = 16 Proof From P Dx = λx and P T Dy = λy, we have x T DP T Dy = λx T Dy, which implies that λx T Dy = for λ 1, 1} The following theorem is a useful tool for proving the final theorem Theorem 49 For every positive integer n, 16
17 (a the columns of (P +D n are invariant sequences of the first ind, and the columns of (P D n are inverse invariant sequences of the first ind 165 (b the columns of (P T + D n are invariant sequences of the second ind, and the columns of (P T D n are inverse invariant sequences of the second ind Proof Let n be a positive integer Then and P D(P + D n = P D(P + D(P + D n 1 = (P + D n P T D(P T D n = P T D(P T D(P T D n 1 = (P T D n, which respectively imply that each column of (P +D n is an invariant sequence of the first ind, and each column of (P T D n is an inverse invariant sequence of the second ind The other assertions of the theorem follow similarly 17 The following theorem is a form of converse of Theorem 48 Theorem 41 Let x i R (i =, 1,,, y R \ }, and let X = [x, x 1, x, ] Then for each i =, 1,,, 175 (a if x i T Dy = and X = P + D (X = P D, then y is an inverse invariant (invariant sequence of the second ind, and x i is an invariant (inverse invariant sequence of the first ind (b if x i T Dy = and X = P T + D (X = P T D, then y is an inverse invariant (invariant sequence of the first ind and x i is an inverse invariant (invariant sequence of the second ind 18 Proof If x T i Dy = and X = P + D, then since X T Dy = (P T + DDy =, we have P T Dy = y and by Theorem 49, P Dx i = x i Thus y is an inverse invariant sequence of the second ind and x i is an invariant sequence of the first ind This proves the first case of part (a For the case of X = P D and part (b, the results can be shown similarly References 185 [1] B Acharya, SK Panigrahy, SK Patra, G Panda, Image encryption using advanced Hill Ciper Algorithm, Int J Recent Trends in Engineering, 1(1 ( [] AL Andrew, Eigenvectors of certain matrices, Linear Algebra Appl, 7 ( [3] RA Brualdi, Introductory combinatorics, Fourth ed, Prentice-Hall, Upper Saddle River, 4 17
18 19 [4] NT Cameron, A Nwanta, On some(pseudo involutions in the Riordan group, J Integer Seq 8, (5 Article 537 [5] G-S Cheon, H Kim, Simple proofs of open problems about the structure of involutions in the Riordan group, Linear Algebra Appl, 48 ( [6] G-S Choi, S-G Hwang, I-P Kim BL Shader, (±1-Invariant sequences and truncated Fibonacci sequences, Linear Algebra Appl, 395 ( [7] L Comtet, Advanced Combinatorics, Reidel, Dordrecht, 1974 [8] R Donaghey, Binomial self-inverse sequences and tangent coefficient, J Combin Theory Ser A, 1 (1976, [9] RA Horn, CR Johnson, Matrix Analysis Second ed, Cambridge University Press, New Yor, 13 [1] H Prodinger, Some information about the binomial transform, Fibonacci Quarterly, 3 ( [11] J Riordan, Combinatorial Identities, Wiley, New Yor, 1968 [1] LW Shapiro, Some open questions about random wals, involutions, limiting distributions and generating functions, Adv Appl Math, 7 ( [13] LW Shapiro, S Getu, W-J Woan, LC Woodson, The Riordan group, Discrete Appl Math, 34 ( [14] ZH Sun, Invariant sequences under binomial transformation, Fibonacci Quarterly, 39 ( [15] Y Wang, Self-inverse sequences related to a binomial inverse pair, Fibonacci Quarterly, 43 ( [16] K Woci, Binomial transform and Dole sequences, J Integer Seq, 18 (15 Article
Applications of Riordan matrix functions to Bernoulli and Euler polynomials
Illinois Wesleyan University From the SelectedWors of Tian-Xiao He 06 Applications of Riordan matrix functions to Bernoulli and Euler polynomials Tian-Xiao He Available at: https://worsbepresscom/tian_xiao_he/78/
More informationDiscrete Applied Mathematics
Discrete Applied Mathematics 157 (2009 1696 1701 Contents lists available at ScienceDirect Discrete Applied Mathematics journal homepage: www.elsevier.com/locate/dam Riordan group involutions and the -sequence
More informationMAPPING AND PRESERVER PROPERTIES OF THE PRINCIPAL PIVOT TRANSFORM
MAPPING AND PRESERVER PROPERTIES OF THE PRINCIPAL PIVOT TRANSFORM OLGA SLYUSAREVA AND MICHAEL TSATSOMEROS Abstract. The principal pivot transform (PPT) is a transformation of a matrix A tantamount to exchanging
More informationYi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian , China (Submitted June 2002)
SELF-INVERSE SEQUENCES RELATED TO A BINOMIAL INVERSE PAIR Yi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian 116024, China (Submitted June 2002) 1 INTRODUCTION Pairs of
More informationBounds for Levinger s function of nonnegative almost skew-symmetric matrices
Linear Algebra and its Applications 416 006 759 77 www.elsevier.com/locate/laa Bounds for Levinger s function of nonnegative almost skew-symmetric matrices Panayiotis J. Psarrakos a, Michael J. Tsatsomeros
More informationGeometric Mapping Properties of Semipositive Matrices
Geometric Mapping Properties of Semipositive Matrices M. J. Tsatsomeros Mathematics Department Washington State University Pullman, WA 99164 (tsat@wsu.edu) July 14, 2015 Abstract Semipositive matrices
More informationThe generalized order-k Fibonacci Pell sequence by matrix methods
Journal of Computational and Applied Mathematics 09 (007) 33 45 wwwelseviercom/locate/cam The generalized order- Fibonacci Pell sequence by matrix methods Emrah Kilic Mathematics Department, TOBB University
More informationOn identities with multinomial coefficients for Fibonacci-Narayana sequence
Annales Mathematicae et Informaticae 49 08 pp 75 84 doi: 009/ami080900 http://amiuni-eszterhazyhu On identities with multinomial coefficients for Fibonacci-Narayana sequence Taras Goy Vasyl Stefany Precarpathian
More informationarxiv: v1 [math.co] 11 Jul 2015
arxiv:50703055v [mathco] Jul 05 Duals of Bernoulli Numbers and Polynomials and Euler Number and Polynomials Tian-Xiao He and Jinze Zheng Department of Mathematics Illinois Wesleyan University Bloomington,
More informationMath 113 Homework 5. Bowei Liu, Chao Li. Fall 2013
Math 113 Homework 5 Bowei Liu, Chao Li Fall 2013 This homework is due Thursday November 7th at the start of class. Remember to write clearly, and justify your solutions. Please make sure to put your name
More informationCentral Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J
Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J Frank Curtis, John Drew, Chi-Kwong Li, and Daniel Pragel September 25, 2003 Abstract We study central groupoids, central
More informationIntrinsic products and factorizations of matrices
Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 5 3 www.elsevier.com/locate/laa Intrinsic products and factorizations of matrices Miroslav Fiedler Academy of Sciences
More informationThe k-fibonacci matrix and the Pascal matrix
Cent Eur J Math 9(6 0 403-40 DOI: 0478/s533-0-0089-9 Central European Journal of Mathematics The -Fibonacci matrix and the Pascal matrix Research Article Sergio Falcon Department of Mathematics and Institute
More informationOn the adjacency matrix of a block graph
On the adjacency matrix of a block graph R. B. Bapat Stat-Math Unit Indian Statistical Institute, Delhi 7-SJSS Marg, New Delhi 110 016, India. email: rbb@isid.ac.in Souvik Roy Economics and Planning Unit
More informationON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS
Hacettepe Journal of Mathematics and Statistics Volume 8() (009), 65 75 ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS Dursun Tascı Received 09:0 :009 : Accepted 04 :05 :009 Abstract In this paper we
More informationMath 3191 Applied Linear Algebra
Math 9 Applied Linear Algebra Lecture 9: Diagonalization Stephen Billups University of Colorado at Denver Math 9Applied Linear Algebra p./9 Section. Diagonalization The goal here is to develop a useful
More informationDefinition (T -invariant subspace) Example. Example
Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin
More informationRESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices
Linear and Multilinear Algebra Vol. 00, No. 00, Month 200x, 1 15 RESEARCH ARTICLE An extension of the polytope of doubly stochastic matrices Richard A. Brualdi a and Geir Dahl b a Department of Mathematics,
More informationOn sum of powers of the Laplacian and signless Laplacian eigenvalues of graphs
On sum of powers of the Laplacian and signless Laplacian eigenvalues of graphs Saieed Akbari 1,2 Ebrahim Ghorbani 1,2 Jacobus H. Koolen 3,4 Mohammad Reza Oboudi 1,2 1 Department of Mathematical Sciences
More informationThe Perron eigenspace of nonnegative almost skew-symmetric matrices and Levinger s transformation
Linear Algebra and its Applications 360 (003) 43 57 www.elsevier.com/locate/laa The Perron eigenspace of nonnegative almost skew-symmetric matrices and Levinger s transformation Panayiotis J. Psarrakos
More informationSOME RESULTS ON q-analogue OF THE BERNOULLI, EULER AND FIBONACCI MATRICES
SOME RESULTS ON -ANALOGUE OF THE BERNOULLI, EULER AND FIBONACCI MATRICES GERALDINE M. INFANTE, JOSÉ L. RAMÍREZ and ADEM ŞAHİN Communicated by Alexandru Zaharescu In this article, we study -analogues of
More informationEVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS
EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS CHARLES HELOU AND JAMES A SELLERS Abstract Motivated by a recent work about finite sequences where the n-th term is bounded by n, we evaluate some classes
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 informationOn the closures of orbits of fourth order matrix pencils
On the closures of orbits of fourth order matrix pencils Dmitri D. Pervouchine Abstract In this work we state a simple criterion for nilpotentness of a square n n matrix pencil with respect to the action
More informationRepresentation of doubly infinite matrices as non-commutative Laurent series
Spec. Matrices 217; 5:25 257 Research Article Open Access María Ivonne Arenas-Herrera and Luis Verde-Star* Representation of doubly infinite matrices as non-commutative Laurent series https://doi.org/1.1515/spma-217-18
More informationON THE MATRIX EQUATION XA AX = X P
ON THE MATRIX EQUATION XA AX = X P DIETRICH BURDE Abstract We study the matrix equation XA AX = X p in M n (K) for 1 < p < n It is shown that every matrix solution X is nilpotent and that the generalized
More informationSome results on the existence of t-all-or-nothing transforms over arbitrary alphabets
Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets Navid Nasr Esfahani, Ian Goldberg and Douglas R. Stinson David R. Cheriton School of Computer Science University of
More informationLinear transformations preserving the strong q-log-convexity of polynomials
Linear transformations preserving the strong q-log-convexity of polynomials Bao-Xuan Zhu School of Mathematics and Statistics Jiangsu Normal University Xuzhou, PR China bxzhu@jsnueducn Hua Sun College
More informationMatrix functions that preserve the strong Perron- Frobenius property
Electronic Journal of Linear Algebra Volume 30 Volume 30 (2015) Article 18 2015 Matrix functions that preserve the strong Perron- Frobenius property Pietro Paparella University of Washington, pietrop@uw.edu
More informationAnimals and 2-Motzkin Paths
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 8 (2005), Article 0556 Animals and 2-Motzkin Paths Wen-jin Woan 1 Department of Mathematics Howard University Washington, DC 20059 USA wwoan@howardedu
More informationThe Jordan Canonical Form
The Jordan Canonical Form The Jordan canonical form describes the structure of an arbitrary linear transformation on a finite-dimensional vector space over an algebraically closed field. Here we develop
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 informationHomework sheet 4: EIGENVALUES AND EIGENVECTORS. DIAGONALIZATION (with solutions) Year ? Why or why not? 6 9
Bachelor in Statistics and Business Universidad Carlos III de Madrid Mathematical Methods II María Barbero Liñán Homework sheet 4: EIGENVALUES AND EIGENVECTORS DIAGONALIZATION (with solutions) Year - Is
More informationInvestigating Geometric and Exponential Polynomials with Euler-Seidel Matrices
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 14 (2011), Article 11.4.6 Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices Ayhan Dil and Veli Kurt Department of Mathematics
More informationFrame 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 informationDeterminants of Partition Matrices
journal of number theory 56, 283297 (1996) article no. 0018 Determinants of Partition Matrices Georg Martin Reinhart Wellesley College Communicated by A. Hildebrand Received February 14, 1994; revised
More informationTakao Komatsu School of Mathematics and Statistics, Wuhan University
Degenerate Bernoulli polynomials and poly-cauchy polynomials Takao Komatsu School of Mathematics and Statistics, Wuhan University 1 Introduction Carlitz [6, 7] defined the degenerate Bernoulli polynomials
More informationRANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA
Discussiones Mathematicae General Algebra and Applications 23 (2003 ) 125 137 RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Seok-Zun Song and Kyung-Tae Kang Department of Mathematics,
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More information(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true.
1 Which of the following statements is always true? I The null space of an m n matrix is a subspace of R m II If the set B = {v 1,, v n } spans a vector space V and dimv = n, then B is a basis for V III
More informationMaximizing the numerical radii of matrices by permuting their entries
Maximizing the numerical radii of matrices by permuting their entries Wai-Shun Cheung and Chi-Kwong Li Dedicated to Professor Pei Yuan Wu. Abstract Let A be an n n complex matrix such that every row and
More informationA q-analogue OF THE GENERALIZED FACTORIAL NUMBERS
J. Korean Math. Soc. 47 (2010), No. 3, pp. 645 657 DOI 10.4134/JKMS.2010.47.3.645 A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS Seok-Zun Song, Gi-Sang Cheon, Young-Bae Jun, and LeRoy B. Beasley Abstract.
More informationHessenberg Pairs of Linear Transformations
Hessenberg Pairs of Linear Transformations Ali Godjali November 21, 2008 arxiv:0812.0019v1 [math.ra] 28 Nov 2008 Abstract Let K denote a field and V denote a nonzero finite-dimensional vector space over
More informationTwo applications of the theory of primary matrix functions
Linear Algebra and its Applications 361 (2003) 99 106 wwwelseviercom/locate/laa Two applications of the theory of primary matrix functions Roger A Horn, Gregory G Piepmeyer Department of Mathematics, University
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More informationA Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices
University of Wyoming Wyoming Scholars Repository Mathematics Faculty Publications Mathematics 9-1-1997 A Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices Bryan L. Shader University
More informationA LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS
Applicable Analysis and Discrete Mathematics, 1 (27, 371 385. Available electronically at http://pefmath.etf.bg.ac.yu A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS H. W. Gould, Jocelyn
More information2 b 3 b 4. c c 2 c 3 c 4
OHSx XM511 Linear Algebra: Multiple Choice Questions for Chapter 4 a a 2 a 3 a 4 b b 1. What is the determinant of 2 b 3 b 4 c c 2 c 3 c 4? d d 2 d 3 d 4 (a) abcd (b) abcd(a b)(b c)(c d)(d a) (c) abcd(a
More informationImproved Upper Bounds for the Laplacian Spectral Radius of a Graph
Improved Upper Bounds for the Laplacian Spectral Radius of a Graph Tianfei Wang 1 1 School of Mathematics and Information Science Leshan Normal University, Leshan 614004, P.R. China 1 wangtf818@sina.com
More informationCatalan triangle numbers and binomial coefficients
Contemporary Mathematics Volume 73, 208 https://doiorg/0090/conm/73/435 Catalan triangle numbers and binomial coefficients Kyu-Hwan Lee and Se-jin Oh Abstract The binomial coefficients and Catalan triangle
More informationCATALAN TRIANGLE NUMBERS AND BINOMIAL COEFFICIENTS arxiv:6006685v2 [mathco] 7 Oct 207 KYU-HWAN LEE AND SE-JIN OH Abstract The binomial coefficients and Catalan triangle numbers appear as weight multiplicities
More informationKernels of Directed Graph Laplacians. J. S. Caughman and J.J.P. Veerman
Kernels of Directed Graph Laplacians J. S. Caughman and J.J.P. Veerman Department of Mathematics and Statistics Portland State University PO Box 751, Portland, OR 97207. caughman@pdx.edu, veerman@pdx.edu
More informationarxiv: v1 [math.ra] 8 Aug 2016
Each n-by-n matrix with n > 1 is a sum of 5 coninvolutory matrices arxiv:1608.02503v1 math.ra 8 Aug 2016 Ma. Nerissa M. Abara b, Dennis I. Merino 1,, Vyacheslav I. Rabanovich c, Vladimir V. Sergeichuk
More informationSpectrally arbitrary star sign patterns
Linear Algebra and its Applications 400 (2005) 99 119 wwwelseviercom/locate/laa Spectrally arbitrary star sign patterns G MacGillivray, RM Tifenbach, P van den Driessche Department of Mathematics and Statistics,
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 informationAN ALGEBRAIC INTERPRETATION OF THE q-binomial COEFFICIENTS. Michael Braun
International Electronic Journal of Algebra Volume 6 (2009) 23-30 AN ALGEBRAIC INTERPRETATION OF THE -BINOMIAL COEFFICIENTS Michael Braun Received: 28 February 2008; Revised: 2 March 2009 Communicated
More informationA hyperfactorial divisibility Darij Grinberg *long version*
A hyperfactorial divisibility Darij Grinberg *long version* Let us define a function H : N N by n H n k! for every n N Our goal is to prove the following theorem: Theorem 0 MacMahon We have H b + c H c
More informationResearch Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation
Applied Mathematics Volume 20, Article ID 423163, 14 pages doi:101155/20/423163 Research Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation
More informationEigenvectors and Reconstruction
Eigenvectors and Reconstruction Hongyu He Department of Mathematics Louisiana State University, Baton Rouge, USA hongyu@mathlsuedu Submitted: Jul 6, 2006; Accepted: Jun 14, 2007; Published: Jul 5, 2007
More informationLecture 15, 16: Diagonalization
Lecture 15, 16: Diagonalization Motivation: Eigenvalues and Eigenvectors are easy to compute for diagonal matrices. Hence, we would like (if possible) to convert matrix A into a diagonal matrix. Suppose
More informationGeneralized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este
More informationCounting Matrices Over a Finite Field With All Eigenvalues in the Field
Counting Matrices Over a Finite Field With All Eigenvalues in the Field Lisa Kaylor David Offner Department of Mathematics and Computer Science Westminster College, Pennsylvania, USA kaylorlm@wclive.westminster.edu
More informationSome Formulas for the Principal Matrix pth Root
Int. J. Contemp. Math. Sciences Vol. 9 014 no. 3 141-15 HIKARI Ltd www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.4110 Some Formulas for the Principal Matrix pth Root R. Ben Taher Y. El Khatabi and
More informationA Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.2.3 A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! Ira M. Gessel 1 and Guoce Xin Department of Mathematics Brandeis
More informationDefinition: An n x n matrix, "A", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X 1 A X
DIGONLIZTION Definition: n n x n matrix, "", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X X D. Theorem: n n x n matrix, "", is diagonalizable
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 informationWavelets and Linear Algebra
Wavelets and Linear Algebra () (05) 49-54 Wavelets and Linear Algebra http://wala.vru.ac.ir Vali-e-Asr University of Rafsanjan Schur multiplier norm of product of matrices M. Khosravia,, A. Sheikhhosseinia
More informationELA ON A SCHUR COMPLEMENT INEQUALITY FOR THE HADAMARD PRODUCT OF CERTAIN TOTALLY NONNEGATIVE MATRICES
ON A SCHUR COMPLEMENT INEQUALITY FOR THE HADAMARD PRODUCT OF CERTAIN TOTALLY NONNEGATIVE MATRICES ZHONGPENG YANG AND XIAOXIA FENG Abstract. Under the entrywise dominance partial ordering, T.L. Markham
More informationThe Jordan Normal Form and its Applications
The and its Applications Jeremy IMPACT Brigham Young University A square matrix A is a linear operator on {R, C} n. A is diagonalizable if and only if it has n linearly independent eigenvectors. What happens
More informationm We can similarly replace any pair of complex conjugate eigenvalues with 2 2 real blocks. = R
1 RODICA D. COSTIN Suppose that some eigenvalues are not real. Then the matrices S and are not real either, and the diagonalization of M must be done in C n. Suppose that we want to work in R n only. Recall
More informationMATH JORDAN FORM
MATH 53 JORDAN FORM Let A,, A k be square matrices of size n,, n k, respectively with entries in a field F We define the matrix A A k of size n = n + + n k as the block matrix A 0 0 0 0 A 0 0 0 0 A k It
More informationLemma 8: Suppose the N by N matrix A has the following block upper triangular form:
17 4 Determinants and the Inverse of a Square Matrix In this section, we are going to use our knowledge of determinants and their properties to derive an explicit formula for the inverse of a square matrix
More informationarxiv: v1 [math.nt] 9 May 2017
The spectral norm of a Horadam circulant matrix Jorma K Merikoski a, Pentti Haukkanen a, Mika Mattila b, Timo Tossavainen c, arxiv:170503494v1 [mathnt] 9 May 2017 a Faculty of Natural Sciences, FI-33014
More informationJORDAN NORMAL FORM. Contents Introduction 1 Jordan Normal Form 1 Conclusion 5 References 5
JORDAN NORMAL FORM KATAYUN KAMDIN Abstract. This paper outlines a proof of the Jordan Normal Form Theorem. First we show that a complex, finite dimensional vector space can be decomposed into a direct
More informatione j = Ad(f i ) 1 2a ij/a ii
A characterization of generalized Kac-Moody algebras. J. Algebra 174, 1073-1079 (1995). Richard E. Borcherds, D.P.M.M.S., 16 Mill Lane, Cambridge CB2 1SB, England. Generalized Kac-Moody algebras can be
More informationIr 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 informationAn Analytic Approach to the Problem of Matroid Representibility: Summer REU 2015
An Analytic Approach to the Problem of Matroid Representibility: Summer REU 2015 D. Capodilupo 1, S. Freedman 1, M. Hua 1, and J. Sun 1 1 Department of Mathematics, University of Michigan Abstract A central
More informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 430 (2009) 532 543 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa Computing tight upper bounds
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationarxiv: v1 [cs.sc] 29 Jul 2009
An Explicit Construction of Gauss-Jordan Elimination Matrix arxiv:09075038v1 [cssc] 29 Jul 2009 Yi Li Laboratory of Computer Reasoning and Trustworthy Computing, University of Electronic Science and Technology
More informationTwo finite forms of Watson s quintuple product identity and matrix inversion
Two finite forms of Watson s uintuple product identity and matrix inversion X. Ma Department of Mathematics SuZhou University, SuZhou 215006, P.R.China Submitted: Jan 24, 2006; Accepted: May 27, 2006;
More informationComputationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:
Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication
More information0.1 Tangent Spaces and Lagrange Multipliers
01 TANGENT SPACES AND LAGRANGE MULTIPLIERS 1 01 Tangent Spaces and Lagrange Multipliers If a differentiable function G = (G 1,, G k ) : E n+k E k then the surface S defined by S = { x G( x) = v} is called
More informationA property concerning the Hadamard powers of inverse M-matrices
Linear Algebra and its Applications 381 (2004 53 60 www.elsevier.com/locate/laa A property concerning the Hadamard powers of inverse M-matrices Shencan Chen Department of Mathematics, Fuzhou University,
More informationThe Fibonacci Identities of Orthogonality
The Fibonacci Identities of Orthogonality Kyle Hawins, Ursula Hebert-Johnson and Ben Mathes January 14, 015 Abstract In even dimensions, the orthogonal projection onto the two dimensional space of second
More informationSome bounds for the spectral radius of the Hadamard product of matrices
Some bounds for the spectral radius of the Hadamard product of matrices Guang-Hui Cheng, Xiao-Yu Cheng, Ting-Zhu Huang, Tin-Yau Tam. June 1, 2004 Abstract Some bounds for the spectral radius of the Hadamard
More informationOn the spectra of striped sign patterns
On the spectra of striped sign patterns J J McDonald D D Olesky 2 M J Tsatsomeros P van den Driessche 3 June 7, 22 Abstract Sign patterns consisting of some positive and some negative columns, with at
More informationA note on restricted Lie supertriple systems
Advances in Theoretical and Applied Mathematics. ISSN 0973-4554 Volume 10, Number 1 (2015), pp. 1 13 c Research India Publications http://www.ripublication.com/atam.htm A note on restricted Lie supertriple
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 informationWAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME)
WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME ITAI BENJAMINI, GADY KOZMA, LÁSZLÓ LOVÁSZ, DAN ROMIK, AND GÁBOR TARDOS Abstract. We observe returns of a simple random wal on a finite graph to a fixed node,
More informationRecall : Eigenvalues and Eigenvectors
Recall : Eigenvalues and Eigenvectors Let A be an n n matrix. If a nonzero vector x in R n satisfies Ax λx for a scalar λ, then : The scalar λ is called an eigenvalue of A. The vector x is called an eigenvector
More informationL(C G (x) 0 ) c g (x). Proof. Recall C G (x) = {g G xgx 1 = g} and c g (x) = {X g Ad xx = X}. In general, it is obvious that
ALGEBRAIC GROUPS 61 5. Root systems and semisimple Lie algebras 5.1. Characteristic 0 theory. Assume in this subsection that chark = 0. Let me recall a couple of definitions made earlier: G is called reductive
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationKey words. Strongly eventually nonnegative matrix, eventually nonnegative matrix, eventually r-cyclic matrix, Perron-Frobenius.
May 7, DETERMINING WHETHER A MATRIX IS STRONGLY EVENTUALLY NONNEGATIVE LESLIE HOGBEN 3 5 6 7 8 9 Abstract. A matrix A can be tested to determine whether it is eventually positive by examination of its
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 informationMULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng
Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationElementary Row Operations on Matrices
King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix
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 information