Abstract
The Structured Total Least Squares (STLS) problem is a natural extension of the Total Least Squares (TLS) problem when constraints on the matrix structure need to be imposed. Similar to the ordinary TLS approach, the STLS approach can be used to determine the parameter vector of a linear model, given some noisy measurements. In many signal processing applications, the imposition of this matrix structure constraint is necessary for obtaining Maximum Likelihood (ML) estimates of the parameter vector. In this paper we consider the Toeplitz (Hankel) STLS problem (i.e., an STLS problem in which the Toeplitz (Hankel) structure needs to be preserved). A fast implementation of an algorithm for solving this frequently occurring STLS problem is proposed. The increased efficiency is obtained by exploiting the low displacement rank of the involved matrices and the sparsity of the associated generators.
The fast implementation is compared to two other implementations of algorithms for solving the Toeplitz (Hankel) STLS problem. The comparison is carried out on a recently proposed speech compression scheme. The numerical results confirm the high efficiency of the newly proposed fast implementation: the straightforward implementations have a complexity of O((m+n)3) and O(m3) whereas the proposed implementation has a complexity of O(mn+n2).
Similar content being viewed by others
References
T.J. Abatzoglou and R.H. Hackman, Constrained total least squares based extraction of acoustic resonances of elastic targets, in: Proc. of IEEE Internat. Conf. on Acoustics, Speech and Signal Processing, Toronto, ON, Canada (May 1991) pp. 3385-3388.
T.J. Abatzoglou, J.M. Mendel and G.A. Harada, The constrained total least squares technique and its applications to harmonic superresolution, IEEE Trans. Signal Processing 39 (1991) 1070-1086.
E. Anderson, Z. Bai, C. Bischof, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, S. Ostrouchov and D. Sorensen, LAPACK Users' Guide (SIAM, Philadelphia, 1995).
E. Anderson, Z. Bai and J.J. Dongarra, Generalized QR factorization and its applications, Technical Report CS-91-131 (LAPACK Working Note 31), Computer Science Department, University of Tennessee, Knoxville (1991).
A.W. Bojanczyk, R.P. Brent, P. Van Dooren and F.R. De Hoog, A note on downdating the Cholesky factorization, SIAM J. Sci. Statist. Comput. 1 (1980) 210-220.
S. Chandrasekaran and A.H. Sayed, Stabilizing the generalized Schur algorithm.
R. Fletcher, Practical Methods of Optimization (Wiley, New York, 1987).
P.E. Gill, W. Murray and M.H. Wright, Practical Optimization (Academic Press, London, 1981).
G.H. Golub and C.F. Van Loan, An analysis of the total least squares problem, SIAM J. Numer. Anal. (1980).
G.H. Golub and C.F. Van Loan, Matrix Computations (Johns Hopkins Univ. Press, Baltimore, MD, 1989).
T. Kailath, Displacement structure and array algorithms, in: Fast Reliable Algorithms for Matrices with Structure, eds. T. Kailath and A.H. Sayed (SIAM, Philadelphia, PA, 1999).
T. Kailath and A.H. Sayed, Displacement structure: Theory and applications, SIAM Rev. 37 (1995) 297-386.
P. Lemmerling, Structured total least squares: analysis, algorithms and applications, Ph.D. thesis, Department of Electrical Engineering, ESAT-SISTA/COSIC, Katholieke Universiteit Leuven (1999).
P. Lemmerling, I. Dologlou and S. Van Huffel, Variable rate speech compression based on exact modeling and waveform vector quantization, in: Signal Processing Sympos. (SPS 98), Katholieke Universiteit Leuven, Leuven (March 1998), IEEE Benelux Signal Processing Chapter, pp. 127-130.
P. Lemmerling and S. Van Huffel, Analysis of the structured total least squares problem for Hankel/Toeplitz matrices, Technical Report/Internal Report 98-100, ESAT-SISTA/COSIC, Katholieke Universiteit Leuven (1998).
L. Ljung, System Identification: Theory for the User (Prentice-Hall, Englewood Cliffs, NJ, 1987).
N. Mastronardi, P. Van Dooren and S. Van Huffel, Stability of the generalized Schur algorithm, Technical Report/Internal Report 99-82, ESAT-SISTA/COSIC, Katholieke Universiteit Leuven (1999).
J.B. Rosen, H. Park and J. Glick, Total least norm formulation and solution for structured problems, SIAM J. Matrix Anal. Appl. 17(1) (1996) 110-128.
M. Stewart and P. Van Dooren, Stability issues in the factorization of structured matrices, SIAM J. Matrix Anal. Appl. 18 (1997) 104-118.
S. Van Huffel, C. Decanniere, H. Chen and P. Van Hecke, Algorithm for time-domain nmr data fitting based on total least squares, J. Magn. Reson. A 110 (1994) 228-237.
S. Van Huffel and J. Vandewalle, The Total Least Squares Problem: Computational Aspects and Analysis, Vol. 9 (SIAM, Philadelphia, PA, 1991).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Lemmerling, P., Mastronardi, N. & Van Huffel, S. Fast algorithm for solving the Hankel/Toeplitz Structured Total Least Squares problem. Numerical Algorithms 23, 371–392 (2000). https://doi.org/10.1023/A:1019116520737
Issue Date:
DOI: https://doi.org/10.1023/A:1019116520737