FR 3141581 A1 20240503 FR 2211177 A 20221027 L’invention concerne un procédé de transfert de temps vers un récepteur en lien radiofréquence avec au moins trois balises terrestres de transfert de temps synchronisées diffusant chacune un message horaire. Le procédé décrit comprend : une étape (201) de détermination d’une mesure de pseudo-distance Pdi associée à chaque balise de transfert de temps, une étape (202) de détermination d’une position initiale et d’un décalage temporel initial dudit récepteur à partir des mesures de pseudo-distances, une étape (203) de calcul d’un décalage temporel fin du récepteur, par un algorithme de résolution optimisant un critère de maximum de vraisemblance résolvant conjointement les positions des balises terrestres de transfert de temps, la position et le décalage temporel du récepteur. L’invention concerne également un équipement de radiocommunications mettant en œuvre le procédé, ainsi qu’un programme d’ordinateur et un support informatique contenant ledit programme. Figure pour l’abrégé : Fig. 2 Procédé de radiodiffusion de temps précis par un réseau de balises de positions imprécises L’invention se situe dans le domaine du transfert d’une information de temps par voie radiofréquence. Elle porte en particulier sur un procédé permettant de transmettre une information de temps précise à une plateforme mobile ou semi-mobile par des balises de radiodiffusion dont les positions sont connues avec une faible précision. Le partage, entre plusieurs équipements, d’une référence de temps commune de bonne précision peut s’avérer indispensable. C’est le cas de manière générale dès lors qu’il est nécessaire de synchroniser des équipements, comme par exemple pour les transmissions radio multiplexées en temps (TDMA, sigle anglais pour Time Division Multiple Access , ou accès multiple à répartition dans le temps) ou étalées en code (CDMA, sigle anglais pour Code Division Multiple Access , ou accès multiple à répartition en code), ou utilisant l’étalement par saut de fréquence (FHSS, sigle anglais pour Frequency Hopping Spread Spectrum ). La précision recherchée dépend de l’application visée, mais un ordre de grandeur standard de la précision recherché est la microseconde ou moins. En règle générale, les équipements cherchant à déterminer un temps précis font appel au positionnement GNSS (sigle anglais pour Global Navigation Satellite Systems , ou système global de navigation par satellite), qui permet de déterminer précisément une position et une heure. Cependant, dans certains cas, la diffusion d’une heure par un réseau de satellites n’est pas possible, par exemple lorsque les conditions de visibilité des satellites sont mauvaises, que l’on ne dispose pas de moyens permettant de recevoir des signaux des satellites, ou lorsque l’on ne souhaite pas s’appuyer sur un réseau GNSS. Il est alors possible de diffuser une information de temps précise par des voies filaires, mais cette solution n’est pas applicable dès lors que les équipements sont mobiles. Une autre solution consiste à transmettre un marquant horaire par trois balises terrestres de transfert de temps ou plus sur une voie radio. Dès lors que les balises sont précisément synchronisées (par des voies filaires ou des horloges de haute précision) et que leurs positions sont transmises ou connues de l’équipement cherchant à déterminer l’heure du réseau, il est possible de calculer une distance avec chaque balise et d’en déduire la position de l’équipement. Cette détermination se fait de manière comparable à un positionnement GNSS, par l’estimation de pseudo-distances avec chacune des balises, obtenues par mesure du temps de propagation, puis par trilatération à partir des pseudo-distances. Cette solution fonctionne dans le cas d’un réseau fixe de balises terrestres de transfert de temps dont la position est parfaitement connue. Ce n’est plus le cas lorsqu’il existe une imprécision élevée sur la position des balises de transfert de temps. Ce cas d’application est cependant fréquent lorsqu’il s’agit de déployer un réseau de manière provisoire sur un site en l’absence de service GNSS, ou dans le cas de réseaux dits « ad-hoc » (réseau sans fil décentralisé ne s’appuyant pas sur une architecture fixe de points d’accès). Dans ces réseaux, chaque nœud peut jouer le rôle de balise de transfert de temps, et contribue ainsi à la synchronisation des autres nœuds du réseau. Les nœuds pouvant être mobiles, leur position n’est pas nécessairement connue de manière très précise. La précision de l’heure diffusée peut alors ne pas atteindre les niveaux d’exigence souhaités, en particulier si les nœuds sont distants, même lorsque les balises de transfert de temps disposent d’horloges très précises. Jusqu’à maintenant, les mesures mises en œuvre pour pallier ce problème consistent à ajouter un poste d’erreur sur la position des balises de transfert de temps dans le calcul de la variance de l’estimateur de temps, à la manière d’un calcul de DOP (sigle anglais pour Dilution Of Precision) pour des liaisons GNSS, ce qui permet de prendre en compte une incertitude sur l’heure, mais ne permet pas d’améliorer la précision de l’heure diffusée. Un objet de l’invention est donc de décrire un procédé de transfert d’une information précise de temps dans un réseau de radiocommunications depuis des balises terrestres en prenant en compte le problème de la précision de leur positionnement A cet effet, la présente invention décrit un procédé de transfert de temps vers un récepteur en lien radiofréquence avec au moins trois balises terrestres de transfert de temps synchronisées diffusant chacune un message horaire comprenant un marquant horaire et une heure d’émission. Le procédé décrit comprend : une première étape de calcul d’un instant d’arrivée des messages horaires transmis par lesdites balises terrestres de transfert de temps, les positions desdites balises terrestres de transfert de temps étant connues par le récepteur ou transmises au récepteur, et de détermination d’une mesure de pseudo-distance Pd i associée à chaque balise de transfert de temps, une deuxième étape de détermination d’une position initiale et d’un décalage temporel initial dudit récepteur à partir des mesures de pseudo-distances calculées lors de la première étape, une troisième étape de calcul d’un décalage temporel fin dudit récepteur, par la mise en œuvre d’un algorithme de résolution optimisant un critère de maximum de vraisemblance utilisant les valeurs calculées à la deuxième étape et résolvant conjointement les positions des balises terrestres de transfert de temps, la position du récepteur et le décalage temporel du récepteur. Selon un mode de réalisation, la troisième étape comprend : la construction d’un vecteur d’état comprenant la position estimée du récepteur, le décalage temporel estimé du récepteur, et les positions estimées des balises terrestres de transfert de temps, la construction d’un vecteur d’observation y comprenant les positions connues ou transmises desdites au moins trois balises terrestres de transfert de temps et les pseudo-distances calculées lors de la première étape, la construction d’un vecteur d’observation linéarisé comprenant un écart entre les positions connues ou transmises desdites au moins trois balises terrestres de transfert de temps et leurs positions estimées, et un écart entre les pseudo-distances calculées lors de la première étape (201) et des pseudo-distances construites à partir des valeurs du vecteur d’état , la construction d’une matrice Jacobienne à partir du vecteur d’état et du vecteur d’observation tel que , et la résolution du système, en calculant . Avantageusement, le calcul de est itéré plusieurs fois, les valeurs du vecteur et de la matrice étant recalculées à chaque itération en fonction du vecteur déterminé à l’itération précédente. Selon différents modes de réalisation, l’algorithme utilisant un critère de maximum de vraisemblance mis en œuvre lors de la troisième étape est un algorithme de Gauss-Newton ou un algorithme de Levenberg-Marquardt. Avantageusement, la deuxième étape est mise en œuvre par un procédé d’estimation pseudo-linéaire de la position initiale et du décalage temporel initial dudit récepteur. Selon un mode de réalisation, l’estimation pseudo-linéaire comprend : le calcul d’une position initiale du récepteur, avec , où est une matrice construite à partir des mesures de pseudo-distances Pd i calculées lors de la première étape du procédé et des positions connues ou transmises des balises terrestres de transfert de temps, est un vecteur d’observation comprenant les mesures de pseudo-distances Pd i faites lors de la première étape du procédé, et est une pseudo-mesure construite à partir des mesures de pseudo-distance Pd i faites lors de la première étape du procédé et des positions connues ou transmises des balises, le calcul d’un décalage temporel initial tel que , avec N le nombre de balises terrestres de transfert de temps et un vecteur comprenant les positions connues ou transmises des balises terrestres de transfert de temps. Le procédé selon l’invention est particulièrement adapté à des cas d’usage où la position des balises terrestres de transfert de temps est connue de manière imparfaite. L’invention porte également sur un équipement de radiocommunications (ou récepteur RF) en lien radiofréquence avec au moins trois balises terrestres de transfert de temps synchronisées diffusant chacune un message horaire comprenant un marquant horaire et une heure d’émission. L’équipement de radiocommunications selon l’invention comprend une chaîne radio configurée pour recevoir et interpréter le message horaire transmis par les balises. Il comprend en outre des moyens de calculs configurés pour mettre en œuvre un procédé de transfert de temps selon l’une des revendications précédentes. Enfin, l’invention porte également sur un programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé de transfert de temps selon l’invention lorsque ledit programme est exécuté sur un ordinateur, et un support d’enregistrement lisible par un ordinateur sur lequel est enregistré un tel programme d’ordinateur. L’invention sera mieux comprise et d’autres caractéristiques, détails et avantages apparaîtront mieux à la lecture de la description qui suit, donnée à titre non limitatif, et grâce aux figures annexées, données à titre d’exemple. La représente une configuration d’équipements dans laquelle un procédé de transfert de temps d’un équipement radio selon l’invention peut être mis en œuvre. La représente schématiquement les étapes d’un procédé de transfert de temps selon un mode de réalisation de l’invention. La représente une configuration opérationnelle dans laquelle est mis en œuvre un procédé de transfert de temps selon l’invention. La représente une configuration d’équipements dans laquelle un procédé de transfert de temps selon l’invention peut être mis en œuvre. Il comprend un équipement 101 cherchant à se synchroniser temporellement avec au moins trois balises terrestres de transfert de temps 102 à 104 avec lesquelles il est en lien radio. Les balises terrestres de transfert de temps peuvent être mobiles ou temporairement fixes mais de positionnements approximatifs, de sorte que leurs positions ne sont connues que de manière imparfaite. Les balises sont synchronisées très précisément, par exemple en embarquant des horloges de haute précision comme des horloges atomiques ou en étant raccordées à un réseau filaire ou un réseau radio dédié implémentant un mécanisme de synchronisation précis. Par la suite, on considèrera que le réseau est plan, et que les inconnues portent uniquement sur des positions orthogonales x et y dans un plan horizontal. Seules trois balises de transfert de temps sont donc nécessaires pour permettre à un équipement radio de récupérer une information de temps précise. Cependant, l’invention s’applique également en dehors d’un plan, en ajoutant la dimension z orthogonale à x et y, auquel cas un minimum de quatre stations est requis. L’ajout de stations de transfert de temps supplémentaires permet d’améliorer la qualité de la détermination de l’heure, en lissant le bruit de mesure. Le procédé de transfert de temps selon l’invention consiste à prendre en compte explicitement les positions des balises de transfert de temps dans un estimateur mis en œuvre par le récepteur, afin d’estimer conjointement les positions précises des balises, la position du récepteur et son biais temporel par rapport au réseau de balises. Il suppose l’existence d’une forme d’onde radio adaptée au besoin de transfert de temps par radiodiffusion, c’est-à-dire la transmission par les balises d’un message horaire et d’un top horloge périodique, que le récepteur radio à synchroniser sait recevoir et mesurer précisément. Le message horaire contient l’instant d’émission du top horloge, défini selon le temps de référence du réseau de balises. La position de chaque balise (potentiellement imprécise) peut être transmise dans le message horaire, mémorisée dans le récepteur, ou transmise au récepteur par une liaison de données. Le procédé selon l’invention nécessite la mesure des instants d’arrivée (en anglais TOA pour Time of Arrival ) des marquants horaires par le récepteur. La précision de cette mesure dépend de la largeur de bande de la forme d’onde utilisée pour transmettre l’information de temps. Elle est fonction du produit BT, avec B la largeur de bande et T la durée d’un symbole, et du SNR (sigle anglais pour Signal to Noise Ratio, ou rapport signal à bruit). La valeur théorique maximale de sa variance, donnée par la borne de Cramer-Rao, est : avec N 0 la densité spectrale de bruit, et la variance de la mesure de temps d’arrivée. Les paramètres de la forme d’onde utilisée pour la transmission radio du message horaire et du top horloge doivent donc être dimensionnés de manière à offrir la précision recherchée. Par exemple, pour une largeur de bande B de 1 MHz avec une durée symbole T de 1 ms et un SNR de 10 dB, la variance de la mesure de temps sera de 1,2 ns. En règle générale, on vise des variances inférieures à 30 ns, voire inférieures à 10 ns. Avantageusement, la forme d’onde choisie peut être une forme d’onde par étalement de spectre, pour lesquelles la largeur de bande est généralement importante et permet d’obtenir des précisions de l’ordre de la nanoseconde pour des SNR très faibles. Par exemple, pour B égal à 20 MHz, T égal à 1 ms et un SNR de -30 dB, est égal à 1,3 ns. La représente schématiquement les étapes d’un procédé de transfert de temps selon un mode de réalisation de l’invention. Il comprend les étapes suivantes : une étape 201 de mesures de pseudo-distances avec au moins trois balises terrestres de transfert de temps, une étape 202 de calcul d’une position initiale du récepteur, à partir des pseudo-distances. Cette étape peut avantageusement mettre en œuvre une résolution initiale efficace reposant sur une forme clause, c’est-à-dire une forme donnant directement la solution recherchée en l’absence de bruit, une étape 203 d’estimation conjointe d’un biais d’horloge du récepteur, par la mise en œuvre d’un algorithme de résolution optimisant un critère de maximum de vraisemblance résolvant conjointement la position des balises terrestres de transfert de temps, la position du récepteur et le décalage temporel du récepteur par rapport aux balises de transfert de temps. L’étape 201 est une étape habituelle pour l’homme du métier, qui consiste à détecter les tops de synchronisation dans les signaux transmis par les balises pour déterminer un instant d’arrivée du signal, et à les confronter avec leurs instants d’émission, transmis dans le message horaire contenu dans le signal, afin de calculer une différence de temps, et par conséquent une pseudo-distance entre le récepteur et les balises de transfert de temps avec lesquelles il est en lien radio. La pseudo-distance Pd i mesurée entre la balise de transfert de temps i et le récepteur vaut Pd i = c(TOA i – Te i + Δt + e i ) = d i + c. Δt + c.e i , avec : TOA i l’heure d’arrivée du top de synchronisation dans le signal provenant de la balise i, mesuré par le récepteur à partir de son horloge locale, Te i l’heure d’émission par la balise du top de synchronisation, heure mesurée selon la référence temps du système de balises et transmise dans le message horaire, Δt le biais de l’horloge du récepteur, soit la différence entre le temps du système de balises et le temps de l’horloge locale du récepteur, e i un bruit blanc gaussien affectant la mesure, c la vitesse de la lumière (vitesse de propagation du signal radio), et d i la distance entre la balise i et le récepteur, Pd le vecteur formé à partir de l’ensemble des valeurs des Pd i . L’étape 202 consiste, à partir de N pseudo-distances Pd i obtenu à la première étape 201, à calculer une position initiale du récepteur et une estimée initiale du biais d’horloge de cet équipement, N étant supérieur ou égal à 3. Par exemple, cette position et ce biais d’horloge peuvent être obtenus par trilatération, en appliquant des techniques bien connues de l’homme du métier. A l’issue de l’étape 202, on dispose d’une position initiale du récepteur et du biais d’horloge. Ces valeurs sont peu précises puisque les mesures sont entachées des erreurs de positionnement des balises de transfert de temps (potentiellement importantes) et des bruits blancs gaussiens e i . L’étape 203 consiste à estimer le biais d’horloge Δt par un algorithme optimisant un critère de maximum de vraisemblance résolvant conjointement les positions des balises terrestres de transfert de temps, le biais d’horloge Δt et la position du récepteur. Cette étape peut être mise en œuvre par exemple par un algorithme de Gauss-Newton, bien connu de l’homme du métier, ou par un algorithme de résolution des moindres carrés équivalent, par exemple un algorithme de Levenberg-Marquardt. La suite de la description sera donnée en utilisant un algorithme de Gauss-Newton, de manière illustrative seulement. Les notations suivantes seront adoptées : désigne la position exacte de la balise de transfert de temps i, avec i [1, N], désigne la position connue de la balise de transfert de temps i, avec avec son erreur de positionnement, et est le vecteur formé à partir de l’ensemble des valeurs des , désigne la position estimée de la balise de transfert de temps i, initialisée à , et est le vecteur formé à partir de l’ensemble des valeurs des , désigne l’estimation du biais d’horloge du récepteur. Toutes les notations de position peuvent être décomposées selon un premier axe x et un deuxième axe y orthogonal à x dans un plan horizontal, par exemple et . Selon un mode de réalisation, l’étape 203 comprend la construction d’un vecteur d’état comprenant la position estimée du récepteur, initialisée à l’étape 202, le biais d’horloge estimé et les positions estimées des balises de transfert de temps, avec par exemple : L’étape 203 comprend également la construction d’un vecteur d’observation à partir des positions connues ou transmises des balises et des pseudodistances calculées lors de la première étape 201, et de son pendant linéarisé , qui comprend l’écart entre les positions connues des balises de transfert de temps et leurs positions estimées, et l’écart entre les pseudo-distances calculées lors de la première étape 201 et des pseudo-distances construites à partir des positions estimées du récepteur et des balises terrestres de transfert de temps, et du décalage temporel entre le récepteur et les balises : avec la distance entre la position estimée du récepteur et la position estimée de la balise terrestre de transfert de temps d’indice , et le vecteur construit à partir des . Lors de la première itération, les positions estimées des balises terrestres de transfert de temps sont égales à leurs positions connues ou transmises , et le décalage temporel estimé est nul. L’étape 203 comprend également la construction de la matrice Jacobienne à partir des vecteurs et , avec : Par construction, et . Par conséquent : Par construction toujours : de sorte que : La solution de ce système est calculée par un estimateur du maximum de vraisemblance mis en œuvre par exemple avec un algorithme de Gauss-Newton initialisé avec les valeurs calculées lors de la deuxième étape 202 du procédé selon l’invention, par la formule : avec la matrice de covariance d’observation, qui permet de caractériser la précision de l’estimateur. L’estimée du vecteur d’état , qui contient en particulier une estimation du biais d’horloge du récepteur peut donc être obtenue très simplement à partir de la matrice et du vecteur . Avantageusement, ce calcul est réalisé itérativement, en mettant à jour à chaque itération le vecteur d’observation et la matrice avec les valeurs du vecteur d’état calculé à l’itération précédente, afin de linéariser et améliorer l’estimation de à chaque itération. Les itérations peuvent être stoppées par exemple en évaluant un critère d’arrêt sur la norme , avec le vecteur d’état obtenu à l’itération , et le vecteur d’état obtenu à l’itération . Avantageusement, des mesures successives du vecteur peuvent adresser un filtre de lissage récursif afin notamment de prendre en compte la mobilité du récepteur et/ou des balises de transfert de temps. Le calcul décrit ci-dessus réalise l’optimisation conjointe de la position du récepteur, de son biais d’horloge, et des positions des balises terrestres de transfert de temps. Cette optimisation conjointe permet d’améliorer la qualité de la mesure du biais d’horloge lorsque la position des balises de transfert de temps est imprécise. La matrice de covariance d’observation vaut : où : est la matrice de covariance d’erreur sur la connaissance des positions des balises. Les mesures de position des balises sont généralement supposées indépendantes entre elles. Cette matrice est donc constituée d’une diagonale de matrices carrées des variances de position de chaque balise i. Le bruit de mesure est supposé connu, gaussien et centré, de sorte que , avec E() l’espérance mathématique et le vecteur bruit gaussien sur la mesure de position Mi de la balise i ( , est la matrice d’observation des pseudo-distances Pdi. Le bruit de mesure est supposé connu, gaussien et centré, et , avec la valeur de à la ligne i et la colonne i. La matrice de covariance de l’estimateur vaut alors : Les équations données ci-dessus sont formées à partir de matrices agencées selon une organisation donnée à titre d’illustration seulement. Les mêmes résultats peuvent être atteints en agençant les matrices de manière différente, par exemple en permutant leurs lignes, sans que cela n’affecte la précision des calculs et les résultats obtenus. Selon un mode de réalisation avantageux, la deuxième étape 202 du procédé, qui vise à déterminer une position initiale et un décalage temporel initial du récepteur, est mise en œuvre par une estimation pseudo-linéaire reposant sur une forme clause. Cette solution consiste à parcourir les triplets (i, j, k) de balises de transfert de temps, avec . Pour chaque triplet d’indice q allant de 1 à Q, on calcule alors une pseudo-mesure f q telle que : avec O le centre du repère, et la différence entre la pseudo-distance Pd i et la pseudo-distance Pd j , soit : On calcule également le vecteur à partir des mesures de pseudo-distances et de la position connue des balises, soit : Enfin, on construit le vecteur d’observation à partir des mesures de pseudo distances, soit : Une estimation pseudo-linéaire de la position initiale du récepteur et de son biais d’horloge Δt est obtenue en minimisant les résidus du système linéaire , où est un vecteur constitué de l’ensemble des pseudo mesures , et une matrice construite par concaténation des vecteurs . Pour cela, il est possible de poser : Les valeurs de la position du récepteur et de son biais d’horloge sont ensuite utilisées pour initialiser le vecteur d’état lors de la troisième étape 203 du procédé de transfert de temps selon l’invention. Ce mode de réalisation est particulièrement intéressant, puisqu’il permet d’obtenir une première approximation très précise de la position et du biais d’horloge du récepteur, ce qui garantit le bon fonctionnement de l’algorithme optimisant un critère de maximum de vraisemblance mis en œuvre lors de la troisième étape 203 du procédé selon l’invention. En effet, compte tenu de la forte non linéarité du problème d’estimation, l’algorithme de Gauss-Newton présente une forte probabilité de divergence en l’absence d’une initialise précise. La représente une configuration opérationnelle dans laquelle peut être mis en œuvre un procédé de transfert de temps selon l’invention, afin de présenter les résultats qu’il permet d’atteindre. Le récepteur 301 est en lien radio avec quatre balises terrestres de transfert de temps 302 à 305, équiréparties sur un cercle centré autour de l’équipement 301. Il s’agit là d’un cas favorable. La position des balises est connue avec un CEP de 50% (acronyme anglais pour Circle Error Probable, ou erreur circulaire probable) de 10 mètres. L’équipement 301 est capable de mesurer les temps d’arrivée des messages horaires avec une variance de 2 ns. La mise en œuvre du procédé de transfert de temps selon l’invention permet de synchroniser l’équipement 301 avec une précision de ±40 ns, pour un intervalle de confiance à 99%. L’invention porte donc sur un procédé de transfert de temps d’un terminal en réception radiofréquence avec au moins trois balises terrestres de transfert de temps synchronisées diffusant chacune un message horaire comprenant un marquant horaire et une heure d’émission, mais également sur le terminal lui-même. Un tel terminal doit disposer d’une chaîne radio lui permettant de recevoir et d’interpréter le message horaire diffusé par les balises pour déterminer un instant d’arrivée de marquants horaires, et de moyens de calculs lui permettant de mettre en œuvre le procédé, comme par exemple un microprocesseur, un DSP (sigle anglais pour Digital Signal Processor , ou processeur de signal numérique), un FPGA (sigle anglais pour Field Programmable Gate Array , ou réseau de portes programmable), un ASIC (acronyme anglais pour Application- Specific Integrated Circuit , ou circuit intégré propre à une application), n’importe quelle association de ces moyens, ou n'importe quel composant matériel permettant d'exécuter les fonctions précitées. Enfin, il porte sur un programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé de transfert de temps selon l’invention lorsque ledit programme est exécuté sur un ordinateur, ainsi que sur un support d’enregistrement lisible par un ordinateur sur lequel est enregistré un tel programme d’ordinateur. Procédé de transfert de temps vers un récepteur (101) en lien radiofréquence avec au moins trois balises terrestres de transfert de temps (102, 103, 104) synchronisées diffusant chacune un message horaire comprenant un marquant horaire et une heure d’émission, le procédé de transfert de temps étant caractérisé en ce qu’il comprend : - une première étape (201) de calcul d’un instant d’arrivée des messages horaires transmis par lesdites balises terrestres de transfert de temps, les positions desdites balises terrestres de transfert de temps étant connues par le récepteur ou transmises au récepteur, et de détermination d’une mesure de pseudo-distance Pd i associée à chaque balise de transfert de temps, - une deuxième étape (202) de détermination d’une position initiale et d’un décalage temporel initial dudit récepteur à partir des mesures de pseudo-distances calculées lors de la première étape (201), - une troisième étape (203) de calcul d’un décalage temporel fin dudit récepteur, par la mise en œuvre d’un algorithme de résolution optimisant un critère de maximum de vraisemblance utilisant les valeurs calculées à la deuxième étape (202) et résolvant conjointement les positions des balises terrestres de transfert de temps, la position du récepteur et le décalage temporel du récepteur. Procédé de transfert de temps selon la revendication 1, dans lequel ladite troisième étape (203) comprend : - la construction d’un vecteur d’état comprenant la position estimée du récepteur, le décalage temporel estimé du récepteur, et les positions estimées desdites au moins trois balises terrestres de transfert de temps, - la construction d’un vecteur d’observation comprenant les positions connues ou transmises desdites au moins trois balises terrestres de transfert de temps et les pseudo-distances calculées lors de la première étape (201), la construction d’un vecteur d’observation linéarisé comprenant un écart entre les positions connues ou transmises desdites au moins trois balises terrestres de transfert de temps et leurs positions estimées, et un écart entre les pseudo-distances calculées lors de la première étape (201) et des pseudo-distances construites à partir des valeurs du vecteur d’état , -la construction d’une matrice Jacobienne à partir du vecteur d’état et du vecteur d’observation tel que , et -la résolution du système, en calculant . Procédé de transfert de temps selon la revendication 2, dans lequel le calcul de est itéré plusieurs fois, les valeurs du vecteur et de la matrice étant recalculées à chaque itération en fonction du vecteur déterminé à l’itération précédente. Procédé de transfert de temps selon l’une des revendications précédentes, dans lequel l’algorithme utilisant un critère de maximum de vraisemblance mis en œuvre lors de la troisième étape (203) est un algorithme de Gauss-Newton ou un algorithme de Levenberg-Marquardt. Procédé de transfert de temps selon l’une des revendications précédentes, dans lequel la deuxième étape (202) est mise en œuvre par un procédé d’estimation pseudo-linéaire de la position initiale et du décalage temporel initial dudit récepteur. Procédé de transfert de temps selon la revendication 5, dans lequel le procédé d’estimation pseudo-linéaire comprend : - le calcul d’une position initiale du récepteur, avec , où est une matrice construite à partir des mesures de pseudo-distances Pd i calculées lors de la première étape 201 du procédé et des positions connues ou transmises des balises terrestres de transfert de temps, est un vecteur d’observation comprenant les mesures de pseudo-distances Pd i faites lors de la première étape (201) du procédé, et est une pseudo-mesure construite à partir des mesures de pseudo-distance Pd i faites lors de la première étape 201 du procédé et des positions connues ou transmises des balises, - le calcul d’un décalage temporel initial tel que , avec N le nombre de balises terrestres de transfert de temps et un vecteur comprenant les positions connues ou transmises des balises terrestres de transfert de temps. Procédé de transfert de temps selon l’une des revendications précédentes, dans lequel la position des balises terrestres de transfert de temps est connue de manière imparfaite. Equipement de radiocommunications (101) en lien radiofréquence avec au moins trois balises terrestres de transfert de temps (102, 103, 104) synchronisées diffusant chacune un message horaire comprenant un marquant horaire et une heure d’émission, l’équipement de radiocommunications comprenant une chaîne radio configurée pour recevoir et interpréter le message horaire transmis par les balises, l’équipement de radiocommunications étant caractérisé en ce qu’il comprend en outre des moyens de calculs configurés pour mettre en œuvre un procédé de transfert de temps selon l’une des revendications précédentes. Programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé de transfert de temps selon l’une des revendications 1 à 7 lorsque ledit programme est exécuté sur un ordinateur. Support d’enregistrement lisible par un ordinateur sur lequel est enregistré un programme d’ordinateur selon la revendication 9.