The newly discovered inertial modes in the Sun offer the opportunity to probe the solar convective zone down to the tachocline. While linear analysis predicts the frequencies and eigenfunctions of the modes, it gives no information about their excitation or their amplitudes. We present here a theoretical formalism for the stochastic excitation of the solar inertial modes by turbulent convection. The amplitudes predicted by our model are in complete agreement with observations, thus supporting the assumption that they are stochastically excited. Our work also uncovers a qualitative transition in the shape of the inertial mode spectrum, between m ≲ 5 where the modes are clearly resolved in frequency, and m ≳ 5 where the modes overlap. This complicates the interpretation of the high-m data, and suggests that a model for the whole shape of the power spectrum is necessary to exploit the full seismic potential of solar inertial modes.