Simulation of a stationary autoregression: A characterization of the normal distribution

Simulating a stationary AR(p), X_t = ∑^p_i=1α_iX_t−i + Z_t, when the innovations {Z_t} are assumed to be i.i.d. is straightforward. Starting the process in the stationary state, however, requires generation of (X₁,X₂,…,X_p) from the stationary p-dimensional distribution. When Z_t is normal this may be achieved by generating X_i as a linear function of X₁,X₂,…,X_i−1 and an independent normal variate for i = 2,3,…, p. It is shown that the ability to initialize a stationary AR(p) in this way characterizes the normal distribution.