Thermostated time evolutions are on a firm ground and widely used in classical molecular dynamics (MD) simulations. Hamilton´s equations of motion are supplemented by time-dependent pseudofriction terms that convert the microcanonical isoenergetic time evolution into a canonical isothermal time evolution, thus permitting the calculation of canonical ensemble averages by time averaging. However, similar methods for quantum MD schemes are still lacking. Given the rich dynamical behavior of ultracold trapped quantum gases depending on the value of the s-wave scattering length, it is timely to investigate how classical thermostating methods can be combined with powerful approximate quantum dynamics schemes to deal with interacting quantum systems at finite temperature. In this work, the popular method of Nose and Hoover to create canonically distributed positions and momenta in classical MD simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum harmonic oscillator, the equations of motion in terms of coherent states may be modified in a Nose-Hoover manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble. The method is developed initially for a single particle and then generalized to the case of an arbitrary number of identical quantum particles, involving entangled distribution functions. The resulting isothermal equations of motion for bosons and fermions contain additional terms leading to Bose-attraction and Pauli-blocking, respectively. Questions of ergodicity are discussed for different coupling schemes. In the many-particle case, the superiority of the Nose-Hoover technique to a Langevin approach is demonstrated. In addition, the work contains an investigation of the Grilli-Tosatti thermostating method applied to the harmonic oscillator, and calculations for quantum wavefunctions moving with a time-invariant shape in a harmonic potential.