Title: Orthogonal basis functions in discrete least-squares rational approximation
Authors: Bultheel, Adhemar ×
Van Barel, Marc
Van gucht, Patrick #
Issue Date: Mar-2004
Publisher: Elsevier
Series Title: Journal of Computational and Applied Mathematics vol:164-165 pages:175-194
Abstract: We consider a problem that arises in the field of frequency domain system identification. If a discrete-time system has an input-output relation Y(z) = G(z)U(z), with transfer function G, then the problem is to find a rational approximation G^_n for G. The data given are measurements of input and output spectra in the frequency points z_k: {U(z_k), Y(z_k): k=1...N} together with some weight. The approximation criterion is to minimize the weighted discrete least squares norm of the vector obtained by evaluating G - G^_n in the measurement points.
If the poles of the system are fixed, then the problem reduces to a linear least-squares problem in two possible ways: by multiplying out the denominators and hide these in the weight, which leads to the construction of orthogonal vector polynomials, or the problem can be solved directly using an orthogonal basis of rational functions. The orthogonality of the basis is important because if the transfer function G^_n is represented with respect to a nonorthogonal basis, then this least-squares problem can be very ill conditioned. Even if an orthogonal basis is used, but with respect to the wrong inner product (e.g., the Lebesgue measure on the unit circle) numerical instability can be fatal in practice.

We show that both approaches lead to an inverse eigenvalue problem, which forms the common framework in which fast and numerically stable algorithms can be designed for the computation of the orthonormal basis.
ISSN: 0377-0427
Publication status: published
KU Leuven publication type: IT
Appears in Collections:NUMA, Numerical Analysis and Applied Mathematics Section
× corresponding author
# (joint) last author

Files in This Item:
File Status SizeFormat
iccam02.pdf Published 291KbAdobe PDFView/Open


All items in Lirias are protected by copyright, with all rights reserved.

© Web of science