FR 3088784 A1 20200522 FR 1871833 A 20181115 PROCEDE ET DISPOSITIF POUR CONSERVER LE POINTAGE D'UNE ANTENNE AVEC UN SATELLITE L'invention concerne un procédé et un dispositif permettant de conserver la direction de pointage d'une antenne avec un satellite. Elle est utilisée pour des communications satellite (Satcom) où les antennes peuvent être montées sur des véhicules terrestres, des aéronefs, des navires. Pour assurer des communications par satellite sur des véhicules en mouvement, les antennes équipant ces véhicules doivent être dirigées to vers un satellite. Pour cela, les antennes sont montées sur des plateformes stabilisées désignées sous l'acronyme APS (Antenne Positioning Systems). Les exigences en termes de précision de pointage dépendent de la bande de fréquence de la télécommunication, de la taille de l'antenne ainsi que de sa technologie. 15 Plusieurs techniques de positionnement d'antenne sont connues de l'art antérieur. Les problèmes techniques à résoudre sont d'assurer un débit de liaison satellite maximal et de ne pas émettre en direction d'un satellite voisin. Une première technique ou « boucle ouverte » demande 20 l'utilisation d'une centrale de navigation de précision, appelée INS (Inertial Navigation System) ou AHRS (Attitude and Heading Reference System), afin de mesurer une attitude, un cap et la position GPS. Ce procédé tient compte des positions absolues. Une deuxième technique ou pointage en boucle fermée utilise le 25 principe d'écartométrie, c'est-à-dire le calcul des écarts angulaires de l'antenne par rapport au satellite. Ce calcul utilise un signal qui est l'image de la qualité de la liaison, appelé signal de tracking, issu par exemple d'un récepteur de balise satellite ou d'un modem. Cette technique utilise une figure d'écadométrie, c'est-à-dire un dépointage volontairement induit à la 30 ligne de visée de l'antenne pour faire varier le signal de tracking et ainsi estimer l'écart au satellite. La méthode en boucle fermée utilisée prend le 2 nom de sa figure d'écartométrie. Il existe deux familles de figures : le « Step Tracking », qui injecte des mouvements intermittents et le « Conical Scan », qui injecte des mouvements continus. Une troisième méthode est la stratégie hybride, qui juxtapose une 5 stratégie en boucle ouverte et en boucle fermée. Une quatrième technique décrite dans la publication de E. Broussard et al, intitulée « Hybrid satellite position estimator using self-induced dither for Satcom on the Move » présente une technique qui utilise des mesures fournies par au moins deux gyromètres positionnés sur l'embase du positionneur. Le balayage ou « scan » nécessaire pour l'estimation de la direction de pointage se fait sans mouvements volontaires de l'antenne. Les mouvements de scan sont obtenus notamment avec le bruit de mesure. Les techniques mises en oeuvre dans l'art antérieur présentent différents avantages, il reste toutefois quelques axes à améliorer pour optimiser la conservation de direction de pointage d'une antenne. L'idée de la présente invention repose sur une nouvelle approche positionnant un dispositif configuré pour mesurer des vitesses angulaires situé au niveau de l'antenne même et d'un traitement de données approprié. L'invention concerne un procédé pour conserver une direction de pointage d'une antenne dans l'axe d'un satellite caractérisé en ce qu'il exécute au moins les étapes suivantes : Pour i variant de 1 à n, n étant un nombre de mesures données, Mesurer une valeur du signal de tracking SNR, reçu depuis le satellite pour un instant donné, Exprimer la valeur du signal de tracking SNR, en fonction des paramètres KI et K2 liés au diagramme de l'antenne, d'une période de temps j, et d'une valeur de rapport signal à bruit SNR0 obtenue sans dépointage volontaire de l'antenne, ( SNR, = 1(1. ez - drIszi + K2. Ey - drI E 3; 1 + SNR0 ri ri 2 2 i-1 3 avec : ez et ey des valeurs d'écarts au satellite, i variant de 1 à n, indice d'une mesure de SNR dans une fenêtre de temps donnée, 5 j indice de la période de temps précédent d'un indice i, on somme sur toutes les mesures précédentes, ez* et e); j les vitesses angulaires pour au moins deux axes y et z, différents de l'axe de pointage de l'antenne vers le satellite, dr une période d'échantillonnage, 10 Définir une fonction de coût à partir de la valeur de SNR, et de sa valeur estimée SNRi Zi-i )2 SNRi = 1(1.( ez - dr1 21) + K 2.(ry - d r lEy. j ± SNRO Er ri = SNRi - SNR j = I(p,.Erri)2 J =1 pi est un facteur de pondération choisi par l'opérateur, Trouver les paramètres optimaux ry et ez au sens des moindres carrés, c'est-à-dire minimisant la fonction J à chaque pas de temps, 15 Utiliser le résultat de la comparaison des valeurs estimées à la valeur de référence pour conserver la direction de pointage. Selon une variante de réalisation, on exprime la valeur SNR sous la forme d'un système linéaire en regroupant les résultats des n mesures dans une première matrice A exprimant les Eij11E;, et les E 4, et une 20 deuxième matrice X regroupant les valeurs d'écart au satellite et la valeur de rapport signal à bruit SNR0, on estime un vecteur X à partir du système linéaire YrAX en utilisant une méthode des moindres carrés ordinaire, pour chaque nouvelle valeur de signal de balise SNR, en déduit une valeur d'écart 4 angulaire ry selon l'axe y et une deuxième valeur angulaire estimée ez selon l'axe z: = -X i-1 e ; il J=1 -IE. J=1 X = [2. Ki d - Ez 2 . K2 . Ey (SN RO + K2 . + Ez2)1 Yi 2 2 -= SNRi- Ki.( dl-lez., - K2-(cl-cle;,j ) avec - Évé3, : les vitesses angulaires mesurées par des gyromètres situés au 5 niveau de l'antenne, - t: la période d'échantillonnage, : l'indice de la mesure dans la fenêtre de temps, i e - j: l'indice de la mesure précédant la mesure courante i, n: le nombre de mesures utilisées dans la fenêtre de temps, 10 les n = -tir mesures de IV, et Y, sont regroupées dans deux matrices A et Y: - 0 0 1 2 2 1 Yj j=1 j=1 -XE" -X zi Yi j=1 j=1 n-1 - re-1 . , 1 j=1 j=1 { - y On estime le vecteur X en utilisant une méthode des moindres carrés ordinaire : Î = (AT . A)- 1. AT. A Y, AT étant la transposée de la matrice A, 5 On obtient une estimation ey, É; pour chacune des vitesses angulaires Éz, à partir de l'estimation de X: 2. Kr dr. ez -= k21 = 2. K2. CIT.ty ri; KiZ2 + K2.-42 ± SNR0 = et r = 2.K . 2.K2.dr 1.dr 5 On prend en compte les biais de mesure des vitesses angulaires et on exécute les étapes suivantes : 2 i-1 SN Ri =- dr1(éz,,,.) + (i - 1). dr. bz + 1=1 j ( K2 Ey i drI(Ey .) i 1). dr. by + SN Ro m I avec : E : les vitesses angulaires mesurées par les gyromètres de 10 l'antenne, - bz, by les biais des gyromètres de l'antenne, - dr : la période d'échantillonnage, - I: l'indice de la mesure courante dans la fenêtre de temps, i e [1...n], - j: l'indice de la mesure précédant la mesure courante 15 - n: le nombre de mesures utilisées dans la fenêtre de temps, On détermine les estimations optimales ez , , h; et by et SNR0 des paramètres s, sy, bz, by et SNR0 au sens des moindres carrés. L'estimation des biais se présente sous la forme d'un système linéaire dans le cas où les estimés des écarts C;, ri, sont issues d'un module 20 de traitement de données et utilisés comme entrées de l'estimateur de biais : Ybi = Wbe Xb avec : i-1 6 2 2 i-1 i-1 Ybi = SA Ri - Ki - ch-1 Ézmi - K2 ey - dr L'y mi J=1 1=1 Wb; = F 2 1 ( G - drI 4, . (i - 1). ch 2. K2. Érz,- drlé), mi . . (i - 1).dr (i - 1)2. dr2 11 4 est la matrice regroupant n lignes de Wbi. Yb est la matrice regroupant n lignes de L'estimation g;', = (Ab T. Ab)-1 . Ab T. Yb permet d'estimer les biais 17, et by : = [6; h.-3, (Ki. bz2 + K2. by2) SNR01T On estime les biais g et ii;, puis on calcule les vitesses angulaires, ty, é; 5=E - b Yi Ymi Y = e; - bz. Les deux axes y et z pour déterminer des vitesses angulaires sont choisis orthogonaux. On peut utiliser la valeur zéro comme valeur de référence pour 10 l'asservissement. Selon une variante, on utilise pour l'estimation un mouvement de l'antenne en continu, généré par n'importe quel phénomène faisant bouger le faisceau de l'antenne par rapport au satellite, tel qu'un bruit de mesure, une flexibilité mécanique, une erreur d'asservissement, ou un mouvement 15 volontaire réalisé par un opérateur. Selon une autre variante, on estime en continu la valeur des écarts angulaires et on utilise une boucle d'asservissement pour maintenir en position l'antenne dans une direction de pointage optimale. L'invention concerne aussi un dispositif pour conserver une 20 direction de pointage d'une antenne dans l'axe d'un satellite caractérisé en ce qu'il comporte au moins les éléments suivants : Un module de mesure d'une valeur d'un signal SNRÉ émis par le satellite pour un instant i, 7 Un dispositif de détermination d'au moins une première vitesse angulaire et d'au moins une deuxième vitesse angulaire au niveau de l'antenne, pour au moins deux axes y et z, non orthogonaux et différents de l'axe de pointage de l'antenne vers le satellite, 5 Un module de traitement des données configuré pour exécuter les étapes du procédé selon l'invention, recevant les valeurs mesurées des vitesses angulaires, la valeur du rapport signal sur bruit mesurée SNR et générant un signal de commande à un dispositif de correction des écarts angulaires. Le dispositif utilisé pour déterminer les vitesses angulaires est un 10 gyromètre ou comprend un estimateur de biais des mesures. Le dispositif peut comporter une boucle d'asservissement configurée pour corriger en continu un écart angulaire. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description d'exemples de réalisation 15 donnés à titre illustratif et nullement limitatif, annexée des figures qui représentent : - Figure 1, un exemple de système comprenant un aéronef en communication avec un satellite, - Figure 2, un schéma représentant le traitement des données pour 20 conserver la direction de pointage, et - Figure 3, une variante de réalisation permettant de calculer la valeur du biais d'une mesure. Afin de bien faire comprendre le procédé et son dispositif associé, la description qui suit est donnée, à titre d'exemple, pour un avion équipé 25 d'une antenne de communication utilisée pour la liaison Satcom (figure 1). Le satellite 1 émet un signal. L'avion 2 est équipé d'une antenne 3 montée par exemple sur son fuselage 4 via un support d'antenne 5. Le montage de l'antenne et de son support est connu de l'art antérieur et ne sera pas détaillé. L'antenne est équipée d'un récepteur 6 de signal émis par le satellite. 8 Un dispositif 7 configuré pour mesurer des vitesses angulaires comprenant, par exemple, au moins deux gyromètres 7y, 7z, qui sont associés directement à l'antenne, et qui sont fixés, par exemple, sur l'antenne. Les deux gyromètres 7y, 7z, sont en liaison avec un processeur 8 5 pourvu d'un module de traitement des données 9 qui va recevoir les valeurs des vitesses angulaires mesurées par les gyromètres et une mesure du facteur de qualité de la liaison satellite, ou signal de tracking, à maximiser, pour les traiter afin de fournir un signal représentatif des écarts mesurés. Le facteur de qualité peut être n'importe quel signal représentant l'image du 10 débit de la liaison satellite. Pour une description plus concrète, l'exemple est donné pour un rapport signal à bruit SNR (Signal over Noise Ratio) issu d'un récepteur de balise satellite. Les mesures et les résultats du traitement des données seront par exemple stockés dans une mémoire 15 de l'avion. Les valeurs des écarts angulaires sont comparées à une valeur de référence, en 15 général la valeur zéro ou une valeur voisine de zéro, le résultat étant transmis à un dispositif de correction ou de compensation d'une erreur de position angulaire 10. Le dispositif de correction va agir sur le support de l'antenne afin de maintenir le pointage de l'antenne ou le ramener dans une position optimale afin de conserver une valeur de rapport signal sur bruit 20 SNR maximale. Le processeur 8 peut faire partie des équipements de bords de l'appareil. L'antenne est reliée à un dispositif 11 de mesure de signal balise. Les gyromètres permettant les mesures sont positionnés directement au niveau de l'antenne et non au niveau du support d'antenne 25 comme il est décrit dans l'art antérieur. La figure 2 schématise un exemple d'enchaînement des étapes du procédé selon l'invention. La valeur du SNR est mesurée par le récepteur 6 de signaux et transmise au processeur 8 et au module de traitement des données 9 30 configuré pour exécuter les étapes du procédé selon l'invention. La valeur du 9 SNR dépend notamment d'une valeur SNRo (valeur maximale du rapport signal à bruit disponible pour le système lorsqu'il n'y a pas de dépointage). Le module de traitement des données 9 reçoit aussi des valeurs de vitesses angulaires estimées par les deux gyromètres, au moins pour la 5 vitesse angulaire selon l'axe y et pour la vitesse angulaire selon l'axe z. Les axes y et z ne sont pas nécessairement orthogonaux. Les axes choisis pour mesurer les écarts au satellite (angles exprimés en degrés) sont différents de l'axe de pointage de l'antenne aligné dans la direction du satellite. Le diagramme 12 de la figure 2 illustre la perte de SNR résultant 113 d'un dépointage de l'antenne. Il dépend de la taille, la forme et la technologie de l'antenne ainsi que de la fréquence de fonctionnement. Le procédé va rechercher au moins deux valeurs d'écarts au satellite, es,, ez, respectivement selon l'axe y de l'antenne et selon l'axe z de l'antenne. 15 Le SNR s'exprime en fonction des écarts, £3, et sz, des coefficients négatifs Ki et K2 liés à la réalisation l'antenne et d'une valeur de SNRo qui est inconnue : SNR = Ki.a+ K2.4 + sNR0 (1). La vitesse angulaire de l'antenne fournit une mesure de la dérivée 20 temporelle de ces écarts éz et 4. Ki et K2 sont des coefficients ou des constantes liées au diagramme de l'antenne qui sont connues ou estimées, SNR0 correspond à la valeur maximale du rapport signal à bruit disponible pour le système lorsqu'il n'y a pas de dépointage. La mise en oeuvre de ces calculs (équ.3) 25 requiert des mouvements de l'antenne et donc un dépointage, pour avoir une estimation correcte. Le dépointage peut être naturel, lié par exemple au bruit de mesure dans l'estimation naturellement présent dans le système, ou volontaire réalisé par un opérateur. Selon une première variante de réalisation qui ne tient pas compte 30 de la valeur du biais des mesures, la première étape de traitement exécutée 10 au niveau du module de traitement consiste à exprimer le diagramme d'antenne en fonction du dépointage de l'antenne. Pour un indice i d'une fenêtre de temps, la valeur du SNR, mesurée s'exprime de la manière suivante : SNRi = Ki.(ez - drIE; I) + K2 (Ey - 5 (3) avec : e" s" correspond à la vitesse angulaire, mesurée par les gyromètres et fait référence à la dérivée temporelle des valeurs d'écarts au satellite recherchés E E 10 i, variant de 1 à n, est l'indice d'une mesure dans la fenêtre de temps, composée de n mesures. n est un paramètre de réglage lié à la fenêtre de temps T = n. dr. Les grandeurs indicées 1 correspondent au pas de temps le plus récent, celles indicées n sont les plus anciennes. j, variant de 1 à est l'indice d'une mesure précédant la mesure i, 15 ez'j et Ej,j sont les vitesses angulaires mesurées par les gyromètres de l'antenne pour l'indice j, dr est la période d'échantillonnage, par exemple choisie en fonction de l'architecture du système, du capteur de réception des signaux, de manière à obtenir suffisamment d'informations pour obtenir une valeur représentative 20 sans saturer le module de traitement, elle sera notamment fonction de la fréquence d'échantillonnage du signal balise. La taille de la fenêtre temporelle est choisie de façon à avoir suffisament d'échantillons pour dissocier le signal utile du bruit, ceci selon des principes connus de l'homme du métier dans les traitements statistiques. 25 On cherche à estimer les paramètres ez, ey et SNR0 optimaux au sens des moindres carrés. C'est-à-dire les valeurs ez, ry et SNR0 qui minimisent la fonction de coût J. La fonction de coût représente la somme sur 2 i-1 i-1 2 j=1 T E) y] J=1 SNR0 11 toute la fenêtre d'observation des erreurs au carré entre la mesure SN RE et la valeur estimée SNI?, : 2 2 1-1 t-1 SN R, = - dr £.2j) + K2. (ry d-rIE; I) + SNR0 Erri = SN Ri - SNRi j = (p,.Erri)2 J=1 pi est un facteur de pondération choisi par l'opérateur. Il est possible d'utiliser toute méthode de résolution des moindres carrés 5 connue de l'homme du métier. Pour l'illustration de cet exemple, une approche linéaire est suivie. L'équation (3) peut se mettre sous la forme : SN Ri - Ki.( el% - K2. cirIe;,j = -2. Ez. Ki. cillez] - 2. Ey. K2. cit-IE; i + SN RO + K2. e, +K1. (4) j=1 j=1 L'étape suivante consiste à exprimer l'équation précédente (4) 10 sous la forme d'un système linéaire en ses paramètres : 2 1-1 i-1 2 Yi = SN Ri - Ki.(drIezi) - K2 (dit E) Yi = 1-1 Wi 4 F-1 1=1 1-1 -1 6 1 .1=1 YI X r= [2. Ki. citez .2 K2. dr. Ey (SNRO + K2- Kl- Ez2 )](5) avec : les vitesses angulaires mesurées par les gyromètres situés au niveau de l'antenne, 2 i-1 2 12 - dr : la période d'échantillonnage, - : l'indice de la mesure dans la fenêtre de temps, i E n], - j: l'indice de la mesure précédant la mesure courante i, j e [1 (i - 1)1, n : le nombre de mesures de vitesses angulaires utilisées dans la 5 fenêtre de temps. Avec une fenêtre de temps de T secondes, et en exécutant les étapes précédemment décrites sur les n mesures, on obtient n = Lar mesures de Wi et Yi, qui sont regroupées dans deux matrices A et Y s'exprimant de la manière suivante : 1 2 2 1 j=1 j=1 1 - s;.1 j=1 j=1 n-1 - ezl _IEy" 10 1 lu Y = };71 On obtient le système d'équation surdéterminé suivant : Y = A. X. On estime le vecteur X en utilisant la méthode des moindres carrés ordinaires : X = (AT. A)-1. AT. Y. 15 AT est la transposée de la matrice A. On obtient une estimation des écarts angulaires ry, ez à chaque nouvelle valeur de signal balise SNR à partir de l'estimation de X: 13 2. K1. dr. = = 2. K2. dr. Éj, I (7) 11-3- Kirez2 K2.-E-y2 Arirto ey 2. K2. dr e = 2 K dr A la fin de l'exécution des étapes décrites précédemment, le processeur dispose d'une valeur estimée pour l'écart au satellite suivant l'axe 5 y et une valeur estimée pour l'écart au satellite selon l'axe z. Les valeurs estimées ty, rz, pour les écarts au satellite sont ensuite comparées au niveau du module de traitement des données à une valeur de référence, la valeur « 0 » qui est la consigne d'écart nul ou une valeur proche de zéro. Si les valeurs des écarts au satellite estimées sont 10 trop éloignées de cette valeur de référence ou prouvent l'existence d'un dépointage trop important, alors le module de traitement 9 va émettre un signal de commande Sr vers le contrôleur de position 10 de l'antenne afin de la repositionner (de la repointer) pour obtenir une valeur de SNR optimal. La correction à apporter est réalisée en continu en utilisant une 15 boucle d'asservissement qui permet à tout moment de faire tendre cette valeur d'écart vers zéro. Lorsque le véhicule sur lequel est montée la station Satcom n'a pas un mouvement trop rapide, il est possible de déterminer l'écart non plus en continu, mais à des instants donnés. 20 La figure 3 illustre une variante de réalisation qui tient compte du biais présent dans les vitesses angulaires mesurées. La vitesse angulaire mesurée ém est égale â: Èrn = É + b avec - É : la vitesse angulaire vraie, non biaisée, 25 - b: le biais de mesure à identifier. 14 L'objectif est de prendre en compte le biais de mesure afin de le supprimer et d'améliorer le calcul. L'hypothèse suivante est posée : - Le biais est constant (varie lentement dans le temps par rapport aux 5 autres signaux en jeux). Pour un indice i, la valeur du SNRi mesurée s'exprime de la manière suivante : SNRi = 1(1. i-1 2 i -1 2 ( Ez - dr I E.2 + K2 + SNR0 (8) ) ) E y - drIs; j ) j=1 j=1 i-1 2 SNRi = Ki.( Ez - dr»zryt ) + (1 - 1). dr. bz + ' é-1 2 K2 (Ey - dr I (Éynt,) + (i - 1). dr. by + SNR0 f=i avec : 10 - E Ymi : les vitesses angulaires mesurées par les gyromètres de l'antenne, - bz, by : biais des gyromètres de l'antenne, - tir: période d'échantillonnage, - : indice de la mesure courante dans la fenêtre de temps, i e [1 n], 15 - j: indice de la mesure précédant la mesure courante, j c [1 (i - 1)], - n: nombre de mesure utilisée dans la fenêtre de temps. On détermine les estimés optimales ey, ez, 6, by et SNR0 au sens des moindres carrés des paramètres 22, ey, bz, by et SNR0. Il est possible d'utiliser toute méthode de résolution des moindres carrés connue de 20 l'homme du métier. 15 Dans le cas où les paramètres Ey, tz sont préalablement connus, un estimateur de biais 20 est disposé entre le récepteur de signaux 6 et le module de traitement des données 9 (figure 1). L'estimateur de biais 20 reçoit en entrée : 5 - La mesure de vitesse angulaire des gyromètres : - La mesure de signal balise : SNRi, - L'estimation des écarts au satellite, c'est-à-dire la sortie du module de traitement des données 9 : ey, ez. Dans ce cas où les valeurs estimées des écarts É,Ç sont issues 10 d'un module de traitement de données et utilisées comme entrées de l'estimateur de biais, une résolution par les moindres carrés ordinaires (linéaire) est utilisée pour calculer /5, by et SNR0, Ybi= WIDE- Xb (9) 15 avec Y 2 é-1 yrri2 b1 =SNR -Ki(-drI4'I-K2 i-drIé j Wb i = F 2. K1 ij - dr 1 Ê . (i - 1). dr 2. K2 . -5, - dr Iéy . .(i -1).dr (i - 1)2 . dr2 Ah est la matrice regroupant n lignes de Wbi. Yb est la matrice regroupant n lignes de L'estimation Xb = (AbT.Ab)_l.Abr .Yb permet d'estimer g, et r) ; et .5W/1-0: = [ 13; 6-3, (K,. bz2 + K2. by2) SNR0 Ir A partir de l'estimation du biais, fi; et 1-5;' on calcule les vitesses angulaires, 20 =E - b Yi Ami Y iJ 16 E- = E - y, Yin, Y - bz. ..mj Les vitesses angulaires estimées ryi et s",1 sont alors utilisées comme entrées de Y1 et 'Nt; (4) pour l'estimation suivante des écarts 5 angulaires ey et É. Si les valeurs des écarts au satellite estimées sont trop éloignées de cette valeur de référence ou prouvent l'existence d'un dépointage trop important, alors le module de traitement 9 va émettre un signal de commande Sr vers le contrôleur de position 10 de l'antenne afin de la 10 repositionner (de la repointer) pour obtenir une valeur de SNR optimale. Sans sortir du cadre de l'invention, il est possible d'utiliser n'importe quelle méthode de résolution linéaire (méthodes récursives des moindres carrés, moindres carrés pondérés, moindres carrés totaux etc.) ou non-linéaire (optimisation numérique). 15 Le procédé et le dispositif selon l'invention présentent, en plus d'assurer un pointage optimal de l'antenne dans la direction du satellite afin de conserver une valeur maximale de signal de suivi, les avantages suivants : - Moins de capteurs nécessaires pour le suivi des satellites, 20 - Les valeurs des écarts estimées sont indépendantes de l'architecture du positionneur, - Les mesures et les calculs en temps réel permettent une estimation des biais par le pointage, REVENDICATIONS1 - Procédé pour conserver une direction de pointage d'une antenne dans l'axe d'un satellite caractérisé en ce qu'il exécute au moins les étapes suivantes : Pour i variant de 1 à n, n étant un nombre de mesures données, Mesurer une valeur d'un signal de tracking SNRi reçu depuis le satellite pour un instant donné, Exprimer la valeur du signal de tracking SNRi en fonction des paramètres K, 10 et K2 liés au diagramme de l'antenne, d'une période de temps j, et d'une valeur de rapport signal à bruit SNR0 obtenue sans dépointage volontaire de l'antenne, 2 2 t-1 t-1 SNR i Ki.(Ez - ciTIE; + K2. (Ey - dT EY)) SNR0 j=1 j=1 avec : Ez et Ey des valeurs d'écarts au satellite, 15 i variant de 1 à n, indice d'une mesure de SNR dans une fenêtre de temps donnée, j indice de la période de temps précédent d'un indice i, on somme sur toutes les mesures précédentes, Ez..) et syj les vitesses angulaires pour au moins deux axes y et z, différents 20 de l'axe de pointage de l'antenne vers le satellite, dr une période d'échantillonnage, Définir une fonction de coût à partir de la valeur de SNRi et de sa valeur estimée SNRi SNRi = Ki.(Ez - drIE; j) + K2. (ty - EYJ SNR0 j=1 j=1 Er n = SNRi - SNR 18 = I(pi. Erri)2 pi est un facteur de pondération choisi par l'opérateur, Trouver les paramètres optimaux ey et é; au sens des moindres carrés, c'est-à-dire minimisant la fonction J à chaque pas de temps, Utiliser le résultat de la comparaison des valeurs estimées à la valeur de 5 référence pour conserver la direction de pointage. 2 - Procédé selon la revendication 1 caractérisé en ce que : on exprime la valeur SNR sous la forme d'un système linéaire en regroupant les résultats des n mesures dans une première matrice A exprimant les 10 EfiClisz'j et les Eij4 Ey. j et une deuxième matrice X regroupant les valeurs d'écart au satellite et la valeur de rapport signal à bruit SNRD, on estime un vecteur X à partir du système linéaire Y=AX en utilisant une méthode des moindres carrés ordinaire, pour chaque nouvelle valeur de signal de balise SNR, on déduit une valeur d'écart angulaire ey, selon l'axe y 15 et une deuxième valeur angulaire estimée É; selon l'axe z: 2 i-1 i-1 )2 =-- SNRI - - K2.(cir E. -1E); j j=1 X = [2. Ki. dr. Ez 2 . K2. dr E y (SNR0 + K2. 4 + ez2 )] avec e" : les vitesses angulaires mesurées par des gyromètres situés au niveau de l'antenne, dr : la période d'échantillonnage, 20 j : l'indice de la mesure dans la fenêtre de temps, i E [1..7/], j: l'indice de la mesure précédant la mesure courante i, 19 - n: le nombre de mesures utilisées dans la fenêtre de temps, les n = - mesures de Wi et Yi, sont regroupées dans deux matrices A et Y: dr 0 0 1- 2 2 Y, 1 A -1Ey. , 1 wn- j=1 j=1 171. On estime le vecteur X en utilisant une méthode des moindres carrés ordinaire :1? = (Ar. A)-H1. -T Y , AT étant la transposée de la matrice A, 5 On obtient une estimation ry, ez, pour chacune des vitesses angulaires , E.z E. à partir de l'estimation de X: ) fri 2. Ki. dr. é; ii = Figil = 2. K2. dr. ry Ili K1.-2 + K2.-4-2 + SN Ro - et ez = ii . Y 2.K2.dr 2.K1.dt 3 - Procédé selon la revendication 1 caractérisé en ce que l'on prend en compte les biais de mesure des vitesses angulaires et on exécute les étapes suivantes : On exprime : ( i-i SN Ri = Ki. ez - drl(iz) ( )2 - 1). dr . bz + 1=1 20 avec : K2 (EY 2 Ch (É y ) (i 1) - ch. by + SNRo 771j j=1 les vitesses angulaires mesurées par les gyromètres de l'antenne, 5 - bz, by : les biais des gyromètres de l'antenne, dr : la période d'échantillonnage, : l'indice de la mesure courante dans la fenêtre de temps, i E [1... n], j: l'indice de la mesure précédant la mesure courante i, n: le nombre de mesures utilisées dans la fenêtre de temps, 10 On détermine les estimations optimales ez, g, g, et by et SNR0 des paramètres E, E, bz, by et SNR0 au sens des moindres carrés. 4 - Procédé selon la revendication 3 caractérisé en ce que l'on exprime le problème d'estimation des biais sous la forme d'un système linéaire dans le 15 cas où les estimées des écarts (,,r,y sont issues d'un module de traitement de données et utilisés comme entrées de l'estimateur de biais Yb1 = Wbi- XI) avec ;1 -= SNRi - Ki tz - drIèz',, - K2(e; - drI é ) ( y ',j i--1 t--1 [2.1(1 (fz - dr X Ézm ). (i - 1).dr 2. K2. e7,- drI É y m (i - 1). dr (i - 1)z di' 1=1 Ab est la matrice regroupant n lignes de Wb,. 4 est la matrice regroupant n 20 lignes de Yb,. -1 L'estimation firi, =- (AbT. Ab) .AbT.Yb permet d'estimer les biais g, et Ç: 2 2 j=1 J j=1 11 21 = [6; -ff; (K1.b,2+ K2. b3,2) SNRo]T On estime les biais 6 et by, puis on calcule les vitesses angulaires, Éy, E. --=-E - b Yrnj Y E.z) = E;111) S;- 5 5 - Procédé selon l'une des revendications 1 à 3 caractérisé en ce que les deux axes y et z pour déterminer des vitesses angulaires sont orthogonaux. 6 - Procédé selon l'une des revendications 1 à 4 caractérisé en ce que l'on utilise la valeur zéro comme valeur de référence pour l'asservissement. 10 7 - Procédé selon l'une des revendications 1 à 4 caractérisé en ce que l'on utilise pour l'estimation un mouvement de l'antenne en continu, généré par n'importe quel phénomène faisant bouger le faisceau de l'antenne par rapport au satellite (1), tel qu'un bruit de mesure, une flexibilité mécanique, une erreur d'asservissement, ou un mouvement volontaire réalisé par un opérateur. 8 - Procédé selon l'une des revendications précédentes caractérisé en ce que l'on estime en continu la valeur des écarts angulaires et on utilise une 20 boucle d'asservissement pour maintenir en position l'antenne dans une direction de pointage optimale. 9 - Dispositif pour conserver une direction de pointage d'une antenne dans l'axe d'un satellite caractérisé en ce qu'il comporte au moins les éléments 25 suivants : Un module de mesure d'une valeur d'un signal SNR, émis par le satellite pour un instant i, 22 Un dispositif de détermination d'au moins une première vitesse angulaire et d'au moins une deuxième vitesse angulaire au niveau de l'antenne, pour au moins deux axes y et z, non orthogonaux et différents de l'axe de pointage de l'antenne vers le satellite, 5 Un module de traitement des données (9) configuré pour exécuter les étapes du procédé selon l'une des revendications 1 à 8, recevant les valeurs mesurées des vitesses angulaires, la valeur du rapport signal sur bruit mesurée SNR et générant un signal de commande à un dispositif de correction des écarts angulaires (10). 10 10 - Dispositif selon la revendication 9 caractérisé en ce que le dispositif permettant de déterminer les vitesses angulaires est un gyromètre (7y, 7z). 11 - Dispositif selon la revendication 9 caractérisé en ce que le dispositif 15 permettant de déterminer les vitesses angulaires comprend un estimateur de biais des mesures (20). 12 - Dispositif selon l'une des revendications 9 ou 10 caractérisé en ce qu'il comporte une boucle d'asservissement configurée pour corriger en continu un écart angulaire.