On the use of the Kalman filter in timescales

The Kalman filter is a very useful tool of estimation theory, successfully adopted in a wide variety of problems. As a recursive and optimal estimation technique, the Kalman filter seems to be the correct tool also for building precise timescales, and various attempts have been made in the past giving rise, for example, to the TA(NIST) timescale. Despite the promising expectations, a completely satisfactory implementation has never been found, due to the intrinsic non-observability of the clock time readings, which makes the clock estimation problem underdetermined. However, the case of the Kalman filter applied to the estimation of the difference between two clocks is different. In this case the problem is observable and the Kalman filter has proved to be a powerful tool. A new proposal with interesting results, concerning the definition of an independent timescale, came with the GPS composite clock, which is based on the Kalman filter and has been in use since 1990 in the GPS system. In the composite clock the indefinite growth of the covariance matrix due to the non-observability is controlled by the so-called `transparent variations'—squeezing operations on the covariance matrix that do not interfere with the estimation algorithm. A useful quantity, the implicit ensemble mean, is defined and the `corrected clocks' (physical clocks minus their predicted bias) are shown to be observable with respect to this quantity. We have implemented the full composite clock and we discuss some of its advantages and criticalities. More recently, the Kalman filter is generating new interest, and a few groups are proposing new implementations. This paper gives an overview of what has been done and of what is currently under investigation, pointing out the peculiar advantages and the open questions in the application of this attractive technique to the generation of a timescale.

[1]  R. E. Kalman,et al.  New Results in Linear Filtering and Prediction Theory , 1961 .

[2]  H. Sorenson Least-squares estimation: from Gauss to Kalman , 1970, IEEE Spectrum.

[3]  Arthur Gelb,et al.  Applied Optimal Estimation , 1974 .

[4]  Peter V Tryon,et al.  Estimating Time From Atomic Clocks. , 1983, Journal of research of the National Bureau of Standards.

[5]  A J Van Dierendonck,et al.  Relationship between Allan variances and Kalman Filter parameters , 1984 .

[6]  Peter V. Tryon,et al.  Continuous time series models for unequally spaced data applied to modeling atomic clocks , 1987 .

[7]  S. R. Stein Kalman filter analysis of precision clocks with real-time parameter estimation , 1989, Proceedings of the 43rd Annual Symposium on Frequency Control.

[8]  Patrizia Tavella,et al.  Time Scale Algorithm: Definition of Ensemble Time and Possible Uses of the Kalman Filter , 1990 .

[9]  C. T. Leondes,et al.  Ensembling clocks of Global Positioning System , 1990 .

[10]  Patrizia Tavella,et al.  Comparative Study of Time Scale Algorithms , 1991 .

[11]  T. Weissert,et al.  AT2, a new time scale algorithm : AT1 plus frequency variance , 1991 .

[12]  W. Bich Bias and Optimal Linear Estimation in Comparison Calibrations , 1992 .

[13]  J. W. Chaffee Observability, ensemble averaging and GPS time , 1992 .

[14]  Claudine Thomas Real-time Restitution of GPS Time Through a Kalman Estimation , 1993 .

[15]  S. R. Stein Advances in time-scale algorithms , 1993 .

[16]  R. J. Douglas,et al.  Time scale algorithms for an inhomogeneous group of atomic clocks , 1993 .

[17]  Walter Bich,et al.  Calibrations by comparison in metrology: A survey , 1994 .

[18]  Charles A Greenhall Kalman Plus Weights: A Time Scale Algorithm , 2001 .

[20]  K. Senior,et al.  Developing an IGS time scale , 2001, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.