We consider the anisotropic Shubin operators (−∆) m + |x| 2k acting on the space L 2 (R n ), with k, m ⩾ 1 some positive integers. We provide sharp quantitative estimates in Gelfand-Shilov spaces for the eigenfunctions of these selfadjoint differential operators with a strategy based on the classical approach to obtain Agmon estimates in spectral theory. By using a Weyl law for the eigenvalues of the anisotropic Shubin operators, we also describe the smoothing properties of the semigroups generated by the fractional powers of these operators, with precise estimates in short times. This description allows us to prove positive null-controllability results for the associated evolution equations posed on the whole space R n , from control supports which are thick with respect to densities and in any positive time. We generalize in particular results known for the evolution equations associated with fractional harmonic oscillators.Résumé. -On considère les opérateurs de Shubin anisotropes (−∆) m + |x| 2k agissant sur l'espace L 2 (R n ), avec k, m ⩾ 1 des entiers strictement positifs. On démontre des inégalités quantitatives et précises dans des espaces de Gelfand-Shilov pour les fonctions propres de ces opérateurs différentiels autoadjoints, grâce à une stratégie basée sur l'approche classique pour obtenir des estimées d'Agmon en théorie spectrale. En utilisant une loi de Weyl pour les valeurs propres des opérateurs de Shubin anisotropes, on décrit également les effets régularisants des semi-groupes engendrés par les puissances fractionnaires de ces opérateurs, et on donne des estimations précises en temps courts. Cette description permet de démontrer des résultats positifs de contrôlabilité à zéro pour les équations d'évolution associées posées sur tout l'espace R n , depuis des supports de contrôle épais par rapport à des densités, et en tout temps strictement positif. On généralise en particulier des résultats connus pour les équations d'évolution associées aux oscillateurs harmoniques fractionnaires.