Optimally regularized inverse of singular value decomposition and application to signal extrapolation

Abstract The singular value decomposition (SVD) is an excellent tool to attain reliable and stable least squares parameters estimation in ill-conditioned cases, by discarding small singular values adequately. However, large singular values and small singular values do not always separate neatly. In the present paper, we introduce multiple regularization parameters to modify the Moore-Penrose pseudo-inverse matrix for the purpose of stabilization of ill-posed least squares problems. Optimal values of the regularization parameters can be determined so as to minimize an estimated mean squares error (EMSE) criterion calculated by using only accessible signal data. Thus, the proposed scheme can successfully give threshold conditions whether smaller singular values should be adopted or discarded. The relationship with the optimal truncation of the singular values is also investigated analytically. Effectiveness of the proposed method is discussed in applications to optimal extrapolation of band-limited signals.

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

[2]  D. Slepian Prolate spheroidal wave functions, fourier analysis, and uncertainty — V: the discrete case , 1978, The Bell System Technical Journal.

[3]  D. Youla,et al.  Image Restoration by the Method of Convex Projections: Part 1ߞTheory , 1982, IEEE Transactions on Medical Imaging.

[4]  A. Papoulis A new algorithm in spectral analysis and band-limited extrapolation. , 1975 .

[5]  Bede Liu,et al.  On the use of singular value decomposition and decimation in discrete-time band-limited signal extrapolation , 1984 .

[6]  J. Cadzow,et al.  An extrapolation procedure for band-limited signals , 1979 .

[7]  W. Hemmerle An Explicit Solution for Generalized Ridge Regression , 1975 .

[8]  D. Luenberger Optimization by Vector Space Methods , 1968 .

[9]  Louis L. Scharf,et al.  Reduced rank methods for solving least squares problems , 1989, Twenty-Third Asilomar Conference on Signals, Systems and Computers, 1989..

[10]  H. Tsuji,et al.  Simultaneous optimization method of regularization and singular value decomposition in least squares parameter identification , 1989, International Conference on Acoustics, Speech, and Signal Processing,.

[11]  Hrishikesh D. Vinod,et al.  Recent Advances in Regression Methods. , 1983 .

[12]  Louis L. Scharf,et al.  Rank reduction for modeling stationary signals , 1987, IEEE Trans. Acoust. Speech Signal Process..

[13]  Harbans S. Dhadwal,et al.  Regularized Iterative and Non-iterative Procedures for Object Restoration from Experimental Data , 1983 .

[14]  A. Jain,et al.  Extrapolation algorithms for discrete signals with application in spectral estimation , 1981 .