We study the effect of Gaussian perturbations on a class of model hyperbolic partial differential equations with double symplectic characteristics in low spatial dimensions, extending some recent work in [5]. The coefficients of our partial differential operators contain harmonic oscillators in the space variables, while the noise is additive, white in time and colored in space. We provide sufficient conditions on the spectral measure of the covariance functional describing the noise that allows for the existence of a random field solution for the resulting stochastic partial differential equation. Furthermore we show how the symplectic structure of the set of multiple points affects the regularity of the noise needed to build a measurable process solution. Our approach is based on some explicit computations for the fundamental solutions of several model partial differential operators together with their explicit Fourier transforms.