FR 3050830 A1 20171103 FR 1653763 A 20160427 1 Domaine Technique [0001] De manière générale, l'invention concerne un simulateur de constellation de satellites de radionavigation, plus particulièrement un procédé et un dispositif de calibration d'un tel simulateur. Arrière-plan technologique [0002] Les signaux de radionavigation (GNSS) émis par les satellites d'un système de positionnement par satellites (GNSS, de l'anglais : « Global Navigation Satellite System ») se présentent sous la forme d'une porteuse radiofréquence (RF) modulée par une forme d'onde d'étalement contenant un code binaire pseudo-aléatoire (« code d'étalement »). Accessoirement, les signaux GNSS transportent des données utiles (notamment le message de navigation contenant les données d'éphémérides) sous forme d'une séquence numérique (à un rythme nettement moins élevé que le code pseudo-aléatoire) modulée en plus sur la porteuse. De manière simplifiée, un récepteur de signaux de radionavigation (récepteur GNSS) détermine sa position par trilatération en mesurant le temps (et donc les distances) de parcours de plusieurs signaux GNSS en provenance de différents satellites GNSS. [0003] Un simulateur de constellation de satellites de radionavigation (simulateur GNSS) émule l'environnement radiofréquence tel qu'il se présenterait à un récepteur GNSS à un endroit, à une date et à une heure donnés. En d'autres mots, un simulateur GNSS produit, en fonction de paramètres d'entrée incluant le lieu, la date et l'heure de la réception, l'ensemble des signaux radiofréquence (RF) qu'un récepteur GNSS recevrait de la constellation de satellites de radionavigation réelle s'il était placé dans ledit lieu auxdites date et heure. L'endroit de réception simulé n'a pas besoin d'être invariable mais peut changer dans le temps dans le contexte d'un scénario avec un récepteur mobile. Le simulateur modélise notamment le mouvement des satellites de radionavigation (satellites GNSS) et le temps de parcours des signaux entre les différents satellites et l'endroit de réception. Selon le perfectionnement du simulateur, il peut également simuler des dégradations des signaux GNSS, comme, par exemple, des perturbations des signaux dues à l'ionosphère, la troposphère, des dérives d'horloges des satellites GNSS, des effets multi-trajets, etc. 3050830 2 [0004] Les simulateurs GNSS sont utilisés pour tester des récepteurs GNSS, des prototypes et des composants de récepteurs dans des conditions parfaitement maîtrisées. [0005] Dans un simulateur, avant la combinaison des signaux sur la sortie RF, chaque 5 signal GNSS est généré dans un canal séparé. Ces canaux comprennent des composants électroniques et RF qui sont physiquement indépendants. De la génération séparée des signaux GNSS peuvent résulter de légères différences entre les signaux GNSS produits. En particulier, les signaux GNSS peuvent être affectés de biais (temporels) différents (biais inter-canaux) non maîtrisés. Ces biais inconnus 10 limitent la précision de la simulation, car ils se répercutent dans les mesures réalisées par le récepteur GNSS à tester. [0006] Le document GB 2524340 A décrit un procédé pour aligner entre eux des signaux de radionavigation d'un simulateur GNSS. Dans le procédé, un premier biais est déterminé entre un premier signal de code délivré par un premier canal et un signal 15 de code de référence. D'autres biais sont déterminés entre des signaux de code délivrés par d'autres canaux et la continuation du signal de code de référence. Les biais sont ensuite appliqués aux différents canaux par une logique de contrôle de sorte à aligner les canaux entre eux. Problème technique 20 [0007] Un objectif d'un aspect de la présente invention est d'améliorer les solutions existantes pour calibrer les simulateurs GNSS. Description générale de l'invention [0008] Un premier aspect de l'invention se rapporte à un procédé (de calibration) pour déterminer un biais interne d'un simulateur de constellation de satellites de 25 radionavigation (simulateur GNSS), par rapport à un signal d'horloge de référence, le simulateur GNSS mettant à disposition en sortie un signal de radionavigation (signal GNSS) en bande RF. Le procédé comprend l'exécution d'un scénario par le simulateur dans lequel le biais théorique entre le signal de radionavigation et le signal d'horloge de référence est connu. Le procédé comprend aussi la numérisation du signal de 30 radionavigation (éventuellement après transposition vers une fréquence 3050830 3 intermédiaire), la génération d'une réplique du code du signal de radionavigation, l'évaluation d'une fonction de corrélation entre le signal de radionavigation et la réplique de code, l'identification d'un pic de corrélation de la fonction de corrélation et la déduction du biais interne du simulateur. L'évaluation de la fonction de corrélation 5 est réalisée par transformation de Fourier et transformation de Fourier inverse, c.-à-d. par détermination des transformées de Fourier du signal de radionavigation et de la réplique de code et par transformation de Fourier inverse du produit des transformées de Fourier. La déduction du biais interne du simulateur se fait sur base du pic de corrélation identifié et du biais théorique. 10 [0009] Dans le contexte de ce document, on entend par « scénario » l'ensemble des paramètres, décrivant une situation et éventuellement son évolution dans le temps, fournis au simulateur GNSS afin que celui-ci puisse calculer et générer l'environnement RF du récepteur comme s'il se trouvait dans la situation décrite. Les options de paramétrage disponibles peuvent dépendre du type de simulateur. 15 Toutefois, comme il s'agit d'équipements dédiés aux spécialistes, la liberté de définir des scénarios est très large. [0010] On notera que le scénario de la simulation peut être choisi de sorte à ce que le biais théorique soit égal à une certaine valeur, p.ex. 0, un nombre entier de périodes de code ou un nombre entier de bribes (en anglais : « chips ») de code. Dans ce 20 contexte, il convient toutefois de noter que le biais théorique n'a pas besoin d'être connu d'avance et/ou au moment de l'exécution du scénario. Pour l'invention, il est en effet suffisant que le scénario de simulation permette de déduire le biais théorique par des considérations purement théoriques et que le biais théorique soit connu au moment du calcul du biais interne. 25 [0011] Comme énoncé plus haut, le simulateur permet de modéliser des environnements radiofréquences, de manière contrôlée et retraçable, pour un scénario donné, au niveau d'un récepteur. Pour un scénario donné, il est donc possible de déterminer et/ou connaître le biais engendré par la propagation d'un signal GNSS vers le récepteur. Ce biais représente le « biais théorique » mentionné ci-dessus. Le biais 30 mesurable dans un test de récepteur se compose du biais théorique, du biais interne du simulateur et du biais interne du récepteur. Dès lors, afin de pouvoir calibrer le 3050830 4 récepteur avec précision, il faut connaître non seulement le biais théorique (dépendant des paramètres de la simulation) mais également les biais internes du simulateur. Une application intéressante est dès lors d'utiliser un simulateur GNSS calibré selon le procédé de la présente invention pour calibrer un récepteur GNSS. Un tel procédé de 5 calibrage serait grandement facilité par rapport aux méthodes de calibrage actuellement employées. L'utilisation d'un simulateur pour calibrer un récepteur offre des avantages par rapport à une calibration en conditions réelles, c'est-à-dire sur des signaux de radionavigation provenant d'une ou plusieurs constellations de satellites de radionavigation. Ces avantages sont, entre autres, la répétabilité, le très bon 10 contrôle des conditions de test, la très bonne précision de la calibration en situation de laboratoire, un coût réduit pour la phase de test du récepteur. [0012] La numérisation du signal GNSS fourni par le simulateur GNSS se fait préférablement à une fréquence d'échantillonnage comprise entre 0,1 GHz et 1 GHz, plus préférablement à une fréquence d'échantillonnage comprise entre 0,1 GHz et 15 0,5 GHz, encore plus préférablement à une fréquence d'échantillonnage comprise entre 0,1 GHz et 0,2 GHz. Il est à noter que même si de telles fréquences d'échantillonnage ne permettent pas d'échantillonner correctement la porteuse du signal GNSS au sens du théorème d'échantillonnage de Nyquist-Shannon, elles permettent par contre d'échantillonner correctement le code d'étalement ou les 20 données compris dans le signal GNSS. [0013] Préférablement, après lancement du scénario de simulation, le signal GNSS est numérisé pendant un temps inférieur ou égal à une période du code compris dans le signal GNSS. De cette manière, on s'assure que la fonction de corrélation entre le signal GNSS et la réplique de code possède un pic de corrélation unique (un maximum 25 global) sur le temps de corrélation. Le temps d'échantillonnage ainsi que la fréquence d'échantillonnage peuvent varier dynamiquement (c.-à-d. au cours du temps) selon les besoins. À titre d'exemple, lorsqu'un signal de radionavigation comprend un code GPS C/A (comprenant une séquence de 1023 bribes de code à une cadence de 1,023 Mcps (1 Mcps = 1 « mega-chips per second » = 106 bribes par seconde)), le temps de 30 numérisation du signal de radionavigation est inférieur ou égal à 1 ms. A une fréquence d'échantillonnage de 1 GHz, le code C/A est échantillonné 1 million de fois par période de code soit environ 1000 fois par bribe. 3050830 5 [0014] Le signal de référence peut être un signal fourni par le simulateur, p.ex. un signal 1PPS (pour « 1 pulse per second » en anglais, signal à une impulsion par seconde). L'instant de référence pour le calcul d'un biais par rapport à un signal de référence peut être un flanc montant ou un flanc descendant du signal de référence. 5 Alternativement, le signal de référence peut être fourni par un système extérieur, par exemple, par une horloge externe. L'horloge du simulateur GNSS est de préférence asservie à une horloge externe, employée également par le dispositif de calibration utilisé pour exécuter le procédé. [0015] Une réplique (locale) du code compris dans le signal GNSS est générée et 10 échantillonnée (ou alternativement directement fournie sous forme d'échantillons) en vue de son utilisation pour calculer la fonction de corrélation. La réplique peut, par exemple, être échantillonnée à la même fréquence d'échantillonnage que le signal GNSS. Le temps d'échantillonnage de la réplique est, de préférence, égal à la période du code compris dans le signal GNSS. 15 [0016] L'identification d'un pic de corrélation de la fonction de corrélation comprend, de préférence, l'évaluation du module de la fonction de corrélation entre la réplique de code et le signal GNSS provenant du simulateur. Le calcul du module de la fonction de corrélation est recommandé si des mesures spécifiques n'ont pas été prises pour faire en sorte que la fonction de corrélation soit réelle (ou purement imaginaire). Si on 20 note C(r) la fonction de corrélation entre la réplique de code r(t) et le signal GNSS s(t) provenant du simulateur, où t est le temps, on peut écrire : C(r) = fI r(t - x)s(t)dt = FTI( FT[r(t)]conj(FT[s(t)])), (1) où I est l'intervalle d'intégration, la réplique de code r(t) est décalée d'un temps 'r (biais) par rapport à s(t), FT[f(t)] représente la transformation de Fourier d'une fonction f(t) et 25 FTI[f(t)] représente la transformation de Fourier inverse d'une fonction f(t). Chaque point de la fonction de corrélation correspond donc à l'intégration de la fonction résultante de la multiplication du signal GNSS avec une version décalée en temps de la réplique de code. De façon imagée, la fonction de corrélation C(r) représente la similarité entre la réplique de code r(t) et le signal GNSS s(t) en fonction du décalage 30 en temps -c. Le module de la fonction de corrélation est maximal lorsque le code du signal GNSS se superpose à la réplique de code. La zone autour du maximum du 3050830 6 module de la fonction de corrélation a, généralement, une forme de triangle. Selon le type signal GNSS le module de la fonction de corrélation peut comprendre une succession de triangles. C'est par exemple le cas si la modulation est de type forme d'onde binaire décalée (« Binary Offset Carrier » (BOC) en anglais). Le maximum du 5 module de la fonction d'onde a éventuellement une pointe écrêtée due aux imperfections du signal GNSS. Ces imperfections peuvent être dues à, par exemple, un temps de montée ou un temps de descente non négligeable par rapport à une bribe du signal de code, des suroscillations (« overshoot » en anglais), etc. [0017] Le maximum du module de la fonction de corrélation ou pic de corrélation peut 10 être identifié par détection du maximum des valeurs échantillonnant cette fonction respectivement son module. Le décalage en temps, pour lequel la valeur du module de la fonction de corrélation est maximale, peut ainsi être identifié. Optionnellement, afin d'augmenter la précision de l'identification du pic de corrélation, une des méthodes décrites dans la suite du document peut être utilisée pour détecter le maximum de la 15 fonction de corrélation ou de son module. Toutes ces méthodes opèrent directement sur la fonction de corrélation ou son module. Il sera apprécié qu'il est possible combiner les résultats de ces méthodes (par exemple par calcul de la moyenne de ceux-ci). [0018] Une première méthode comprend un ajustement (en anglais « fit ») par une 20 droite de chaque flanc du pic de corrélation. La zone d'ajustement de la fonction de corrélation ou de son module (c.-à-d. les points de la fonction de corrélation ou de son module considérés pour l'ajustement) de chaque droite est configurée pour contenir uniquement la zone linéaire (les flancs) autour du pic de corrélation. L'intersection de ces deux droites identifie à la fois la valeur du maximum théorique (c.-à-d. sans 25 écrêtage) de la fonction de corrélation ou de son module et la valeur du décalage en temps (biais) associé. Le pic de corrélation est ainsi localisé précisément. Optionnellement des contraintes peuvent être imposées lors de l'ajustement. Ces contraintes peuvent comprendre, par exemple une contrainte de symétrie des flancs du pic de corrélation. 30 [0019] Une seconde méthode comprend l'ajustement de la fonction de corrélation ou de son module par la fonction théorique (c.-à-d. la fonction de corrélation attendue 3050830 7 pour la modulation en question). La localisation du pic de corrélation peut être aidée par la définition d'une zone de recherche du maximum de la fonction de corrélation ou de son module ou par la définition de valeurs initiales proches des valeurs réelles du pic de corrélation pour l'ajustement de la fonction de corrélation. L'ajustement fournit 5 le décalage en temps associé au maximum de la fonction de corrélation respectivement de son module. [0020] Une troisième méthode d'identification du pic de corrélation comprend la localisation du maximum par une approche du type avance-retard (en anglais « earlylate »). Dans ce cas, on place deux points séparés d'une distance donnée de part et 10 d'autre de la zone culminante de la fonction de corrélation de sorte à ce qu'ils soient à la même hauteur (valeur sur l'axe des C(T)) et on détermine l'abscisse (valeur sur l'axe des -c) du milieu du segment reliant les deux points. On répète ce procédé pour différentes distances entre les deux points. En fin de compte, on identifie comme l'abscisse du pic de corrélation la valeur moyenne des abscisses de tous les milieux 15 trouvés. [0021] Une quatrième méthode comprend le calcul des coefficients de Fourier intervenant dans la décomposition de Fourier de la fonction de corrélation ou de son module. La localisation du pic de corrélation comprend la détermination du maximum, et le décalage associé, de la fonction de corrélation résultante de la décomposition de 20 Fourier. [0022] Selon un mode de réalisation de l'invention, le décalage Doppler éventuel du signal GNSS est déterminé. Le signal GNSS peut effectivement comprendre un décalage en fréquence dû, par exemple, à un mouvement relatif entre le satellite (simulé) émettant des signaux GNSS et le récepteur, à une dérive de l'horloge interne 25 du satellite (simulé) ou à une dérive entre l'horloge interne du simulateur par rapport à l'horloge interne du dispositif de calibration utilisé pour exécuter le procédé. Le décalage Doppler du signal GNSS en sortie du simulateur peut, par exemple, être déterminé itérativement. Alternativement, pour éviter une méthode itérative, l'espace fréquentiel peut être directement découpé finement. 30 [0023] La déduction du biais interne du simulateur se fait sur base du pic de corrélation identifié ainsi que du biais théorique. Le décalage dans le temps entre le 3050830 8 signal GNSS en sortie du simulateur et la réplique de code correspond à la somme du biais théorique et du biais interne du simulateur. Le biais interne du simulateur peut donc être trouvé par la différence entre le décalage mesuré et le biais théorique. [0024] Selon un mode de réalisation de l'invention, le simulateur modélise plusieurs 5 signaux GNSS. Ces signaux GNSS sont générés sur des canaux séparés du simulateur. Chaque canal a son propre biais interne. Ces biais internes peuvent être déterminés de manière séquentielle ou en parallèle. Les plusieurs signaux GNSS générés par le simulateur GNSS peuvent alors être alignés sur base des plusieurs biais internes. Si la détermination des différents biais internes repose sur des 10 scénarios faisant intervenir plusieurs biais théoriques, lors de l'alignement des signaux GNSS, on en tiendra compte. [0025] Un second aspect de l'invention se rapporte à un dispositif (de calibration) pour déterminer un biais interne, par rapport à un signal d'horloge de référence, d'un simulateur GNSS mettant à disposition, en sortie, un signal GNSS résultant de 15 l'exécution un scénario dans lequel le biais théorique entre le signal de radionavigation et le signal d'horloge de référence est connu. Le dispositif comprend une entrée pour le signal de radionavigation issu par le simulateur GNSS, un numériseur pour numériser le signal GNSS, un générateur de répliques de code configuré pour produire une réplique du code du signal GNSS, un processeur, une sortie configurée pour 20 fournir le biais interne déduit par le processeur. Le processeur est configuré pour évaluer une fonction de corrélation entre le signal de radionavigation et la réplique de code par transformation de Fourier et transformation de Fourier inverse, identifier le pic de corrélation de la fonction de corrélation et déduire du pic de corrélation identifié et du biais théorique le biais interne entre le signal de radionavigation et le signal 25 d'horloge de référence. Le processeur peut être configuré pour lesdites tâches par programmation ou du fait qu'il s'agit d'un matériel spécialement conçu, p.ex. d'un circuit électronique spécialement adapté (« application-specific integrated circuit » en anglais ou ASIC) ou d'un circuit logique programmable (p.ex. un « field-programmable gate array » ou FPGA (réseau de portes programmables in situ)). Au cas où le 30 processeur est configuré pour certaines tâches par programmation, il comprend une mémoire ou il a accès à une mémoire dans laquelle le programme correspondant est enregistré de manière interprétable ou exécutable par le processeur. 3050830 9 [0026] Préférablement, le générateur de réplique de code ainsi que le processeur sont configuré pour déterminer un décalage Doppler du signal de radionavigation. [0027] Le dispositif pour déterminer un biais entre un signal de référence et un signal de radionavigation provenant d'un simulateur de constellation de satellites de 5 radionavigation est préférablement un dispositif en boucle ouverte, c.-à-d. un dispositif sans boucle de rétroaction afin d'aligner, sur base des un ou plusieurs biais, les un ou plusieurs signaux de radionavigation provenant du simulateur. Brève description des dessins [0028] D'autres particularités et caractéristiques de l'invention ressortiront de la 10 description détaillée de certains modes de réalisation avantageux présentés ci- dessous, à titre d'illustration, avec référence aux dessins annexés qui montrent : Fig. 1: un schéma d'un montage de calibration d'un simulateur GNSS à l'aide d'un dispositif de calibration ; Fig. 2: un chronogramme d'un signal d'horloge de référence 1PPS et d'un signal GPS 15 L1 C/A généré par un simulateur GNSS ; Fig. 3: un chronogramme d'un signal 1PPS et d'un code d'étalement d'un signal GNSS généré en sortie par un simulateur GNSS ; Fig. 4: un ordinogramme comprenant les étapes d'un procédé de détermination d'un biais interne d'un simulateur GNSS ; 20 Fig. 5: une illustration de la détermination d'une fonction de corrélation entre un signal GNSS généré par le simulateur GNSS et une réplique locale de code pseudoaléatoire ; Fig. 6: le graphe d'une fonction de corrélation obtenue ; Fig. 7: une illustration d'un ajustement d'une fonction de corrélation par deux droites ; 25 Fig. 8: une illustration d'un ajustement d'une fonction de corrélation par une fonction triangle ; Fig. 9: une illustration de la détermination du pic de corrélation par une approche de type avance-retard ; 3050830 10 Fig. 10: une illustration des différentes fonctions de corrélation pour différents décalages Doppler ; Fig. 11: une partie d'un signal GNSS en sortie du simulateur GNSS observé par un oscilloscope après amplification. 5 Description détaillée de modes de réalisation de l'invention [0029] La Fig. 1 illustre un dispositif de calibration 10 pour déterminer un biais interne 12 d'un simulateur GNSS 16 mettant à disposition en sortie un signal GNSS 18 et un signal d'horloge de référence 14. Le signal d'horloge de référence est, dans le cas décrit, un signal 1PPS. Le dispositif de calibration 10 possède une première entrée 10 pour le signal d'horloge de référence 14 et une deuxième entrée pour le signal GNSS 18. Le dispositif de calibration 10 possède aussi une sortie afin de mettre à disposition le biais interne déterminé. Afin d'éviter ou au moins limiter des erreurs liées, par exemple, à une éventuelle dérive entre les horloges internes du simulateur GNSS 16 et du dispositif de calibration 10, un oscillateur 20 fournit une fréquence de référence 15 ou un signal d'horloge de référence communs au simulateur GNSS 16 et au dispositif de calibration 10. [0030] Le simulateur GNSS 16 comprend différents canaux, réalisés par des composants électriques et/ou électroniques séparés, qui sont responsables de la génération de signaux GNSS. La génération des signaux GNSS par des canaux 20 indépendants peut résulter dans des différences entre les signaux GNSS produits. En particulier, les signaux GNSS peuvent être affectés de biais différents qui limitent, s'ils restent inconnus, la précision des simulations réalisées par le simulateur GNSS 16. [0031] Les signaux GNSS sont générés en accord avec un scénario de simulation comprenant, entre autres, des trajectoires pour les satellites GNSS et le récepteur 25 GNSS. Le scénario exécuté par le simulateur GNSS 16 étant configurable virtuellement à souhait, il est possible de modéliser des situations réalistes ou même des situations impossibles à réaliser, p.ex. un scénario dans lequel les satellites GNSS se trouvent sur d'autres que leurs orbites réelles. Le simulateur GNSS est configuré pour envoyer, lors de l'exécution d'une simulation, et en fonction du scénario choisi 30 par l'utilisateur, des instructions aux canaux afin que ceux-ci génèrent les signaux GNSS correspondants. 3050830 11 [0032] La Fig. 2 montre un exemple de signaux produits en sortie du simulateur GNSS 16. Le signal d'horloge de référence (14, voir Fig. 1) prend la forme d'un signal 1PPS 22, c.-à-d. un signal possédant une impulsion toutes les secondes mais constant par ailleurs. Idéalement, le flanc montant du signal 1PPS est aligné au début d'une 5 période de code pour chaque signal GNSS pour une distance de propagation nulle. Or, pour chaque canal du simulateur GNSS 16, il existe un biais entre le flanc montant de l'impulsion du signal 1PPS et la transition de code et ce biais n'est pas identique (sauf après une éventuelle calibration préalable et sauf accident) pour tous les signaux GNSS. Pour déterminer ce biais, on pourrait penser à simplement mesurer le temps entre le flanc montant de l'impulsion du signal 1PPS et la transition de code à l'aide d'un oscilloscope. Or, pour ce faire, il faudrait d'abord amplifier le signal GNSS concerné, ce qui nécessiterait l'introduction d'un amplificateur à haut gain. Comme le code pseudo-aléatoire module une porteuse, l'observation qu'on peut faire sur un oscilloscope est une forme d'onde oscillant rapidement, dont l'amplitude est 15 localement réduite à l'instant de la transition de code (voir Fig. 11). Cette approche ne permet pas de déterminer le biais concerné à moins de 200 ps près, précision considérée insuffisante pour des besoins de calibration. [0033] Le procédé selon l'invention s'appuie sur une technique de corrélation pour détecter le biais interne. Cette technique permet de s'affranchir de l'amplification du 20 signal GNSS ou au moins d'utiliser un amplificateur moins puissant. [0034] Pour l'illustration d'un mode de réalisation préféré de l'invention, le signal GNSS en sortie du simulateur est représenté par un signal GPS C/A 24 émis dans la bande de fréquences L1. On comprendra toutefois que le procédé s'applique à n'importe quel autre signal GNSS. La période du code C/A 24 étant de 1 ms, un pic du 25 signal 1PPS 22 interviendra toutes les mille périodes du code. [0035] La Fig. 3 illustre un biais entre le signal 1PPS 22 (graphe (a)) et le code pseudo-aléatoire 26 d'un signal GNSS en sortie du simulateur (graphe (b)). Le biais 28 entre le flanc montant du signal 1PPS et le début du code 26 se compose d'un biais théorique 32, dû à la propagation (simulée) du signal GNSS entre le satellite GNSS et 30 le récepteur, et un biais interne 30, spécifique au canal du simulateur ayant généré le signal. Les durées représentées sur la Fig. 3 ne sont pas nécessairement 3050830 12 représentatives des proportions entre le bais interne 30 et le biais théorique 32 formant le biais total 28. Le biais théorique 32 peut être déterminé sur base du scénario utilisé. Il est avantageux, lors de la calibration du simulateur, de sélectionner un scénario, dans lequel le biais théorique est constant (p.ex. 0). On utilise donc de préférence un 5 scénario dans lequel le mouvement relatif entre le satellite GNSS et le récepteur est nul et dans lequel les satellites GNSS sont tous sur une orbite géostationnaire, etc. [0036] Pour pouvoir visualiser plus facilement le biais interne, on peut concevoir un scénario dans lequel, en l'absence de biais interne, le flanc montant du signal 1PPS serait aligné avec le début d'une bribe du code (p.ex. la première bribe de la période 10 de code). Un tel scénario est maintenant décrit plus en détail. Afin de rendre le scénario le plus simple possible, on considère une propagation des signaux GNSS dans le vide. On note d la distance entre le satellite GNSS et le récepteur. Pour que, en absence de biais interne, le flanc montant du signal 1PPS soit aligné avec le début d'une bribe de code, il faut que la condition 15 d/c = N.Tcode k.Tbribe, avec N, k des entiers 0 soit satisfaite (c représentant la vitesse de la lumière dans le vide, Tcode la période de code et Tbribe la durée d'une bribe de code.) Pour remplir cette condition dans le scénario de la simulation, peut placer le satellite GNSS sur une orbite géostationnaire (p.ex. à 0°E à une altitude altGE0) et le récepteur à une distance d juste en-dessous du 20 satellite (p.ex. à la position 0°E, 0°N et à l'altitude altGEo-d). Les valeurs de N et de k peuvent être choisies librement ; p.ex. N = 0 et k = 0 place le récepteur à l'endroit de l'émetteur. Avec un tel choix de la distance d, tout décalage entre le flanc montant du signal 1PPS et le début d'une bribe du code est directement imputable au biais interne du simulateur (voir graphe (c) de la Fig. 3). 25 [0037] La Fig. 4 illustre le procédé selon un aspect de l'invention pour déterminer le biais interne d'un simulateur GNSS sous forme d'un ordinogramme. La première étape S10 comprend la calibration (en phase, en amplitude, en TPG, etc.) d'éléments électriques et/ou électroniques entre le simulateur et le dispositif de calibration (par exemple des câbles électriques 32, voir Fig. 1) et éventuellement d'éléments 30 électriques et/ou électroniques du dispositif de calibration (par exemple d'amplificateurs, etc.). Cette calibration est réalisée préférentiellement grâce à un 3050830 13 analyseur de réseau vectoriel (« vector network analyser » en anglais). À l'étape suivante S12, le scénario à utiliser pour la détermination du biais interne est généré et encodé dans le simulateur. Avantageusement, le scénario décrit précédemment est utilisé. 5 [0038] Le scénario est ensuite exécuté par le simulateur GNSS à l'étape S14 et le dispositif de calibration numérise le signal en sortie du simulateur. La numérisation peut se faire, p.ex., à des fréquences d'échantillonnage comprise entre 0,1 GHz et 1 GHz. La porteuse du signal GNSS sera ainsi sous-échantillonnée et le code du signal GNSS sera, par contre, correctement échantillonné (c.-à-d. à une fréquence 10 suffisamment élevée pour respecter les conditions du théorème d'échantillonnage de Nyquist-Shannon). En vue d'une détermination aisée du biais interne par la corrélation entre le signal GNSS et le signal 1PPS, une réplique du code comprise dans le signal GNSS est générée et échantillonnée à la même fréquence par le dispositif à étape S16. La génération de la réplique de code est faite en accord avec la valeur de k 15 choisie à l'étape S12. De cette façon, un pic de corrélation situé en 0 sur l'axe des abscisses signifie que le code du signal GNSS est aligné avec la réplique de code générée. Le signal GNSS provenant du simulateur GNSS est échantillonné pendant un temps inférieur ou égal à une période du signal afin qu'un seul pic de corrélation soit présent dans la fonction de corrélation. Le signal GNSS numérisé et échantillonné 20 est complété par des valeurs zéros (« zero-padding » en anglais) afin d'obtenir le même nombre de points entre la réplique du code échantillonnée et le signal GNSS échantillonné. [0039] Une fonction de corrélation est calculée à l'étape S18 selon l'Eq. 1. Afin de s'affranchir de la phase de la porteuse, le module de la fonction de corrélation sera 25 considéré dans la suite. La Fig. 5 représente la corrélation entre code pseudo-aléatoire 26 du signal GNSS et la réplique de code 31 du signal GNSS. Pour chaque valeur de décalage temporel T de la réplique de code 31 du signal GNSS, une valeur de la fonction de corrélation est associée. Autrement dit, pour chaque translation donnée (décalage temporel -r) de la réplique de code (désignée par la flèche pleine à la Fig. 5), 30 une valeur de la fonction de corrélation est déterminée. La valeur du maximum de la fonction de corrélation donne le décalage temporel entre ceux-ci et est liée à la valeur recherchée du biais interne 30 du simulateur GNSS. Il est à noter que la largeur du pic 3050830 14 de la fonction de corrélation représentée à la Fig. 5 n'est pas représentative d'un cas réel. [0040] Un exemple de fonction de corrélation est donné à la Fig. 6. Idéalement, la fonction de corrélation est de forme triangulaire 36. Cependant, les signaux GNSS en 5 sortie du simulateur GNSS sont entachés d'imperfections induites principalement par la traversée des signaux dans les composants constitutifs du simulateur GNSS, du dispositif de calibration et des éléments de connexion. Ces imperfections provoquent un écrêtement 34 du pic de corrélation visible à la Fig. 6. Les flancs 38 du pic de corrélation ne sont, par contre, pas (ou très peu) affectés par ces imperfections. Io [0041] Selon un mode de réalisation préféré de l'invention et avec référence à la Fig. 7, l'identification du pic de corrélation 40 de la fonction de corrélation 44 à l'étape S20 est avantageusement réalisée par un ajustement linéaire (c.-à-d. par une droite) de chaque flanc du pic de corrélation. Seuls les points compris dans la zone d'ajustement sont pris en compte pour réaliser l'ajustement. L'ajustement du pic de 15 corrélation peut être réalisé avec une contrainte de symétrie sur les deux droites (c.- à-d. mêmes valeurs absolues des pentes des droites) vu qu'il n'y a aucune source d'asymétrie, telle que des interférences de type trajets multiples (ou multi-trajets, « multipath » en anglais), dans le scénario considéré. L'intersection 40' de ces deux droites 42 identifie le pic de corrélation 40 et en particulier le décalage temporel 48 20 recherché. [0042] Selon un autre mode de réalisation de l'invention et avec référence à la Fig. 8, l'identification du pic de corrélation 50 de la fonction de corrélation 44 à l'étape S20 est réalisée par un ajustement du pic de corrélation par une fonction triangle 46. Afin d'identifier le pic de corrélation 50, une zone de recherche 49 est définie lors de 25 l'ajustement. Le pic du triangle 50' (le maximum de la fonction triangle résultante de l'ajustement) identifie le pic de corrélation 40 et son abscisse le décalage temporel 52. Optionnellement, lors de l'ajustement, un poids plus élevé peut être affecté aux points compris sur les flancs du pic de corrélation. [0043] Selon un troisième mode de réalisation de l'invention et avec référence à la 30 Fig. 9, l'identification du pic de corrélation 54 de la fonction de corrélation 56 à l'étape S20 est réalisée par une approche du type avance-retard. Selon cette approche, deux 3050830 15 points 58, 60 sont placés à gauche et à droite du pic de corrélation 54 à une même hauteur (c.-à-d. les deux points possèdent la même ordonnée). Le milieu 62 du segment 64 entre les deux points 58, 60 est déterminé et sauvegardé. Le procédé est répété pour un nombre prédéterminé de paires de points 66, 68, 70, 72 dans des zones 5 de sélection à gauche 74 et à droite 76 du pic de corrélation 54. Les milieux 62, 78, 80 des segments créés par les pairs de points et sauvegardés sont reportés dans un histogramme et la moyenne de ces milieux identifie le décalage temporel 82 recherché. [0044] Selon un quatrième mode de réalisation de l'invention, le pic de corrélation est identifié par décomposition de Fourier de la fonction de corrélation calculée. Après détermination des coefficients de Fourier associés à la décomposition de Fourier, le maximum de la décomposition est recherché, ce qui fournit ainsi le décalage temporel. [0045] Avec référence à la Fig. 4, l'étape suivante S22 du procédé déduit la valeur du biais interne à partir du décalage temporel trouvé à l'étape S20. Vu la définition de la 15 réplique de code dans l'exemple discuté en détail (biais théorique = 0), le décalage temporel est dans ce cas identique au biais interne recherché. [0046] Optionnellement, à l'étape S24, une autre détermination du biais interne suivant le même ou un autre scénario encodé dans le simulateur GNSS est possible. Le procédé entre ainsi dans une boucle pour déterminer le biais interne et peut, par 20 exemple, réaliser une moyenne des biais interne ainsi déterminés. [0047] Un fois le biais théorique et le biais interne du simulateur GNSS connus, un alignement du signal GNSS avec le signal 1PPS peut être réalisé grâce à, par exemple, une boucle de rétroaction ou, de préférence, par une correction précise en boucle ouverte. 25 [0048] Dans le cas où le simulateur GNSS modélise plusieurs signaux GNSS, plusieurs canaux de génération de signaux GNSS sont utilisés. Le procédé de détermination d'un biais interne d'un canal du simulateur GNSS peut être utilisé pour déterminer le biais interne de chaque canal du simulateur. Après la détermination des biais internes, un alignement des signaux GNSS au signal 1PPS est réalisé. 30 [0049] Optionnellement, le décalage Doppler éventuel du signal GNSS est déterminé. Le signal GNSS peut effectivement comprendre un décalage en fréquence dû à, par 3050830 16 exemple, un scénario comprenant un mouvement relatif entre le satellite GNSS et le récepteur ou à une dérive entre l'horloge interne du simulateur GNSS par rapport à l'horloge interne du dispositif de calibration utilisé pour exécuter le procédé. Le décalage Doppler du signal GNSS peut, par exemple, être déterminé itérativement. 5 Dans ce cas et avec référence à la Fig. 10, l'espace fréquentiel (ou espace Doppler), correspondant à l'espace des décalages Doppler possibles du signal GNSS, est découpé en plusieurs classes Doppler 84 (« Doppler bins » en anglais). Dans un premier temps, l'espace fréquentiel est découpé grossièrement (c.-à-d. des classes Doppler assez larges telles que représentées à la Fig. 10). Pour chaque classe Doppler 84, - le signal GNSS est multiplié par une porteuse dont la fréquence Doppler est celle associée à la classe Doppler 84 ; - la fonction de corrélation entre le signal GNSS multiplié et la réplique du code compris dans le signal est déterminée selon l'Eq. 1 ; 15 - le module de la fonction de corrélation 86 est calculé ; - le maximum global 88 de la fonction de corrélation (c.-à-d. selon le décalage temporel et selon la fréquence Doppler) de toutes les fonctions de corrélations est identifié Le procédé itératif sauvegarde le décalage en temps -cm= et la fréquence Doppler 20 fD,max associés au maximum global. A la prochaine itération, la classe Doppler pour laquelle le maximum global a été obtenu est redécoupée en classes plus fines et le procédé itératif recommence. Le procédé itératif est répété jusqu'à l'obtention d'une précision voulue pour le décalage Doppler et le décalage en temps. Alternativement à la méthode itérative, un découpage fin de l'espace fréquentiel peut être envisagé dès 25 le départ afin de ne pas recourir à une méthode itérative. [0050] Alors que des modes de réalisation particuliers viennent d'être décrits en détail, l'homme du métier appréciera que diverses modifications et alternatives à ceux-là puissent être développées à la lumière de l'enseignement global apporté par la présente divulgation de l'invention. Par conséquent, les agencements et/ou procédés 30 spécifiques décrits ci-dedans sont censés être donnés uniquement à titre d'illustration, sans intention de limiter la portée de l'invention. REVENDICATIONS1. Un procédé pour déterminer un biais interne, par rapport à un signal d'horloge de référence, d'un simulateur de constellation de satellites de radionavigation mettant à disposition en sortie un signal de radionavigation, caractérisé par l'exécution d'un scénario par ledit simulateur dans lequel le biais théorique entre le signal de radionavigation et le signal d'horloge de référence est connu ; la numérisation dudit signal de radionavigation ; la génération d'une réplique du code du signal de radionavigation ; l'évaluation d'une fonction de corrélation entre le signal de radionavigation et ladite réplique de code par transformation de Fourier et transformation de Fourier inverse ; l'identification d'un pic de corrélation de la fonction de corrélation ; et la déduction, sur base du pic de corrélation identifié et du biais théorique, dudit biais interne du simulateur. 2. Le procédé selon la revendication 1, dans lequel la numérisation dudit signal de radionavigation se fait à une fréquence d'échantillonnage comprise entre 0,1 GHz et 1 GHz. 3. Le procédé selon l'une quelconque des revendications 1 à 2, dans lequel ledit signal de radionavigation est numérisé pendant un temps inférieur ou égal à une période du signal de code dudit signal de radionavigation. 4. Le procédé selon l'une quelconque des revendications 1 à 3, comprenant la détermination d'un décalage Doppler dudit signal de radionavigation. 5. Le procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'identification d'un pic de corrélation comprend l'ajustement par une droite de chaque flanc du pic de corrélation et la localisation de l'intersection de ces deux droites. 6. Le procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'identification d'un pic de corrélation comprend l'ajustement du pic de corrélation par une fonction théorique et la localisation du maximum global de la fonction théorique résultante de l'ajustement. 3050830 18 Le procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'identification d'un pic de corrélation comprend la localisation du maximum par une approche du type avance-retard. Le procédé selon l'une quelconque des revendications 1 à 4, dans lequel l'identification d'un pic de corrélation comprend la décomposition de Fourier de la fonction de corrélation et la localisation du maximum de la fonction résultante de la décomposition de Fourier. Le procédé selon l'une quelconque des revendications 1 à 8, dans lequel le simulateur de constellation de satellites de radionavigation met à disposition en sortie plusieurs signaux de radionavigation, le procédé comprenant la détermination de plusieurs biais internes du simulateur de constellation de satellites de radionavigation. Le procédé selon la revendication 9, dans lequel les plusieurs signaux de radionavigation provenant du simulateur de constellation de satellites de radionavigation sont alignés sur base des plusieurs biais internes. Un dispositif pour déterminer un biais interne, par rapport à un signal d'horloge de référence, d'un simulateur de constellation de satellites de radionavigation mettant à disposition, en sortie, un signal de radionavigation, le simulateur configuré pour exécuter un scénario dans lequel le biais théorique entre le signal 20 de radionavigation et le signal d'horloge de référence est connu, le dispositif étant caractérisé par une entrée pour ledit signal de radionavigation ; un numériseur pour numériser ledit signal de radionavigation ; un générateur de répliques de code configuré pour produire une réplique du code 25 dudit signal de radionavigation ; un processeur configuré pour o évaluer une fonction de corrélation entre le signal de radionavigation et ladite réplique de code par transformée de Fourier et transformée de Fourier inverse ; 30 o identifier le pic de corrélation du module de la fonction de corrélation ;7. 8. 5 9. 10 10. 15 11. 3050830 19 o déduire du pic de corrélation identifié et du biais théorique ledit biais interne entre ledit signal de radionavigation et le signal d'horloge de référence ; et une sortie configurée pour fournir ledit biais interne déduit par le processeur. 5 12. Le dispositif selon la revendication 11, dans lequel le générateur de répliques de code et le processeur sont configurés pour déterminer un décalage Doppler du signal de radionavigation. 13. Le dispositif selon l'une quelconque des revendications 10 à 12, dans lequel le dispositif est un dispositif en boucle ouverte. 10