Brevet FR2913160A1 — publié le 20080829

DÉCODEUR À MAXIMUM DE VRAISEMBLANCE POUR SYSTÈME MULTI- SOURCE À MODULATION DE POSITION D'IMPULSION DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne le domaine du décodage à maximum de vraisemblance, plus particulièrement pour un système multi-source utilisant une modulation de position d'impulsion. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Les récepteurs utilisant un critère de maximum de vraisemblance, dénommés aussi récepteurs ML (Maximum Likelihood) sont bien connus dans le domaine des télécommunications pour être optimaux lorsque le canal de transmission est gaussien. On trouvera par exemple une description de ces récepteurs dans l'ouvrage de J.G. Proakis intitulé Digital communications , [dème édition, pages 242-247. Les récepteurs à maximum de vraisemblance ont été notamment envisagés dans le domaine des télécommunications mobiles. Afin d'éliminer l'interférence multi-accès ou MAI (Mufti Access Interference), on peut recourir à un récepteur ML capable de décoder simultanément les symboles émis par les différents utilisateurs sur le canal de transmission (récepteur ML multi-utilisateur). On peut montrer que l'estimation des symboles transmis par ces utilisateurs selon le critère du maximum de vraisemblance revient à rechercher parmi les points d'un réseau celui qui est le plus proche d'un point représentatif du signal reçu dans un espace de dimension MK où M est la dimensionnalité de la modulation utilisée par les K utilisateurs. Le réseau de point est généré par les constellations de modulation des différents utilisateurs. Cette méthode s'avère rapidement complexe pour K élevé aussi recourt-on classiquement à un décodage dit par sphère qui limite la recherche du plus proche voisin aux points du réseau appartenant à une boule de bruit centrée sur le point reçu. On trouvera une description du décodeur par sphère dans l'article d'E. Viterbo et al. intitulé A universal lattice code decoder for fading channels publié dans IEEE Transactions on Information Theory, vol. 45, pages 1639-1642, Juillet 1999. Le décodage par sphère a été appliqué à des systèmes utilisant des constellations de modulation de type PAM ou QAM. Plus récemment, le décodage par sphère a été proposé pour la réalisation de récepteurs de systèmes MIMO (Multiple Input Multiple Output). On entend par système MIMO un système de télécommunication dans lequel au moins un émetteur transmet des symboles d'information au moyen d'une pluralité d'antennes. Le récepteur peut posséder une seule antenne (on parle alors plus particulièrement de système MISO, acronyme de Multiple Input Single Output) ou une pluralité d'antennes, le terme MIMO étant choisi ici pour désigner indifféremment ces deux configurations. Dans le cas d'un système MIMO, le réseau de points est généré par les constellations de modulation utilisées pour transmettre les symboles par les différentes antennes. On trouvera un exemple de réalisation de décodage par sphère pour système MIMO dans l'article de M.O. Damen et al. intitulé On Maximum-Likelihood detection and the search for the closest lattice point publié dans IEEE Trans. on Information Theory, Vol. 49, N 10, Oct. 2003, pages 2389-2402. Dans ce décodeur par sphère, seules des modulations PAM et QAM sont envisagées. Dans la suite nous désignerons par le terme générique multi-source aussi bien une configuration multi-utilisateur qu'une configuration MIMO. On comprendra que dans le premier cas les sources représentent les flux de symboles des ou à destination des différents utilisateurs et, dans le second cas, les flux de symboles émis par les différentes antennes. Bien entendu ces deux cas peuvent être combinés lorsque les terminaux des utilisateurs sont de type multiantenne. Nous supposerons par ailleurs que les flux de symboles sont synchrones. Le récepteur ML à décodage par sphère sera décrit dans le cas de K utilisateurs, chaque utilisateur 25 kE {1,..,K} transmettant à l'aide de ik antennes à destination du récepteur, soit un nombre total de K P = Eik sources. Le signal reçu peut être exprimé sous k=1 la forme vectorielle : 30 x=Ha+n (1) où . x est un vecteur de variables de décision de dimension P' où P' est égal au produit du nombre d'antennes du récepteur avec le nombre de variables de décision observées par antenne de réception, par exemple le nombre de trajets résolus par antenne. a est un vecteur de taille PM, obtenu par la concaténation de vecteurs aw,W4,..,WP), chaque vecteur W0, pE {1,..,P} étant la représentation vectorielle du symbole d'information transmis par la pi ème source et étant de dimension M égale à la dimensionnalité de la modulation utilisée ; H est la matrice de taille P'xPM représentant le canal de transmission, elle décrit notamment les interférences entre utilisateurs et entre trajets issus des différentes antennes ; n est un vecteur de dimension P' dont les composantes sont des échantillons de bruit blanc centré additif gaussien affectant le signal reçu. Le récepteur à maximum de vraisemblance estime le vecteur â minimisant l'écart quadratique x ù Hat avec le vecteur reçu, soit : â=argmm xùHa2 (2) aE CP où CP est la constellation produit des constellations respectives des P sources. On peut montrer que l'expression (2) peut s'écrire de manière équivalente : â = arg mm z ù Ra r aE cP avec z = QTx et où Q et R sont respectivement une matrice unitaire de taille P'xP et une matrice triangulaire supérieure de taille PxP obtenue par décomposition QR de la matrice H, c'est-à-dire RTR=HTH. Pour des constellations PAM ou QAM, il est toujours possible, au moyen d'une opération linéaire élémentaire f, dont nous ferons par la suite abstraction, de se ramener au cas où les points de la constellation produit sont de éléments de ZMP où Z est l'ensemble des entiers relatifs et M la dimensionnalité des constellations de modulation . Les vecteurs Ra peuvent alors être représentés comme les points d'un réseau A de matrice génératrice R. Le décodage par sphère consiste à faire une recherche de plus proche voisin au sein d'une boule (dans un espace de dimension MP) centrée sur le point représentatif du signal reçu z. La recherche est effectuée parmi les points de A compris au sein de cette boule. La manière dont les candidats sont successivement sélectionnés (ou énumérés) est critique pour les performances de l'algorithme de décodage. On connaît (3) essentiellement deux techniques d'énumération, la première dite de Pohst, la seconde dite de Schnorr-Euchner. Nous décrirons tout d'abord un décodage par sphère utilisant une énumération de Pohst. On supposera pour les besoins d'illustration que la modulation est de type PAM, c'est-à-dire que l'espace considéré est de dimension P. Il s'applique également à une modulation QAM, la dimension de l'espace étant alors 2P. La recherche est effectuée, dimension par dimension, ou selon la terminologie usuelle, couche par couche, en sélectionnant dans chaque couche la coordonnée d'un point candidat de A. La contribution de la couche P à la distance quadratique (3) vaut simplement : (zp û rppap) ( 4) De même la contribution de la couche Pû1 à cette distance quadratique s'écrit : (Zp û rp_ ,p_ ap_ û rp_ ,pap (5) et, de manière générale, la contribution de la couche i s'exprime par: (zi ûri,i ai ûEi )2 (6) P avec Ei défini par Ei = Si l'on désigne par Ti, j=i+1 la contribution des couches successives i+1,...,P, celle-ci peut être obtenue par récurrence descendante : 5 T_1=T+ziùEiùriiai 2 (7) avec Tp =O. Avec les conventions précédentes la recherche dans la sphère de rayon quadratique d revient à déterminer 10 pour chacune des couches les valeurs ai de Z vérifiant . z ûE~ +T 20 où M' est une puissance de 2 et si la transformation linéaire précitée est (m')=2m'û1ûM' pour m'E {1,..,M'}, la recherche peut être réduite aux entiers impairs appartenant à l'intervalle [Ai,Bi] où : (9) Ai = max 1ùM', où [xi est le plus petit entier supérieur à x et Lx] est le plus grand entier inférieur à x. Cette réduction de l'intervalle de recherche correspond en fait à une intersection de l'intervalle tel que défini en (9), avec la projection de la constellation sur la dimension en question. Un intervalle de validité est ainsi défini pour chaque couche. Il donne l'excursion maximale admissible dans cette couche, compte tenu du choix déjà fait pour les valeurs ai dans les couches supérieures. La Fig. 1 représente le principe du décodeur par sphère dans le cas simple où P=2 et M'=4. La sphère est centrée sur le point z représenté par une croix. Les points de la constellation produit C2 sont représentés par des cercles pleins. L'intervalle de validité [A2,B2] de la couche 2 n'est autre que l'intersection de la projection de la sphère avec celle de la constellation, ici [A2,B2]= [-3,3] . L'intervalle de validité [A1,B1] est ensuite déterminé par le choix d'une valeur a2 dans [A2,B2]. On voit ainsi dans l'exemple illustré que pour a2=1, l'intervalle de validité est [A1, B1 ] _ [ù 3,1] . 8 ( r l'i,i / Bi = min M'ù1, zi ùEi+V dùTi l'i,i et (10) La Fig. 2 représente un organigramme de la méthode de décodage par sphère utilisant une énumération de Pohst. A l'étape 210, on procède à l'initialisation de l'indice de couche, des variables intermédiaires et du rayon de recherche de la sphère, soit : i=P ; Tp =0 ; Ep =0 ; d =D où D est le rayon quadratique initial de la sphère de recherche. D est choisi en fonction d'une estimation de la puissance de bruit. A l'étape 220, on teste si Ti >d, c'est-à-dire si l'intervalle se réduit à l'ensemble vide. Si c'est le cas, on teste en 223 si i=P et dans l'affirmative l'algorithme se termine en 225. Dans la négative, on incrémente i de 1 et l'on passe à l'étape 230. Si T L'algorithme exposé ci-dessus effectue une recherche en balayant les intervalles de validité couche par couche à partir de leurs bornes inférieures vers leurs bornes supérieures, le balayage le plus rapide étant effectué dans la couche la plus haute (i= P) et le plus lent dans la couche la plus basse (i =1) . A chaque sélection d'une nouvelle coordonnée de candidat dans une couche, les bornes des intervalles de validité des couches supérieures sont recalculées. Par ailleurs, dès que l'on améliore la distance au point reçu, on met à jour le rayon de la sphère, ce qui permet d'accélérer la recherche. Enfin, le vecteur de sortie de l'algorithme est le point de la constellation produit le plus proche du point reçu z. Le décodage Pohst est assez des intervalles inférieures. Le décodage Schnorr-Euchnerpar sphère utilisant l'énumération de lent du fait du balayage systématique de validité à partir de leurs bornes par sphère utilisant l'énumération de est sensiblement plus rapide que le précédent. Il de faire une constellation part du principe qu'il est plus efficace recherche en partant d'un point ê de la produit correspondant à une première estimation du signal transmis. Plus précisément, ce point est obtenu au moyen d'une égalisation de type ZF- DFE: on soustrait l'interférence due à une couche (un utilisateur) aux couches (inférieures) suivantes, la décision à chaque couche étant prise par le décodeur par sphère. Soit, en commençant par la couche P : êi = round(ei) avec ep =zplrPp p e i = (zi û i i eti = ~,JêJ pour i = Pû1,...,1 (11) J=i+1 où round(y) représente l'entier relatif impair (avec la forme précédente de transformation linéaire L) le plus proche de y. Chaque couche i contribue à la distance quadratique du point a du réseau au point z reçu pour di = ril(ai û 2 La recherche du plus proche voisin s'effectue dans chaque couche i en considérant successivement les coordonnées ê+26, ê+46, etc. où bi =sgn(ei ùêi) . On comprendra que l'on part ainsi de ei et que l'on progresse en zigzag de part et d'autre de cette estimation tout en s'assurant de rester dans l'intervalle de validité. La Fig. 3 représente un organigramme de la méthode de décodage par sphère utilisant une énumération de Schnorr-Euchner. A l'étape 310, on initialise l'indice de couche, les variables intermédiaires et le rayon de recherche de la sphère, soit : i=P; p=0 ;TP=0 ; d=D où D est le rayon quadratique initial de la sphère de recherche. A l'étape 320, on procède à l'égalisation ZF-DFE : êi_ round(ei) avec ei_ (ziù i)/il, et l'on détermine la valeur de coordonnée initiale al=êl ainsi que la direction de recherche initiale : 8i=sgn(eiùêi). A l'étape 330, on teste si le point est dans la Y voisin â=a en 365. Dans tous les cas on passe par 370 pour incrémenter i de 1. Enfin, à l'étape 380, on met à jour al par al =al+BZ avec 8z = -8z -2sgn(8Z) et l'on retourne à l'étape 330. On a représenté schématiquement, en Fig. 4A et 4B, le balayage des points selon l'énumération de Pohst et celle de Schnorr-Euchner, respectivement. On a supposé que la constellation était 4-PAM et que le nombre de sources P était de 3. Les différentes lignes correspondent aux projections selon les P=3 dimensions. On a représenté par des croix les composantes de z et par des cercles les points des constellations de modulation PAM. Dans les deux cas, les branches successivement parcourues par l'algorithme ont été indiquées par (1), (2) etc. Alors que le décodage par sphère a fait l'objet de nombreuses études pour les modulations PAM et QAM, il n'a été envisagé que récemment pour les modulations PPM (Pulse Position Modulation). Un décodage par sphère pour ce type de modulation a été proposé pour un système MIMO UWB de type impulsionnel dans l'article de C. Abou-Rjeily et al. intitulé MIMO UWB communications using modified Hermite pulses , publié dans Proc. of the 17th Annual IEEE Internal Symposium on Personal, Indoor and Mobile Communications (PIRMC'06). On rappelle qu'un système UWB impulsionnel est un système utilisant la transmission en bande de base de trames d'impulsions très brèves (de l'ordre de quelques picosecondes). Les symboles transmis appartiennent à un alphabet de modulation M -PPM, ou de manière plus générale à une modulation hybride AI-PPMAl'-PAM. Le cardinal de l'alphabet de modulation est MM'. Pour chacune des M positions temporelles, M' amplitudes de modulation sont possibles. Un symbole a de cet alphabet peut être représenté par une séquence am, m=O,..,Mû1 avec am =ô(mû )a où est une position de la modulation AI-PPM, î une amplitude de la modulation PAM et 8(.) la distribution de Dirac. Chaque symbole PPM ou PPM-PAM module les positions temporelles et, le cas échéant, les positions temporelles et les amplitudes des impulsions élémentaires d'une trame. On peut considérer un symbole AI-PPM (resp. At-PPM-Al'-PAM) comme un vecteur de dimension Al, dont une seule des composantes est égale à 1 (resp. î) et les autres sont nulles. A la différence d'une modulation QAM, les composantes d'un symbole PPM ou PPM-PAM ne sont donc pas indépendantes. Si l'on suppose que P sources transmettent des symboles PPM, l'espace des signaux à considérer est de dimension MP. Dans cet espace, chaque source est représentée par une couche, elle-même composée de M positions temporelles (ou sous-couches) d'un symbole PPM ou PPM-PAM. Le décodage par sphère décrit dans l'article précité propose de traiter de manière conjointe M couches consécutives en raison de leur interdépendance. La relation (8) est alors remplacée par une contrainte liant M sous-couches successives, soit pour la couche d'ordre p .30 PM z(p) ù r a m f (p,m),j j j=f(p,m) d-T M t m=1 2 (12) où . z(P) est un sous vecteur de z, représentatif du signal reçu, autrement dit le vecteur z est formé par concaténation de P vecteurs z(1) , .. , z(P) correspondant chacun à une source, z(p) étant la m1ème composante du symbole PPM ou PPM-PAM émis par cette source ; la fonction f est donnée par f(p,m)=(pû1)M+m ; Tp est la somme des contributions à la distance quadratique des couches p+1,...,P. Etant donné que, pour chaque source p, la transmission n'a lieu que pour une seule position c(p)E{1,..,M}, l'expression (12) peut être réécrite sous la forme : M 2 z mp) ùrf(P,m),f(p,c(p))af(P,c(P)) ùE'P,m d ùTp (13) m=1 PM avec E'pm= Lrf(pm) jaj , ou bien encore sous la forme plus j=f (p+1,1) compacte: Z(P) _a(P) R(P,P) ùE Z d ûT c(P) •,c(P) P (14) où le vecteur a, comme le vecteur z, est une concaténation de P vecteurs au),a(2),...,a(P), le vecteur a(P) étant de composantes amp) = af(pm) , m =1,..,M ; R(`'') est une sous-matrice M xM de R constituée des 5 éléments rmm, avec m=(iû1)M+1,..,iM et m'=(jûl)M+1,..,jM ; pour une matrice SZ donnée, SZk désigne conventionnellement la kieme colonne de cette matrice ; E.4) est la pierre colonne d'une matrice E et est définie par P E û a(p') R(P,P') c(p) •,c(Pt) p'=p+1,1 10 On comprendra de l'expression (14) que les M sous-couches correspondant à la source p sont traitées conjointement. En particulier, les bornes des intervalles de validité de chacune des composantes a(.''), 15 m =1,..,M sont déterminées en même temps. Ces bornes dépendent des candidats a(Pn),..a(P) déjà choisis pour les couches supérieures. La technique d'énumération utilisée dans ce décodage par sphère peut être considérée comme une 20 extension d'une énumération de Pohst. Les positions temporelles de chaque couche sont balayées une à une, à partir de leurs bornes inférieures, en prenant la précaution de ne sélectionner qu'une position par couche. La convergence du décodage est donc 25 relativement lente. D'autre part, l'application directe d'une technique d'énumération de type Schnorr-Euchner est loin d'être triviale puisque les sous-couches, c'est-à-dire les positions PPM d'une même couche ne sont pas indépendantes. Le but de la présente invention est de proposer un procédé de décodage par sphère pour une pluralité de sources PPM, qui converge plus rapidement que celui connu de l'état de la technique. EXPOSÉ DE L'INVENTION La présente invention est définie par un procédé de décodage par sphère pour récepteur à maximum de vraisemblance destiné à recevoir des symboles PPM d'une pluralité P de sources, chaque source émettant un flux de symboles PPM à M positions de modulation, P symboles PPM émis simultanément par les P sources étant représentés par un point de la constellation de modulation produit dans un espace des signaux transmis de dimension MP décomposé en P couches, chaque couche représentant les M positions de modulation possibles d'un symbole PPM émis par cette source, le signal reçu par ledit récepteur étant transformé en un point représentatif de ce signal, dit point reçu, dans l'espace des signaux transmis, ledit procédé déterminant le point de la constellation produit le plus proche du point reçu à l'intérieur d'une sphère de rayon quadratique donné. Selon ledit procédé, pour chaque couche de rang p . (a) on effectue une égalisation ZF-DFE du signal reçu dans ladite couche compte tenu des symboles PPM estimés dans les P û p couches précédentes, dites couches supérieures à la couche p; (b) on classe dans une liste les M symboles PPM de ladite couche en fonction des contributions qu'ils apporteraient à la distance quadratique au point reçu ; (c) on sélectionne le symbole PPM réalisant la plus faible contribution et l'on ajoute cette contribution à celles obtenues pour les couches précédentes pour obtenir une somme de contributions ; les étapes (a),(b),(c) étant répétées jusqu'à ce que l'on atteigne la couche la plus basse ; et l'on met à jour le rayon quadratique de la sphère et ledit point le plus proche si ladite somme de contributions est inférieure au rayon quadratique de la sphère. Avantageusement, si pour une couche donnée et un symbole PPM sélectionné dans cette couche, ladite somme des contributions excède le rayon quadratique de la sphère, on passe à la couche supérieure et on sélectionne dans cette couche le symbole suivant de la liste qui lui est associée. Si tous les symboles de ladite liste associée ont déjà fait l'objet d'une sélection, on passe à la couche encore supérieure pour y sélectionner dans cette dernière le symbole suivant dans la liste qui lui est associée. Si l'on a atteint la couche la plus haute et que, soit l'on a sélectionné tous les symboles PPM de ladite couche, soit la contribution de ladite couche calculée pour le symbole sélectionné excède le rayon quadratique de la sphère, le procédé de décodage se termine en fournissant ledit point le plus proche. Selon un exemple de réalisation, les symboles PPM des sources de rangs respectifs 1,..,P, estimés au sens du maximum de vraisemblance, sont obtenus comme sous- vecteurs (aopt , âopt , .. âop)t) à M composantes du vecteur à 5 MP composantes représentant ledit point le plus proche (âopt) . La contribution dp d'une couche de rang p , pE {1,..,P} est avantageusement calculée par : R(pmp) _el)) 2 avec E(p) = R(p'1 où 10 R(p'p) est la pieme sous-matrice de taille MxM sur la diagonale de la matrice triangulaire supérieure R de taille MPxMP, obtenue par transformation QR d'une 15 matrice H représentative du canal de transmission entre les P sources et ledit récepteur ; 'rIp est le résultat de l'égalisation ZF-DFE à la couche p ; R(p,p) est la mième colonne de la matrice R(p'p), où 20 mE {1,..,M} est la position de modulation du symbole PPM T sélectionné dans la couche p ; et il 2 =il . Le classement des positions mE {1,..,M} des symboles PPM de la couche p peut être réalisé en ordonnant ces symboles selon les valeurs : 2(R(pmp) )TE(p) + (RP) )TR(pmp). m =1,..,M 25 Le symbole PPM de la couche p réalisant la plus faible contribution correspond à la position de modulation pos(p) obtenue de préférence par : pos(p) = arg min 21R(P' )T E(P) + (RP)T R(PmP) m / L'invention est également définie par un récepteur à maximum de vraisemblance destiné à recevoir des 10 symboles PPM d'une pluralité P de sources, comprenant un filtre adapté au canal de transmission entre les sources et le récepteur et, le cas échéant, au codage spatio-temporel utilisé à l'émission, ledit récepteur comprenant en outre un décodeur par sphère comprenant 15 des moyens adaptés à exécuter les étapes du procédé de décodage par sphère selon l'une des revendications précédentes, le décodeur par sphère recevant en entrée la sortie dudit filtre adapté. 20 BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention fait en 25 référence aux figures jointes parmi lesquelles La Fig. 1 illustre le principe d'un décodage par sphère pour un système à deux sources PAM ; La Fig. 2 représente l'organigramme d'un premier algorithme de décodage par sphère connu de l'état de la 30 technique ;5 La Fig. 3 représente l'organigramme d'un second algorithme de décodage par sphère connu de l'état de la technique ; Les Figs. 4A et 4B illustrent le balayage de points 5 utilisé respectivement par le premier et le second algorithme de décodage par sphère ; La Fig. 5 représente l'organigramme d'un algorithme de décodage par sphère selon un mode de réalisation de l'invention ; 10 La Fig. 6 représente un exemple de progression à travers la constellation produit de l'algorithme de décodage par sphère selon l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS 15 Nous considérerons ci-après un système comprenant P sources de symboles PPM. L'alphabet de modulation comprend M positions temporelles. Comme indiqué plus haut les sources de symboles peuvent être des flux de symboles émis par les antennes d'un terminal MIMO, flux 20 de symboles en provenance ou à destination d'une pluralité d'utilisateurs, voire une combinaison de ces deux situations si chaque terminal utilisateur est équipé d'une pluralité d'antennes. Dans le cas général, K P = Eik où ik est le nombre d'antennes du terminal de k=1 25 l'utilisateur k et K le nombre d'utilisateurs. On supposera dans la suite que les flux de symboles sont émis de manière synchrone. Ces flux de symboles modulent avantageusement mais non nécessairement des trames de signaux UWB impulsionnels. Le récepteur est adapté à calculer P'M variables de décision avec P'> P. Par exemple, chaque antenne de réception est équipée d'un récepteur Rake adapté à séparer L trajets, et pour chaque trajet les M positions temporelles de la modulation PPM. Les variables de décision sont alors les P'M sorties de ces récepteurs Rake avec P'= L.Pr où Pr est le nombre d'antennes de réception. Si un codage spatio-temporel est utilisé à l'émission, un filtrage adapté à ce code est utilisé en réception pour séparer les différents utilisateurs et dans ce cas P'=L.Pr.P. Dans tous les cas, le signal reçu par le récepteur peut être représenté par un vecteur x de dimension P'M correspondant aux P'M valeurs des variables de décision. Le vecteur x peut s'exprimer, en adoptant les mêmes notations que pour l'équation (1), sous la forme x = Ha+n , à la différence près que le vecteur n est ici de dimension PMI, le vecteur a de dimension PM et la matrice H de taille P'MxPM . Le vecteur a peut être considéré comme la concaténation de P sous-vecteurs a(l),a(2),...,a(P) chaque sous-vecteur étant relatif à une source. La matrice H représente le canal de transmission et prend en compte les interférences multi-utilisateur et multi-trajet. Si R est la matrice de taille PMxPM obtenue par transformation QR de la matrice H, z=QTx est un résumé exhaustif du signal reçu appartenant à l'espace des signaux transmis. Le décodage par sphère consiste àrechercher le point â, appartenant à l'intersection de la constellation produit et d'une boule de bruit centrée sur z et pour lequel la distance à z est minimale. On appelle constellation produit la constellation générée par toutes les positions des constellations PPM élémentaires. Le récepteur détermine parmi les MP points de la constellation produit, c'est-à-dire parmi les MP combinaisons de symboles PPM possibles émis par les P sources, celle qui satisfait au critère de maximum de vraisemblance. On notera, comme précédemment, pour une matrice il quelconque de taille PMxPM, fl(''j) la sous-matrice de taille MxM de R constituée des éléments com,m, avec m = (iûl)M+1,..,iM et m'=(j-l)M+l,..,jM, avec 1 20 eP_1 = V(z û Rêp) , (17) ou encore, puisque R et V sont triangulaires supérieures . 25 e(P 1) = 0-1,p-11 ((') û R(P-1,P)êp(P) )+ 0-4) ((P) _ R(P,P)êPP) (18) et l'estimation ZF-DFE n'est autre que ap_1 =epP 1) =1.pos(P-1) où pos(Pù1) est l'indice de la première composante après l'opération de tri, comme précédemment indiqué. On procède ainsi de couche en couche, en éliminant à chaque couche p, l'interférence due aux P ù p couches supérieures. Lorsque les P couches ont été soumises à l'égalisation ZF-DFE, on dispose d'un vecteur â défini par la concaténation des â(P) pour p=1,..,P. Ce vecteur peut être représenté par un point de la constellation produit. On part de ce point pour rechercher le plus proche voisin de z. La recherche se poursuit ensuite au sein de ladite constellation, selon une méthode d'énumération exposée ci-après. On a vu plus haut que, pour chaque couche p, on déterminait le symbole â(P) le plus proche de 4P). Afin d'alléger les notations, nous noterons ici aP=â(P) et 11P =epP) . La contribution de cette couche à la distance quadratique entre z et â peut s'écrire : LR(P'P, ( m,m R(P,P) (al) ,qP 2 m=1 gym'=m (19) puisque R est triangulaire supérieure. al) étant un symbole PPM, a ,=8(m'ùm0) , mp E {1,..M} c'estûà-dire que 25 l'on suppose que la p1eme source transmet un symbole PPM en position m0 La contribution à la distance quadratique s'exprime alors : R(P'P) (û11 R(P'P) _ E(P) .,mo 2 2 (20) avec E(P) = R(P'P)rIP . On notera qu'en l'absence d'interférence entre les différentes positions de modulation, c'est-à-dire si la réponse impulsionnelle du canal de transmission est plus courte que l'écart temporel entre ces positions, la matrice R(P'P) est diagonale et que la distance quadratique dp se confond avec l'écart quadratique entre le symbole PPM testé (m = mp) et le signal égalisé par source. Pour chaque couche, on peut classer les symboles PPM (autrement dit les positions PPM) à partir des contributions respectives qu'ils apporteraient à la distance quadratique. Le premier choix dans l'énumération des candidats dans cette couche porte comme on l'a vu plus haut, sur : 2' (21) i pos(p) = arg min mR(PmP) _ E(P) puis l'on passe successivement aux positions réalisant des contributions de plus en plus élevées. Etant donné que E(P) ne dépend pas de m, le classement des distances quadratiques dp est équivalent à celui des produits scalaires : {_2(R,tTE(t ) )+(R(PmP))TR(PmP), m=1,..,M} (22) et par conséquent : pos(p) = arg min - 2(R(Pr f )T E(P) + (RP) )T R(P,P) (23) m / Le décodage par sphère commence par la dernière couche (couche P). On ordonne la liste des symboles PPM de cette couche en fonction des contributions respectives qu'ils apporteraient à la distance quadratique au point reçu. On sélectionne d'abord celui qui apporterait la plus faible contribution, ce qui donne une estimation du symbole émis par la source P. On élimine l'interférence générée par ce symbole sur toutes les couches inférieures. On passe ensuite à la couche immédiatement inférieure et on recommence ainsi pour chaque couche. Plus précisément, pour une couche courante donnée on ordonne la liste des positions PPM de cette couche en fonction des contributions qu'apporteraient les symboles correspondants à la distance quadratique au point reçu. On sélectionne celle réalisant la plus faible contribution en utilisant (21) ou (23). On élimine ensuite des couches inférieures l'interférence générée par le symbole PPM ainsi sélectionné. On progresse ainsi de couche en couche, les symboles sélectionnés dans les différentes couches forment une branche. Si la somme des contributions des couches supérieures et de la couche courante est supérieure au rayon courant de la sphère, on abandonne la branche en question et on passe à la couche supérieure à la couche courante pour sélectionner une nouvelle position PPM. A l'inverse, on poursuit la branche jusqu'à ce que l'on atteigne la couche de rang 1. La distance quadratique entre le point â de la constellation produit au point reçu est obtenue comme la somme des contributions des différentes couches. Si la distance est plus faible que le rayon courant de la sphère, le rayon ainsi que le meilleur point courant âopt sont mis à jour. Plus précisément, le nouveau rayon de la sphère prend la valeur de ladite distance et âopt=a. On retourne ensuite à la couche immédiatement supérieure pour sélectionner le second symbole PPM de la liste. On élimine des couches inférieures (en l'occurrence de la couche de rang 1) l'interférence générée par ce symbole et l'on procède comme précédemment vers les couches inférieures. A chaque remontée vers une couche supérieure, on sélectionne successivement les différents symboles de cette couche dans l'ordre de la liste en commençant par la contribution la plus faible. Lorsque la liste est épuisée, c'est-à-dire lorsque l'on a déjà précédemment testé les M symboles de cette liste, on passe à la couche supérieure suivante. On progresse ainsi, par remontées successives vers des couches de rangs de plus en plus élevés et avec des intervalles de validité de plus en plus étroits (étant donné que le rayon de la sphère est mis à jour), jusqu'à ce que l'on ait épuisé la liste de la couche racine (couche P), auquel cas l'algorithme s'arrête. En sortie de décodage, le meilleur point en sortie âopt est celui qui réalise le minimum de distance à z. Le point âopt fournit l'estimation ML des symboles PPM transmis par les P sources, soit âôpt , âopt , .. , âopt . La Fig. 5 représente de manière schématique l'organigramme du procédé de décodage par sphère d'une pluralité de sources PPM selon l'invention. A l'étape 510, on initialise l'indice de couche par celui de la plus haute, p = P. On initialise les variables intermédiaires, notamment le rayon quadratique d de la sphère, d = D et P variables représentant respectivement la somme des contributions à la distance quadratique jusqu'aux P couches p respectives 6p = Edp , soit 6p =O pour p =1,..,P. k=p A l'étape 515, on effectue l'égalisation ZF-DFE de la couche courante p compte tenu des symboles â(P) , â(P-1) , .. , â(p+l) sélectionnés dans les couches supérieures précédentes. A l'étape 520, on classe les symboles PPM de cette couche selon la relation d'ordre (20) ou (22). On dispose ainsi d'une liste ordonnée L(p) et on choisit le premier symbole PPM de cette liste. A l'étape 525, on met à jour 6p, la somme des contributions des couches p,..,P à la distance euclidienne à z, soit 6p =6p+1 +dp . On teste ensuite si la somme de ces contributions excède le rayon quadratique courant d de la sphère. Si c'est le cas, il est inutile de conserver la branche courante â(P),â(P-1),...,â(t) puisque elle ne pourra, quels que soient les symboles que l'on sélectionnera dans les couches inférieures suivantes, conduire au plus proche voisin. On passe alors au test 535. Par contre si 6p est inférieure au dit rayon quadratique on passe au test 540. En 535, on teste si p = P. Si c'est le cas l'algorithme s'arrête en 537 et le plus proche voisin de z appartenant au réseau est âopt• Si p =P , on incrémente p et pos(p) de 1 en 555, c'est-à-dire que l'on passe à la couche supérieure et on sélectionne le symbole PPM suivant dans la liste L(p), s'il existe. On teste en 560 si l'on a épuisé les symboles de la liste L(p). Dans l'affirmative (pos(p)>M), on retourne au test 535. Dans la négative, on met à jour la valeur de dp correspondant au nouveau symbole sélectionné et on retourne à l'étape 525. En pratique, si l'on a 25 préalablement calculé et stocké les valeurs dp pour les différentes valeurs de m =1,..,M , il n'est pas nécessaire de recalculer dp. Si l'on décide de conserver la branche courante â(P),â(P1),...,e) en 530, on teste en 540 si cette branche a atteint la couche la plus basse Q,=1), auquel cas on dispose d'un point du réseau améliorant la distance courante à z et l'on passe à l'étape 550. Si, en revanche ce n'est pas le cas (p #1), on passe à la couche inférieure suivante en 545, puis on retourne à l'étape d'égalisation ZF-DFE. A l'étape 550, on met à jour le rayon de la sphère, soit d =G1 et le meilleur candidat âopt = â , puis on poursuit la recherche en incrémentant p à l'étape 555. Le meilleur candidat donne une estimation ML conjointe âopt âopt , .. âopt des symboles PPM émis par les P sources. On a illustré en Fig. 6 la progression du décodage par sphère au sein de la constellation produit. L'espace des symboles transmis est divisé en P couches, chacune étant de dimension A T. On a représenté en ordonnées pour chacune des positions PPM de chaque couche p les M composantes de epp), obtenues comme on l'a vu par égalisation ZF-DFE. La branche notée (1) est la première branche testée. Dans le cas présent, la branche progresse jusqu'à la couche la plus basse et améliore la distance quadratique à z. Le rayon quadratique d de la sphère est alors mis à jour. On passe à la couche 2 et le test d'une nouvelle branche (2). Etant donné que l'on a sélectionné un nouveau candidat dans la couche (2), on procède préalablement à une nouvelle égalisation de la couche 1. On a représenté par (e) une branche en cours de test. On remarque que lorsque cette branche atteint la couche p, la somme 6p des contributions des couches supérieures et de la couche courante P,Pû1,..,p, excède le rayon quadratique courant d de la sphère. La recherche se poursuit en passant à la couche supérieure p +l et en sélectionnant dans cette couche le symbole PPM suivant de la liste 41+1) (cf. étape 555 de la Fig. 5). On actualise alors la somme des contributions, 6 p+l Si cette somme est inférieure à d on poursuit avec la nouvelle branche (+1), sinon on remonte d'une couche. Lorsque l'on atteint la couche la plus haute (cf. étape 535 de la Fig. 5) et que soit on a épuisé les symboles de la liste L(P), soit la contribution de la couche P excède d, l'algorithme de décodage se termine en fournissant âopt On a donné en annexe le pseudo-code de la méthode de décodage par sphère selon une réalisation possible de l'invention. Les notations suivantes ont été adoptées . OM est le vecteur nul de taille M ; p est un vecteur de taille M ; e et E sont des matrices de taille MP x P, les colonnes de e donnent le résultat de l'égalisation ZFDFE ; A matrice de taille MxP dont les colonnes sont les symboles PPM sélectionnés pour les différentes couches ; Order matrice de taille MxP dont les colonnes donnent les listes d'indices L(p) précitées pour chacune des couches ; Diag(SL) est la matrice diagonale composée des éléments diagonaux de la matrice il ; sort(0) est un vecteur constitué des composantes du vecteur 0 triées par ordre croissant. La méthode de décodage par sphère exposée ci-dessus est destinée à être implémentée grâce à des moyens logiciels ou matériels dans un récepteur ML opérant dans un système multi-source, par exemple un système MIMO et/ou un système multi-utilisateur. Avantageusement, les symboles PPM d'une source servent à moduler un signal impulsionnel du type TH-UWB (Time Hopped UWB), DS-UWB (Direct Spread UWB) ou TH-DS-UWB. Les signaux ainsi modulés sont transmis par exemple grâce à des antennes UWB conventionnelles ou des diodes laser. Annexe Décodage par sphère (input: z , R , P ,M ,D; output: âopt 5 Step 1 (Initialisation): Set k = P +1 ; dist(k) = 0 ; V = R1; p = 0M ; bestdist = D. Step 2 : newdist = dist(k)+ pTp; if ( newdist M , go to Step 4 endif, pos(k) = Orderstep(k),k , A•,k = I.,pos(k) 25 p = E(k) ù R(pos(k), go to Step 2. 10 REVENDICATIONS 1. Procédé de décodage par sphère pour récepteur à maximum de vraisemblance destiné à recevoir des symboles PPM d'une pluralité P de sources, chaque source émettant un flux de symboles PPM à M positions de modulation, P symboles PPM émis simultanément par les P sources étant représentés par un point (a) de la constellation de modulation produit dans un espace des signaux transmis de dimension MP décomposé en P couches, chaque couche représentant les M positions de modulation possibles d'un symbole PPM émis par cette source, le signal reçu (x) par ledit récepteur étant transformé en un point représentatif (z) de ce signal, dit point reçu, dans l'espace des signaux transmis, ledit procédé déterminant le point de la constellation produit le plus proche (âopt) du point reçu à l'intérieur d'une sphère de rayon quadratique donné, caractérisé en ce que, pour chaque couche de rang p . (a) on effectue une égalisation ZF-DFE du signal reçu dans ladite couche compte tenu des symboles PPM estimés dans les P û p couches précédentes, dites couches supérieures à la couche p; (b) on classe dans une liste les M symboles PPM de ladite couche en fonction des contributions qu'ils apporteraient à la distance quadratique au point reçu ; (c) on sélectionne le symbole PPM réalisant la plus faible contribution (dp) et l'on ajoute cette contribution à celles obtenues pour les couchesprécédentes pour obtenir une somme de contributions (6p); les étapes (a),(b),(c) étant répétées jusqu'à ce que l'on atteigne la couche la plus basse ; et l'on met à jour le rayon quadratique de la sphère et ledit point le plus proche si ladite somme de contributions est inférieure au rayon quadratique de la sphère. 2. Procédé de décodage par sphère selon la revendication 1, caractérisé en ce que, si pour une couche donnée et un symbole PPM sélectionné dans cette couche, ladite somme des contributions excède le rayon quadratique de la sphère, on passe à la couche supérieure et on sélectionne dans cette couche le symbole suivant de la liste qui lui est associée. 3. Procédé de décodage par sphère selon la revendication 2, caractérisé en ce que si tous les symboles de ladite liste associée ont déjà fait l'objet d'une sélection, on passe à la couche encore supérieure pour y sélectionner dans cette dernière le symbole suivant dans la liste qui lui est associée. 4. Procédé de décodage par sphère selon la revendication 2 ou 3, caractérisé en ce que si l'on a atteint la couche la plus haute et que, soit l'on a sélectionné tous les symboles PPM de ladite couche, soit la contribution de ladite couche calculée pour le symbole sélectionné excède le rayon quadratique de lasphère, le procédé de décodage se termine en fournissant ledit point le plus proche (âopt)• 5. Procédé de décodage par sphère selon la revendication 4, caractérisé en ce que les symboles PPM des sources de rangs respectifs 1,..,P, estimés au sens du maximum de vraisemblance, sont obtenus comme sous- vecteurs ( aopt , âopt , .. âop)t) à M composantes du vecteur à MP composantes représentant ledit point le plus proche ( aopt) . 6. Procédé de décodage par sphère selon l'une des revendications précédentes, caractérisé en ce que la contribution dp d'une couche de rang p , pE {1,..,P} est calculée par : R(pmp) _el)) 2 avec E(p) = R(p'1 où dp = R(p'p) est la p1eme sous-matrice de taille MxM sur la diagonale de la matrice triangulaire supérieure R de taille MPxMP, obtenue par transformation QR d'une matrice H représentative du canal de transmission entre les P sources et ledit récepteur ; 'rIp est le résultat de l'égalisation ZF-DFE à la couche p R(P'P) est la m1ème colonne de la matrice R(P'P), où mE {1,..,M} est la position de modulation du symbole PPM sélectionné dans la couche p ; et S 2= n . 7. Procédé de décodage par sphère selon la revendication 6, caractérisé en ce que le classement des positions mE {1,..,M} des symboles PPM de la couche p est réalisé en ordonnant ces symboles selon les valeurs . 1.2(R(P,P) )TE(P) + (RP))TR(PmP), =1,..,M 8. Procédé de décodage par sphère selon la revendication 7, caractérisé en ce que le symbole PPM de la couche p réalisant la plus faible contribution correspond à la position de modulation pos(p) obtenue par : 2(R(Pmp) )T E(p) + (R(pmP) )T R(pmp) i 9. Récepteur à maximum de vraisemblance destiné à recevoir des symboles PPM d'une pluralité P de sources, comprenant un filtre adapté au canal de transmission entre les sources et le récepteur et, le cas échéant, au codage spatio-temporel utilisé à l'émission, caractérisé en ce qu'il comprend en outre un décodeur par sphère comprenant des moyens adaptés à exécuter les étapes du procédé de décodage par sphère pos(p) = arg min mselon l'une des revendications précédentes, le décodeur par sphère recevant en entrée la sortie dudit filtre adapté.