Context. The development of space-borne missions such as CoRoT or Kepler now provides us with numerous and precise asteroseismic measurements, which allows us to put better constraints on our theoretical knowledge of the physics of stellar interiors. In order to utilise the full potential of these measurements, however, we need a better theoretical understanding of the coupling between stellar oscillations and turbulent convection. Aims. This series of papers aims at building a new formalism specifically tailored to study the impact of turbulence on the global modes of oscillation in solar-like stars. In building this formalism, we circumvent some fundamental limitations inherent to the more traditional approaches, in particular the need for separate equations for turbulence and oscillations, and the reduction of the turbulent cascade to a unique length and time scale. This first paper aims at deriving a linear wave equation that directly and consistently contains the turbulence as an input to the model, and therefore naturally contains the information on the coupling between the turbulence and the modes, through the stochasticity of the equations. Methods. We use a Lagrangian stochastic model of turbulence based on Probability Density Function methods to describe the evolution of the properties of individual fluid particles through stochastic differential equations. We then transcribe these stochastic differential equations from a Lagrangian frame to an Eulerian frame, more adapted to the analysis of stellar oscillations. We combine this method with Smoothed Particle Hydrodynamics, where all the mean fields appearing in the Lagrangian stochastic model are estimated directly from the set of fluid particles themselves, through the use of a weighting kernel function allowing to filter the particles present in any given vicinity. The resulting stochastic differential equations on Eulerian variables are then linearised. As a first step, the gas is considered to follow a polytropic relation, and the turbulence is assumed anelastic. Results. We obtain a stochastic, linear wave equation governing the time evolution of the relevant wave variables, while at the same time containing the effect of turbulence. The wave equation generalises the classical, unperturbed propagation of acoustic waves in a stratified medium (which reduces to the exact deterministic wave equation in the absence of turbulence) to a form that, by construction, accounts for the impact of turbulence on the mode in a consistent way. The effect of turbulence consists in a non-homogeneous forcing term, responsible for the stochastic driving of the mode, and a stochastic perturbation to the homogeneous part of the wave equation, responsible for both the damping of the mode and the modal surface effects. Conclusions. The stochastic wave equation obtained here represents our baseline framework to properly infer properties of turbulenceoscillation coupling, and can therefore be used to constrain the properties of the turbulence itself with the help of asteros...