Numerical solution of the eigenvalue problem for Hermitian Toeplitz-like matrices
|
|
- Erick Melton
- 5 years ago
- Views:
Transcription
1 TR-CS-9-1 Numerical solution of the eigenvalue problem for Hermitian Toeplitz-like matrices Michael K Ng and William F Trench July 199 Joint Computer Science Technical Report Series Department of Computer Science Faculty of Engineering and Information Technology Computer Sciences Laboratory Research School of Information Sciences and Engineering
2 This technical report series is published jointly by the Department of Computer Science, Faculty of Engineering and Information Technology, and the Computer Sciences Laboratory, Research School of Information Sciences and Engineering, The Australian National University Please direct correspondence regarding this series to: Technical Reports Department of Computer Science Faculty of Engineering and Information Technology The Australian National University Canberra ACT Australia or send to: A list of technical reports, including some abstracts and copies of some full reports may be found at: Recent reports in this series: TR-CS-9-1 Michael K Ng Blind channel identification and the eigenvalue problem of structured matrices July 199 TR-CS-9-1 Michael K Ng Preconditioning of elliptic problems by approximation in the transform domain July 199 TR-CS-9-11 Richard P Brent, Richard E Crandall, and Karl Dilcher Two new factors of Fermat numbers May 199 TR-CS-9-1 Andrew Tridgell, Richard Brent, and Brendan McKay Parallel integer sorting May 199 TR-CS-9-9 M Manzur Murshed and Richard P Brent Constant time algorithms for computing the contour of maximal elements on the Reconfigurable Mesh May 199 TR-CS-9-8 Xun Qu, Jeffrey Xu Yu, and Richard P Brent socket April 199 A mobile TCP
3 NUMERICAL SOLUTION OF THE EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES MICHAEL K NG AND WILLIAM F TRENCH y Abstract An iterative method based on displacement structure is proposed for computing eigenvalues and eigenvectors of a class of Hermitian Toeplitz{like matrices which includes matrices of the form T T where T is arbitrary Toeplitz matrix, Toeplitz{block matrices and block{toeplitz matrices The method obtains a specic individual eigenvalue (ie, the i-th smallest, where i is a specied integer in [1 :::n]) of an n n matrix at a computational cost of O(n ) operations An associated eigenvector is obtained as a byproduct The method is more ecient than general purpose methods such as the QR algorithm for obtaining a small number (compared to n) of eigenvalues Moreover, since the computation of each eigenvalue is independent of the computation of all other eigenvalues, the method is highly parallelizable Numerical results illustrate the eectiveness of the method Key words Toeplitz matrix, displacement structure, Toeplitz{like matrix, eigenvalue, eigenvector, root{nding AMS subject classications 15A18, 15A5, 5F15 ACSys, Computer Sciences Laboratory, Research School of Information Sciences and Engineering, The Australian National University, Canberra, ACT, Australia mng@cslabanueduau y Department of Mathematics, Trinity University, 15 Stadium Drive, San Antonio, Texas, 81, USA wtrench@trinityedu Research supported by National Science Foundation Grant DMS
4 Michael K Ng and William F Trench 1 Introduction In this paper we consider the eigenvalue problem for an n n Hermitian matrix A n which has displacement structure in the sense that (1) A n Z n, Z n A n = G n H T n where Z n is the shift matrix Z n = G n and H n are in C n, and is small compared to n (For discussions of other types of displacement structure see [8, 1, 11]) The smallest integer for which (1) holds with some G n and H n in C n, is called the fz n Z n g-displacement rank of A n we will call it simply the displacement rank of A n Henceforth we will say that a matrix which satises (1) with small compared to n is a Toeplitz{like matrix There are many ecient direct methods that exploit displacement structure to invert Toeplitz{like matrices, or to solve Toeplitz{like systems A n x = b [, 8, 11] There are also preconditioned conjugate gradient methods for solving Toeplitz-like systems with O(n log n) operations [, ] However, numerical solution of the Toeplitz eigenvalue problem has only recently received attention [5, 9, 15, 1] In particular, Cybenko and Van Loan [5] presented a method for using Levinson's algorithm [1] to nd the smallest eigenvalue of an n n Hermitian Toeplitz matrix with O(n ) operations In [15], Trench extended their method and gave an iterative method for computing arbitrary eigenvalues and associated eigenvectors of Hermitian Toeplitz and Toeplitz{plus{Hankel matrices at a cost of O(n ) per eigenvalue The purpose of this paper is to use Trench's method to compute the eigenvalues and eigenvectors of Hermitian Toeplitz{like matrices In xwe propose an algorithm for nding individual eigenvalues of an nn matrix with displacement rank not greater than at a computation cost of O(n ) each In x we give examples of Hermitian matrices with displacement structure (1), along with specic formulas for the associated matrices G n and H n In x we discuss an application to signal processing In x5 we describe the results of numerical experiments with the algorithm The algorithm The following theorem from [1] provides the motivation and the theoretical basis for the method Part of this theorem goes back at least to Wilkinson [1] (For the statement concerning the inertia of A n, I n, see also Browne [1]) Theorem 1 Let A n =[a ij ] n ij=1 be a Hermitian matrix, and dene Let p () =1, A m =[a ij ] m ij=1 1 m n: p m () = det (A m, I m ) 1 m n and q m () = p m() p m () 1 m n:
5 EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES Dene v m = a 1m+1 a m m n, 1: a mm+1 Let S m be the spectrum of A m and S n = [ n m=1 S m If is real let Neg n () be the number (counting multiplicities) of eigenvalues of A n less than For each = S n let w () =and w m () = be the solutions of the systems w 1m () w m () w mm () 5 1 m n, 1 () (A m, I m )w m () =v m 1 m n, 1: Dene () Then y m () = wm () m n: () (A m, I m )y m () =,q m ()e m m n where e m =[ 1] T is the last column of I m Moreover, q m () =a mm,, v m w m() 1 m n q m() =,kw m ()k and Neg n () equals the number of negative quantities in fq 1 ()q ():::q n ()g Finally, if is an eigenvalue of A n, then y n () is an associated eigenvector Theorem 1 provides a way to compute p n ()=p n () and the inertia of A n,i n Therefore, in principle it can be used in conjunction with a root{nding procedure to determine a given eigenvalue i of A n, provided that i is not \too close" to an eigenvalue of one of the principal submatrices A 1 A :::A n of A n This method is not practical for general Hermitian matrices, because in general O(n ) operations are required to solve the systems () for each value of However, Theorem 1 can be useful if A n is structured so that this computational cost is O(n ) In [15], Trench described a computational strategy combining Theorem 1 with bisection and the Pegasus root{nding method for computing individual eigenvalues and eigenvectors of Hermitian Toeplitz matrices with O(n ) operations In [1] he applied the same strategy to Hermitian Toeplitz{plus{Hankel matrices We will now show that a similar approach can be used to compute individual eigenvalues and eigenvectors of Toeplitz{ like matrices
6 Michael K Ng and William F Trench Henceforth G m and H m (1 m n) are the m matrices obtained by dropping rows m +1:::n from G n and H n thus (5) G m = U mn G n and H m = U mn H n where U mn is the m n matrix obtained by dropping the same rows from I m We denote the jth column of G m by thus, g (m) j = g 1j g mj 5 G m =[g (m) 1 g (m) g (m) ]: The following result of Heinig and Rost [8, p11] is crucial to our approach Lemma If A n satises (1) then () A m Z m, Z m A m = G m H T m, v me T m 1 m n, 1: Proof It is easily veried that U mn A n Z n U T mn = A m Z m + v m e T m and U mn Z n A n U T mn = Z m A m : Therefore we can obtain () by multiplying (1) on the left by U mn and on the right by U T mn, and invoking (5) The following algorithm provides an O(n ) method for solving the linear systems () if A n satises (1) with G n H n C n The algorithm is an adaptation of a recursion formula given in [8, p11] for solving systems with Toeplitz{like matrices Algorithm If =S n then q 1 ():::q n () can be computed as follows: q 1 () =a 11, w 1 () = a 1 q 1 () f (1) j () = g 1j q 1 () 1 j and for m n, () (8) and (9) f (m) j () = w m () = f (m) j () q m () =a mm,, v m w m() w m () y m () = wm () (g mj, v m f (m) j ()), y m () 1 j q m () h, f (m) 1 () f (m) () f (m) () i H T m y m():
7 EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES 5 (1) Proof Adding and subtracting Z m on the left side of () yields (A m, I m )Z m, Z m (A m, I m )=G m H T m, v me T m 1 m n, 1: From () and (), for m n, e T m y m() = Z m (A m, I m )y m () = and Z m ()y m () = Therefore, multiplying (1) on the right by y m () yields (A m, I m ) w m () = G m H T m y m()+v m m n: w m () Multiplying by (A m, I m ) and recalling () shows that this is equivalent to where w m () = w m () which we write in terms of its columns as These columns are the solutions of (11) Since, F m ()H T m y m() F m () =(A m, I m ) G m F m () =[f (m) 1 () f (m) () f (m) ()]: (A m, I m )f (m) j () =g (m) j = g (m) j g mj 1 j : (A m, I m )f (m) j () =g (m) j and (A m, I m )y m () =,q m ()e m it follows that the solutions of (11) are given by (8) Examples The following are examples of Hermitian matrices with the kind of displacement structure indicated in (1) (i) A Hermitian Toeplitz matrix T n =[t i,j ] n ij=1 (where t,r = t r ) has displacement rank at most, since T n Z n, Z n T n = 1 t n t 1 5 t1 t n (ii) If T n =[t i,j ] n ij=1 is an arbitrary (not necessarily Hermitian) n n Toeplitz matrix, then A n = Tn T n has displacement rank at most (see [8, p1] and [1]) It can be shown that A n satises (1) with G n = 1 t t n b 1 5 and H n = t,n+1 t 1 b n t,t n c 1 t,n+1,t 1 c n 5 : :
8 Michael K Ng and William F Trench where b i = nx k=1 t k,i+1 t k,n and c j = n X k=1 t k t k,j : P k (iii) A matrix of the form i=1 T n (i) T n (i), where T n (1) :::T n (k) are arbitrary Toeplitz matrices, has displacement rank at most k [] Matrices like this arise in solving the normal equations of Toeplitz least squares problems in signal and image processing [] (iv) A Hermitian Toeplitz{block matrix of the form (1) A n = T m (11) T m (1) T m (1s) T m (1) T m () T m (s) T m (1s) T m (1s) T m (s) T m (ss) T m (s) T m (ss) where n = sm and ft m (ij) g s ij=1 are Toeplitz matrices given by [T m (ij) ] m = kl=1 t(ij) k,l, has displacement rank at most s [8, p1] For example, if s = it can be shown that (1) holds with n =m, G n = 1 t (1) t (1) 1, t (11),m+1,t (1),m+1 t (1) m, t(11),t (1) 1 t () t () 1, t (1),m+1,t (),m+1 t () m, t(1),t () 5 and H n = 5 t (11) t (1), t (11) m t (11),m+1 t (1),m+1, t(11) 1,t (11) 1 t (1) t (), t (1) m t (1),m+1 t (),m+1, t(1) 1,t (1) 1 (v) For an example closely related to (iv) let B n be the block{toeplitz matrix B n = where each block is an s s matrix thus, h C (r) s = C (p,q) s h c (r) ij i m pq=1 i s : ij=1 Now let P n be the n n permutation matrix dened as follows: for k =1 :::m, rows (k, 1)s + 1 through ks of P n are rows k k + m:::k+(s, 1)m of I n Then A n = P n B n Pn T is the Toeplitz{block matrix where A n = h h Tm ij = T (ij) m c (p,q) ij i s ij=1 i m : pq=1 Moreover, if B n is Hermitian then so is A n that is, A n is of the form (1) Finally, if is an eigenvalue and x is an associated eigenvector of A n, then is an eigenvalue and P T n x is an associated eigenvector of B n 5 :
9 EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES An application to signal processing The input fx k g and the output fy k g of a transversal lter of order n are related by y r = n X k= w k x r,k : In signal processing problems it is often necessary to estimate the lter coecients fw w 1 :::w n g given observed values fx 1 x :::x m g and fy 1 y :::y m g of the input and output, where m>n One way to do this is to choose fw w 1 :::w n g so as to minimize (w w 1 :::w n )= mx r=1 y r, n X k= w k x r,k! where it is assumed that x j = if j An elementary argument shows that fw w 1 :::w n g should be chosen so that nx j=1 a ij w j = mx r=1 y r x r,i+1 1 i n where a ij = mx r=1 x r,i+1 x r,j+1 : The matrix A n =[a ij ] n ij=1 is given by A n = X T X, where X is the m n Toeplitz matrix (1) X = x 1 x x 1 x x 1 x n x1 x m xm,n x m x m x m,n+ x m,n+1 The matrix X T X is called the normal equations matrix or the information matrix of the corresponding least squares problem [1, 1] It is an approximation to the the correlation matrix of the input signal data We are interesting in computing the eigenvalues of X T X because, for example, the smallest and the largest eigenvalues of X T X are related to the accuracy of the least squares computations and the stability of least squares algorithms [1, 1] In [] it was shown that the lter coecients that maximize the output signal-to-noise ratio can be obtained from the eigenvector of X T X associated with its largest eigenvalue 5 :
10 8 Michael K Ng and William F Trench It can be shown that A n = X T X satises (1) with G n = 1 x m u 1 5 and H n = x m,n+ u n,x m v 1,x m,n+ v n 5 where u i = m,n+1 X l=1 x l x l+n,i and v j = mx l=j+1 x l x l,j : Therefore each iteration of Algorithm requires O(n ) operations 5 Numerical results We tried Algorithm on Toeplitz{block matrices (with s =)asmentioned in x and on matrices of the form T n T n where T n is an arbitrary real Toeplitz matrix The elements of these matrices are randomly generated with a uniform distribution in [ 1] All computations were done with Matlab in double precision Let 1 n be the eigenvalues of a Toeplitz-like matrix A n, and suppose we wish to nd i, where i is a specied integer in [1:::n] We assume that i is not an eigenvalue of any of the principal submatrices A 1 :::A n We rst nd an interval ( ) containing i but not any other eigenvalues of A n,orany eigenvalues of A n On such aninterval q n is continuous In [15] itwas shown that and satisfy this requirement if and only if Neg n () =i, 1 Neg n () =i q n () > and q n () < and a strategy was given for obtaining ( ) by means of bisection After ( ) is determined, we use the Matlab M-le \fzero" to nd i as a root of the function q n () (This root{nding algorithm was originated by T Dekker and further improved by R Brent see Matlab on-line documentation) We stop the iteration for i when the dierence between successive iterates k and k obtained by the root nder satises the inequality j k, k j 1 1 maxfj k j 1g: To check the accuracy of the individual eigenvalues and associated eigenvectors of the randomly generated Toeplitz{like matrices, we computed the residual norms i = ka ny n ( ~ i ), ~ i y n ( ~ i )k ky n ( ~ i )k where ~ i is the approximate ith eigenvalue and y n ( ~ i ) (as dened in () with = ~ i ) is an approximate i {eigenvector Tables 1 and show the distribution of f i g for 5 randomly generated matrices of order 1, 5 of order 5, and 5 of order 1, for two types of Toeplitz{like matrices Table lists the average number of iterations per eigenvalue for two types of Toeplitz{like matrices
11 EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES 9 For each randomly generated Toeplitz{like matrix of order n we formed the diagonal matrix D n consisting of the computed eigenvalues and the matrix n whose columns are the corresponding computed eigenvectors For each matrix we computed the reconstruction and orthogonality errors = ka n, n D n T n k F ka n k F and = ki n, n T n k F p n where kk F is the Frobenius norm The results are shown in Tables and 5 We also tried Algorithm on matrices of the form X T X where X is as in (1), with m = 1 and n = 18 We considered 5 cases with fx 1 :::x 1 g generated by the second-order autoregressive (AR) process x k, 1:x k +:5x k, = k and 5 cases with fx 1 :::x 1 g generated by the second-order moving-average (MA) process x k = k +:5 k +:5 k, : In all instances f k g is a Gaussian process with mean zero and variance one, and E( j k )= jk Tables and show the distribution of the residual norm i and the relative error between the eigenvalues computed by Algorithm and those computed by the QR method, respectively Table 8 shows the values of and for these two input processes The average numbers of iterations per eigenvalue for the AR and MA processes were 1 and 15 respectively Table 1 Distribution of errors f i g for 5 matrices A n = T n Tn, where Tn are randomly generated nonsymmetric n n Toeplitz matrices Number of errors Interval n = 1 n = 5 n = 1 [1, 1 ) [1, 1, ) 1 [1, 1, ) 1 [1,5 1, ) 1 1 [1, 1,5 ) 5 1 [1, 1, ) [1,8 1, ) [1,9 1,8 ) [1 1,9 ) [1 1 1 ) [1 1 1 ) 5 [1 1 )
12 1 Michael K Ng and William F Trench Table Distribution of errors f i g for 5 randomly generated Toeplitz-block matrices with s =and n =m Number of errors Interval n = 1 n = 5 n = 1 [1, 1 ) [1, 1, ) [1, 1, ) 1 [1,5 1, ) 1 15 [1, 1,5 ) [1, 1, ) 91 9 [1,8 1, ) [1,9 1,8 ) [1 1,9 ) [1 1 1 ) [1 1 1 ) [1 1 ) 1 Table Average number of iterations per eigenvalue for computations summarized in Tables 1 and Number of iterations Type n = 1 n = 5 n = 1 Tn T n where T n are nonsymmetric Toeplitz matrices Toeplitz{block matrices Table Reconstruction and orthogonality errors for 5 matrices A n = T ntn where Tn are randomly generated nonsymmetric Toeplitz matrices n = 1 n = 5 n = 1 Interval [1,5 1, ) [1, 1,5 ) 1 1 [1, 1, ) 11 1 [1,8 1, ) [1,9 1,8 ) [1 1,9 ) 8 [1 1 1 ) Table 5 Reconstruction and orthogonality errors for 5 randomly generated Toeplitz-block matrices with s =and n =m n = 1 n = 5 n = 1 Interval [1, 1, ) 1 8 [1,8 1, ) [1,9 1,8 ) [1 1,9 ) [1 1 1 ) 5 [1 1 1 ) 1
13 EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ{LIKE MATRICES 11 Table Distribution of errors f i g for 5 matrices X T X with m = 1 and n = 18 Number of errors Interval AR Process MA Process [1, 1, ) 1 [1,8 1, ) 8 1 [1,9 1,8 ) [1 1,9 ) [1 1 1 ) Table Distribution of the relative error between the eigenvalues computed by Algorithm : method and those computed by QR method for 5 matrices X T X with m = 1 and n = 18 Number of errors Interval AR Process MA Process [1, 1, ) [1,8 1, ) 1 1 [1,9 1,8 ) [1 1,9 ) 1 [1 1 1 ) 5 5 Table 8 Reconstruction and orthogonality errors for 5 matrices for X T X with m = 1 and n = 18 AR Process MA Process Interval [1,8 1, ) [1,9 1,8 ) [1 1,9 ) 5 8 [1 1 1 ) 1 1
14 1 Michael K Ng and William F Trench Summary The experimental results reported here show that Algorithm is an ecient and eective method for computing individual eigenvalues of Hermitian Toeplitz-like matrices For an n n Toeplitz{like matrix, the computational cost of each eigenvalue and an associated eigenvector is O(n ) operations The method is more ecient than general purpose methods such as the QR algorithm for obtaining a small number (compared to n) of eigenvalues (See [15]) Since the computation of each eigenvalue is independent of the computation of all others, the method is highly parallelizable Moreover, if q 1 (),, q n () are computed with a parallel processing machine utilizing as many processors as necessary to exploit the full parallelism in the algorithm, the multiplications as well as additions required to compute in (), (8) and (9) can be carried out simultaneously The inner products in (), (8) and (9) can also be computed simultaneously by employing parallel processors in O(log n) time units Therefore, the computations of fq 1 (),,q n ()g when performed by O(n) parallel processors, can be accomplished in O(n log n) time units Hence the computations of each eigenvalue, when performed by O(n) processors, can be accomplished in O(n log n) time REFERENCES [1] E T Browne, On the separation property of the roots of the secular equation, Amer J Math, 5 (19), pp 81{85 [] R Chan, J Nagy, and R Plemmons, Circulant Preconditioned Toeplitz Least Squares Iterations, SIAM J Matrix Anal Appl, 15 (199), pp 8{9 [] R Chan, J Nagy, and R Plemmons, Displacement Preconditioner for Toeplitz Least Squares Iterations, Elec Trans Numer Anal, (199), pp {5 [] R Chan and M Ng, Conjugate Gradient Methods for Solving Toeplitz Systems, SIAM Review, 8 (199), pp {8 [5] G Cybenko and C Van Loan, Computing the Minimum Eigenvalue of a Symmetric Positive Denite Toeplitz Matrix, SIAM J Sci Stat Comput, (198), pp 1{11 [] I Gohberg, T Kailath and V Olshevsky, Fast Gaussian Elimination with Partial Pivoting for Matrices with Displacement Structure, Math Comp, (1995), pp 155{15 [] S Haykin, Adaptive Filter Theory, nd ed, Prentice-Hall, Englewood Clis, NJ, 1991 [8] G Heinig and K Rost, Algebraic Methods for Toeplitz-like Matrices and Operators, Operator Theory: Advances and Applications V1, Birkhauser Verlag, Stuttgart, 198 [9] Y H Hu and S Y Kung, Toeplitz Eigensystem Solver, IEEE Trans Acoustics, Speech, Sig Proc, (1985), pp 1{11 [1] T Kailath, S Kung and M Morf, Displacement Ranks of Matrices and Linear Equations, J Math Anal Appl, 8 (199), pp 95{ [11] T Kailath and A Sayed, Displacement Structure: Theory and Applications, SIAM Review, (1995), pp 9{8 [1] N Levinson, The Wiener RMS (Root Mean Square) Error Criterion in Filter Design and Prediction, J Math and Phys, 5 (19), pp 1{8 [1] W Ferng, G Golub and R Plemmons, Adaptive Lanczos Methods for Recursive Condition Estimation, Numer Algo, 1 (1991), pp 1{19 [1] D Pierce and R Plemmons, Tracking the Condition Number for RLS in Signal Processing, Math Control Signals Systems, 5 (199), pp {9 [15] W F Trench, Numerical Solution of the Eigenvalue Problem for Hermitian Toeplitz Matrices, SIAM J Matrix Th Appl, 9 (1988), pp 91{ [1] W F Trench, Numerical Solution of the Eigenvalue Problem for Eciently Structured Hermitian Matrices, Linear Algebra Appl, 15{15 (1991), pp 15{ [1] J H Wilkinson, The Algebraic Eigenvalue Problem, Clarendon Press, 195
Preconditioning of elliptic problems by approximation in the transform domain
TR-CS-97-2 Preconditioning of elliptic problems by approximation in the transform domain Michael K. Ng July 997 Joint Computer Science Technical Report Series Department of Computer Science Faculty of
More informationNUMERICAL SOLUTION OF THE EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ MATRICES. William F. Trench* SIAM J. Matrix Anal. Appl.
NUMERICAL SOLUTION OF THE EIGENVALUE PROBLEM FOR HERMITIAN TOEPLITZ MATRICES William F. Trench* SIAM J. Matrix Anal. Appl. 10 (1989) 135-156 Abstract. An iterative procedure is proposed for computing the
More informationToeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y
Toeplitz-circulant Preconditioners for Toeplitz Systems and Their Applications to Queueing Networks with Batch Arrivals Raymond H. Chan Wai-Ki Ching y November 4, 994 Abstract The preconditioned conjugate
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More information1 Introduction The conjugate gradient (CG) method is an iterative method for solving Hermitian positive denite systems Ax = b, see for instance Golub
Circulant Preconditioned Toeplitz Least Squares Iterations Raymond H. Chan, James G. Nagy y and Robert J. Plemmons z November 1991, Revised April 1992 Abstract We consider the solution of least squares
More informationThe restarted QR-algorithm for eigenvalue computation of structured matrices
Journal of Computational and Applied Mathematics 149 (2002) 415 422 www.elsevier.com/locate/cam The restarted QR-algorithm for eigenvalue computation of structured matrices Daniela Calvetti a; 1, Sun-Mi
More informationHeinrich Voss. Dedicated to Richard S. Varga on the occasion of his 70th birthday. Abstract
A Symmetry Exploiting Lanczos Method for Symmetric Toeplitz Matrices Heinrich Voss Technical University Hamburg{Harburg, Section of Mathematics, D{27 Hamburg, Federal Republic of Germany, e-mail: voss
More informationOUTLINE 1. Introduction 1.1 Notation 1.2 Special matrices 2. Gaussian Elimination 2.1 Vector and matrix norms 2.2 Finite precision arithmetic 2.3 Fact
Computational Linear Algebra Course: (MATH: 6800, CSCI: 6800) Semester: Fall 1998 Instructors: { Joseph E. Flaherty, aherje@cs.rpi.edu { Franklin T. Luk, luk@cs.rpi.edu { Wesley Turner, turnerw@cs.rpi.edu
More informationSplit Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hankel Matrices
Split Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hanel Matrices Georg Heinig Dept of Math& CompSci Kuwait University POBox 5969 Safat 13060, Kuwait Karla Rost Faultät für Mathemati
More informationParallel Computation of the Eigenstructure of Toeplitz-plus-Hankel matrices on Multicomputers
Parallel Computation of the Eigenstructure of Toeplitz-plus-Hankel matrices on Multicomputers José M. Badía * and Antonio M. Vidal * Departamento de Sistemas Informáticos y Computación Universidad Politécnica
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationFraction-free Row Reduction of Matrices of Skew Polynomials
Fraction-free Row Reduction of Matrices of Skew Polynomials Bernhard Beckermann Laboratoire d Analyse Numérique et d Optimisation Université des Sciences et Technologies de Lille France bbecker@ano.univ-lille1.fr
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationApproximate Inverse-Free Preconditioners for Toeplitz Matrices
Approximate Inverse-Free Preconditioners for Toeplitz Matrices You-Wei Wen Wai-Ki Ching Michael K. Ng Abstract In this paper, we propose approximate inverse-free preconditioners for solving Toeplitz systems.
More informationRITZ VALUE BOUNDS THAT EXPLOIT QUASI-SPARSITY
RITZ VALUE BOUNDS THAT EXPLOIT QUASI-SPARSITY ILSE C.F. IPSEN Abstract. Absolute and relative perturbation bounds for Ritz values of complex square matrices are presented. The bounds exploit quasi-sparsity
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 informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
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 informationOld and new algorithms for Toeplitz systems
Old and new algorithms for Toeplitz systems Richard P. Brent Computer Sciences Laboratory Australian National University Canberra, Australia ABSTRACT Toeplitz linear systems and Toeplitz least squares
More informationCS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3
CS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3 Felix Kwok February 27, 2004 Written Problems 1. (Heath E3.10) Let B be an n n matrix, and assume that B is both
More informationON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH
ON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH V. FABER, J. LIESEN, AND P. TICHÝ Abstract. Numerous algorithms in numerical linear algebra are based on the reduction of a given matrix
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 informationAutomatica, 33(9): , September 1997.
A Parallel Algorithm for Principal nth Roots of Matrices C. K. Koc and M. _ Inceoglu Abstract An iterative algorithm for computing the principal nth root of a positive denite matrix is presented. The algorithm
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationNP-hardness of the stable matrix in unit interval family problem in discrete time
Systems & Control Letters 42 21 261 265 www.elsevier.com/locate/sysconle NP-hardness of the stable matrix in unit interval family problem in discrete time Alejandra Mercado, K.J. Ray Liu Electrical and
More informationSubset selection for matrices
Linear Algebra its Applications 422 (2007) 349 359 www.elsevier.com/locate/laa Subset selection for matrices F.R. de Hoog a, R.M.M. Mattheij b, a CSIRO Mathematical Information Sciences, P.O. ox 664, Canberra,
More informationA Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations
A Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations Jin Yun Yuan Plamen Y. Yalamov Abstract A method is presented to make a given matrix strictly diagonally dominant
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis Volume 8, pp 127-137, 1999 Copyright 1999, ISSN 1068-9613 ETNA etna@mcskentedu BOUNDS FOR THE MINIMUM EIGENVALUE OF A SYMMETRIC TOEPLITZ MATRIX HEINRICH VOSS
More informationG1110 & 852G1 Numerical Linear Algebra
The University of Sussex Department of Mathematics G & 85G Numerical Linear Algebra Lecture Notes Autumn Term Kerstin Hesse (w aw S w a w w (w aw H(wa = (w aw + w Figure : Geometric explanation of the
More informationAlgorithms to solve block Toeplitz systems and. least-squares problems by transforming to Cauchy-like. matrices
Algorithms to solve block Toeplitz systems and least-squares problems by transforming to Cauchy-like matrices K. Gallivan S. Thirumalai P. Van Dooren 1 Introduction Fast algorithms to factor Toeplitz matrices
More informationOutline Introduction: Problem Description Diculties Algebraic Structure: Algebraic Varieties Rank Decient Toeplitz Matrices Constructing Lower Rank St
Structured Lower Rank Approximation by Moody T. Chu (NCSU) joint with Robert E. Funderlic (NCSU) and Robert J. Plemmons (Wake Forest) March 5, 1998 Outline Introduction: Problem Description Diculties Algebraic
More informationReduced Synchronization Overhead on. December 3, Abstract. The standard formulation of the conjugate gradient algorithm involves
Lapack Working Note 56 Conjugate Gradient Algorithms with Reduced Synchronization Overhead on Distributed Memory Multiprocessors E. F. D'Azevedo y, V.L. Eijkhout z, C. H. Romine y December 3, 1999 Abstract
More informationNumerical solution of the inverse eigenvalue problem for real symmetric Toeplitz matrices
Trinity University From the SelectedWorks of William F. Trench 1997 Numerical solution of the inverse eigenvalue problem for real symmetric Toeplitz matrices William F. Trench, Trinity University Available
More informationPROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS
PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix
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 informationChapter 3 Least Squares Solution of y = A x 3.1 Introduction We turn to a problem that is dual to the overconstrained estimation problems considered s
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
More informationOn the Convergence of the Inverses of Toeplitz Matrices and Its Applications
180 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 49, NO. 1, JANUARY 2003 On the Convergence of the Inverses of Toeplitz Matrices Its Applications Feng-Wen Sun, Member, IEEE, Yimin Jiang, Member, IEEE,
More informationA Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Squares Problem
A Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Suares Problem Hongguo Xu Dedicated to Professor Erxiong Jiang on the occasion of his 7th birthday. Abstract We present
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More information0 Introduction Let the matrix A C nn be given by all its n entries The problem is to compute the product Ab of a matrix A by input vector b C n Using
Fast algorithms with preprocessing for matrix{vector multiplication problems IGohberg and VOlshevsky School of Mathematical Sciences Raymond and Beverly Sackler Faculty of Exact Sciences Tel Aviv University,
More informationNumerical Linear Algebra
Numerical Linear Algebra The two principal problems in linear algebra are: Linear system Given an n n matrix A and an n-vector b, determine x IR n such that A x = b Eigenvalue problem Given an n n matrix
More informationANONSINGULAR tridiagonal linear system of the form
Generalized Diagonal Pivoting Methods for Tridiagonal Systems without Interchanges Jennifer B. Erway, Roummel F. Marcia, and Joseph A. Tyson Abstract It has been shown that a nonsingular symmetric tridiagonal
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationDense LU factorization and its error analysis
Dense LU factorization and its error analysis Laura Grigori INRIA and LJLL, UPMC February 2016 Plan Basis of floating point arithmetic and stability analysis Notation, results, proofs taken from [N.J.Higham,
More informationarxiv: v1 [math.ra] 13 Jan 2009
A CONCISE PROOF OF KRUSKAL S THEOREM ON TENSOR DECOMPOSITION arxiv:0901.1796v1 [math.ra] 13 Jan 2009 JOHN A. RHODES Abstract. A theorem of J. Kruskal from 1977, motivated by a latent-class statistical
More information14.2 QR Factorization with Column Pivoting
page 531 Chapter 14 Special Topics Background Material Needed Vector and Matrix Norms (Section 25) Rounding Errors in Basic Floating Point Operations (Section 33 37) Forward Elimination and Back Substitution
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More information12 CHAPTER 1. PRELIMINARIES Lemma 1.3 (Cauchy-Schwarz inequality) Let (; ) be an inner product in < n. Then for all x; y 2 < n we have j(x; y)j (x; x)
1.4. INNER PRODUCTS,VECTOR NORMS, AND MATRIX NORMS 11 The estimate ^ is unbiased, but E(^ 2 ) = n?1 n 2 and is thus biased. An unbiased estimate is ^ 2 = 1 (x i? ^) 2 : n? 1 In x?? we show that the linear
More informationInverses of regular Hessenberg matrices
Proceedings of the 10th International Conference on Computational and Mathematical Methods in Science and Engineering, CMMSE 2010 27 30 June 2010. Inverses of regular Hessenberg matrices J. Abderramán
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationNORTHERN ILLINOIS UNIVERSITY
ABSTRACT Name: Santosh Kumar Mohanty Department: Mathematical Sciences Title: Ecient Algorithms for Eigenspace Decompositions of Toeplitz Matrices Major: Mathematical Sciences Degree: Doctor of Philosophy
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationJordan Journal of Mathematics and Statistics (JJMS) 5(3), 2012, pp A NEW ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS OF EQUATIONS
Jordan Journal of Mathematics and Statistics JJMS) 53), 2012, pp.169-184 A NEW ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS OF EQUATIONS ADEL H. AL-RABTAH Abstract. The Jacobi and Gauss-Seidel iterative
More informationPart I: Preliminary Results. Pak K. Chan, Martine Schlag and Jason Zien. Computer Engineering Board of Studies. University of California, Santa Cruz
Spectral K-Way Ratio-Cut Partitioning Part I: Preliminary Results Pak K. Chan, Martine Schlag and Jason Zien Computer Engineering Board of Studies University of California, Santa Cruz May, 99 Abstract
More informationConvexity of the Joint Numerical Range
Convexity of the Joint Numerical Range Chi-Kwong Li and Yiu-Tung Poon October 26, 2004 Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement. Abstract Let A = (A 1,..., A m ) be an
More informationThe geometric mean algorithm
The geometric mean algorithm Rui Ralha Centro de Matemática Universidade do Minho 4710-057 Braga, Portugal email: r ralha@math.uminho.pt Abstract Bisection (of a real interval) is a well known algorithm
More informationOn the stability of the Bareiss and related Toeplitz factorization algorithms
On the stability of the Bareiss and related oeplitz factorization algorithms A. W. Bojanczyk School of Electrical Engineering Cornell University Ithaca, NY 14853-5401, USA adamb@toeplitz.ee.cornell.edu
More informationDirect Methods for Solving Linear Systems. Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le
Direct Methods for Solving Linear Systems Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le 1 Overview General Linear Systems Gaussian Elimination Triangular Systems The LU Factorization
More informationIntroduction to Linear Algebra. Tyrone L. Vincent
Introduction to Linear Algebra Tyrone L. Vincent Engineering Division, Colorado School of Mines, Golden, CO E-mail address: tvincent@mines.edu URL: http://egweb.mines.edu/~tvincent Contents Chapter. Revew
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More information1. Introduction Let the least value of an objective function F (x), x2r n, be required, where F (x) can be calculated for any vector of variables x2r
DAMTP 2002/NA08 Least Frobenius norm updating of quadratic models that satisfy interpolation conditions 1 M.J.D. Powell Abstract: Quadratic models of objective functions are highly useful in many optimization
More informationUMIACS-TR July CS-TR 2494 Revised January An Updating Algorithm for. Subspace Tracking. G. W. Stewart. abstract
UMIACS-TR-9-86 July 199 CS-TR 2494 Revised January 1991 An Updating Algorithm for Subspace Tracking G. W. Stewart abstract In certain signal processing applications it is required to compute the null space
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 informationLecture 5 Singular value decomposition
Lecture 5 Singular value decomposition Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationFast Transforms: Banded Matrices with Banded Inverses
Fast Transforms: Banded Matrices with Banded Inverses 1. Introduction Gilbert Strang, MIT An invertible transform y = Ax expresses the vector x in a new basis. The inverse transform x = A 1 y reconstructs
More informationJae Heon Yun and Yu Du Han
Bull. Korean Math. Soc. 39 (2002), No. 3, pp. 495 509 MODIFIED INCOMPLETE CHOLESKY FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC POSITIVE DEFINITE MATRIX Jae Heon Yun and Yu Du Han Abstract. We propose
More informationInterlacing Inequalities for Totally Nonnegative Matrices
Interlacing Inequalities for Totally Nonnegative Matrices Chi-Kwong Li and Roy Mathias October 26, 2004 Dedicated to Professor T. Ando on the occasion of his 70th birthday. Abstract Suppose λ 1 λ n 0 are
More information1182 L. B. Beasley, S. Z. Song, ands. G. Lee matrix all of whose entries are 1 and =fe ij j1 i m 1 j ng denote the set of cells. The zero-term rank [5
J. Korean Math. Soc. 36 (1999), No. 6, pp. 1181{1190 LINEAR OPERATORS THAT PRESERVE ZERO-TERM RANK OF BOOLEAN MATRICES LeRoy. B. Beasley, Seok-Zun Song, and Sang-Gu Lee Abstract. Zero-term rank of a matrix
More informationA Note on Eigenvalues of Perturbed Hermitian Matrices
A Note on Eigenvalues of Perturbed Hermitian Matrices Chi-Kwong Li Ren-Cang Li July 2004 Let ( H1 E A = E H 2 Abstract and à = ( H1 H 2 be Hermitian matrices with eigenvalues λ 1 λ k and λ 1 λ k, respectively.
More informationPerturbation results for nearly uncoupled Markov. chains with applications to iterative methods. Jesse L. Barlow. December 9, 1992.
Perturbation results for nearly uncoupled Markov chains with applications to iterative methods Jesse L. Barlow December 9, 992 Abstract The standard perturbation theory for linear equations states that
More informationAN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES
AN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES ZHONGYUN LIU AND HEIKE FAßBENDER Abstract: A partially described inverse eigenvalue problem
More informationA PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
IJMMS 25:6 2001) 397 409 PII. S0161171201002290 http://ijmms.hindawi.com Hindawi Publishing Corp. A PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
More informationDirect Methods for Solving Linear Systems. Matrix Factorization
Direct Methods for Solving Linear Systems Matrix Factorization Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationExponential Decomposition and Hankel Matrix
Exponential Decomposition and Hankel Matrix Franklin T Luk Department of Computer Science and Engineering, Chinese University of Hong Kong, Shatin, NT, Hong Kong luk@csecuhkeduhk Sanzheng Qiao Department
More informationthen kaxk 1 = j a ij x j j ja ij jjx j j: Changing the order of summation, we can separate the summands, kaxk 1 ja ij jjx j j: let then c = max 1jn ja
Homework Haimanot Kassa, Jeremy Morris & Isaac Ben Jeppsen October 7, 004 Exercise 1 : We can say that kxk = kx y + yk And likewise So we get kxk kx yk + kyk kxk kyk kx yk kyk = ky x + xk kyk ky xk + kxk
More informationTotal least squares. Gérard MEURANT. October, 2008
Total least squares Gérard MEURANT October, 2008 1 Introduction to total least squares 2 Approximation of the TLS secular equation 3 Numerical experiments Introduction to total least squares In least squares
More informationPerformance Comparison of Two Implementations of the Leaky. LMS Adaptive Filter. Scott C. Douglas. University of Utah. Salt Lake City, Utah 84112
Performance Comparison of Two Implementations of the Leaky LMS Adaptive Filter Scott C. Douglas Department of Electrical Engineering University of Utah Salt Lake City, Utah 8411 Abstract{ The leaky LMS
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More informationIterative Algorithm for Computing the Eigenvalues
Iterative Algorithm for Computing the Eigenvalues LILJANA FERBAR Faculty of Economics University of Ljubljana Kardeljeva pl. 17, 1000 Ljubljana SLOVENIA Abstract: - We consider the eigenvalue problem Hx
More informationProblem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class).
Math 3, Fall Jerry L. Kazdan Problem Set 9 Due In class Tuesday, Nov. 7 Late papers will be accepted until on Thursday (at the beginning of class).. Suppose that is an eigenvalue of an n n matrix A and
More informationIs Every Matrix Similar to a Toeplitz Matrix?
Is Every Matrix Similar to a Toeplitz Matrix? D. Steven Mackey Niloufer Mackey Srdjan Petrovic 21 February 1999 Submitted to Linear Algebra & its Applications Abstract We show that every n n complex nonderogatory
More informationA TOUR OF LINEAR ALGEBRA FOR JDEP 384H
A TOUR OF LINEAR ALGEBRA FOR JDEP 384H Contents Solving Systems 1 Matrix Arithmetic 3 The Basic Rules of Matrix Arithmetic 4 Norms and Dot Products 5 Norms 5 Dot Products 6 Linear Programming 7 Eigenvectors
More information6.4 Krylov Subspaces and Conjugate Gradients
6.4 Krylov Subspaces and Conjugate Gradients Our original equation is Ax = b. The preconditioned equation is P Ax = P b. When we write P, we never intend that an inverse will be explicitly computed. P
More informationDepartment of. Computer Science. Functional Implementations of. Eigensolver. December 15, Colorado State University
Department of Computer Science Analysis of Non-Strict Functional Implementations of the Dongarra-Sorensen Eigensolver S. Sur and W. Bohm Technical Report CS-9- December, 99 Colorado State University Analysis
More informationRank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col
Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an
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 informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS. 1. Linear Equations and Matrices
ELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS KOLMAN & HILL NOTES BY OTTO MUTZBAUER 11 Systems of Linear Equations 1 Linear Equations and Matrices Numbers in our context are either real numbers or complex
More informationarxiv: v1 [math.ra] 11 Aug 2014
Double B-tensors and quasi-double B-tensors Chaoqian Li, Yaotang Li arxiv:1408.2299v1 [math.ra] 11 Aug 2014 a School of Mathematics and Statistics, Yunnan University, Kunming, Yunnan, P. R. China 650091
More informationEfficient and Accurate Rectangular Window Subspace Tracking
Efficient and Accurate Rectangular Window Subspace Tracking Timothy M. Toolan and Donald W. Tufts Dept. of Electrical Engineering, University of Rhode Island, Kingston, RI 88 USA toolan@ele.uri.edu, tufts@ele.uri.edu
More informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More information