Abstract
We propose a new class ofhyperbolic Gram-Schmidt methods to simultaneously update and downdate the Cholesky factor of a sample covariance matrix efficiently with applications to sliding window recursive least squares (RLS) filtering problems. Several vectorized versions of this Gram-Schmidt approach are introduced, which include conventional column-updating, modified row/column-updating, and square-root-free methods. Comparisons to the existing known methods, such as Householder transformation and Givens rotation, are also given. Upon further reformulating these algorithms, a systolic triarray structure is proposed to facilitate VLSI implementations.
Similar content being viewed by others
References
M.G. Bellanger,Adaptive digital filters and signal analysis, New York and Basel: Marcel Dekker 1987.
J.M. Cioffi, “The fast adaptive ROTOR's RLS algorithm,”IEEE Trans. Acoust., Speech, Signal Processing, vol. 38, Apr. 1990, pp. 631–653.
R.T. Compton Jr.,Adaptive antennas: Concepts and performance, Englewood Cliffs, NJ: Prentice Hall, 1988.
S. Haykin,Adaptive filter theory, Englewood Cliffs, NJ: Prentice-Hall, 1986.
S.F. Hsieh and K. Yao, “Hyperbolic Gram-Schmidt pseudo orthogonalization with applications to sliding window RLS filtering,”24-th Annual Conference on Information Science and System, Princeton: Princeton University, 1990.
S. Kalson and K. Yao, “Systolic array procesing for order and time recursive generalized least-squares estimation,”Proc. SPIE 564, Real-Time Signal Processing VIII, 1985, pp. 28–38. A more detailed version appeared as S. Kalson and K. Yao, “A Class of Least-Squares Filtering and Identification Algorithms with Systolic Array Architectures,”IEEE Trans. on Information Theory, vol. 37, pp. 43–52, 1991.
K.J.R. Liu, S.F. Hsieh, and K. Yao, “Two-level pipelined implementation of systolic block Householder transformations with application to RLS algorithm,”Proc. Int'l Conf. on Application-Specific Array processors, Princeton, 1990, pp. 758–769.
F. Ling, D. Manolakis, and J.G. Proakis, “A recursive modified Gram-Schmidt algorithm for least-squares estimation,”IEEE Trans. on Acous., Speech, and Signal Processing, vol. ASSP-34, 1986, pp. 829–836.
J.G. McWhirter, “Recursive least-squares minimization using a systolic array,”Proc. SPIE 431, Real-time signal processing VI, 1983, pp. 105–112.
C.M. Rader and A.O. Steinhardt, “Hyperbolic Householder transformations,”IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP-34, 1986, pp. 1589–1602.
R. Schreiber, “Implementation of adaptive array algorithms,”IEEE Trans. on Acoust., Speech, Signal Processing, vol. ASSP-34, 1986, pp. 338–346.
P. Strobach,Linear prediction theory: A mathematical basis for adaptive systems, New York: Springer-Verlag, 1990.
M.L. Honig and D.G. Messerschmitt,Adaptive filters, Boston: Kluwer Academic Publishers, 1984.
S.T. Alexander, C.T. Pan, and R.J. Plemmons, “Numerical properties of a hyperbolic rotation method for windowed RLS filtering,”IEEE ICASSP, 1987, pp. 423–426.
S.F. Hsieh and K. Yao, “Systolic implementation of windowed recursive LS estimation,”Proc. of IEEE Int'l Symp. on CAS, 1990, pp. 1931–1934.
G.H. Golub and C.F. Van Loan,Matrix computations, 2nd ed., Baltimore, MD, Johns Hopkins Press, 1989.
W. Murray, P.E. Gill, G.H. Golub, and M.A. Saunders, “Methods for modifying matrix factorizations,”Mathematics of Computation, vol. 28, 1974, pp. 505–535.
N.-K. Tsao, “A note on implementing the Householder transformation,”SIAM J. Numer. Anal., vol. 12, 1975, pp. 53–58.
Å. Björck, “Solving least squares problems by Gram-Schmidt orthogonalization,”BIT 7, 1967, pp. 1–21.
Å. Björck, “Least Squares Methods” in Handbook of Numerical Analysis Vol. II: Finite difference methods—Solution of equations in Rn, North Holland: Elsevier, 1989.
J.R. Rice, “Experiments on Gram-Schmidt orthogonalization,”Math. Comp. 20, 1966, pp. 325–328.
W. Hoffmann, “Iterative algorithms for Gram-Schimdt orthogonalization,”Computing, vol. 41, 1989, pp. 335–348.
W.M. Gentleman, “Least squares computations by Givens transformations without square roots,”J. Inst. Maths Applics, vol. 12, 1973, pp. 329–336.
S.F. Hsieh, K.J.R. Liu, and K. Yao, “A unified sqrt-free rank-1 up/down-dating approach for recursive least-squares problems,” to be presented inIEEE Int'l Conf. on ASSP, Toronto, 1991.
H.T. Kung, “Why systolic architectures?”IEEE Computer, 1982.
S.Y. Kung,VLSI array processors, Englewood Cliffs, NJ: Prentice-Hall, 1988.
W.M. Gentelman and H.T. Kung, “Matrix triangularization by systolic array,”Proc. SPIE, vol. 298:Real-time signal processing IV, 1981, pp. 19–26.
C.L. Lawson and R.J. Hanson,Solving least squares problems, Englewood Cliffs, NJ: Prentice-Hall, 1974.
Author information
Authors and Affiliations
Additional information
This work is partially supported by a UC MICRO grant and the NSF grant NCR-8814407. It is also partially supported by NSF grant ECD-8803012-06.
Rights and permissions
About this article
Cite this article
Hsieh, S.F., Liu, K.J.R. & Yao, K. Systolic implementations of up/down-dating cholesky factorization using vectorized Gram-Schmidt pseudo orthoganalization. J VLSI Sign Process Syst Sign Image Video Technol 3, 151–161 (1991). https://doi.org/10.1007/BF00925826
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/BF00925826