FR 3047568 A1 20170811 FR 1600197 A 20160205 1 Méthode de calibrage d'un récepteur de radio-navigation par satellites L'invention concerne le domaine des récepteurs de radio-navigation par satellites, encore appelés récepteurs GNSS de l'acronyme anglo-saxon couramment utilisé « Global Navigation Satellite Systems », et plus précisément les récepteurs multi-antennes comprenant plusieurs voies de réception. L'invention concerne également le domaine du traitement des signaux à spectre étalé qui englobent les signaux de radio-navigation par satellite. L'objet de l'invention est une méthode de calibrage d'un récepteur GNSS multi-voies conduisant à la détermination d'un filtre d'égalisation large bande et/ou d'un filtre d'égalisation bande étroite, ainsi qu'un récepteur comprenant de tels filtres. Dans un récepteur multi-voies, un problème général concerne la calibrage des différentes voies de réception afin de compenser les différences en temps, phase et amplitude qui peuvent impacter le signal reçu sur les différentes voies. En effet, chaque voie de réception comprend une antenne ainsi qu'une chaine de réception radiofréquence qui comporte notamment un ou plusieurs filtre(s) de mise en forme et un ou plusieurs amplificateur(s). Ces différents composants peuvent ne pas présenter exactement les mêmes caractéristiques, en particulier, les temps de groupe, les gains en amplitude et en phase peuvent être sensiblement différents. Ces différences entraînent un désappairage des voies de réception qui doivent, en théorie, toutes avoir la même réponse en fréquence. Les signaux reçus sur les différentes voies de réception peuvent ainsi être impactés par ces différences et les traitements appliqués aux signaux en vue de déterminer des informations de positionnement du récepteur ou de direction d'arrivée des signaux peuvent être perturbés. Dans un contexte d'applications qui nécessitent une grande fiabilité de l'information de positionnement, il est donc nécessaire de calibrer très 3047568 2 précisément les récepteurs afin d'éviter au maximum les distorsions entre les différentes voies de réception. Cette nécessité est notamment importante pour des applications aéronautiques ou liées aux transports. 5 Les solutions connues permettant de résoudre le problème de calibrage d'un récepteur GNSS multi-voies peuvent être basées sur un pré-calibrage en usine effectué sur chaque équipement et utilisant des mesures en laboratoire ou en chambre anéchoïque. La réponse des voies de réception RF étant sensible à la température, il est nécessaire d'utiliser des 10 composants thermostatés ou des senseurs thermiques. D'autres solutions utilisent un signal radiofréquence dédié au calibrage qui est généré par un dispositif spécifique interne ou externe au récepteur. De manière générale, les solutions existantes présentent l'inconvénient de nécessiter des dispositifs ou équipements supplémentaires. 15 On connait par ailleurs d'autres solutions décrites dans les documents EP805510, US5525997 et US7304605. L'invention propose une méthode de calibrage d'un récepteur multivoies GNSS qui ne nécessite pas l'utilisation d'un générateur de signal 20 spécifique et qui peut être mise en oeuvre directement à partir de simples mesures réalisées sur un récepteur en fonctionnement. L'invention consiste notamment à déterminer un premier filtre d'égalisation large bande qui peut être positionné en sortie des voies de réception RF et en entrée des corrélateurs pour corriger le désappairage 25 entre les différentes voies de réception RF. L'invention consiste également à déterminer un second filtre d'égalisation bande étroite pour corriger les erreurs de gain et phase résiduels. L'invention s'applique à tout système de positionnement par satellite, notamment à des récepteurs compatibles des systèmes GNSS GPS, 30 GALILEO ou GLONASS. 3047568 3 L'invention a ainsi pour objet une méthode de calibrage d'un récepteur de signaux de navigation par satellites, ledit récepteur comprenant une pluralité de voies de réception comprenant chacune une antenne de réception et une pluralité de corrélateurs pour corréler le signal reçu par 5 l'antenne avec un code local retardé d'un retard paramétrable, ladite méthode de calibrage comprenant la détermination d'au moins un filtre égaliseur large bande destiné à être appliqué entre l'antenne de réception et les corrélateurs de chaque voie de réception, lesdits filtres égaliseur large bande étant obtenus par les étapes suivantes : 10 - Construire, pour chaque voie de réception, un vecteur de corrélation comprenant les points de corrélation obtenus simultanément en sortir desdits corrélateurs, - Déterminer, à partir de la fonction de corrélation du code local avec le signal reçu, estimée, pour chaque voie de réception, à partir 15 des vecteurs de corrélation, la voie de réception, appelée voie de référence, pour laquelle le signal reçu est le plus en retard par rapport au code local, - Déterminer, pour chaque voie de réception sauf la voie de référence, le filtre égaliseur large bande en calculant la 20 déconvolution du vecteur de corrélation de ladite voie de réception par le vecteur de corrélation de la voie de référence. Selon une variante de réalisation, la méthode selon l'invention comprend au préalable la recherche et la retenue, parmi les signaux reçus par ledit récepteur et provenant de plusieurs satellites, de ceux qui ne sont 25 pas affectés par la présence de trajets multiples du signal et la sélection d'un signal parmi ceux retenus. Selon une variante de réalisation, la méthode selon l'invention comprend la recherche et la sélection, parmi les signaux retenus, du signal émis par le satellite ayant l'élévation la plus élevée parmi les satellites 30 émetteurs des signaux reçus par ledit récepteur. 3047568 4 Selon une variante de réalisation, la méthode selon l'invention comprend la recherche et la sélection, parmi les signaux retenus, du signal qui présente le rapport signal à bruit le plus élevé. Selon une variante de réalisation, la méthode selon l'invention 5 comprend une étape de filtrage des sorties desdits corrélateurs pour diminuer le bruit de mesure. Selon une variante de réalisation, la méthode selon l'invention comprend la détermination d'un filtre égaliseur bande étroite comprenant un coefficient complexe destiné à être appliqué au signal filtré par le filtre 10 égaliseur large bande sur chaque voie de réception, ledit filtre égaliseur bande étroite étant obtenu par les étapes suivantes : - Pour chaque signal reçu provenant d'un satellite distinct et pour chaque voie de réception, filtrer le vecteur de corrélation par le filtre égaliseur large bande, 15 - Déterminer, pour chaque signal reçu provenant d'un satellite distinct, une matrice d'intercorrélation spatiale Rk à partir des vecteurs de corrélation de chaque voie de réception, - Déterminer, pour chaque signal reçu provenant d'un satellite distinct, la direction d'arrivée théorique du signal reçu à partir de la 20 position du satellite, - Déterminer une fonction de coût Q représentative de l'erreur globale entre les directions d'arrivée estimées et les directions d'arrivée théoriques des signaux reçus, une direction d'arrivée estimée étant déterminée à partir de la matrice d'intercorrélation 25 spatiale Rk et du coefficient complexe du filtre égaliseur bande étroite appliqué au signal reçu, - Rechercher, pour chaque voie de réception, les valeurs des coefficients complexes qui permettent de minimiser la fonction de coût Q. 30 3047568 5 L'invention a aussi pour objet une méthode de calibrage d'un récepteur de signaux de navigation par satellite, ledit récepteur comprenant une pluralité de voies de réception comprenant chacune une antenne de réception, une pluralité de corrélateurs pour corréler le signal reçu par 5 l'antenne avec un code local retardé d'un retard paramétrable et un filtre égaliseur large bande, obtenu par l'exécution de la méthode de calibrage selon l'invention et appliqué entre l'antenne de réception et les corrélateurs de chaque voie de réception, ladite méthode de calibrage comprenant la détermination d'un filtre égaliseur bande étroite comprenant un coefficient 10 complexe destiné à être appliqué au signal filtré par le filtre égaliseur large bande sur chaque voie de réception, ledit filtre égaliseur bande étroite étant obtenu par les étapes suivantes : - Déterminer, pour chaque signal reçu provenant d'un satellite distinct, une matrice d'intercorrélation spatiale Rk à partir des 15 vecteurs de corrélation de chaque voie de réception, - Déterminer, pour chaque signal reçu provenant d'un satellite distinct, la direction d'arrivée théorique du signal reçu à partir de la position du satellite, - Déterminer une fonction de coût Q représentative de l'erreur 20 globale entre les directions d'arrivée estimées et les directions d'arrivée théoriques des signaux reçus, une direction d'arrivée estimée étant déterminée à partir de la matrice d'intercorrélation spatiale Rk et du coefficient complexe du filtre égaliseur bande étroite appliqué au signal reçu, 25 - Rechercher, pour chaque voie de réception, les valeurs des coefficients complexes qui permettent de minimiser la fonction de coût Q. Selon un aspect particulier de l'invention, la matrice d'intercorrélation spatiale Rk est moyennée sur plusieurs instants consécutifs. 3047568 6 Selon un aspect particulier de l'invention, les directions d'arrivée des signaux reçus sont transposées dans un repère lié aux antennes de réception à partir de l'attitude du récepteur. Selon un aspect particulier de l'invention, l'attitude du récepteur est 5 intégrée comme paramètre de la fonction de coût Q et la méthode comprend en outre la recherche de la valeur de l'attitude du récepteur qui permet de minimiser la fonction de coût Q. Selon un aspect particulier de l'invention, la fonction de coût Q dépend des vecteurs directionnels attjk,Ok) des antennes de réception du récepteur 10 calculés dans les directions d'arrivée théoriques des signaux reçus. Selon un aspect particulier de l'invention, la fonction de coût Q est obtenue en sommant plusieurs fonctions de coûts chacune représentative de l'erreur entre la direction d'arrivée estimée et la direction d'arrivée théorique d'un signal reçu provenant d'un satellite. 15 Selon un aspect particulier de l'invention, la fonction de coût Q est Ns donnée par la relation Q=EagdokrelRilT.a(9k,Ok) avec Ns le nombre k=1 de signaux reçus provenant de satellites distincts, Rk1 est l'inverse de la matrice d'intercorrélation spatiale, F est une matrice diagonale dont les composantes sont les coefficients complexes des filtres d'égalisation bande 20 étroite sur chaque voie de réception et H est l'opérateur transposé conjugué. L'invention a aussi pour objet un récepteur de signaux de navigation par satellite comprenant une pluralité de voies de réception comprenant chacune une antenne de réception, une pluralité de corrélateurs pour 25 corréler le signal reçu par l'antenne avec un code local retardé d'un retard paramétrable et un filtre égaliseur large bande, disposé entre l'antenne de réception et lesdits corrélateurs et obtenu par l'exécution de la méthode de calibrage selon un premier mode de réalisation de l'invention. Dans une variante particulière de réalisation, le récepteur selon 30 l'invention comprend en outre, sur chaque voie de réception, un filtre 3047568 7 égaliseur bande étroite à un coefficient complexe, disposé entre la sortie du filtre égaliseur large bande et l'entrée desdits corrélateurs et obtenu par l'exécution de la méthode de calibrage selon un second mode de réalisation de l'invention. 5 Dans une variante particulière de réalisation, le récepteur selon l'invention comprend en outre, sur chaque voie de réception, au moins une ligne à retard connectée en sortie desdits corrélateurs. L'invention a encore pour objet un programme d'ordinateur comportant des instructions pour l'exécution de la méthode de calibrage d'un 10 récepteur de signaux de navigation par satellite selon l'invention, lorsque le programme est exécuté par un processeur. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux 15 dessins annexés qui représentent : La figure 1, un schéma d'un récepteur GNSS multi-voies, La figure 2, un organigramme décrivant les étapes de mise en oeuvre d'une méthode de calibrage selon un premier mode de réalisation de l'invention, 20 - La figure 3, un organigramme décrivant les étapes de mise en oeuvre d'une méthode de calibrage selon un second mode de réalisation de l'invention, La figure 4, un schéma d'un récepteur GNSS multi-voies calibré selon l'invention. 25 La figure 1 représente un schéma d'un exemple de récepteur GNSS multi-voies pour lequel l'invention s'applique. Un tel récepteur est décrit en détail dans la demande internationale de brevet publiée sous le numéro WO 2012/025306. 30 Il comprend plusieurs voies de réception comprenant chacune une antenne élémentaire Al ,A2,...An conçue pour capter un ou plusieurs signaux 3047568 8 de radio-navigation émis par un ou plusieurs satellites. Chaque voie de réception comporte également une chaine de traitement radio-fréquence (non représentée sur la figure 1) comprenant notamment un ou plusieurs filtres de mise en forme, un ou plusieurs amplificateurs ainsi qu'un 5 convertisseur analogique numérique. Les signaux analogiques échantillonnés sont ensuite fournis à plusieurs corrélateur(s) C11,C12,..Cn1,Cn2. Un corrélateur a pour rôle de réaliser une corrélation entre le signal reçu et une version locale du code d'étalement utilisé pour émettre le signal de radio-navigation. On rappelle que les signaux de radio- 10 navigation sont des signaux à étalement de spectre, par exemple des signaux utilisant la technique d'étalement de spectre CDMA (Code Division Multiple Access) en ce qui concerne les systèmes GPS et Galileo. Un code spécifique est associé à chaque satellite de la constellation. Plusieurs corrélateurs C11,C12 sont utilisés pour former un banc de 15 corrélateurs sur chaque voie de réception. Chaque corrélateur du banc utilise une version retardée du code local, le retard étant multiple de la période d'échantillonnage du signal. Un corrélateur fournit en sortie des points de corrélation et les différents points de corrélation fournis par le banc de corrélateurs permettent 20 de former une fonction de corrélation. La recherche du pic (ou maximum) de cette fonction de corrélation permet de synchroniser le signal reçu sur le code local et de déterminer ainsi une information de positionnement par le biais de plusieurs mesures de pseudo-distances entre plusieurs satellites et le récepteur. 25 Chaque corrélateur peut également être relié en sortie à une ou plusieurs lignes à retard R11m,R111,R12m,R121,.. Rn1m,Rn11,Rn2m. L'écart introduit par chaque ligne à retard correspond à un multiple de la durée d'intégration d'un corrélateur. Les lignes à retard permettent de 30 sauvegarder les résultats de corrélation sur une durée égale au nombre de lignes à retard multiplié par la durée d'intégration. 3047568 9 Les différents résultats de corrélation sont enfin fournis à un module de calcul 11P qui implémente un algorithme de traitement de signal permettant de déterminer une information de positionnement. 5 Un objectif de l'invention est de proposer une méthode de calibrage d'un récepteur multi-voies du type de celui décrit à la figure 1. La figure 2 représente un organigramme détaillant les étapes de mise en oeuvre de la méthode de calibrage selon un premier mode de réalisation de l'invention qui comprend la détermination d'un filtre égaliseur large bande 10 pour chaque voie de réception. Un filtre égaliseur est déterminé pour chaque voie de réception pour être positionné derrière chaque antenne et derrière la chaine de réception RF comme indiqué sur la figure 4. Le filtre égaliseur large bande déterminé selon l'invention a pour 15 objectif d'égaliser les fluctuations et différences de temps de groupe, de gain et de phase, dans la bande des signaux utiles, entre les différentes voies de réception. Selon une première étape 201, on construit, pour chaque voie de réception, un vecteur de corrélation dont chaque composante est égale à un 20 point de corrélation obtenu en sorti d'un corrélateur d'un banc de corrélateurs. Le vecteur de corrélation comprend ainsi les sorties des corrélateurs à un même instant. On rappelle que chaque corrélateur du banc réalise une corrélation avec un code local retardé d'un multiple de la période d'échantillonnage. 25 La première étape 201 est réalisée pour chaque voie de réception en considérant les signaux reçus par le récepteur pour un ou plusieurs satellites. Dans une étape optionnelle 202, les vecteurs de corrélation sont analysés pour identifier les signaux qui sont affectés par des multitrajets 30 provenant de réflexions du signal émis par un satellite sur des obstacles, par exemple des obstacles urbains. 304 756 8 10 Plus précisément, si n signaux émis par n satellites différents sont reçus par le récepteur, on analyse les vecteurs de corrélation calculés pour ces n signaux afin de déterminer si certain de ces signaux sont affectés par des multitrajets. 5 Différentes méthodes sont envisageables pour identifier la présence de multitrajets. De façon générale, il est possible de comparer la forme de la fonction de corrélation (déterminée à partir d'un ou plusieurs vecteurs de corrélation) à un modèle théorique. Si la forme de la fonction de corrélation mesurée diffère significativement du modèle théorique, cela peut signifier 10 que le signal est affecté par des multitrajets. A l'issue de l'étape 202, les n signaux reçus sont classés en deux groupes, un premier groupe de m signaux (m inférieur ou égal à n) non affectés par des multitrajets et un second groupe de n-m signaux affectés par des multitrajets. Les signaux du second groupe sont par la suite éliminés et 15 ne sont plus pris en compte pour les traitements suivants. L'étape 202 est optionnelle dans la mesure où il est possible de mettre en oeuvre l'invention en considérant qu'aucun signal reçu n'est affecté par la présence de multitrajets, par exemple parce que l'environnement du récepteur n'est pas susceptible de provoquer des réflexions du signal sur des 20 obstacles. La méthode se poursuit ensuite avec deux étapes alternatives 203,204 optionnelles qui ont pour objectif de ne retenir qu'un seul signal parmi les m signaux sélectionnés à l'étape 202. 25 L'étape 203 consiste à déterminer, pour chacun des m signaux retenus, l'angle d'élévation du satellite émetteur puis à retenir uniquement le signal émis par le satellite ayant l'élévation la plus élevée. L'angle d'élévation peut être déterminé à partir de la position du satellite obtenue à l'aide d'éphémérides. 30 Alternativement l'étape 204 consiste à retenir le signal qui présente le rapport signal à bruit le plus élevé parmi les m signaux. Le rapport signal à 3047568 11 bruit peut être estimé par différentes méthodes de traitement de signal, par exemple à partir de la méthode dite des moments ou de la méthode « narrow wideband power ratio » appliquées en sortie du corrélateur « prompt ». Le document « Carrier-to-Noise Algorithms », M. Petovello, pp.20-27, Inside 5 GNSS Jan/Feb 2010 donne plusieurs exemples d'algorithmes de mesure du rapport signal à bruit dans un récepteur GNSS. Si les étapes 202,203,204 ne sont pas exécutées, la solution la plus simple consiste à sélectionner un signal parmi les n signaux disponibles de 10 façon aléatoire. Dans une étape 205, appliquée au signal retenu parmi les n signaux disponibles, on évalue, à partir des différents vecteurs de corrélation déterminés à l'étape 201, quelle est la voie de réception qui présente le 15 temps de propagation de groupe le plus important, autrement dit la voie de réception sur laquelle le signal reçu est le plus retardé. Pour cela, on compare les retards estimés sur les différentes voies de réception entre le signal reçu et le code généré localement. Une façon d'estimer ce retard consiste à mesurer, sur chaque voie de réception, l'instant temporel associé 20 au pic (ou maximum) de la fonction de corrélation puis à comparer les différentes mesures afin d'identifier celle qui correspond au retard le plus important. La voie de réception ainsi identifiée est désignée par la suite voie de référence. 25 Dans une étape 206, on détermine ensuite, pour chaque voie de réception, sauf la voie de référence, un filtre égaliseur large bande. Ce filtre a pour fonction de compenser les différences de retard, gain et phase, dans la bande du signal utile, entre chaque voie de réception et la voie de référence qui correspond à la référence de temps de propagation de 30 groupe, gain et phase. Ainsi, dans le récepteur calibré selon l'invention tel que décrit à la figure 4, la voie de référence (voie 2 sur la figure 4) ne 3047568 12 comporte pas de filtre égaliseur tandis que toutes les autres voies comportent un filtre égaliseur positionné entre la chaine RF et les corrélateurs. Les coefficients du filtre égaliseur sont déterminés en réalisant une 5 opération de déconvolution du vecteur de corrélation de la voie de réception considérée par le vecteur de corrélation de la voie de référence. En appliquant les filtres égaliseurs ainsi déterminés en amont des corrélateurs, on obtient, en sortie des corrélateurs, des vecteurs de corrélation sensiblement identiques sur chaque voie de réception. 10 Si on note h() la réponse du filtre égaliseur à déterminer, C() la réponse d'un banc de corrélateurs et Sr le signal reçu pris en sortie de chaine RF, alors le vecteur de corrélation Vc obtenu en appliquant le filtre égaliseur entre la chaine RF et le banc de corrélateurs est Vc= C(h(Sr)) = h(C(Sr)) = h(V'e) où V'c est le vecteur de corrélation obtenu sans appliquer le 15 filtre égaliseur. Si h(X)=X* V'c -1 Vcref ou * désigne une opération de convolution, alors on obtient Vc= Vcref sur chaque voie de réception égalisée. Dans une variante de réalisation, les signaux reçus sur chaque voie peuvent être filtrés par un filtre supplémentaire de sorte à augmenter le 20 rapport signal à bruit. La figure 3 représente un organigramme détaillant les étapes de mise en oeuvre de la méthode de calibrage selon un second mode de réalisation de l'invention qui comprend la détermination d'un filtre égaliseur bande 25 étroite pour chaque voie de réception. L'application de la méthode de calibrage large bande décrite ci-dessus permet d'égaliser les désappairages entre les voies. L'invention propose également une méthode de calibrage bande étroite qui est appliquée en supplément de la méthode de calibrage large bande et qui permet d'égaliser 30 les résidus de gain et phase statique entre les différentes voies de réception. 304 756 8 13 La méthode de calibrage bande étroite selon l'invention s'applique à des signaux mesurés sur un récepteur qui a été, au préalable, calibré par la méthode de calibrage large bande selon l'invention. Ainsi, la méthode décrite à l'appui de la figure 3 peut être appliquée à 5 des mesures de corrélation réalisées sur un récepteur du type de celui décrit à la figure 1 mais comprenant en outre les filtres égaliseurs obtenus par application de la méthode de calibrage large bande décrite à la figure 2. Alternativement, il est possible d'utiliser les mêmes mesures de corrélation exploitées pour la méthode de calibrage large bande, c'est-à-dire 10 des mesures réalisées sur un récepteur du type de la figure 1 qui ne comporte aucun filtre égaliseur, puis de simuler l'effet de ces filtres égaliseurs sur les mesures de corrélation afin de générer de nouvelles mesures de corrélation « filtrées » qui peuvent être utilisées comme valeurs d'entrée de la méthode de calibrage bande étroite. 15 La méthode de calibrage bande étroite selon l'invention comprend une première étape 301 identique à la première étape 201 de la méthode de calibrage large bande décrite à la figure 2. Autrement dit, cette étape 301 consiste à construire, pour chaque voie 20 de réception, un vecteur de corrélation dont chaque composante est égale à un point de corrélation obtenu en sortie d'un corrélateur d'un banc de corrélateurs. La première étape 301 est réalisée pour chaque voie de réception en considérant les signaux reçus par le récepteur pour un ou plusieurs satellites. Comme indiqué, les vecteurs de corrélation sont obtenus 25 par mesure sur un récepteur comprenant des filtres égaliseurs large bande ou à partir de mesures de corrélation réalisées sur un récepteur non calibré sur lesquelles l'influence des filtres égaliseurs large bande est simulée. Une étape 302, identique à l'étape 202 de la méthode de calibrage 30 large bande, peut ensuite être appliquée optionnellement. 3047568 14 Selon cette étape 302, les vecteurs de corrélation sont analysés pour identifier les signaux qui sont affectés par des multitrajets provenant de réflexions du signal émis par un satellite sur des obstacles, par exemple des obstacles urbains. 5 Différentes méthodes sont envisageables pour identifier la présence de multitrajets. De façon générale, il est possible de comparer la forme de la fonction de corrélation (déterminée à partir d'un ou plusieurs vecteurs de corrélation) à un modèle théorique. Si la forme de la fonction de corrélation mesurée diffère significativement du modèle théorique, cela peut signifier 10 que le signal est affecté par des multitrajets. A l'issue de l'étape 302, on retient les signaux qui ne sont pas affectés par des multitrajets. On note Ns le nombre de ces signaux retenus. Ce nombre est au moins égal à deux. 15 Dans une étape 303, on construit, pour chaque signal retenu indicé par k, une matrice d'intercorrélation spatiale Rk du récepteur. La matrice d'intercorrélation spatiale Rk peut être déterminée à partir des vecteurs de corrélation des différentes voies de réception. Plus précisément, chaque composante de la matrice d'intercorrélation peut 20 correspondre au produit scalaire entre deux vecteurs de corrélation (ou entre un sous-ensemble des composantes des vecteurs de corrélation) associés à deux voies de réception différentes. Afin de diminuer le niveau de bruit impactant les valeurs de la matrice d'intercorrélation et d'obtenir une matrice de rang plein, on peut exploiter les 25 lignes à retard du récepteur GNSS situées derrière les corrélateurs (cf figure 1) pour calculer plusieurs matrices d'intercorrélation à différents instants successifs puis calculer la moyenne des matrices d'intercorrélation déterminées. 30 Dans une étape 304, on détermine, pour chaque signal associé à un satellite émetteur différent, la direction théorique d'arrivée du signal, 3047568 15 autrement dit la direction théorique entre le repère attaché au réseau d'antenne du récepteur et le satellite. Cette direction théorique peut être déterminée par un couple d'angle d'azimut et d'angle d'élévation (ijk,Ok ) ou par tout autre couple d'angles permettant d'identifier une direction d'arrivée 5 d'un signal. La direction théorique peut être déterminée à partir des coordonnées du satellite qui sont obtenues grâce, par exemple, à des éphémérides. La détermination de la direction théorique d'arrivée du signal nécessite également une transposition de la direction déterminée à partir des coordonnées du satellite, qui sont exprimées dans un repère fixe par rapport 10 à la Terre, vers un repère lié au récepteur et plus précisément au réseau d'antennes que comporte le récepteur. Ce changement de repère nécessite une estimée de l'attitude du réseau d'antennes. L'attitude peut être exprimée à l'aide de trois angles. Selon une variante de réalisation de l'invention, dans le cas où 15 l'attitude du réseau d'antennes n'est pas connue, il est possible de prendre plusieurs hypothèses a priori sur la valeur de l'attitude et d'intégrer cette donnée comme une variable de la fonction de coût déterminée à l'étape suivante 305. Ce point sera explicité plus en détail par la suite. 20 Dans une étape 305, on détermine une fonction de coût Q représentative de l'erreur entre les directions estimées d'arrivée des signaux et les directions théoriques d'arrivée des signaux. La fonction de coût Q peut être déterminée comme la somme de fonctions de coûts Qk associées à chaque signal provenant d'un satellite 25 d'indice k. Une fonction de coût Qk dépend de la matrice d'intercorrélation spatiale calculée à l'étape 304 et plus précisément de l'inverse de cette matrice. Elle dépend également des coefficients des filtres égaliseurs bande étroite déterminés in fine par la méthode de calibrage bande étroite. Ces 30 coefficients peuvent être regroupés dans une matrice diagonale I" . Un coefficient complexe d'égalisation est ainsi destiné à être appliqué sur 3047568 16 chaque voie de réception, préférentiellement en amont des corrélateurs. Les coefficients complexes d'égalisation, formant le filtre d'égalisation bande étroite selon l'invention, visent à corriger les erreurs commises sur l'estimation des directions d'arrivées des signaux du fait de défauts résiduels 5 dans l'appairage des voies de réception. Enfin, la fonction de coût Qk dépend du vecteur directionnel a(ê k , k) du réseau antennaire du récepteur GNSS qui comprend les amplitudes et les déphasages relatifs des signaux reçus par chaque antenne du réseau par rapport à une antenne du réseau prise comme référence, en fonction de la direction du signal reçu. 10 La fonction de coût Qk peut prendre la forme suivante : Qk = k Ok r eiRkina(êk/Ok où l'exposant H désigne l'opérateur transposé conjugué. La fonction de coût globale Q est donnée par la relation Ns Q= , -\ll 1'111E1;1145k Ok)' k=1 Toute variante de calcul de la fonction de coût introduite ci-dessus 15 entre dans le périmètre de l'invention dans la mesure où celle-ci est bien représentative de l'erreur entre les directions estimées d'arrivée des signaux et les directions théoriques d'arrivée des signaux. Les directions estimées d'arrivée des signaux sont déterminées en prenant en compte un filtre d'égalisation bande étroite appliqué à chaque voie de réception et 20 comprenant un seul coefficient complexe permettant de compenser le gain et la phase du signal reçu sur chaque voie de sorte à corriger la direction estimée d'arrivée du signal par rapport à la direction théorique déterminée à partir des positions connues des satellites. La fonction de coût Q déterminée à l'issue de l'étape 305 a pour 25 inconnues les coefficients d'égalisation bande étroite appliqués à chaque voie de réception et formalisés par la matrice diagonale F . Selon une variante de réalisation de l'invention introduite précédemment, si l'attitude du réseau antennaire du récepteur GNSS n'est pas connue et qu'il n'est pas possible de la déterminer, alors le vecteur 3047568 17 directionnel a(êk,Ok) peut inclure l'attitude du réseau antennaire en tant que paramètre. Dans une étape 306, on recherche les valeurs des coefficients 5 d'égalisation bande étroite qui permettent de minimiser la fonction de coût Q déterminée à l'issue de l'étape 305. En réécrivant la matrice diagonale F sous la forme d'un vecteur colonne K, la fonction de coût peut se mettre sous la forme plus avantageuse : 10 Ns Q = KH[ EdiaeaWk , k\\H Ri-idiag(a(ijk ,Ok)) ic= K H ZK k=1 Le problème de minimisation de la fonction de coût Q revient donc à minimiser le terme KHZK sous la contrainte K1 =1 où KI désigne le premier 15 terme du vecteur K. Cette contrainte est nécessaire pour éviter la solution triviale K=0. En résolvant ce problème d'optimisation quadratique à contrainte linéaire, on obtient la solution : û = avec w = [1 0 - - - Of wHz-lw Qmin = 1 20 wHz-`w Selon une variante de réalisation de l'invention introduite précédemment, si l'attitude du réseau antennaire du récepteur n'est pas connue, elle est ajoutée aux inconnues de la fonction de coût et est considérée comme paramètre supplémentaire dans la minimisation de la 25 fonction de coût. Par exemple, l'étape 306 décrite ci-dessus peut être répétée pour plusieurs hypothèses d'attitude du réseau antennaire prises dans une plage de variation d'attitude prédéterminée. Pour chaque hypothèse d'attitude, le minimum Qmin de la fonction de coût est déterminé 3047568 18 ainsi que le vecteur d'égalisation K qui permet d'obtenir ce minimum. Ensuite, une étape supplémentaire de recherche de minimum parmi les minima Qmin est effectuée et on retient in fine le minimum Qmin de la fonction de coût parmi les minima calculés pour chaque hypothèse d'attitude et le 5 vecteur d'égalisation associé. Le procédé de calibrage selon le premier ou le second mode de réalisation de l'invention peut être mis en oeuvre en tant que programme d'ordinateur comportant des instructions pour son exécution. Le programme 10 d'ordinateur selon l'invention prend en entrée des points de mesures réalisées en sorties des corrélateurs d'un récepteur GNSS tel qu'explicité ci-dessus. Le programme d'ordinateur selon l'invention fournit en sortie les coefficients des filtres d'égalisation large bande et/ou des filtres d'égalisation bande étroite. Le programme d'ordinateur peut être enregistré sur un support 15 d'enregistrement lisible par un processeur. Le support peut être électronique, magnétique, optique ou électromagnétique. La figure 4 représente un schéma d'un récepteur GNSS multi-voies 20 modifié selon l'invention. Le récepteur décrit à la figure 4 comporte les mêmes éléments que le récepteur de la figure 1. Selon une première variante de l'invention, il comporte en outre, sur chaque voie de réception sauf une, un filtre d'égalisation large bande FIRi,FIRn. Sur l'exemple de la figure 4, la voie de 25 réception numéro deux est la voie de référence. Les filtres d'égalisation large bande FIRi,FIRn sont déterminés par la méthode de calibrage large bande selon l'invention telle que décrite à l'appui de la figure 2. Les filtres d'égalisation large bande FIRi,FIRn sont positionnés, comme indiqué sur la figure 4, entre les chaines RF et les corrélateurs de 30 chaque voie de réception. 3047568 19 Les filtres d'égalisation large bande FIRi ,FIRn peuvent être mis en oeuvre au moyen de filtres à réponse impulsionnelle finie comportant un nombre de coefficients inférieur au nombre de corrélateurs. Selon une seconde variante de l'invention, le récepteur GNSS multi5 voies comporte en outre, non seulement les filtres d'égalisation large bande FIRi,FIRn décrits ci-dessus mais également un filtre d'égalisation bande étroite G1,G2,Gn sur chaque voie de réception (y compris la voie de référence). Les filtres d'égalisation bande étroite Gi ,G2,Gn sont déterminés par la méthode de calibrage bande étroite selon l'invention telle que décrite à 10 l'appui de la figure 3. Les filtres d'égalisation bande étroite Gi ,G2,Gil peuvent être positionnés, comme indiqué sur la figure 4, entre les filtres d'égalisation large bande et les corrélateurs de chaque voie de réception. Alternativement, ils peuvent également être positionnés en sortie des bancs de corrélateurs de 15 chaque voie. Les filtres d'égalisation bande étroite Gi ,G2,Gr, sont mis en oeuvre au moyen d'un multiplieur complexe, c'est-à-dire qui affecte un gain d'amplitude et un déphasage aux signaux reçus sur chaque voie. Les différents filtres d'égalisation ainsi que les différents éléments qui 20 composent le récepteur GNSS selon l'invention peuvent être mis en oeuvre par des moyens matériels et/ou logiciels. En particulier, le récepteur GNSS dans son ensemble ou chaque élément du récepteur peut être implémenté par un dispositif comprenant un processeur et une mémoire. Le processeur peut être un processeur générique, un processeur spécifique, un circuit 25 intégré propre à une application (connu aussi sous le nom anglais d'ASIC pour « Application-Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). 30 REVENDICATIONS1. Méthode de calibrage d'un récepteur de signaux de navigation par satellites, ledit récepteur comprenant une pluralité de voies de réception comprenant chacune une antenne de réception (A1,A2..An) et une pluralité de corrélateurs (C11,C12,..Cn1,Cn2) pour corréler le signal reçu par l'antenne avec un code local retardé d'un retard paramétrable, ladite méthode de calibrage comprenant la détermination d'au moins un filtre égaliseur large bande (FIR1,FIR2,..FIRn) destiné à être appliqué entre l'antenne de réception et les corrélateurs de chaque voie de réception, lesdits filtres égaliseur large bande étant obtenus par les étapes suivantes : - Construire (201), pour chaque voie de réception, un vecteur de corrélation comprenant les points de corrélation obtenus simultanément en sortir desdits corrélateurs, - Déterminer (205), à partir de la fonction de corrélation du code local avec le signal reçu, estimée, pour chaque voie de réception, à partir des vecteurs de corrélation, la voie de réception, appelée voie de référence, pour laquelle le signal reçu est le plus en retard par rapport au code local, - Déterminer (206), pour chaque voie de réception sauf la voie de référence, le filtre égaliseur large bande en calculant la déconvolution du vecteur de corrélation de ladite voie de réception par le vecteur de corrélation de la voie de référence. 2. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon la revendication 1 comprenant au préalable la recherche et la retenue (202), parmi les signaux reçus par ledit récepteur et provenant de plusieurs satellites, de ceux qui ne sont pas affectés par la présence de trajets multiples du signal et la sélection d'un signal parmi ceux retenus. 3047568 21 3. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon la revendication 2 comprenant la recherche et la sélection (203), parmi les signaux retenus, du signal émis par le satellite ayant 5 l'élévation la plus élevée parmi les satellites émetteurs des signaux reçus par ledit récepteur. 4. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon la revendication 2 comprenant la recherche et la sélection 10 (204), parmi les signaux retenus, du signal qui présente le rapport signal à bruit le plus élevé. 5. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une des revendications précédentes comprenant une 15 étape de filtrage des sorties desdits corrélateurs pour diminuer le bruit de mesure. 6. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une des revendications précédentes comprenant la 20 détermination d'un filtre égaliseur bande étroite (G1,G2,...Gn) comprenant un coefficient complexe destiné à être appliqué au signal filtré par le filtre égaliseur large bande (FIR1,FIR2,..FIRn) sur chaque voie de réception, ledit filtre égaliseur bande étroite étant obtenu par les étapes suivantes : - Pour chaque signal reçu provenant d'un satellite distinct et pour 25 chaque voie de réception, filtrer le vecteur de corrélation par le filtre égaliseur large bande, - Déterminer (303), pour chaque signal reçu provenant d'un satellite distinct, une matrice d'intercorrélation spatiale Rk à partir des vecteurs de corrélation de chaque voie de réception, 30 - Déterminer (304), pour chaque signal reçu provenant d'un satellite distinct, la direction d'arrivée théorique du signal reçu à partir de la position du satellite, 3047568 22 - Déterminer (305) une fonction de coût Q représentative de l'erreur globale entre les directions d'arrivée estimées et les directions d'arrivée théoriques des signaux reçus, une direction d'arrivée estimée étant déterminée à partir de la matrice d'intercorrélation 5 spatiale Rk et du coefficient complexe du filtre égaliseur bande étroite appliqué au signal reçu, - Rechercher (306), pour chaque voie de réception, les valeurs des coefficients complexes qui permettent de minimiser la fonction de coût Q. 10 7. Méthode de calibrage d'un récepteur de signaux de navigation par satellite, ledit récepteur comprenant une pluralité de voies de réception comprenant chacune une antenne de réception (Al ,A2..An), une pluralité de corrélateurs (C11,C12,..Cn1,Cn2) pour corréler le signal reçu par 15 l'antenne avec un code local retardé d'un retard paramétrable et un filtre égaliseur large bande (FIR1,FIR2,..FIRn), obtenu par l'exécution de la méthode de calibrage selon l'une des revendications 1 à 5 et appliqué entre l'antenne de réception et les corrélateurs de chaque voie de réception, ladite méthode de calibrage comprenant la détermination d'un 20 filtre égaliseur bande étroite (G1,G2,...Gn) comprenant un coefficient complexe destiné à être appliqué au signal filtré par le filtre égaliseur large bande (FIR1,FIR2,..FIRn) sur chaque voie de réception, ledit filtre égaliseur bande étroite étant obtenu par les étapes suivantes : - Déterminer (303), pour chaque signal reçu provenant d'un satellite 25 distinct, une matrice d'intercorrélation spatiale Rk à partir des vecteurs de corrélation de chaque voie de réception, - Déterminer (304), pour chaque signal reçu provenant d'un satellite distinct, la direction d'arrivée théorique du signal reçu à partir de la position du satellite, 30 - Déterminer (305) une fonction de coût Q représentative de l'erreur globale entre les directions d'arrivée estimées et les directions 3047568 23 d'arrivée théoriques des signaux reçus, une direction d'arrivée estimée étant déterminée à partir de la matrice d'intercorrélation spatiale Rk et du coefficient complexe du filtre égaliseur bande étroite appliqué au signal reçu, 5 - Rechercher (306), pour chaque voie de réception, les valeurs des coefficients complexes qui permettent de minimiser la fonction de coût Q. 8. Méthode de calibrage d'un récepteur de signaux de navigation par 10 satellite selon l'une des revendications 6 ou 7 dans laquelle la matrice d'intercorrélation spatiale Rk est moyennée sur plusieurs instants consécutifs. 9. Méthode de calibrage d'un récepteur de signaux de navigation par 15 satellite selon l'une des revendications 6 à 8 dans laquelle les directions d'arrivée des signaux reçus sont transposées dans un repère lié aux antennes de réception à partir de l'attitude du récepteur. 10.Méthode de calibrage d'un récepteur de signaux de navigation par 20 satellite selon la revendication 9 dans laquelle l'attitude du récepteur est intégrée comme paramètre de la fonction de coût Q et la méthode comprend en outre la recherche (306) de la valeur de l'attitude du récepteur qui permet de minimiser la fonction de coût Q. 25 11.Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une des revendications 6 à 10 dans laquelle la fonction de coût Q dépend des vecteurs directionnels a(Bk , y'à k) des antennes de réception du récepteur calculés dans les directions d'arrivée théoriques des signaux reçus. 30 12. Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une des revendications 6 à 11 dans laquelle la fonction de 3047568 24 coût Q est obtenue en sommant plusieurs fonctions de coûts chacune représentative de l'erreur entre la direction d'arrivée estimée et la direction d'arrivée théorique d'un signal reçu provenant d'un satellite. 5 13.Méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une des revendications 6 à 12 dans laquelle la fonction de Ns coût Q est donnée par la relation Q = Zek,OkrTHR-kir.a(Bk,Ok) avec k=1 NS le nombre de signaux reçus provenant de satellites distincts, Rk-1 est l'inverse de la matrice d'intercorrélation spatiale, T est une matrice io diagonale dont les composantes sont les coefficients complexes des filtres d'égalisation bande étroite sur chaque voie de réception et H est l'opérateur transposé conjugué. 14. Récepteur de signaux de navigation par satellite comprenant une pluralité 15 de voies de réception comprenant chacune une antenne de réception (A1,A2..An), une pluralité de corrélateurs (C11,C12,..Cn1,Cn2) pour corréler le signal reçu par l'antenne avec un code local retardé d'un retard paramétrable et un filtre égaliseur large bande (FIR1,FIR2,..FIRn), disposé entre l'antenne de réception et lesdits corrélateurs et obtenu par 20 l'exécution de la méthode de calibrage selon l'une des revendications 1 à 5. 15. Récepteur de signaux de navigation par satellite selon la revendication 14 comprenant en outre, sur chaque voie de réception, un filtre égaliseur 25 bande étroite (G1,G2...Gn) à un coefficient complexe, disposé entre la sortie du filtre égaliseur large bande et l'entrée desdits corrélateurs et obtenu par l'exécution de la méthode de calibrage selon l'une des revendications 6 à 13. 30 16. Récepteur de signaux de navigation par satellite selon la revendication 15 comprenant en outre, sur chaque voie de réception, au moins une ligne à 3047568 25 retard (R11m,R111,R12m,R121...Rn1m,Rn11,Rn2m,Rn21) connectée en sortie desdits corrélateurs. 17. Programme d'ordinateur comportant des instructions pour l'exécution de 5 la méthode de calibrage d'un récepteur de signaux de navigation par satellite selon l'une quelconque des revendications 1 à 13, lorsque le programme est exécuté par un processeur.