A Levinson-type algorithm for modeling fast-sampled data

The standard discrete-time autoregressive model is poorly suited for modeling series obtained by sampling continuous-time processes at fairly rapid rates. Large computational errors can occur when the Levinson algorithm is used to estimate the parameters of this model, because the Toeplitz covariance matrix is ill-suited for inversion. An alternative model is developed based on an incremental difference operator, rather than the shift operator. It is shown that, as the sampling period goes to zero, unlike the standard autoregressive parameters, the coefficients of this model converge to certain parameters that depend directly on the statistics of the continuous-time process. A Levinson-type algorithm for efficiently estimating the parameters of this model is derived. Numerical examples are given to show that when the sampling interval is small this algorithm is considerably less sensitivity to arithmetic roundoff errors than the Levinson algorithm