Characterization of vector valued,gaussian, stationary,markov processes

We extend an old result by Doob characterizing real-valued, Gaussian, stationary, Markov processes to the vector case. In this case a deterministic component appears that consists of a system of harmonic oscillators while the random part is a collection of independent oscillator processes, modulo linear changes of coordinates.