De nos jours, les centrales nucléaires produisent de l'énergie par des réactions de fission (division d'un noyau atomique lourd en plusieurs noyaux atomiques légers et neutrons). Une alternative serait d'utiliser plutôt la réaction de fusion de noyaux légers de deutérium et de tritium, isotopes de l'hydrogène. Toutefois, cette technique reste encore du domaine de la recherche en physique des plasmas. Les expériences effectuées dans ce domaine ont révélé que les réacteurs à configuration magnétique toroïdale, dite tokamak, sont les plus efficaces. Un mélange gazeux d'isotopes de l'hydrogène appelé plasma est confiné grâce à un champ magnétique produit par des bobines. Ce plasma doit être chauffé à une température très élevée afin que les réactions de fusion aboutissent. De même, un courant intense doit être maintenu dans le plasma afin d'obtenir une configuration magnétique qui permet de le confiner. Une des méthodes les plus attrayantes parmi les techniques connues pour générer du courant est basée sur l'injection d'ondes électromagnétiques dans le plasma à la fréquence proche de la résonance hybride. Cette méthode offre la possibilité de contrôler le profil de densité dans le plasma. Une analyse de type Full-Wave permet alors de modéliser la propagation et l'absorption de l'onde hybride à partir des équations de Maxwell. Le but de cette thèse est de développer une méthode numérique pour cette simulation Full-Wave. Le chapitre 2 présente les équations de propagation d'ondes en mettant en évidence les caractères physiques du plasma. Une approche variationnelle de type mixte augmentée est développée et une analyse mathématique de cette dernière est effectuée dans le chapitre 3. Dans le contexte de la géométrie d'un tokamak, le problème Full-Wave dépendant de trois paramètres peut être réduit en une série de problèmes à deux variables à l'aide de la transformation de Fourier, ce sera l'objet du chapitre 4. Dans le chapitre 5, la formulation variationnelle obtenue à partir du problème mode par mode est discrétisée en utilisant des éléments finis nodaux de type Taylor-Hood. Le chapitre 6 concerne les méthodes de résolution du système linéaire après discrétisation. À l'aide de différents diagnostics physiques présentés dans le chapitre 7, des résultats de la simulation Full-Wave obtenues à partir d'un code MATLAB sont présentées dans le chapitre 8. Enfin, dans le but de développer une version parallèle de la simulation, le chapitre 9 est consacré à une méthode de décomposition de domaine sans recouvrement associé au système Full-Wave.