VLSI algorithms and architectures for complex householder transformation with applications to array processing

The Householder transformation is considered to be desirable among various unitary transformations due to its superior computational efficiency and robust numerical stability. Specifically, the Householder transformation outperforms the Givens rotation and the modified Gram-Schmidt methods in numerical stability under finite-precision implementations, as well as requiring fewer arithmetical operations. Consequently, the QR decomposition based on the Householder transformation is promising for VLSI implementation and real-time high throughput modern signal processing. In this paper, a recursive complex Householder transformation (CHT) with a fast initialization algorithm is proposed and its associated parallel/pipelined architecture is also considered. Then, a CHT based recursive least-squares algorithm with a fast initialization is presented. Its associated systolic array processing architecture is also considered.

[1]  Seth Z. Kalson,et al.  A class of least-squares filtering and identification algorithms with systolic array architectures , 1991, IEEE Trans. Inf. Theory.

[2]  John M. Cioffi.  The fast Householder filters-RLS adaptive filter , 1990, International Conference on Acoustics, Speech, and Signal Processing.

[3]  Che-Fen Tom Tang.  Adaptive Array Systems Using QR-Based RLS and CRLS Techniques with Systolic Array Architectures , 1991 .

[4]  L. Johnsson.  A Computational Array for the QR-Method , 1982 .

[5]  J. G. McWhirter.  Algorithmic Engineering-An Emerging Discipline , 1989, Optics & Photonics.

[6]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[7]  Kui Liu,et al.  On systolic arrays for recursive complex Householder transformations with applications to array processing , 1991, [Proceedings] ICASSP 91: 1991 International Conference on Acoustics, Speech, and Signal Processing.

[8]  Charles M. Rader,et al.  Hyperbolic householder transformations , 1986, IEEE Trans. Acoust. Speech Signal Process..

[9]  S. M. Yuen.  Algorithmic, architectural, and beam pattern issues of sidelobe cancellation , 1989 .

[10]  James W. Longley.  Least squares computations using orthogonalization methods , 1984 .

[11]  Nai–Kuan Tsao.  A Note on Implementing the Householder Transformation , 1975 .

[12]  J. G. McWhirter,et al.  Application of a systolic array to adaptive beamforming , 1984 .

[13]  J. G. McWhirter,et al.  A novel algorithm and architecture for adaptive digital beamforming , 1986 .

[14]  J. Greg Nash,et al.  Modified Faddeeva Algorithm for Concurrent Execution of Linear Algebraic Operations , 1988, IEEE Trans. Computers.

[15]  Gene H. Golub,et al.  Matrix computations , 1983 .

[16]  A. Steinhardt,et al.  Householder transforms in signal processing , 1988, IEEE ASSP Magazine.

[17]  Robert Schreiber,et al.  Implementation of adaptive array algorithms , 1986, IEEE Trans. Acoust. Speech Signal Process..

[18]  Robert Schreiber,et al.  On systolic arrays for updating the Cholesky factorization , 1986, BIT Comput. Sci. Sect..

[19]  T. J. Shepherd,et al.  A systolic array processor for MVDR beamforming , 1989 .

[20]  K.J.R. Liu,et al.  Recursive LS filtering using block Householder transformations , 1990, International Conference on Acoustics, Speech, and Signal Processing.

[21]  J. G. McWhirter,et al.  Recursive Least-Squares Minimization Using A Systolic Array , 1983, Optics & Photonics.

[22]  Fuyun Ling,et al.  A recursive modified Gram-Schmidt algorithm for least- squares estimation , 1986, IEEE Trans. Acoust. Speech Signal Process..

[23]  K. Yao,et al.  Systolic array processing for order and time recursive generalized least-squares estimation , 1986, Optics & Photonics.

[24]  J. E. Hudson,et al.  Parallel Weight Extraction By A Systolic Least Squares Algorithm , 1989, Optics & Photonics.

[25]  H. T. Kung,et al.  Matrix Triangularization By Systolic Arrays , 1982, Optics & Photonics.

[26]  K. J. Ray Liu,et al.  Systolic block Householder transformation for RLS algorithm with two-level pipelined implementation , 1992, IEEE Trans. Signal Process..