Let {Xt} follow a discrete Gaussian vector autoregression with deterministic components. We derive the exact finite-sample joint moment generating function (MGF) of the quadratic forms that form the basis for the sufficient statistic. The formula is then specialized to the limiting MGF of functionals involving multivariate and univariate Ornstein–Uhlenbeck processes, drifts, and time trends. Such processes arise asymptotically from more general non-Gaussian processes and also from the Gaussian {Xt} and have also been used in areas other than time series, such as the “goodness of fit” literature.