  The flux density spectra of radio pulsars are known to be steep and, to first
order, described by a power-law relationship of the form S_{\nu} \propto
\nu^{\alpha}, where S_{\nu} is the flux density at some frequency \nu and
\alpha is the spectral index. Although measurements of \alpha have been made
over the years for several hundred pulsars, a study of the intrinsic
distribution of pulsar spectra has not been carried out. From the result of
pulsar surveys carried out at three different radio frequencies, we use
population synthesis techniques and a likelihood analysis to deduce what
underlying spectral index distribution is required to replicate the results of
these surveys. We find that in general the results of the surveys can be
modelled by a Gaussian distribution of spectral indices with a mean of -1.4 and
unit standard deviation. We also consider the impact of the so-called
"Gigahertz-peaked spectrum" pulsars. The fraction of peaked spectrum sources in
the population with significant turn-over at low frequencies appears to be at
most 10%. We demonstrate that high-frequency (>2 GHz) surveys preferentially
select flatter-spectrum pulsars and the converse is true for lower-frequency
(<1 GHz) surveys. This implies that any correlations between \alpha and other
pulsar parameters (for example age or magnetic field) need to carefully account
for selection biases in pulsar surveys. We also expect that many known pulsars
which have been detected at high frequencies will have shallow, or positive,
spectral indices. The majority of pulsars do not have recorded flux density
measurements over a wide frequency range, making it impossible to constrain
their spectral shapes. We also suggest that such measurements would allow an
improved description of any populations of pulsars with 'non-standard' spectra.
