FR 3041103 A1 20170317 FR 1558621 A 20150915 1 DOMAINE TECHNIQUE La présente invention appartient au domaine du contrôle d'attitude de satellites en orbite terrestre, et concerne plus particulièrement un procédé et un système d'estimation d'attitude d'un satellite de télécommunications. ÉTAT DE LA TECHNIQUE Dans un système de télécommunications par satellite, une station sol émet un signal utile à destination d'un satellite en orbite terrestre, par exemple géostationnaire (« Geostationary Orbit » ou GEO dans la littérature anglo-saxonne). Puis le satellite retransmet le signal utile vers la Terre, généralement après l'avoir décalé en fréquences, à destination d'une autre station sol. Par exemple, dans un système de diffusion par satellite (« broadcast » ou « broadband »), le satellite reçoit d'une station sol un signal utile à diffuser, tel qu'un signal DVB-S, et le retransmet à destination d'une pluralité de stations sol, telles que des terminaux utilisateurs. Un tel satellite est prévu pour desservir, sur le lien descendant à destination des terminaux utilisateurs, une zone de couverture prédéterminée. A cet effet, le satellite comporte une antenne d'émission, pour la diffusion du signal utile, dont le diagramme de rayonnement par rapport audit satellite, lorsque ledit satellite est placé dans une attitude de mission prédéterminée, présente, dans la direction de chaque point à l'intérieur ladite zone de couverture, un gain supérieur à celui en direction d'un point à l'extérieur de ladite zone de couverture. On comprend donc qu'il est nécessaire, pour assurer que les terminaux utilisateurs qui se trouvent dans la zone de couverture prévue puissent recevoir le signal utile, de contrôler avec précision l'attitude du satellite, afin que celle-ci soit la plus proche possible de l'attitude de mission. A cet effet, les satellites actuels sont équipés d'un ensemble de capteurs permettant d'estimer avec précision l'attitude dudit satellite, et également un ensemble d'actionneurs permettant de modifier l'attitude du satellite pour réduire l'écart entre l'attitude estimée et l'attitude de mission. Notamment, il est connu d'utiliser des capteurs optiques (capteur d'horizon terrestre ou solaire, capteur stellaire, etc.) pour estimer l'attitude de manière absolue, par opposition à d'autres types de capteurs qui ne permettent 3 0 4 1 1 0 3 2 que d'estimer les variations d'attitude entre deux instants. Les capteurs optiques présentent l'avantage d'une précision élevée, mais présentent l'inconvénient d'être aveugles lorsqu'ils sont face à un trou d'étoiles ou éblouis par un astre constituant une source de lumière intense. En 5 outre, l'estimation d'attitude n'est plus possible en cas de panne temporaire ou définitive d'un ou de plusieurs capteurs. Le brevet américain US 4630058 décrit un procédé d'estimation d'attitude d'un satellite, dans lequel ledit satellite est configuré pour générer à bord deux signaux qui sont émis par des antennes d'émission respectives 10 différentes, dont l'antenne d'émission des signaux utiles. Les signaux ainsi émis sont reçus par quatre stations sol, et les signaux reçus sont comparés afin d'estimer l'attitude dudit satellite. Une telle solution est avantageuse dans la mesure où elle ne met pas en oeuvre de capteurs à bord du satellite. 15 Toutefois, la solution décrite par le brevet US 4630058 repose sur un satellite ayant une architecture très spécifique, et présente par conséquent l'inconvénient de ne pas être compatible avec les satellites actuellement à poste en orbite GEO. EXPOSÉ DE L'INVENTION 20 La présente invention a pour objectif de remédier à tout ou partie des limitations des solutions de l'art antérieur, notamment celles exposées ci-avant, en proposant une solution pour l'estimation d'attitude d'un satellite qui ne nécessite pas d'utiliser des mesures de capteurs à bord du satellite, et qui soit en outre compatible avec les satellites actuellement à poste en orbite terrestre. 25 A cet effet, et selon un premier aspect, l'invention concerne un procédé d'estimation d'attitude d'un satellite, ledit satellite étant configuré pour générer et pour émettre un premier signal de référence sur un lien descendant au moyen d'une première antenne d'émission, ledit satellite étant en outre configuré pour recevoir des signaux utiles sur un lien montant, et pour 30 retransmettre lesdits signaux utiles sur le lien descendant au moyen d'une seconde antenne d'émission, la première antenne d'émission desservant une zone de couverture englobant une zone de couverture de la seconde antenne d'émission. Le procédé d'estimation comporte : 3041103 3 - l'émission sur le lien montant d'un second signal de référence, qui est reçu par ledit satellite au moyen d'une antenne de réception des signaux utiles, - la retransmission par ledit satellite du second signal de référence 5 au moyen de la seconde antenne d'émission, - la réception du premier signal de référence et du second signal de référence par au moins trois stations sol se trouvant dans la zone de couverture de la seconde antenne d'émission, dites « stations de mesure », 10 - le calcul, pour chaque station de mesure, d'une grandeur représentative d'un rapport entre les puissances respectives du premier signal de référence et du second signal de référence reçus, dite « mesure différentielle de puissance », - l'estimation de l'attitude du satellite en fonction des mesures 15 différentielles de puissance calculées. Ainsi, le procédé d'estimation d'attitude repose également sur la réception de signaux émis par le satellite, et ne nécessite pas d'utiliser des mesures de capteurs à bord du satellite. En outre, le procédé d'estimation est compatible avec les satellites 20 actuellement à poste en orbite terrestre. En effet, seul l'un des signaux émis par le satellite est généré à bord, en l'occurrence le premier signal de référence. De tels signaux existent déjà dans tous les satellites actuels. On peut citer notamment le signal de télémesure ou le signal balise, qui sont générés à bord de chaque satellite et émis sur le lien descendant au moyen 25 d'antennes d'émission peu directives. Le signal de télémesure ou le signal balise peuvent être utilisés en tant que premier signal de référence dans le procédé d'estimation selon l'invention. L'autre signal, c'est-à-dire le second signal de référence, est lui généré au sol et est reçu par le satellite sur le lien montant, avec des signaux utiles. Il 30 est retransmis sur le lien descendant avec lesdits signaux utiles. Par conséquent, un tel second signal de référence peut être reçu et retransmis par les satellites actuels, par les mêmes moyens que ceux mis en oeuvre pour la réception et la retransmission des signaux utiles. 3041103 4 Les puissances du premier signal de référence et du second signal de référence reçus dépendent notamment du gain de la première antenne d'émission et du gain de la seconde antenne d'émission. Elles dépendent également des pertes de propagation entre le satellite et la Terre, qui ne sont 5 pas connues a priori mais affectent sensiblement de la même façon le premier signal de référence et le second signal de référence s'ils sont émis, par le satellite, suffisamment proches dans les domaines fréquentiel et temporel. Par conséquent, si les pertes de propagation sont sensiblement les mêmes pour le premier signal de référence et pour le second signal de référence, alors le 10 calcul des mesures différentielles de puissance permet d'annuler lesdites pertes de propagation sans avoir à les estimer. En outre, du fait que le premier signal de référence et le second signal de référence sont émis par des antennes d'émission respectives différentes, les gains respectifs desdites première antenne d'émission et seconde antenne d'émission dans la direction 15 d'une station de mesure sont différents et ne s'annulent donc pas en calculant lesdites mesures différentielles de puissance. Ensuite, l'attitude du satellite est estimée en fonction desdites mesures différentielles de puissance. En effet, une variation de l'attitude du satellite entraîne un déplacement des diagrammes de rayonnement respectifs de la 20 première antenne d'émission et de la seconde antenne d'émission, ce qui entraine une variation des gains desdites première antenne d'émission et seconde antenne d'émission dans les directions des stations de mesure. Les positions respectives du satellite et des stations de mesure sont généralement connues, de même que les diagrammes de rayonnement 25 respectifs de la première antenne d'émission et de la seconde antenne d'émission. Ainsi, on peut déterminer, pour chaque attitude possible du satellite, des gains théoriques desdites première antenne d'émission et seconde antenne d'émission dans les directions des stations de mesure, et en déduire des mesures différentielles de puissance théoriques associées 30 respectivement auxdites stations de mesure. L'attitude du satellite peut alors être estimée, par exemple, en comparant les mesures différentielles de puissance réelles, c'est-à-dire celles calculées à partir du premier signal de référence et du second signal de référence effectivement reçus par les stations 3041103 5 de mesure, aux mesures différentielles de puissance théoriques. Dans des modes particuliers de mise en oeuvre, le procédé d'estimation d'attitude peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons 5 techniquement possibles. Dans des modes particuliers de mise en oeuvre, ledit procédé comporte le calcul, pour au moins deux paires de stations de mesure, d'une grandeur représentative d'un rapport entre les mesures différentielles de puissance respectives des stations de mesure de ladite paire, dite « mesure 10 différentielle inter-station », et l'attitude du satellite est estimée en fonction desdites mesures différentielles inter-station calculées. Les puissances du premier signal de référence et du second signal de référence reçus dépendent également de leurs puissances d'émission respectives par le satellite, qui ne sont pas forcément connues a priori et/ou 15 sont connues mais avec une précision qui peut s'avérer insuffisante. Par contre, lesdites puissances d'émission affectent de la même façon les puissances reçues par les différentes stations de mesure. Par conséquent, le calcul des mesures différentielles inter-station permet d'annuler les incertitudes liées à une méconnaissance totale ou partielle desdites puissances d'émission. 20 Dans des modes particuliers de mise en oeuvre, ledit procédé comporte la synchronisation temporelle des stations de mesure entre elles, et les mesures différentielles inter-station sont calculées à partir de mesures différentielles de puissance calculées pour des signaux reçus en même temps par les stations de mesure. 25 Les puissances d'émission par le satellite peuvent éventuellement varier au cours du temps. Si ces variations sont trop rapides, alors il peut s'avérer avantageux de synchroniser les stations de mesure afin de s'assurer de comparer des mesures différentielles de puissance pour lesquelles les puissances d'émission par le satellite sont effectivement les mêmes d'une 30 station de mesure à une autre. De telles dispositions permettent par conséquent d'améliorer la précision de l'estimation de l'attitude du satellite. Dans des modes particuliers de mise en oeuvre, les stations de mesure comportant des chaînes de réception respectives, ledit procédé 3041103 6 comporte, pour chaque station de mesure, l'estimation d'une erreur introduite par ladite chaîne de réception sur le calcul de la mesure différentielle de puissance, dite « erreur différentielle de puissance », et la compensation de l'erreur différentielle de puissance sur la mesure différentielle de puissance. 5 Dans des modes particuliers de mise en oeuvre, l'estimation de l'erreur différentielle de puissance comporte, pour l'une au moins des stations de mesure, l'injection dans la chaîne de réception de ladite station de mesure d'un premier signal de calibration et d'un second signal de calibration, et le calcul d'une grandeur représentative d'un rapport entre les puissances respectives du 10 premier signal de calibration et du second signal de calibration mesurés en sortie de ladite chaîne de réception. Dans des modes particuliers de mise en oeuvre, le premier signal de calibration et le second signal de calibration sont injectés sur des fréquences centrales respectives F'1 et F'2 différentes. 15 Dans des modes particuliers de mise en oeuvre, ledit procédé comporte la génération du premier signal de calibration et du second signal de calibration par mélange non linéaire de deux signaux de fréquences centrales respectives Fg1 et Fg2 telles que : - IFg1 - Fg2I = F'1 et IFg1 + Fg2I = F'2, ou 20 - IFg1 - Fg2I = F'2 et IFg1 + Fg2I = F'1. Dans des modes particuliers de mise en oeuvre, le premier signal de référence et le second signal de référence sont émis sur des fréquences centrales respectives F1 et F2 différentes sur le lien descendant, et de préférence simultanément. 25 Dans des modes particuliers de mise en oeuvre, le premier signal de référence et/ou le second signal de référence sont des signaux sinusoïdaux. Selon un second aspect, la présente invention concerne un système d'estimation d'attitude d'un satellite, ledit satellite comportant des moyens configurés pour générer et pour émettre un premier signal de référence sur un 30 lien descendant au moyen d'une première antenne d'émission, recevoir des signaux utiles sur un lien montant et retransmettre lesdits signaux utiles sur le lien descendant au moyen d'une seconde antenne d'émission, ladite première antenne d'émission desservant une zone de couverture englobant une zone de 3041103 7 couverture de ladite seconde antenne d'émission. Le système d'estimation comporte : - une station sol, dite « station d'émission de référence », configurée pour émettre un second signal de référence sur le lien montant, 5 adapté à être reçu par le satellite avec les signaux utiles, - au moins trois stations sol adaptées à recevoir le premier signal de référence émis au moyen de la première antenne d'émission et le second signal de référence retransmis au moyen de la seconde antenne d'émission, dites « stations de mesure », 10 - des moyens configurés pour calculer, pour chaque station de mesure, une grandeur représentative d'un rapport entre les puissances respectives du premier signal de référence et du second signal de référence reçus, dite « mesure différentielle de puissance », 15 - des moyens configurés pour estimer l'attitude du satellite en fonction des mesures différentielles de puissance calculées. Dans des modes particuliers de réalisation, le système d'estimation d'attitude peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement 20 possibles. Dans des modes particuliers de réalisation, ledit système comporte des moyens configurés pour calculer, pour au moins deux paires de stations de mesure, une grandeur représentative d'un rapport entre les mesures différentielles de puissance respectives des stations de mesure de ladite paire, 25 dite « mesure différentielle inter-station », et l'attitude du satellite est estimée en fonction desdites mesures différentielles inter-station calculées. Dans des modes particuliers de réalisation, les stations de mesure sont synchronisées temporellement entre elles. Dans des modes particuliers de réalisation, les trois stations de 30 mesure considérées sont situées, dans la zone de couverture de la seconde antenne d'émission du satellite placé dans une attitude de mission, de telle sorte que la valeur du gradient du diagramme de rayonnement de ladite seconde antenne d'émission dans la direction d'une desdites trois stations de 3041103 8 mesure, dite « station de mesure de référence », est inférieure aux valeurs dudit gradient du diagramme de rayonnement dans les directions respectives des deux autres stations de mesure. Dans des modes particuliers de réalisation, chaque station de mesure 5 comportant une chaîne de réception, ledit système comporte des moyens configurés pour estimer, pour chaque station de mesure, une erreur introduite par ladite chaîne de réception sur le calcul de la mesure différentielle de puissance, dite « erreur différentielle de puissance », et pour compenser ladite erreur différentielle de puissance sur la mesure différentielle de puissance. 10 Dans des modes particuliers de réalisation, l'une au moins des stations de mesure comporte des moyens configurés pour générer et pour injecter un premier signal de calibration et un second signal de calibration dans la chaîne de réception de ladite station de mesure, et des moyens configurés pour calculer une grandeur représentative d'un rapport entre les puissances 15 respectives du premier signal de calibration et du second signal de calibration mesurés en sortie de ladite chaîne de réception. Dans des modes particuliers de réalisation, le premier signal de calibration et le second signal de calibration étant injectés sur des fréquences centrales respectives F'1 et F'2 différentes, les moyens de génération dudit 20 premier signal de calibration et dudit second signal de calibration sont configurés pour mélanger de manière non linéaire deux signaux de fréquences centrales respectives Fg1 et Fg2 telles que : - IFg1 - Fg2I = F'1 et IFg1 + Fg2I = F'2, ou - IFg1 - Fg2I = F'2 et IFg1 + Fg2I = F'1. 25 PRÉSENTATION DES FIGURES L'invention sera mieux comprise à la lecture de la description suivante, donnée à titre d'exemple nullement limitatif, et faite en se référant aux figures qui représentent : - Figure 1 : une représentation schématique d'un système 30 d'estimation d'attitude d'un satellite, - Figure 2 : une représentation schématique d'un exemple de zone de couverture d'un satellite et de positionnement de stations sol du système d'estimation d'attitude, 3041103 9 - Figure 3: un diagramme illustrant les principales étapes d'un procédé d'estimation d'attitude, - Figure 4 : une représentation schématique d'un mode préféré de réalisation d'une station sol du système d'estimation d'attitude. 5 Dans ces figures, des références identiques d'une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté, les éléments représentés ne sont pas à l'échelle, sauf mention contraire. DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION La figure 1 représente schématiquement un système 10 d'estimation fo d'attitude d'un satellite 20 en orbite autour de la Terre T. Dans la suite de la description, on se place de manière non limitative dans le cas où le satellite 20 est en orbite GEO. Rien n'exclut cependant, suivant d'autres exemples, de considérer un satellite en orbite non-géostationnaire, par exemple une orbite géosynchrone quelconque, une orbite 15 basse (« Low Earth Orbit » ou LEO dans la littérature anglo-saxonne), une orbite moyenne (« Medium Earth Orbit » ou MEO dans la littérature anglo-saxonne), etc., sous réserve de connaître la position dudit satellite lorsque son attitude doit être estimée. Tel qu'illustré par la figure 1, le satellite 20 comporte une première 20 antenne d'émission 21 d'un premier signal de référence, généré à bord dudit satellite 20, sur un lien descendant à destination de la Terre T. En outre, le satellite 20 comporte une charge utile comportant notamment une antenne de réception 22 pour la réception de signaux utiles sur un lien montant entre la Terre T et ledit satellite 20. La charge utile comporte 25 également une seconde antenne d'émission 23 au moyen de laquelle les signaux utiles reçus sont retransmis sur le lien descendant. Les traitements effectués, avant retransmission, sur les signaux utiles reçus sortent du cadre de l'invention. Il est à noter que l'invention, dans son principe, ne nécessite pas de traitements spécifiques à bord et trouve une 30 application particulièrement avantageuse dans le cas de satellites dits « transparents » (c'est-à-dire dont les traitements effectués à bord sur les signaux utiles comportent principalement l'introduction d'un décalage fréquentiel et une amplification) actuellement à poste en orbite terrestre. Il est 3 0 4 1 1 0 3 10 cependant également possible de considérer des satellites dits « régénératifs », en particulier pour des satellites qui n'ont pas encore été lancés. Par contre, la compatibilité avec des satellites régénératifs actuellement à poste n'est pas forcément assurée, et dépend du type de traitements 5 spécifiques effectués à bord. Afin d'effectuer sa mission de collecte et de retransmission de signaux utiles, le satellite 20 est placé dans une attitude de mission prédéterminée. Lorsque le satellite 20 est placé dans l'attitude de mission, la seconde antenne d'émission 23 dessert une zone de couverture Zc2. 10 Dans l'exemple illustré par la figure 1, il est considéré de manière non limitative que l'antenne de réception 22 dessert la même zone de couverture Zc2 que la seconde antenne d'émission 23. Dans un tel cas, l'antenne de réception 22 et la seconde antenne d'émission 23 peuvent être, dans des variantes de réalisation, une seule et même antenne utilisée à la fois pour la 15 réception et la retransmission de signaux utiles. La première antenne d'émission 21, mise en oeuvre pour l'émission du premier signal de référence, dessert une zone de couverture Zc1. Tel qu'illustré par la figure 1, la zone de couverture Zc1 englobe la zone de couverture Zc2 de la seconde antenne d'émission 23. 20 Le système 10 d'estimation comporte en outre une station sol, dite « station d'émission de référence » 30, configurée pour émettre un second signal de référence sur le lien montant, adapté à être reçu par le satellite 20 avec les signaux utiles. En particulier, le satellite 20 est configuré pour recevoir les signaux utiles émis, sur le lien montant, dans une bande fréquentielle 25 prédéfinie, dite « bande fréquentielle montante ». Par conséquent, ledit second signal de référence est également émis dans ladite bande fréquentielle montante, sur une fréquence centrale F0, afin de pouvoir être reçu par ledit satellite 20. Le second signal de référence est de préférence émis de sorte à 30 limiter les interférences sur les signaux utiles. Par exemple, le second signal de référence est un signal de faible largeur spectrale, de préférence un signal sensiblement sinusoïdal, afin de pouvoir être émis dans une partie de la bande fréquentielle montante dans 3041103 11 laquelle aucun signal utile n'est émis, mais que le satellite 20 peut néanmoins recevoir, par exemple en bordure de ladite bande fréquentielle montante. Alternativement, le second signal de référence peut être un signal à spectre étalé au moyen d'une séquence d'étalement prédéfinie. En effet, la 5 densité spectrale de puissance d'un signal à spectre étalé peut être rendue très faible, limitant ainsi les interférences sur les signaux utiles. En outre, la corrélation par la séquence d'étalement permet d'introduire un gain de traitement qui permet d'extraire ledit second signal de référence malgré sa faible densité spectrale de puissance. 10 Généralement, la charge utile d'un satellite 20 de télécommunications introduit un décalage fréquentiel de sorte qu'un signal utile est retransmis sur le lien descendant sur une fréquence centrale différente de celle sur laquelle il a été reçu sur le lien montant. Ainsi, la bande fréquentielle dans laquelle les signaux utiles sont retransmis sur le lien descendant, dite « bande fréquentielle 15 descendante », est généralement différente de la bande fréquentielle montante, et le second signal de référence est retransmis sur une fréquence centrale F2 différente de la fréquence centrale F0. De préférence, la fréquence centrale FO du second signal de référence est choisie de telle sorte que la fréquence centrale F2 est proche d'une 20 fréquence centrale F1 sur laquelle est émis le premier signal de référence, lequel n'est pas nécessairement émis à l'intérieur de la bande fréquentielle descendante des signaux utiles. En effet, les pertes de propagation dépendent de la fréquence utilisée et il est par conséquent avantageux de considérer des fréquences centrales F1 et F2 proches afin d'assurer que les pertes de 25 propagation subies, sur le lien descendant, par le premier signal de référence et par le second signal de référence sont sensiblement les mêmes. Par exemple, l'écart fréquentiel entre les fréquences centrales F1 et F2 est de l'ordre de quelques dizaines de mégahertz (MHz), voire inférieur à 10 MHz. Dans le cas d'une mise en oeuvre pour estimer l'attitude d'un satellite 30 20 existant, déjà à poste, le premier signal de référence est choisi parmi les signaux, autres que les signaux utiles, émis par le satellite 20 au moyen d'une antenne d'émission différente de celle utilisée pour la retransmission des signaux utiles. Par exemple, les satellites actuels émettent généralement un 3041103 12 signal de télémesure et un signal balise au moyen d'une antenne d'émission peu directive, typiquement une antenne cornet, dont la zone de couverture englobe la zone de couverture de l'antenne d'émission des signaux utiles. En outre, dans les satellites actuels, le signal de télémesure et le signal balise sont 5 émis sur des fréquences centrales proches de la bande fréquentielle descendante des signaux utiles. Le signal de télémesure et/ou le signal balise peuvent donc être utilisés en tant que premier signal de référence. Rien n'exclut cependant, dans le cas notamment d'un satellite non encore à poste en orbite terrestre, 10 d'embarquer dans ledit satellite des moyens spécifiques pour l'émission du premier signal de référence, distincts des moyens d'émission du signal de télémesure et du signal balise. Dans la suite de la description, on se place de manière non limitative dans le cas où les fréquences centrales F1 et F2 sont distinctes. Rien n'exclut 15 cependant, suivant d'autres exemples, de considérer la fréquence centrale F1 égale à la fréquence centrale F2, par exemple si le premier signal de référence et le second signal de référence sont multiplexés temporellement (c'est-à-dire qu'ils ne sont pas émis simultanément) et/ou multiplexés par étalement de spectre au moyen de séquences d'étalement orthogonales, etc. 20 Dans des modes particuliers de mise en oeuvre, le premier signal de référence et/ou le second signal de référence sont des signaux à faible largeur spectrale, par exemple des signaux sinusoïdaux. De telles dispositions permettent, lorsque les fréquences centrales F1 et F2 sont proches, d'assurer une cohérence fréquentielle des pertes de propagation. En outre, l'estimation 25 de l'attitude du satellite 20 s'appuie sur des mesures des puissances avec lesquelles ledit premier signal de référence et ledit second signal de référence sont reçus par des stations sol, et des signaux sinusoïdaux permettent de simplifier lesdites mesures de puissance. Dans la suite de la description, on se place de manière non limitative 30 dans le cas où le premier signal de référence et le second signal de référence, émis sur le lien descendant, sont des signaux sinusoïdaux de fréquences centrales respectives F1 et F2. En outre, le premier signal de référence et le second signal de référence sont avantageusement émis simultanément, afin 3041103 13 d'assurer une meilleure cohérence temporelle des pertes de propagation. Par exemple, si la station d'émission de référence 30 émet le second signal de référence en continu, alors il sera également retransmis en continu par le satellite 20, il se trouvera toujours des intervalles de temps au cours desquels 5 le premier signal de référence et le second signal de référence sont émis simultanément par ledit satellite 20. Tel qu'illustré par la figure 1, le système 10 de contrôle d'attitude comporte également au moins trois stations sol, dites « stations de mesure » 31, adaptées à recevoir le premier signal de référence et le second signal de 10 référence émis par le satellite 20. Dans l'exemple illustré par la figure 1, la station d'émission de référence 30 est distincte des stations de mesure 31. Rien n'exclut cependant, suivant d'autres exemples d'avoir une station sol qui soit à la fois une station d'émission de référence et une station de mesure. Lesdites trois stations de mesure 31 ne sont pas alignées avec le 15 satellite 20. En d'autres termes, lesdites trois stations de mesure 31 et le satellite 20 ne se trouvent tous dans le même plan. De telles dispositions permettent d'estimer l'attitude du satellite 20 suivant trois axes, désignés généralement par axes de roulis, de tangage et de lacet. En effet, l'erreur d'attitude en roulis d'un satellite 20 en orbite GEO 20 entraîne un déplacement de la zone de couverture Zc2 de la seconde antenne d'émission 23 suivant l'axe Nord-Sud, alors que l'erreur d'attitude en tangage entraîne un déplacement suivant l'axe Est-Ouest. En outre, l'erreur d'attitude d'un satellite 20 suivant le troisième axe, désigné généralement par axe de lacet et qui relie le satellite 20 au centre de la Terre T, entraîne une rotation de 25 ladite zone de couverture Zc2 autour dudit axe de lacet. De tels déplacements peuvent être détectés au moyen de trois stations de mesure 31 qui ne sont pas alignées avec le satellite 20. Toutefois, en pratique, l'erreur d'attitude en lacet d'un satellite 20 en orbite GEO est généralement négligeable. La figure 2 représente schématiquement un exemple de zone de 30 couverture Zc2 de la seconde antenne d'émission 23 du satellite 20. Plus particulièrement, la figure 2 représente la projection du diagramme de rayonnement de ladite seconde antenne d'émission 23 sur la surface de la Terre T, le satellite 20 étant dans son attitude de mission. 3041103 14 Tel qu'illustré par la figure 2, la projection du diagramme de rayonnement de la seconde antenne d'émission 23 peut être représentée sous la forme de courbes de niveau qui joignent les points à la surface de la Terre T pour lesquels le gain du diagramme de rayonnement est le même. 5 Dans l'exemple illustré par la figure 2, plusieurs courbes de niveau sont représentées, à savoir des courbes de niveau correspondant respectivement à des atténuations de 1 décibel (dB), 2 dB, 3dB, 4 dB et 5 dB par rapport au gain maximal G mAx du diagramme de rayonnement de la seconde antenne d'émission 23. Dans l'exemple non limitatif illustré par la 10 figure 2, la zone de couverture Zc2 est définie comme étant la zone pour laquelle l'atténuation par rapport au gain maximal G mAx n'est pas supérieure à trois décibels, c'est-à-dire la zone pour laquelle le gain est égal ou supérieur à (G mAx - 3 dB). Les trois stations de mesure 31 sont par exemple situées, dans la 15 zone de couverture Zc2 de la seconde antenne d'émission 23, de telle sorte que la valeur du gradient du diagramme de rayonnement de ladite seconde antenne d'émission 23 dans la direction d'une desdites trois stations de mesure, dite « station de mesure de référence », est inférieure aux valeurs dudit gradient du diagramme de rayonnement dans les directions respectives 20 des deux autres stations de mesure, voire très inférieure à celles-ci. En d'autres termes, la station de mesure de référence est placée dans une zone où les courbes de niveau sont moins resserrées que dans les zones où sont placées les deux autres stations de mesure. Dans l'exemple illustré par la figure 2, la station de mesure de 25 référence, désignée par 31a, est située sensiblement au centre de la zone de couverture Zc2, et les deux autres stations de mesure, désignées respectivement par 31b et 31c, sont situées sensiblement en bordure de la zone de couverture Zc2, dans des zones où les courbes de niveau sont resserrées. Plus particulièrement, dans l'exemple illustré par la figure 2, la 30 station de mesure 31b se trouve dans une zone dans laquelle les courbes de niveau sont resserrées suivant l'axe Est-Ouest, et la station de mesure 31c se trouve dans une zone dans laquelle les courbes de niveau sont resserrées suivant l'axe Nord-Sud. Ainsi, une erreur sur l'attitude en roulis et/ou en 3041103 15 tangage du satellite 20, par rapport à l'attitude de mission : - modifie peu le gain de la seconde antenne d'émission 23 dans la direction de la station de mesure de référence 31a par rapport au gain théorique attendu lorsque le satellite 20 est dans son attitude 5 de mission, - modifie fortement le gain de la seconde antenne d'émission 23 dans la direction de la station de mesure 31b et/ou dans la direction de la station de mesure 31c par rapport aux gains théoriques attendus lorsque le satellite 20 est dans son attitude de 10 mission. Le système 10 d'estimation d'attitude comporte également un ensemble de moyens configurés pour estimer l'attitude du satellite 20 en fonction du premier signal de référence et du second signal de référence reçus par les différentes stations de mesure 31. Ces moyens peuvent être intégrés 15 dans un même équipement ou distribués sur plusieurs équipements. Au niveau de chaque équipement, lesdits moyens se présentent par exemple sous la forme d'un module de traitement (non représenté sur les figures) comportant par exemple un ou plusieurs processeurs et des moyens de mémorisation (disque dur magnétique, mémoire électronique, disque 20 optique, etc.) dans lesquels est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter pour mettre en oeuvre une ou plusieurs étapes d'estimation de l'attitude du satellite 20. Alternativement ou en complément, le module de traitement peut comporter un ou des circuits logiques programmables (FPGA, PLD, etc.), et/ou 25 un ou des circuits intégrés spécialisés (ASIC), et/ou un ensemble de composants électroniques discrets, etc., adaptés à mettre en oeuvre une ou plusieurs desdites étapes d'estimation de l'attitude du satellite 20. En d'autres termes, lesdits moyens sont configurés de façon logicielle (produit programme d'ordinateur spécifique) et/ou matérielle (FPGA, PLD, 30 ASIC, composants électroniques discrets, etc.) pour estimer l'attitude du satellite 20 en fonction du premier signal de référence et du second signal de référence reçus par les différentes stations de mesure 31. Dans la suite de la description, on se place de manière non limitative 3041103 16 dans le cas où ces moyens sont distribués sur les stations de mesure 31 et sur une station centrale (non représentée sur les figures). La figure 3 représente schématiquement les principales étapes d'un procédé 50 d'estimation de l'attitude d'un satellite 20. Tel qu'illustré par la 5 figure 3, le procédé 50 d'estimation d'attitude comporte tout d'abord des étapes de : - 51 émission sur le lien montant, par la station d'émission de référence 30, du second signal de référence, qui est reçu par le satellite 20 au moyen de l'antenne de réception 22 des signaux 10 utiles, - 52 retransmission par ledit satellite 20 du second signal de référence sur le lien descendant au moyen de la seconde antenne d'émission 23, - 53 émission par ledit satellite 20 du premier signal de référence sur 15 le lien descendant au moyen de la première antenne d'émission 21 - 54 réception du premier signal de référence et du second signal de référence par les au moins trois stations de mesure 31. Le procédé 50 d'estimation d'attitude comporte ensuite une étape 55 20 de calcul, pour chaque station de mesure 31, d'une grandeur représentative d'un rapport entre les puissances respectives du premier signal de référence et du second signal de référence reçus, dite « mesure différentielle de puissance ». Le calcul de telles mesures différentielles de puissance permet d'annuler, notamment, les pertes de propagation sur le lien descendant. 25 L'étape 55 de calcul de mesure différentielle de puissance est par exemple exécutée par chaque station de mesure 31. Le procédé 50 d'estimation d'attitude comporte ensuite une étape 56 d'estimation de l'attitude du satellite en fonction des mesures différentielles de puissance calculées. L'étape 56 d'estimation d'attitude est par exemple 30 exécutée par la station centrale, en fonction des mesures différentielles de puissance reçues des différentes stations de mesure 31. Par exemple, l'attitude du satellite 20 est estimée en fonction d'une connaissance a priori des positions respectives du satellite 20 et des stations 3041103 17 de mesure 31, et des diagrammes de rayonnement respectifs de la première antenne d'émission 21 et de la seconde antenne d'émission 23 par rapport audit satellite 20. Ainsi, on peut déterminer, pour chaque attitude possible du satellite 20, des gains théoriques desdites première antenne d'émission 21 et 5 seconde antenne d'émission dans les directions des stations de mesure 31, et en déduire des mesures différentielles de puissance théoriques associées respectivement auxdites stations de mesure 31. L'attitude du satellite 20 peut alors être estimée, par exemple, en comparant les mesures différentielles de puissance réelles, c'est-à-dire celles calculées à partir du premier signal de 10 référence et du second signal de référence effectivement reçus par les stations de mesure 31, aux différentes mesures différentielles de puissance théoriques. Dans des modes préférés de mise en oeuvre, l'étape 56 d'estimation d'attitude comporte le calcul, pour au moins deux paires de stations de mesure 31, d'une grandeur représentative d'un rapport entre les mesures différentielles 15 de puissance respectives des stations de mesure 31 de ladite paire, dite « mesure différentielle inter-station ». Les paires de stations de mesure 31 considérées sont différentes, c'est-à-dire qu'au moins une station de mesure 31 diffère d'une paire à une autre. Par contre, il est possible de considérer des paires comportant toutes 20 une même station de mesure 31, par exemple comportant toutes la station de mesure de référence 31a dans l'exemple illustré par la figure 2. L'attitude du satellite 20 est ensuite avantageusement estimée en fonction desdites mesures différentielles inter-station. De telles dispositions sont avantageuses en ce qu'elles permettent 25 d'annuler les éventuelles incertitudes sur les puissances d'émission du premier signal de référence et du second signal de référence par ledit satellite 20. En effet, la précision avec laquelle lesdites puissances d'émission peuvent être estimées peut s'avérer insuffisante et limiter les performances de l'estimation d'attitude du satellite 20. 30 Notamment, lesdites puissances d'émission sont susceptibles de fluctuer sensiblement au cours du temps, par exemple du fait des variations de température subies par le satellite 20. Par conséquent, la précision avec laquelle lesdites puissances d'émission peuvent être estimées est limitée par 3041103 18 l'amplitude de ces fluctuations. En outre, la puissance d'émission, par le satellite 20, du second signal de référence peut s'avérer difficile à estimer dans la mesure où elle dépend de plusieurs paramètres qui peuvent être eux-mêmes difficiles à estimer avec une 5 précision suffisante, dont : - la puissance avec laquelle le second signal de référence a été émis par la station d'émission de référence 30, - les gains respectifs de l'antenne d'émission de ladite station d'émission de référence 30 et de l'antenne de réception 22 du 10 satellite 20, - les pertes de propagation sur le lien montant. Par conséquent, le calcul des mesures différentielles inter-station permet d'améliorer la précision d'estimation de l'attitude du satellite 20. Par exemple, la puissance Pri[m] avec laquelle le premier signal de 15 référence est reçu par une station de mesure de rang m, et la puissance Pr2[m] avec laquelle le second signal de référence est reçu par ladite station de mesure de rang m, peuvent être exprimées sous la forme suivante en échelle logarithmique (dB) : Pri[m] = Pt1 + Gt1[m] + Gr[m] - Ld[m] 20 Pr2[m] = Pt2 + Gt2[m] + Gr[m] - Ld[m] expressions dans lesquelles : - Pt1 correspond à la puissance d'émission du premier signal de référence par le satellite 20, - Gt1[m] correspond au gain de la première antenne d'émission 21 25 dans la direction de la station de mesure de rang m, - Gr[m] correspond au gain de l'antenne de réception de la station de mesure de rang m dans la direction du satellite 20, - Ld[m] correspond aux pertes de propagation sur le lien descendant entre le satellite 20 et la station de mesure de rang m, 30 - Pt2 correspond à la puissance d'émission du second signal de référence par le satellite 20, - Gt2[m] correspond au gain de la seconde antenne d'émission 23 dans la direction de la station de mesure de rang m. 3041103 19 Les puissances Pri[m] et Pr2[m] sont par exemple mesurées dans le domaine fréquentiel, après avoir appliqué une fenêtre temporelle de type « flat-top » sur le signal reçu par la station de mesure 31. Une telle fenêtre temporelle permet de concentrer la puissance d'un signal non modulé sur une 5 fréquence, et est considérée comme optimale pour détecter et mesurer la puissance du premier signal de référence et du second signal de référence, en particulier dans le cas où ceux-ci sont des signaux sinusoïdaux. La mesure différentielle de puissance U[m] est par exemple calculée, pour la station de mesure de rang m, selon l'expression suivante : 10 EP[m] = Pt1 - Pt2 = Pt1 + Gt1[m] - Pt2 - Gt2[m] La mesure différentielle inter-station AP[m,n], entre la station de mesure de rang m et la station de mesure de rang n (m étant différent de n), est par exemple calculée selon l'expression suivante : AP[m,n] = EP[m] - EP[n] = Gt1[m] - Gt2[m] - (Gt1[n] - Gt2[n]) 15 On constate donc que les mesures différentielles inter-station ne dépendent plus que des différentes valeurs des gains de la première antenne d'émission 21 et de la seconde antenne d'émission 23 dans les directions respectives des stations de mesure de rang m et n. Par conséquent, il est possible de calculer des mesures différentielles 20 inter-station théoriques APT[m,n](Y) pour chaque attitude Y possible du satellite 20, en fonction des positions respectives du satellite 20 et des stations de mesure 31, et en fonction des diagrammes de rayonnement respectifs de la première antenne d'émission 21 et de la seconde antenne d'émission 23 par rapport audit satellite 20. 25 L'attitude du satellite 20 peut être estimée en optimisant une fonction de coût prédéterminée. Par exemple, l'attitude estimée YEST du satellite 20 est déterminée, en fonction des mesures différentielles inter-station OP[k] calculées pour Np paires de stations de mesure 31 (1 k Np), de sorte à minimiser l'erreur quadratique entre les mesures différentielles inter-station 30 A P[k] calculées et les mesures différentielles inter-station théoriques APT[k], par exemple pondérée selon l'expression suivante : 3041103 20 1\1,) { ( YEST = arg min 1 -,8,P[k] - APT [k](Y Y k=1 ak expression dans laquelle o-k correspond à la variance du bruit de la mesure différentielle inter-station AP[k]. Tel qu'indiqué précédemment, l'erreur d'attitude en lacet est 5 généralement négligeable pour un satellite 20 en orbite GEO. Par conséquent, il est possible de ne considérer que des attitudes Y possibles pour lesquelles l'attitude en lacet correspond à l'attitude de mission. Alternativement, il est possible de considérer plusieurs valeurs possibles de l'attitude en lacet lors du calcul des mesures différentielles inter-station théoriques A P-r[k]. Il est 10 également possible d'estimer l'attitude en lacet par d'autres moyens, par exemple au moyen de mesures fournies par des capteurs embarqués à bord du satellite 20 (l'attitude en lacet estimée étant par exemple émise par ledit satellite 20 jusqu'à la station centrale), et de ne calculer des mesures différentielles inter-station théoriques APT[k] que pour des attitudes Y possibles 15 pour lesquelles l'attitude en lacet est égale à cette estimation. Dans des modes préférés de mise en oeuvre, les stations de mesure 31 sont synchronisées temporellement entre elles. Par « synchronisées temporellement entre elles », on entend que le système 10 d'estimation est capable de recaler temporellement entre elles des mesures différentielles de 20 puissance calculées pour deux stations de mesure 31 différentes. De telles dispositions permettent par conséquent d'assurer une meilleure cohérence temporelle des mesures différentielles de puissance prises en compte pour calculer les mesures différentielles inter-station. En effet, il est alors possible d'identifier des mesures différentielles de puissance 25 calculées pour des signaux reçus en même temps par les différentes stations de mesure 31, éventuellement aux différences de temps de propagation entre le satellite 20 et les différentes stations de mesure 31 près. On assure alors que les signaux pris en compte dans les différentes stations de mesure 31 ont été émis avec sensiblement les mêmes puissances d'émission Pt1, pour le 30 premier signal de référence, et Pt2 pour le second signal de référence. Le procédé 50 d'estimation d'attitude ainsi obtenu est par conséquent 12 1 3041103 21 plus robuste aux variations temporelles des puissances d'émission par le satellite 20 du premier signal de référence et du second signal de référence. Les stations de mesure 31 sont de préférence synchronisées temporellement par GPS (« Global Positioning System »). Toutefois, toute 5 méthode de synchronisation temporelle peut être mise en oeuvre, et le choix d'une méthode particulière ne constitue qu'une variante d'implémentation de l'invention. La précision avec laquelle les stations de mesure 31 sont synchronisées temporellement entre elles est de préférence de l'ordre de 10 millisecondes (ms), voire inférieure. 10 Dans des modes préférés de mise en oeuvre, les stations de mesure 31 comportant des chaînes de réception respectives, le procédé 50 d'estimation d'attitude comporte, pour chaque station de mesure, l'estimation d'une erreur introduite par ladite chaîne de réception sur le calcul de la mesure différentielle de puissance, dite « erreur différentielle de puissance », et la 15 compensation de ladite erreur différentielle de puissance sur la mesure différentielle de puissance. En effet, les chaînes de réception des stations de mesure 31 sont constituées d'équipements dont les caractéristiques peuvent varier avec la fréquence. Ainsi, le gain introduit, par une chaîne de réception, sur le premier 20 signal de référence (reçu sur la fréquence centrale F1) peut être différent du gain introduit sur le second signal de référence (reçu sur la fréquence centrale F2), et ces différences de gain sont susceptibles d'introduire un biais sur les mesures différentielles de puissance et, in fine, sur les mesures différentielles inter-station. L'estimation des erreurs différentielles de puissance vise 25 précisément à calibrer ces différences de gain et à les compenser. Par exemple, l'estimation de l'erreur différentielle de puissance comporte, pour l'une au moins des stations de mesure 31, l'injection dans la chaîne de réception de ladite station de mesure d'un premier signal de calibration et d'un second signal de calibration. 30 De préférence, le premier signal de calibration et le second signal de calibration sont injectés le plus en amont possible de la chaîne de réception, afin que ceux-ci traversent essentiellement les mêmes équipements que le premier signal de référence et le second signal de référence. Par exemple, le 3041103 22 premier signal de calibration et le second signal de calibration sont présentés en entrée du LNB (« Low Noise Bock converter »), c'est-à-dire immédiatement en sortie de l'antenne de réception de la station de mesure 31. De préférence, le premier signal de calibration et le second signal de 5 calibration sont injectés dans la chaîne de réception simultanément à la réception du premier signal de référence et du second signal de référence. De telles dispositions sont avantageuses en ce que la différence de gain introduite par une chaîne de réception peut varier au cours du temps, et on assure de la sorte une meilleure cohérence temporelle entre l'erreur différentielle de 10 puissance estimée et la mesure différentielle de puissance calculée. Par exemple, il est possible d'injecter, au cours des opérations d'estimation d'attitude, le premier signal de calibration et le second signal de calibration de manière continue dans la chaîne de réception, ce qui permet d'assurer qu'ils seront bien reçus en même temps que le premier signal de référence et le 15 second signal de référence. Dans des modes préférés de mise en oeuvre, le premier signal de calibration et le second signal de calibration sont injectés sur des fréquences centrales respectives F'1 et F'2 différentes, et sont par exemple des signaux sinusoïdaux. Les fréquences centrales F'1 et F'2 sont de préférence choisies 20 proches des fréquences centrales respectivement F1 et F2, afin d'assurer que le premier signal de calibration et le second signal de calibration soient soumis sensiblement aux mêmes différences de gains que le premier signal de référence et le second signal de référence. Par exemple, l'écart entre les fréquences centrales F'1 et F1 (respectivement F'2 et F2) est inférieur à 100 25 kilohertz (kHz). Le premier signal de calibration et le second signal de calibration sont injectés, dans la chaîne de réception, avec la même puissance ou avec un rapport de puissances prédéterminé. L'erreur différentielle de puissance est ensuite calculée, comme la mesure différentielle de puissance, sous la forme 30 d'une grandeur représentative d'un rapport entre les puissances respectives du premier signal de calibration et du second signal de calibration mesurés en sortie de ladite chaîne de réception (éventuellement corrigé par ledit rapport des puissances d'injection dudit premier signal de calibration et dudit second 3041103 23 signal de calibration dans la chaîne de réception, si lesdites puissances d'injection ne sont pas égales). La figure 4 représente schématiquement un mode préféré de réalisation d'une station de mesure 31. 5 Tel qu'illustré par la figure 4, la station de mesure 31 comporte notamment une antenne de réception, par exemple du type comportant un réflecteur 310 et une source 311, en sortie de laquelle se trouve une chaîne de réception 312. L'entrée de la chaîne de réception 312 comporte un LNB 313, et la sortie de ladite chaîne de réception 312 est reliée à un module de traitement 10 314 comportant des moyens configurés pour numériser et traiter les signaux en sortie de ladite chaîne de réception 312. En outre, la station de mesure 31 comporte un module de génération 315 du premier signal de calibration et du second signal de calibration dont la sortie est reliée à un coupleur 316 en entrée du LNB 313. 15 Plus particulièrement, dans l'exemple illustré par la figure 4, le module de génération 315 comporte un premier générateur 317 d'un signal sur une fréquence centrale Fg1 et un second générateur 318 d'un signal sur une fréquence centrale Fg2. Les fréquences centrales Fg1 et Fg2 sont avantageusement choisies de telle sorte que : 20 - si F'1 F'2 : IFg1 - Fg21= F'2 et IFg1 + Fg21= F'1. Le premier générateur 317 et le second générateur 318 sont reliés à un mélangeur non linéaire 319 dont la sortie est reliée au coupleur 316. De manière connue, un tel mélangeur non linéaire 319 (« mixer » dans la 25 littérature anglo-saxonne) réalise le produit des signaux en entrée. Ce produit fournit en sortie deux signaux de même amplitude sur des fréquences centrales respectives égales à IFg1 - Fg2I et IFg1 + Fg21. Par conséquent, les signaux en sortie du mélangeur non linéaire 319 sont centrés sur les fréquences centrales F'1 et F'2, et correspondent respectivement au premier 30 signal de calibration et au second signal de calibration. De telles dispositions sont avantageuses en ce qu'elles permettent d'assurer un meilleur contrôle des puissances d'injection du premier signal de calibration et du second signal de calibration dans la chaîne de réception 312. 3041103 24 En effet, les signaux générés par le premier générateur 317 et le second générateur 318 sont susceptibles de subir des atténuations différentes, du fait notamment qu'ils sont centrés sur des fréquences centrales respectives Fg1 et Fg2 différentes et que les moyens de connexion jusqu'au mélangeur non 5 linéaire 319 peuvent présenter une longueur importante. Dans le cas illustré par la figure 4, le mélangeur non linéaire 319 peut être agencé le plus près possible du coupleur 316, afin de limiter la longueur des moyens de connexion entre ledit mélangeur non linéaire 319 et ledit coupleur 316. Grâce au mélangeur non linéaire 319, le premier signal de calibration et le second signal 10 de calibration sont injectés avec sensiblement la même puissance, et ce quelles que soient les puissances respectives des signaux en entrée dudit mélangeur non linéaire 319. De manière plus générale, il est à noter que les modes de mise en oeuvre et de réalisation considérés ci-dessus ont été décrits à titre d'exemples 15 non limitatifs, et que d'autres variantes sont par conséquent envisageables. En particulier, l'invention a été décrite en considérant une estimation d'attitude d'un satellite 20 basée principalement sur des mesures différentielles de puissance entre des signaux émis par ledit satellite 20 et reçus par différentes stations de mesure 31. Rien n'exclut, suivant d'autres exemples, de 20 considérer également d'autres mesures. Notamment, l'estimation d'attitude peut également prendre en compte les mesures des capteurs à bord du satellite 20. Par exemple, les mesures desdits capteurs peuvent être transmises au sol et le procédé 50 peut estimer l'attitude du satellite 20 en fonction à la fois des mesures différentielles de 25 puissance et des mesures desdits capteurs. Il est également possible, notamment, d'estimer au sol l'attitude du satellite 20 en fonction des mesures différentielles de puissance, et d'envoyer l'attitude estimée au satellite 20. Le satellite 20 peut alors à son tour estimer son attitude en fonction des mesures des capteurs et en fonction de l'attitude estimée reçue du sol. 30 REVENDICATIONS1 Procédé (50) d'estimation d'attitude d'un satellite (20), ledit satellite étant configuré pour générer et pour émettre un premier signal de référence sur un lien descendant au moyen d'une première antenne d'émission (21), ledit satellite étant en outre configuré pour recevoir des signaux utiles sur un lien montant, et pour retransmettre lesdits signaux utiles sur le lien descendant au moyen d'une seconde antenne d'émission (23), la première antenne d'émission desservant une zone de couverture (Zc1) englobant une zone de couverture (Zc2) de la seconde antenne d'émission, caractérisé en ce qu'il comporte : - l'émission (51) sur le lien montant d'un second signal de référence, qui est reçu par ledit satellite au moyen d'une antenne de réception (22) des signaux utiles, - la retransmission (52) par ledit satellite du second signal de référence au moyen de la seconde antenne d'émission, - la réception (54) du premier signal de référence et du second signal de référence par au moins trois stations sol se trouvant dans la zone de couverture de la seconde antenne d'émission, dites « stations de mesure » (31), - le calcul (55), pour chaque station de mesure, d'une grandeur représentative d'un rapport entre les puissances respectives du premier signal de référence et du second signal de référence reçus, dite « mesure différentielle de puissance », - l'estimation (56) de l'attitude du satellite en fonction des mesures différentielles de puissance calculées. 2 - Procédé (50) selon la revendication 1, comportant le calcul, pour au moins deux paires de stations de mesure, d'une grandeur représentative d'un rapport entre les mesures différentielles de puissance respectives des stations de mesure de ladite paire, dite « mesure différentielle inter- station », et dans lequel l'attitude du satellite (20) est estimée en fonction desdites mesures différentielles inter-station calculées. 3 Procédé (50) selon la revendication 2, comportant la synchronisation temporelle des stations de mesure entre elles, et dans lequel les mesures 3041103 26 différentielles inter-station sont calculées à partir de mesures différentielles de puissance calculées pour des signaux reçus en même temps par les stations de mesure. 4 Procédé (50) selon l'une des revendications précédentes, caractérisé en 5 ce que, les stations de mesure comportant des chaînes de réception respectives, ledit procédé comporte, pour chaque station de mesure, l'estimation d'une erreur introduite par ladite chaîne de réception sur le calcul de la mesure différentielle de puissance, dite « erreur différentielle de puissance », et la compensation de ladite erreur différentielle de 10 puissance sur la mesure différentielle de puissance. 5 Procédé (50) selon la revendication 4, dans lequel l'estimation de l'erreur différentielle de puissance comporte, pour l'une au moins des stations de mesure, l'injection dans la chaîne de réception de ladite station de mesure d'un premier signal de calibration et d'un second signal de calibration, et le 15 calcul d'une grandeur représentative d'un rapport entre les puissances respectives du premier signal de calibration et du second signal de calibration mesurés en sortie de ladite chaîne de réception. 6 Procédé (50) selon la revendication 5, dans lequel le premier signal de calibration et le second signal de calibration sont injectés sur des 20 fréquences centrales respectives F'1 et F'2 différentes. 7 Procédé (50) selon la revendication 6, comportant la génération du premier signal de calibration et du second signal de calibration par mélange non linéaire de deux signaux de fréquences centrales respectives Fg1 et Fg2 telles que : 25 - IFg1 - Fg2I = F'1 et IFg1 + Fg2I = F'2, ou - IFg1 - Fg21= F'2 et IFg1 + Fg21= F'1. 8 Procédé (50) selon l'une des revendications précédentes, dans lequel le premier signal de référence et le second signal de référence sont émis sur des fréquences centrales respectives F1 et F2 différentes sur le lien 30 descendant. 9 - Système (10) d'estimation d'attitude d'un satellite (20), ledit satellite comportant des moyens configurés pour émettre un premier signal de référence sur un lien descendant au moyen d'une première antenne 3041103 27 d'émission (21), recevoir des signaux utiles sur un lien montant et retransmettre lesdits signaux utiles sur le lien descendant au moyen d'une seconde antenne d'émission (23), ladite première antenne d'émission desservant une zone de couverture (Zc1) englobant une zone de 5 couverture (Zc2) de ladite seconde antenne d'émission, caractérisé en ce qu'il comporte : - une station sol, dite « station d'émission de référence » (30), configurée pour émettre un second signal de référence sur le lien montant, adapté à être reçu par le satellite avec les signaux utiles, 10 - au moins trois stations sol adaptées à recevoir le premier signal de référence émis au moyen de la première antenne d'émission et le second signal de référence retransmis au moyen de la seconde antenne d'émission, dites « stations de mesure » (31), - des moyens configurés pour calculer, pour chaque station de 15 mesure, une grandeur représentative d'un rapport entre les puissances respectives du premier signal de référence et du second signal de référence reçus, dite « mesure différentielle de puissance », - des moyens configurés pour estimer l'attitude du satellite en 20 fonction des mesures différentielles de puissance calculées. 10 - Système (10) selon la revendication 9, comportant des moyens configurés pour calculer, pour au moins deux paires de stations de mesure, une grandeur représentative d'un rapport entre les mesures différentielles de puissance respectives des stations de mesure de ladite paire, dite 25 « mesure différentielle inter-station », et dans lequel l'attitude du satellite (20) est estimée en fonction desdites mesures différentielles inter-station calculées. 11 - Système (10) selon la revendication 10, dans lequel les stations de mesure sont synchronisées temporellement entre elles. 30 12 - Système (10) selon l'une des revendications 9 à 11, dans lequel les trois stations de mesure considérées sont situées, dans la zone de couverture de la seconde antenne d'émission du satellite placé dans une attitude de mission, de telle sorte que la valeur du gradient du diagramme de 3041103 28 rayonnement de ladite seconde antenne d'émission dans la direction d'une desdites trois stations de mesure, dite « station de mesure de référence » (31a), est inférieure aux valeurs dudit gradient du diagramme de rayonnement dans les directions respectives des deux autres stations 5 de mesure (31b, 31c). 13 - Système (10) selon l'une des revendications 9 à 12, dans lequel, chaque station de mesure comportant une chaîne de réception, ledit système comporte des moyens configurés pour estimer, pour chaque station de mesure, une erreur introduite par ladite chaîne de réception sur le calcul 10 de la mesure différentielle de puissance, dite « erreur différentielle de puissance », et pour compenser ladite erreur différentielle de puissance sur la mesure différentielle de puissance. 14 - Système (10) selon la revendication 13, dans lequel l'une au moins des stations de mesure comporte des moyens configurés pour générer et pour 15 injecter un premier signal de calibration et un second signal de calibration dans la chaîne de réception de ladite station de mesure, et des moyens configurés pour calculer une grandeur représentative d'un rapport entre les puissances respectives du premier signal de calibration et du second signal de calibration mesurés en sortie de ladite chaîne de réception. 20 15 - Système (10) selon la revendication 14, dans lequel, le premier signal de calibration et le second signal de calibration étant injectés sur des fréquences centrales respectives F'1 et F'2 différentes, les moyens de génération dudit premier signal de calibration et dudit second signal de calibration sont configurés pour mélanger de manière non linéaire deux 25 signaux de fréquences centrales respectives Fg1 et Fg2 telles que : - IFg1 - Fg2I = F'1 et IFg1 + Fg2I = F'2, ou - IFg1 - Fg21= F'2 et IFg1 + Fg21= F'1.