In this paper, we show how to carry out a relatively more realistic and complete reconstruction of supernova neutrino spectra in the future large liquid-scintillator detectors, by implementing the method of singular value decomposition with a proper regularization. For a core-collapse supernova at a distance of 10 kpc in the Milky Way, itsν e spectrum can be precisely determined from the inverse beta-decay processν e þ p → e þ þ n, for which a 20 kiloton liquid-scintillator detector with the resolution similar to the Jiangmen Underground Neutrino Observatory may register more than 5000 events. We have to rely predominantly on the elastic neutrino-electron scattering ν þ e − → ν þ e − and the elastic neutrino-proton scattering ν þ p → ν þ p for the spectra of ν e and ν x , where ν denotes collectively neutrinos and antineutrinos of all three flavors and ν x for ν μ and ν τ as well as their antiparticles. To demonstrate the validity of our approach, we also attempt to reconstruct the neutrino spectra by using the time-integrated neutrino data from the latest numerical simulations of delayed neutrino-driven supernova explosions.