In the paper we consider higher-order partial differential equations from the class of linear dispersive equations. We investigate solutions to these equations subject to random initial conditions given by harmonizable φ-sub-Gaussian processes. The main results are the bounds for the distributions of the suprema for solutions. We present the examples of processes for which the assumptions of the general result are verified and bounds are written in the explicit form. The main result is also specified for the case of Gaussian initial condition.