Abstract
In earlier work we described a circulant embedding approach for simulating scalar-valued stationary Gaussian random fields on a finite rectangular grid, with the covariance function prescribed. Here, we explain how the circulant embedding approach can be used to simulate Gaussian vector fields. As in the scalar case, the simulation procedure is theoretically exact if a certain non-negativity condition is satisfied. In the vector setting, this exactness condition takes the form of a nonnegative definiteness condition on a certain set of Hermitian matrices. The main computational tool used is the Fast Fourier Transform. Consequently, when implemented appropriately, the procedure is highly efficient, in terms of both CPU time and storage.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Chan, G. and Wood, A. T. A. (1997) An algorithm for simulating stationary Gaussian random fields. Applied Statistics, Algorithm Section, 46, 171–181.
Chan and Wood, A. T. A. (1998) Simulation of multifractional Brownian motion, pages 233–238, Proceedings in Computational Statistics 1998. R. Payne and P. J. Green (Eds). New York: Physica–Verlag.
Cressie, N. A. C. (1991) Statistics for Spatial Data. New York: John Wiley.
Davies, R. and Harte, D. S. (1987) Tests for Hurst effect. Biometrika, 74, 95–101.
Dembo, A., Mallows, C. L. and Shepp, L. A. (1989) Embedding nonnegative definite Toeplitz matrices in nonnegative definite circulant matrices, with application to covariance estimation. IEEE Trans. Inform. Theory, 35, 1206–1212.
Dietrich, C. and Newsam, G. (1993) A fast and exact method for multidimensional Gaussian stochastic simulations. Water Resources Research, 29, 2861–2869.
Feuerverger, A., Hall, P. and Wood, A. T. A. (1994) Estimation of fractal index and fractal dimension of a Gaussian process by counting the number of level crossings. J. Time Series Analysis, 15, 587–606.
NAG (1991) NAG Fortran Library Manual Mark 15, Numerical Algorithms Limited, Oxford, U.K.
Ripley, B. D. (1987) Stochastic Simulation. New York: John Wiley.
Wood, A. T. A. and Chan, G. (1994) Simulation of stationary Gaussian process in [0; 1]d. J. Comp. Graph. Stat., 3, 409–432.
Rights and permissions
About this article
Cite this article
Chan, G., Wood, A.T.A. Simulation of stationary Gaussian vector fields. Statistics and Computing 9, 265–268 (1999). https://doi.org/10.1023/A:1008903804954
Issue Date:
DOI: https://doi.org/10.1023/A:1008903804954