Infrared and Raman spectroscopies are ubiquitous techniques employed in many experimental laboratories, thanks to their fast and non-destructive nature able to capture materials' features as spectroscopic fingerprints. Nevertheless, these measurements frequently need theoretical support in order to unambiguously decipher and assign complex spectra. Linear-response theory provides an effective way to obtain the higher-order derivatives needed, but its applicability to modern exchange-correlation functionals remains limited. Here, we devise an automated, open-source, user-friendly approach based on ground-state density-functional theory and the electric enthalpy functional to allow seamless calculations of first-principles infrared and Raman spectra. By employing a finite-displacement and finite-field approach, we allow for the use of any functional, as well as an efficient treatment of large low-symmetry structures. Additionally, we propose a simple scheme for efficiently sampling the Brillouin zone with different electric fields. To demonstrate the capabilities of our approach, we provide illustrations using the ferroelectric LiNbO₃ crystal as a paradigmatic example. We predict infrared and Raman spectra using various (semi)local, Hubbard corrected, and hybrid functionals. Our results also show how PBE0 and extended Hubbard functionals yield in this case the best match in term of peak positions and intensities, respectively.