Les satellites artificiels baignent dans un environnement radiatif hostile qui conditionne en partie leur fiabilité et leur durée de vie en opération : les ceintures de Van Allen. Afin de les protéger, il est nécessaire de caractériser la dynamique des électrons énergétiques piégés dans les ceintures radiatives. Elle est déterminée essentiellement par les interactions entre les électrons énergétiques et les ondes électromagnétiques existantes. <br /><br />Le travail de cette thèse a consisté à concevoir un schéma numérique original pour la résolution du système d'équations modélisant ces interactions : les équations de Vlasov-Maxwell relativistes. Notre choix s'est orienté vers des méthodes d'intégration directe. Nous proposons trois nouvelles méthodes spectrales pour discrétiser en impulsion les équations : une méthode de Galerkin et deux méthodes de type collocation. Ces approches sont basées sur des fonctions de Hermite qui ont la particularité de dépendre d'un facteur d'échelle permettant d'obtenir une bonne résolution en vitesse. <br /><br />Nous présentons dans ce manuscript les calculs conduisant à la discrétisation et à la résolution du problème de Vlasov-Poisson monodimensionnel ainsi que les résultats numériques obtenus. Puis nous étudions les extensions possibles des méthodes au problème complet relativiste. Afin de réduire les temps de calcul, une parallélisation et une optimisation des algorithmes ont été mises en \oe uvre. Enfin, les calculs de validation du code 1Dx-3Dv, à partir d'instabilités de types Weibel et whistlers, à une ou deux espèces d'électrons, sont détaillés.