Abstract
We apply a fast adaptive condition estimation scheme, calledACE, to recursive least squares (RLS) computations in signal processing.ACE is fast in the sense that onlyO(n) operations are required forn parameter problems, and is adaptive over time, i.e., estimates at timet are used to produce estimates at timet + 1. RLS algorithms for linear prediction of time series are applied in various fields of signal processing: identification, estimation, and control. However, RLS algorithms are known to suffer from numerical instability problems under finite word-length conditions, due to ill-conditioning. We apply adaptive procedures, linear in the order of the problem, for accurately tracking relevant extreme eigen-values or singular values and the associated condition numbers over timet. In this paper exponentially weighted data windows are considered. The sliding data window case, which involves downdating as well as updating, is considered else-where. Numerical experiments indicate thatACE yields an accurate, yet inexpensive, RLS condition estimator for signal processing applications.
Similar content being viewed by others
References
S. T. Alexander,Adaptive Signal Processing, Springer-Verlag, New York, 1986.
S. T. Alexander, C.-T. Pan, and R. J. Plemmons, Analysis of a recursive least squares hyperbolic rotation scheme for signal processing,Linear Algebra Appl.,98 (1988), 3–40.
K. J., Åström and B. Wittenmark,Adaptive Control, Addison-Wesley, Reading, MA, 1989.
G. J. Bierman,Factorization Methods for Discrete Sequential Estimation, Academic Press, New York, 1977.
C. H. Bischof, Incremental condition estimation,SIAM J. Matrix Anal. Appl.,11 (1990), 312–322.
C. H. Bischof and G. M. Shroff, On updating signal subspaces, Technical Report, Argonne National Lab., 1989;IEEE Trans. Signal Process (to appear).
A. Björck, Least squares methods, inHandbook of Numerical Analysis. Volume 1:Solution of Equations (P. G. Ciarlet and J. L. Lions, eds.), Elsevier/North Holland, Amsterdam, 1990.
A. Björck, Error analysis of least squares algorithms, inNumerical Linear Algebra, Digital Signal Processing and Parallel Algorithms (G. Golub and P. Van Dooren, eds.), pp. 41–73, NATO ASI Series F, vol. 70, Springer-Verlag, Heidelberg, 1991.
P. Bolzern, A multiple forgetting algorithm for robust adaptive estimation,Internat. J. Adaptive Control Signal Process.,5 (1991) (to appear).
J. M. Cioffi, Limited precision effects in adaptive filtering,IEEE Trans. Circuits and Systems,34 (1987), 821–833.
J. M. Cioffi and T. Kailath, Windowed fast transversal filters, adaptive algorithms with normalization,IEEE Trans. Acoust. Speech Signal Process.,33 (1985), 607–625.
P. Comon and G. H. Golub, Tracking a few extreme singular values and vectors in signal processing, Stanford Univ. Technical Report NA-89-01,1989.
R. Ferng, G. H. Golub, and R. J. Plemmons, Adaptive Lanczos methods for recursive condition estimation,Numerical Algorithms,1 (1991), 1–20.
P. E. Gill, G. H. Golub, W. Murray, and M. A. Saunders, Methods for modifying matrix factorizations,Math. Comp.,28 (1974), 505–535.
G. H. Golub and C. Van Loan,Matrix Computations, 2nd edn., Johns Hopkins University Press, Baltimore, MD, 1989.
S. Haykin,Adaptive Filter Theory, Prentice-Hall, Englewood Cliffs, NJ, 1986.
C. S. Henkel, M. T. Heath, and R. J. Plemmons, Chokesky downdating on a hypercube, inHypercube Concurrent Computers and Applications (G. Fox, ed.), pp. 1592–1598, ACM Press, 1988.
C. S. Henkel and R. J. Plemmons, Recursive least squares computations on a hypercube multiprocessor using the covariance factorization,SIAM J. Sci. Statist. Comput.,12 (1990), 95–106.
N. J. Higham, A survey of condition estimators for triangular matrices,SIAM Rev.,29 (1987), 575–596.
Y. H. Hu, Adaptive methods for real time Pisarenko spectrum estimation,Proceedings of the International Conference on Acoustics, Speech, and Signal Processing, Tampa, FL, 1985, pp. 105–109.
F. Milinazzo, C. Zala, and I. Barrodale, On the rate of growth of condition numbers for convolution matrices,IEEE Trans. Acoust. Speech Signal Process.,35 (1987), 471–475.
C. Moler, J. Little, and S. Bangert,Pro-Matlab Users Guide, The Mathworks, Sherborn, MA, 1987.
M. Moonen, P. Van Dooren, and J. Vanderwalle, An SVD updating algorithm for subspace tracking, preprint, 1989.
M. Morf and T. Kailath, Square root algorithms and least squares estimation,IEEE Trans. Automat. Control,20 (1975) 487–497.
C.-T. Pan and R. J. Plemmons, Least squares modifications with inverse factorizations: parallel implications,Comp. Appl. Math.,27 (1989), 109–127.
D. J. Pierce and R. J. Plemmons, Fast adaptive condition estimation,SIAM J. Matrix Anal. Appl.,13 (1992) (to appear).
R. J. Plemmons, Recursive least squares computations, inSignal Processing, Scattering and Operator Theory, and Numerical Methods M. A. Kaashoek, J. H. van Schuppen, and A. C. M. Ran, eds.), pp. 495–502, Birkhäuser, Cambridge, MA, 1990.
S. Sastry and M. Bodson,Adaptive Control: Stability, Convergence, and Robustness, Prentice-Hall, Englewood Cliffs, NJ, 1989.
G. M. Shroff and C. H. Bischof, Adaptive condition estimation for rank-one updates of QR factorizations, Technical Report, Argonne National Lab., 1989;SIAM J. Matrix Anal. Appl. (to appear).
M. Verhaegen and P. Van Dooren, Numerical aspects of different Kalman filter implementations,IEEE Trans. Automat. Control,31 (1986), 907–917.
Author information
Authors and Affiliations
Additional information
Research supported by the Air Force under Grant No. AFOSR-88-0285 and by the National Science Foundation under Grant No. DMS-89-02121.
Rights and permissions
About this article
Cite this article
Pierce, D.J., Plemmons, R.J. Tracking the condition number for RLS in signal processing. Math. Control Signal Systems 5, 23–39 (1992). https://doi.org/10.1007/BF01211974
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01211974