A comparison of the evolution of density fields in perturbation theory and numerical simulations – I. Non-linear evolution of the power spectrum