Pour mener à bien leur mission, les satellites de télécommunications doivent rester à la verticale d'un même point de la Terre, sur une orbite dite géostationnaire, pour laquelle la période de révolution des satellites sur leur orbite est identique à la période de rotation de la Terre sur elle-même. Cependant, à cause des perturbations orbitales, les satellites tendent à s'en éloigner, et il est alors nécessaire de concevoir des stratégies de commande pour les maintenir dans un voisinage de cette position de référence. Du fait de leur grande valeur de poussée, les systèmes à propulsion chimique ont largement été utilisés, mais aujourd'hui les systèmes à propulsion électrique avec leur grande impulsion spécifique sont des alternatives viables pour réduire la masse d'ergols du satellite, et ainsi le coût au lancement, ou allonger la durée de vie du satellite, ce qui permettrait de limiter l'encombrement dans l'espace. Cependant, l'utilisation d'un tel système propulsif induit des contraintes opérationnelles issues en partie du caractère limité de la puissance électrique disponible à bord. Ces contraintes sont difficiles à prendre en compte dans la transcription du problème de maintien à poste en un problème de contrôle optimal à consommation minimale avec contraintes sur l'état et le contrôle. Ce manuscrit propose deux approches pour résoudre ce problème de commande optimale. La première, basée sur le développement et l'exploitation de conditions nécessaires d'optimalité, consiste à découper le problème initial en trois sous-problèmes pour former une méthode de résolution à trois étapes. La première étape permet de résoudre un problème de maintien à poste expurgé des contraintes opérationnelles, tandis que la deuxième, initialisée par le résultat de la première, produit une solution assurant le respect de ces dernières contraintes. La troisième étape permet d'optimiser la valeur des instants d'allumage et d'extinction des propulseurs dans le cadre du formalisme des systèmes à commutation. La seconde approche, dite " directe ", consiste à paramétrer le profil de commande par une fonction binaire et à le discrétiser sur l'horizon temporel de résolution. Les contraintes opérationnelles sont ainsi facilement transcrites en contraintes linéaires en nombres entiers. Après l'intégration numérique de la dynamique, le problème de contrôle optimal se résume à un problème linéaire en nombres entiers. Après la résolution du problème de maintien à poste sur un horizon court d'une semaine, le problème est résolu sur un horizon long d'un an par résolutions successives sur des horizons courts d'une durée de l'ordre de la semaine. Des contraintes de fin d'horizon court doivent alors être ajoutées afin d'assurer la faisabilité de l'enchaînement des problèmes sur l'horizon court constituant le problème sur l'horizon long. / Geostationary spacecraft have to stay above a fixed point of the Earth, on a so-called geostationary Earth orbit. For this orbit, the orbital period of the spacecraft is equal to the rotation period of the Earth. Because of orbital disturbances, spacecraft drift away their station keeping position. It is therefore mandatory to create control strategies in order to make the spacecraft stay in the vicinity of the station keeping position. Due to their high thrust capabilities, chemical thrusters have been widely used. However nowadays electric propulsion based thrusters with their high specific impulse are viable alternative in order to decrease the spacecraft mass or increase its longevity. The use of such a system induce the necessity to handle operational constraints because of the limited on-board power. These operational constraints are difficult to take into account in the mathematical transcription of the station keeping problem in an optimal control problem with control and state constraints. This thesis proposed two techniques in order to solve this optimal control problem. The first one is based on the computation of first order necessary conditions and consists in decomposing the overall problem in three sub-problems, leading to a three-step decomposition method. The first step solves an optimal control problem without the operational constraints. The second steps enforces these operational constraints thanks to dedicated equivalence schemes and the third one optimises the switching times of the control profile thanks to a method borrowed from the switched systems theory. The second proposed method consists in parametrising the on-off control profile with binary functions. After a time discretisation of the station keeping horizons, the operational constraints are easily recast as linear constraints on integer variables, the dynamics is numerically integrated and the station keeping problem is recast as a mixed integer linear programming problem. After the resolution of the problem over a short time horizon of one week, the station keeping problem is solved over a long time horizon of one year. To this end, the long time horizon is split in shorter horizons over which the problem is successively solved. End-of-cycle constraints have been set up in order to ensure the feasibility of the solution one short horizon after another.