ON THE COMPUTATION OF THE MINIMAL POLYNOMIAL OF A POLYNOMIAL MATRIX
|
|
- Nancy Stewart
- 6 years ago
- Views:
Transcription
1 Int J Appl Math Comput Sci, 25, Vol 5, No 3, ON THE COMPUTATION OF THE MINIMAL POLYNOMIAL OF A POLYNOMIAL MATRIX NICHOLAS P KARAMPETAKIS, PANAGIOTIS TZEKIS Department of Mathematics, Aristotle University of Thessaloniki Thessaloniki 546, Greece karampet@mathauthgr The main contribution of this work is to provide two algorithms for the computation of the minimal polynomial of univariate polynomial matrices The first algorithm is based on the solution of linear matrix equations while the second one employs DFT techniques The whole theory is illustrated with examples Keywords: minimal polynomial, discrete Fourier transform, polynomial matrix, linear matrix equations Introduction It is well known from the Cayley Hamilton theorem that every matrix A R r r satisfies its characteristic equation (Gantmacher, 959), ie, if p (s) :det(si r A) s r +p s r + +p r,then p (A) The Cayley Hamilton theorem is still valid for all cases of matrices over a commutative ring (Atiyah and McDonald, 964), and thus for multivariable polynomial matrices Another form of the Cayley-Hamilton theorem, also known as the relative Caley-Hamilton theorem, is given in terms of the fundamental matrix sequence of the resolvent of the matrix, ie, if (si r A) i Φ is i then Φ k +p Φ k + + p r Φ k r The Cayley-Hamilton theorem was investigated for the matrix pencil case A (s) A + A s in (Mertzios and Christodoulou, 986), and the respective relative Cayley-Hamilton theorem in (Lewis, 986) The Cayley-Hamilton theorem was extended to matrix polynomials (Fragulis, 995; Kitamoto, 999; Yu and Kitamoto, 2), to standard and singular bivariate matrix pencils (Givone and Roesser, 973; Ciftcibaci and Yuksel, 982; Kaczorek, 995a; 989; Vilfan, 973), M-D matrix pencils in (Gałkowski, 996; Theodorou, 989) and n-d polynomial matrices (Kaczorek, 25) The Cayley-Hamilton theorem was also extended to non-square matrices, nonsquare block matrices and singular 2D linear systems with non-square matrices (Kaczorek, 995b; 995c; 995d) The reason behind the interest in the Caley-Hamilton theorem is its applications in control systems, ie, the calculation of controllability and observability grammians and the state-transition matrix, electrical circuits, systems with delays, singular systems, 2-D linear systems, the calculation of the powers of matrices and inverses, etc Of particular importance for the determination of the characteristic polynomial of a polynomial matrix A (s) A + A s + + A q s q R r r [s] are: (a) the Faddeev- Leverrier algorithm (Faddeev and Faddeeva, 963; Helmberg et al, 993) which is fraction free and needs r 3 (r ) polynomial multiplications, (b) the CHTB method presented in (Kitamoto, 999), which needs r 3 (q +) polynomial multiplications (its shortcomings are that it cannot be used for a polynomial matrix A (s) when A has multiple eigenvalues, and it needs to compute first the eigenvalues and eigenvectors of A ), and (c) the CHACM method presented in (Yu and Kitamoto, 2), which needs 7 2 r4 +O ( r 3) polynomial multiplications (a CHTB method given with an artifical constant matrix in order to release the restrictions of the CHTB method, which needs no condition on the given matrix, does not have to solve any eigenvalue problem and is fraction free) Except for the characteristic polynomial of a constant matrix, say p (s), with the nice property p (A), there is also another polynomial, known as the minimal polynomial, say m (s), which is the least degree monic polynomial that satisfies the equation m (A) (Gantmacher, 959) Since the minimal polynomial has a lower degree than the characteristic polynomial, it helps us to solve faster problems such as the computation of the inverse or power of a matrix A number of algorithms have been proposed for the computation of the minimal polynomial of a constant matrix (Augot and Camion, 997), but there is not much interest in polynomial matrices of one or more variables Therefore, the aim of this work is to propose two algorithms for the computation of the minimal polynomial of univariate polynomial matrices The first one is presented
2 34 in Section 2 and is based on the solution of linear matrix equations, while the second one is based on Discrete Fourier Transform (DFT) techniques and is presented in Section 3 The proposed algorithms are illustrated via examples 2 Computation of the Minimal Polynomial of Univariate Polynomial Matrices Consider the polynomial matrix q A(s) A i s i R r r [s], () i q is the greatest power of s in A (s) Definition Every polynomial p (z,s)z p + p (s) z p + + p p (s) for which p ( A(s),s ) A (s) p + p (s) A (s) p + + p p (s) I r (2) is called the annihilating polynomial for the polynomial matrix A (s) R r r [s] The monic annihilating polynomial with a lower degree in z is called the minimal polynomial It is well known that the characteristic polynomial p (z,s) det(zi r A (s)) is an annihilating polynomial, but not necessarily a minimal polynomial Example Let A (s) s s s Then p (z,s) det(zi 3 A (s)) det and z s z s z s (z s) 3 z 3 3sz 2 +3s 2 z s 3 p ( A ( s ),s ) (A (s) si 3 ) 3 3 3,3 NP Karampetakis and P Tzekis The coefficients of the characteristic polynomial can be computed in a recursive way by an algorithm presented in (Fragulis et al, 99; Kitamoto, 999; Yu and Kitamoto, 2) As we shall see below, the characteristic polynomial of the above example is not the only polynomial of the third order that satisfies (2), and does not coincide with the minimal polynomial Let now B(s) p B i s i R r r [s], i p is the greatest power of s in B (s) Then the productof B(s)A(s) is given by ( p+q l ) B(s)A(s) B i A l i s l l i Note that the coefficientmatriceswith indicesgreater than p (resp q) for B (s) (resp for A (s)) aretakento be zero If B (s) Φ, : I r,then ( q q l ) A(s) : Φ,l s l Φ, A(s) Φ,i A l i s l, l Φ,l, l, and thus Φ,l l l Φ,i A l i A l i i Similarly, if we set B (s) q i Φ,is i A (s), Φ,i A i,i,,,q,then 2q A 2 (s) : Φ 2,l s l A(s)A(s) and thus l 2q l Φ 2,l ( l ) Φ,i A l i s l, i l Φ,i A l i i In the general case, Φ k,i is the matrix coefficient of s i in the matrix A (s) k,wehave I r if k, A k (s) kq (3) Φ k,l s l if k, Φ k,l l l Φ k,i A l i with l,,,kq (4) i and Φ,l A l, Φ, I r
3 On the computation of the minimal polynomial of a polynomial matrix 34 Let now the minimal polynomial of A (s) be of the form p (z,s)z m + p m (s)z m + + p (s)z + p (s), m r, with p i (s) (m i)q k Then (2) can be rewritten as p i,k s k, p i,k R (5) p (A(s),s)A(s) m + p m (s)a(s) m or, equivalently, + + p (s)a(s)+p (s)i r r,r, (6) p m (s)a(s) m + + p (s)a(s) Equation (7) can be rewritten as + p (s)i r A(s) m (7) mq mq f i s i Φ m,i s i (8) i i Using (3), (4) and (7) in (8), we get the formula f k m max(,k iq) i jmin(k,(m i)q) Φ m i,j p m i,k j (9) for k,,,mqdefine now the matrices f Φ m, f F m, Φ Φ m, m, f mq P m p m, I r p m, I r p m,q I r p m 2, I r p m 2,2q I r p, I r p,mq I r Φ m,mq, and Φ m, cf Eqn (), n r (mq +), m r m (iq +) i ( ) qm(m +)+m r, 2 and Φ, I r, Φ,i A i From (9) and () we have F m Φ m P m Φ m () Let Φ m i be the matrix that contains i mod r columns of the matrix Φ m and Ki m be the matrix that contains i columns of the matrix Φ m Then () can be rewritten as p m, Φ m Φ m 2 Φ m r {{ F m p m, p m,q p m 2, p m 2,2q p, p,mq {{ P m K m K2 m K m r {{ K m, (2) F m R n2 m2,n 2 r 2 (qm +),m 2 qm(m +)+m, with m r Note that 2 λ n 2 m 2 2 r 2 (qm +) qm(m +)+m r(qm +) q(m +)+ r(qm +) 2 (mq +)+(q +) rm, m 2 r(qm +) 2m 2(mq +) and thus, in general, the number of rows of F m is greater than or equal to the number of its columns, with equality in the case r m Note that the relation (9) and the matrices presented in (), and therefore in (2), are used in (Fragulis, 995) for the computation of the characteristic polynomial of a polynomial matrix, but in a wrong form By using known numerical procedures, such as the Gauss elimination method, the QR factorization or the Cholesky factorization, we can easily determine the values of p i,j and therefore the polynomials p i (s),i,,,m It is easily checked that the upper bound for m is r, ie, the degree of the characteristic polynomial An algorithm for the computation of the
4 342 NP Karampetakis and P Tzekis Φ m Φ m, Φ m, Φ m, Φ m,2 Φ m, Φ m, Φ m,q Φ m,q Φ m,q 2 Φ m,(m )q Φ m,(m )q Φ m,(m )q 2 Φ m,(m 2)q Φ m,(m )q Φ m,(m )q Φ m,(m 2)q+ Φ m,(m )q Φ m,(m 2)q+2 {{ Φ m,(m )q q+ Φ 2, Φ 2, Φ 2, Φ 2,2q Φ 2,2q Φ 2,2q Φ 2, Φ 2,2q Φ, Φ, Φ, Φ,q Φ,q Φ,q Φ, Φ,q Φ 2,2q Φ,q {{{{ (m 2)q+ (m )q+ Φ, Φ, Φ, R n m, () Φ, {{ Φ, mq+
5 On the computation of the minimal polynomial of a polynomial matrix 343 minimal polynomial or otherwise for the coefficients p i,j for i,,,m and j,,,mq is given below Algorithm Computation of the minimal polynomial of a polynomial matrix Step x Step 2 Do x x + Define the matrix Φ x R n m (see ()), with n r (xq +)and m r x i (iq +) Rewrite the equations Φ x P x Φ x as F x P x K x (see (2)) While NOT F x P x K x has a solution Step 3 The coefficients of the minimal polynomial are given by the solution of the system F x P x K x The above algorithm can be easily modified in order to find an annihilating polynomial of the same degree with the characteristic polynomial of the polynomial matrix A (s) as follows: Algorithm 2 Computation of an annihilating polynomial of the same degree with the characteristic polynomial of a polynomial matrix Step Construct the matrices Φ i,, Φ i,,, Φ i,iq, i,,,r Step 2 Define the matrix Φ r (see ()) Step 3 Construct the matrices F r and K r that contain i mod r columns of Φ r and Φ, respectively Step 4 The coefficientsofthe annihilatingpolynomialare given by the solution of the system (2) with m r, ie, F r P r K r The main advantage of the above method of determining of the coefficients p i,j is the use of numerically stable procedures for the solution of the linear system (2), while the main disadvantage is the use of largescale matrices The upper bound for the complexity of the algorithm for the computation of the minimal polynomial is O ( q3 r 9) (q is the greatest power of s in A (s), r stands for the dimension of the matrix A (s) ) and this is the case the minimal polynomial coincides with the characteristic polynomial The lower bound for the complexity of the above algorithm is O ( 2 q3 r 4) The complexity can also be reduced by using fast matrix multiplication techniques (Coppersmith and Winograd, 99) with the complexity O ( n 2376) and fast linear matrix solvers that exploit the sparsity of the matrices, ie, a conjugate-gradient linear system solver, with the complexity O ( n 2) instead of O ( n 3), which we have used for the computation of the complexity of Algorithm 2 Now, if we take into account the fact that the upper bound for the complexity of polynomial multiplication of polynomials with degree at most rq is O (rq log (rq)) (by using FFT transforms), then the CHACM method for the computation of the characteristic polynomial (Yu and Kitamoto, 2) needs r 3 (q +) polynomial multiplications or otherwise the upper bound for its complexity is O ( r 4 q (q +)log(rq) ), which is better than those in Algorithm 2 However, Algorithm 2 may gives better results than those in (Yu and Kitamoto, 2) in the case the minimal polynomial has a much smaller degree in z than the characteristic polynomial Example 2 Let A (s) s s s + {{ A Step x Step 2 x Define the compound matrices [ I3 Φ P {{ A I 3 [ ] p, I 3, Φ p, I 3 ] R 6 6, s R [s] 3 3 [ ] Φ,, Φ, Φ, A, Φ, A We can easily check that the matrix equation Φ P Φ or equivalently P Φ has no solution and thus we proceed to the next step ie x x + Define the compound matrices Φ 2 Φ, I 3 Φ, Φ, I 3 Φ, I 3 R 9 5,
6 344 P 2 p, I 3 p, I 3 p, I 3 p, I 3 p,2 I 3, Φ 2 Φ 2, Φ 2,, Φ 2,2 Φ 2,, Φ 2, 2, Φ 2,2 Rewrite the equations Φ 2 P 2 Φ 2 as F 2 P 2 K 2 (see (2)) with p, p, F 2, P 2 p, p,, K 2 2, p,2 and get the unique solution p, p, p, p, p,2 2 Step 3 The coefficients of the minimal polynomial are given by the solution of the system F 2 P 2 K 2 p (z,s) (z s) 2 z 2 2sz + s 2 NP Karampetakis and P Tzekis In case we would like to apply the Algorithm 2, then we need to proceed one step more ie x r 3Define the compound matrices Φ 2, Φ, I 3 Φ 3 Φ 2, Φ 2, Φ, Φ, I 3 Φ 2,2 Φ 2, Φ, Φ, I 3 Φ 2,2 Φ, I p,3 I 3 I I I I p 2, I 3 p 2, I 3 p, I 3 p, I 3 P 3 p,2 I 3, Φ 3 p, I 3 p, I 3 p,2 I 3 Φ 3,3 Φ 3,, Φ 3,2, Φ 3, Φ 3, Φ 3, Φ 3,2 Φ 3,3 3,
7 On the computation of the minimal polynomial of a polynomial matrix 345 Then we construct the matrices F 3, K 3 and P 3 as defined above: F P 3 p 2, p 2, p, p, p,2 p, p, p,2 p,3,, K 3 3, and solve the equations F 3 P 3 K 3 The coefficients of the annihilating polynomial of degree 3 in z are given by the solution of the above system of equations: 6 z + 3 z 4 6 z 8 p 2, p 4 2, 3 6 z z 5 6 z 9 p, p, 3 z 2 3 z z 8 P 3 p,2 p, z z z 9 p, p,2 p,3 6 z + 3 z 4 6 z z z 5 6 z 9 Thus the annihilating polynomial of order 3 in z will have the following form: [( p (z,s) z z z 5 ) 6 z 9 s ( + 6 z + 3 z 4 )] 6 z 8 z 2 [( z z 5 + ) 3 z 9 s 2 ( + 3 z 2 3 z 4 + ) ] 3 z 8 s z [( z z 5 ) 6 z 9 s 3 ( + 6 z + 3 z 4 ) 6 z 8 s 2], or, equivalently, if we set x z 2 2z 5 + z 9 and y z 2z 4 + z 8, p ( z,s ) (( z ) ( 6 x s + )) 6 y z 2 (( ) ( ) ) 3 x s y s z (( + ) ( 6 x s 3 + )s 6 y 2) (z s) 2( (6 + x)s + y ) z 6 Note that (a) p (A (s),s) and (b) the characteristic polynomial is a special case of p(z,s) for x, y Any other polynomial of the form p (z,s) (z s) 2 (z a (x,x 2,,x n,s)) is also an annihilating polynomial of degree 3 in z since, as we have shown above, the minimal polynomial of the matrix A (s) is (s z) 2
8 346 3 DFT Calculation of a Minimal Polynomial The main disadvantage of the algorithm presented in the previous section is its complexity In order to overcome this difficulty, we canuseothertechniquessuchas interpolation methods Schuster and Hippe (992) use interpolation techniques to find the inverse of a polynomial matrix The speed of interpolation algorithms can be increased by using Discrete Fourier Transforms (DFT) techniques or better Fast Fourier Transforms (FFT) Some of the advantages of DFT-based algorithms are that there are very efficient algorithms available both in software and hardware, and that they are parallel in nature (through symmetric multiprocessing or other techniques) Paccagnella and Pierobon (976) use FFT methods for the computation of the determinant of a polynomial matrix In this section we present an algorithm based on the Discrete Fourier Transform (DFT) which is by an order of magnitude faster than the algorithm presented in the previous section Multidimensional Fourier Transforms arise very frequently in many scientific fields such as image processing, statistics, etc Let us now present the strict definition of a DFT pair Consider the finite sequences X(k,k 2 ) and X(r,r 2 ), k i,r i,,,m i In order for the sequence X(k,k 2 ) and X(r,r 2 ) to constitute a DFT pair, the following relations should hold: X(r,r 2 ) X(k,k 2 ) R M M 2 k k 2 M r r 2 X(k,k 2 )W kr W k2r2, (3) M 2 X(r,r 2 )W kr W k2r2, (4) W i e 2πj M i +, i, 2, (5) R (M +)(M 2 +), (6) and X, X are discrete argument matrix-valued functions The relation (3) is the forward Fourier transform of X(k,k 2 ), while (4) is the inverse Fourier transform of X(r,r 2 ) The great advantage of FFT methods is their reduced complexity The complexity of D DFT on a matrix M R R is O(R 2 ), while the FFT has a complexity of O(R log R) Similarly, the complexity of the DFT of a matrix M R m m2 is O( 2 i m2 i ) which, using the FFT, reduces to O(( 2 i m i)( 2 i log m i) Theinverse DFT is of the same complexity as the forward one In the following, we propose a new algorithm for the calculation of the minimal polynomial of A (s) using discrete Fourier transforms From (5) it is easily seen that the NP Karampetakis and P Tzekis greatest powers of the variables s and z in the minimal polynomial p(s, z) are deg z (p(z,s)) b : m ( r), deg s (p(z,s)) b : mq ( rq) Thus, the polynomial p(z,s) can be written as p(z,s) b b k k (7) (p kk ) ( z k s k) (8) and numerically computed via interpolation using the following R points: u i (r j )W rj i, i, and r j,,,b i, (9) W i e 2πj b i +, (2) R (b +)(b +) (2) In order to evaluate the coefficients p kk,define p rr p (u (r ),u (r )), (22) weusean O ( r 3) algorithm for the computation of the minimal polynomial of the above constant matrix A(u (r )) (Augot and Camion, 997) From (8), (2) and (22) we get p rr b b l l ( ) (p ll ) W rl W rl Notice that [p ll ] and [ p rr ] form a DFT pair and thus using (4) we derive the coefficients of (8), ie, p ll R b b r r p rr W rl W rl, (23) l i,,b i and i, Having in mind the above theoretical deliberations, we will continue by describing the algorithm as an outline for computation Algorithm 3 DFT computation of the minimal polynomial Step Calculate the number of interpolation points b i using (7) Step 2 Compute R points u i (r j ) for i, and r j,,,b i in (9) Step 3 Determine the values at u (r ) of the minimal polynomials of the constant matrices A (u (r )) and thus construct the values p rr in (22) Step 4 Use the inverse DFT (23) for the points p rr in order to construct the values p ll
9 On the computation of the minimal polynomial of a polynomial matrix 347 The above algorithm can also be used for the computation of the characteristic polynomial of a matrix polynomial by making necessary changes in Step 3 (the computation of characteristic polynomial of A (u (r )) instead of the minimal polynomial) The upper bound for the complexity of the above algorithm is O ( r 4 q 2) if we use DFT techniques or O ( r 4 q log (q) ) if we use FFT techniques, and is better than the CHACM method for the characteristic polynomial of A (s) while being comparable to Algorithm 2 when the minimal polynomial has a much smaller degree in z than the characteristic polynomial Example 3 Consider the polynomial matrix A (s) of Example 2 Then by applying Algorithm 3 we have the following results: Step Calculate the number of interpolation points b i by (7) Step 2 Compute R b deg z p(z,s) r 3, b deg s a(z,s) rq 3 (b i +)(3+)(3+)6 i points u i (r j )W rj i,w i e 2πj b i +,i, and r j,,,b j in (9) We get u () W, u () W e 2πj 3+ e πj 2, u (2) W 2 2πj 2 e 3+ e πj, u (3) W 3 2πj 3 e 3+ e 3πj 2, u () W, u () W e 2πj 3+ e 2πj 4, u (2) W 2 2πj 2 e 3+ e 4πj 4, u (3) W 3 2πj 3 e 3+ e 6πj 4 Step 3 Determine the minimal polynomials of the constant matrices A (u (r )): p (z,u ()) z 2 2z +, p (z,u ()) z 2 +2jz, p (z,u (2)) z 2 +2z +, p (z,u (3)) z 2 2jz, and then the values of each polynomial at u (r ), p, p (u (),u ()) ( z 2 2z + ) z, p, p (u (),u ()) ( z 2 2z + ) 2j, ze πj/2 p 2, p (u (2),u ()) ( z 2 2z + ) 4, ze πj p 3, p (u (3),u ()) ( z 2 2z + ) 2j, ze 3πj/2 p, p (u (),u ()) ( z 2 +2jz ) z 2j, p, p (u (),u ()) ( z 2 +2jz ), ze πj/2 p 2, p (u (2),u ()) ( z 2 +2jz ) 2j, ze πj p 3, p (u (3),u ()) ( z 2 +2jz ) 4, ze 3πj/2 p,2 p (u (),u (2)) ( z 2 +2z + ) 4, z p,2 p (u (),u (2)) ( z 2 +2z + ) 2j, ze πj/2 p 2,2 p (u (2),u (2)) ( z 2 +2z + ), ze πj p 3,2 p (u (3),u (2)) ( z 2 +2z + ) 2j, ze 3πj/2 p,3 p (u (),u (3)) ( z 2 2jz ) 2j, z p,3 p (u (),u (3)) ( z 2 2jz ) 4, ze πj/2 p 2,3 p (u (2),u (3)) ( z 2 2jz ) 2j, ze πj p 3,3 p (u (3),u (3)) ( z 2 2jz ), ze 3πj/2 and thus construct the values p rr in (22)
10 348 Step 4 Use the inverse DFT (23) for the points p rr in order to construct the values p ll 3 r p r r W rl W rl 6 3 r p,, p,, p,2, p,3, p,, p, 2, p,2, p,3, p 2,, p 2,, p 2,2, p 2,3, p 3,, p 3,, p 3,2, p 3,3, and thus the minimal polynomial is p (z,s) s 2 2sz + z 2 (z s) 2 In the case when we are interested in the characteristic polynomial of the matrix A (s), we have to change Steps 3 and 4 as follows: Step 3a Determine the characteristic polynomials of the constant matrices A (u (r )): p (z,u ()) z 3 3z 2 +3z, p (z,u ()) z 3 3z 2 e 2 πj +3ze πj e 3 2 πj, p (z,u (2)) z 3 3z 2 e πj +3ze 2πj e 3πj, p (z,u (3)) z 3 3z 2 e 3 2 πj +3ze 3πj e 9 2 πj, and then the values of each polynomial at u (r ), 6 3 r p,, p, 2 2j, p 2, 8, p 3, 2+2j, p, 2+2j, p,, p 2, 2+2j, p 3, 8j, p,2 8, p,2 2 2j, p 2,2, p,3 2 2j, p 3,2 2+2j, p,3 8j, p 2,3 2 2j, p 3,3, and thus construct the values p rr in (22) Step 4a Use the inverse DFT (23) for the points p rr in order to construct the values p ll 3 r p r r W rl W rl, p,, p,, p,2, p,3, p,, p,, p,2 3, p,3, p 2,, p 2, 3, p 2,2, p 2,3, p 3,, p 3,, p 3,2, p 3,3, and thus the characteristic polynomial is NP Karampetakis and P Tzekis p (z,s)z 3 3z 2 s +3zs 2 s 3 Note that, when using the DFT method or the CHACM method, it is not easy to find the family of annihilating polynomials of degree 3 in z, as we have already done in Example 2 4 Conclusions Two algorithms for the computation of the minimal polynomial and the characteristic polynomial of univariate matrices have been developed The proposed algorithms are easily implemented on a digital computer and are very useful in many problems, such as the computation of the powers of polynomial matrices, the evaluation of the controllability and observability grammians of polynomial matrix descriptions, the calculation of the state-transition matrix, which is used for the evaluation of the solution of homogeneous matrix differential equations of the form A (ρ) β (t) with A (ρ) R [ρ] r r The algorithm presented in Section 2 is based on the solution of linear matrix equations Its main advantage is that (a) it creates a family of annihilating polynomials of the same degree with the characteristic polynomial, and (b) has a complexity comparable with the DFT method presented in Section 3 in the case the degree of the minimal polynomial in z is lower enough in contrast to the degree of the characteristic polynomial Its main disadvantage is that the upper bound for its complexity is big enough since it uses large-scale matrices The algorithm presented in Section 3 is based on DFT techniques and has a lower upper bound complexity than the ones in Section 2 An extension of these algorithms to the n-variable case can be easily done (Tzekis and Karampetakis, 25) Acknowledgments The authors would like to thank the reviewers for their valuable comments References Atiyah MF and McDonald IG (964): Introduction to Commutative Algebra Reading, MA: Addison-Wesley Augot D and Camion P (997): On the computation of minimal polynomials, cyclic vectors, and Frobenius forms Lin Alg App, Vol 26, pp 6 94 Ciftcibaci T and Yuksel O (982): On the Cayley-Hamllton theorem for two-dimensional systems IEEE Trans Automat Contr, Vol 27, No, pp 93 94
11 On the computation of the minimal polynomial of a polynomial matrix 349 Coppersmith D and Winograd S (99): Matrix multiplication via arithmetic progressions J Symb Comput, Vol 9, No3, pp Faddeev DK and Faddeeva VN (963): Computational Methods of Linear Algebra San Francisco: Freeman Fragulis GF, Mertzios B and Vardulakis AIG (99): Computation of the inverse of a polynomial matrix and evaluation of its Laurent expansion Int J Contr, Vol 53, No 2, pp Fragulis GF (995): Generalized Cayley-Hamilton theorem for polynomial matrices with arbitrary degree IntJ Contr, Vol 62, No 6, pp Gałkowski K (996): Matrix description of polynomials multivariable Lin Alg Applic, Vol 234, pp Gantmacher FR (959): The Theory of Matrices New York: Chelsea Givone DD and Roesser RP (973): Minimization of multidimensional linear iterative circuits IEEE Trans Comput, Vol C-22 pp Helmberg G, Wagner P and Veltkamp G (993): On Faddeev- Leverrier s method for the computation of the characteristic polynomial of a matrix and of eigenvectors Lin Alg Applic, Vol 85, pp Kaczorek T (989): Existence and uniqueness of solutions and Cayley-Hamilton theorem for a general singular model of 2-D systems Bull Polish Acad Sci Techn Sci, Vol 37, No 5-6, pp Kaczorek T (995a): An extension of the Cayley-Hamilton theorem for 2-D continuous-discrete linear systems Appl Math Comput Sci, Vol 4, No 4, pp Kaczorek T (995b): Generalization of the Cayley Hamilton theorem for nonsquare matrices Proc Int Conf Fundamentals of Electrotechnics and Circuit Theory XVIII- SPETO, Gliwice-Ustroń, Poland, pp Kaczorek T (995c): An extension of the Cayley-Hamilton theorem for nonsquare block matrices and computation of the left and right inverses of matrices Bull Polish Acad Sci Techn Sci, Vol 43, No, pp Kaczorek T (995d): An extension of the Cayley-Hamilton theorem for singular 2D linear systems with non-square matrices Bull Polish Acad Sci, Techn Sci, Vol 43, No, pp Kaczorek T (998): An extension of the Cayley-Hamilton theorem for a standard pair of block matrices Appl Math Comput Sci, Vol 8, No 3, pp 5 56 Kaczorek T (25): Generalization of Cayley-Hamilton theorem for n-d polynomial matrices IEEE Trans Automat Contr, Vol 5, No 5, pp Kitamoto T (999): Efficient computation of the characteristic polynomial of a polynomial matrix IEICE Trans Fundament, Vol E83-A, No5, pp Lewis FL (986): Further remarks on the Cayley-Hamilton theorem and Leverrier s method for the matrix pencil (se A) IEEE Trans Automat Contr, Vol 3, No 9, pp Mertzios B and Christodoulou M (986): On the generalized Cayley-Hamilton theorem IEEE Trans Automat Contr, Vol 3, No 2, pp Paccagnella L E and Pierobon G L (976): FFT calculation of a determinantal polynomial IEEE Trans Automat Contr, Vol 2, No 3, pp 4 42 Schuster A and Hippe P (992): Inversion of polynomial matrices by interpolation IEEE Trans Automat Contr, Vol 37, No 3, pp Theodorou NJ (989): M-dimensional Cayley-Hamilton theorem IEEE Trans Automat Contr, Vol 34, No 5, pp Tzekis P and Karampetakis N (25): On the computation of the minimal polynomial of a two-variable polynomial matrix Proc 4th Int Workshop Multidimensional (nd) Systems, Wuppertal, Germany Victoria J (982): A block Cayley-Hamilton theorem Bull Math Soc Sci Math Roum, Vol 26, No, pp Vilfan B (973): Another proof of the two-dimensional Cayley Hamilton theorem IEEE Trans Comput, Vol C-22, pp 4 Yu B and Kitamoto T (2): The CHACM method for computing the characteristic polynomial of a polynomial matrix IEICE Trans Fundament, E83-A, No7, pp 45 4 Received: 8 March 25 Revised: 9 May 25
Infinite elementary divisor structure-preserving transformations for polynomial matrices
Infinite elementary divisor structure-preserving transformations for polynomial matrices N P Karampetakis and S Vologiannidis Aristotle University of Thessaloniki, Department of Mathematics, Thessaloniki
More informationA q x k+q + A q 1 x k+q A 0 x k = 0 (1.1) where k = 0, 1, 2,..., N q, or equivalently. A(σ)x k = 0, k = 0, 1, 2,..., N q (1.
A SPECTRAL CHARACTERIZATION OF THE BEHAVIOR OF DISCRETE TIME AR-REPRESENTATIONS OVER A FINITE TIME INTERVAL E.N.Antoniou, A.I.G.Vardulakis, N.P.Karampetakis Aristotle University of Thessaloniki Faculty
More informationPOSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS
Int J Appl Math Comput Sci 7 Vol 7 No 65 7 DOI: 478/v6-7-5- POSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical
More informationMINIMAL POLYNOMIALS AND CHARACTERISTIC POLYNOMIALS OVER RINGS
JP Journal of Algebra, Number Theory and Applications Volume 0, Number 1, 011, Pages 49-60 Published Online: March, 011 This paper is available online at http://pphmj.com/journals/jpanta.htm 011 Pushpa
More informationComputing Minimal Polynomial of Matrices over Algebraic Extension Fields
Bull. Math. Soc. Sci. Math. Roumanie Tome 56(104) No. 2, 2013, 217 228 Computing Minimal Polynomial of Matrices over Algebraic Extension Fields by Amir Hashemi and Benyamin M.-Alizadeh Abstract In this
More informationSYMBOLIC IMPLEMENTATION OF LEVERRIER-FADDEEV ALGORITHM AND APPLICATIONS
SYMBOLIC IMPLEMENTATION OF LEVERRIER-FADDEEV ALGORITHM AND APPLICATIONS Predrag S. Stanimirovic* andn. P. Karampetakis** *University of of Niš, Faculty of Philosophy, Department of Mathematics, Ćirila
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 informationCONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2002, Vol.2, No.2, 73 80 CONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS JERZY KLAMKA Institute of Automatic Control, Silesian University of Technology ul. Akademicka 6,
More informationComputing generalized inverse of polynomial matrices by interpolation
Computing generalized inverse of polynomial matrices by interpolation Predrag S. Stanimirović, Marko D. Petković, University of Niš, Department of Mathematics, Faculty of Science, Višegradska 33, 18000
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationNotes on n-d Polynomial Matrix Factorizations
Multidimensional Systems and Signal Processing, 10, 379 393 (1999) c 1999 Kluwer Academic Publishers, Boston. Manufactured in The Netherlands. Notes on n-d Polynomial Matrix Factorizations ZHIPING LIN
More informationAlgebro-geometric aspects of Heine-Stieltjes theory
Dedicated to Heinrich Eduard Heine and his 4 years old riddle Algebro-geometric aspects of Heine-Stieltjes theory Boris Shapiro, Stockholm University, shapiro@math.su.se February, 9 Introduction and main
More informationCAIPORXI INSTITUTh OF TZCHNOLOGY
'754 MA THEMA TICS: E. T. BELL PiLoc. N. A. S. ;It is evident that if M is cyclicly connected, then Lim M(x) = M for every point x of M and conversely. TH1oRZM 13.-If M is a continuous curve, either (1)
More informationFast reversion of formal power series
Fast reversion of formal power series Fredrik Johansson LFANT, INRIA Bordeaux RAIM, 2016-06-29, Banyuls-sur-mer 1 / 30 Reversion of power series F = exp(x) 1 = x + x 2 2! + x 3 3! + x 4 G = log(1 + x)
More informationEXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci., 1, Vol.11, No.4, 957 964 EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS Tadeusz KACZOREK The notions of externally and internally positive time-varying
More informationCharacteristic Polynomial
Linear Algebra Massoud Malek Characteristic Polynomial Preleminary Results Let A = (a ij ) be an n n matrix If Au = λu, then λ and u are called the eigenvalue and eigenvector of A, respectively The eigenvalues
More informationOn Böttcher s mysterious identity
AUSTRALASIAN JOURNAL OF COBINATORICS Volume 43 (2009), Pages 307 316 On Böttcher s mysterious identity Ömer Eğecioğlu Department of Computer Science University of California Santa Barbara, CA 93106 U.S.A.
More informationTopics in linear algebra
Chapter 6 Topics in linear algebra 6.1 Change of basis I want to remind you of one of the basic ideas in linear algebra: change of basis. Let F be a field, V and W be finite dimensional vector spaces over
More informationSOLVING SOLVABLE QUINTICS. D. S. Dummit
D. S. Dummit Abstract. Let f(x) = x 5 + px 3 + qx + rx + s be an irreducible polynomial of degree 5 with rational coefficients. An explicit resolvent sextic is constructed which has a rational root if
More informationThe System of Linear Equations. Direct Methods. Xiaozhou Li.
1/16 The Direct Methods xiaozhouli@uestc.edu.cn http://xiaozhouli.com School of Mathematical Sciences University of Electronic Science and Technology of China Chengdu, China Does the LU factorization always
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 2: Direct Methods PD Dr.
More informationOn Systems of Diagonal Forms II
On Systems of Diagonal Forms II Michael P Knapp 1 Introduction In a recent paper [8], we considered the system F of homogeneous additive forms F 1 (x) = a 11 x k 1 1 + + a 1s x k 1 s F R (x) = a R1 x k
More informationRevision Problems for Examination 1 in Algebra 1
Centre for Mathematical Sciences Mathematics, Faculty of Science Revision Problems for Examination 1 in Algebra 1 Arithmetics 1 Determine a greatest common divisor to the integers a) 5431 and 1345, b)
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationGRE Subject test preparation Spring 2016 Topic: Abstract Algebra, Linear Algebra, Number Theory.
GRE Subject test preparation Spring 2016 Topic: Abstract Algebra, Linear Algebra, Number Theory. Linear Algebra Standard matrix manipulation to compute the kernel, intersection of subspaces, column spaces,
More informationTHE CAYLEY-HAMILTON THEOREM AND INVERSE PROBLEMS FOR MULTIPARAMETER SYSTEMS
THE CAYLEY-HAMILTON THEOREM AND INVERSE PROBLEMS FOR MULTIPARAMETER SYSTEMS TOMAŽ KOŠIR Abstract. We review some of the current research in multiparameter spectral theory. We prove a version of the Cayley-Hamilton
More informationQUADRATIC RESIDUE CODES OVER Z 9
J. Korean Math. Soc. 46 (009), No. 1, pp. 13 30 QUADRATIC RESIDUE CODES OVER Z 9 Bijan Taeri Abstract. A subset of n tuples of elements of Z 9 is said to be a code over Z 9 if it is a Z 9 -module. In this
More informationarxiv:math/ v1 [math.co] 13 Jul 2005
A NEW PROOF OF THE REFINED ALTERNATING SIGN MATRIX THEOREM arxiv:math/0507270v1 [math.co] 13 Jul 2005 Ilse Fischer Fakultät für Mathematik, Universität Wien Nordbergstrasse 15, A-1090 Wien, Austria E-mail:
More informationLOCALLY POSITIVE NONLINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci. 3 Vol. 3 No. 4 55 59 LOCALLY POSITIVE NONLINEAR SYSTEMS TADEUSZ KACZOREK Institute of Control Industrial Electronics Warsaw University of Technology ul. Koszykowa 75 66
More informationConjugate Gradient Method
Conjugate Gradient Method direct and indirect methods positive definite linear systems Krylov sequence spectral analysis of Krylov sequence preconditioning Prof. S. Boyd, EE364b, Stanford University Three
More informationMath 489AB Exercises for Chapter 2 Fall Section 2.3
Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationInvariance properties in the root sensitivity of time-delay systems with double imaginary roots
Invariance properties in the root sensitivity of time-delay systems with double imaginary roots Elias Jarlebring, Wim Michiels Department of Computer Science, KU Leuven, Celestijnenlaan A, 31 Heverlee,
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 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 informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationStability of interval positive continuous-time linear systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 66, No. 1, 2018 DOI: 10.24425/119056 Stability of interval positive continuous-time linear systems T. KACZOREK Białystok University of
More informationLinear algebra and applications to graphs Part 1
Linear algebra and applications to graphs Part 1 Written up by Mikhail Belkin and Moon Duchin Instructor: Laszlo Babai June 17, 2001 1 Basic Linear Algebra Exercise 1.1 Let V and W be linear subspaces
More informationAn Efficient Algorithm for Deriving Summation Identities from Mutual Recurrences
Discrete Mathematics, Algorithms and Applications c World Scientific Publishing Company An Efficient Algorithm for Deriving Summation Identities from Mutual Recurrences BERKELEY R CHURCHILL Department
More informationSECOND ORDER TWO DIMENSIONAL SYSTEMS: COMPUTING THE TRANSFER FUNCTION
Journal of ELECTRICAL ENGINEERING, VOL. 55, NO. 11-12, 2004, 296 300 SECOND ORDER TWO DIMENSIONAL SYSTEMS: COMPUTING THE TRANSFER FUNCTION Gege E. Antoniou Marinos T. Michael In this paper the discrete
More informationProduct distance matrix of a tree with matrix weights
Product distance matrix of a tree with matrix weights R B Bapat Stat-Math Unit Indian Statistical Institute, Delhi 7-SJSS Marg, New Delhi 110 016, India email: rbb@isidacin Sivaramakrishnan Sivasubramanian
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationOn the Faddeev-Sominsky s algorithm
Available online at www.worldscientificnews.com WSN 106 (2018) 238-244 EISSN 2392-2192 SHORT COMMUNICATION On the Faddeev-Sominsky s algorithm H. Torres-Silva 1, J. López-Bonilla 2, *, S. Vidal-Beltrán
More informationControl for stability and Positivity of 2-D linear discrete-time systems
Manuscript received Nov. 2, 27; revised Dec. 2, 27 Control for stability and Positivity of 2-D linear discrete-time systems MOHAMMED ALFIDI and ABDELAZIZ HMAMED LESSI, Département de Physique Faculté des
More informationA new approach to the realization problem for fractional discrete-time linear systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol 64 No 26 DOI: 55/bpasts-26-2 A new approach to the realization problem for fractional discrete-time linear systems T KACZOREK Faculty of
More informationMATHEMATICS 217 NOTES
MATHEMATICS 27 NOTES PART I THE JORDAN CANONICAL FORM The characteristic polynomial of an n n matrix A is the polynomial χ A (λ) = det(λi A), a monic polynomial of degree n; a monic polynomial in the variable
More informationA new parallel polynomial division by a separable polynomial via hermite interpolation with applications, pp
Electronic Journal of Linear Algebra Volume 23 Volume 23 (2012 Article 54 2012 A new parallel polynomial division by a separable polynomial via hermite interpolation with applications, pp 770-781 Aristides
More information-dimensional space. d 1. d k (1 ρ 2 j,j+k;j+1...j+k 1). Γ(2m) m=1. 2 ) 2 4d)/4. Γ d ( d 2 ) (d 2)/2. Γ d 1 (d)
Generating random correlation matrices based on partial correlation vines and the onion method Harry Joe Department of Statistics University of British Columbia Abstract Partial correlation vines and the
More informationNew algebraic decoding method for the (41, 21,9) quadratic residue code
New algebraic decoding method for the (41, 21,9) quadratic residue code Mohammed M. Al-Ashker a, Ramez Al.Shorbassi b a Department of Mathematics Islamic University of Gaza, Palestine b Ministry of education,
More informationCourse Notes: Week 1
Course Notes: Week 1 Math 270C: Applied Numerical Linear Algebra 1 Lecture 1: Introduction (3/28/11) We will focus on iterative methods for solving linear systems of equations (and some discussion of eigenvalues
More informationPolynomial Multiplication over Finite Fields using Field Extensions and Interpolation
009 19th IEEE International Symposium on Computer Arithmetic Polynomial Multiplication over Finite Fields using Field Extensions and Interpolation Murat Cenk Department of Mathematics and Computer Science
More informationCHAPTER XV. Exponential algebra
CHAPTER XV Exponential algebra 15.1. Introduction. We develop some elementary features of the intricate and hyperintricate exponential algebra and in particular the bosonic and fermionic Euler relations
More informationDESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS
Int. J. Appl. Math. Comput. Sci., 3, Vol. 3, No., 39 35 DOI:.478/amcs-3-3 DESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical
More informationGenerating Functions of Partitions
CHAPTER B Generating Functions of Partitions For a complex sequence {α n n 0,, 2, }, its generating function with a complex variable q is defined by A(q) : α n q n α n [q n ] A(q). When the sequence has
More information1 Introduction It will be convenient to use the inx operators a b and a b to stand for maximum (least upper bound) and minimum (greatest lower bound)
Cycle times and xed points of min-max functions Jeremy Gunawardena, Department of Computer Science, Stanford University, Stanford, CA 94305, USA. jeremy@cs.stanford.edu October 11, 1993 to appear in the
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 informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
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 informationarxiv: v1 [math.dg] 1 Jul 2014
Constrained matrix Li-Yau-Hamilton estimates on Kähler manifolds arxiv:1407.0099v1 [math.dg] 1 Jul 014 Xin-An Ren Sha Yao Li-Ju Shen Guang-Ying Zhang Department of Mathematics, China University of Mining
More informationON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS
Journal of Circuits Systems and Computers Vol. 3 No. 5 (2004) 5 c World Scientific Publishing Company ON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS GEORGE E. ANTONIOU Department of Computer
More informationALGORITHMS FOR CENTROSYMMETRIC AND SKEW-CENTROSYMMETRIC MATRICES. Iyad T. Abu-Jeib
ALGORITHMS FOR CENTROSYMMETRIC AND SKEW-CENTROSYMMETRIC MATRICES Iyad T. Abu-Jeib Abstract. We present a simple algorithm that reduces the time complexity of solving the linear system Gx = b where G is
More informationAlgebraic Algorithm for 2D Stability Test Based on a Lyapunov Equation. Abstract
Algebraic Algorithm for 2D Stability Test Based on a Lyapunov Equation Minoru Yamada Li Xu Osami Saito Abstract Some improvements have been proposed for the algorithm of Agathoklis such that 2D stability
More informationa c b b a + c b + c. c b a + c
TRACE AND NORM KEITH CONRAD 1. Introduction Let L/K be a finite extension of fields, with n = [L : K]. We will associate to this extension two important functions L K, called the trace and the norm. For
More informationOn the computation of the Jordan canonical form of regular matrix polynomials
On the computation of the Jordan canonical form of regular matrix polynomials G Kalogeropoulos, P Psarrakos 2 and N Karcanias 3 Dedicated to Professor Peter Lancaster on the occasion of his 75th birthday
More informationNewcriteria for H-tensors and an application
Wang and Sun Journal of Inequalities and Applications (016) 016:96 DOI 10.1186/s13660-016-1041-0 R E S E A R C H Open Access Newcriteria for H-tensors and an application Feng Wang * and De-shu Sun * Correspondence:
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
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 informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationMultivariable ARMA Systems Making a Polynomial Matrix Proper
Technical Report TR2009/240 Multivariable ARMA Systems Making a Polynomial Matrix Proper Andrew P. Papliński Clayton School of Information Technology Monash University, Clayton 3800, Australia Andrew.Paplinski@infotech.monash.edu.au
More informationTHE REAL AND THE SYMMETRIC NONNEGATIVE INVERSE EIGENVALUE PROBLEMS ARE DIFFERENT
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 12, December 1996, Pages 3647 3651 S 0002-9939(96)03587-3 THE REAL AND THE SYMMETRIC NONNEGATIVE INVERSE EIGENVALUE PROBLEMS ARE DIFFERENT
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationOn Connections between the Cauchy Index, the Sylvester Matrix, Continued Fraction Expansions, and Circuit Synthesis
2st International Symposium on Mathematical Theory of Networks Systems July 7-, 24 On Connections between the Cauchy Index, the Sylvester Matrix, Continued Fraction Expansions, Circuit Synthesis Timothy
More informationDynamics of finite linear cellular automata over Z N
Dynamics of finite linear cellular automata over Z N F. Mendivil, D. Patterson September 9, 2009 Abstract We investigate the behaviour of linear cellular automata with state space Z N and only finitely
More informationKnowledge Discovery and Data Mining 1 (VO) ( )
Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory
More informationScientific Computing WS 2018/2019. Lecture 9. Jürgen Fuhrmann Lecture 9 Slide 1
Scientific Computing WS 2018/2019 Lecture 9 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 9 Slide 1 Lecture 9 Slide 2 Simple iteration with preconditioning Idea: Aû = b iterative scheme û = û
More informationCOMMUTATIVE/NONCOMMUTATIVE RANK OF LINEAR MATRICES AND SUBSPACES OF MATRICES OF LOW RANK
Séminaire Lotharingien de Combinatoire 52 (2004), Article B52f COMMUTATIVE/NONCOMMUTATIVE RANK OF LINEAR MATRICES AND SUBSPACES OF MATRICES OF LOW RANK MARC FORTIN AND CHRISTOPHE REUTENAUER Dédié à notre
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Eigenvalue Problems; Similarity Transformations Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Eigenvalue
More informationOn the Cross-Correlation of a p-ary m-sequence of Period p 2m 1 and Its Decimated
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 58, NO 3, MARCH 01 1873 On the Cross-Correlation of a p-ary m-sequence of Period p m 1 Its Decimated Sequences by (p m +1) =(p +1) Sung-Tai Choi, Taehyung Lim,
More informationCLUSTER ALGEBRA ALFREDO NÁJERA CHÁVEZ
Séminaire Lotharingien de Combinatoire 69 (203), Article B69d ON THE c-vectors AND g-vectors OF THE MARKOV CLUSTER ALGEBRA ALFREDO NÁJERA CHÁVEZ Abstract. We describe the c-vectors and g-vectors of the
More informationA Note on the Transcendence of Zeros of a Certain Family of Weakly Holomorphic Forms
A Note on the Transcendence of Zeros of a Certain Family of Weakly Holomorphic Forms Jennings-Shaffer C. & Swisher H. (014). A Note on the Transcendence of Zeros of a Certain Family of Weakly Holomorphic
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More informationRadix-4 Factorizations for the FFT with Ordered Input and Output
Radix-4 Factorizations for the FFT with Ordered Input and Output Vikrant 1, Ritesh Vyas 2, Sandeep Goyat 3, Jitender Kumar 4, Sandeep Kaushal 5 YMCA University of Science & Technology, Faridabad (Haryana),
More informationDiameter of the Zero Divisor Graph of Semiring of Matrices over Boolean Semiring
International Mathematical Forum, Vol. 9, 2014, no. 29, 1369-1375 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.47131 Diameter of the Zero Divisor Graph of Semiring of Matrices over
More informationFundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved
Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural
More informationAn Interpolation Algorithm for List Decoding of Reed-Solomon Codes
An Interpolation Algorithm for List Decoding of Reed-Solomon Codes Kwankyu Lee Department of Mathematics San Diego State University San Diego, USA Email: kwankyu@sogangackr Michael E O Sullivan Department
More informationCHARACTERISTIC AND MINIMAL POLYNOMIALS OF LINEAR CELLULAR AUTOMATA
CHARACTERISTIC AND MINIMAL POLYNOMIALS OF LINEAR CELLULAR AUTOMATA STEVEN LETTIERI, JOHN G. STEVENS AND DIANA M. THOMAS DEPARTMENT OF MATHEMATICAL SCIENCES MONTCLAIR STATE UNIVERSITY UPPER MONTCLAIR, NJ
More informationCONTINUITY OF SOLUTIONS OF RICCATI EQUATIONS FOR THE DISCRETE TIME JLQP
Int. J. Appl. Math. Comput. Sci., 2002, Vol.12, No.4, 539 543 CONTINUITY O SOLUTIONS O RICCATI EQUATIONS OR THE DISCRETE TIME JLQP ADAM CZORNIK, ANDRZEJ ŚWIERNIAK Department of Automatic Control Silesian
More information= AΦ, Φ(0) = I, = k! tk A k = I + ta t2 A t3 A t4 A 4 +, (3) e ta 1
Matrix Exponentials Math 26, Fall 200, Professor David Levermore We now consider the homogeneous constant coefficient, vector-valued initial-value problem dx (1) = Ax, x(t I) = x I, where A is a constant
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 informationPolynomial functions over nite commutative rings
Polynomial functions over nite commutative rings Balázs Bulyovszky a, Gábor Horváth a, a Institute of Mathematics, University of Debrecen, Pf. 400, Debrecen, 4002, Hungary Abstract We prove a necessary
More informationFast reversion of power series
Fast reversion of power series Fredrik Johansson November 2011 Overview Fast power series arithmetic Fast composition and reversion (Brent and Kung, 1978) A new algorithm for reversion Implementation results
More informationQuasi-primary submodules satisfying the primeful property II
Hacettepe Journal of Mathematics and Statistics Volume 44 (4) (2015), 801 811 Quasi-primary submodules satisfying the primeful property II Hosein Fazaeli Moghimi and Mahdi Samiei Abstract In this paper
More informationMatrix Factorization and Analysis
Chapter 7 Matrix Factorization and Analysis Matrix factorizations are an important part of the practice and analysis of signal processing. They are at the heart of many signal-processing algorithms. Their
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 informationStabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks
Control and Cybernetics vol. 39 (2010) No. 3 Stabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks by Tadeusz Kaczorek
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 informationSums of diagonalizable matrices
Linear Algebra and its Applications 315 (2000) 1 23 www.elsevier.com/locate/laa Sums of diagonalizable matrices J.D. Botha Department of Mathematics University of South Africa P.O. Box 392 Pretoria 0003
More informationLow complexity bit-parallel GF (2 m ) multiplier for all-one polynomials
Low complexity bit-parallel GF (2 m ) multiplier for all-one polynomials Yin Li 1, Gong-liang Chen 2, and Xiao-ning Xie 1 Xinyang local taxation bureau, Henan, China. Email:yunfeiyangli@gmail.com, 2 School
More informationHonors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35
Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35 1. Let R be a commutative ring with 1 0. (a) Prove that the nilradical of R is equal to the intersection of the prime
More information