Brevet FR2737362A1 — publié le 19970131

PROCEDE DE SELECTION DES RETARDS DE PROPAGATION RETENUS POUR RECEVOIR DES MESSAGES TRANSMIS PAR RADIOCOMMUNICATION A ETALEMENT DE SPECTRE La présente invention concerne le domaine des communications radio à étalement de spectre par séquence directe (DSSS: "Direct Sequence Spread Spectrum Communications"). Elle s'applique notamment aux systèmes de radiotéléphonie à accès multiple à répartition par codes (AMRC ou, en anglais, CDMA pour "Code Division Multiple Access"). L'AMRC est un mode de communication numérique à spectre étalé dans lequel on crée différents canaux de transmission en utilisant, pour chaque canal, des séquences d'étalement modulant les bits d'information à transmettre. Les séquences d'étalement ont une cadence supérieure à celle des bits d'information pour réaliser l'étalement du spectre du signal radio. Leurs propriétés d'auto- et d'inter- corrélation sont adaptées pour permettre le multiplexage des différents canaux: ce sont en général des séquences pseudo- aléatoires, orthogonales ou quasi-orthogonales entre elles, prenant les valeurs -1 ou +1. Ces séquences sont souvent complexes, c'est-à-dire qu'elles comprennent deux séquences réelles modulant le même flux binaire et modulées par deux ondes porteuses en quadrature. L'utilisation de l'AMRC dans le domaine de la radiotéléphonie cellulaire est décrite dans le chapitre I de l'ouvrage "Mobile radio communications de Raymond Steele, Pentech Press, Londres, 1992, ainsi que dans l'article "On the system design aspects of code division multiple access (CDMA) applied to digital cellular and personal communications networks", de A. Salmasi et K.S. Gilhousen, Proc. of the 41st IEEE Vehicular Technology Conference, St Louis, MO, 19-22 mai 1991, et dans le brevet US 5 103 459. Pour la liaison descendante, les canaux de transmission multiplexés sont formés au niveau de la station de base de chaque cellule du réseau. Chaque station mobile située dans la cellule utilise une séquence d'étalement particulière pour retrouver, dans le signal radio global émis par la station de base, les bits d'information qui lui sont destinés. L'un des canaux AMRC de la liaison descendante est un canal pilote sur lequel la station de base émet une séquence de référence sans modulation de bits d'information. La réception sur le canal pilote permet aux stations mobiles de se synchroniser et d'évaluer la réponse du canal de propagation pour assurer une démodulation cohérente sur les autres canaux. Pour la liaison montante, chaque station mobile émet un signal radio formé en utilisant une séquence d'étalement particulière. Une corrélation avec cette séquence permet à la station de base de retrouver les bits d'information transmis par la station mobile, parmi les signaux radio reçus de différentes stations mobiles dans la cellule. Il n'y a pas de canal pilote dans le signal émis par une station mobile. Pour la plupart des systèmes de radiocommunication, les trajets de propagation multiples entre les stations émettrice et réceptrice constituent un problème délicat, en particulier parce qu'ils provoquent des évanouissements du signal (fading de Rayleigh). Un avantage important de l'AMRC est de réduire les inconvénients liés aux trajets multiples, grâce à l'étalement du spectre du signal. Avec 1'AMRC, on peut même profiter de l'existence des trajets multiples pour améliorer les performances de réception, par des techniques de diversité spatiale. On utilise alors un récepteur en râteau (rake receiver), dont chaque bras assure la réception des mêmes signaux suivant un trajet de propagation sélectionné. Chaque trajet est identifié par un retard qu'on applique à la séquence d'étalement pour déterminer une corrélation entre le signal reçu et la séquence retardée. Les corrélations obtenues dans différents bras du récepteur peuvent ensuite être combinées pour restituer les bits d'information transmis. Cette technique de diversité spatiale est expliquée dans le brevet US 5109390. Dans le système décrit dans ce document, l'un des bras du récepteur en râteau sert à balayer le domaine temporel pour rechercher de nouveaux trajets de propagation. Cette recherche a pour but de sélectionner pour la réception les trajets/retards pour lesquels l'énergie du signal reçu après corrélation est la plus importante. Les retards sont testés séquentiellement dans une fenêtre temporelle prédéfinie, en déterminant leur énergie associée. Chaque retard testé correspond à un décalage d'un nombre entier ou fractionnaire d'échantillons ou "chips", de la séquence d'étalement qu'on applique à celle-ci pour le calcul de la corrélation. La fenêtre temporelle des retards testés a typiquement une durée de l'ordre de 100ps, ce qui représente une propagation sur environ 30 km. Une fois qu'un retard a été sélectionné et affecté à un bras du récepteur en râteau, le trajet correspondant est poursuivi au moyen d'une boucle de poursuite, du type boucle de Costas, incluse dans ce bras. Un exemple d'une telle boucle est décrit dans l'article "Theory of Spread- Spectrum Communications- A Tutorial" de R.L. Pickholtz et al., IEEE Trans. on Communications, Vol. COM-30, N 5, mai 1982. Dans l'histogramme des énergies évaluées en fonction des retards testés dans la fenêtre temporelle, chaque trajet de propagation apparaît comme une impulsion dont la largeur temporelle est proportionnelle à l'inverse de la largeur de la bande passante de transmission. Pour améliorer la résolution dans l'estimation des retards, le signal reçu en bande de base est suréchantillonné, c'est-à-dire qu'il est numérisé à une cadence fs = 1/Ts multiple de la cadence fc = 1/Tc des séquences d'étalement: fs = M.fc Le lobe principal de la réponse impulsionnelle du système de transmission global, en supposant un trajet de propagation unique, a alors une largeur de plusieurs échantillons. Lorsqu'un premier retard est sélectionné à partir du lobe principal, il y a une probabilité assez forte de sélectionner aussi un ou plusieurs retards supplémentaires parmi les M retards précédant et les M retards suivant ce premier retard dans la fenêtre temporelle. Ces retards supplémentaires n'apportent un gain en diversité spatiale que si les trajets multiples de propagation ont une distribution concentrée. Sinon, ces retards supplémentaires contribuent certes à augmenter l'énergie reçue, mais ils n'apportent pas de gain en diversité. Un but de la présente invention est de proposer un procédé de sélection des retards pour un récepteur en râteau qui permette de séparer efficacement les contributions des trajets de propagation et celles du bruit de manière à exploiter pleinement la diversité spatiale/temporelle sans augmenter le bruit dans la sortie démodulée. L'invention propose ainsi un procédé de sélection de retards de propagation multiples pour la réception de messages transmis dans un système de radiocommunication à étalement de spectre, chaque retard sélectionné étant appliqué à au moins une séquence d'étalement pour déterminer une corrélation entre un signal reçu et ladite séquence d'étalement dans un bras de réception d'un récepteur en râteau à K bras. On échantillonne le signal reçu à une cadence d'échantillonnage supérieure à la cadence des séquences d'étalement. On teste des retards consécutifs multiples de l'inverse de la cadence d'échantillonnage dans une fenêtre temporelle, en évaluant, pour chaque retard testé, une énergie de la corrélation entre le signal reçu et une séquence d'étalement de référence à laquelle ledit retard testé est appliqué, et on sélectionne les retards sur la base de la distribution des énergies évaluées pour les retards testés. Le procédé comprend les étapes /a/ à /d/ ci-après /a/ déterminer un seuil de sélection à partir de la plus grande des énergies évaluées pour les retards testés dans la fenêtre temporelle; /b/ construire une première liste de retards et une seconde liste de retards, la première liste comportant des retards de la fenêtre pour lesquels l'histogramme des énergies évaluées en fonction des retards présente un maximum local supérieur au seuil de sélection, et la seconde liste comportant des retards qui sont chacun voisins d'au moins un retard de la première liste; /c/ si le nombre de retards de la première liste est égal ou supérieur au nombre K de bras du récepteur en râteau, sélectionner les K retards de la première liste pour lesquels les énergies évaluées sont les plus grandes; /d/ si le nombre de retards de la première liste est inférieur au nombre de bras du récepteur en râteau, sélectionner les retards de la première liste et au moins un retard de la seconde liste pour lequel l'énergie évaluée est supérieure au seuil de sélection. De cette façon, le processus de sélection favorise les retards de la première liste, qui correspondent en général à des trajets de propagation distincts. La seconde liste n'est finalement explorée, pour augmenter la puissance reçue sans nécessairement contribuer à la diversité de réception, que si le nombre de retards sélectionnés dans la première liste ne suffit pas aux K bras du récepteur en râteau. On peut ainsi mieux exploiter la diversité des trajets multiples, particulièrement lorsque le nombre de bras du récepteur en râteau est limité en raison de restrictions sur la complexité matérielle et/ou sur la consommation énergétique du récepteur. D'autres particularités et avantages de l'invention apparaîtront à la lecture de la description ci-après d'un exemple de réalisation non limitatif, en référence aux dessins annexés dans lesquels: - la figure 1 est un schéma d'un récepteur AMRC mettant en oeuvre l'invention; - la figure 2 est un schéma d'une unité d'estimation du canal de propagation du récepteur de la figure 1; - les figures 3 à 10 sont des histogrammes illustrant un procédé de sélection de retards selon l'invention; et - la figure 11 est un organigramme illustrant une façon d'exécuter le procédé. Dans un système de radiotéléphonie cellulaire, des stations de base sont réparties sur le territoire à couvrir. Dans la zone, ou cellule, couverte par une station de base, plusieurs stations mobiles peuvent communiquer simultanément avec la station de base. Le procédé selon l'invention est décrit ci-après dans son application aux liaisons descendantes, c'est-à-dire à la transmission des signaux depuis une station de base vers des stations mobiles, pour lesquelles on utilise un accès multiple à répartition par codes (AMRC). Chaque station de base établit des canaux de transmission AMRC définis chacun par une séquence d' étalement respective CP, Cl,..., Cx. Les séquences d'étalement sont par exemple complexes, avec des parties réelle et imaginaire prenant chacune les valeurs +1 ou -1. Ces séquences sont périodiques, pseudo-aléatoires et pratiquement décorrélées. Elles ont une cadence fc = 1/Tc supérieure à celle fb = 1/Tb des bits d'information à transmettre, par exemple une cadence de 5 MHz. Sur chaque canal AMRC de communication vocale, les bits d'information à transmettre, qui représentent des signaux vocaux ou de données préalablement codés par des techniques classiques, sont multipliés par la séquence d'étalement C1,..., Cx produite par un générateur approprié. Chacune des séquences Cl,..., Cx est caractéristique de l'utilisateur (station mobile) avec lequel la communication s'effectue sur le canal correspondant. L'un des canaux AMRC formés par la station de base est un canal pilote, sur lequel est émise une séquence d'étalement de référence CP, sans modulation de bits d'information. Les signaux à spectre étalé en bande de base formés dans les différents canaux AMRC sont combinés puis modulés sur une fréquence porteuse. La combinaison consiste en une sommation, éventuellement pondérée, des signaux en bande de base. La modulation peut être une modulation par déplacement de phase à deux ou à quatre états (MDP-2 ou MDP-4). La fréquence porteuse est par exemple de 2,4 GHz. Le signal radio obtenu après modulation est émis par la station de base à destination des stations mobiles situées dans la cellule. Un récepteur d'une station mobile est représenté sur la figure 1. Le signal radio recu sur l'antenne 6 est d'abord traité par l'étage radio 7 pour fournir un signal complexe en bande de base. Ce signal complexe est soumis à un filtrage passe-bas 8 puis numérisé par le convertisseur analogique- numérique 9 à une cadence d'échantillonnage fs supérieure à la cadence fc des séquences d'étalement: fs = M'fc, avec M > 1, M étant généralement un entier, typiquement M = 2 ou 4. Dans le signal de sortie r(t) du convertisseur 9 sont incluses les composantes AMRC émises par la station de base sur différents canaux AMRC. La station mobile peut retrouver les bits d'informations qui lui sont destinés en calculant les corrélations entre le signal r est la séquence d'étalement Cx qui lui est propre. Le récepteur de la station mobile comprend un récepteur en râteau 10. Dans l'exemple représenté, le récepteur en râteau 10 comporte K = 2 bras de réception de données 1, 2. Chaque bras de réception 1, 2 assure la réception du signal suivant un trajet de propagation identifié par un retard d(l), d(2). Les retards d(l), d(2) sont affectés aux bras de réception 1, 2 par une unité 40 d'estimation du canal de propagation. Le récepteur en râteau 10 comporte en outre un générateur pseudo- aléatoire 4 qui délivre aux bras 1, 2 la séquence d'étalement Cx propre au récepteur, et une unité de combinaison 5 qui reçoit les signaux B1, B2 délivrés par les bras de réception pour en déduire une estimation pondérée B des bits d'information reçus. Une base de temps 30 fournit aux éléments du récepteur différents signaux d'horloge obtenus par division de fréquence à partir d'un oscillateur à cristal, à savoir: - CKS, horloge d'échantillonnage de période Ts fournie au convertisseur analogique-numérique 9; - CKc, horloge chip donnant la cadence des séquences d'étalement, de période Tc = M.TS; - CKb, horloge bit donnant la cadence des bits d'information, de période Tb = Mb.Ts; CKi, horloge fournie à l'unité 40 d'estimation du canal de propagation, de période Ti = Mi.Ts, fixant la durée d'intégration Ti utilisée pour estimer la réponse impulsionnelle complexe du canal de propagation, et la période de mise à jour des amplitudes complexes A(k)ej propagation d(k) fournis aux bras du récepteur en râteau. Chaque bras 1, 2 du récepteur en râteau comprend un corrélateur formé, dans l'exemple représenté, par une mémoire tampon 11, 21, un multiplieur complexe 12, 22 et un sommateur-accumulateur 13, 23. La mémoire tampon 11, 21 est dimensionnée pour recevoir P échantillons successifs du signal reçu r(t), P.Ts étant la plage de variation possible des retards de propagation, définie en fonction de la taille maximale des cellules du réseau. Par exemple, pour une fréquence d'échantillonnage fs de 10 MHz, on peut prendre typiquement P = 1024, soit un retard maximal de l'ordre de is correspondant à un trajet de propagation de l'ordre de 30 km. En écriture, la mémoire 11, 21 se comporte comme un registre à décalage. Une lecture est effectuée à chaque cycle de l'horloge CKS, à une adresse dépendant de la valeur de retard d(k) (exprimée en unités de la période Ts). A chaque période d'échantillonnage, la mémoire 11, 21 produit la valeur du signal reçu retardé de P-d(k) durées d'échantillon Ts. Cette valeur est multipliée par l'échantillon de la séquence d'étalement Cx délivré par le générateur 4. La position initiale du générateur 4 est calée en fonction de données de synchronisation reçues de la station de base sur un canal de synchronisation dual du canal pilote (voir par exemple EP-A-0 630 120), et en prenant en compte un retard constant de P/M chips, de sorte que le retard entre la séquence d'étalement et le signal avec lequel elle est corrélée dans le bras k est bien de (P/M)Tc - [P- d(k)]Ts = d(k).Ts. La sortie du multiplieur 12, 22 est intégrée par le sommateuraccumulateur 13, 23 à la cadence d'échantillonnage fs La durée d'intégration correspond à la durée Tb d'un bit d'information, fournie par l'horloge bit CKb. Un multiplieur 14, 24 multiplie la sortie du sommateur 13, 23 par le complexe conjugué, fourni par l'unité 40, de l'amplitude complexe A(k)ejl(k) de la réponse du canal de propagation selon le trajet correspondant au retard d(k). A chaque temps bit, la sortie du multiplieur 14, 24 représente une estimation B1, B2 de la valeur de bit reçue selon le trajet considéré. Ces estimations sont combinées par l'unité 5 pour restituer l'estimation globale B. La figure 2 montre une unité 40 d'estimation du canal de propagation, utilisable dans le récepteur de la figure 1. Cette unité 40 comprend un corrélateur à glissement 50 ("sliding correlator") qui détermine les corrélations entre le signal reçu r(t) et la séquence de référence CP affectée d'un retard variant dans une fenêtre temporelle de recherche de N retards incluse dans la plage de P retards. Dans l'exemple représenté, le corrélateur à glissement 50 comporte une mémoire tampon 51 identique aux mémoires 11, 21 des bras du récepteur en râteau 10, et un filtre adapté 53 d'une longueur de Mi échantillons. Un compteur 43 sur log2N bits (on suppose ici que N est une puissance de 2, par exemple N = 128) délivre l'indice n du retard testé, incrémenté d'une unité à chaque période de l'horloge CKs(Oen obtenir la somme d = no + n représentant le retard testé. L'entier no est fourni à la mémoire 51 pour qu'elle délivre l'échantillon du signal reçu r(t) avec un retard de P-no échantillons. La séquence de référence CP est fournie au filtre adapté 53 par le générateur 45 à la cadence chip donnée par l'horloge CKc, avec un retard fixe de P/M chips. A chaque débordement du compteur 43, celui- ci délivre au filtre adapté 53 une impulsion CKS/N pour mettre à jour les coefficients du filitre adapté 53. Ainsi, lors du i-ième cycle du signal CKS/N de période NTS, c'est-à-dire entre les instants t=iNTs et t=(i+l)NTs, les Mi coefficients du filtre adapté 53 sont les échantillons de la séquence de référence CP aux instants d'échantillonnage (iN-P-Mi+l)Ts à (iN-P)TS. A l'instant t=(iN+n)Ts, le filtre adapté délivre la corrélation: Mi -1 _ r [ t-(P-n0+j) Ts)] *CP[ (iN-P-j) Ts] j=O0 Mi-1 = r ( t'-jT) *CP( t/'-jT -dT) j=0 avec t'=T-(P-no)Ts. Le retard d=no+n est donc bien testé lorsque le compteur 43 affiche la valeur n, en calculant une corrélation entre le signal d'entrée r et la séquence pilote CP sur une durée d'intégration Ti = Mi.Ts. La corrélation délivrée par le corrélateur à glissement 50 lorsque le compteur 43 affiche la valeur n est une estimation de l'amplitude complexe A(n)ej donnant la durée d'intégration Ti. La sortie du multiplieur 53 est en outre fournie à un élément 47 qui calcule son module carré et le fournit à une unité 48 de calcul de valeurs moyennes. A chaque période de l'horloge CKS, l'unité 48 reçoit la valeur relative n du retard testé et une estimation w(n) de l'énergie de la corrélation entre le signal reçu et la séquence de référence affectée de ce retard. Cette estimation w(n) est le module carré fourni par l'élément 47, la durée d'estimation de l'énergie étant la période Ti de l'horloge CKi. Pour chaque valeur de n (0 s n évaluation de l'énergie revue pour le retard relatif n. obtenue en moyennant Ma/Mi valeurs successives de l'estimation W(n). C'est sur la base des énergies évaluées W(n) que l'unité 49 sélectionne les retards d(k). Une fois que ces retards sont sélectionnés, l'unité 49 les fournit aux bras du récepteur en râteau 10 et à l'unité 46 de traitement de la RIC. Les retards de propagation d(k) sont en général stables sur des durées plus longues que les amplitudes complexes A(k)ej en compte pour chaque énergie évaluée W(n). Les figures 3 à 9 illustrent une procédure de sélection des retards à partir des énergies W(n) évaluées lors d'une période de mise à jour de durée Ta dans un cas particulier o K = 4, M = 2, N = 128. La figure 3 montre un exemple d'histogramme de W(n) fourni à l'unité 49 de sélection de trajets par l'unité 48 de calcul de moyennes. Un seuil de sélection TH est d'abord défini en retirant une marge de la plus grande valeur de W(n), exprimée en décibels. Le retard n = 25 pour lequel l'énergie W(n) est la plus grande est inclus dans une première liste Ll initialement vide. M = 2 retards voisins (n = 24 et n = 26) sont inclus dans une seconde liste L2 initialement vide. Un trajet de propagation se manifeste par un pic Xl- X5 dans l'histogramme des énergies, dont la largeur correspond à M + 1 échantillons, ou à M échantillons dans le cas particulier o le retard est très proche d'un multiple entier et impair de la moitié de la période d'échantillonnage Ts (cas des pics X3 et X4 sur la figure 3). Dans le cas particulier o M = 2, un pic de largeur M + 1 = 3 présente en général un maximum pour l'échantillon central n (coincidence entre les deux valeurs de chip sur toute la durée de l'échantillon) et une valeur encore relativement élevée pour les échantillons n - 1 et n + 1 (coincidence entre les deux valeurs de chip sur une partie seulement de la durée de l'échantillon). Il convient donc de prendre les retards n - 1 et n + 1 comme les deux "voisins" à inclure dans la liste L2. Si le pic a une largeur M = 2 seulement, l'un des deux retards inclus dans la liste L2 représente une corrélation faiblement énergétique, mais ce retard aura finalement une faible probabilité d'être sélectionné. Si M > 2, les valeurs de chip coincident sur M - 1 Z 2 durées d'échantillons consécutifs, de sorte que le maximum de l'énergie ne correspond pas nécessairement au centre du pic. Une possibilité est de retenir comme "voisins" les M/2 retards précédant et les M/2 retards suivant le retard pour lequel l'énergie est maximale. On peut également analyser la forme du pic pour éventuellement prendre davantage de voisins d'un côté du maximum que de l'autre pour tenir compte des cas o le maximum du pic ne serait pas central. Dans tous les cas, un retard "voisin" d'un autre retard est distant d'au plus M-1 échantillons de cet autre retard. Dans une variante de réalisation, on n'inclut un retard voisin dans la seconde liste L2 qu'à condition que l'énergie évaluée pour ce retard soit supérieure au seuil TH. Ceci ne change pas les retards qui seront finalement sélectionnés puisque l'énergie évaluée pour un retard sélectionné doit toujours être supérieure au seuil. Une fois que le premier retard a été inclus dans la liste Ll et ses voisins dans la liste L2, ces M + 1 retards sont éliminés de l'histogramme n,W(n) (figure 4). On recherche alors le retard restant pour lequel l'énergie W(n) est maximum (n = 47 dans l'exemple représenté). Si ce maximum est supérieur au seuil TH, ce retard est inclus dans la liste Ll et ses M = 2 voisins (n = 46 et n = 48) sont inclus dans la liste L2. Le processus de construction des listes mentionné au paragraphe précédent est réitéré jusqu'à ce que l'histogramme n. W(n) ne présente plus d'énergie supérieure au seuil TH (figure 7). Les figures 5,6 et 7 montrent l'histogramme n. W(n) après 2, 4 et 5 itérations, respectivement. La figure 8 montre le contenu final de la seconde liste L2, les 10 retards étant rangés en ordre croissant. La sélection des retards est d'abord effectuée dans la première liste Ll. Si cette liste Ll ne contient pas suffisamment de retards pour les K bras du récepteur en rateau, la seconde liste L2 est examinée pour y sélectionner d'autres retards. Dans le cas de la figure 9, le récepteur en râteau comporte K = 4 bras. La liste L1 contenant 5 retards, on sélectionne les 4 premiers n(1) à n(4), c'est-à-dire ceux pour lesquels les énergies W(n) sont les plus grandes. Dans le cas de la figure 10, le récepteur en râteau comporte K = 8 bras. Les 5 retards n(l) à n(5) de la liste Ll sont sélectionnés, de même que les 3 retards n(6) à n(8) de la liste L2 pour lesquels les énergies W(n) sont les plus grandes. Les retards absolus d(k) finalement affectés par l'unité 49 aux bras du récepteur sont les retards relatifs n(k)augmentés de la position no du premier échantillon de la fenêtre de recherche. On note que, dans le cas de la figure 10, les retards sélectionnés n(1) à n(8) correspondent aux 8 retards pour lesquels les énergies sont les plus grandes au départ (figure 3). Dans ce cas particulier le processus de sélection conduit au même résultat qu'un processus classique basé seulement sur les niveaux énergétiques. Mais dans le cas de la figure 9, o le nombre de bras est moindre, on sélectionne 4 retards correspondant à 4 trajets distincts, tandis que le processus classique conduirait à sélectionner les retards n = 24, 25, 26 et 47, incluant 3 retards voisins probablement dûs au même trajet de propagation. L'invention permet donc de mieux tirer parti de la diversité en réception lorsque le nombre de bras est limité. La figure 11 montre un exemple d'organigramme de la procédure de sélection des retards à partir de l'histogramme n,W(n) obtenu sur une période Ta. A la première étape 60, l'unité 49 calcule le seuil TH en retirant la marge A de la plus grande des énergies W(n) sur la fenêtre, et initialise à 0 l'entier i. Une boucle 62, 64, 66 est alors exécutée plusieurs fois pour construire les listes Ll et L2 initialement vides. A l'étape 62, ondétermine le retard j pour lequel l'énergie évaluée est la plus grande dans la fenêtre de recherche à l'exclusion des retards des listes Ll et L2. Cette plus grande énergie est comparée au seuil TH à l'étape 64. Si la comparaison 64 montre que W(j) > TH, on exécute l'étape 66 avant de revenir à l'étape 62. A l'étape 66, l'entier i est incrémenté d'une unité et le retard n(i) est pris égal à l'entier j qui a précédemment maximisé l'énergie. Le retard n(i) est inclus dans la liste Li, et les retards voisins (n(i)-l et n(i)+l dans le cas o M = 2) sont inclus dans la liste L2. Lorsque la comparaison 64 montre que W(j) 5 TH, on sort de la boucle 62, 64, 66, les listes Li et L2 étant alors complètes. Le nombre i de retards de la première liste est comparé, en 70, au nombre K de bras du récepteur en râteau. Si i 2 K, les K premiers retards de la première liste Ll sont sélectionnés à l'étape 72, les retards absolus d(k) étant obtenus pour 1 e k eK par d(k) = no + n(k). Si la première comparaison 70 montre que i TH, on exécute l'étape 78 avant de revenir à l'étape de comparaison 70. A l'étape 78, l'entier i est incrémenté d'une unité, et le retard n(i) est pris égal à l'entier j qui a précédemment maximisé l'énergie. Le retard n(i) est inclus dans la liste L3. Lorsque la comparaison 70 montre que i = K, K retards sont sélectionnés à l'étape 72, à savoir tous les retards de la liste Ll et un ou plusieurs retards de la liste L2 (ceux pour lesquels les énergies sont les plus grandes). Si une occurrence de la comparaison 76 montre que W(j) 5 TH, la seconde liste L2 ne contient pas suffisamment de retards supérieurs au seuil TH pour fournir un total de K retards avec ceux de la première liste. Dans ce cas, les retards sont sélectionnés à l'étape 80, en nombre égal à la valeur courante de l'entier i. Une fois que les retards ont été sélectionnés à l'étape 72 ou 80, l'unité 49 calcule à l'étape 82 deux quantités Es et Io, puis détermine, à l'étape 84, les paramètres A et no pour la période Ta suivante. La nouvelle position no de la fenêtre de recherche est déterminée à partir de la distribution des retards qui viennent d'être sélectionnés. Une possibilité simple est de prendre no égal au plus petit des retards d(k) diminué d'une valeur constante. Des algorithmes plus sophistiqués peuvent néanmoins être appliqués pour positionner la fenêtre. La marge A est déterminée sur la base d' estimations du rapport signal-à-interféreurs (SIR: 'Signal-to- Interferers Ratio"). L'énergie globale reçue Es a été déterminée à l'étape 82 en faisant la somme des énergies W[n(k)] évaluées pour les retards n(k) retenus à l'étape 72 ou 80, tandis que la densité de puissance Io des interféreurs a été déterminée en faisant la somme des énergies W(n) évaluées pour les retards de la fenêtre de recherche non inclus dans les listes Ll et L2. Le quotient Es/Io représente une estimation du rapport signal-à-interféreurs à chaque période Ta. Ces estimations sont soumises à un filtrage passe-bas, et la marge A est prise égale à l'estimation Es/Io filtrée multipliée par un coefficient prédéterminé. A titre indicatif, la demanderesse a testé le procédé dans un cas particulier o la cadence d'échantillonnage fs était de 10 MHz, pour une cadence chip de 5 Mchips/s (M = 2). La fenêtre de recherche couvrait N = 128 échantillons soit 12,8 ps. Dans le corrélateur à glissement, la durée d'intégration Ti correspondait à Mi = 2560 échantillons ou 1280 chips. La période de mise à jour des amplitudes complexes était égale à la durée d'intégration Ti, soit 256gs. La période Ta de mise à jour des retards était de ,24ms (Ma = 102400). Ma/Mi = 40 estimations w(n) étaient prises en compte dans l'évaluation des énergies W(n). Avec un récepteur à K = 4 bras, la demanderesse a observé une meilleure exploitation des retards correspondant à des trajets distincts du plus énergétique, contrairement aux récepteurs classiques qui ont trop tendance à sélectionner simultanément plusieurs retards voisins correspondant au trajet le plus énergétique. L'invention a été décrite ci-dessus dans un exemple d'application aux liaisons descendantes d'un système radiotéléphonique. Elle peut bien entendu s'appliquer de façon semblable aux liaisons montantes lorsque le récepteur de la station de base comporte un récepteur à diversité (voir US 5109390). Dans ce cas, la séquence de référence utilisée pour le calcul des énergies de réception n'est plus une séquence émise sur un canal pilote, mais la séquence d'étalement servant aux communications. REVENDICATIONS 1. Procédé de sélection de retards de propagation multiples pour la réception de messages transmis dans un système de radiocommunication à étalement de spectre, chaque retard sélectionné étant appliqué à au moins une séquence d'étalement pour déterminer une corrélation entre un signal reçu et ladite séquence d'étalement dans un bras de réception d'un récepteur en râteau à K bras, procédé dans lequel on échantillonne le signal reçu à une cadence d'échantillonnage supérieure à la cadence des séquences d'étalement, on teste des retards consécutifs multiples de l'inverse de la cadence d'échantillonnage dans une fenêtre temporelle, en évaluant, pour chaque retard testé, une énergie de la corrélation entre le signal reçu et une séquence d'étalement de référence à laquelle ledit retard testé est appliqué, et on sélectionne les retards sur la base de la distribution des énergies évaluées pour les retards testés, caractérisé en ce que la sélection des retards comprend les étapes suivantes: /a/ déterminer un seuil de sélection à partir de la plus grande des énergies évaluées pour les retards testés dans la fenêtre temporelle; /b/ construire une première liste de retards et une seconde liste de retards, la première liste comportant des retards de la fenêtre pour lesquels l'histogramme des énergies évaluées en fonction des retards présente un maximum local supérieur au seuil de sélection, et la seconde liste comportant des retards qui sont chacun voisins d'au moins un retard de la première liste; /c/ si le nombre de retards de la première liste est égal ou supérieur au nombre K de bras du récepteur en râteau, sélectionner les K retards de la première liste pour lesquels les énergies évaluées sont les plus grandes; /d/ si le nombre de retards de la première liste est inférieur au nombre de bras du récepteur en râteau, sélectionner les retards de la première liste et au moins un retard de la seconde liste pour lequel l'énergie évaluée est supérieure au seuil de sélection. 2. Procédé selon la revendication 1, caractérisé en ce que la cadence d'échantillonnage est égale à M fois la cadence des séquences d'étalement, M étant un entier plus grand que 1, et en ce que chaque retard de la seconde liste est distant d'au plus M-1 échantillons d'un retard de la première liste. 3. Procédé selon la revendicaiton 2, caractérisé en ce que l'étape /b/ de construction des listes comprend, à partir de listes initialement vides, plusieurs itérations de la succession des sous-étapes suivantes: /bl/ déterminer le retard testé, parmi les retards de la fenêtre temporelle à l'exclusion des retards des première et seconde listes, pour lequel l'énergie évaluée est la plus grande; /b2/ comparer l'énergie évaluée pour le retard déterminé à l'étape /bl/ au seuil de sélection; /b3/ si l'étape /b2/ montre une énergie supérieure au seuil de sélection, inclure dans la première liste le retard déterminé à l'étape /bl/, et inclure dans la seconde liste au moins chaque retard, parmi M retards voisins du retard déterminé à l'étape /bl/, pour lequel l'énergie évaluée est supérieure au seuil de sélection. 4. Procédé selon la revendication 3, caractérisé en ce que, M étant un entier pair, lesdits M retards voisins d'un retard déterminé à l'étape /bl/, parmi lesquels au moins un retard est inclus dans la seconde liste à l'étape /b3/ suivante, sont les M/2 retards précédant et les M/2 retards suivant ledit retard déterminé à l'étape /bl/. 5. Procédé selon l'énergie des revendications 1 à 4, caractérisé en ce que les listes sont construites à l'étape /b/ de façon qu'aucune des énergies évaluées pour les retards de la fenêtre temporelle à l'exclusion des retards des première et seconde listes ne soit supérieure au seuil de sélection, le procédé comprenant en outre une étape d'estimation d'un rapport signal-à-interféreurs basée sur le quotient entre d'une part la somme des énergies évaluées pour les retards sélectionnés à l'étape /c/ ou à l'étape /d/ et d'autre part la somme des énergies évaluées pour les retards de la fenêtre temporelle à l'exclusion des retards des première et seconde listes. 6. Procédé selon la revendication 5, caractérisé en ce que les retards sont sélectionnés par périodes de mise à jour successives de Ma échantillons, et en ce que le seuil de sélection déterminé à l'étape /a/ pour une période de mise à jour est calculé en fonction de la plus grande des énergies évaluées pour les retards testés dans la fenêtre temporelle lors de ladite période de mise à jour et du rapport signal-à- interféreurs précédemment estimé. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que chaque bras du récepteur en râteau est agencé pour effectuer une démodulation cohérente pour le retard de propagation (d(k)) qui lui est affecté, sur la base d'une estimation d'une amplitude complexe de réception (A(k)eji(k)) associée audit retard, et en ce que les estimations des amplitudes complexes (A(k)eje(k)) sont mises à jour à une fréquence plus élevée que la fréquence de mise à jour des retards.