A class of weakly stationary non-Gaussian models

Two models are developed for real-valued, weakly stationary non-Gaussian processes using the spectral representation X(t)=∫0∞[cos(νt)dM1(ν)+sin(νt)dM2(ν)]. The processes MiMi, i=1,2i=1,2, in this representation are (1) square integrable martingales defined as functions Hn(Bi(ν),ν)Hn(Bi(ν),ν), ν≥0ν≥0, depending on Hermite polynomials of independent Brownian motions BiBi, i=1,2i=1,2, or (2) have increments dMi(ν)=Gi(ν)dBi(ν), ν≥0ν≥0, where GiGi are processes with finite variance independent of BiBi, i=1,2i=1,2. It is shown that the models can be calibrated to any target spectral density. Efficient methods are developed for calculating their higher order statistics. Features and limitations of the proposed models are illustrated by numerical examples. Monte Carlo algorithms for generating non-Gaussian samples based on the models in the paper resemble classical spectral representation-based algorithms for generating stationary Gaussian samples.