We present in this contribution the predictions on the multiwavelength spectral energy distribution of our evolutionary population synthesis models including single and binary stellar systems. The high energy computations include the emission associated to X-ray binaries and supernovae remnants, as well as the mechanical energy released into the interstellar medium, which can be partially reprocessed into thermal X-rays. With these components we compute the spectral energy distribution of starburst galaxies from X-ray to radio ranges, and analyze finally the effects of the high energy emission on the H and He ionizing continuum.