FR 2988860 A1 20131004 FR 1200948 A 20120329 DOMAINE TECHNIQUE La présente invention concerne l'estimation d'attitude d'un engin aérien ou spatial en fonction de mesures de phase de signaux radioélectriques reçus par cet engin. Plus particulièrement, la présente invention concerne la calibration d'erreurs, dites « erreurs spatiales », induites par des biais de phase affectant les mesures de phase, lesdits biais de phase dépendant d'une direction d'incidence des signaux radioélectriques par rapport à l'engin. ÉTAT DE LA TECHNIQUE Il est en effet connu, pour estimer l'attitude d'un engin spatial de type satellite, de mesurer les phases de signaux radioélectriques reçus d'un émetteur par deux bases de réception de l'engin, une base de réception consistant en deux antennes de réception, et les antennes de réception desdites deux bases de réception n'étant pas toutes alignées. En effet, il est possible de déterminer une différence de marche des signaux radioélectriques reçus de l'émetteur par les deux antennes de réception d'une base de réception en comparant les phases mesurées sur lesdites deux antennes de réception. La direction d'incidence de l'émetteur par rapport à l'engin, c'est-à-dire la direction de l'émetteur dans un repère fixe par rapport à l'engin, peut ensuite être déterminée en fonction de deux différences de marche déterminées pour les deux bases de réception. En estimant les directions d'incidence de deux émetteurs, il est possible d'estimer l'attitude de l'engin relativement à ces émetteurs. L'attitude absolue de l'engin peut en outre être estimée si les positions absolues respectives de l'engin et des émetteurs sont connues. De manière générale, la précision de l'estimation de l'attitude de l'engin est améliorée en considérant un nombre d'émetteurs supérieur à deux. En considérant des émetteurs tels que des satellites d'un système global de navigation par satellites (GPS, Galileo, etc.), l'estimation de l'attitude de l'engin sera grandement facilitée. En effet, des satellites GPS, par exemple, émettent des informations de navigation qui permettent de déterminer la position de l'engin à partir d'informations de navigation reçues par des récepteurs GPS montés à bord dudit engin. Toutefois, l'estimation d'attitude en fonction de mesures de phase repose généralement sur l'hypothèse que les signaux radioélectriques reçus ont suivi un trajet unique direct depuis l'émetteur jusqu'au centre de phase de chacune des antennes de réception. Or, il existe souvent d'autres trajets, indirects, entre l'émetteur et chacune des antennes de réception. Ces trajets indirects sont liés à la présence d'éléments réfléchissants et/ou diffractants à proximité des antennes de réception de l'engin. Ainsi, un signal radioélectrique reçu de l'émetteur par une antenne de réception sera la somme d'une composante reçue via le trajet direct, dite « composante directe », et d'une composante reçue via des trajets indirects, dite « composante indirecte ». 10 La composante indirecte mesurée pour une antenne de réception peut généralement être assimilée à un biais de phase dont la valeur dépend essentiellement de la direction d'incidence du signal radioélectrique reçu via le trajet direct, et éventuellement de la position et de l'orientation des éléments réfléchissants et/ou diffractants lorsque ceux-ci sont mobiles. 15 Toutefois, pour une même direction d'incidence, les biais de phase dus aux multi-trajets varient d'une antenne de réception à une autre, de sorte que ces biais de phase ne se compensent pas dans le calcul des différences de phases entre paires d'antennes de réception et introduisent des erreurs sur les différences de phases utilisées pour l'estimation de l'attitude de l'engin. Ces 20 erreurs sur les différences de phase sont dites « erreurs spatiales » car les biais de phase qui en sont à l'origine dépendent de la direction d'incidence considérée. Pour calibrer les erreurs spatiales, il est possible de comparer l'attitude, estimée en fonction des différences de phases de signaux 25 radioélectriques reçus par les trois antennes de réception, à une autre estimation de l'attitude fournie par un senseur d'attitude, distinct des trois antennes de réception, tel qu'un senseur stellaire. On comprend cependant qu'une telle solution nécessite d'embarquer dans l'engin deux systèmes distincts pour estimer l'attitude de l'engin. 30 EXPOSÉ DE L'INVENTION La présente invention a pour objectif de remédier à tout ou partie des limitations des solutions de l'art antérieur, notamment celle exposées ci-avant, en proposant une solution qui permette de calibrer les erreurs spatiales, induites par la réception sur les antennes de réception de l'engin de composantes indirectes, qui ne nécessite pas d'utiliser un senseur d'attitude distinct desdites antennes de réception. A cet effet, et selon un premier aspect, l'invention concerne un 5 procédé de calibration d'erreurs spatiales induites par des biais de phase affectant des mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception non alignées d'un engin, lesdits biais de phase dépendant d'une direction d'incidence des signaux radioélectriques par rapport à l'engin. Le procédé de calibration comporte des 10 étapes de : - obtention de positions respectives de l'engin et d'au moins une paire de satellites se trouvant dans un champ de vue radioélectrique de l'engin, estimation de l'écart angulaire inter-satellites de la paire de 15 satellites en fonction des positions respectives de l'engin et desdits satellites de la paire, de sorte à obtenir un écart angulaire inter-satellites dit « théorique », - estimation de directions d'incidence respectives des satellites de la paire par rapport à l'engin en fonction de mesures de phase de 20 signaux radioélectriques reçus desdits satellites de la paire par les trois antennes de réception de l'engin, - estimation de l'écart angulaire inter-satellite de la paire de satellites en fonction des directions d'incidence respectives des satellites de la paire par rapport à l'engin, de sorte à obtenir un écart angulaire 25 inter-satellites dit « estimé », estimation des erreurs spatiales pour chacune des directions d'incidence respectives des satellites de la paire par rapport à l'engin en fonction de l'écart angulaire inter-satellites théorique et de l'écart angulaire inter-satellites estimé déterminé pour ladite 30 paire de satellites. En pratique, plusieurs paires de satellites seront considérées pour la calibration des erreurs spatiales, lesdites paires de satellites pouvant par exemple être considérées simultanément et/ou successivement. Par « écart angulaire inter-satellites » d'une paire de satellites, on entend l'angle entre les directions respectives des deux satellites de cette paire par rapport à l'engin. L'intérêt d'estimer des écarts angulaires inter-satellites réside dans le fait que ceux-ci sont indépendants de l'attitude de l'engin, mais 5 dépendent uniquement des positions respectives des satellites et de l'engin. La calibration des erreurs spatiales est donc réalisée en estimant les écarts angulaires inter-satellites de paires de satellites de deux façons différentes : en fonction des positions respectives des satellites et de l'engin, de 10 sorte à obtenir des écarts angulaires inter-satellites théoriques, en fonction des directions d'incidence respectives des satellites, estimées en fonction de mesures de phase affectées par les biais de phase à calibrer, de sorte à obtenir des écarts angulaires inter-satellites estimés. 15 En considérant les écarts angulaires inter-satellites théoriques comme dépourvus d'erreurs, et en comparant les écarts angulaires inter-satellites estimés, affectés par les biais de phase, auxdits écarts angulaires inter-satellites théoriques, il sera possible de calibrer les erreurs spatiales dans la limite de la précision des positions respectives des satellites et de l'engin. 20 Il est cependant à noter que, du fait que les écarts angulaires inter- satellites considérés ne dépendent pas de l'attitude de l'engin, les erreurs spatiales pourront dans certains cas n'être calibrées qu'à un biais d'attitude près. En d'autres termes, l'attitude de l'engin, éventuellement estimée après calibration et correction des erreurs spatiales, pourra être affectée d'un biais. 25 Toutefois, ce biais d'attitude sera le même quelle que soit l'attitude réelle de l'engin et pourra, le cas échéant, être aisément détecté et à son tour calibré en mettant en oeuvre des moyens considérés comme connus de l'homme de l'art. Dans des modes particuliers de mise en oeuvre, le procédé de calibration comporte l'une ou plusieurs des caractéristiques suivantes, prises 30 isolément ou suivant toutes les combinaisons techniquement possibles. Dans un mode particulier de mise en oeuvre, le champ de vue radioélectrique de l'engin est préalablement partitionné en une pluralité de cônes distincts d'incidence de signaux radioélectriques par rapport à l'engin, et les erreurs spatiales sont considérées comme étant les mêmes pour toute direction d'incidence à l'intérieur d'un même cône d'incidence. De telles dispositions permettent de réduire le nombre de directions d'incidence différentes pour lesquelles des erreurs spatiales doivent être 5 calibrées, et donc de réduire le nombre de calculs à effectuer ainsi que la quantité d'erreurs spatiales à mémoriser. Dans un mode particulier de mise en oeuvre, les cônes d'incidence d'élévation faible par rapport à un plan déterminé par les centres de phase des trois antennes de réception de l'engin sont d'angles solides respectifs inférieurs 10 à ceux de cônes d'incidence d'élévation élevée par rapport audit plan. Il a en effet été constaté que la composante indirecte est plus faible pour des directions d'incidence sensiblement normales (élévation élevée), par rapport au plan déterminé par les trois antennes de réception, que pour les directions d'incidence rasantes (élévation faible) par rapport audit plan. 15 De telles dispositions permettent donc d'optimiser la calibration des erreurs spatiales en estimant davantage d'erreurs spatiales pour les directions d'incidence pour lesquelles les erreurs spatiales varient beaucoup que pour les directions d'incidence pour lesquelles les erreurs spatiales varient peu. Dans un mode particulier de mise en oeuvre, les satellites 20 appartiennent à un système global de navigation par satellites et les trois antennes de réception sont des antennes de récepteurs dudit système global de navigation par satellites. En outre, les positions respectives des satellites et de l'engin sont déterminées à partir d'informations de navigation incluses dans les signaux radioélectriques reçus des satellites. 25 De telles dispositions sont avantageuses en ce que la calibration des erreurs spatiales peut être réalisée exclusivement à partir des signaux radioélectriques reçus des satellites : - les écarts angulaires inter-satellites théoriques sont déterminés à partir des informations de navigation incluses dans les signaux 30 radioélectriques reçus des satellites, - les écarts angulaires inter-satellites estimés sont déterminés à partir des mesures de phase desdits signaux radioélectriques reçus des satellites. La calibration des erreurs spatiales peut ainsi être réalisée sans utiliser d'autre senseur d'attitude que les trois antennes de réception. Il est à noter que, dans ce mode particulier de mise en oeuvre, les écarts angulaires inter-satellites peuvent être considérés comme dépourvus d'erreurs. En effet, les positions de tels satellites sont couramment estimées à quelques mètres près. D'autre part, la position de l'engin peut être estimée à quelques centaines de mètres près par les récepteurs du système global de navigation par satellites. Du fait de l'éloignement entre les satellites et l'engin qui, dans le cas d'un engin en orbite basse, est très supérieur à 1000 kilomètres, il s'en suit que l'estimation des écarts angulaires inter-satellites théoriques à partir des positions des satellites et de l'engin est extrêmement précise et peut être considérée au premier ordre comme dépourvue d'erreurs. Dans un mode particulier de mise en oeuvre, les étapes du procédé sont itérées pour chaque paire de satellites considérée. Dans un mode particulier de mise en oeuvre, lorsque Ns satellites organisés en Np paires se trouvent dans le champ de vue radioélectrique de l'engin, les erreurs spatiales pour chacune des directions d'incidence respectives desdits Ns satellites par rapport à l'engin sont estimées simultanément en fonction des écarts angulaires inter-satellites théoriques et des écarts angulaires inter-satellites estimés déterminés pour lesdites Ns paires de satellites. Selon un second aspect, l'invention concerne un procédé d'estimation d'attitude d'un engin en fonction de mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception non alignées de l'engin. Ledit procédé d'estimation d'attitude comporte également la calibration d'erreurs spatiales conformément à l'invention. Selon un troisième aspect, l'invention concerne un produit programme d'ordinateur comportant un ensemble d'instructions de code de programme qui, lorsqu'elles sont exécutées par un ou plusieurs processeur, mettent en oeuvre un procédé de calibration d'erreurs spatiales ou un procédé d'estimation d'attitude conformément à l'invention. Selon un quatrième aspect, l'invention concerne un système d'estimation d'attitude d'un engin en fonction de mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception non alignées de l'engin. Ledit système d'estimation d'attitude comporte en outre des moyens configurés pour calibrer les erreurs spatiales conformément à l'invention. Dans un mode particulier de réalisation du système d'estimation d'attitude, les trois antennes de réception sont des antennes de récepteurs d'un système global de navigation par satellites, par exemple des récepteurs GPS ou des récepteurs Galileo. Selon un cinquième aspect, l'invention concerne un engin spatial 10 comportant un système d'estimation d'attitude conforme à l'invention. Dans un mode particulier de réalisation, l'engin spatial comporte un corps comportant plusieurs faces externes sensiblement planes, et les trois antennes de réception sont agencées sur une même face externe dudit corps. PRÉSENTATION DES FIGURES 15 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 engin spatial et de satellites en orbite autour de la Terre, 20 Figure 2 : une représentation schématique d'une vue rapprochée de l'engin spatial de la figure 1, Figure 3: un diagramme illustrant les principales étapes d'un procédé de calibration d'erreurs spatiales, Figure 4 : une représentation schématique d'une vue rapprochée 25 de l'engin spatial illustrant l'expression d'une direction d'un satellite par rapport audit engin spatial, dans un repère fixe par rapport audit engin spatial, Figure 5: une représentation schématique d'un champ de vue radioélectrique de l'engin spatial, 30 Figure 6: une représentation schématique d'une partition particulière du champ de vue radioélectrique de l'engin spatial. 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 engin 10 spatial, dont on cherche à estimer l'attitude, et des satellites 20 en orbite autour de la Terre. Par exemple, l'engin 10 spatial est en orbite basse (LEO pour « Low Earth Orbit )) autour de la Terre. Dans la suite de la description, on se place de manière non limitative dans le cas où les satellites 20 sont des satellites d'un système global de navigation par satellites (GNSS pour « Global Navigation Satellite System »), 10 tel qu'un système GPS et/ou Galileo, etc. Plus particulièrement, on considère de manière non limitative des satellites 20 GPS. De tels satellites 20 GPS sont en orbite circulaire autour de la Terre, à une altitude approximative de 20000 kilomètres. La figure 1 représente également les directions respectives de deux 15 satellites 20 GPS par rapport à l'engin 10 spatial, lesquelles sont représentées par des vecteurs respectivement um et un. L'angle emn entre lesdites directions un, et un correspond à l'écart angulaire inter-satellites de la paire formée par lesdits deux satellites 20 GPS. La figure 2 représente schématiquement une vue rapprochée de 20 l'engin 10 spatial. Tel qu'illustré par cette figure, l'engin 10 spatial comporte un corps 12 comportant plusieurs faces externes sensiblement planes. Par exemple, le corps 12 est sensiblement en forme de cube. L'engin 10 spatial comporte également trois antennes de réception non alignées, respectivement A1, A2 et A3, adaptées à recevoir des signaux 25 radioélectriques émis par les satellites 20 GPS. Rien n'exclut, suivant d'autres exemples, de considérer un nombre plus important d'antennes de réception A1, A2, A3 pour estimer l'attitude dudit engin 10 spatial. Toutefois, il est nécessaire d'avoir au moins trois antennes de réception non alignées pour pouvoir estimer l'attitude dudit engin 10 spatial suivant trois axes. 30 Les trois antennes de réception A1, A2, A3 appartiennent à des récepteurs GPS. Ces récepteurs GPS peuvent être complètement distincts ou partager certains moyens, comme par exemple des moyens de calcul adaptés à déterminer la position de l'engin 10 en fonction des signaux radioélectriques reçus sur l'une quelconque desdites trois antennes de réception A1, A2, A3. De préférence, et tel qu'illustré par la figure 2, les trois antennes de réception A1, A2, A3 sont agencées sur une même face externe dudit corps. De la sorte, il est possible de maximiser un champ de vue radioélectrique 5 commun auxdites trois antennes de réception. Ceci permet de maximiser le nombre de satellites 20 GPS se trouvant simultanément dans les champs de vues radioélectriques desdites trois antennes de réception. Dans la suite, on désigne par « champ de vue radioélectrique de l'engin spatial » le champ de vue radioélectrique commun aux trois antennes de réception A1, A2, A3. 10 Les trois antennes de réception A1, A2, A3 sont organisées en deux bases de réception : une première base de réception B1 formée par les antennes de réception Al et A2, et une seconde base de réception B2 formée par les antennes de réception A2 et A3. Il est à noter que rien n'exclut, suivant d'autres exemples, de considérer également une troisième base de réception 15 formée par les antennes de réception Al et A3. De manière nullement limitative, on se place dans le cas où les deux bases de réception B1, B2 sont sensiblement orthogonales. Dans tout le contexte de la présente demande, on désigne par « direction d'incidence » d'un satellite 20 GPS l'expression de la direction de ce 20 satellite GPS par rapport à l'engin 10 spatial dans un repère fixe par rapport audit engin spatial. Les directions un, et un peuvent en effet être exprimées dans un repère quelconque ; lorsqu'elles sont exprimées dans un repère fixe par rapport à l'engin 10 spatial, elles sont désignées respectivement par vm et vn. On comprend que seules les directions d'incidence vm et vn des satellites 25 20 GPS par rapport à l'engin 10 spatial englobent une information sur l'orientation de l'engin 10 spatial par rapport aux satellites 20 GPS, information qui permet, combinée aux positions absolues respectives de l'engin 10 et des satellites 20 GPS, d'estimer l'attitude de l'engin. La figure 2 illustre un exemple de repère fixe par rapport à l'engin 10 30 spatial, dit « repère engin », dans lequel les directions d'incidence des différents satellites 20 GPS seront exprimées dans la suite de la description. Tel qu'illustré par la figure 2, le repère engin est sensiblement centré sur l'antenne de réception A2, et est défini par trois vecteurs unitaires x, y et z orthogonaux entre eux. Le vecteur x est sensiblement colinéaire avec les centres de phase des antennes de réception A1, A2 de la base de réception B1, et le vecteur y est sensiblement colinéaire avec les centres de phase des antennes de réception A2, A3 de la base de réception B2. Tel qu'indiqué précédemment, la présence d'éléments réfléchissants et/ou diffractants à proximité des trois antennes de réception A1, A2, A3 fait que les signaux radioélectriques reçus sur chacune desdites trois antennes de réception présenteront des composantes indirectes plus ou moins importantes. Lesdites composantes indirectes sont à l'origine de biais de phase qui dépendent d'une direction d'incidence des signaux radioélectriques par rapport à l'engin 10 spatial, et lesdits biais de phase introduisent des erreurs spatiales dans l'estimation de l'attitude dudit engin 10 spatial. La figure 3 représente les principales étapes d'un procédé 50 de calibration de ces erreurs spatiales, lesquelles sont : 51 obtention de positions respectives de l'engin 10 spatial et de Ns satellites 20 GPS se trouvant dans le champ de vue radioélectrique de l'engin spatial, les Ns satellites GPS étant organisés en Np paires, 52 estimation d'écarts angulaires inter-satellites respectifs des Np paires de satellites 20 GPS en fonction des positions respectives de l'engin 10 et des Ns satellites GPS, de sorte à obtenir des écarts angulaires inter-satellites dits « théoriques », 53 estimation de directions d'incidence respectives des Ns satellites 20 GPS par rapport à l'engin 10 spatial en fonction de mesures de phase de signaux radioélectriques reçus desdits satellites GPS par les trois antennes de réception A1, A2, A3 de l'engin 10 spatial, 54 estimation des écarts angulaires inter-satellites respectifs des Np paires de satellites 20 GPS en fonction des directions d'incidence respectives des Ns satellites GPS par rapport à l'engin 10 spatial, de sorte à obtenir des écarts angulaires inter-satellites dits « estimés », 55 estimation des erreurs spatiales pour chacune des directions d'incidence respectives des Ns satellites 20 GPS par rapport à l'engin 10 spatial en fonction desdits écarts angulaires inter-satellites théoriques et desdits écarts angulaires inter-satellites estimés. Ainsi, la calibration des erreurs spatiales repose sur l'estimation des écarts angulaires inter-satellites respectifs des Np paires de satellites 20 GPS de deux façons différentes de sorte à obtenir : des écarts angulaires inter-satellites théoriques indépendants des erreurs spatiales, des écarts angulaires inter-satellites estimés affectés par lesdites erreurs spatiales. Obtention des positions des satellites GPS et de l'engin spatial De manière connue, les satellites 20 GPS incorporent des informations de navigation dans les signaux radioélectriques qu'ils émettent. Lesdites informations de navigation comprennent notamment, pour chaque signal radioélectrique émis par un satellite 20 GPS, la position de ce satellite 20 GPS ainsi que l'heure d'émission dudit signal radioélectrique. Ainsi, l'engin 10 spatial, équipé d'un récepteur GPS, obtient directement les positions respectives des Ns satellites 20 GPS dans lesdites 20 informations de navigation, et peut déterminer de manière conventionnelle sa position à partir desdites informations de navigation. Détermination des écarts angulaires inter-satellites théoriques Dans la suite de la description, on suppose, afin de simplifier les équations, que les positions respectives des Ns satellites 20 GPS et de l'engin 25 10 spatial, en fonction desquelles les écarts angulaires inter-satellites théoriques sont déterminés, sont obtenues sans erreurs. En pratique, il a été vérifié que la précision desdites positions permet d'avoir en théorie une erreur sur l'estimation des écarts angulaires inter-satellites de l'ordre de 10e-04 degrés, ce qui est négligeable par rapport à la 30 précision visée pour l'estimation de l'attitude de l'engin 10 spatial, qui est de l'ordre de 0.1 degrés. Tel qu'illustré par la figure 2, l'écart angulaire inter-satellites entre deux satellites de directions respectivement um et un correspond à un angle emn entre lesdites directions um et u,'. Par exemple, les directions um et un sont estimées en fonction des positions des satellites 20 GPS considérés et de l'engin 10 spatial. Dans la suite de la description, on se place de manière non limitative dans le cas où l'écart angulaire inter-satellites emn est estimé sous la forme du produit scalaire um.un, ce qui revient à estimer en pratique cos(emn) lorsque les directions um et un sont des vecteurs unitaires. Les positions respectives de l'engin 10 spatial et des satellites 20 GPS sont par exemple obtenues dans un repère fixe centré sur la Terre. Dans ce repère fixe centré sur la Terre, on désigne par : - (xR, yR, zR) les coordonnées de l'engin 10 spatial, - (xm, ym, zm) les coordonnées du satellite 20 GPS de direction um, - (xn, yn, zn) les coordonnées du satellite 20 GPS de direction un. Dans ce cas, les directions um et un sont par exemple déterminées selon les expressions suivantes : 1 Ym YR, Zm - ZR) U = - - XR, m pm Un = 1 pn - ( xr., - xR, Yn - YR, zn zR) expressions dans lesquelles Pm et P,, sont des coefficients de normalisation. Puis l'écart angulaire inter-satellites théorique entre les satellites 20 20 GPS de directions un, et un est déterminé sous la forme du produit scalaire um.un qui est égal à : (xm xR) - (xn xR) + (Ym YR)- (yn - YR) + (Zm zR) - (zn zR) um un = Pm - Pn Estimation des directions d'incidence en fonction de mesures de phase Tel qu'indiqué précédemment, on entend par « estimer la direction 25 d'incidence » d'un satellite 20 GPS le fait d'estimer la direction de ce satellite GPS par rapport à l'engin 10 spatial dans un repère fixe par rapport audit engin spatial, tel que le repère engin. La figure 4 représente de manière plus détaillée la direction d'incidence vm dans le repère engin, afin d'expliquer comment ladite direction 30 d'incidence vm peut être estimée. Tel qu'illustré par la figure 4, la direction d'incidence vm fait un angle al m avec le vecteur x du repère engin, et un angle a2m avec le vecteur y dudit repère engin. On se place en outre de manière nullement limitative dans le cas où les deux bases de réception B1, B2 sont de même longueur b (non illustrée par les figures). Par conséquent, on a cos(al m) = Arl m/b et cos(a2m) = Ar2m/b, expressions dans lesquelles : Art m est la différence de marche des signaux radioélectriques reçus dans la direction d'incidence vm par les antennes de réception A1, A2 de la base de réception B1, Ar2m est la différence de marche des signaux radioélectriques reçus dans la direction d'incidence vm par les antennes de réception A2, A3 de la base de réception B2. Il est à noter que, dans ces expressions, Arl m (respectivement Ar2m) 15 est négatif dès lors que al m (respectivement a2m) est supérieur à Tr/2. La direction d'incidence vm s'exprime par exemple : ( Arim Ar2m b b (Arlm12 (Ar2m) b ) b Dans cette expression, la longueur b est connue. En outre, et de manière connue de l'homme de l'art, il est possible de 20 déterminer une différence de marche des signaux radioélectriques reçus par les deux antennes de réception d'une base de réception B1, B2 en fonction de la différence des phases mesurées sur lesdites deux antennes de réception. Il est à noter qu'une ambiguïté peut se produire du fait que, la phase étant mesurée modulo 2u, la différence de marche est estimée modulo la 25 longueur d'onde des signaux radioélectriques reçus. Toutefois, il existe des procédés, considérés comme connus de l'homme de l'art, pour lever l'ambiguïté sur la différence de marche (voir par exemple FR 2926891). Par conséquent, une estimation vem de la direction d'incidence vm peut être obtenue en estimant les différences de marche Art m et Ar2m par des 30 procédés considérés comme connus de l'homme de l'art. Détermination des écarts angulaires inter-satellites estimés Vm = Les écarts angulaires inter-satellites estimés sont déterminés en fonction des estimations vek (1 É k Ns) des directions d'incidence vk des Ns satellites 20 GPS par rapport à l'engin 10 spatial. Par exemple, des estimations vem et ven des directions d'incidence 5 respectivement vm et vn sont déterminées, par exemple tel que décrit précédemment, et l'écart angulaire inter-satellite estimé correspondant est déterminé sous la forme du produit scalaire vem-ven. Estimation des erreurs spatiales Dans la suite de la description, on se place dans le cas où les erreurs 10 spatiales sont estimées, pour chaque base de réception B1, B2, sous la forme d'une erreur sur une différence de marche des signaux radioélectriques. En effet, bien que les biais de phase varient d'une antenne de réception à une autre, c'est bien la différence de phase entre deux antennes de réception qui est utilisée pour estimer la différence de marche, elle-même 15 utilisée pour estimer l'attitude de l'engin 10 spatial. Ainsi, les erreurs spatiales induites par lesdits biais de phase peuvent tout à fait être estimées, pour chaque base de réception B1, B2, sous la forme d'une erreur sur une différence de marche ou sur une différence de phase entre les signaux radioélectriques reçus par les deux antennes de réception de la base de 20 réception considérée. Rien n'exclut cependant, suivant d'autres exemples, d'estimer les erreurs spatiales sous d'autres formes, y compris sous la forme de biais de phase pour chacune des antennes de réception A1, A2, A3. On établit ci-après un exemple de relation linéaire pouvant être utilisée 25 pour estimer les erreurs spatiales en fonction des écarts angulaires inter- satellites théoriques et des écarts angulaires inter-satellites estimés. Plus particulièrement, on considère dans cet exemple une paire de satellites 20 GPS de directions respectives um et un (directions d'incidence vm et vn dans le repère engin), et on établit une relation linéaire entre les 30 paramètres suivants : - une différence AYmn entre l'écart angulaire inter-satellites théorique et l'écart angulaire inter-satellites estimé déterminés pour ladite paire de satellites, en d'autres termes : AYmn = Vem'Ven - Um'Un = Vem'Ven - Vm.Vn un vecteur 5mn comportant les erreurs spatiales pour chaque base de réception B1, B2 et chaque direction d'incidence vm et vn. Plus particulièrement, on se place dans cet exemple dans le cas où le vecteur 6mn est exprimé sous la forme suivante : 5mn = (6rnE312 Eire ôrr132 5 expression dans laquelle : 6rB2m est l'erreur sur la différence de marche pour la base de réception B2 dans la direction d'incidence vm, ôrBlm est l'erreur sur la différence de marche pour la base de réception B1 dans la direction d'incidence vm, 6rB2n est l'erreur sur la différence de marche pour la base de réception B2 dans la direction d'incidence vn, ein est l'erreur sur la différence de marche pour la base de réception B1 dans la direction d'incidence vn. On note à ce stade que : une erreur uniquement sur la différence de marche de la base de réception B2 (c'est-à-dire 6rBlk = 0) correspond à une rotation autour de l'axe du repère engin de vecteur x, une erreur uniquement sur la différence de marche de la base de réception B1 (c'est-à-dire 6rB2k = 0) correspond à une rotation autour de l'axe du repère engin de vecteur y. Par conséquent, il est possible d'exprimer les directions d'incidences vm et vn comme étant obtenues, à partir de leurs estimations respectivement 25 vem et ven, par rotations successives autour des axes du repère engin de vecteurs x et y. Pour la direction d'incidence vk et son estimée vek, k égal à m ou n, cela revient à avoir vk = Mk.vek, avec : (1 0 0 r 1 0 Eyk Mk = 0 1 - Exk 0 1 0 Exk 1 Eyk 0 1 30 expression dans laquelle Exk et Eyk (k égal à m ou n) correspondent aux angles de rotation autour des axes du repère engin de vecteurs respectivement x et y, angles de rotation qui sont considérés comme faibles de sorte que cos(Exk) et cos(Eyk) sont sensiblement égaux à respectivement Exk et Eyk, et de sorte que sin(Exk) et sin(Eyk) sont sensiblement égaux à 1. On peut également négliger les termes de second ordre de sorte à obtenir la relation suivante : Mk ( 1 0 1 Eyk EXk En remplaçant, dans l'expression de la différence AYmn, les directions d'incidence vm et vn par respectivement Mm.vem et Mn.ven, on obtient une relation du type : ( vemx vent' - vemy - venz \T ( EXm\ tEXm\ - vemz - Venx + Vemx - VenzEYm EYn1 AYmn + o(E) = a(vem, ven ) - + o(E) - vemz - veny + vemy - venz EXn EXn vemz ve -nxve-m - x venz cYn j \eYni expression dans laquelle o(E) correspond aux termes de second ordre. En négligeant ces termes de second ordre, on obtient : (EXm\ AYmn a(vern, yen) - CYm EXn EYni Pour la base de réception B2, si on désigne par Are2k l'estimée de la différence de marche Ar2k, on a : ôrB2k = Are2k - Ar2k = b.vek.y - b.vk.y Du fait que l'erreur 5rB2k ne dépend pas de l'angle de rotation Eyk, on peut écrire : ( Par conséquent, on obtient : _ 6r-2k B = ID'exk.vekz. 5e2 = (1 0 vek - 0 1 Exk 0 (0 0 -Exk -vek -b-y= 0 0 1 \O - b - Exk 0 \ b - Exk - vek - y 0 17 De manière analogue, on obtient pour 5rBik l'expression suivante : k = _b.-- Ey Vekz- D'après l'expression ci-dessus de la différence AYmn, on a la relation suivante : ( 6re \ vemz \ " Vemx - Veny - Vemy - Venz Eire 1 - Vemz - Venx + Vemx - Venz Venu AYmn =7' D - Venu - Veny + Vemy - venz 6rnB2 vemz Venx - vemx - venz ) venz 6rnB1 venz c'est-à-dire : 1 AYmn = - ( ( B2\ Eirm 6rmB1 6r1E132 6rE31 n \T vemx - veny - vemy - venz vemz vemz - venx vemx - venz vemz - vemz - veny + vemy - venz venz -Vemz-Venx+Vemx-Venz Venz Sur la base de la relation précédente, donnée à titre d'exemple non limitatif, on comprend que tout algorithme connu de l'homme de l'art peut être mis en oeuvre pour estimer les erreurs spatiales à partir de la différence entre d'une part les écarts angulaires inter-satellites théoriques et, d'autre part, les écarts angulaires inter-satellites estimés. La relation précédente relie une observation scalaire (différence AYmn) à quatre paramètres à estimer (erreurs spatiales 6rB2m, Eer.B2n, &Bi n). En considérant par exemple un algorithme d'estimation du type moindres carrés ordinaires, on a une solution unique dès lors que l'on considère cinq satellites 20 GPS (Ns = 5) dans le champ de vue radioélectrique de l'engin 10 spatial. En effet, on a alors au plus dix paires de satellites GPS (au maximum Np = 10) pour lesquelles on peut avoir dix observations. On a en outre cinq directions d'incidence différentes et donc cinq erreurs spatiales à estimer pour chaque base de réception B1, B2, soit dix paramètres à estimer. D'autres algorithmes d'estimation peuvent être mis en oeuvre. En particulier, si on considère un algorithme d'estimation récursif (moindres carrés récursifs, filtre de Kalman, etc.), il n'est pas nécessaire d'avoir une solution unique à chaque itération car l'effet de mémoire inhérent à de tels algorithmes d'estimation permet d'affiner progressivement l'estimation des erreurs spatiales. Par exemple, dans le cas d'un algorithme d'estimation récursif, il est possible d'exécuter l'étape 55 d'estimation des erreurs spatiales paire de satellites 20 GPS par paire de satellites GPS (ou de considérer à chaque itération des groupes différents de paires de satellites). Il est à noter que, au cours du temps, des satellites 20 GPS peuvent entrer ou sortir du champ de vue radioélectrique de l'engin 10. Ainsi, à mesure que des satellites 20 GPS entrent dans le champ de vue radioélectrique de l'engin 20, ceux-ci peuvent être considérés au sein de paires pour la calibration des erreurs spatiales pour les directions d'incidence de ces satellites 20 GPS. En outre, les directions d'incidence d'une même paire de satellites 20 GPS se trouvant dans le champ de vue radioélectrique de l'engin 10 vont varier au cours du temps, du fait du défilement de l'engin 10 en orbite LEO. Une même paire de satellites 20 GPS peut donc être considérée plusieurs fois au cours du temps pour estimer les erreurs spatiales pour des directions d'incidence différentes. De manière plus générale, on comprend donc que le procédé 50 de calibration réalise avantageusement, dans des modes particuliers de mise en oeuvre, un filtrage temporel pour estimer les erreurs spatiales en fonction non seulement des observations à un instant donné, mais également en fonction d'observations antérieures. Un tel filtrage temporel est notamment inhérent aux algorithmes d'estimation récursifs. Plusieurs implémentations sont possibles pour la prise en compte de telles observations réalisées à des instants différents, et on comprend que le choix d'une implémentation particulière ne constitue qu'une variante de l'invention parmi d'autres. Suivant un exemple non limitatif, les étapes du procédé 50 de calibration d'erreurs spatiales peuvent être itérées à chaque fois qu'une nouvelle paire de satellites 20 GPS a été détectée. De préférence, l'étape 55 d'estimation des erreurs spatiales met dans ce cas en oeuvre un filtre de Kalman. Il est à noter que rien n'exclut d'estimer des erreurs spatiales pour chaque direction d'incidence des satellites 20 GPS se trouvant, au fil du temps et des défilements respectifs desdits satellites GPS et de l'engin spatial au-5 dessus de la surface de la Terre, dans le champ de vue radioélectrique de l'engin 10 spatial. Toutefois, dans un mode particulier de mise en oeuvre de procédé de calibration d'erreurs spatiales, le champ de vue radioélectrique de l'engin 10 spatial est préalablement partitionné en une pluralité de cônes distincts 10 d'incidence de signaux radioélectriques par rapport à l'engin 10. En outre, les erreurs spatiales sont considérées comme étant les mêmes pour toute direction d'incidence d'un même cône d'incidence. De la sorte, on peut n'estimer qu'une erreur spatiale par base de réception B1, B2 pour chaque cône d'incidence, ce qui permet de réduire le 15 nombre de calculs à effectuer et/ou la quantité d'erreurs spatiales à mémoriser. La figure 5 représente schématiquement le champ de vue radioélectrique (désigné par « CV ») de l'engin 10 spatial comme étant sensiblement semi-hémisphérique au-dessus d'un plan déterminé par les trois antennes de réception A1, A2, A3. La figure 5 représente également la 20 convention adoptée de manière non limitative pour exprimer une direction d'incidence vm sous la forme de deux angles d'incidence : un angle cp d'azimut, choisi comme étant l'angle entre le vecteur x et la projection de la direction d'incidence vm sur la plan déterminé par les trois antennes de réception A1, A2, A3, 25 un angle 8 d'élévation, choisi comme étant l'angle entre la direction d'incidence vm et le plan déterminé par lesdites trois antennes de réception A1, A2, A3. On peut par exemple définir un ensemble de cônes d'incidence adjacents, chacun desdits cônes d'incidence correspondant à une largeur de 5 30 degrés en azimut et 5 degrés en élévation. Dans une partition particulièrement avantageuse du champ de vue radioélectrique CV de l'engin 10 spatial, les cônes d'incidence d'élévation faible par rapport à un plan déterminé par les trois antennes de réception A1, A2, A3 de l'engin 10 spatial sont d'angles solides respectifs inférieurs à ceux de cônes d'incidence d'élévation élevée par rapport audit plan. Il a en effet été constaté que la composante indirecte est plus faible pour des directions d'incidence sensiblement normales (élévation élevée), par 5 rapport au plan déterminé par les trois antennes de réception, que pour les directions d'incidence rasantes (élévation faible) par rapport audit plan. De telles dispositions permettent d'avoir un pas d'échantillonnage spatial plus faible dans les zones où les erreurs spatiales varient beaucoup que dans les zones où lesdites erreurs spatiales varient peu. 10 La figure 6 représente schématiquement un exemple d'un tel champ de vue radioélectrique CV, dans lequel plus l'élévation 0 est faible et plus la largeur des cônes d'incidence, à la fois en azimut et en élévation, est faible. Le procédé 50 de calibration d'erreurs spatiales est avantageusement utilisé dans un procédé d'estimation d'attitude de l'engin 10 spatial. Par 15 exemple, les erreurs spatiales sont calibrées sous la forme d'erreurs sur une différence de marche des signaux radioélectriques pour chaque base de réception B1, B2. L'estimation d'attitude s'effectue par exemple en estimant des différences de marche en fonction de mesures de phase sur chacune des antennes de réception A1, A2, A3. Ensuite, les erreurs spatiales sont 20 compensées sur lesdites différences de marche estimées, et l'attitude de l'engin 10 spatial est estimée en fonction des différences de marche obtenues après compensation, par exemple au moyen d'un filtre de Kalman. La présente invention concerne également un système d'estimation d'attitude de l'engin 10 spatial en fonction de mesures de phase de signaux 25 radioélectriques reçus sur les trois antennes de réception A1, A2, A3. Le système d'estimation d'attitude comporte notamment, outre des moyens conventionnels, des moyens configurés pour calibrer les erreurs spatiales conformément à l'invention. Ces moyens se présentent par exemple sous la forme d'une unité de traitement comportant au moins un processeur et 30 une mémoire électronique dans laquelle est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter par le processeur pour calibrer les erreurs spatiales et estimer l'attitude de l'engin 10. Dans une variante, l'unité de traitement du système d'estimation d'attitude comporte des circuits logiques programmables, de type FPGA, PLD, etc., et/ou circuits intégrés spécialisés (ASIC), configurés pour effectuer tout ou partie de la calibration des erreurs spatiales. Il est à noter que le système d'estimation d'attitude est soit distribué entre l'engin 10 spatial et une station sol, soit entièrement embarqué dans ledit engin 10 spatial. De manière plus générale, la portée de la présente invention ne se limite pas aux modes de mise en oeuvre et de réalisation décrits ci-dessus à titre d'exemples non limitatifs. Par exemple, il a été considéré ci-avant que les Ns satellites étaient des satellites d'un système global de navigation par satellites. Rien n'exclut, suivant d'autres exemples, de considérer d'autres types de satellites. Dans un tel cas, les positions respectives des Ns satellites et de l'engin peuvent être obtenues par tout moyen connu de l'homme de l'art. Par exemple, lesdites positions peuvent être estimées par une station sol et transmises à l'engin pour estimation des écarts angulaires inter-satellites théoriques. En outre, il a été considéré ci-avant un engin 10 spatial. Rien n'exclut, suivant d'autres exemples de considérer un engin aérien. De plus, l'invention a été décrite en considérant des signaux 20 radioélectriques émis par des satellites. Rien n'exclut suivant d'autres exemples de considérer, en plus desdits satellites ou à la place desdits satellites, des émetteurs aériens ou terrestres. Par exemple, il est possible de considérer des signaux radioélectriques émis par des pseudolites terrestres. La description ci-avant illustre clairement que par ses différentes 25 caractéristiques et leurs avantages, la présente invention atteint les objectifs qu'elle s'était fixés. En particulier, la présente invention permet de calibrer, sans nécessiter de senseur d'attitude autre que les antennes de réception, les erreurs spatiales induites par des biais de phase affectant des mesures de phase de signaux radioélectriques reçus par lesdites antennes de réception. REVENDICATIONS1 - Procédé (50) de calibration d'erreurs, dites « erreurs spatiales », induites par des biais de phase affectant des mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception (A1, A2, A3) non alignées d'un engin (10), lesdits biais de phase dépendant d'une direction d'incidence des signaux radioélectriques par rapport à l'engin (10), caractérisé en ce qu'il comporte des étapes de : (51) obtention de positions respectives de l'engin (10) et d'au moins une paire de satellites (20) se trouvant dans un champ de vue radioélectrique de l'engin (10), (52) estimation de l'écart angulaire inter-satellites de la paire de satellites (20) en fonction des positions respectives de l'engin (10) et desdits satellites de la paire, de sorte à obtenir un écart angulaire inter-satellites dit « théorique », (53) estimation de directions d'incidence respectives des satellites (20) de la paire par rapport à l'engin (10) en fonction de mesures de phase de signaux radioélectriques reçus desdits satellites de la paire par les trois antennes de réception de l'engin (10), (54) estimation de l'écart angulaire inter-satellite de la paire de satellites (20) en fonction des directions d'incidence respectives des satellites de la paire par rapport à l'engin (10), de sorte à obtenir un écart angulaire inter-satellites dit « estimé (55) estimation des erreurs spatiales pour chacune des directions d'incidence respectives des satellites (20) de la paire par rapport à l'engin (10) en fonction de l'écart angulaire inter-satellites théorique et de l'écart angulaire inter-satellites estimé déterminé pour ladite paire de satellites (20). 2 - Procédé (50) selon la revendication 1, caractérisé en ce que le champ de vue radioélectrique de l'engin (10) est préalablement partitionné en une pluralité de cônes distincts d'incidence de signaux radioélectriques par rapport à l'engin (10), et en ce que les erreurs spatiales sont considérées comme étant les mêmes pour toute direction d'incidence à l'intérieur d'un même cône d'incidence.3 Procédé (50) selon la revendication 2, caractérisé en ce que les cônes d'incidence d'élévation faible par rapport à un plan déterminé par les centres de phase des trois antennes de réception (A1, A2, A3) de l'engin (10) sont d'angles solides respectifs inférieurs à ceux de cônes d'incidence d'élévation élevée par rapport audit plan. 4 Procédé (50) selon l'une des revendications précédentes, caractérisé en ce que les satellites de la paire appartiennent à un système global de navigation par satellites et les trois antennes de réception (A1, A2, A3) sont des antennes de récepteurs dudit système global de navigation par satellites, et en ce que les positions respectives des satellites (20) de la paire et de l'engin (10) sont déterminées à partir d'informations de navigation incluses dans les signaux radioélectriques reçus desdits satellites (20). 5 Procédé (50) selon l'une des revendications précédentes, caractérisé en ce que plusieurs paires différentes de satellites sont considérées au cours du temps, et en ce que les étapes du procédé sont itérées pour chaque paire de satellites considérée. 6 - Procédé (50) selon l'une des revendications 1 à 4, caractérisé en ce que, Ns satellites organisé en Np paires se trouvant dans un champ de vue radioélectrique de l'engin (10), les erreurs spatiales pour chacune des directions d'incidence respectives des Ns satellites (20) par rapport à l'engin (10) sont estimées simultanément en fonction des écarts angulaires inter-satellites théoriques et des écarts angulaires inter-satellites estimés déterminés pour lesdites Ns paires de satellites (20). 7 Procédé d'estimation d'attitude d'un engin (10) en fonction de mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception (A1, A2, A3) non alignées de l'engin (10), caractérisé en ce que ledit procédé comporte la calibration d'erreurs spatiales conformément à l'une des revendications précédentes, et la compensation desdites erreurs spatiales estimées lors de l'estimation de l'attitude de l'engin (10). 8 Produit programme d'ordinateur, caractérisé en ce qu'il comporte un ensemble d'instructions de code de programme qui, lorsqu'elles sontexécutées par un processeur, mettent en oeuvre un procédé selon l'une des revendications précédentes. 9 - Système d'estimation d'attitude d'un engin (10) en fonction de mesures de différences de phases de signaux radioélectriques reçus par au moins trois antennes de réception (A1, A2, A3) non alignées de l'engin (10), caractérisé en ce qu'il comporte des moyens configurés pour calibrer des erreurs spatiales conformément à l'une des revendications 1 à 6. - Système d'estimation d'attitude selon la revendication 9, caractérisé en ce que les trois antennes de réception (A1, A2, A3) sont des antennes de 10 récepteurs d'un système global de navigation par satellites. 11 - Engin (10) spatial, caractérisé en ce qu'il comporte un système d'estimation d'attitude selon l'une des revendications 9 à 10. 12 - Engin (10) spatial selon la revendication 11, caractérisé en ce que, ledit engin spatial comportant un corps (12) comportant plusieurs faces externes sensiblement planes, les trois antennes de réception (A1, A2, A3) sont agencées sur une même face externe dudit corps.