FR 3085568 A1 20200306 FR 1857873 A 20180831 PROCEDE DE DATATION DE SIGNAUX DE TELEMESURE DOMAINE DE L'INVENTION L'invention concerne un procédé de datation de réception de s données numériques d'un signal modulé. Par ailleurs, l'invention concerne un récepteur de télécommunication recevant un signal radio modulé. ETAT DE LA TECHNIQUE to Dans le contexte d'aéronefs ou d'engins spatiaux, il est nécessaire de pouvoir, depuis le sol, localiser précisément dans l'espace et le temps ces derniers. Typiquement, dans la technique de télémétrie (dite de « ranging » en anglais), il est procédé à la mesure des temps de trajet d'un signal émis par l'engin ou depuis le sol et reçu dans plusieurs stations 15 distantes au sol. De nombreuses techniques existent : télémétrie par impulsions, télémétrie semi-active, télémétrie passive par corrélation. L'exigence de précision nécessite souvent d'émettre un signal ad-hoc, de largeur de bande suffisante, sans gigue (fluctuation du signal) et incluant des zo marqueurs caractéristiques (télémétrie active ou semi active). Cette solution présente un coût intrinsèque en termes de ressources comme le dispositif d'émission spécifique, de la consommation spectrale ... Dans le cas d'une technique de télémétrie passive, il est utilisé de façon opportuniste un signal télécom, mais ce signal n'est pas toujours 25 présent (satellite d'observation ou LEO) et nécessite de recevoir et traiter séparément des volumes de données importants. Egalement, il est possible d'utiliser un signal de télémesure TM toujours présent sur les engins inhabités 10 (satellites, drones ...) qui est reçu dans les stations de contrôle, et de le dater au niveau de chacune de 30 ces stations, comme illustré en figure 1. Ainsi, les différentes dates de réception du signal TM, TOA1 au niveau de la station 20A, TOA2 au niveau de la station 20B, permettent de calculer les écarts de temps de trajet 2 entre un engin et les différentes stations, et de localiser ensuite ledit engin par trilatération hyperbolique par exemple. La figure 2 illustre la structure typique d'un récepteur de signal de télémesure Stm d'une station de contrôle, ledit signal Stm étant obtenu s après échantillonnage et numérisation asynchrone du signal TM. Ledit récepteur comporte classiquement un bloc 40 démodulateur/démappeur permettant d'extraire les symboles du train d'échantillons, un bloc d'estimation/correction 50 d'écart de phase et de fréquence entre un oscillateur local et la porteuse, un bloc 60 d'estimation 10 du rythme et phase symbole. Ce démodulateur est précédé par un ré-échantillonneur 30 capable de synthétiser des échantillons représentatifs du signal Stm entre ceux échantillonnés à l'entrée dudit récepteur. Ainsi, le ré-échantillonneur 30 peut ajuster la fréquence d'échantillonnage et la phase d'interpolation efficaces servant à re-cadencer de manière 15 synchrone avec les symboles les échantillons intermédiaires. En effet, il est important d'échantillonner le signal de télémesure exactement à l'instant de l'ouverture maximale du diagramme de l'oeil pour obtenir le symbole avec le rapport signal sur bruit le plus élevé. Un décalage de cet instant réduira de manière significative le rapport signal à 20 bruit (SNR) effectif vu par les blocs récepteurs suivants ainsi que l'interférence inter-symbole (ISI). Cette structure classique permet moyennant un asservissement (par exemple du type boucle de Costas) du récepteur sur une horloge de référence externe (type GPS) d'obtenir des précisions de datation du signal 25 reçu de l'ordre de la microseconde. Cependant, avec un tel récepteur, il n'est pas possible d'améliorer la précision de la datation. En effet, comme vu précédemment, le ré-échantillonneur procède au recalage à un instant optimal vis-à-vis du rythme symbole, ce qui introduit une gigue temporelle liée à la 30 caractéristique aléatoire de la séquence de symbole et au bruit reçu dans la bande passante. Un tel traitement génère ainsi une gigue importante 3 proportionnelle à la période d'échantillonnage d'entrée, de l'ordre de la microseconde dans le cas de signaux de télémesure satellite. Dans le contexte de la localisation d'aéronefs ou d'engins spatiaux, un tel ordre de précision de la datation se révèle insuffisant pour permettre s une précision de localisation appropriée (erreur de l'ordre d'un ou plusieurs kilomètres). Il existe donc un besoin d'améliorer les récepteurs de signaux de télémesure pour pouvoir proposer une datation plus précise, de la réception de tels signaux. 10 PRESENTATION DE L'INVENTION L'invention a pour but de pallier les inconvénients de l'art antérieur décrits ci-avant. En particulier, un but de l'invention est de proposer un procédé de 15 datation de données numériques d'un signal modulant de manière cohérente la porteuse ou la sous-porteuse, et de tirer avantage des propriétés d'une telle modulation pour en améliorer la datation. Un tel procédé permet avantageusement de corriger avec une grande précision la datation d'une séquence pilote (pattern connu), et zo réduit considérablement la variance sur la datation des transitions des bits pour permettre entre autres une meilleure localisation des engins mobiles émettant un tel signal de télémesure Un autre but de l'invention est de permettre au procédé d'exploiter un signal de télémesure existant et ainsi d'éviter l'émission d'autres 25 signaux pour faire la mesure de distance. Ce procédé permet également d'augmenter la précision de la resynchronisation à bord de satellites ne disposant pas d'une référence temporelle d'une grande précision, lorsque celle-ci est faite à partir de la date de réception d'une séquence donnée reçue au sol, et dont la date de 30 réception est renvoyée au bord. 4 A cet égard, l'invention a pour objet un procédé de datation de réception de données numériques d'un signal modulé, ledit signal résultant de la modulation d'une porteuse ou sous-porteuse par un signal numérique, le débit de symboles du signal numérique étant un sous- multiple entier N de la fréquence de la porteuse ou sous-porteuse, ledit procédé comprenant les étapes mises en oeuvre par un processeur d'un récepteur de télécommunication, consistant à, après réception et échantillonnage du signal modulé : démodulation cohérente au moyen d'une boucle de phase d'une io pluralité d'échantillons sm(k) du signal modulé et obtention d'une pluralité d'échantillons démodulés sdm(k) et de la phase de la porteuse reconstituée ; association de la pluralité échantillons démodulés sdm(k) avec la phase simultanée ça(k) de la porteuse reconstituée 15 détermination d'une pluralité de dates de passage à une valeur déterminée p0 de la phase de la porteuse reconstituée, relativement à la date de réception d'au moins un échantillon sm(k), selon un décalage de temps pendant lequel la porteuse tourne de p(k) à (p0. 20 Avantageusement, mais facultativement, le procédé selon l'invention peut en outre comprendre au moins l'une des caractéristiques suivantes : l'étape de démodulation cohérente comprend les sous-étapes de : 25 a. détection d'une erreur de phase par un détecteur de phase qui efface la modulation du signal reçu; b. filtrage passe-bas de l'erreur de phase ; dans l'étape de détermination des dates de passage à (p0, cette date est calculée relativement aux échantillons sm(k) du signal modulé 30 qui précédent ou suivent immédiatement un passage à rp0 de la phase de la porteuse reconstituée ; le procédé comporte en outre les étapes de : 5 o sélection d'une date de transition symbole parmi la pluralité de dates de passage à p0, ladite date de transition symbole correspondant à un instant de transition de symboles du signal numérique ; 5 o décimation d'un facteur N des dates de passage à p0, à partir de ladite date de transition symbole détectée ; et o détermination des dates de réception d'une pluralité de symboles consécutifs correspondant aux dates de passage à p0 résultantes de la décimation ; l'étape de sélection de la date de transition symbole, comprend une sous-étape d'intégration I(n) des échantillons démodulés sdm (k) sur une fenêtre de N échantillons consécutifs du signal ; - la sous-étape d'intégration I(n) est effectuée sur des durées d'un symbole démarrant à chaque échantillon sdm(k) du signal démodulé correspondant à un passage à p0 ; - la sous-étape d'intégration comporte également : c. le calcul pour une pluralité d'indices consécutifs p de la valeur moyenne M(p) de la valeur absolue de ladite intégration I(m), pour m prend Ns valeurs qui valent toutes p modulo N, et d. la détermination de la date de transition symbole en sélectionnant le passage à p0 d'indice p qui maximise M(p) ; - la valeur de l'intégration I(n) démarrant à la date de transition symbole fournit le symbole correspondant du signal numérique ; la comparaison des valeurs de M calculées aux passages à p0 précédents M(n-1) et suivants M(n+1) la transition de symbole détectée permet de caractériser une ambiguïté de phase liée à la modulation lorsque ces valeurs ne sont pas sensiblement égales mais que l'une est sensiblement égale à M(n) ; l'ambiguïté de phase est résolue en modifiant la phase de passage prédéterminée (p0 à une autre valeur ; le procédé comporte également le décodage des symboles en données binaires ; 6 le procédé comporte également une étape de synchronisation de trames qui comprend une seconde décimation d'un facteur Q égal à une taille prédéterminée d'une trame de données binaires, à partir d'une date de transition symbole, et une étape de détermination des s dates de réception d'une pluralité de trames consécutives correspondant aux dates de passage à (p0 résultantes de la décimation ; et l'étape de synchronisation de trames comporte en outre la sélection d'une trame ayant la meilleure corrélation avec une séquence pilote 10 déterminée, parmi une pluralité de trames consécutives et l'étape de détermination de la date de réception de ladite trame sélectionnée. Par ailleurs, l'invention concerne un récepteur de télécommunication 15 recevant un signal radio modulé, ledit signal résultant de la modulation d'une porteuse ou sous-porteuse par un signal numérique, le débit de symboles du signal numérique étant un sous-multiple entier N de la fréquence de la porteuse ou sous-porteuse, ledit récepteur comprenant : une unité de démodulation configurée pour, après réception et 20 échantillonnage du signal modulé, la : o démodulation cohérente au moyen d'une boucle de phase d'une pluralité d'échantillons sm(k) du signal modulé et obtention d'une pluralité d'échantillons démodulés sdm(k) et de la phase de la porteuse reconstituée ; 25 o association de la pluralité échantillons démodulées sdm(k) avec la phase simultanée 4)(0 de la porteuse reconstituée o détermination d'une pluralité de dates de passage à une valeur déterminée (p0 de la phase de la porteuse reconstruite, relativement à la date de réception d'au moins un échantillon 30 sm(k). 7 Avantageusement, mais facultativement, le récepteur de télécommunication selon l'invention peut en outre comprendre au moins l'une des caractéristiques suivantes : une unité (130) de synchronisation configurée pour la : 5 o sélection d'une date de transition symbole parmi la pluralité de dates de passage à (p0, ladite date de transition symbole correspondant à un instant de transition de symboles du signal numérique ; o décimation d'un facteur N des dates de passage à TO, à partir 10 de ladite date de transition symbole détectée ; et o détermination des dates de réception d'une pluralité de symboles consécutifs correspondant aux dates de passage à (p0 résultantes de la décimation. 15 En outre, l'invention concerne un système de localisation d'engin spatial ou aéronef mettant en oeuvre un procédé et/ou un récepteur de télécommunication selon l'une des caractéristiques précédemment décrites DESCRIPTION DES FIGURES 20 D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre, au regard des figures annexées, données à titre d'exemples non limitatifs et sur lesquelles : la figure 1 déjà présentée, représente un système de télémesure 25 satellite selon l'état de l'art ; la figure 2 déjà présentée, représente un récepteur radio selon l'état de l'art ; la figure 3A représente les principales étapes d'un procédé de datation mise en oeuvre par un récepteur radio selon l'invention ; 30 la figure 3B illustre schématiquement un récepteur radio 100 selon l'invention ; 8 la figure 4 représente schématiquement une unité 120 de démodulation d'un récepteur radio 100 selon l'invention ; la figure 5 représente schématiquement une unité 130 de synchronisation d'un récepteur radio 100 selon l'invention ; 5 la figure 6 représente schématiquement une étape d'intégration mise en oeuvre par l'unité 120 de démodulation selon l'invention ; et la figure 7 représente schématiquement une unité 140 de détection d'un récepteur radio 100 selon l'invention. 10 DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE REALISATION DE L'INVENTION 15 La télémesure satellite, conformément aux normes en vigueur (définies selon le Comité Consultatif pour les Systèmes de Données Spatiales : CCSDS), transmet des données numériques, le plus souvent binaires ou « bits » qui peuvent moduler une porteuse selon plusieurs méthodes, comme la modulation de phase (PM), le changement de phase 20 (PSK), et/ou l'impulsion et codage (PCM). Dans tous les cas, elle présente une modulation des bits cohérente de la porteuse PSK, ce qui revient à dire que le rythme symbole Rs est un sous-multiple entier N de la fréquence porteuse (ou sous-porteuse) PSK, Fp. Ainsi, dans ce contexte Fp = N x Rs. Un signal de télémétrie typique comporte alors un débit binaire 25 compris entre 1400 et 31200 bits/seconde. Le ratio entre la fréquence Fp de la porteuse et le débit binaire est typiquement compris entre 4 et 16. Dans la suite de la description, le terme « porteuse » peut être également remplacé par le terme « sous-porteuse ». En effet, vu du signal modulant il n'y a pas lieu de différencier ces termes, la porteuse PSK 30 pouvant également moduler un autre signal, ainsi le procédé décrit ci-après s'applique indifféremment à une porteuse ou sous-porteuse. 9 La figure 3A illustre un procédé de datation de signaux mis en oeuvre par un récepteur 100 de communication radio recevant en entrée un signal radio, tel qu'un signal de télémesure, provenant d'un satellite par exemple. 5 En référence à la figure 3B, il est illustré un tel récepteur 100 radio. Ce dernier comprend : - une unité 120 de démodulation ; - une unité 130 de synchronisation ; et - une unité 140 de détection de séquence pilote. 10 Unité de démodulation Le récepteur 100 reçoit un signal Stm, correspondant à un signal de télémesure TM préalablement échantillonné à un rythme d'échantillonnage déterminé Fe. 15 La fréquence d'échantillonnage Fe peut être verrouillée sur la base locale de temps ce qui signifie que la période Te= 1/Fe est parfaitement connue sans erreur ni dérive, les dates d'échantillonnages sont donc exactement à k*Te avec k entier. De plus, il peut être admis sans limitation que l'échantillonnage 20 respecte la condition de Nyquist et donc en particulier que Fe > 2 Fp ce qui implique qu'on a au moins deux échantillons par tour de phase de la sous-porteuse ou porteuse. Une unité de démodulation 120 du récepteur, illustrée en figure 4, dans une étape E11, réceptionne ainsi une pluralité d'échantillons Sm(k), 25 où k est un indice compris entre 0 et Ne-1 où Ne désigne le nombre d'échantillons de la portion du signal radio, et constitue des blocs d'échantillons. L'unité 120 date un échantillon k (par exemple le premier) d'un bloc d'échantillons avec une date donnée Tref. Cette date Tref peut être définie par rapport à une base locale de temps du récepteur ou 30 asservie à une référence externe au récepteur, par exemple un décalage par rapport à un pulse GPS (PPS), un compteur d'une horloge synchronisée sur une référence externe (exemple 10MHz GPS), etc. 10 L'unité de démodulation 120 du récepteur 100 procède, dans une étape E21, à la démodulation du signal réceptionné Sm(k) en Sdm(k). L'unité de démodulation 120 peut comprendre un bloc de filtrage adapté 5 121 qui maximise le rapport signal sur bruit dudit signal. L'unité de démodulation 120 comprend également un bloc décimateur 122 de l'unité de démodulation 120 qui permet de réduire réduit le nombre d'échantillons à traiter par la suite. L'unité 120 comprend également une boucle à verrouillage de phase, 10 mise en oeuvre sur le signal reçu Sm(k) par un bloc détecteur de phase 123, et un bloc filtre de boucle 124, ladite boucle étant asservie par un signal de référence généré par un bloc oscillateur 125, de type NCO (oscillateur contrôlé numériquement) par exemple. Le bloc détecteur de phase 123 génère un signal d'erreur qui 15 commande la boucle à verrouillage de phase. Ainsi placé après le bloc décimateur 122, le bloc détecteur de phase 123 mesure la différence entre la phase du signal généré par un oscillateur local 125 et la phase de la porteuse du signal reçu Sm(k). Dans le cas où la modulation du signal modulé Sm(k) possède une 20 ambiguïté de phase d'ordre Netat, alors le détecteur de phase 123 doit fournir une erreur nulle pour les Netat valeurs de phase ( 21rp Netatpour p entier) que peut prendre la modulation. La fonction du détecteur de phase 123 présente alors une périodicité d'ordre Netat ce qui est généralement obtenu en multipliant la phase par Netat. 25 Le bloc filtre de boucle 124 (souvent d'ordre 2 pour annuler le biais constant de phase) filtre le signal d'erreur de phase afin de fournir un meilleur signal au bloc oscillateur 125. En référence de nouveau à la figure 3A, dans une étape E22, l'unité 120 transmet à une unité de synchronisation 130, la date Tref, et une 30 pluralité de couples associant un échantillon démodulé Sdm(k) avec la phase simultanée (p(k) de la porteuse reconstruite par l'oscillateur 125. 11 Avantageusement, cette fourniture de la phase reconstruite associée à chaque échantillon permet d'affiner la datation. Dans un mode de réalisation particulier, à l'étape E22, une zone mémoire 126, telle qu'une zone tampon, stocke la date Tref associée à la 5 pluralité des couples Sdm(k), go(k)ainsi décrits. Avantageusement, la connaissance de Tref permet de connaître la date précise de chaque couple Sdm(k), cp(k) au sein de cette zone de stockage puisque que les échantillons sont espacés d'une période Te qui est parfaitement connue. 10 Ensuite, ladite zone mémoire 126 est transmise à une unité de synchronisation 130. A ce stade, le récepteur 100 dispose donc d'échantillons datés dont la phase porteuse est parfaitement renseignée en regard du signal transmis par l'émetteur. 15 Unité de synchronisation (datation) Une unité de synchronisation 130 reçoit en entrée la zone mémoire 126 sortant de l'unité 120. Ainsi, comme illustré en figure 5, un bloc d'intégration 131 reçoit sur son 20 entrée les échantillons démodulés Sdm(k) accompagnés de leur phase porteuse reconstruite co(k) et de la datation Tref, par exemple via la zone mémoire 126. Par la suite, il est considéré que Tref est l'instant défini pour k=0, k référençant l'indice de l'échantillon dans la zone mémoire En référence à la figure 6, il est illustré un symbole, échantillonné 25 en une pluralité d'échantillons k. Il y a plusieurs périodes porteuses Scf par symbole et plusieurs échantillons par période porteuse Scf. Ainsi, Scf(n) représente la n-ième période de la porteuse. La modulation des symboles étant cohérente de la porteuse PSK, le rythme symbole Rs est un sous-multiple entier de la fréquence porteuse PSK notée Fp = N x Rs, et donc 30 un symbole module la porteuse sur N périodes de celle-ci. Ainsi, chaque symbole commence sur un des N passages de la phase à une même valeur déterminée (,00. Cette valeur déterminée (,00 est, en 12 principe (mais pas nécessairement) proche de zéro, et dans le reste de la description sera appelée « passage à zéro / 0 » de la phase w(k) de la porteuse reconstruite. Dans une étape E30, la pluralité d'échantillons démodulés Sdm(k) 5 accompagnés de leur phase porteuse reconstruite w(k) est traitée pour déterminer un ensemble d'échantillons k reçus en fonction de la période de la porteuse. Ainsi, dans une sous-étape E31 de E30, Le bloc de synchronisation 131 détecte un passage à zéro du signal de porteuse reconstitué en 10 détectant le passage à 0 (4)0 modulo 27r) de la phase ço(k). Un instant T(n) est celui du nième passage de la porteuse à la phase 0, la période n de la porteuse sur laquelle portera une intégration du signal sdm(k) est définie entre les instants T(n) et T(n+1). A chaque passage à 0 (rp0 modulo 2n), le bloc 131 d'intégration 15 procède à l'intégration du signal Sdm(k) sur une période n de la porteuse. Au passage à 0 suivant, débutant la période n+1, le bloc 131 d'intégration réinitialise le calcul d'intégration qui est remis à 0 et le résultat S(n) de l'intégration précédente est délivré simultanément en sortie à une pluralité de blocs 132 de sommation qui opèrent en parallèle. 20 Dans une sous-étape E31 de E30, chaque sortie S(n) est par ailleurs associée à une datation obtenue par extrapolation de la date du kième échantillon, par exemple, selon la formule suivante valable pour des phases calculées en radians: T(n) = Tk +AçOk/27CFp = Tref + k(n)*Te - ko(k(n)) - ço0 I / 27cF p , k(n) 25 correspondant à l'indice de l'échantillon k désigné comme succédant au passage à 0 de la phase reconstruite. Le terme correctif de T(n) qui vaut frp (k) - 4)0 1127cFp calcule le délai exact qui sépare le passage à zéro de la phase de la porteuse avec l'échantillon d'indice k et T(n) date donc précisément le passage à zéro qui 30 correspond possiblement à l'instant de début du symbole. De façon alternative, k pourrait désigner l'indice de l'échantillon précédant le passage à zéro de la phase reconstruite et la datation de T(n) 13 serait alors : T(n) = Tk +(27r- 4400/27(Fp = Tref + k(n)*Te + [27c + p0 - p(k(n))1/27cFp. On remarque ici que le passage à zéro de la phase correspond à cp(k(n)) E Fru-Fp°, TO] et v(k(n) + [(p0,+7c+(p0[ qui a forcément une s solution à chaque tour puisqu'il y a au moins deux échantillons par tour. Les valeurs de phase prises dans cet intervalle sont donc exemptes d'ambiguïté à 2n près. Également, de façon alternative, la formule suivante d'interpolation entre les échantillons peut être utilisée : 10 T(n) = Tk +A(pk/(cpk+i -pk)Te = Tref + [k(n) +(po - pk)/(pk+i - pk) ]*Te. Cette dernière formule est un peu plus précise que les précédentes si on y utilise la valeur nominale de Fp, en revanche elles sont équivalentes si Fp prend la valeur réelle de commande instantanée du NCO puisque 15 ( = 2nFp*Te.k+t - Pk) Un dispositif de démodulation classique n'aurait pas utilisé la valeur de la phase reconstituée et n'aurait donc pas daté finement le passage à 0 mais aurait synchronisé le rythme par une boucle de récupération de zo rythme en rajoutant une gigue (bruit de phase) inhérente à ce type de dispositif. Le dispositif proposé procède à l'intégration du signal pour détecter les symboles, mais à la différence d'un démodulateur classique le fait avant même d'avoir déterminé les instants de transition entre symboles. 25 En effet, dans un mode de réalisation de la sous-étape E32 de E30, N blocs 132 de sommation réalisent en parallèle l'intégration sur une fenêtre de N échantillons (soit un symbole) S(n) sortant du bloc 131 par la fonction d'intégration symbole suivante [p(m) = E',llig±±Pp++1*Z+N-15(j) avec o hz par Netat ce qui crée l'ambiguïté modulo 20 2n/Netat. Ainsi, dans le cas d'une modulation à deux états de phase, si la valeur M(p) définie plus haut a été calculée à partir de la bonne demi-période sur une période porteuse, alors M(p-1) est environ égal à M(p+1) au bruit près et M(p) est strictement supérieur à ces deux valeurs. 25 Dans le cas contraire, si les intégrations se sont produites sur la mauvaise demi-période alors M(p) est sensiblement égal à M(p-1) ou M(p+1) car le pic est alors partagé sur la position d'une demi-période avant et d'une demi-période après la bonne position de phase de la porteuse. Pour lever l'ambiguïté de demi-période, dans une réalisation 30 préférée, le module 133 calcule un seuil SO = (M(p modulo N) + min(M((p- 1) modulo N ),M((p+1) modulo N)))/2,. Ensuite, le module 133 compte le 16 nombre de valeurs au-dessus de ce seuil parmi M(p modulo N),M((p-1) modulo N), M((p+1) modulo N). Dans le cas où il n'y a qu'une valeur au-dessus du seuil SO, alors l'intégrale du symbole est calculée sur la bonne demi-période et est s d'amplitude maximale. Dans le cas contraire, les intégrations se sont produites sur la mauvaise demi-période, il faut alors décaler l'intégration d'une demi-période. Le bloc 133 détermine ainsi la demi-période correcte et fournit cette 10 information au bloc 131 qui détecte en conséquence le passage à 0 ou à TC (en fait ço0+7c) de la phase co(k) en fonction du sélecteur de demi-période de démarrage. Sur sa sortie, le bloc 133 fournit à chaque symbole d'indice i, à un rythme correspondant au rythme de la porteuse décimé par N, un couple 15 de données (I(i), Tm(i)). I(i) étant l'intégrale d'amplitude maximale I(popt) sélectionnée à l'étape précédente et Tm(i) étant T(popt). De façon alternative, T(i) peut être une date moyennée sur un intervalle de durée fixe autour de popt et décalée d'un délai fixe pour indiquer le milieu ou un autre endroit prédéfini dans le symbole selon la convention de datation zo adoptée. Ce couple de données est fourni à une unité de détection 140. Les informations de décalage de cycle et demi-cycle de la porteuse vis-à-vis des transitions de symboles sont déterminées par le bloc 133 sélecteur et retournées au bloc 131 dans un objectif d'initialisation des 25 calculs de couple de données par l'unité de synchronisation 130. Lesdites informations peuvent être également déterminées à nouveau et à intervalle régulier dans un objectif de vérification des calculs de couple de données par l'unité de synchronisation 130. Dans le cas où l'ambiguïté de phase est d'un ordre supérieur à deux, 30 il convient de modifier quelque peu l'algorithme présenté ci-dessus. La méthode générale consiste à calculer le critère de synchronisation symbole utilisé dans la boucle de récupération de rythme standard pour la 17 modulation considérée, ceci aux N positions possibles de passage à zéro, puis à retenir celle qui maximise ce critère. Par exemple, dans la boucle de Costas pour 4 états, on effectuera une conversion de fréquence en nombre complexe dans l'unité de démodulation 120 et le critère devient le s maximum de M(n) calculé à la fois sur la voie I en phase et la voie Q en quadrature. A ce stade, le récepteur 100 dispose donc d'une valeur intégrale I(p) correspondant à la valeur d'un symbole brut (avant décodage et décision) à déterminer, ladite valeur étant associée à une datation de l'instant de 10 passage à zéro de la phase de la porteuse qui débute le symbole. Avantageusement, cette datation est plus précise que celle que l'on aurait obtenue par un algorithme de synchronisation de rythme classique, car celui-ci introduit une gigue liée à l'asynchronisme intrinsèque entre le rythme d'échantillonnage du signal et le rythme symbole. 15 Unité de détection 140 L'unité 140 de détection de séquence pilote, illustrée en figure 7, permet ensuite de détecter et de dater le premier symbole ou tout autre symbole prédéterminé, de ladite séquence. Cette séquence pilote peut être 20 une séquence de symboles encodés ou non, communément appelé mot de synchronisation, ou par exemple un ton pur. Dans une étape E40, le couple de données (IO), Tm(i)) reçu de l'unité 130 est traité pour décoder chaque somme IO) sélectionnée, par décodage souple/dur pour estimer au moins une valeur binaire la plus 25 probable. Un décodage dur est utilisé pour traiter des données qui prennent un ensemble fixe et dénombrable de valeurs possibles (par exemple 0 ou 1), tandis qu'un décodage souple traite des données définies dans une gamme de valeurs représentant généralement une vraisemblance. 30 Ainsi, dans une sous-étape E41 de E40, le couple de données (I(i), Tm(i)) reçu de l'unité 130 est traité par un bloc 141 d'estimation de probabilité de l'unité 140, qui présente sur sa sortie la probabilité que I(n) 18 soit égal à 1 ou 0, dans une étape E41. Le bloc 141 d'estimation de probabilité est par exemple de type circuit limiteur (« slicer » en anglais). La probabilité est ensuite traitée, dans une sous-étape E42 de E40 par un bloc 142 de décodage en fonction du codage de canal qui a été utilisé 5 (turbo codage, LDPC, Viterbi, Reed-Solomon, etc.), L'unité 140 de détection permet également, dans une sous-étape E43 de E40, de traiter I(i) directement dans le cas où la séquence de symboles modulant le signal Stm n'est pas issue d'un codage correcteur d'erreur pour le canal de transmission. 10 Dans ce contexte, un bloc 143 de décision détermine la valeur du symbole émis en comparant I(i) à un seuil. Par exemple, pour un codage binaire où les symboles sont des bits, si I(i) est négatif, le bit associé est égal à 0, et dans le cas contraire égal à 1 (« hard bit decoding »). En fonction de la présence de ce codage de canal ou non, un bloc de 15 sélection 144, dans une étape E44, présente soit le signal décodé, soit le signal non décodé à un bloc 145 de synchronisation de trame. Dans une étape E45, le bloc 145 de synchronisation de trames est également configuré pour détecter une séquence pilote prédéfinie à détecter, dans une trame de longueur donnée Q bits. De façon classique, 20 cette détection se fait en calculant de façon glissante la corrélation de la séquence de bits avec la séquence pilote attendue, le critère de détection étant celui qui maximise la valeur absolue de cette corrélation Le bloc 145 de synchronisation de trames détecte donc tous les Q bits, la séquence pilote. En cas de détection, l'unité 140 peut présenter sur 25 sa sortie Thead(j) la date Tm(qt) du premier bit B(q) de la séquence détectée. Dans une étape E46, un bloc 146 de génération de paquet réceptionne ces données et peut générer alors un message contenant la date Thead(j), ainsi qu'un identificateur de paquet Id(j) (par exemple, un 30 numéro de compteur de trame). L'identificateur de paquet Id(j) sert de référence commune entre les paquets générés par des stations de réceptions distinctes. 19 Ce paquet est envoyé au centre de contrôle qui peut, à partir de ces informations, déterminer la position du mobile par exemple par trilatération, en comparant les dates des paquets, ayant le même identificateur Id(n), envoyés par au moins deux stations de réception de s signaux de télémesure. Le récepteur 100 de signaux radio, et le procédé proposé, en tirant avantage de la modulation cohérente, permettent de ne plus ré-échantillonner le signal à l'aide par exemple d'une boucle de récupération 10 de rythme pour conserver tout au long du traitement la cohérence entre échantillon et phase porteuse. Avantageusement, un tel récepteur 100/procédé permet de mesurer avec une grande précision la datation d'une séquence pilote, et améliore considérablement la variance sur la datation des transitions des symboles pour permettre entre autres une 15 meilleure localisation des engins mobiles émettant un tel signal de télémesure. Ainsi, la mise en oeuvre d'un tel récepteur 100/procédé permet d'atteindre un écart type de l'ordre de la nanoseconde (en considérant la base de temps comme parfaite), REVENDICATIONS1. Procédé de datation de réception de données numériques d'un signal modulé, ledit signal résultant de la modulation d'une porteuse ou s sous-porteuse par un signal numérique, le débit de symboles du signal numérique étant un sous-multiple entier N de la fréquence de la porteuse ou sous-porteuse, ledit procédé comprenant les étapes mises en oeuvre par un processeur d'un récepteur de télécommunication, consistant à, après réception et échantillonnage du signal modulé : démodulation cohérente au moyen d'une boucle de phase d'une pluralité d'échantillons sm(k) du signal modulé et obtention d'une pluralité d'échantillons démodulés sdm(k) et de la phase de la porteuse reconstituée ; 15 association de la pluralité échantillons démodulés sdm(k) avec la phase simultanée (p(k) de la porteuse reconstituée détermination d'une pluralité de dates de passage à une valeur déterminée cp0 de la phase de la porteuse reconstituée, relativement à la date de réception d'au moins un échantillon sm(k), 20 selon un décalage de temps pendant lequel la porteuse tourne de cp(k) à cp0. 2. Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel l'étape de démodulation 25 cohérente comprend les sous-étapes de : a. détection d'une erreur de phase par un détecteur de phase qui efface la modulation du signal reçu; b. filtrage passe-bas de l'erreur de phase. 30 3. Procédé de datation de réception de données numériques selon l'une des revendications précédentes pour lequel, dans l'étape de détermination des dates de passage à cp0, cette date est calculée 21 relativement aux échantillons sm(k) du signal modulé qui précédent ou suivent immédiatement un passage à p0 de la phase de la porteuse reconstituée. 4. Procédé de datation de réception de données numériques selon l'une des revendications précédentes comprenant en outre les étapes de : sélection d'une date de transition symbole parmi la pluralité de dates de passage à rp0, ladite date de transition symbole correspondant à un instant de transition de symboles du signal numérique ; décimation d'un facteur N des dates de passage à (p0, à partir de ladite date de transition symbole détectée ; et détermination des dates de réception d'une pluralité de symboles consécutifs correspondant aux dates de passage à rp() résultantes de la décimation. 5. Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel l'étape de sélection de la date de transition symbole, comprend une sous-étape d'intégration I(n) des échantillons démodulés sdm (k) sur une fenêtre de N échantillons consécutifs du signal. 6. Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel la sous-étape d'intégration I(n) est effectuée sur des durées d'un symbole démarrant à chaque échantillon sdm(k) du signal démodulé correspondant à un passage à (p0. 7. Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel la sous-étape d'intégration comporte également : 22 e. le calcul pour une pluralité d'indices consécutifs p de la valeur moyenne M(p) de la valeur absolue de ladite intégration I(m), pour m prend Ns valeurs qui valent toutes p modulo N, et f. la détermination de la date de transition symbole en 5 sélectionnant le passage à p0 d'indice p qui maximise M(p). 8. Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel la valeur de l'intégration I(n) démarrant à la date de transition symbole fournit le symbole 10 correspondant du signal numérique. 9. Procédé de datation de la réception de données numériques selon la revendication 7 dans lequel la comparaison des valeurs de M calculées aux passages à (p0 précédents M(n-1) et suivants M(n+1) 15 la transition de symbole détectée permet de caractériser une ambiguïté de phase liée à la modulation lorsque ces valeurs ne sont pas sensiblement égales mais que l'une est sensiblement égale à M(n). 20 10.Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel l'ambiguïté de phase est résolue en modifiant la phase de passage prédéterminée p0 à une autre valeur. 25 11.Procédé de datation de la réception de données numériques selon l'une des revendications 4 à 9 dans lequel le procédé comporte également le décodage des symboles en données binaires. 12.Procédé de datation de la réception de données numériques selon la 30 revendication précédente comportant en outre une étape de synchronisation de trames qui comprend une seconde décimation d'un facteur Q égal à une taille prédéterminée d'une trame de 23 données binaires, à partir d'une date de transition symbole, et une étape de détermination des dates de réception d'une pluralité de trames consécutives correspondant aux dates de passage à (p0 résultantes de la décimation. 5 13.Procédé de datation de la réception de données numériques selon la revendication précédente dans lequel l'étape de synchronisation de trames comporte en outre la sélection d'une trame ayant la meilleure corrélation avec une séquence pilote déterminée, parmi 10 une pluralité de trames consécutives et l'étape de détermination de la date de réception de ladite trame sélectionnée. 14.Récepteur (100) de télécommunication recevant un signal radio modulé, ledit signal résultant de la modulation d'une porteuse ou 15 sous-porteuse par un signal numérique, le débit de symboles du signal numérique étant un sous-multiple entier N de la fréquence de la porteuse ou sous-porteuse, ledit récepteur comprenant : une unité (120) de démodulation configurée pour, après réception et échantillonnage du signal modulé, la : 20 o démodulation cohérente au moyen d'une boucle de phase d'une pluralité d'échantillons sm(k) du signal modulé et obtention d'une pluralité d'échantillons démodulés sdm(k) et de la phase de la porteuse reconstituée ; o association de la pluralité échantillons démodulées sdm(k) 25 avec la phase simultanée cp(k) de la porteuse reconstituée o détermination d'une pluralité de dates de passage à une valeur déterminée (1)0 de la phase de la porteuse reconstruite, relativement à la date de réception d'au moins un échantillon sm(k). 30 24 15.Récepteur (100) de télécommunication recevant un signal radio modulé selon la revendication précédente, ledit récepteur comprenant également : - une unité (130) de synchronisation configurée pour la : s o sélection d'une date de transition symbole parmi la pluralité de dates de passage à (p0, ladite date de transition symbole correspondant à un instant de transition de symboles du signal numérique ; o décimation d'un facteur N des dates de passage à (p0, à partir 10 de ladite date de transition symbole détectée ; et o détermination des dates de réception d'une pluralité de symboles consécutifs correspondant aux dates de passage à p0 résultantes de la décimation.