Brevet FR2845546A1 — publié le 20040409

La présente invention concerne un procédé et un dispositif destinés à un estimateur de décalages de phase et de fréquence combiné pour une transmission par modulation par déplacement de phase multiple (appelée par la 5 suite "MPSK"), qui est capable de générer un plus petit retard par comparaison à un procédé classique estimant séquentiellement les décalages de phase et de fréquence en effectuant simultanément des estimations de décalages de phase et de fréquence en utilisant des estimations de 10 phases calculées à partir des échantillons dans chaque intervalle d'observation découpé, et qui est capable de réduire la complexité en partageant une pluralité de composants matériels. Ces dernières années, une technique de synchro15 nisation est devenue de manière variée nécessaire dans les traitements de transmission et de réception de signal d'une communication numérique. Les techniques de synchronisation peuvent être classées en une synchronisation en salve et une synchronisation continue. De nombreux systè20 mes de communication numérique modernes utilisent une transmission en mode de salve dans laquelle un nombre relativement petit de bits d'informations sont transmis en envoyant une séquence d'un petit nombre de formes d'onde. Dans une transmission par paquet en salve générale, les 25 données de transmission sont divisées en paquets, une séquence de symboles d'un préambule court est annexée au début de chaque paquet, et ensuite un récepteur démodule les symboles reçus à l'aide de la réalisation d'une synchronisation en utilisant le préambule court. Les données sont transmises en utilisant un mode de salve dans une transmission en mode duplex par répartition dans le temps (appelée par la suite '"TDD"') d'un système d'accès sans fil large bande fixe (appelé par la suite "BWA"). Du fait que le système BWA fixe uti35 lise une communication cohérente, il doit compenser un décalage de phase et un décalage de fréquence. Le décalage de phase est provoqué par la disparité des phases de référence d'un émetteur distant et d'un récepteur local. Le décalage de fréquence est provoqué par la disparité 5 des fréquences de l'émetteur distant et du récepteur local. En conséquence, le décalage de fréquence et le décalage de phase doivent être nécessairement compensés. Dans la communication en mode de salve, un court préambule est annexé dans la partie avant d'une partie de début de cha10 que transmission de salve et la salve ayant le court préambule est transmise. Afin d'obtenir un décalage de fréquence et un décalage de phase, un procédé synchrone anticipatif est en général utilisé. Dans un système typique, le décalage de phase 15 et le décalage de fréquence sont estimés d'une manière indépendante. A savoir, deux matériels exclusifs provoquant une complexité structurelle sont utilisés pour estimer le décalage de phase et le décalage de fréquence, c'est-àdire qu'un estimateur de phase de vraisemblance 20 maximale est utilisé pour estimer l'estimation de décalage de phase et un estimateur de décalage de fréquence assisté par des données est utilisé pour estimer l'estimation de décalage de fréquence. Du fait qu'une estimation de décalage de phase commence après une compensation 25 d'un décalage de fréquence, le temps d'attente croît dans un démodulateur. Le dispositif et le procédé classiques destinés à estimer une fréquence et une phase calculent les phases sur chacun des symboles reçus, et estiment ensuite un dé30 calage de phase et un décalage de fréquence en utilisant les phases calculées. Pour obtenir un décalage de phase, les phases calculées sont additionnées, et la somme des phases calculées est divisée par une valeur prédéterminée. Pour obtenir une estimation de décalage de fré35 quence, un intervalle d'observation de N symboles au to- tal est divisé en trois intervalles de longueur égale. Les phases du premier intervalle d'observation sont multipliées par une constante négative, les phases du deuxième sont multipliées par zéro, et les phases du troisième sont multipliées par une constante positive ayant la même quantité que la constante négative dans le premier intervalle d'observation de symboles. Les valeurs obtenues en résultat des multiplications dans les premier, deuxième et troisième intervalles sont additionnées 10 et le résultat est divisé par une valeur prédéterminée, estimant ainsi un décalage de fréquence. La présente invention est destiné à résoudre les problèmes précédents. C'est un but de la présente invention de fournir un procédé et un dispositif pour un 15 estimateur de décalages de phase et de fréquence combiné pour une transmission MPSK, o l'intervalle d'observation global est divisé en sous-intervalles adjacents et une estimation de phase séparée repliée est indépendamment obtenue à partir de chacun des sousintervalles et en20 suite est dépliée pour supprimer la discontinuité dans l'estimation de phase calculée. Le décalage de phase est ensuite calculé en mettant en moyenne les estimations de phases dépliées tandis que le décalage de fréquence est calculé en mettant en moyenne les-différences entre les 25 estimations de phases dépliées adjacentes. On fournit également un dispositif pour un estimateur de décalage de phase et de décalage de fréquence combiné pour la transmission MPSK, le dispositif étant caractérisé en ce qu'il comporte: un multiplicateur complexe pour multiplier un signal reçu par une conjuguée complexe d'un préambule, un accumulateur pour accumuler les sorties du multiplicateur complexe pour un intervalle d'observation divisé d'une longueur MT, un calculateur de phases pour calculer les phases des sorties de l'accumulateur pour obtenir N phases (o N est un nombre entier), une section de dépliement pour mettre en conti5 nu les N phases consécutives obtenues par le calculateur de phases, un estimateur de décalage de phase pour calculer une valeur moyenne des N phases dépliées provenant de la section de dépliement et estimer un décalage de phase 10 sur la base de la valeur moyenne après l'exécution d'une fonction (mod2z), et un estimateur de décalage de fréquence pour calculer une valeur moyenne de (N-1) différences de phases entre les estimations de phases dépliées adjacentes 15 provenant de la section de dépliement, et diviser la valeur moyenne de (N-l) différences de phases dépliées par 2nM afin d'estimer un décalage de fréquence. Les buts mentionnés ci-dessus ainsi que d'autres caractéristiques et avantages de la présente inven20 tion vont mieux apparaître dans la description en détail d'un mode préféré de réalisation de celle-ci en se reportant aux dessins annexés, sur lesquels: - la figure 1 est un schéma fonctionnel représentant une configuration d'un dispositif pour un estima25 teur de décalages de phase et de fréquence combiné pour une transmission MPSK selon le mode préféré de réalisation de la présente invention, - la figure 2 est un ordinogramme permettant d'illustrer un procédé pour un estimateur de décalages de 30 phase et de fréquence combiné pour une transmission MPSK selon le mode préféré de réalisation de la présente invention, - la figure 3 est une vue permettant de représenter un résultat qui compare un intervalle d'observa35 tion d'une estimation de décalage de phase selon la pré- sente invention à un intervalle d'observation classique d'une estimation de décalage de phase, - la figure 4 est un schéma fonctionnel illustrant un fonctionnement du dispositif représenté sur la figure 1, - la figure 5 est un schéma fonctionnel illustrant un fonctionnement d'un bloc d'estimateur de phase représenté sur la figure 4, et - la figure 6 est un ordinogramme permettant 10 d'illustrer un algorithme d'un bloc de dépliement représenté sur la figure 4. Par la suite, on va décrire plus en détail un procédé et un dispositif pour un estimateur de décalage de phase et de décalage de fréquence combiné pour une 15 transmission MPSK selon le mode préféré de réalisation de la présente invention en se reportant aux dessins annexés. Un dispositif pour un estimateur de décalage de phase et de décalage de fréquence combiné pour une trans20 mission MPSK sera décrit en se reportant aux figures 1, 2, 4, 5 et 6. La figure 1 est un schéma fonctionnel représentant une configuration d'un dispositif pour un estimateur de décalages de phase et de fréquence combiné pour une 25 transmission MPSK selon le mode préféré de réalisation de la présente invention. Le dispositif pour un estimateur de décalages de phase et de fréquence combiné pour une transmission MPSK met en oeuvre simultanément un estimateur de déca30 lage de phase 5 et.un estimateur de décalage de fréquence 6 en modifiant un estimateur de décalage de phase de vraisemblance maximale classique. En résultat, la présente invention génère un plus petit retard par comparaison à un procédé classique qui estime séquentiellement un 35 décalage de fréquence et un décalage de phase. Du fait que la présente invention partage une pluralité de composants matériels, elle réduit la complexité structurelle. Le dispositif inclut un multiplicateur complexe 1, un accumulateur 2, un calculateur de phases 3, une 5 section de dépliement 4, un estimateur de décalage de phase 5, et un estimateur de décalage de fréquence 6. Le multiplicateur complexe 1 multiplie un signal reçu r(kT) par une conjuguée complexe d'un préambule, ck. L'accumulateur 2 accumule une pluralité de 10 sorties du multiplicateur complexe 1 pour un intervalle d'observation divisé d'une longueur MT. Le calculateur de phases 3 calcule les phases des sorties provenant de l'accumulateur 2 pour obtenir N phases, o N est un nombre entier. La section de dépliement 4 prolonge les N 15 phases obtenues par le calculateur.de.phases 3 pour avoir des valeurs continues entre les estimations de phases adjacentes en utilisant un algorithme d'un ordinogramme d'un bloc de dépliement représenté sur la figure 6. L'estimateur de décalage de phase 5 calcule une 20 valeur moyenne des N phases de dépliement provenant de la section de dépliement 4, estime un décalage de phase en soumettant à la fonction (x mod27c) le résultat de la moyenne. L'estimateur de décalage de fréquence 6 calcule une valeur moyenne d'une différence de (N-1) estimations 25 de phases dépliées adjacentes provenant de la section de dépliement 4. L'estimateur de décalage de fréquence 6 divise la valeur moyenne de (N-1) différences de phases dépliées par 2nM afin d'estimer un décalage de fréquence. En d'autres termes, la présente invention es30 time un décalage de phase et un décalage de fréquence sur la base d'un décalage de phase dans un intervalle d'observation divisé. Dans un mode de réalisation. de la présente invention, on suppose qu'un cadencement de symboles est 35 idéal, qu'un décalage de fréquence normalisé est d'une manière significative inférieur à 1 (c'est-à-dire que fdT " 1), et la séquence de symboles reçue est une séquence connue. En conséquence, le k-ième signal reçu dans le multiplicateur complexe 1 est donné par r(kT)=c kej(2ndT+eo) + n(kT) (1) o ck est un signal MPSK reçu dans un k-ième intervalle de symboles, il est normalisé pour que IckI2=i, 00 est un décalage de phase, fd est un petit décalage de fréquence, 10 T est une période de symbole, n(kT) est un bruit de Gauss complexe blanc additif, dont les composantes réelle et imaginaire sont des variables aléatoires de Gauss indépendantes et réparties d'une manière identique ayant un écart nul et une variante de No 2E8 respectivement. Les composantes du signal et les composantes du bruit sont statistiquement indépendantes. Comme représenté sur la figure 2, on suppose qu'un intervalle d'observation total est LT. L'intervalle 20 d'observation total LT est divisé en N intervalles d'observation ayant chacun une longueur MT. Une estimation de phase pour chaque intervalle d'observation divisé est effectuée indépendamment. Un procédé pour un estimateur de décalages de 25 phase et de fréquence combiné pour une transmission MPSK selon le mode préféré de réalisation de la présente invention sera décrit en se reportant aux figures 2 et 3. La figure 2 montre un ordinogramme permettant d'illustrer un procédé pour un estimateur de décalages de phase et de 30 fréquence combiné pour une transmission MPSK selon le mode préféré de réalisation de la présente invention. La figure 3 est une vue permettant de comparer un intervalle d'observation d'une estimation de décalage de phase selon la présente invention à un intervalle d'observation classique d'une estimation de décalage de phase. Un intervalle d'observation total LT est donné. L'intervalle d'observation total LT est divisé en N in5 tervalles d'observation chacun ayant une longueur prédéterminée MT. Ici, N est un nombre entier. Une pluralité d'estimations de phases pour les N intervalles d'observation divisés sont indépendamment effectuées pour obtenir les estimations de phases. Les estimations de phase sont dépliées pour avoir une phase continue entre les estimations adjacentes. Un décalage de phase et un décalage de fréquence sont simultanément estimés en mettant en moyenne les estimations de phases dépliées et mettant en moyenne les 15 différences entres les estimations de phases dépliées adjacentes, respectivement. La figure 3(a) représente un intervalle d'observation d'une estimation de décalage de phase classique. Mais la figure 3(b) représente un procédé de divi20 sion de l'intervalle d'observation de l'estimation de décalages de phase et de fréquence selon la présente invention, c'est-à-dire que l'intervalle d'observation global LT est divisé en N intervalles d'observation de manière à ne pas être en chevauchement les uns avec les autres. Les 25 N intervalles d'observation divisés ont la même longueur, à savoir MT. Une estimation de phase est effectuée indépendamment sur chaque intervalle d'observation. La i-ième estimation de phase est indiquée par OM(i), o i est compris entre 0 et (N-1). La figure 4 est un schéma fonctionnel illustrant le fonctionnement du dispositif représenté sur la figure 1. Un bloc d'estimateur de phase 10 effectue indépendamment une pluralité d'estimations de phases pour les N intervalles d'observation divisés pour obtenir une plu35 ralité d'estimations de phases. L'estimation du i-ième intervalle d'observation divisé est indiquée par OM(i). La référence numérique 11 représente un bloc de dépliement. Le fonctionnement détaillé du bloc de dépliement 11 est décrit sur la figure 6. La référence numérique 12 re5 présente un bloc d'estimateur de décalage de phase. La référence numérique 13 représente un bloc d'estimateur de décalage de fréquence. Une estimation de phase pour un i-ième intervalle d'observation est exprimée par: f(i+l)M-1] OM(i) = arg E ckr(kT) (2) k=iM) o c est une conjuguée complexe de c, arg(x) est une fonction qui calcule une phase d'une valeur complexe de "x" et dont la valeur obtenue en résultat est dans la plage de [-s, 7]. La figure 5 est un schéma fonctionnel détaillé illustrant un fonctionnement d'un estimateur de phase représenté sur la figure 4. Comme représenté sur la figure 5, un multiplicateur complexe 1, un accumulateur 2, et un calculateur de phase 3 calculent séquentiellement chaque 20 partie de l'équation (2). En d'autres termes, le multiplicateur complexe 1 multiplie un signal reçu par une conjuguée complexe d'un k-ième signal de transmission connu. L'accumulateur 2 accumule les valeurs obtenues en résultat des multiplications provenant du multiplicateur 25 complexe 1 pour une période de temps MT, c'est-à-dire qu'il calcule (i+l)M-1] {Mc*r(kT)}. k=iM Le calculateur de phase 3 calcule les phases de la sortie provenant de l'accumulateur 2 pour obtenir une estimation 30 de phase, c'est-à-dire qu'il calcule [(i+l)M-Z] ar Ckr À argt E c*r(kT)f. k=,iM L'accumulateur de phase 2 est initialisé à une valeur complexe (0+ jO) au début de chaque accumulation pour le i-ième intervalle d'observation divisé. Ensuite, une accumulation des résultats de multiplication complexe pour 5 une période de temps MT et l'exécution du calcul de arg() permettent d'obtenir une estimation de phase. Ces procédés sont répétés jusqu'à ce que N estimations de phases soient calculées. En se reportant maintenant à la figure 4, une 10 phase OM(i) estimée par le bloc d'estimateur de phase 10 est limitée à une estimation de [-s, n] par le calcul de arg(). Dans une estimation de prolongation, lorsqu'une iième phase d'estimation se trouve très proche de _, une différence AOM(i) entre une i-ième et une (i-1)-ième 15 estimations de phases est en réalité très faible, mais devient d'environ 2n. Ceci fait qu'une importante erreur apparaît dans une estimation de décalage de phase et une estimation de décalage de fréquence. Afin d'éviter l'apparition d'une grande erreur, 20 une phase estimée courante calculée dans un intervalle d'observation courant doit avoir une valeur continue sur la base d'une estimation de phase précédemment calculée dans un intervalle d'observation divisé juste avant. Le bloc de dépliement 11 effectue une telle fonction. La 25 i-ième sortie du bloc de dépliement 11 est indiquée par AOM(i) À Une différence entre une i-ième et une (i-1)-ième phase d'estimation est calculée par l'équation (3) suivante: AOM(i) a O(i) - OM(i-1) (3) o une valeur initiale de OM(i) est Om(-1) et la valeur de OM(-1) est initialisée à zéro. C'est-à-dire que: OM(-1) = 0 (4) Afin d'avoir une estimation de phase continue 35 d'une valeur continue à une proximité proche de -n ou a, une différence de phase estimée est compensée selon l'équation (5) suivante: _AG (i), pour IAoM(i)Il(( A (i) {AO(i) - 21 sgn(AO.(i)), sinon 1, pour x 2 0 o sgn(x) = - 1, pour x(0 (6) Une phase dépliée dans un (i-l)-ième intervalle d'observation OM(i-1) est ajoutée à la i-ième différence d'estimation de phase compensée AOM(i) pour obtenir une phase dépliée dans un i-ième intervalle d'observation. La phase dépliée dans un i-ième intervalle d'observation 10 GM(i) est indiquée par l'équation suivante: OM(i) = OM(i-l)+A âM(i) (7) o une valeur initiale OM(-1) d'une phase de prolongation AOM(i) est initialisée à zéro. C'est-à-dire, ôM(-l) = 0 (8) Une sortie du bloc de dépliement 11 a la iième phase de prolongation ÈM(i) et une différence d'estimation de phase compensée AâM(i) calculées dans un procédé de calcul afin d'obtenir la iième phase de prolongation 20 OM(i). La figure 6 est un ordinogramme permettant d'illustrer un algorithme d'un bloc de dépliement représenté sur la figure 4. Le bloc d'estimateur de décalage de phase 12 calcule une valeur moyenne d'une pluralité 25 d'estimations de phases estimées dépliées pour obtenir une estimation de décalage de phase. L'estimation de décalage de phase est calculée par l'équation (9) suivante: 0O= {-: OM(i)}mod27r (9) Ni=O o y = xmod2z effectue l'addition d'un multiple entier de 2n à 'x' ou une soustraction de ce multiple à partir de x' afin d'obtenir 'y' dans une région de valeurs [-r, n]. Un bloc d'estimateur de décalage de fréquence 13 représenté sur la figure 4 calcule un décalage de fré5 quence estimé fdT en utilisant une différence d'estimation de phase de prolongation par l'équation (10) suivante. N-l _ i y. e(i) "dT = 2=M(N- 1) =(10) Lorsque L et M satisfont à l'équation (11) à un 10 rapport signal-bruit élevé, une estimation de décalage de fréquence ayant une variance minimale est implémentée, L =M + 2M fdT tan-'(MOfdT) (11) 2Moic tan-'(MorfdT)- 1 o M0 est une valeur optimale M et a un indice inférieur "0". Du fait que tan (x) donne approximativement 15 tan-l(x) x-1 pour une valeur de x petite, l'équation (11) est simplifiée en équation (12) suivante: L Mapp =- (12) o du fait que M est un nombre entier et est un facteur de L, M est sélectionné comme étant une valeur la plus 20 proche de L parmi les facteurs de L. La présente invention modifie le schéma d'estimation de phase classique nécessaire pour une transmis25 sion cohérente afin d'estimer le décalage de fréquence ainsi que le décalage de phase, de sorte qu'un calcul nécessaire supplémentaire est minimisé. De même, selon la présente invention, le dispositif pour un estimateur de décalage de phase et de décalage de fréquence combiné 30 pour une transmission MPSK est utilisé pour un estimateur de décalage de fréquence et un estimateur de décalage de phase précis en tant que composants principaux d'un récepteur numérique. Bien que la présente invention ait été particulièrement représentée et décrite en se reportant à un 5 mode particulier de réalisation de celleci, l'homme du métier doit comprendre que divers changements en termes de forme et de détail peuvent y être effectués sans s'écarter du domaine et de la portée de la présente invention telle que définie par les revendications an10 nexées. REVENDICATIONS 1. Procédé pour un estimateur de décalage de phase (5) et de décalage de fréquence (6) combiné pour une transmission par modulation par déplacement de phase 5 multiple, caractérisé en ce que le procédé comporte les étapes consistant à: (i) diviser l'intervalle d'observation total LT en N (o N est un nombre entier) intervalles d'observation chacun ayant une longueur prédéterminée MT (o M 10 est un nombre entier), (ii) effectuer indépendamment des estimations de phases pour les N intervalles d'observation divisés pour obtenir une pluralité d'estimations de phases, (iii) déplier la pluralité d'estimations de 15 phases obtenues à l'étape (ii) pour avoir des valeurs continues entre les estimations de phases adjacentes, et (iv) estimer simultanément un décalage de phase et un décalage de fréquence en mettant en moyenne les estimations de phases dépliées et en mettant en 20 moyenne les différences de phase entre les estimations de phases dépliées adjacentes, respectivement. 2. Procédé selon la revendication 1, caractérisé en ce que l'étape (ii) inclut les sous-étapes consistant à: diviser l'intervalle d'observation total LT divisé en N intervalles d'observation de la même longueur MT de manière à ne pas être en chevauchement les uns avec les autres, et effectuer indépendamment une estimation de 30 phase pour le i-ième intervalle d'observation en utilisant l'équation suivante (ig k)- ON(i) = arg{ E c*r(kT), k=iM o N est un nombre entier, T est une période de symboles, c est une conjuguée complexe de c, arg(x) est une fonc- tion qui calcule une phase d'une valeur complexe de "x" et a une valeur résultante de [-sr,]. 3. Procédé selon la revendication 1, caractérisé en ce que l'étape (iii) inclut les sous-étapes consis5 tant à: initialiser OM(-1) à zéro, initialiser AOM(-1) à zéro, calculer une différence entre la i-ième phase estimée et la (i-l)-ième phase estimée par l'équation 10 suivante: AOM(i) - OM(i) - om(i-l), compenser la différence de phase conformément à IAOM(i) I par l'équation suivante: ( AeK(i), pour |AeM(i)1 compensée pour obtenir une phase dépliée OM(i) dans un i-iéme intervalle d'observation qui est indiquée par l'équation suivante: OM(i) = OM(i-1) + A0M(i). 4. Procédé selon la revendication 1, caractéri20 se en ce que l'estimation d'un décalage de phase de l'étape (iv) calcule une valeur moyenne des estimations de phases dépliées pour obtenir une estimation de décalage de phase qui est indiquée par l'équation suivante: o = {N ê(i)} mod2ir, 0 8i ood Ni=0 o y = xmod2n signifie que y devient une valeur comprise dans [-sr,] en ajoutant un multiple entier de 2c à x ou en soustrayant celui-ci à partir de x. 5. Procédé selon la revendication 1, caractérisé en ce que l'estimation d'un décalage de fréquence de 30 l'étape (iv) calcule un décalage de fréquence fT en utilisant les différences de phase entre les estimations de phases dépliées adjacentes par l'équation suivante: N-1 2fdT 2M(N - l) i=1 6. Procédé selon les revendications i ou 2, caractérisé en ce que le procédé utilise un multiplicateur complexe (1) pour multiplier un signal reçu par une 5 conjuguée complexe d'un k-ième signal de transmission, un accumulateur (2) étant initialisé à une valeur complexe de (O+jO) au début d'une borne d'observation divisée et pour accumuler les valeurs résultant des multiplications complexes du multiplicateur complexe (1) pour la période 10 de temps MT, et un calculateur arg() (3), qui fonctionne une fois par intervalle d'observation divisé, pour calculer la phase d'une sortie complexe de l'accumulateur (2). 7. Procédé selon les revendications i ou 2, caractérisé en ce qu'une valeur optimale de M, représentée 15 par Mo, est calculée par l'équation: Mo + 2M27tfdT tan'(MO7rfdT) L = + 2Mo0r tan'(MnfdT - il est déterminé que M = Mo lorsque Mo est un nombre entier et un facteur de L, et sinon, M est déterminé comme étant une valeur la plus proche de Mo parmi les facteurs 20 de L. 8. Procédé selon les revendications 1 ou 2, caractérisé en ce qu'un intervalle d'observation optimal simplifié Mapp est calculé par l'équation suivante: L Mapp = 3 il est déterminé que M = Mapp lorsque Mapp est un nombre entier et un facteur de L, et sinon, M est déterminé comme étant une valeur la plus proche de Mapp parmi les facteurs de L. 9. Dispositif pour un estimateur de décalage de 30 phase (5) et de décalage de fréquence (6) combiné pour une transmission par modulation par déplacement de phase multiple, le dispositif étant caractérisé en ce qu'il comporte: un multiplicateur complexe (1) pour multiplier un signal reçu par une conjuguée complexe d'un symbole connu, un accumulateur (2) pour accumuler une plurali5 té de sorties du multiplicateur complexe (1) pour un intervalle d'observation divisé d'une longueur MT, un calculateur de phases (3) pour calculer les phases de la pluralité de sorties provenant de l'accumulateur (2) pour obtenir N phases, une section de dépliement (11) pour mettre en continu les N phases consécutives obtenues par le calculateur de phases (3) pour avoir des valeurs continues, un estimateur de décalage de phase (5) pour calculer une valeur moyenne des N phases dépliées par la 15 section de dépliement (11), et estimer une phase par une opération (mod2z) de la moyenne des N phases dépliées, et un estimateur de décalage de fréquence (6) pour calculer une valeur moyenne de (N-1) différences de phase entre les estimations de phases dépliées adjacentes pro20 venant de la section de dépliement (11), et diviser la valeur moyenne de (N-1) différences de phase par 2zM afin d'estimer un décalage de fréquence. 10. Dispositif selon la revendication 9, caractérisé en ce que l'estimateur de décalage de phase (5) 25 déplie une estimation de phase pour chaque intervalle d'observation divisé pour avoir des valeurs continues, calcule une valeur moyenne des estimations de phases dépliées et estime une phase par une opération (mod27t) de la moyenne des N phases dépliées. 30. 11. Dispositif selon la revendication 9, caractérisé en ce que l'estimateur de décalage de fréquence (6) déplie les estimations de phases pour chaque intervalle d'observation divisé pour avoir des valeurs continues, calcule une valeur moyenne des différences de phase 35 entre les estimations de phases adjacentes, et divise la valeur moyenne des différences de phase entre les estimations de phases dépliées adjacentes par 2nM afin d'estimer un décalage de fréquence.