FR 3063405 A1 20180831 FR 1751626 A 20170228 1 L'invention concerne un procédé d'estimation d'un écart d'horloge entre un émetteur et un récepteur durant une transmission sur un canal de communication d'un deuxième signal de communication issu d'une application à un premier signal de communication d'une technique d'étalement de spectre utilisant un signal, dit signal d'étalement, pour obtenir le deuxième signal de communication et un dispositif mettant en oeuvre le procédé. Au cours de la seconde guerre mondiale, des méthodes de communication sans fils utilisant une technologie dite par étalement de spectre, sont apparues. La technologie d'étalement de spectre avait deux motivations : résister aux efforts des ennemis pour brouiller un signal de communication, puis dissimuler la communication elle-même. La technologie d'étalement de spectre permet en effet d'obtenir des communications bas débit ayant une bonne robustesse. Les Figs. 4A, 4B, 4C, 5A, 5B et 5C illustrent la technologie d'étalement de spectre appliquée à un signal de données binaires. La Fig. 4A illustre schématiquement un signal 11(t) de données binaires sur lequel doit être appliqué un étalement de spectre. Dans l'exemple de la Fig. 4A, le signal 11(t) prend des valeurs de tension « 1» ou « -1 », la valeur de tension « 1» indiquant un « 1» binaire, la valeur de tension « -1» indiquant un « 0 » binaire. Dans le signal b(t), chaque bit a une durée Tb. La Fig. 4B illustre schématiquement un signal c(t) binaire pseudo aléatoire utilisé dans l'étalement de spectre. Dans l'exemple de la Fig. 4B, le signal c(t) prend des valeurs de tension « 1» ou « -1 », la valeur de tension « 1» indiquant un « 1» binaire, la valeur de tension « -1 » indiquant un « o» binaire. Dans le signal c(t) chaque bit a une durée Tc. avec la durée Tc. très inférieure à la duréeTb Tb). On note que le signal c(t) est parfois appelé signal d'étalement. On définit un facteur d'étalement FE tel que : Tb FE = - T, Les signaux d'étalement les plus courants sont les suivants : - les signaux basés sur des séquences pseudo-aléatoires de type séquences à longueur maximale ; 3063405 2 - les signaux basés sur des séquences de Walsh-Hadamard ; - les signaux basés sur des codes de Gold ; - les signaux basés sur des codes de Kasami. La Fig. 4C illustre schématiquement un signal obtenu après application d'un 5 étalement de spectre au signal de données binaires b(t). Lors de l'étalement de spectre, le signal b(t) est combiné au signal d'étalement c(t) : une combinaison d'une valeur de tension « 1» avec une valeur de tension « 1» donne une valeur de tension « 1 » ; une combinaison d'une valeur de tension « 1» avec une valeur de tension « -1» donne une valeur de tension « -1 » ; une combinaison d'une valeur de tension « -1 » avec une 10 valeur de tension « -1 » donne une valeur de tension « 1». Le signal résultant de l'étalement de spectre, ou signal étalé, avant modulation est noté e(t)=11(t)-c(t) On peut écrire : n=+co ( k=+co e(t) = b(nTb)8(t - nTb)) - ( c(kT,)8(t - kTe)) n=-co k=-co où n et k sont deux variables entières et Enn:tro 8(t - nTb) (respectivement 15 retro 8(t - kTe)) est un peigne de Dirac de période Tb (respectivement Te). La Fig. 5A illustre schématiquement une densité spectrale de puissance du signal de données binaires b(t). La densité spectrale de puissance, notée DSP(b(t)), prend la forme d'un lobe principal 50 centré sur une fréquence centrale et ayant une largeur de bande de fréquences de -2 entouré de lobes secondaires (51, 52) ayant une largeur de Tb 1 20 bande de fréquences de -Tb. La Fig. 5B illustre schématiquement une densité spectrale de puissance du signal binaire pseudo-aléatoire c(t). La densité spectrale de puissance, notée DSP(c(t)), prend la forme d'un lobe principal centré sur une fréquence centrale et ayant une largeur de bande de fréquences de -T2 entouré de lobes secondaires ayant une largeur de bande de 25 fréquences de La Fig. 5C illustre schématiquement une densité spectrale de puissance du signal étalé e(t). La densité spectrale de puissance, notée DSP(e(t)), prend la forme d'un lobe principal centré sur une fréquence centrale et ayant une largeur de bande de fréquences de -2 entouré de lobes secondaires ayant une largeur de bande de fréquences de -. On Tc Tc 30 note qu'après application de l'étalement de spectre sur le signal de données binaires b(t), le signal étalé e(t) conserve la même énergie que le signal b(t) mais répartie sur une bande de fréquence plus grande. 3063405 3 Nous supposons que le signal étalé est modulé à l'aide d'un schéma de modulation BPSK (modulation par changement de phase binaire, « Binary Phase Shift Keying » en terminologie anglo-saxonne). On note : tx(t) = BPSK(e(t)) 5 où BPSK() représente une fonction de modulation BPSK. On remarque qu'ici le signal e(t) est implicitement modulé suivant une modulation BPSK puisque les valeurs prises par e(t) sont « -1» et « +1 » ce qui correspond à une définition habituellement choisie pour la modulation BPSK. On pourrait donc écrire : tx(t) = e(t) 10 Suite à une transmission du signal étalé modulé tx(t) sur un canal de communication sans fils par un dispositif émetteur, un signal rx(t) est reçu par un dispositif récepteur. Lorsque le canal de communication est gaussien, le signal reçu rx(t) peut s'écrire Eb rx(t) = tx(t) + (-No) N(t) 15 où n (3E. ) désigne un coefficient de pondération fonction du rapport énergie par No bit sur densité spectrale de puissance du bruit -Eb (« energy per bit to noise power spectral 0 density ratio » en terminologie anglo-saxonne). Le rapport -Eb est donc une métrique No représentative du bruit ayant affecté le signal tx(t) lors de sa transmission. N(t) désigne une variable aléatoire gaussienne complexe d'espérance nulle. 20 Lorsque les dispositifs émetteur et récepteur sont synchronisés en fréquence et en temps, l'opération inverse de l'étalement, appelée étalement inverse est appliquée au signal reçu comme suit : k=+00 d(t) = rx(t) - c(kT,)8(t - kn)) Le signal b(t) est estimé de la manière suivante : k=(n+1)FE-1 25 b(nTb) = BPSK-1 rx(kn).c(kTc) ,VnEZ k=n.FE où BPSK-10 représente l'opération de démodulation BPSK et le symbole - signifie estimée. Le signal b(t) est par exemple une estimation du signal b(t). L'estimation du signal b(t) est donc calculée à partir d'un signal bruité ce qui peut entraîner une mauvaise estimation du signal b(t). ). Par ailleurs, l'hypothèse de 30 synchronisation en temps et fréquences des dispositifs émetteur et récepteur n'est pas 3063405 4 forcément vérifiée. En effet, chaque dispositif possède une horloge interne permettant audit dispositif d'effectuer des mesures de temps nécessaires au fonctionnement dudit dispositif. Les horloges internes sont en générale très précises, mais on observe des fluctuations de cette précision dues à de nombreux facteurs tels que des variations de 5 température, des variations d'altitude ou dues à la mobilité des dispositifs provoquant un effet Doppler. Ces fluctuations de la précision des horloges n'affectant pas de manière identique tous les dispositifs d'un système de communication, elles peuvent provoquer des écarts d'horloge entre lesdits dispositifs et donc des désynchronisations. Il est souhaitable de pallier ces inconvénients de l'état de la technique. Il est 10 notamment souhaitable de proposer une méthode qui permette de resynchroniser temporellement les dispositifs émetteur et récepteur. Il est par ailleurs souhaitable de proposer une méthode qui soit simple à mettre en oeuvre et à faible coût. Selon un premier aspect de la présente invention, la présente invention concerne 15 un procédé d'estimation d'un écart d'horloge entre un émetteur et un récepteur durant une transmission sur un canal de communication d'un deuxième signal de communication résultant d'une application à un premier signal de communication d'une technique d'étalement de spectre utilisant un signal, dit signal d'étalement. Le procédé comprend suite à la réception du second signal par le récepteur : obtenir au moins un 20 signal, dit signal d'étalement orthogonal, orthogonal au signal d'étalement; pour chaque signal d'étalement orthogonal obtenu, appliquer ledit signal d'étalement orthogonal obtenu au second signal reçu par le récepteur pour obtenir un troisième signal ; obtenir, une estimation d'un rapport énergie par bit sur densité spectrale de puissance d'un bruit de communication ayant affecté le deuxième signal durant sa transmission à partir de 25 chaque troisième signal ; et, estimer une information représentative de l'écart d'horloge entre l'émetteur et le récepteur en fonction de l'estimation dudit rapport. L'estimation de l'écart d'horloge permet de resynchroniser temporellement les dispositifs émetteur et récepteur. Selon un mode de réalisation, lorsqu'un seul signal d'étalement orthogonal est 30 déterminé, ledit rapport est estimé à partir d'une fonction prédéterminée d'une moyenne arithmétique d'informations élevées au carré, lesdites informations étant représentatives du bruit de communication ayant affecté le deuxième signal durant sa transmission et ayant été obtenues à partir du troisième signal. 3063405 5 Selon un mode de réalisation, lorsque une pluralité de signaux d'étalement orthogonaux est obtenue, pour chaque signal d'étalement orthogonal de la pluralité, une première moyenne arithmétique est calculée sur des informations élevées au carré, lesdites informations étant représentatives du bruit de communication ayant affecté le 5 deuxième signal durant sa transmission et ayant été obtenues à partir du troisième signal obtenu pour ledit signal d'étalement orthogonal, ledit rapport étant estimé à partir d'une fonction prédéterminée d'une seconde moyenne arithmétique calculée à partir des premières moyennes arithmétiques calculées pour chaque signal d'étalement orthogonal. 10 L'estimation de l'écart d'horloge est ainsi fiabilisée. Selon un mode de réalisation, l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur comprend une estimation de la valeur absolue dudit écart d'horloge. Selon un deuxième aspect de l'invention, l'invention concerne un procédé itératif 15 d'estimation d'un écart d'horloge entre un émetteur et un récepteur. Le procédé selon le premier aspect est mis en oeuvre régulièrement sur le deuxième signal, l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur étant remise à jour à chaque itération. De cette manière, l'écart d'horloge entre l'émetteur et le récepteur est 20 constamment corrigé. Selon un mode de réalisation, l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur comprend en outre un signe de l'écart d'horloge, ledit signe étant déterminé en comparant un écart d'horloge d'une itération courante avec un écart d'horloge d'une itération précédente, le signe de l'écart d'horloge de l'itération 25 courante étant l'inverse du signe de l'écart d'horloge de l'itération précédente lorsque l'écart d'horloge augmente et le signe de l'écart d'horloge de l'itération courante étant identique au signe de l'écart d'horloge de l'itération précédente lorsque l'écart d'horloge n'augmente pas. Selon un troisième aspect de l'invention, l'invention concerne un dispositif 30 comprenant des moyens pour mettre en oeuvre le procédé selon le premier ou le deuxième aspect. Selon un quatrième aspect de l'invention, l'invention concerne un programme d'ordinateur, comprenant des instructions pour mettre en oeuvre, par un dispositif, le 3063405 6 procédé selon le premier ou le deuxième aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif Selon un cinquième aspect de l'invention, l'invention concerne des moyens de stockage stockant un programme d'ordinateur comprenant des instructions pour mettre 5 en oeuvre, par un dispositif, le procédé selon le premier ou le deuxième aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description suivante d'un exemple de réalisation, ladite description étant faite en relation avec les dessins joints, parmi 10 lesquels : - La Fig. 1 illustre schématiquement un système de communication dans lequel est mise en oeuvre l'invention; - la Fig. 2 illustre schématiquement un module de traitement compris dans un dispositif de réception ; 15 - la Fig. 3A illustre schématiquement un procédé selon l'invention d'estimation d'un écart d'horloge entre un émetteur et un récepteur; - la Fig. 3B illustre schématiquement un procédé selon l'invention d'estimation d'une information représentative d'un bruit de communication ayant affecté un signal transmis en utilisant une technologie d'étalement de spectre; 20 - la Fig. 4A illustre schématiquement un signal de données binaires sur lequel doit être appliqué un étalement de spectre; - la Fig. 4B illustre schématiquement un signal binaire pseudo aléatoire utilisé dans l'étalement de spectre; - la Fig. 4C illustre schématiquement un signal obtenu après application d'un 25 étalement de spectre au signal de données binaires; - la Fig. 5A illustre schématiquement une densité spectrale de puissance du signal de données binaires ; - la Fig. 5B illustre schématiquement une densité spectrale de puissance du signal binaire pseudo-aléatoire ; et, 30 - la Fig. 5C illustre schématiquement une densité spectrale de puissance du signal de données binaires sur lequel a été appliqué un étalement de spectre. La Fig. 1 illustre schématiquement un système de communication 1 dans lequel est mise en oeuvre l'invention. 3063405 7 Le système de communication 1 comprend un dispositif émetteur 10 et un dispositif récepteur 11. Le dispositif émetteur 10 transmet le signal tx(0 tel que décrit plus haut, i.e. le signal tx(0 résultant d'une modulation BPSK du signal étalé e(t, le signal étalé e(t résultant de l'application du signal d'étalement c(t au signal de données 5 binaires b(t. Le signal tx(0 est transmis sur un canal de communication sans fils 12 utilisant par exemple les bandes de fréquences connues sous l'appellation « Bande ISM» (Industrie, Science et Médical) comprenant des bandes de fréquences pouvant être utilisées librement pour des applications industrielles, scientifiques et médicales. Le dispositif récepteur 11 reçoit le signal rx(t, tel que décrit plus haut, résultant 10 de la transmission du signal tx(0 sur le canal de transmission sans fil 12. Le dispositif récepteur 11 comprend un module de traitement 110, décrit en relation avec la Fig. 2, apte notamment à fournir une estimation du signal b(t. La Fig. 2 illustre schématiquement un exemple d'architecture matérielle du module de traitement 110 compris dans le dispositif récepteur 11. 15 Selon l'exemple d'architecture matérielle représenté à la Fig. 2, le module de traitement 110 comprend alors, reliés par un bus de communication 1100 : un processeur ou CPU (« Central Processing Unit » en anglais) 1101 ; une mémoire vive RAM (« Random Access Memory » en anglais) 1102; une mémoire morte ROM (« Read Only Memory » en anglais) 1103 ; une unité de stockage telle qu'un disque dur 20 ou un lecteur de support de stockage, tel qu'un lecteur de cartes SD (« Secure Digital » en anglais) 1104 ; au moins une interface de communication 1105 permettant au module de traitement 110 de communiquer avec d'autres modules ou dispositifs. Par exemple l'interface de communication 1105 permet au module de traitement 110 de communiquer avec le dispositif émetteur 10. 25 Le processeur 1101 est capable d'exécuter des instructions chargées dans la RAM 1002 à partir de la ROM 1003, d'une mémoire externe (non représentée), d'un support de stockage (tel qu'une carte SD), ou d'un réseau de communication. Lorsque le dispositif récepteur 11 est mis sous tension, le processeur 1101 est capable de lire de la RAM 1102 des instructions et de les exécuter. Dans un mode de 30 réalisation, ces instructions forment un programme d'ordinateur causant la mise en oeuvre complète ou partielle, par le processeur 1101, des procédés décrits ci-après en relation avec les Figs. 3A et 3B. Les procédés décrits en relation avec les Figs. 3A et 3B peuvent être implémentés sous forme logicielle par exécution d'un ensemble d'instructions par une machine 3063405 8 programmable, par exemple un DSP (« Digital Signal Processor » en anglais) ou un microcontrôleur, ou être implémentés sous forme matérielle par une machine ou un composant dédié, par exemple un FPGA (« Field-Programmable Gate Array » en anglais) ou un ASIC (« Application-Specific Integrated Circuit » en anglais). 5 La Fig. 3A illustre schématiquement un procédé selon l'invention d'estimation d'un écart d'horloge entre l'émetteur 10 et le récepteur 11. Dans un mode de réalisation, le procédé décrit en relation avec la Fig. 3A est mis en oeuvre à intervalles réguliers de manière à corriger régulièrement les éventuelles dérives d'horloge entre le dispositif émetteur 10 et le dispositif récepteur 11. 10 On suppose ici que le module de traitement 110 du dispositif récepteur reçoit en continu le signal rx(t). Dans une étape 31, le module de traitement 110 obtient une information représentative du bruit de communication ayant affecté le signal tx(t) durant sa transmission. Pour ce faire le module de traitement 110 met en oeuvre un procédé décrit 15 en relation avec la Fig. 3B. La Fig. 3B illustre schématiquement un procédé selon l'invention d'estimation d'une information représentative d'un bruit de communication ayant affecté un signal transmis en utilisant une technologie d'étalement de spectre. Dans une étape 311, le module de traitement 110 obtient un signal 20 ci(t) orthogonal au signal d'étalement c(t) que nous appelons signal d'étalement orthogonal. Dans un mode de réalisation, le signal d'étalement c(t) est obtenu à partir d'une séquence de Walsh-Hadamard, elle-même obtenue en utilisant, par exemple, une matrice de Hadamard, notée H6 4 , de taille 64x64. 25 Une matrice de Hadamard se construit récursivement. Soit H1 une matrice de Hadamard de taille 1 x 1, H2 une matrice de Hadamard de taille 2x2 et H4 une matrice de Hadamard de taille 4x4 : 30 H 4 = H2= -1 H2n--H2n-= [ 1-] -1 On peut donc écrire : H2n' [1_ ri il = H2 g^H2n- H2n = [n 2n- 1 [i 3063405 9 où le symbole 0 désigne le produit de Kronecker. Le signal d'étalement c(t) est obtenu en faisant un OU exclusif entre une séquence de Walsh-Hadamard et une séquence PN pseudo aléatoire. Le signal d'étalement c(t) est obtenu, par exemple, à partir d'une séquence de Walsh-Hadamard 5 donnée par une ligne de la matrice de Hadamard H64. Dans ce cas, le signal d'étalement orthogonal ci(t) est obtenu à partir d'une séquence de Walsh-Hadamard donnée par une des « 63 » lignes restantes de la matrice de Hadamard H64. Dans une étape 312, le module de traitement 110 applique le signal d'étalement orthogonal ci (t) au signal reçu rx(t). =+ co 10 di(t) = rx(t) - ci(kT,)8(t - kT,)) k=-co Le module de traitement 110 obtient alors un signal di(t). Il est possible d'introduire le signal reçu rx(t) dans la somme. On obtient alors : (k=+co di (t) = rx(kTc)c±(kT,)8(t - kT,)) /c. Dans une étape 313, le module de traitement 110 obtient une estimation du rapport Eb 15 ' qui est une information représentative du bruit de communication ayant affecté le No signal tx(t) durant sa transmission, à partir du signal cl-'-(t). Soit une quantité Q±(t) définie de la manière suivante : k=(n+1)FE-1 QI (nTb) = rx(kTc)c±(kTc) k=nFE k=(n+1)FE-1 k=nFE 20 Le coefficient de pondération n ) peut être extrait de la somme puisqu'il ne No dépend pas de la variable k. On obtient alors : k=(n+1)FE-1 Eb Eb (2±(nTo) = (-No) - N (kTc). CI (kTc) = 11 ()s'(nTb)vn E k=nFE Z où S-L(nTb) = Ekkin(nF+E1)F E -1 N(kTc)c±(kTc),Vn E E. Comme nous l'avons vu plus haut, N(t) (et donc N(kTc)) désigne une variable 25 aléatoire gaussienne complexe d'espérance nulle. S±(nTo) est aussi une variable aléatoire gaussienne complexe d'espérance nulle de sorte que la quantité (2' (nTo) peut être considérée comme une quantité (i.e. une information) représentative du bruit de (tx(kTc) + N(kTc)) ci(kTc) 3063405 10 communication ayant affecté le signal tx(t) durant sa transmission. Soit une quantité e(P) correspondant à une moyenne arithmétique de la quantité (2±(nTb) élevée au carré sur P dernières valeurs disponibles de (2±(nTb), chaque valeur disponible de la quantité (2±(nTb) correspondant à un symbole reçu par le module de traitement 110 et 5 issu du signal rx( t). p=n 1 1(2±(PTb)12 p=n-P+1 En développant la moyenne arithmétique (1)(P) on obtient On note donc qu'il est possible d'obtenir une relation entre la moyenne 10 arithmétique e(P) et le rapport énergie par bit sur densité spectrale de puissance du bruit Lb. Il est alors possible d'obtenir une estimation () dudit rapport en utilisant No No une fonction rl de la moyenne arithmétique e(P). Dans un mode de réalisation, la fonction rl est une fonction du logarithme népérien de e(P), noté log(P(P)), ce qui permet d'obtenir une fonction rl représentative d'une droite affine dans un espace à 15 deux dimensions. (1) = i (i0g(e(P))) = a1i0g(e(P)) b1 No al et b1 sont des coefficients prédéterminés de la fonction Ti et l'estimation (Eb(1) est donnée en dB. 40) E (1) Lors de l'étape 313, le module de traitement 110 obtient l'estimation () du No 20 rapport -Eb en introduisant le logarithme népérien de la moyenne arithmétique e(P) dans No la fonction Ti dont les paramètres al et b1 ont été prédéterminés. Comme décrit plus haut, la moyenne arithmétique e(P) est obtenue en appliquant le signal d'étalement (t) aux P derniers symboles reçus issus du signal rx(t). Dans un mode de réalisation, les coefficients al et b1 de la fonction rl du 25 logarithme népérien de e(P) sont prédéterminés par simulation. Dans ce mode de réalisation, une séquence PN de type MILS (Séquence à longueur maximale : «Maximum Length Sequence » en terminologie anglo-saxonne) obtenue à l'aide d'un registre à décalage à rétroaction linéaire, est utilisée. La séquence PN est produite par exemple à partir d'un polynôme primitif de degré p, le degré p étant tel que la longueur 2 p=n «ro =1 (Eb)) P No ISI(gb)12 p=n-P+1 3063405 11 de la séquence MLS produite (i.e. 2P -1), soit très supérieure à la longueur de la séquence de Walsh-Hadamard choisie. Une transmission sur un canal Gaussien d'un signal modulé obtenu en appliquant le signal d'étalement c(t) au signal b(t) et en appliquant une modulation BPSK au signal étalé obtenu est simulée. Une démodulation 5 et un étalement inverse par un signal d'étalement orthogonal c.1(0 orthogonal à c(t) sont appliqués au signal issu de la transmission simulée. Par exemple, on obtient par simulation al = - -133 et b1 = O. Dans un mode de réalisation, l'estimation du rapport -Eb est améliorée en itérant plusieurs fois les étapes 311 et 312 du procédé décrit en relation avec la Fig. 3B. A 10 chaque itération un signal d'étalement orthogonal différent est utilisé. Dans le mode de réalisation où le signal d'étalement c(t) est obtenu à partir d'une séquence de Walsh-Hadamard donnée par une ligne de la matrice de Hadamard H64, chaque signal d'étalement orthogonal est obtenu à partir d'une séquence de Walsh-Hadamard donnée par une des « 63 » lignes restantes de la matrice de Hadamard H64, différente de la ligne 15 de la matrice de Hadamard H64 utilisée pour générer le signal d'étalement c(t). Dans ce cas de V=1 à V=63 signaux d'étalement orthogonaux Cd(t) différents (u E [1; V]) peuvent être obtenus à partir de la matrice de Hadamard H64. A chaque itération des étapes 311 et 312 du procédé décrit en relation avec la Fig. 3B, i.e. pour chaque signal d'étalement orthogonal, une quantité Q± (nTb), 20 représentative du bruit de communication ayant affecté le signal tx(t) durant sa transmission, est obtenue par le module de traitement 110 : k=(n+1)FE-1 (2,±(nTb) = rx(kTc)c'±(kTc) k=nFE Suite à chaque itération des étapes 311 et 312, lors de l'étape 313, le module de traitement 110 calcule une moyenne arithmétique eu (P) p=n 1 25 (P) = - iQ(P71)12 p=n-P+1 Lors de l'étape 313, le module de traitement 110 calcule une quantité tp(P) qui est une moyenne arithmétique calculée sur les V moyennes arithmétiques eu(P) obtenues à chaque itération des étapes 311 et 312 en appliquant les V signaux d'étalement orthogonaux Cd (t) différents aux P symboles issus du signal reçu rx(t). 1 30 V ci-(t)E[C1} 3063405 12 où {C±} est l'ensemble de V signaux d'étalement orthogonaux Cu-L(t). Il est alors possible d'obtenir une estimation () du rapport E) en utilisant une No fonction r2 de la moyenne arithmétique p(P). Dans un mode de réalisation, la fonction r2 est une fonction du logarithme népérien de tP(P) ce qui permet d'obtenir une fonction 5 r2 représentative d'une droite affine dans un espace à deux dimensions. (Eb -No) = r2(10g(1P(P)))= a210g(e(P))+ b2 a2 et b2 sont des coefficients prédéterminés de la fonction r2 et l'estimation (2) 10 (E'b (2) 40) est donnée VAR ((Ei) -(2)) No (1) (El)) No en dB. L'estimation est associée à une variance No inférieure à une variance VAR (()) associée à l'estimation --E-(1)) No (2) (1) une information VAR(() - I .VAR(rb) ) \No/ V \No/ Dans une étape 33, le module de traitement 110 estime représentative d'un écart d'horloge E entre le dispositif émetteur 10 et le dispositif récepteur 11. 15 Soit 12(nTb)c et 12(nTb)[c1l deux quantités utilisées dans la détermination de l'écart d'horloge E: k=(n+1)FE-1 n,(nTb)c rx(kTc + E)c(kTc) k=nFE 20 çl(nTb)tcll = k=(n+1)FE-1 rX(kTc E)C7±(kTc) ci-Efc+} k=nFE En appliquant un développement de Taylor au signal rx(t) au voisinage de kTc, on obtient : 25 rX(kTc E) = rX(kTc) E - rx'(kTc) Les quantités f/(nTb)c et f/(nTb)[cil s'écrivent alors : 3063405 13 k=(n+1)FE-1 n,(nTb) = c(kTc)(rx(kTc) + E - rXi (kTc)) k=nFE k=(n+1)FE-1 l/(nTb)[cil = ci±(kTc) (rx(kTc) + E - rXi (OC)) CjjE[C-} k=nFE 5 Avec : rx(kTc) = (tx * hRCS)(kTc) )(N * hRcs)(kTc) in=-FMRcs rX(kTc) = ilcs(0)BPSekTc) + 11 (E) hRcs(MTc) N ((k - M)Tc) No in=-FMcs rXi (kTc) = h'c-s(MTc)BPSK((k - M)Tc) m=-Mcs rn=-EMRcs + 7 1 11;?c-s (MT c) N ((k - Mec) 0 RCS 10 Où: - * : désigne l'opération de convolution - hRcs : désigne un filtre en racine de cosinus surélevé (comprenant 2MRc's + 1 coefficients) --1?' CS = dhRcs 15 - h dt (i.e. dérivée première de hRcs) - h -cs = hRcs * hRcs : désigne un filtre en cosinus surélevé (comprenant 2Mc5 + 1 coefficients) obtenu à partir du filtre en racine de cosinus surélevé hRcs - hc' s = -dhcs (i.e. dérivée première de hcs) dt 20 Nous avons alors : Eb S1(11,Tb) = him,(0)BPSK(11,71)FE + (ATD)S(nTb) E Eb Eb Eb - (D(nTb) + (-No)S'(nTb)) = Ac + E - Bc IVO -MR CS 3063405 14 Eb Eb n(nTbhcil = (ATD)s±(01) + E - n (k)s_UnTb) = Atc11 (Z) + E - B[c11 ()-No Avec : k=(n+1)FE-1 m=+Mcs D (nT b) = c(kTc) h:cs(mTc)BPSK((k - m)Tc) k=nFE m-Mcs 5 k=(n+1)FE-1 m=+MRcs S(nTb) = c (k Tc) hRcs(mTc) N ((k - m)Tc) k=nFE in=-MRCS k=(n+1)FE-1 rn=+mRcs S'(nTb) = c(kTc) h'Rcs(mTc) N ((k - m)Tc) k=nFE m=-MRcs k=(n+1)FE-1 rn=-EMRcs S _L(nT b) = ci±(kTc) hRcs(mTc) N ((k - m)Tc) ci-Efc-L} k=nFE m--mRcs k=(n+1)FE-1 in-+MRcs S 1(nT b) = c%+ Y h.scs(mTc)N((k -m)Tc) k=nFE m--mEcs Nous constatons que les quantités 12(nTb)c et (nTb)[cil dépendent de l'écart E. Pour déterminer l'écart E, nous formons les quantités auxiliaires suivantes : k=n 1 Inc(01)12 If2I= P p k=n-P+1 k=n 2 1 2 Infcl} IF = Ifitcil (01) k=n-P+1 où P est un nombre entier positif prédéfini. On effectue donc une somme des 20 valeurs absolues élevées au carré des P dernières valeurs des quantités l/.(nTb)c et l/.(nTb)[cil. Pour un rapport (-Eb) donné, les quantités auxiliaires peuvent être No modélisées à l'aide d'une fonction polynômiale du second degré de la manière suivante : 10 15 3063405 15 112c 1r,2 = ctE2 + aïE + ag = f (E) 12 [cl} [cl} [cl} 112fci}ip = a2 £2 +a1 £ ao =9(E) 5 f() et g0 sont des fonctions polynômiales paires. Par conséquent, seule la valeur [} absolue de l'écart 1E1 peut être mesurée et aï = a1 = 0. ag, ao[cll [cl} Eb et a2 sont des coefficients fonction du rapport (ou d'un No S R Eb rapport -s puisque - = - - où R et W désignent respectivement un taux binaire et une N W No [cl} [} bande passante du récepteur). Chaque coefficient ag, ao et ac1 2 est décrit par un 10 modèle mathématique prédéterminé à l'aide d'un simulateur. ,c _ (E b), ac _ hc (E b), atc11 _ E b c at = c 1 Eb cl L'e0 11.0 No 2 2 No 1\i0) 2 2 No Le simulateur comprend un émetteur, un canal (gaussien) et un récepteur. Le signal reçu, sur-échantillonné avec une période de sur-échantillonnage donnée, est supposé parfaitement synchronisé en fréquence (c'est-à-dire que l'écart de fréquence 15 entre le signal émis et le signal reçu est nul). Une première plage de variation est définie pour le rapport -Eb (ou -s) et une deuxième plage de variation est définie pour l'écart No temporel £ (exprimé en multiple entier, éventuellement nul, de la période de sur-échantillonnage). Des simulations sont effectuées à l'aide de deux boucles imbriquées : une première boucle permet de faire varier le rapport -Eb (ou -s) sur la première plage de No 20 variation et une deuxième boucle permet de faire varier l'écart temporel £ sur la deuxième plage de variation. L'objectif est de déterminer les quantités auxiliaires inc 2 et 112fc1} I _. Pour chaque valeur du rapport (ou -s) de la première plage de variation, No 2 deux ensembles finis de valeurs des quantités auxiliaires Ificir,2 et 112fc1} IF sont obtenus en faisant varier l'écart temporel £ sur la deuxième plage de variation. Le cardinal de 25 ces deux ensembles est donc égal au nombre d'écarts temporels c simulés. Les deux ensembles sont ensuite interpolés en utilisant les fonctions polynômiales f et g, pour obtenir pour chaque valeur du rapport -Eb (ou -s) deux couples de coefficients: No N (a,c, (e), a (Eb)) , 2 No et anIc1} ()E ), a.,1c1} (El') . L'ensemble de ces couples forme deux No No tableaux à deux colonnes dont le nombre de lignes correspond au nombre de valeurs 30 simulées du rapport -Eb (ou -s). No 3063405 16 L'interpolation des valeurs des colonnes de ces tableaux à l'aide d'autres fonctions polynômiales h (de degré > 2) permet de modéliser respectivement les coefficients aoc, a2c, [cl} et a2 en fonction du rapport signal sur bruit et d'obtenir les modèles mathématiques suivants : ag c11 = = hç (Eb) at = htc 1 (Eb) et alcil = ° ) 2 0 » 0 No 2 5 1112c1} (ENb3 Soit un estimateur «E) : intA2 (E) ht2c1} GvEbc) E2 + 7 n ,}(Eb) No .(E) = P ' (Eb)E2 _L hc (Eb) No E bE b[i E b[Ci (E b 10 Une fois les fonctions hg () () hCj déterminées, il est No 2 No 0 (No) ' hj 2 No) possible de calculer des valeurs numériques de l'estimateur «E) ce qui permet de déduire VI: (E\DP) - hg e,(,) (NEbo) (E) hr} (NEb) Afin de diminuer la complexité de l'estimation, il est possible de simplifier l'estimateur «E) ci-dessus en ne considérant que l'une ou l'autre des quantités Sl(nTb)c et l/(nTb)[cil 2 Dans un mode de réalisation, la quantité 112fci} IF, déterminée à l'aide d'un 20 ensemble {c-L} réduit est utilisée : .(E) = In fc1}1132 = 9(E) = (E) E 2 + hIcil (4) 2 No ° \Nof Comme précédemment, la connaissance de la valeur numérique .(E) permet de 25 déduire VI: VI 15 3063405 17 VI «E) - hfc1} (E o N)o hfcl} (E) b\ No Dans une étape 34, le module de traitement 110 vérifie si l'estimation de la valeur absolue de l'écart d'horloge VI est une première estimation. 5 S'il s'agit d'une première estimation, le module de traitement 110 considère que l'écart d'horloge est de signe positif lors d'une étape 35. Dans une étape 36, le module de traitement 110 sauvegarde la valeur de l'estimation de l'écart d'horloge ê dans une variable Ep'. Dans une étape 41, le module de traitement 110 corrige son horloge interne en 10 utilisant la valeur de l'estimation de l'écart d'horloge ê. S'il ne s'agit pas d'une première estimation, lors d'une étape 37, le module de traitement 110 compare l'estimation de l'écart d'horloge ê avec la variable sauvegardée dans la variable Ep' qui correspond à l'estimation précédente de l'écart d'horloge ê. Si ê > a lEp' I, avec a désignant un coefficient 1, i.e. si l'écart d'horloge augmente 15 significativement, le module de traitement considère que le signe de l'écart d'horloge a mal été déterminé lors de l'estimation précédente de l'écart d'horloge. Dans ce cas, lors d'une étape 38, le module de traitement 110 considère que le signe de l'estimation de l'écart d'horloge ê est l'inverse du signe de l'estimation précédente de l'écart Ep'. Sinon, lors d'une étape 39, le module de traitement 110 considère que le signe de 20 l'estimation de l'écart d'horloge ê est le même que le signe de l'estimation précédente de l'écart Ep'. Suite aux étapes 38 et 39, le module de traitement retourne à l'étape 36. Suite à l'étape 40, le module de traitement 110 se met en attente pendant une durée prédéfinie avant de procéder à une nouvelle estimation de l'écart d'horloge ê. 25 REVENDICATIONS1) Procédé d'estimation d'un écart d'horloge entre un émetteur et un récepteur durant une transmission sur un canal de communication d'un deuxième signal de communication résultant d'une application à un premier signal de communication d'une technique d'étalement de spectre utilisant un signal, dit signal d'étalement, caractérisé en ce que le procédé comprend suite à la réception du second signal par le récepteur : obtenir (311) au moins un signal, dit signal d'étalement orthogonal, orthogonal au signal d'étalement; pour chaque signal d'étalement orthogonal obtenu, appliquer (312) ledit signal d'étalement orthogonal obtenu au second signal reçu par le récepteur pour obtenir un troisième signal ; obtenir (313) une estimation d'un rapport énergie par bit sur densité spectrale de puissance d'un bruit de communication ayant affecté le deuxième signal durant sa transmission à partir de chaque troisième signal ; et, estimer (33) une information représentative de l'écart d'horloge entre l'émetteur et le récepteur en fonction de l'estimation dudit rapport. 2) Procédé selon la revendication 1, caractérisé en ce que, lorsqu'un seul signal d'étalement orthogonal est déterminé, ledit rapport est estimé à partir d'une fonction prédéterminée d'une moyenne arithmétique d'informations élevées au carré, lesdites informations étant représentatives du bruit de communication ayant affecté le deuxième signal durant sa transmission et ayant été obtenues à partir du troisième signal. 3) Procédé selon la revendication 1, caractérisé en ce que, lorsque une pluralité de signaux d'étalement orthogonaux est obtenue, pour chaque signal d'étalement orthogonal de la pluralité, une première moyenne arithmétique est calculée sur des informations élevées au carré, lesdites informations étant représentatives du bruit de communication ayant affecté le deuxième signal durant sa transmission et ayant été obtenues à partir du troisième signal obtenu pour ledit signal d'étalement orthogonal, ledit rapport étant estimé à partir d'une fonction prédéterminée d'une seconde moyenne arithmétique calculée à partir des premières moyennes arithmétiques calculées pour chaque signal d'étalement orthogonal. 3063405 19 4) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur comprend une estimation de la valeur absolue dudit écart d'horloge. 5 5) Procédé itératif d'estimation d'un écart d'horloge entre un émetteur et un récepteur, caractérisé en ce que le procédé selon l'une quelconque des revendications 1 à 4 est mis en oeuvre régulièrement sur le deuxième signal, l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur étant remise à jour à chaque itération. 10 6) Procédé selon la revendication 5, caractérisé en ce que l'information représentative de l'écart d'horloge entre l'émetteur et le récepteur comprend en outre un signe de l'écart d'horloge, ledit signe étant déterminé en comparant un écart d'horloge d'une itération courante avec un écart d'horloge d'une itération précédente, 15 le signe de l'écart d'horloge de l'itération courante étant l'inverse du signe de l'écart d'horloge de l'itération précédente lorsque l'écart d'horloge augmente et le signe de l'écart d'horloge de l'itération courante étant identique au signe de l'écart d'horloge de l'itération précédente lorsque l'écart d'horloge n'augmente pas. 20 7) Dispositif (110) comprenant des moyens pour mettre en oeuvre le procédé selon l'une quelconque des revendications 1 à 6. 8) Programme d'ordinateur, caractérisé en ce qu'il comprend des instructions pour mettre en oeuvre, par un dispositif (110), le procédé selon l'une quelconque des 25 revendications 1 à 6, lorsque ledit programme est exécuté par un processeur dudit dispositif (110). 9) Moyens de stockage, caractérisés en ce qu'ils stockent un programme d'ordinateur comprenant des instructions pour mettre en oeuvre, par un dispositif (110), 30 le procédé selon l'une quelconque des revendications 1 à 6, lorsque ledit programme est exécuté par un processeur dudit dispositif (110).