Fast and exact synthesis of stationary multivariate Gaussian time series using circulant embedding

A fast and exact procedure for the numerical synthesis of stationary multivariate Gaussian time series with a priori prescribed and well controlled auto- and cross-covariance functions is proposed. It is based on extending the circulant embedding technique to the multivariate case and can be viewed as a modification and variation around the Chan and Wood algorithm proposed earlier to solve the same problem. The procedure is shown to yield time series possessing exactly the desired covariance structure, when sufficient conditions are satisfied. Such conditions are discussed theoretically and examined on several examples of multivariate time series models. Issues related to prescribing a priori the spectral structure rather than the covariance one are also discussed. Matlab routines implementing this procedure are publicly available at http://www.hermir.org.

[1]  Howell Tong,et al.  ON TIME-REVERSIBILITY OF MULTIVARIATE LINEAR PROCESSES , 2005 .

[2]  E. Bullmore,et al.  Fractal connectivity of long-memory networks. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  Katsumi Shimotsu,et al.  Gaussian semiparametric estimation of multivariate fractionally integrated processes , 2007 .

[4]  Balachander Krishnamurthy,et al.  Sketch-based change detection: methods, evaluation, and applications , 2003, IMC '03.

[5]  C. R. Dietrich,et al.  Fast and Exact Simulation of Stationary Gaussian Processes through Circulant Embedding of the Covariance Matrix , 1997, SIAM J. Sci. Comput..

[6]  William H. Press,et al.  Numerical recipes in C , 2002 .

[7]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

[8]  Hiroshi Esaki,et al.  Live E! Project: Establishment of Infrastructure Sharing Environmental Information , 2007, 2007 International Symposium on Applications and the Internet Workshops.

[9]  Ching-Fan Chung.  SAMPLE MEANS, SAMPLE AUTOCOVARIANCES, AND LINEAR REGRESSION OF STATIONARY MULTIVARIATE LONG MEMORY PROCESSES , 2002, Econometric Theory.

[10]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[11]  Richard G. Baraniuk,et al.  Multiscale nature of network traffic , 2002, IEEE Signal Process. Mag..

[12]  P. Craigmile.  Simulating a class of stationary Gaussian processes using the Davies–Harte algorithm, with application to long memory processes , 2003 .

[13]  Donald B. Percival,et al.  Simulating Gaussian Random Processes with Specified Spectra , 2001 .

[14]  Helmut Ltkepohl,et al.  New Introduction to Multiple Time Series Analysis , 2007 .

[15]  I. M. Pyshik,et al.  Table of integrals, series, and products , 1965 .

[16]  Joseph Lipka,et al.  A Table of Integrals , 2010 .

[17]  J. Goldberger,et al.  Sympathovagal balance: how should we measure it? , 1999, The American journal of physiology.

[18]  A. T. A. Wood,et al.  Simulation of stationary Gaussian vector fields , 1999, Stat. Comput..

[19]  Colin L. Mallows,et al.  Embedding nonnegative definite Toeplitz matrices in nonnegative definite circulant matrices, with application to covariance estimation , 1989, IEEE Trans. Inf. Theory.

[20]  Patrice Abry,et al.  Testing fractal connectivity in multivariate long memory processes , 2009, 2009 IEEE International Conference on Acoustics, Speech and Signal Processing.

[21]  A. Philippe,et al.  Generators of long-range dependent processes: A survey , 2003 .

[22]  G. Didier,et al.  Integral representations and properties of operator fractional Brownian motions , 2011, 1102.1822.

[23]  Murad S. Taqqu,et al.  Theory and applications of long-range dependence , 2003 .

[24]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[25]  Donald B. Percival.  Exact simulation of complex-valued Gaussian stationary processes via circulant embedding , 2006, Signal Process..

[26]  François Ingelrest,et al.  The hitchhiker's guide to successful wireless sensor network deployments , 2008, SenSys '08.

[27]  Zhi-Quan Luo,et al.  Distributed signal processing in sensor networks , 2006 .

[28]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[29]  Tilmann Gneiting,et al.  Power-law correlations, related models for long-range dependence and their simulation , 2000, Journal of Applied Probability.

[30]  R. Davies,et al.  Tests for Hurst effect , 1987 .

[31]  Qiansheng Cheng.  Miscellanea. On time-reversibility of linear processes , 1999 .

[32]  M. Vetterli.  Distributed signal processing for sensor networks , 2005, IEEE 6th Workshop on Signal Processing Advances in Wireless Communications, 2005..

[33]  Marcus J. Chambers.  The simulation of random vector time series with given spectrum , 1995 .

[34]  Patrice Abry,et al.  Log Wavelet Leaders Cumulant Based Multifractal Analysis of EVI fMRI Time Series: Evidence of Scaling in Ongoing and Evoked Brain Activity , 2008, IEEE Journal of Selected Topics in Signal Processing.

[35]  Catalin Starica,et al.  Gaussian and Non-Gaussian Linear Time Series and Random Fields , 2001 .

[36]  A. Wood,et al.  Simulation of Stationary Gaussian Processes in [0, 1] d , 1994 .