Rapid singular value decomposition for time-domain analysis of magnetic resonance signals by use of the lanczos algorithm

Abstract

Time-series analysis has become an integral part of signal interpretation in magnetic resonance and other experiments. As an example, linear prediction with singular value decomposition (LPSVD) (I, 2) and the related method, Hankel SVD (HSVD) (3)) have been recently applied to problems in NMR (2)) 2D NMR (4), 2D electron spin-echo spectroscopy (2D ESE) (5,6), Fourier transform ESR correlation spectroscopy ( 7)) time-resolved femtosecond spectroscopy (8)) and nonstationary neurological currents ( 9). The minimal goal of time-series analysis is to separate signal from noise. Spectral transform methods, such as FFT filtering (IO), maximum entropy ( MEM ) ( I I ) , and LPZ ( 12)) replace the original data with a new time series (or its spectrum) that has an enhanced signal-to-noise ratio. In contrast spectral decomposition methods, such as LPSVD or HSVD, which assume a certain functional form for the time series, provide both a recipe for separating signal from noise and a listing of the harmonic components contained in the original data. This harmonic list simplifies spectral analysis and can be used to (1) extend the noise-reduced time series, (2) selectively reconstruct certain regions of the harmonic spectrum, and (3) data-compress the signal. The substantial benefits of spectral decomposition are often offset by the requirements of considerable computational time and memory storage. For a time series of N data points spectral transform techniques usually require computational time of O(N logzN) to O(N’) and storage of O(N) whereas spectral decomposition, such as HSVD, requires computational time of 0( N3) and storage of O(N2). We have investigated a new approach for solving the spectral decomposition problem which is based on the Lanczos algorithm (LA) (13). In the past the LA has provided an effective means for solving slow-motional magnetic resonance spectra ( 14, 15) as well as other complicated problems in chemical physics ( 16, 17). The LA is an extremely fast method for solving large eigenelement problems when either (1)