FR 3077941 A1 20190816 FR 1800146 A 20180215 PROCEDE POUR CALCULER UNE ESTIMATION D'UN SIGNAL NUMERIQUE MODULE ET DE SA FIABILITE L'invention concerne un procédé pour améliorer le calcul de l'estimation des symboles d'un signal numérique modulé et de leur fiabilité, pour toute tâche qui nécessite une telle estimation. Dans le cadre des communications numériques, elle peut être utilisée pour des récepteurs de communication sans fil, ou pour des récepteurs de communications filaires. Par exemple, l'estimation des symboles présentée ici peut être utilisée pour dériver une estimation du canal (ou de certains de ses paramètres) entre l'émetteur et le récepteur. Elle peut être utilisée pour dériver une estimation des paramètres du signal ou du récepteur comme, par exemple, le décalage d'horloge, ou l'instant de réception du signal. Elle peut être utilisée aussi pour dériver des métriques de qualité du lien, comme par exemple l'information is mutuelle entre les symboles émis et ceux reçus, la fiabilité de l'estimée du signal envoyé, le rapport signal sur bruit ou signal sur bruit plus interférence. Ces métriques sont utilisées avec différentes finalités, comme l'adaptation du lien, par exemple l'adaptation des paramètres d'une chaine de transmission, ou alors à la réception pour optimiser le point de fonctionnement du zo récepteur ou d'un ou plusieurs de ses blocs fonctionnels internes, en choisissant les bons algorithmiques et les paramétrages les plus adaptés, en fonction de la flexibilité du récepteur même. L'invention concerne aussi un procédé pour supprimer des interférences au sein d'un signal reçu sur un récepteur. Elle s'applique 25 notamment à de l'égalisation avec un récepteur adaptatif de type égaliseur à retour de décision ou DFE (Digital Feedback Equalizer) permettant une égalisation des symboles reçus. Le domaine de l'invention est celui des systèmes de communications numériques par voie radio et entre autres les émetteurs et 30 les récepteurs de communications multivoies, c'est-à-dire comprenant plusieurs antennes. 2 L'invention concerne également les systèmes multi-utilisateurs dans lesquels les ressources de communications sont partagées entre plusieurs utilisateurs qui peuvent communiquer simultanément par un partage de bandes fréquentielles ou de créneaux ou slots temporels. 5 Plus généralement, l'invention concerne tous les systèmes de communication multi-utilisateurs dans lesquels de forts niveaux d'interférences sont générés à la fois entre des émetteurs associés à des utilisateurs différents mais également entre les symboles véhiculés par un signal transmis par un utilisateur du fait des perturbations inhérentes au io canal de propagation. Afin de supprimer ou au moins limiter l'interférence générée sur le signal reçu, il est connu de mettre en oeuvre un procédé de suppression d'interférence au sein du récepteur. Cette fonctionnalité a pour but de nettoyer le signal reçu des diverses sources d'interférences avant son décodage. Dans ce contexte, l'invention concerne précisément le domaine de la suppression d'interférence dans un contexte multi-utilisateurs dans le cadre des récepteurs itératifs, qui consiste en l'itération des fonctions de suppression d'interférence et de décodage dans l'objectif final d'améliorer le taux d'erreur bit ou le taux d'erreur paquet sur les symboles décodés. L'invention trouve entre autres son application dans les systèmes de communication cellulaires tels que le système 3GPP LTE. L'invention vise également une méthode d'égalisation dans le domaine fréquentiel et qui soit flexible, adaptée à un contexte soit mono- utilisateur, soit multi-utilisateurs et avec un récepteur mettant en oeuvre un traitement soit non-itératif, soit itératif. Un transfert d'information d'une source à une destination implique la propagation par un canal qui peut être, par exemple, un canal radio, un canal filaire (tel qu'un câble coaxial), etc. Certains moyens de propagation génèrent une interférence dite inter-symbole sur le signal reçu. Autrement dit, le signal reçu échantillonné à un instant donné, après compensation des 3 retards de propagation et de traitement, et ayant une synchronisation correcte, ne contient pas seulement le symbole envoyé (éventuellement amplifié et avec une perturbation de phase) plus du bruit, mais un mélange ou combinaison linéaire de symboles envoyés. 5 Il est usuel dans ce cas d'utiliser des algorithmes d'égalisation pour réduire l'effet parasite de l'interférence inter-symbole (ISI). De nombreuses méthodes d'égalisation sont décrites dans l'art antérieur. L'objectif est d'atteindre des performances optimales données par la borne du filtre adapté (qui est une borne inférieure sur le taux d'erreur paquet) tout io en étant facile d'utilisation. Les algorithmes doivent donc avoir une complexité de calcul, une occupation de mémoire et une latence de traitement compatibles, à la fois aux applications qui utilisent les récepteurs en question et aux contraintes des plateformes matérielles sur lesquelles ils sont implémentés. Le problème technique posé est donc de trouver des 15 algorithmes qui trouvent un bon compromis performances/complexité d'implémentation par rapport aux applications visées. Il existe plusieurs classes d'égaliseurs : les égaliseurs linéaires, les égaliseurs à retour de décision (DFE ou Decision Feedback Equalizers), les égaliseurs à effacement d'interférences (Interference Cancellation) et les 20 détecteurs du maximum a posteriori (Maximum A Posteriori - MAP), ou des détecteurs qui estiment une séquence au maximum de vraisemblance (Maximum Likelihood Sequence Estimation - MLSE). Ces égaliseurs peuvent être déclinés en forme itérative. Des récepteurs itératifs particuliers sont les récepteurs dits « turbo », où il existe un échange réitéré 25 d'information probabiliste extrinsèque entre blocs de traitements, par exemple entre le bloc d'égalisation et le bloc de décodage. Dans le domaine des égaliseurs itératifs, les égaliseurs de l'état de l'art peuvent se distinguer en fonction du type de signal de retour : Retour dur (DFE conventionnel) ou souple à partir du démodulateur 30 ou du décodeur de canal [13], Retour extrinsèque du décodeur de canal [2], 4 - Retour a posteriori du décodeur de canal [14]. Ces types de récepteurs élémentaires ont été étendus à des cas de signaux soumis à interférences inter-symboles, inter-antennes ou inter-utilisateurs [12]. 5 Plus en détail, la structure de l'égaliseur DFE fréquentiel classique a été présentée dans la référence [13] et étendue au cas d'égaliseur fractionnel dans [8]. Ces articles présentent des égaliseurs sans boucle turbo, où le signal de retour est soit de type classique dur, soit souple et calculé avec la probabilité a posteriori sur le signal égalisé. Les égaliseurs de io type fréquentiel ont été ensuite étendus pour prendre en compte un signal de retour calculé à partir de probabilités extrinsèques fournies par le décodeur, il s'agit des égaliseurs de type « turbo », voir par exemple [2]. Dans l'article [14] les auteurs montrent sur un cas de détecteur linéaire pour des systèmes multi-antennaires (on traite ici l'interférence entre antennes) qu'un signal de 15 retour calculé à partir de probabilités a posteriori fournies par le décodeur peut améliorer les performances du système par rapport à un signal de retour calculé sur les probabilités extrinsèque. Les auteurs de l'article [12] appliquent les retours extrinsèques en provenance du décodeur à un cas de transmissions sans fil multiutilisateurs où les émetteurs sont dotés 20 d'antennes multiples et appliquent un code espace-temps. Thomas Minka propose dans son rapport technique [3] la propagation d'espérance (Expectation Propagation - EP - en anglais), une technique de l'inférence bayésienne, qui est un algorithme itératif pour l'estimation de densité de probabilité a posteriori. Cette méthode mathématique générique peut être 25 utilisée pour approximer la solution d'un problème de maximisation de la probabilité a posteriori ou MAP de façon itérative. La propagation d'espérance EP peut aussi être vue comme un algorithme de passage de messages connu sous l'expression anglo-saxonne « message passing », en particulier comme une généralisation de l'algorithme de propagation de 30 croyance (Belief Propagation - BP - en anglais) aux cas de distributions de probabilité non-catégoriques mais appartenant à la famille exponentielle. En 5 particulier, ce concept permet d'effectuer un meilleur « message passing » avec des estimateurs à erreur quadratique moyenne minimum (Minimum Mean Square Error - MMSE). Dans le domaine des communications numériques et en particulier des récepteurs turbo, le concept « d'expectation 5 propagation » permet de calculer un type de retour souple différent. En particulier, la référence [6] utilise ce nouveau retour souple dans le cadre d'un récepteur à entrées multiples sorties multiples, (Multiple-Input Multiple-Output - MIMO), ce qui permet de mieux approcher les performances MAP. Dans le cadre de l'égalisation, la référence [15] étudie un signal de retour io basé sur le concept de l'EP à partir du démodulateur souple, pour les égaliseurs linéaires par bloc dans le domaine temporel. En couplant ce retour avec le retour extrinsèque en provenance du décodeur, dans une structure avec une double boucle (égaliseur/démodulateur, démodulateur/décodeur), les auteurs obtiennent de bonnes performances. 15 La référence [16] traite de l'égalisation fréquentielle avec EP dans le cas multi-utilisateurs où les émetteurs mobiles sont équipés d'une seule antenne de transmission. Dans une des approches proposées, le retour EP est calculé sur une distribution Gaussienne coloré, c'est-à-dire, en laissant chaque symbole estimé au démodulateur avoir une mesure de fiabilité propre 20 à lui. Cela contraint le récepteur fréquentiel à avoir une structure très complexe, nécessitant l'inversion d'une matrice pleine à chaque bloc de données. Un récepteur alternatif moins complexe est aussi dérivé, avec une hypothèse de moindres carrés, mais l'égaliseur dévient alors inutilisable dès lors que le canal de propagation présente des zéros spectraux. 25 Les méthodes connues de l'art antérieur ne permettent pas de prendre en compte un ou plusieurs des problèmes listés ci-après, ou encore conjointement, les aspects suivants : La mise en oeuvre d'une chaîne de réception numérique itérative qui inclut un récepteur adaptatif à suppression d'interférence qui 30 transforme le canal de propagation en un canal équivalent à bruit additif Gaussien, 6 Le récepteur adaptatif à suppression d'interférence est auto-itéré avec un démodulateur qui lui fournit un retour qui est une estimation du signal transmis, Le calcul d'un retour de type propagation d'espérance au sein du 5 démodulateur souple, sur une distribution gaussienne blanche, c'est- à-dire où l'ensemble des symboles estimés du bloc de données est caractérisé par une seule mesure de fiabilité, L'utilisation d'un procédé de lissage adaptatif sur les symboles estimés par le retour EP, et sur la variance caractérisant la fiabilité ce 10 retour. Ce procédé permet de choisir un compromis entre la robustesse des estimations et la vitesse de convergence du procédé itératif, L'utilisation éventuelle d'une technique de protection supplémentaire contre des estimations erronées, en remplaçant le retour EP par un 15 retour a posteriori, suivant la valeur de la fiabilité des estimés EP. La prise en considération conjointe de ces points permet alors l'utilisation d'un démodulateur souple à retour EP blanchi pour adresser directement les problématiques suivantes (entre autres) : La mise en oeuvre d'un récepteur traitant l'interférence entre symboles 20 dans le domaine fréquentiel, La mise en oeuvre d'un récepteur fractionnaire qui traite de l'égalisation fréquentielle du canal avec un facteur de sur-échantillonnage, La mise en oeuvre d'un récepteur fréquentiel avec recouvrement de 25 blocs de données, où le retour EP est par ailleurs utilisé pour la réduction de l'interférence entre blocs, L'utilisation des retours EP au sein d'un détecteur multi-utilisateur et où les utilisateurs peuvent être équipés d'un nombre variable (une ou plusieurs) d'antennes de transmission. Le retour EP blanchi est alors 30 utilisé pour traiter l'interférence entre utilisateurs et l'interférence multi- antennaire, en plus de l'interférence entre symboles. 7 La suite de la description fait appel aux notations suivantes : x dénote un scalaire; x dénote un vecteur; X dénote une matrice; 5 [x1] j1 est le vecteur contenant les coefficients xi; [x Li R est une matrice par bloc qui contient les matrices )(id; j=1...M xT est le vecteur transposé de x (pareil pour les matrices); x est le vecteur conjugué de x (pareil pour les matrices); x" est le vecteur hermitien, c'est-à-dire le vecteur transposé et conjugué, de 10 x (pareil pour les matrices); x dénote la représentation en domaine fréquentielle du vecteur x; X dénote la représentation en domaine fréquentielle de la matrice X; tr(X) est la trace de la matrice X; 0 dénote le produit de Kronecker; 15 Q dénote le produit terme à terme entre deux vecteurs (c =a0b est tel que cri = ab), aussi appelé produit matriciel de Hadamard, ou produit de Schur ; circ(x) est une matrice circulante avec première colonne égale à x; diag(x) est une matrice diagonale avec les coefficients x sur sa diagonale; 20 b : vecteur des bits d'information ; d: vecteur des bits codés et entrelacés ; E[x]: espérance de x selon la distribution de probabilités P, E[x] = EaEX aP(a) ; Fm : la matrice de transformée de Fourier de taille M x M, l'entrée k,1 de la 1 25 matrice est [Fm]kj = 7/Tiexp (-2kl m ), k = 0, ...,M -1 et 1 = 0, ...,M -1. L'inverse de la matrice de transformée de Fourier coïncide avec sa transposée conjugué : ; : matrice identité de taille NxN ; M: taille d'un vecteur de données à égaliser, correspondant aussi à la taille 30 de FFT dans un émetteur de type SC-FDMA ; 8 : fonction d'étiquetage binaire d'un symbole de la constellation ; N: taille de la IFFT dans un émetteur de type SC-FDMA, correspondant au nombre total de sous-porteuses du modulateur/démodulateur. Cela aussi est équivalent au nombre d'échantillons temporels ; 5 R: nombre d'antennes de réception ; s : identifiant entier de l'auto-itération (entre 0 et Si) ; : nombre d'auto-itérations du récepteur (nombre entier égal ou supérieur à zéro) pour la -r-ième turbo-itération ; : identifiant entier de la turbo-itération (entre 0 et T); 10 : nombre de turbo-itérations du récepteur (nombre entier égal ou supérieur à zéro) ; Tu : nombre d'antenne de transmission de l'utilisateur u; U: nombre d'utilisateurs (notamment d'émetteurs) ; : fiabilité globale de l'estimation ï; 15 Van, [X] : variance de x selon la distribution de probabilité P, V arp[x] = Es[lx - ET[x]12] ; X u : constellation des symboles de l'utilisateur u; ï: estimation du symbole transmis x fournie par le démodulateur ; : estimation du symbole transmis x fournie par le récepteur adaptatif à 20 suppression d'interférence ; Ew : la matrice de covariance du bruit dans le domaine fréquentiel L'invention concerne un procédé pour améliorer l'estimation des symboles d'un signal modulé et de la fiabilité de l'estimée, le signal comprenant des blocs de données le procédé étant exécuté au niveau d'un 25 processeur et comportant au moins les étapes suivantes : - Une étape de calcul d'estimations souples du signal par une technique de maximum de vraisemblance exacte ou approximée, en utilisant les symboles estimés X du bloc de données courant, les symboles estimés du bloc étant caractérisés par une mesure de fiabilité moyenne 9 - Une étape de division gaussienne du signal estimé par maximum de vraisemblance, supposé blanc, par le signal estimé X pour générer une estimation souple du signal transmis, l'estimation étant caractérisée par une mesure de fiabilité par bloc de données, 5 - Une étape de lissage adaptatif sur le signal estimé par la division gaussienne et sur les mesures de fiabilités, pour fournir une estimation X du signal transmis et une mesure de fiabilité 17 associée. Selon une variante de réalisation, le procédé comporte en plus : - Une étape de calcul de distributions a priori du signal transmis à partir 10 d'entrées/sorties souples binaires, - Une étape de calcul d'estimations souples du signal transmis, par la technique du maximum a posteriori, en utilisant le signal estimé par le maximum de vraisemblance et les distributions a priori calculées à l'étape précédente, ladite estimation est caractérisée par une mesure 15 de fiabilité moyenne. Le procédé peut aussi comporter : - Une étape de décodage par un décodeur binaire à entrées/sorties souples, - Une étape d'initialisation du démodulateur souple, en calculant une estimation a priori du signal transmis, ladite estimation étant caractérisée par 20 une fiabilité moyenne en amont de l'étape de calcul d'estimations souples du signal transmis. Selon une variante de réalisation, le procédé comprend en outre une étape de désentrelacement et d'entrelacement des informations souples sur le signal transmis, lors de l'échange itératif entre un démodulateur souple 25 et un décodeur binaire à entrées/sorties souples. Le procédé est utilisé pour la suppression d'interférence au sein d'un signal issu d'au moins un émetteur reçu sur un récepteur comprenant une ou plusieurs antennes, ledit signal comprenant un ou plusieurs blocs de données et comporte en outre les étapes suivantes : 30 - Une étape de calcul des paramètres du récepteur, en utilisant la connaissance du canal, 10 - Une étape de traitement du signal reçu par un récepteur à suppression d'interférence qui génère une autre estimation X du signal transmis en utilisant une estimation préalable X du signal transmis x et la connaissance du canal, l'estimation du signal étant effectuée par le s récepteur à suppression d'interférence, - Le calcul de la fiabilité o-, de l'estimation du signal transmis X obtenue par le récepteur à suppression d'interférence, en utilisant les paramètres du récepteur, - Une étape d'échange itératif des estimations du signal transmis, entre 10 le récepteur à suppression d'interférence et le démodulateur souple, - Une étape de calcul d'information souple sur les bits du signal transmis, faite à partir du signal démodulé. Le procédé peut aussi comporter les étapes suivantes : - Une étape de calcul des paramètres du récepteur adaptatif, en 15 utilisant la connaissance du canal et la fiabilité des retours souples du démodulateur, - Une étape d'échange itératif des fiabilités des estimations du signal transmis entre le bloc qui calcule les paramètres du récepteur adaptatif et le démodulateur souple. 20 Selon une variante de réalisation, le procédé pour supprimer l'interférence au niveau d'un signal comporte les étapes suivantes : - Une étape d'estimation des paramètres du canal entre l'émetteur et le récepteur en utilisant le signal reçu, - Une étape de calcul des paramètres du récepteur adaptatif ou non, en 25 utilisant l'estimation des paramètres du canal et la fiabilité des retours souples du démodulateur, - Une étape d'échange itératif des fiabilités des estimations du signal transmis entre le bloc qui calcule les paramètres du récepteur adaptatif et le démodulateur souple. 11 Le procédé peut comprendre en outre une étape itérative d'estimation des paramètres du canal, en utilisant l'estimation souple du signal transmis et sa fiabilité en provenance du démodulateur souple, conjointement avec le signal reçu. s L'étape de traitement du signal reçu par un récepteur à suppression d'interférence et l'étape de calcul de ses paramètres comprennent, par exemple: - Une étape de conversion du signal reçu dans le domaine fréquentiel, - Une étape de calcul des paramètres de l'égaliseur dans le domaine 10 fréquentiel, qui réalise les fonctions du bloc, - Une étape d'égalisation linéaire à suppression d'interférence dans le domaine fréquentiel, paramétrée selon les sorties obtenues lors du calcul des paramètres de l'égaliseur dans le domaine fréquentiel, - Une étape de conversion du signal traité dans le domaine temporel, 15 - Une étape de conversion dans le domaine fréquentiel de l'estimation souple du signal transmis fournie par le démodulateur souple. L'étape de traitement du signal reçu par un récepteur à suppression d'interférence et l'étape de calcul de ses paramètres comprennent, par exemple: 20 - Une étape de séparation du signal reçu y en P segments yp de Np symboles chacun, avec recouvrement, p = 1, ..., P. - Une étape de conversion de chacun de segments du signal reçu dans le domaine fréquentiel, - Une étape du calcul des paramètres des P égaliseurs dans le 25 domaine fréquentiel, qui réalise les fonctions du bloc, - Une étape d'égalisation linéaire à suppression d'interférence dans le domaine fréquentiel de chacun des segments, paramétrée selon les sorties obtenues lors du calcul des paramètres des P égaliseurs dans le domaine fréquentiel, 12 - Une étape de conversion de chacun des signaux traités dans le domaine temporel, et d'extraction des symboles non recouverts pour obtenir une estimation X, - Une étape de conversion dans le domaine fréquentiel de l'estimation 5 souple du signal transmis fournie par le démodulateur souple. Le calcul des paramètres des P égaliseurs dans le domaine fréquentiel peut utiliser la connaissance de la statistique de l'interférence entre segments, estimée à partir de la fiabilité des retours du démodulateur souple. 10 L'interférence résiduelle entre segments est régénérée, par exemple, en utilisant les retours du démodulateur souple, et soustraite de chacun des P segments yp, avant l'étape de conversion dans le domaine fréquentiel. Selon un mode de réalisation l'étape de traitement du signal reçu 15 par un récepteur à suppression d'interférence et l'étape de calcul de ses paramètres comprennent: - Un sur-échantillonnage du signal reçu avec un facteur de sur-échantillonnage o, - Une étape de conversion en fréquence du signal reçu 20 suréchantillonné, - Une étape de calcul du filtre de l'égaliseur dans le domaine fréquentiel suréchantillonné, - Une étape d'égalisation linéaire à suppression d'interférence dans le domaine fréquentiel suréchantillonné selon la sortie obtenue lors du 25 calcul des paramètres de l'égaliseur dans le domaine fréquentiel suréchantillonné, - Une étape de conversion du signal traité dans le domaine temporel, suivi d'un échantillonneur de facteur d'échantillonnage o égal au facteur de sur-échantillonnage du signal reçu. 13 - Une étape d'interpolation idéale de l'estimation souple du signal transmis fournie par le démodulateur souple par un facteur de suréchantillonnage o égal à celui du signal reçu (cela consiste en l'ajout de o - 1 zéros entre deux échantillons), suivi d'une conversion 5 dans le domaine fréquentiel suréchantillonné. Le procédé est utilisé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type SC-FDMA ou SS-SC-FDMA comprenant une étape de mise en trame utilisant respectivement une modulation SC-FDMA ou SS-SC-FDMA et une étape d'exécution du procédé io selon l'invention. Le procédé peut aussi être utilisé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type SC ou SS-SC comprenant une étape de mise en trame utilisant respectivement une modulation SC ou SS-SC et une étape d'exécution du procédé selon 15 l'invention. L'invention concerne aussi un procédé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type OFDM ou OFDMA comprenant une étape de mise en trame utilisant respectivement une modulation OFDM ou OFDMA et une étape d'exécution du procédé 20 selon l'invention dans lequel, l'étape de traitement du signal reçu par un récepteur à suppression d'interférence et l'étape de calcul de ses paramètres comprennent: - Une étape de conversion du signal reçu dans le domaine fréquentiel, - Une étape de calcul des paramètres de l'égaliseur dans le domaine 25 fréquentiel, qui réalise les fonctions du bloc, - Une étape d'égalisation linéaire à suppression d'interférence dans le domaine fréquentiel, paramétrée selon les sorties obtenues lors du calcul des paramètres de l'égaliseur dans le domaine fréquentiel. Selon une variante, le procédé de suppression d'interférence 30 comprend en outre: 14 - Une étape de génération du signal à l'émetteur selon un code espace-temps par blocs de type Alamouti ou autre code espace-temps par bloc de type orthogonal, - Une étape d'estimation des paramètres des canaux entre les 5 antennes de l'émetteur et celles du récepteur, - Une étape du calcul des paramètres du récepteur dans le domaine fréquentiel, qui prennent en compte la structure du code espace-temps par blocs utilisé, - Une étape d'égalisation linéaire à suppression d'interférence dans le 10 domaine fréquentiel en utilisant les filtres obtenus lors du calcul des paramètres des égaliseurs dans le domaine fréquentiel, et qui combine les signaux reçus en entrée selon la méthode associé au code espace-temps par blocs utilisé par l'émetteur. L'invention concerne aussi un procédé pour supprimer des 15 interférences au sein des signaux issus d'un émetteur reçus sur un récepteur, dans un réseau de type MIMO où l'émetteur est doté de plusieurs antennes de transmissions et d'un entrelaceur qui répartit les bits codés sur les antennes disponibles, comprenant une étape d'exécution du procédé selon l'invention comprenant en outre : 20 - Une étape de traitement des signaux reçus en provenance des antennes de l'émetteur par un récepteur à suppression d'interférence, qui en utilisant ses paramètres et une estimation préalable des signaux transmis Xt par les antennes d'émission t = 1, ..., T1, génère une autre estimation des signaux émis .ît, 25 - Une étape du calcul des paramètres du récepteur qui prend en compte la présence de plusieurs antennes d'émission et les estimations des paramètres des canaux entre les antennes de l'émetteur et du récepteur, - Une étape de démodulation souple où l'estimation du signal de 30 chaque antenne d'émission est traitée indépendamment. 15 L'invention concerne aussi un procédé pour supprimer des interférences au sein des signaux issus de plusieurs émetteurs reçus sur un récepteur, dans un réseau de type MU-MIMO où les émetteurs sont dotés d'une ou plusieurs antennes de transmissions et d'un entrelaceur qui répartit 5 les bits codés sur les antennes disponibles, comprenant une étape d'exécution du procédé selon l'invention comprenant: - Une étape de traitement des signaux reçus en provenance des émetteurs par un récepteur à suppression d'interférence, qui en utilisant ses paramètres et une estimation préalable des signaux 10 transmis ït,u par les antennes d'émission t = 1, ..., Tu des utilisateurs= 1, U, génère une autre estimation des signaux émis - Une étape du calcul des paramètres du récepteur qui prend en compte la présence de plusieurs utilisateurs, éventuellement avec plusieurs antennes d'émission, et d'éventuelles estimations des 15 paramètres des canaux entre les antennes des émetteurs et du récepteur, - Une étape de démodulation souple où l'estimation du signal de chaque antenne d'émission de chaque utilisateur est traitée indépendamment, 20 - Une étape d'ordonnancement des utilisateurs à décoder selon un critère basé sur la puissance des signaux reçus, ou selon une métrique dérivée de la puissance des signaux reçus et de la fiabilité de leurs constellations, ou encore selon une métrique basée sur les écarts de puissances et fiabilité entre couples d'utilisateurs. 25 L'invention concerne aussi un récepteur comportant au moins une antenne pour recevoir un signal émis par au moins un émetteur et un processeur configuré pour exécuter le procédé d'amélioration de calcul des symboles et leurs fiabilités, lesdits symboles étant utilisés pour la modulation d'un signal selon les étapes du procédé de l'invention. 16 Le récepteur peut comporter un démodulateur souple comprenant au moins les éléments suivants : - Un module de calcul des distributions a priori, - Un premier module d'estimations souples, ayant une sortie 5 reliée à un module de calcul de distributions a posteriori, - Un module de calcul des vraisemblances gaussiennes ayant une sortie reliée au module de calcul de distributions a posteriori, dont la sortie est reliée à un deuxième module d'estimations souples, - Le deuxième module d'estimations souples transmet ses io résultats directement à un module de division gaussienne et à un dispositif effectuant une moyenne dont la sortie est reliée au module de division gaussienne, - Le module de division gaussienne est relié à un lissage adaptatif fournissant les valeurs des symboles estimées ainsi que leurs 15 variances. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description d'exemples de réalisation donnés à titre illustratif et non limitatif annexée des figures qui représentent : - Figure 1, une vue du système mono-utilisateur avec émetteur et 20 récepteur utilisant une seule antenne, SISO, - Figure 2, un schéma-bloc avec une vue fonctionnelle d'un émetteur mono-antennaire, - Figure 3, un schéma de la structure d'un turbo DFE dans le cas SISO, pour un système SC-FDMA, 25 - Figure 4, une structure d'égaliseur linéaire à suppression d'interférence, - Figure 5, une structure du démodulateur souple, - Figure 6 et figure 7 des courbes de performances obtenues dans le cas SISO, 30 - Figure 8, un schéma généralisant le récepteur turbo DFE à des modulations quelconques, 17 - Figure 9, une vue du système multi-utilisateurs, - Figure 10, une vue d'un système multi-utilisateurs MIMO, - Figure 11, un schéma bloc fonctionnel pour un émetteur multiantennaire, 5 - Figure 12, un schéma de structure d'un turbo DFE dans le cas MU- MIMO, pour un système SC-FDMA, - Figure 13, un démodulateur souple pour le cas MU-MIMO, - Figure 14, le filtrage et le procédé de suppression des interférences dans le cas MU-MIMO, et 10 - Figure 15, un détail d'un bloc pour le démodulateur souple. La technique de turbo égalisation consiste en l'itération entre la fonction d'égalisation, de démodulation souple (ici appelé aussi demapper ou demapping souple) et de décodage, généralement avec la finalité d'améliorer le taux d'erreur bit (Bit Error Rate - BER) ou le taux d'erreur paquet (Packet 15 Error Rate - PER) tout en maitrisant la complexité du récepteur. Afin de bien faire comprendre l'invention, la première variante de réalisation est donnée, à titre d'exemple non limitatif, dans le cas d'un système simple entrée simple sortie ou SISO (Single Input Single Output) illustré à la figure 1. Le système est constitué d'un émetteur 10 et d'un 20 récepteur 20 tous les deux équipés d'une seule antenne à l'émission et à la réception. L'émetteur mono-antennaire générique 10 pour l'application SISO est représenté en figure 2. L'émetteur 10 prend en entrée des bits d'information, b, les bits sont codés dans un codeur 11 avec un codeur 25 correcteur d'erreur qui peut être un code convolutif, un turbo-code, ou un code à contrôle de parité à faible densité (Low-Density Parity-Check - LDPC) ou tout autre code pour lequel il existe un algorithme de décodage qui donne des informations souples. Les bits codés sont entrelacés avec un entrelaceur 12. Les bits entrelacés sont ensuite modulés par un modulateur 13. Le 30 modulateur 13 fournit en sortie des symboles tirés d'une constellation X. Les symboles modulés sont transmis à un bloc de mise en trame 14 qui organise 18 les données par blocs dans une trame et qui peut aussi insérer des séquences pilotes qui serviront, par exemple au récepteur, pour l'estimation du canal. Les séquences pilotes sont générées par un générateur de séquences pilotes 15. De plus, le bloc de mise en trame 14 implémente une 5 méthode de périodisation partielle des blocs de données qui permet, à la réception, d'implémenter un égaliseur dans le domaine fréquentiel. Le signal mis en trame est transmis vers une chaîne radio fréquence RF, 16, pour émission par une antenne Ae. Le bloc de mise en trame 14 peut utiliser une modulation à 10 multiplexage par répartition orthogonale de la fréquence (Orthogonal Frequency Division Multiplexing -OFDM) avec un total de N sous-porteuses dont M sous-porteuses utilisées avec un préfixe cyclique (Cyclic Prefix - CP) et éventuellement un suffixe cyclique (Cyclic Suffix - CS). Le bloc de mise en trame peut aussi mettre en place une modulation OFDM à accès multiple, 15 appelée dans la littérature OFDMA. Le bloc de mise en trame peut aussi mettre en place une modulation mono-porteuse à accès multiple par répartition en fréquence Single Carrier-Frequency Division Multiple Access (SC-FDMA), avec M le nombre de sous-porteuses utilisées pour le précodage avec une transformée de Fourier discrète (Discrete Fourier 20 Transform - DFT). Ici aussi le CP et éventuellement le CS sont utilisés. Dans ce cadre, si N = M, l'émetteur implémente un signal mono-porteuse (Single Carrier - SC). Les CP et CS peuvent être substitués par une séquence constante (par exemple des zéros, ou des séquences pilotes) ou évoluant 25 d'une trame à l'autre (par exemple par un procédé pseudo-aléatoire connu des utilisateurs), ce qui permet d'obtenir encore un signal avec la bonne propriété de périodicité partielle sur les blocs de données étendus avec l'une de ces séquences. Il est aussi possible de considérer l'utilisation d'une modulation monoporteuse avec une mise-en-forme spectrale (Spectrally 30 Shaped - Single Carrier - SS-SC) ou une modulation mono-porteuse avec une mise-en-forme spectrale à accès multiple par répartition en fréquence 19 (Spectrally Shaped - Single Carrier - Frequency Division Multiple Access - SS-SC-FDMA), où le signal pourrait être filtré en temps, ou en fréquence, par un filtre de mise-en-forme, après l'ajout des CP/CS. Sans sortir du cadre de l'invention, il est possible d'utiliser des 5 constellations ayant une propriété de non-circularité (appelée aussi impropriété dans la littérature [1]. La non-circularité s'exprime formellement par le fait que si x(n) est un symbole aléatoire de la constellation émis au temps n, alors l'espérance du symbole au carré est différente de zéro E[x2(n)] # O. La quantité E[x2(n)] est appelée aussi pseudo-covariance dans io la littérature. Cette propriété s'étend aux signaux échantillonnés ou continus aussi. L'invention s'applique donc à des émetteurs utilisant des constellations complexes comme la modulation d'amplitude en quadrature (Quadrature Amplitude Modulation - QAM), la modulation par changement de phase ou Phase Shift Keying (PSK), la modulation par changement de phase et 15 d'amplitude (Amplitude Phase Shift Keying - APSK) ; des constellations réelles, comme la modulation par changement de phase binaire ou « Binary Phase Shift Keying » (BPSK), ou les modulations en amplitude d'impulsion (Pulse Amplitude Modulation - PAM). La technique peut être appliquée aussi à des constellations tournées périodiquement comme la Tr/2-BPSK, où sur 20 les symboles pairs on utilise une constellation BPSK classique {+1, -1} et sur les symboles impairs on utilise une constellation tournée de -rr/2 radians {+j, - j}, et plus généralement les constellations Tr/M-MPSK, où la rotation périodique est de -rr/M radians. Il est aussi possible d'appliquer la technique décrite ici à des constellations appelées quasi-rectilignes, c'est-à-dire des 25 constellations dont les symboles peuvent être obtenus par filtrage complexe d'un signal décrit par les symboles d'une constellation réelle. Des exemples sont la modulation par déplacement de fréquence minimum (Minimum Shift Keying - MSK), la modulation par déplacement de fréquence minimum gaussien (Gaussian Minimum Shift Keying - GMSK), les modulations de 30 phase continues (Continuous Phase Modulation - CPM) avec alphabet 20 binaire ou la modulation d'amplitude en quadrature décalée en quadrature (Offset Quadrature Amplitude Modulation - OQAM). La figure 3 illustre un schéma fonctionnel du récepteur 20 dans le cas où le signal émis est obtenu avec une technique SC-FDMA avec préfixe 5 CP. On suppose ici que le signal de l'émetteur est synchronisé au récepteur avec une précision adéquate (i.e. inférieure à la durée du CP, si ce dernier est présent) et qu'un algorithme de synchronisation a fourni l'instant de synchronisation au récepteur. Le signal fourni en entrée au bloc 101 est donc déjà synchronisé en temps et fréquence et à la bonne fréquence 10 d'échantillonnage. Si dans la figure 3 on suppose qu'à l'émission un signal SS-SC-FDMA a été utilisé, et que le signal en entrée du bloc 101 a été aussi filtré par un filtre de mise en forme adapté à celui utilisé à l'émission (selon les méthodes connues par l'homme du métier), alors le récepteur peut être utilisé pour recevoir un signal de type SS-SC-FDMA. Dans la figure 3, si le 15 bloc de dé-allocation des sous-porteuses 103 coïncide avec l'identité (c'est- à-dire s'il est transparent) et M = N, ce récepteur peut être utilisé dans un système SC avec CP. Dans ce même cas (bloc 103 supprimé et M = N), si le signal émis est un signal SS-SC et si le signal en entrée du bloc 101 a été aussi filtré par un filtre de mise en forme adapté à celui utilisé à l'émission 20 (selon les méthodes connues par l'homme du métier), alors le récepteur peut être utilisé pour recevoir un signal de type SS-SC. Pour les deux derniers cas de signaux émis de type SC et SS-SC, si dans la figure 3 le bloc « suppression CP » est substitué par un bloc « Extraction du bloc de traitement » qui extrait les séquences connues (par exemple pour les utiliser 25 pour estimer le canal) et les blocs de données à traiter, alors le récepteur peut recevoir des signaux SC et SS-SC avec des séquences pilotes connues à la place des CP et CS. De plus, si dans la figure 3 les blocs 109 et 107 sont supprimés (ou équivalemment s'ils implémentent la transformation identité), alors le récepteur est adapté à recevoir un signal de type OFDM ou 30 OFDMA. 21 Le récepteur selon l'invention contient deux boucles de rétroaction. La première boucle est positionnée entre l'égaliseur DFE 108 et le démodulateur souple SISO 105. Dans la suite, le parcours de cette boucle est appelé « auto-itération » et un égaliseur qui utilise cette boucle est un 5 « égaliseur auto-itéré ». La deuxième boucle se situe entre le démodulateur souple SISO 105 et le décodeur bloc « décodage » 111. Dans la suite, le parcours de cette boucle est désigné sous l'expression « turbo-itération)> et un égaliseur qui utilise cette boucle est un « égaliseur turbo-itéré » ou un « turbo-égaliseur ». 10 Le récepteur 20 effectue un nombre 7' > 0 de turbo-itérations. Pour chaque turbo-itération donnée r = 0, le récepteur effectue un nombre ST > 0 d'auto-itérations. Il est possible de configurer le récepteur avec la suite (TA, soit au début de la réception du signal, soit dynamiquement pendant le traitement du signal ; les détails seront donnés 15 plus loin dans la description. Par exemple prenons (7' = 3,S0 = 5,S = 3,S2 = 1,S3 = 1), cela veut dire que le récepteur effectue cinq auto-itérations lorsque r = 0 (sans avoir encore décodé le signal), trois auto-itérations quand r = 1 c'est-à-dire à la première turbo itération (après le premier passage dans le décodeur), et une auto-itération quand = 2 ou 3, Le., à la 20 deuxième et troisième turbo itération. Les nombres d'auto-itérations ou de turbo-itérations règlent le compromis entre performances, complexité de calcul et mémoire nécessaire, au sein du récepteur. Dans la suite de la description, les paramètres en entrée ou en sortie du turbo égaliseur, du calcul égaliseur et les paramètres du 25 démodulateur souple, et du démodulateur souple, à l'exception des observations initiales et estimations initiales y, H, Ew, dépendent de l'identifiant de la turbo-itération courante = 0, ...,T et éventuellement de l'identifiant de l'auto-itération courante s = 0, ...,Sr. Par exemple, le filtre égaliseur es) indique le filtre utilisé à la turbo-itération -r et à l'auto-itération 30 S. Ceci-dit, dans la suite, pour simplifier la notation, l'apex ne sera pas 22 explicité, et nous nous référerons implicitement à la turbo-itération et auto-itération courante. De plus, il convient d'initialiser les blocs au début du traitement de réception et au début de chaque turbo-itération, donc pour chaque -r = 5 0, Dans cet exemple, la description est donnée pour un cas SCFDMA, en sachant que les adaptations pour les autres cas cités précédemment peuvent être déduites par un homme du métier. Après avoir extrait les blocs de données en provenance de io l'antenne de réception A,- du récepteur, à travers la suppression du CP dans le cas des systèmes SC avec CP ou du SC-FDMA, 101, les blocs de données passent à travers une transformation de Fourier rapide 102 (Fast Fourier Transform - FFT) de taille N (passage dans le domaine fréquentiel). Le signal en sortie de la FFT passe ensuite dans un bloc « De-alloc sous-p » 15 103 qui permet d'accéder aux M ressources fréquentielles (sous-porteuses) parmi N, à travers lesquelles le signal a été transmis. Ce bloc sélectionne seulement les sous-porteuses sur lesquelles le signal a été transmis (dé-allocation des sous-porteuses), et en sortie de ce dernier bloc de dé-allocation 103 le signal est représenté par un vecteur y de taille M, 20 regroupant les sous-porteuses utilisés. En parallèle, les séquences pilotes sont extraites des signaux provenant de l'antenne de réception Ar et envoyées dans le bloc « Estimation des canaux et de la variance du bruit » 104. Elles sont utilisées pour calculer une estimation de la réponse fréquentielle du canal entre l'émetteur et le 25 récepteur, sur les M sous-porteuses d'intérêts (celles utilisées par l'utilisateur pour envoyer l'information). La réponse fréquentielle du canal entre l'émetteur et le récepteur est organisée dans un vecteur h de taille M, ou de façon équivalente, dans une matrice diagonale H = diag(h) de taille M x M. Le bloc « Estimation des canaux et de la variance du bruit)> 104 fournit aussi 30 une estimation de la variance du bruit. Ce bruit contient la contribution du bruit thermique sur l'antenne de réception Ar et d'éventuelles interférences 23 dues à plusieurs effets parasites en provenance soit de traitements numériques, soit d'autres signaux extérieurs. Le bloc « Estimation des canaux et de la variance du bruit)> fournit la matrice de covariance du bruit dans le domaine fréquentiel Ew. Typiquement, on a une estimation de canal 5 et de covariance du bruit par bloc de données. Les traitements de la partie du récepteur où il y a la présence des deux boucles précitées sont alors activés. Au tout premier passage (1- = 0), OU à chaque début de turbo-itération (1- = 1, ...,7., si 3' > 0), le bloc 105, démodulateur souple SISO, est io activé pour une phase d'initialisation. Par convention cette phase est associée à l'indice d'auto-itération s = -1. Dans cette phase d'initialisation, le bloc 105 effectue plusieurs opérations : a) il vide ses registres internes ; b) il met à zéro les entrées du vecteur d'entrée X (signal égalisé) ; c) il pose o-, = +o. (ou autre valeur faisant office de l'infini). Puis le démodulateur 15 souple SISO 105 regarde s'il y a une entrée d'information a priori en provenance du décodeur et plus précisément du bloc entrelaceur 112. Ensuite, il génère une l'estimation souple X et sa fiabilité correspondante V. La sortie vers le bloc désentrelaceur 110, n'est pas calculée pendant la phase d'initialisation. A la toute première activation de la phase d'initialisation 20 (1- = 0) l'information a priori en provenance du bloc entrelaceur 112 correspond par convention à un vecteur de zéros. Ainsi, l'estimation souple X sera un vecteur nul et e = cri, où ()- est l'énergie moyenne des symboles des constellations d'origine de l'utilisateur. L'estimation du canal H et l'estimation de la matrice de covariance du bruit 25 sont ensuite transmises au bloc « calcul des paramètres de l'égaliseur » 106. Ce bloc prend aussi en entrée la quantité V en provenance du demapping souple qui mesure la fiabilité moyenne des symboles envoyés par l'utilisateur à l'intérieur du bloc de données qui est en cours de traitement (la moyenne est donc sur la longueur de bloc M dans le cas de l'exemple). Dans le cas 30 plus général, avec s> 0, le bloc de « calcul des paramètres de l'égaliseur » 24 106 fournit en sortie à chaque auto-itération les coefficients de l'égaliseur f de taille M dans le domaine fréquentiel, qui est en général un turbo-égaliseur linéaire MMSE auto-itéré avec retour EP. Dans tous les cas, le bloc de calcul 106 fournit en sortie aussi une estimation du bruit après égalisation pour 5 chaque bloc de données, calculée sous l'hypothèse de bruit gaussien et modèle additif X = x+ y. Ces deux sorties (filtre en domaine fréquentiel et variance du bruit résiduel) sont calculées de la façon suivante : M-1 1 VI ILlm I2 (1) 2, fm = = Imi=Odt, + 171L1m1 1 2 M = 0, , M - 1, (2) 2 (3) av - où est le biais de l'égaliseur fréquentiel MMSE et cr2m pour m = 0, M- 1 sont les covariances sur la diagonale de la matrice Ew (il n'est donc pas io nécessaire d'estimer les entrées de cette matrice hors de la diagonale principale), fm sont les coefficients du filtre f. Le vecteur de signal reçu y de taille M est passé dans le bloc « égaliseur linéaire à suppression d'interférence » 108 qui prend en entrée aussi les coefficients du filtre de l'égaliseur f, la réponse fréquentielle du is canal H et le vecteur de taille M d'estimations souples des symboles émis dans le domaine fréquentiel X, qui est généré par le bloc de FFT normalisée de taille M 107 à partir du vecteur de taille M d'estimations souples dans le domaine temporel X, en provenance du démodulateur souple SISO 105. La fiabilité globale de cette estimation est représentée par la valeur 17, plus elle 20 est proche de zéro plus la fiabilité est grande (i est une variance). Le bloc égaliseur linéaire à suppression d'interférence 108 réalise la suppression d'interférence (Interference Cancellation - IC) pour l'ISI, dans ce cas SISO. Le bloc égaliseur linéaire à suppression d'interférence 108 produit un vecteur 25 de taille M X qui représente une estimation des symboles dans le domaine fréquentiel. Le vecteur X est ensuite repassé dans le domaine temporel à travers une _ IFFT 109 normalisée de taille M, pour obtenir le vecteur X du signal temporel 5 égalisé, qui est envoyé au démodulateur souple SISO 105. A partir de X, de la variance du bruit après égalisation 4 et des informations souples a priori, par exemple sous la forme de métriques en rapport de vraisemblance logarithmique (Log-Likelihood Ratio - LLRs) La en provenance du décodeur, le démodulateur souple produit des informations souples pour chaque bit du io signal d'entrée, par exemple sous la forme de LLR extrinsèques Le. Ce démodulateur prend des formes différentes selon la statistique du signal après égalisation : si la constellation de départ est réelle on peut utiliser un démodulateur pour des statistiques gaussiennes complexes symétriques, autrement un démodulateur pour des statistiques gaussiennes avec pseudo- 15 covariance non nulle est plus adapté. Les métriques souples sont ensuite désentrelacées par le bloc désentrelaceur 110 qui est le bloc inverse du bloc entrelaceur 112. Ensuite, quand tous les bits du paquet sont récupérés, les métriques souples sont envoyées au décodeur 111 qui génère, produit des estimations de 20 probabilités des bits d'information envoyés (par exemple des probabilités a posteriori). Ces estimations de probabilités peuvent être utilisées pour obtenir une estimation dur b' des bits transmis b avec un détecteur à seuil. Si le récepteur utilise des turbo-itérations, le décodeur produit aussi des estimations des bits codés envoyés, c'est-à-dire des probabilités 25 extrinsèques (EXT), par exemple sous forme de LLRs, qui mesurent la probabilité que les bits codés envoyés soient 0 ou 1. Les EXT sont ensuite envoyées à l'entrelaceur 112 pour être entrelacées. L'entrelaceur 112 situé dans le récepteur 20 opère de la même manière que l'entrelaceur 12 se trouvant dans l'émetteur 10, avec la seule différence que ce dernier opère 30 sur des données binaires, alors que l'entrelaceur 112 opère sur des LLRs. Les EXT entrelacées deviennent des informations a priori du point de vue du 26 reste du récepteur, par exemple sous forme de LLRs a priori La et entrent dans le démodulateur souple SISO 105. En fonction de l'ordonnancement de l'activation du démodulateur souple SISO, de l'égaliseur et du décodeur, différentes boucles peuvent être 5 réalisées. Le bloc « égaliseur linéaire à suppression d'interférence » 108 réalise la suppression d'interférence et applique une égalisation linéaire en fréquence sur le signal d'entrée. Dans la littérature il est appelé aussi égaliseur linéaire à suppression d'interférence (Linear Equalizer - 10 lnterference Cancellation - LE-IC) puisqu'il régénère une estimation du signal interférant pour le soustraire. Cette structure est déjà bien connue dans la littérature [2]. La différence par rapport à l'état de l'art existant est dans la valeur de l'entrée X (estimation des symboles dans le domaine _ fréquentiel) fournie au bloc 108, ainsi que le filtre utilisé, ces deux quantités 15 étant calculées de façon nouvelle. La sortie X du bloc égaliseur linéaire à suppression d'interférences 108 sera donc différente que celle qui peut être trouvée dans [2] à égalité des autres hypothèses. La figure 4 donne un exemple détaillé pour le filtrage et le procédé de suppression de l'interférence, qui est réalisé de façon efficace comme il 20 suit. Le bloc « égaliseur linéaire à suppression d'interférence » 108 comprend un soustracteur 201 qui agit sur des vecteurs de taille M et qui retranche du signal réellement reçu y, une estimation de ce même signal reçu (signal utile avec interférences du canal de propagation). Cela permet d'obtenir un signal correctif introduit dans un filtre linéaire 202 de type MMSE 25 de réponse fréquentielle f, calculée par le bloc de calcul des paramètres de l'égaliseur 106. Le filtre linéaire 202 est suivi lui-même d'un sommateur 203 de vecteurs de taille M, qui réintroduit dans le signal correctif filtré la contribution du signal utile déjà estimé. Le bloc de retranchement du signal reçu estimé 201 obtient l'estimation à retrancher par un module de 30 génération du signal reçu estimé (signal utile avec interférences du canal de propagation) 200. 27 Le bloc 200 de génération du signal reçu estimé produit une estimation du signal reçu (signal utile avec interférences du canal de propagation) à partir de l'estimation des symboles dans le domaine fréquentiel X, et l'estimation du canal H. Ceci est obtenu par simple 5 multiplication terme à terme du vecteur X avec l'estimation de la réponse fréquentielle du canal h, tel que h = diag(11). Donc, le signal reçu estimé en fréquence est h C) X ou équivalemment HX en notations matricielles. Quand = S = 0, X est un vecteur de zéros (pas d'estimation disponible). Le sommateur 201 implémente un retranchement dans le domaine fréquentiel du signal reçu estimé au signal reçu, y - h O X. Le résultat est nommé signal correctif. Le bloc 202 applique simplement le filtre f. Il s'agit là encore tout simplement d'une multiplication entrée par entrée du vecteur du filtre de taille M et du vecteur d'entrée de même taille. Cette étape correspond à l'égalisation 15 proprement dite qui réduit d'avantage l'ISI résiduel du signal correctif. Enfin, le sommateur 203 génère une estimation fréquentielle des symboles égalisés en ajoutant au signal correctif filtré l'estimation du signal utile obtenue à partir de l'information souple en provenance du démodulateur souple SISO à l'auto-itération précédente et si s = 0, à la turbo-itération 20 précédente. Ceci est implémenté de manière à améliorer progressivement (itération après itération) l'estimation du signal utile final. La sortie X s'écrit De façon équivalente, cette relation aurait aussi pu s'exprimer par des produits matriciels avec X = ï + diag(f*)(y - HX). _ _ La structure du démodulateur souple SISO 105 est illustrée à la 25 figure 5. Tout d'abord, l'utilisateur choisit les symboles à envoyer parmi une constellation (potentiellement multidimensionnelle, ici l'invention est présentée avec une constellation avec une dimension complexe) X c C, où C est l'espace des nombre complexes, de cardinalité Q = III, tel que 28 q = log2 Q. Nous supposons que les symboles de la constellation ont une moyenne nulle et une puissance statistique égale à cri. Il est aussi possible de normaliser une constellation de façon que 4 = 1. Le modulateur dans l'émetteur 10, à partir de chaque vecteur de q bits codés et entrelacés en s sortie du bloc entrelaceur 12, génère le symbole correspondant de la constellation selon une fonction M, qui est aussi appelée fonction d'étiquetage ou simplement étiquetage, notion connue de l'homme du métier. En se focalisant sur un vecteur de données x d'émission de longueur M, ses entrées xm, m = 0, ...,M - 1 sont obtenues ainsi : xm = .7e(dia), où dm = 10 [drag, ... , d(n+i)a_i] est le m-ième groupe de q éléments du vecteur d qui est le vecteur des bits codés et entrelacés de taille qM correspondant au bloc de données. Le vecteur de bits dm peut être interprété comme l'étiquette binaire du symbole xm. On note traditionnellement aussi dmq+i = .74T1(xm) le j-eme bit (de l'étiquette) du symbole xm, avec j = 0, ..., g - 1. Pour terminer, on 15 pose aussi if = ta E Xlittri-1(a) = b), i.e. l'ensemble de tous les symboles de la constellation qui ont un bit de valeur b (0 ou 1) dans la j-eme position de leur étiquette binaire. Ces notations sont nécessaires pour expliciter facilement les calculs des probabilités des bits et des symboles dans la suite. La figure 5 décrit un exemple détaillé pour la structure du 20 démodulateur souple. L(d) représente le vecteur des LLRs a priori sur les bits codés et entrelacés, en provenance du décodeur et qui sont associés au vecteur de signal estimé après égalisation X(T,$) en train d'être traité. A la turbo-itération I" = 0 (i.e. le décodage n'a pas encore eu lieu), La(d) est un vecteur de zéros (aucune information a priori sur les bits d disponible). 25 A la toute première utilisation du récepteur (1" = 0), ou à chaque début de turbo-itération (1" = 1, ...,T, si / > 0), après avoir activé le décodeur 111 et l'entrelaceur 112, le démodulateur souple SISO 105, est activé pour une phase d'initialisation, avec un indice d'auto-itération conventionnel de s = -1. Dans cette phase, le démodulateur souple SISO 105 exécute 30 plusieurs étapes préparatoires : a) il vide ses registres internes en 29 correspondance du bloc 307; b) il met à zéro les entrées du vecteur d'entrée X (signal égalisé) ; c) il pose 4 = +oe (ou autre valeur faisant office de l'infini). Ensuite, la chaîne de calcul constituée par un bloc de calcul des distributions a priori 300, un bloc d'estimations souples 301, un estimateur de 5 moyenne effectuant une moyenne 302 et un interrupteur 303 est activée. Cette chaîne de calcul n'est plus activée pour les autres indices s > 0, mais uniquement quand une nouvelle initialisation est instruite. Pendant la phase d'initialisation, en fonctionnement normal les autres blocs ne sont pas activés. 10 Lors de la phase d'initialisation, le bloc de calcul des distributions a priori 300 effectue le calcul des probabilités de masse (distributions de probabilité) a priori sur les symboles de la constellation X (ce bloc est connu de l'état de l'art) à partir des La(d). Pour chaque symbole xm, m = 0, ...,M - 1 de chaque bloc de données utilisé pour transmettre le mot de code, on 15 calcule une distribution selon les formules suivantes : ge (a) = riri=-01 exp(-.74)-1 (a) IC) (d qm+ i)) (5) Ise (a) = nel' je) (a) , avec nem, = EaExie (a) (6) Il est possible d'implémenter ces calculs dans le domaine logarithmique (selon des méthodes existantes dans la littérature) pour des questions de réduction de complexité de calcul. Le bloc d'estimations souples 301 calcule une estimation des 20 symboles émis, ainsi qu'une estimation de leur fiabilité, en partant des distributions de probabilité en entrée. Puisque l'ensemble des distributions a priori lipm(t)(a)} (pour tous les symboles de tous les blocs de données composant un même mot de code) change seulement à chaque turbo-itération, le bloc 301 est appliqué aux {9(a)} une seule fois par turbo- 25 itération et produit les estimations suivantes : Xm = E (ro) l,(,)[xml; vm = V ar yeo[xm] (7) 30 avec les définitions standard d'espérance El,(0[.] et de variance V ar 1,(0[]. Quand 1- = 0, i.e., il n'y a pas de turbo-itération et donc pas d'information a priori disponible en provenance du décodeur, 471,(°'°) = 0 et vm(") = cr pour tous les m. Autrement dit, les symboles envoyés ne peuvent pas être estimés 5 et ils ont une variance égale à la puissance statistique de la constellation d'origine (c'est-à-dire, l'imprécision est maximale). Cette information pour r = 0, peut être exploitée pour éviter les calculs du bloc de calcul de distributions a priori 300 et le calcul de la moyenne des variances vmerm), et réduire ainsi la complexité de la solution. 10 Le bloc de moyenne 302 calcule une simple moyenne des valeurs en entrée. Ainsi, quand le bloc moyenneur est appliqué aux vm("), m = 0, ...,M - 1 d'un bloc donné, il calcule tout simplement 1)-(r,o) = 114-1 Emml vin(T,O) et il génère des moyennes, bloc par bloc. Il est aussi possible, dans une autre implémentation, de considérer tous les symboles 15 correspondants à un mot de code et d'effectuer la moyenne sur cette longueur. f(") est la variance moyenne a priori de l'estimation des symboles émis (soit sur un bloc de données, soit sur tout le mot de code). Elle permet à la fois d'obtenir un estimateur plus robuste que les variances instantanées de la fiabilité moyenne des symboles estimés et aussi de pouvoir 20 implémenter une stratégie d'égalisation dans le domaine fréquentiel, en général moins coûteuse d'un point de vue calculs. Les estimations a priori des symboles émis .5e0) ainsi que leur variance moyenne 1-7(T,°) sont acheminées vers l'interrupteur 303. Si l'auto-itération courante est s = -1 (i.e. phase d'initialisation), l'interrupteur 303 25 fournit en sortie es+1) = 4.,(T'°) pour tout m et 17(1,s') = f(T,°), c'est-à-dire les estimations a priori des symboles et leur variance moyenne a priori. Si s > 0, l'interrupteur 303 fournit en sortie les estimations des symboles 5es+1) et leur variance moyenne 17(t's+1) en provenance du bloc de lissage adaptatif 307. Lorsque s = -1, pour toute turbo-itération, la sortie de l'interrupteur 303 31 correspond donc à la sortie d'un mapping souple avec information souple en provenance du décodeur, tel que décrit dans la référence [2]. La sortie 1, e(r) (d) n'est pas calculée en phase d'initialisation. Une fois que les sorties de l'interrupteur 303 sont calculées, on considère que l'indice s d'entrée est 5 incrémenté de 1. Nous décrivons par la suite la chaîne de calcul quand le bloc 105 est activé avec un indice d'auto-itération s = 0, ...,S - 1. Dans ce cas, il n'est pas nécessaire de parcourir la chaîne de calcul constituée par les blocs 300, 301, 302. L'interrupteur 303 laissera passer les sorties du lissage adaptatif 10 307. La figure 5 montre qu'il n'est pas nécessaire non plus de parcourir la chaîne comprenant le calcul LLRs a posteriori 309, le sommateur 310 pour calculer les LLRs extrinsèques de sortie. En effet, l'interrupteur 309 est fermé seulement lorsque s = Sr, c'est-à-dire quand les auto-itérations sont 15 terminées. Selon une variante, il est possible de calculer les LLRs a posteriori pour chaque indice d'auto-itération, pour pouvoir lancer le décodeur. Nous décrivons par la suite la chaîne de traitement constituée par les blocs 304, 305, 301, 302, 306, 307 et 303. Le bloc de calcul des vraisemblances 304 calcule des 20 vraisemblances à partir des estimations Jemer's) (les M symboles égalisés du bloc courant) et de la variance du bruit résiduel après égalisation crv2Cr's), en utilisant un modèle gaussien non-biaisé pour les estimations -.41..$). A noter qu'il est décrit dans l'état de l'art comment utiliser un modèle gaussien biaisé. Ce calcul donne : (r,$)(,.X(r,$) (_ les) _ al2/0.12(r,$)) , Lm m la) = exp m = 0, ...,M - 1, 25 pour chaque a appartenant à la constellation X. Il est possible d'implémenter ce calcul dans le domaine logarithmique (selon des méthodes existantes dans la littérature) pour des questions de réduction de la complexité de calcul. 32 Le bloc de calcul de distributions de probabilité a posteriori 305 calcule les distributions des symboles émis en connaissant l'observation (symboles égalisés), selon une méthode connue de l'homme du métier: dg-,$) (a) ei,$) (je ,$) la )3,41r) ( a ) 73^'s) (a) = n el;rri) (a), avec nd = EaEx ce (a) Il est possible d'implémenter ces calculs dans le domaine logarithmique (selon des méthodes existantes dans la littérature) pour des questions de réduction de complexité computationnelle. L'ensemble 17)'s) (a)} des distributions a posteriori pour les symboles «EX, du bloc de données courant est ensuite traité par le bloc d'estimations souples 301 pour générer des estimations a posteriori des symboles émis it = ED(r,$)[xm] et leur correspondante variance a posteriori ym(t's) = Varip,$)[xm]. La nouvelle estimation des symboles émis, générée en combinant l'information a priori disponible avec l'information contenue dans l'observation selon l'invention présente des caractéristiques supérieures à l'estimation a priori, ou à une estimation dérivée de la seule observation. 15 Dans le cadre de la théorie d' « expectation propagation » (EP), cette opération permet aussi de projeter une description probabiliste mixte sur les symboles envoyés (des variables gaussiennes pour les observations, combinées avec des estimations a priori qui sont des variables catégoriques) sur une description probabiliste unique, selon un modèle gaussien tel que 20 décrit dans la référence [3]. Plus dans le détail, nous utilisons ici un résultat de la théorie de l'information, connu sous le nom de « reverse 1-projection », ou « M-projection », qui stipule que le fait de trouver une distribution de type exponentiel (dans notre cas, gaussienne) qui minimise la divergence de Kullback-Leibler (une mesure de distance probabiliste) par rapport à une 25 distribution quelconque, revient à trouver la distribution exponentielle ayant les mêmes moments statistiques (« moment matching »). 33 Les variances a posteriori sont ensuite moyennées par le bloc moyenneur 302 pour obtenir une variance a posteriori moyenne fCe,$) = M-1 EMmIP yr(r.T.'s) sur le bloc de données en cours de traitement. Comme dans le cas des estimations a priori, la variance moyenne a posteriori permet à la fois d'obtenir un estimateur plus robuste de la fiabilité moyenne des symboles émis estimés du bloc de données et aussi de pouvoir implémenter une stratégie d'égalisation dans le domaine fréquentiel, en général moins coûteuse d'un point de vue des calculs. Les estimations a posteriori min(T's) des symboles émis et leurs variances io moyennes a posteriori es) sont ensuite envoyées dans le bloc de « division gaussienne » 306 qui prend aussi en entrée les observations (symboles égalisé) 277,(T's) et la variance du bruit résiduel après égalisation o-v2er's). En suivant le cadre de la théorie « expectation propagation », ce bloc va générer une information extrinsèque sur les symboles émis, en suivant un modèle 15 gaussien. Les solutions de l'état de l'art utilisent soit une information a priori (catégorique, extrinsèque du décodeur), soit une information a posteriori (catégorique aussi, combinée de l'information extrinsèque du décodeur et de la vraisemblance du démodulateur souple). Cette opération permet donc de mieux respecter le principe turbo en retranchant à l'information a posteriori, 20 l'information a priori disponible au démodulateur, au début du traitement (c'est l'extrinsèque de l'égaliseur) et en permettant ainsi d'améliorer la qualité du signal de retour et, in fine, les performances du récepteur. La division gaussienne 306 va calculer des estimations extrinsèques des symboles émis -,(r,s+1) X m et une variance extrinsèque moyenne -17'er's+1) selon la procédure 25 suivante : (r a ,$) 2(r,$) Y 37 ' 34 { _(T ' s) --er,s+1) 2(T,$) s) V av Y ' 2(r,$) Cr > Y ' 2er,$) ' av i/ Y ' En effet, le cas o-v2(') 0, T. Un exemple d'implémentation est de choisir une décroissance exponentielle Au = cr(r,$)e-Kr,$), où A(.) sera une fonction croissante avec s et/ou T, et a(.) sert de pondération pour s et petits. Dans cet exemple, au fur et à mesure des itérations, nous donnons plus de confiance aux estimations courantes. En gardant la même expression, on io pourrait imaginer f3 comme une fonction nulle, une fonction a(.) qui reste constante, proche de un, pendant les quelques itérations initiales, puis qui décroît linéairement après avoir décidé que les points fixes erronés ne seront pas pris en compte. Encore une autre méthode de lissage adaptatif possible et de pondérer les précisions, c'est-à-dire l'inverse des variances, à la place 15 que les variances elles-mêmes. Par exemple, la formule (14) peut être substituée par : f(T,s1-1) = 1 \ 1 - (1-Au je(T,s+ 1.) +PT,S(r,$) De plus, la formule (13) peut être substituée par : 1.771(r,s+i) = yr,si,, m(r,s+1.) p 1-11-,sA'm OÙ Y r,s r,s n'est pas forcement 1. Par exemple on pourrait choisir Au= ef(r,s+1)/77(T,$) et yr,s = (1 Oi(T,51-1)/et(T,S-1-1). est possible d'implémenter 20 des filtres adaptatifs, par exemple de type Kalman. En général on peut envisager des formules de lissages aussi non-linéaires eCr's+1) = nr,$) (t--(T,$) r,s +1)) et .5-e7;.,r's+1) = ((t,$) ïr(.tir,s -F1) ,e(r,$), f(r,s+1), f(T,s+1). ) où la fonction de lissage sur les fiabilités /vec's) est fonction seulement des fiabilités courantes et passées pour réduire la complexité du calcul. 25 Au final, selon l'une des méthodes décrites précédemment ou évidentes pour l'homme de l'art, le bloc de lissage adaptatif 307 génère des estimations extrinsèques lissées i'm(T's+1) des symboles émis, ainsi que pour chaque 36 vecteur leur variance moyenne extrinsèque lissée f(r,s") qui est une mesure de la fiabilité du vecteur de ces estimations. Ces quantités sont transmises à l'interrupteur 303, elles sont considérées comme des sorties du bloc 105, si s = 0, ..., - 1. 5 Nous décrivons maintenant le fonctionnement du démodulateur souple SISO 105 quand s = S. Une valeur de s = ST indique la dernière auto-itération. Après avoir exécutée la démodulation souple SISO 105, le traitement de décodage sera enclenché. Dans ce cas, le démodulateur souple SISO 105 produit en sortie, pour chaque nouvelle turbo-itération avec le décodeur, des LLRs extrinsèques L(d) sur les bits codés du bloc de données traité. Dans cette phase, le bloc interrupteur 308 est fermé et seuls le bloc 309 et le sommateur 310 sont utilisés. Il n'est pas nécessaire de calculer les autres sorties du bloc 105. Par un procédé appelé dans la littérature « marginalisation » [4], les Le(r)(d) 15 sont générées à partir des LLRs a posteriori calculés par le bloc 309 et en soustrayant les LLRs a priori L(d) à travers le bloc somme 310. Le bloc de calcul des LLRs a posteriori 309 est un bloc connu de la littérature : Ecrei? es)(a) L(r)(dmq+ j) = in (r,$) m = 0, ...,M -1, j = 0, ...,q -1 (17) Ectee Dm (a) D'autres bases de logarithme et des méthodes de simplification des calculs sont aussi applicables. Ensuite les LLRs extrinsèques sont calculés par le 20 bloc somme 310 comme suit : L(er)(dmq+i) = L(r)(dmq+j) - L(j)(dmq+i). Dans des implémentations sur des programmes d'ordinateur ou sur du matériel ,il convient aussi de maitriser la dynamique des LLRs du demapping souple et du décodeur afin de stabiliser, rendre plus robuste le comportement de l'algorithme, et de réduire son occupation mémoire tout en 25 gardant des performances acceptables. Différentes stratégies telles que décrites dans les références [5], [6], [7] sont possibles et peuvent être appliquées indépendamment du procédé décrit ci-dessus. 37 Les figures 6 et 7 illustrent un exemple de performance de l'invention dans un cas de système SC avec synchronisation et estimation de canal parfaites, dans le cadre SISO. La constellation utilisée par l'émetteur est une 8-PSK, le canal de propagation est le canal Proakis C (avec 5 puissance moyenne normalisé à un) avec un profil de puissance de type triangulaire (1, 2, 3, 2, 1). La figure 6 montre comment différentes instanciations de l'invention, dénommées FD S-SILE, avec ST = S pour tout T, et avec S = 1, 2, 5 produisent des meilleures performances que l'égaliseur MMSE linéaire standard et l'égaliseur turbo linéaire MMSE fréquentiel avec 10 IC de [1], dénommé dans les figures « FD ». La figure 6 en particulier montre les performances avec 7' = 0, sans turbo-itération, et avec 3" = 1, en turboitérant une fois avec de décodeur. Dans ces deux cas les gains sont considérables. Dans les mêmes hypothèses, la figure 7 illustre les performances 15 avec T = 10 turbo-itérations. Les propositions considérées arrivent à mieux approcher la borne du filtre adapté que le turbo-égaliseur MMSE linéaire fréquentiel avec IC. Sans sortir du cadre de l'invention, le procédé peut s'appliquer dans différentes structures d'égaliseur, dont quelques exemples sont donnés 20 ci-après à titre illustratif et nullement limitatif. L'invention peut s'appliquer à un égaliseur turbo fractionnairé voir par exemple [8]. Dans ce contexte, l'unique différence est que le signal en provenance de l'antenne est sur-échantillonné avec un facteur avec un facteur o > 1, typiquement entier. En faisant référence à la figure 3, les blocs 25 de données auront donc une taille oN dans la FFT 102, oM dans la transformée de fourrier inverse IFFT 109 et dans la FFT 107. De plus, le bloc IFFT 109 est suivi par un échantillonneur de débit o et le bloc FFT 107 est précédé par un interpolateur idéal de facteur o (qui introduit o - 1 zéro entre deux échantillons entrants). Le bloc d' « estimation de canal et variance du 30 bruit » doit effectuer ses estimations dans le domaine interpolé en incluant les filtres éventuels d'émission et réception, et puis passer éventuellement 38 les estimations en fréquence. Le calcul des paramètres de l'égaliseur va effectuer les calculs du filtre dans le domaine sur-échantillonné, en revanche il va estimer la variance du bruit résiduel d'égalisation au temps symbole (pour la passer au démodulateur souple). Le bloc turbo égalisation 108 va 5 donc travailler dans le domaine sur-échantillonné. En revanche, le bloc de démodulation souple SISO 105 conserve la description précédente. L'invention peut s'appliquer à un égaliseur de type Overlap FDE [9]. Dans ce type de récepteur, le bloc de données entre deux séquences pilotes est plus long du bloc de traitement, et pour des questions de io complexité ou de nécessité (par exemple un canal variable dans le temps qui change dans l'horizon d'un bloc de données), le traitement d'égalisation fréquentiel est décomposé sur plusieurs segments (avec recouvrement) du bloc de données initial. Ces segments souffrent d'une interférence entre eux. Cela implique que la matrice de covariance du bruit en sortie du bloc 15 d' « estimation de canal et de variance du bruit » 104 inclut aussi la variance due à l'interférence, ce qui modifie les valeurs du filtre d'égalisation. L'invention peut être appliquée à la fois dans des stratégies de réjection d'interférence [10] et de suppression d'interférence [11]. Dans ce dernier cas, la variance 17(T,5) en sortie du démodulateur souple SISO 105 20 peut aussi être utilisée dans le bloc 104 (lien en pointillé dans la figure 3) pour fournir une meilleure estimation de la variance du bruit plus interférence résiduelle par exemple selon la formule suivante : Ew = Ew, + 2.17(T,5)FN, Gp F1,4p où Np est la taille du segment, FNp est la matrice de transformation de Fourier discrète de taille N' x N', Ew, est la matrice de covariance du bruit 25 avant ajout du terme relatif à l'interférence entre segments, Gp est une matrice N, x Np avec les L - 1 diagonales supérieures non-nulles et égales à celles d'une matrice circulante générée à partir de l'estimation de la réponse impulsionnelle du canal hp = [ho,p, hl_1317. obtenue en correspondance du segment p de données. De plus, si l'estimation souple du signal envoyé et sa 39 variance correspondante en sortie du démodulateur souple SISO 105 sont fournis au bloc d'estimation de canal et de variance de bruit 104 (liens en pointillé dans la figure 3), ce dernier bloc peut implémenter des procédés d'estimation de canal et de variance de bruit en exploitant ces informations s additionnelles. Les nouvelles estimations de canal et de variance de bruit (possiblement à chaque itération) peuvent être exploitées par le reste du récepteur selon la description précédente. La technique proposée peut s'appliquer aussi à des égaliseurs linéaires au sens large. 10 De manière générale, l'invention peut être utilisée au sein de tout type de récepteur itératif étant structuré comme illustré sur la Figure 8. En effet, on considère disposer d'un émetteur 10 utilisant une modulation 13, une mise en trame 14 et une chaîne RF 15, quelconques, sur un bloc de données x pour le transmettre à travers un canal de propagation 15 quelconque. On considère disposer d'un frontal RF 500 qui extrait l'observation y du signal, reçu à l'antenne, perturbé par le canal de communication, et qui le fournit au bloc d'estimation canal 501. Cet estimateur extrait du signal y, les séquences pilotes, afin d'utiliser une des techniques d'estimations, connues de l'homme de l'art, pour obtenir des 20 paramètres caractérisant le canal. On suppose disposer d'un bloc « récepteur adaptatif à suppression d'interférence » 502, capable de fournir une estimation R = x + v modélisée selon un modèle statistique gaussien additif (v étant le bruit gaussien équivalent ajouté au signal utile). L'estimation R est obtenue grâce à un procédé de suppression d'interférence 25 qui, en utilisant une estimation préalable ,- du signal transmis x, et une estimation des paramètres du canal, génère une estimation de l'interférence et la soustrait au signal reçu. Par ailleurs, on dispose d'un bloc de calcul des paramètres du récepteur 503, qui utilise les estimations du canal et la fiabilité de l'estimation 30 préalable i7, pour calculer les paramètres intrinsèques du récepteur 502, et pour calculer la variance a caractérisant la puissance d'interférence résiduelle y, en sortie du récepteur 502. L'invention s'applique au niveau du démodulateur souple 504, où, pour un bloc de données avec des statistiques supposées blanches, la 5 technique de propagation d'espérance est utilisée pour calculer les estimations souples 5-c qui serviront comme retour de décision pour le récepteur, avec la variance ÇT caractérisant la fiabilité des ces mêmes estimations. Suite à un nombre approprié d'itérations S, le démodulateur souple 504 calcule des LLRs binaires extrinsèques sur les bits déterminant io les symboles qui peuvent servir à estimer les bits transmis à travers une décision dure 505, consistant à comparer chaque LLR binaire à zéro. Si les bits transmis ont été traités au préalable par un encodeur canal (comme par exemple dans l'émetteur 10 détaillé à la figure 2), les modifications nécessaires au récepteur sont représentées en tirets sur la is figure 8. Le bloc décodeur binaire à entrées/sorties souples 507 est utilisé conjointement avec le démodulateur souple 504, et le récepteur adaptatif 502, selon un ordonnancement choisi de manière appropriée, pour aider au décodage des bits d'information de l'émetteur. En particulier, on pourrait utiliser l'ordonnancement à deux boucles, décrite plus haut pour l'application zo SC-FDMA. Si les mots de code en sortie de l'encodeur canal ont été entrelacés par un entrelaceur binaire (voir par exemple émetteur 10 à la figure 2), alors les blocs desentrelaceur 506 et entrelaceur 508 sont inclus dans boucle d'itération entre le démodulateur souple 504 et le décodeur 507. Finalement, en considérant un récepteur itératif qui améliore aussi son 25 estimation du canal au fur et à mesure des itérations, en plus du signal transmis, on peut considérer que les estimations 5-c et la fiabilité en sortie du démodulateur souple 504 sont dirigées aussi vers l'estimateur du canal 501 (flèches en pointillés sur la figure 8). Dans ce cas, l'estimateur 501 peut utiliser, à partir de la première auto-itération du récepteur, l'estimation du 30 signal transmise, conjointement avec l'observation du canal y, pour raffiner les paramètres du canal estimé, par exemple avec une méthode de type 41 « moindres carrés » sur la totalité du bloc de données, au lieu d'utiliser seulement les pilotes. En effet, l'invention porte, notamment, sur la technique de l'estimation du signal transmis, au niveau du démodulateur souple 504 qui s travaille sur un modèle équivalent à un bruit additif gaussien. Elle peut être utilisée au sein de tout récepteur itératif qui passe par une telle étape de démodulation et se distingue des démodulateurs durs, utilisés entre autres dans les égaliseurs à retour de décision conventionnels, qui sont sujets à des problèmes de propagation de l'erreur. La technique considérée ici est io différente aussi des techniques de démodulation souple utilisant une estimation a posteriori du signal, qui n'enlèvent pas de l'estimation du signal de retour l'information déjà connue par le récepteur, induisant ainsi le récepteur à l'erreur en le biaisant avec ses propres informations. Comparée à d'autres techniques d'estimation avec propagation d'espérance utilisées 15 dans les récepteurs, le procédé selon l'invention fournit une fiabilité moyenne par bloc de données, ce qui permet de rendre le récepteur moins complexe et plus robuste face à des erreurs ponctuelles d'estimation. Par ailleurs, l'invention applique un lissage adaptatif sur le signal estimé, afin de fournir un degré de liberté additionnel qui permet de calibrer le rapport entre la 20 robustesse des performances et la vitesse de convergence du récepteur vers les performances limites. Par simplicité, les cas d'utilisations décrits ci-dessus, s'intéressent au traitement d'un signal mono-antennaire provenant d'un seul utilisateur. Nous allons montrer par un exemple d'usage que l'invention peut être utilisée 25 au sein d'un récepteur multiutilisateur itératif. L'invention s'applique aussi au cas mono-utilisateur single input multiple output (SIMO), c'est-à-dire à un système de communications avec un seul émetteur et un seul récepteur utilisant R antennes à la réception. L'invention s'applique aussi au cas mono-utilisateur multiple input 30 multiple output (MIMO) où le seul émetteur est doté de plusieurs, Tu, antennes de transmission. 42 L'invention s'applique aussi à un système de transmission sans fil multi-utilisateurs, où les utilisateurs u (101,..10u) sont dotés d'une antenne (Tu = 1) et le récepteur 20 est doté de R antennes où R peut être égal à un ou plus. Nous supposons que U utilisateurs transmettent vers le récepteur et 5 utilisent les mêmes ressources temporelles et fréquentielles. Le système est représenté en figure 9. Ce type de système est appelé aussi système Multiple Input Multiple Output (MIMO) distribué ou MIMO virtuel, ou simplement MIMO. Si ces utilisateurs utilisent des ressources temporelles et fréquentielles distinctes, permettant au récepteur d'apercevoir leurs signaux io sans interférence (i.e. orthogonaux), le récepteur peut traiter chacun des utilisateurs comme un système SIMO/SISO mono-utilisateur indépendant. L'invention s'applique aussi à un système de transmission sans fil multiutilisateur, où les utilisateurs u (101,..10,) sont dotés d'une ou plusieurs antennes (Tu k 1) et le récepteur 20 est doté de R antennes ou R 1. Nous 15 supposons que U utilisateurs transmettent vers le récepteur et utilisent les mêmes ressources temporelles et fréquentielles. Le système est représenté en figure 10. Ce type de système est appelé aussi système Multi-User MIMO (MU-MIMO). Si ces utilisateurs utilisent des ressources temporelles et fréquentielles distinctes, permettant au récepteur d'apercevoir leurs signaux 20 sans interférence (i.e. orthogonaux), le récepteur peut traiter chacun des utilisateurs comme un système MIMO/MISO mono-utilisateur indépendant. Dans ce cas, seule l'interférence entre les signaux des antennes d'un même émetteur sera présente ou non-négligeable. Nous décrivons l'invention d'abord dans le cas simple d'un 25 récepteur multi-antennaire dans le cadre d'un système de type MU-MIMO. Nous notons que les autres scénarios ne sont que des cas particuliers de ce dernier et peuvent s'en déduire. La figure 11 illustre le cas d'un émetteur multi-antenne générique 10u, et le schéma bloc fonctionnel associé. Le codeur correcteur d'erreur 11u 30 utilisé dans ce cas, n'est pas sujet à des contraintes différentes de celles d'un émetteur mono-antennaire, il prend en entrée des bits d'information, et 43 fournit des bits encodées en sortie. L'entrelaceur spatio-temporel 12u fournit Tu flux de bits entrelacés indépendants en sortie. Chacun des flux de bits entrelacés, est ensuite modulé par un modulateur 13u,1 , 13,0-, dédié à chaque antenne qui fournit en sortie des symboles tirés d'une constellation 5 XU,t qui peut être différente d'une antenne à l'autre et d'un utilisateur à l'autre. Les remarques sur les constellations auxquelles l'invention reste applicable s'étend pour le cas de chacun de T, modulateurs d'un émetteur multi-antenne. Les blocs de mise en trame 14 organisent les symboles en utilisant la même stratégie pour chacun des antennes. Le bloc 15u io générateur de séquences pilotes est étendu au cas multi-antenne, avec une capacité de fournir des séquences pilotes distinctes pour les différentes antennes. Les exemples de modulations données, tel qu'OFDM, (SS)-SCFDMA ou (SS)-SC, restent valables. Ensuite chacun des signaux préparés pour une antenne est passé par une chaîne radio fréquence (RF) 16, 15 identique pour chaque antenne, et émis. Le code correcteur, l'entrelaceur, le modulateur et la génération de séquences pilotes peuvent présenter des caractéristiques différentes choisies en fonction de l'utilisateur. De plus, le modulateur peut présenter des caractéristiques différentes choisies en fonction de l'antenne d'émission 20 pour un utilisateur donné. La figure 12 représente un schéma-bloc fonctionnel pour le récepteur dans le cas où le signal émis est obtenu avec une technique SCFDMA avec CP. Cependant, comme pour le cas SISO, si le bloc de de-allocation des sous-porteuses coïncide avec l'identité (c'est-à-dire s'il est 25 transparent), ce récepteur peut être utilisé dans un système SC ou SS-SC ou SS-SC-FDMA avec CP. En plus de cela, si à la place du bloc « suppression CP » on met un bloc d' « Extraction du bloc de traitement », le récepteur est adapté aux différents types de systèmes SC avec des séquences pilotes ou constantes à la place des CP. On suppose ici que le signal de l'émetteur est 30 synchronisé au récepteur avec une fiabilité adéquate (i.e. inférieure à la 44 durée du CP, si présent) et qu'un algorithme de synchronisation a fourni l'instant de synchronisation au récepteur. Comme pour le cas SISO (figure 3), le récepteur MU-MIMO de la figure 12 contient deux boucles de rétroaction : une boucle d'auto-itération 5 (itération demappeur souple / égaliseur) et une boucle de turbo-itération (itération demappeur souple / décodeur). En plus de cette structure, en présence de plusieurs utilisateurs, le récepteur doit décider, pour chaque boucle, les utilisateurs dont les signaux seront décodés / traités et quels signaux ne seront pas décodés/traités. Ici nous utilisons un procédé de io traitement en parallèle de tous les signaux de tous les utilisateurs dans la boucle d'auto-itération, et une procédé de décodage en série par utilisateur (typiquement à partir de l'utilisateur reçu avec le rapport signal sur interférence plus bruit le plus fort) pour la boucle de turbo-itération. D'autres ordonnancements des traitements des utilisateurs et de leurs antennes sont 15 cependant envisageables. En général, soit U = (1, ..., in l'ensemble des identifiants des utilisateurs et soit Pu = fui n = 0, j # j; Uctti = U} une partition de U. Par exemple, si U = (1,2,3} une partition pourrait être = (1,3), 44,2 = [2]. Un autre exemple d'ordonnancement est, pour chaque turbo-itération, de décoder d'abord les utilisateurs appartenant à puis faire 20 des auto-itérations et puis décoder les utilisateurs appartenant ài4,2 (sans ré- decoder ceux dans 44,1). La composition de la partition des utilisateurs et l'ordre de ses sous-ensembles définissent l'ordonnancement en ce qui concerne le décodage des utilisateurs dans la boucle turbo. Les regroupements d'utilisateurs dans la partition peuvent se faire selon plusieurs 25 métriques possibles, par exemple des utilisateurs ordonnancés selon la puissance de réception de leurs signaux, ou selon une métrique dérivée de la puissance de réception de leurs signaux et de la fiabilité/variance de leurs constellations, ou encore des métriques basées sur les écarts de puissance de réception et fiabilité entre couple d'utilisateurs, etc. 30 Comme pour le cas SISO, dans la suite, les paramètres en entrée ou en sortie des blocs 106, 108, 105T, à l'exception des observations initiales yR sur les différentes antennes, et des estimations initiales H, Ew, dépendent de l'identifiant de la turbo-itération courante -t- = O, ... a', de l'ensemble d'utilisateurs 44,i à décoder et éventuellement de l'identifiant de l'auto-itération courante s = 0, ..., ST. La lettre « T » est utilisée pour désigner 5 que l'on considère T antennes c'est-à-dire tous les utilisateurs confondus. Ce bloc a tous les T signaux estimés en entrée. Pour présenter la stratégie d'ordonnancement plus clairement, supposons que nous avons une partition Pu = -4,11) de "Li et une fonction d'ordre Ki) qui ordonne la partition de façon séquentielle 44,0,(i) du plus petit au plus grand selon un critère à 10 déterminer (des exemples ont été donnés précédemment). Dans la suite, nous supposons que cette fonction d'ordre a été appliquée à la partition pour chaque turbo-itération et nous dénommons la suite ordonnée de la partition (44, par simplicité et sans perte de généralité. A chaque turbo- itération, le récepteur essaye de décoder l'ensemble des utilisateurs -4,41_, puis 15 il fait des auto-itérations éventuelles, puis (au sein de la même turbo- itération) il essaie de décoder les utilisateurs dans 44,2 tout en sachant que des informations souples du décodage des utilisateurs dans 44,1 existent et ainsi de suite. A une turbo-itération donnée, le récepteur essaie donc de décoder 44,i en ayant les informations souples de décodage déjà pour tout 441 20 avec] 103 qui sélectionnent uniquement les éléments des vecteurs correspondants aux sous-porteuses où les signaux ont été envoyés (dé-allocation des sous-porteuses). En sortie de ces derniers blocs, les signaux sont représentés par R vecteurs yr de taille M. Pour faciliter et 5 simplifier la lecture du document par la suite, ces vecteurs peuvent être concaténés l'un à la suite de l'autre dans un vecteur y de taille MR. En parallèle, les séquences pilotes sont extraites des signaux en provenance de l'antenne de réception et envoyées dans le bloc « Estimation des canaux et de la variance du bruit » 104 et utilisées pour calculer une io estimation de la réponse fréquentielle des canaux MIMO entre les émetteurs et le récepteur, sur les M sous-porteuses d'intérêts (celles utilisées par les utilisateurs pour envoyer l'information). La réponse fréquentielle du canal entre l'antenne t de l'utilisateur u et l'antenne r du récepteur est organisée dans un vecteur h,-,t,u de taille M, ou de façon équivalente, dans une matrice 15 diagonale Hau = diag(hr,t,u) de taille M x M. Pour faciliter la lisibilité du document, les matrices correspondantes à une antenne de transmission t d'un utilisateur u peuvent être concaténées dans une matrice bloc-diagonale Htu de taille MR x M. Ces dernières peuvent être concaténées séquentiellement dans une matrice bloc-diagonale H de taille MR x MT, les 20 canaux des antennes des émetteurs sont concaténés sur la même ligne, pour la même antenne de réception. La quantité T =Euu=iTu est le nombre total de toutes les antennes d'émission, tous les émetteurs confondus. Le bloc « Estimation des canaux et de la variance du bruit » 104 fournit aussi une estimation de la variance du bruit. Ce bruit contient la contribution du 25 bruit thermique sur les antennes de réception et d'éventuelles interférences dues à la plusieurs effets en provenance soit de traitement numérique, soit d'autres signaux extérieurs. On suppose que le bloc « Estimation des canaux et de la variance du bruit » fournit la matrice de covariance du bruit dans le domaine fréquentiel Ew de taille MR x MR. Typiquement, on a une estimation 30 de canal et de covariance du bruit par bloc de données. 48 Une fois que les observations et les estimations du canal MUMIMO global H et de la covariance du bruit sont disponibles, la première boucle de traitement (auto-itération) est activée. Au tout premier passage = 0), ou à chaque début de turbo-itération (1- = 1, ...,T, si T > 0) et 5 également après chaque décodage d'un ensemble d'utilisateurs, le bloc 105T, démodulateur souple, est activé pour une phase d'initialisation. Par convention nous associons cette phase à l'indice d'auto-itération s = -1. Au tout début d'une turbo-itération, si aucun paquet d'aucun utilisateur n'a pas encore été décodé, une commande d'initialisation (s = -1) est générée avec io le bloc « Gestionnaire d'ordonnancement et d'initialisation » pour tous les utilisateurs, ce qui fait que le bloc 105T (l'explication est donnée en relation aux figures 12 et 13): a) vide ses registres internes ; b) met à zéro les entrées des U vecteurs (signaux égalisés) 54' de 15 longueur TuM formés par la concaténation des vecteurs 5-et de taille M (signal égalisé de l'antenne t de l'utilisateur u) ; c) pose les entrées des vecteursor?. - v,u = [0.1,2,1,u, - - - , 4,t,u, - - - 19-1/2 à l'infini ou toute autre valeur faisant office de l'infini ; d) à la toute première activation de la phase d'initialisation (1- = 0) les 20 informations a priori en provenance des blocs 112u sont mises par convention égales à des vecteurs de zéros ; e) les blocs 105T, de la figure 13 s'exécutent avec un indice d'auto-itération s = -1, notamment en utilisant d'éventuelles entrées non-nulles d'information a priori en provenance des décodeurs 111,, et plus 25 précisément des blocs 112u (pour tous les u = 1, U). La notation Tu est une convention indiquant à la fois que c'est le démodulateur de l'utilisateur u (en indice) et que l'entrée est égale aux signaux égalisés des Tu antennes de l'utilisateur u. Ainsi, pour tout u, le bloc 105T génère une estimation souple iu de longueur TuM formée par la 30 concaténation des vecteurs Yt,u de taille M (estimation souple du 49 signal de l'antenne t de l'utilisateur u) et le vecteur de leurs fiabilités correspondantes = [171,u, ---, --- 111/4,1,11, qui collecte une fiabilité par antenne de l'utilisateur u. Nous notons qu'à la toute première activation de la phase 5 d'initialisation (I- = 0) les estimations souples Yu seront des vecteurs nuls et 17" = pour tout t et u, où cx est l'énergie moyenne des symboles des constellations d'origine de l'utilisateur (on suppose ici que toutes les constellations ont la même énergie moyenne, ce qui n'est pas restrictif). la f) Les sorties vers les blocs 110u, désentrelaceurs, ne sont pas calculées parce qu'inutiles. Pendant une turbo-itération donnée, quand les paquets de l'ensemble des utilisateurs ui ont été décodés, en ayant ainsi généré des informations souples, plusieurs stratégies d'initialisation sont possibles. La 15 première est que le bloc « Gestionnaire d'ordonnancement et d'initialisation » figure 13 génère une commande d'initialisation (s = -1) pour tous les utilisateurs, ce qui produit une procédure d'initialisation comme celle décrite pour le cas « aucun paquet d'aucun utilisateur n'a pas encore été décodé » précédemment décrit. Une autre stratégie d'initialisation est que le bloc 20 « G estio n na i re d'ordonnancement et d'initialisation )> génère une commande d'initialisation (s = -1) pour tous les utilisateurs u E Uul, c'est-à-dire que le démodulateur souple MU-MIMO 105T utilise l'information a priori en provenance des décodeurs des utilisateurs qui viennent d'être décodés et des utilisateurs qui n'ont pas encore été décodés, alors que pour les autres 25 utilisateurs u E Ui 0, le bloc de calcul égaliseur 106 fournit en sortie à chaque auto- 30 itération les coefficients des égaliseurs ft,u de taille MR dans le domaine 51 fréquentiel, pour l'égalisation du signal en provenance de l'antenne t de l'utilisateur u. Ce filtre agit sur les observations concaténées y dans le domaine fréquentiel. Ceci est en général un turbo-égaliseur linéaire MMSE auto-itéré avec retour (feedback) EP. Les filtres d'un même utilisateur 5 peuvent être regroupés par commodité dans une matrice Fu = rf --- f Tu,u de taille MR x T. Le bloc de calcul des paramètres de l'égaliseur fournit aussi en sortie les vecteurs - v,u = [civ2,1,u, --- (31,t - - - , eV2 U = U qui collectent des estimations du bruit après égalisation pour chaque bloc de données de 10 l'antenne t de l'utilisateur u, calculées sous l'hypothèse de bruit gaussien et modèle additif î --t,u = xt,u + vt,u- Appelons o-,24, = diag(Zw) le vecteur des covariances sur la diagonale de la matrice de covariance du bruit. On calcule d'abord la matrice suivante (représentant la covariance du signal reçu en intégrant l'information souple en provenance du démodulateur souple MU- 15 MIMO souple) : U Tu E d = diag(crl) + (19) u=1 t=1 Cette matrice est diagonale par blocs. Le récepteur doit calculer son inverse et il existe des méthodes efficaces en termes de calculs ou de mémoire pour implémenter ce calcul. La matrice inverse est de taille MR x MR et elle est aussi diagonale par blocs qui ont une taille de M x M correspondante à la 20 longueur du bloc de données : diag(.1") - .- - diag(À1,R) = (20) diag(4,1) diag(* ÀR,R) où la matrice diag(.1,,) capture l'inverse de la cross-covariance du signal reçu entre deux antennes r et r' du récepteur. Les sorties sont calculées de la façon suivante, pour tout u= 1, U et t= 1, ..., Tu : (21) 52 fr,t,u = 1 r m = 0, M - 1 (22) lir ri Ilri,t,u, 1 1 _ ,.2 (23) t,u Le vecteur total ftu résulte de la concaténation des vecteurs des filtres fr tu _ pour r = 1, R. Les vecteurs de signaux reçus sur des antennes et transformés en fréquence yr de taille M sont passés au bloc « égaliseur linéaire à 5 suppression d'interférence » 108 qui reçoit aussi en entrée les coefficients des égaliseurs Fu, pour tous les utilisateurs u (qui sont traités parallèlement), la réponse fréquentielle globale du canal MU-MIMO H et les vecteurs de taille M d'estimations souples des symboles émis dans le domaine fréquentiel Zu, pour tous les utilisateurs u, qui sont générés à partir des io vecteurs de taille M d'estimations souples dans le domaine temporel Yu, en provenance du démodulateur souple MU-MIMO 105T. La fiabilité globale de ces estimations par antenne de l'utilisateur u est représentée par les valeurs de vii = 1,u,--- 115t,UP --- fiTwu] plus elle est proche de 0 plus la fiabilité est grande (ces valeurs représentent des variances). L'égaliseur LE-IC 108 15 réalise la suppression d'interférence (IC) pour l'ISI, le MAI et le MUI dans le cas général MU-MIMO. L'égaliseur 108 produit des vecteurs Xt,u de taille M qui représentent des estimations des symboles de chaque antenne de chaque utilisateur dans le domaine fréquentiel. Les vecteurs Tt,u sont ensuite passés à nouveau dans le domaine 20 temporel à travers Tu IFFT normalisées de taille M (109), pour obtenir des vecteurs représentant les signaux temporels égalisés, qui sont envoyés au démodulateur souple MU-MIMO 105T. Ce dernier, à partir des îct,u, des correspondantes variances du bruit après égalisation uv,t,u et des informations souples a priori, par exemple sous la forme de métriques Log- 25 Likelihood Ratios (LLRs) La(du) en provenance des décodeurs, produit les estimations souples des symboles émis dans le domaine temporel iu, pour 53 tous les utilisateurs u, ainsi que leurs fiabilités correspondantes eu. Si la boucle de turbo-itération est activée après, le demappeur souple produit aussi des informations souples pour chaque bit des signaux des utilisateurs appartenant à par exemple sous la forme de LLR extrinsèques Le(du), 5 U E A ce moment, si l'ordonnancement commande d'exécuter une boucle d'auto-itération, les traitements poursuivent avec une réactivation des blocs 107, 106, 108, 109, 105T. Autrement, si l'ordonnancement commande d'exécuter une io nouvelle boucle d'auto-itération ou de continuer la boucle de turbo-itération courante, les métriques souples Le(du), U E /Li sont ensuite désentrelacées par les blocs 110u pour les utilisateurs u Ensuite, les métriques souples sont envoyées aux décodeurs 111, des utilisateurs u E ui qui produisent des estimations des probabilités des bits d'informations envoyés (par exemple 15 des probabilités a posteriori), si la dernière turbo-itération est atteinte, i.e. T = T. Elles sont utilisées pour décider les bits envoyés avec des détecteurs à seuil. Autrement, les décodeurs 111u des utilisateurs u e u1 produisent aussi des estimations des bits codés envoyés c'est-à-dire des probabilités extrinsèques (EXT), par exemple sous forme de LLRs, qui mesurent la 20 probabilité que les bits codés envoyés soient 0 et 1. Les EXT sont ensuite envoyées aux entrelaceurs 112, des utilisateurs u e Les EXT entrelacées des utilisateurs u e ui deviennent des informations a priori La (du) et entrent dans le démodulateur souple MU-MIMO 105T. Avant de repasser par le démodulateur souple MU-MIMO, l'indice qui parcourt la partition des 25 utilisateurs Pu est incrémenté à i + 1, sauf si i = I. Dans ce dernier cas, on vient de décoder les derniers utilisateurs de la partition, donc l'indice pour parcourir la partition est réaffecté à i = 1 et l'indice de turbo-itération est incrémenté à T + 1 (nouvelle turbo-itération qui commence). Le démodulateur souple MU-MIMO 105T traite ses entrées et la 30 procédure continue en fonction de l'activation d'une auto-itération ou d'une turbo-itération. 54 Le bloc « égaliseur linéaire à suppression d'interférence » 108 réalise la suppression d'interférence (IC) pour l'ISI, le MAI et le MUI et égalise les signaux reçus dans le domaine fréquentiel et spatial. Comme pour le cas SISO, la différence par rapport à l'état de l'art existant réside 5 dans les entrées Xu (estimations des symboles dans le domaine fréquentiel pour tous les utilisateurs), ainsi que les filtres utilisés, ces quantités étant calculées selon le procédé de l'invention. La figure 14 décrit en détail le filtrage et le procédé de suppression de l'interférence. Le bloc de concaténation 204 sert simplement à concaténer io les vecteurs des estimations souples Xi, de taille M; dans un unique vecteur X de taille MT. Le bloc « Génération des signaux reçus estimés » 205 génère une estimation des R signaux reçus (partie utile plus interférences de tout type) au récepteur, concaténés dans le vecteur de sortie HY de taille MR. Nous notons que pour r = s = 0, X est un vecteur de zéros (pas d'estimation 15 disponible). Le séparateur 206 sépare le vecteur Hi en R vecteurs de taille M, un pour chaque voie du récepteur. Ces derniers sont soustraits par les sommes vectorielles 207 aux symboles reçus dans le domaine fréquentiel yr. On voit bien que le procédé est ici en parallèle, dans le sens que toutes les estimations disponibles dans la boucle d'auto-itération des signaux des 20 utilisateurs sont utilisées pour générer les signaux reçus. Ceci permet de prendre en compte (à travers l'estimation de canal global) l'influence que les utilisateurs ont les uns sur les autres, et que les antennes ont les unes sur les autres. Les signaux correctifs en sortie des sommes 207, sont envoyés à un banc de filtres MMSE linéaire 208 qui appliquent singulièrement le même 25 traitement du cas SISO mais en utilisant les filtres fr tu calculés par le bloc 106, pour tout u et t. Cela est fait pour toutes les antennes de réception. Les signaux correctifs égalisés (de taille M) en sortie des filtres (égaliseurs) 208 et qui correspondent au même utilisateur u et à la même antenne t, sont sommés ensemble par les sommes vectorielles 209. Cela revient à combiner 30 ensemble des estimations des mêmes signaux égalisés en provenance d'antennes de réception différentes pour améliorer la fiabilité des estimations. Les blocs de concaténation 210u concatènent ensemble les signaux des Tu antennes d'un même utilisateur u. Enfin, les sommes vectorielles 211u génèrent les estimations fréquentielles (de taille TuM) des 5 symboles égalisés en ajoutant les estimations des signaux utiles aux signaux correctifs filtrés en provenance des bancs de filtres 208, comme une généralisation des sommes issues du bloc 201 au cas SISO. Structure du démodulateur souple MU-MIMO Nous détaillons par la suite la structure du démodulateur souple 10 MU-MIMO 105T. Comme il a été écrit pour la phase d'initialisation, le démodulateur souple MU-MIMO 105T représenté en figure 13, se compose de U démodulateurs souples MIMO séparés 105Tu, qui traitent séparément les signaux égalisés correspondant aux antennes d'émission d'un même utilisateur. Le gestionnaire d'ordonnancement et d'initialisation 311 (figure 15 13) active tous les démodulateurs souples MIMO pendant les auto-itérations en générant ainsi toutes les estimations souples .Xu et leurs fiabilités moyennes correspondantes eu afin que les blocs implémentant les fonctions liées à l'égalisation à suppression d'interférence puissent avoir des estimations des signaux de tous les utilisateurs. En revanche, quand un 20 cycle d'auto-itérations prend fin parce que s = ST, le gestionnaire d'ordonnancement configure seulement les démodulateurs souples MIMO des utilisateurs u E u qui doivent être décodés selon l'ordonnancement courant pour que ces démodulateurs calculent les LLRs extrinsèques Le(du) à envoyer aux décodeurs des utilisateurs intéressés. Quand les nouvelles 25 informations a priori La(du), u E -ut sont disponibles à l'entrée du démodulateur souple MU-MIMO 105T, le gestionnaire 311 va déclencher l'initialisation des démodulateurs souples MIMO 105Tu selon la stratégie choisie (voir description plus haut pour le turbo-récepteur MU-MIMO de figure 12). Si i = I et r = T, la boucle turbo est terminée et la dernière turbo- 30 itération est achevée. Dans ce cas il faut réinitialiser r = 0, S = -1, j = 1 pour le prochain paquet. Si i = I mais r REVENDICATIONS1 - Procédé pour améliorer l'estimation des symboles d'un signal modulé et de la fiabilité de l'estimée, le signal comprenant des blocs de données 5 comportant au moins les étapes suivantes exécutées au niveau d'un processeur: - Une étape de calcul d'estimations souples (301) du signal par une technique de maximum de vraisemblance exacte ou approximée (304), en utilisant les symboles estimés X du bloc de données courant, les symboles io estimés du bloc étant caractérisés par une mesure de fiabilité moyenne (302), - Une étape de division gaussienne (306) du signal estimé par maximum de vraisemblance, supposé blanc, par le signal estimé X pour générer une estimation souple du signal transmis, l'estimation étant caractérisée par une 15 mesure de fiabilité par bloc de données, - Une étape de lissage adaptatif (307) sur le signal estimé par la division gaussienne et sur les mesures de fiabilités, pour fournir une estimation ï du signal transmis et une mesure de fiabilité e associée. 20 2 - Procédé selon la revendication 1 comprenant en outre: - Une étape de calcul de distributions a priori du signal transmis (300) à partir d'entrées/sorties souples binaires, - Une étape de calcul d'estimations souples (301) du signal transmis, par la technique du maximum a posteriori (305), en utilisant le signal 25 estimé par le maximum de vraisemblance (304) et les distributions a priori calculées à l'étape précédente, ladite estimation est caractérisée par une mesure de fiabilité moyenne (302). 3 - Procédé selon l'une des revendications 1 ou 2 caractérisé en ce qu'il 30 comporte en outre : 61 - Une étape de décodage par un décodeur binaire à entrées/sorties souples (507, 111), - Une étape d'initialisation du démodulateur souple (504,105), en calculant une estimation a priori du signal transmis (300,301), ladite estimation étant caractérisée par une fiabilité moyenne (302) en amont de l'étape de calcul d'estimations souples (301) du signal transmis. 4 - Procédé selon la revendication 3 comprenant en outre une étape de désentrelacement (506) et d'entrelacement (508) des informations souples 10 sur le signal transmis, lors de l'échange itératif entre un démodulateur souple (105) et un décodeur binaire à entrées/sorties souples (111). - Procédé selon l'une des revendications 1 à 4 pour la suppression d'interférence au sein d'un signal issu d'au moins un émetteur reçu sur un 15 récepteur comprenant une ou plusieurs antennes, ledit signal comprenant un ou plusieurs blocs de données caractérisé en ce qu'il comporte en outre les étapes suivantes : - Une étape de calcul des paramètres du récepteur (503), en utilisant la connaissance du canal, 20 - Une étape de traitement du signal reçu par un récepteur à suppression d'interférence (502) qui génère une autre estimation X du signal transmis en utilisant une estimation préalable .X du signal transmis x et la connaissance du canal, l'estimation du signal étant effectuée par le récepteur à suppression d'interférence, 25 - Le calcul de la fiabilité cfe de l'estimation du signal transmis X obtenue par le récepteur à suppression d'interférence (502), en utilisant les paramètres du récepteur, - Une étape d'échange itératif des estimations du signal transmis, entre le récepteur à suppression d'interférence (502) et le démodulateur 30 souple (504,105), 62 - Une étape de calcul d'information souple sur les bits du signal transmis, faite à partir du signal démodulé (309-310). 6 - Procédé de suppression d'interférence d'un signal reçu selon la s revendication 5 comprenant en outre: - Une étape de calcul des paramètres du récepteur adaptatif (503), en utilisant la connaissance du canal et la fiabilité des retours souples du démodulateur (504), - Une étape d'échange itératif des fiabilités des estimations du signal transmis entre le bloc qui calcule les paramètres du récepteur adaptatif (503) et le démodulateur souple (504,105). 7 - Procédé de suppression d'interférence d'un signal reçu selon la revendication 5 comportant: - Une étape d'estimation des paramètres du canal entre l'émetteur et le récepteur (501) en utilisant le signal reçu, - Une étape de calcul des paramètres du récepteur adaptatif ou non (503), en utilisant l'estimation des paramètres du canal et la fiabilité des retours souples du démodulateur (504), - Une étape d'échange itératif des fiabilités des estimations du signal transmis entre le bloc qui calcule les paramètres du récepteur adaptatif (503) et le démodulateur souple (504,105). 8 - Procédé de suppression d'interférence d'un signal reçu selon la revendication 7, comprenant en outre une étape itérative d'estimation des paramètres du canal (501), en utilisant l'estimation souple du signal transmis et sa fiabilité en provenance du démodulateur souple (504,105), conjointement avec le signal reçu. 9 - Procédé de suppression d'interférence d'un signal reçu selon l'une des revendications 5 à 8 dans lequel, l'étape de traitement du signal reçu par un 63 récepteur à suppression d'interférence (502) et l'étape de calcul de ses paramètres (503) comprennent: - Une étape de conversion du signal reçu dans le domaine fréquentiel (102), - Une étape de calcul des paramètres de l'égaliseur dans le domaine fréquentiel (106), qui réalise les fonctions du bloc (503). - Une étape d'égalisation linéaire à suppression d'interférence (108) dans le domaine fréquentiel, paramétrée selon les sorties obtenues lors du calcul des paramètres de l'égaliseur dans le domaine 10 fréquentiel (106), - Une étape de conversion du signal traité dans le domaine temporel (109), - Une étape de conversion dans le domaine fréquentiel (107) de l'estimation souple du signal transmis fournie par le démodulateur 15 souple. 10 - Procédé selon l'une des revendications 5 à 8 dans lequel, l'étape de traitement du signal reçu par un récepteur à suppression d'interférence (502) et l'étape de calcul de ses paramètres (503) comprennent: 20 - Une étape de séparation du signal reçu y en P segments y de Np symboles chacun, avec recouvrement, p = 1, ...,P. - Une étape de conversion de chacun des segments du signal reçu dans le domaine fréquentiel (102), - Une étape du calcul des paramètres des P égaliseurs dans le 25 domaine fréquentiel, qui réalise les fonctions du bloc (503), - Une étape d'égalisation linéaire à suppression d'interférence (108) dans le domaine fréquentiel de chacun des segments, paramétrée selon les sorties obtenues lors du calcul des paramètres des P égaliseurs dans le domaine fréquentiel (503), 64 - Une étape de conversion de chacun des signaux traités dans le domaine temporel (109), et d'extraction des symboles non recouverts pour obtenir une estimation 5e, - Une étape de conversion dans le domaine fréquentiel (107) de s l'estimation souple du signal transmis fournie par le démodulateur souple. 11 - Procédé selon la revendication 10 où le calcul des paramètres des P égaliseurs dans le domaine fréquentiel utilise la connaissance de la 10 statistique de l'interférence entre segments, estimée à partir de la fiabilité des retours du démodulateur souple (504). 12 - Procédé selon l'une des revendications 10 ou 11 où l'interférence résiduelle entre segments est régénérée en utilisant les retours du 15 démodulateur souple (504), et soustraite de chacun des P segments yp, avant l'étape de conversion dans le domaine fréquentiel. 13 - Procédé selon l'une des revendications 5 à 8, dans lequel, l'étape de traitement du signal reçu par un récepteur à suppression d'interférence (502) 20 et l'étape de calcul de ses paramètres (503) comprennent: - Un sur-échantillonnage du signal reçu avec un facteur de sur-échantillonnage o, - Une étape de conversion en fréquence (102) du signal reçu suréchantillonné, 25 - Une étape de calcul du filtre de l'égaliseur dans le domaine fréquentiel suréchantillonné, - Une étape d'égalisation linéaire à suppression d'interférence (108) dans le domaine fréquentiel suréchantillonné selon la sortie obtenue lors du calcul des paramètres de l'égaliseur dans le domaine 30 fréquentiel suréchantillonné (503), 65 - Une étape de conversion du signal traité dans le domaine temporel (109), suivi d'un échantillonneur de facteur d'échantillonnage o égal au facteur de sur-échantillonnage du signal reçu. - Une étape d'interpolation idéale de l'estimation souple du signal 5 transmis fournie par le démodulateur souple par un facteur de suréchantillonnage o égal à celui du signal reçu (cela consiste en l'ajout de o - 1 zéros entre deux échantillons), suivi d'une conversion dans le domaine fréquentiel suréchantillonné (107). 10 14 - Procédé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type SC-FDMA ou SS-SC-FDMA comprenant une étape de mise en trame utilisant respectivement une modulation SC-FDMA ou SSSC-FDMA et une étape d'exécution du procédé selon l'une des revendications 9 à 13. 15 15 - Procédé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type SC ou SS-SC comprenant une étape de mise en trame utilisant respectivement une modulation SC ou SS-SC et une étape d'exécution du procédé selon l'une des revendications 9 à 13. 20 16 - Procédé pour supprimer des interférences au sein d'un signal reçu sur un récepteur de type OFDM ou OFDMA comprenant une étape de mise en trame utilisant respectivement une modulation OFDM ou OFDMA et une étape d'exécution du procédé selon l'une des revendications 5 à 8 dans 25 lequel, l'étape de traitement du signal reçu par un récepteur à suppression d'interférence (502) et l'étape de calcul de ses paramètres (503) comprennent: - Une étape de conversion du signal reçu dans le domaine fréquentiel (102), 66 - Une étape de calcul des paramètres de l'égaliseur dans le domaine fréquentiel (106), qui réalise les fonctions du bloc (503), - Une étape d'égalisation linéaire à suppression d'interférence (108) dans le domaine fréquentiel, paramétrée selon les sorties obtenues s lors du calcul des paramètres de l'égaliseur dans le domaine fréquentiel (106). 17 - Procédé de suppression d'interférence d'un signal reçu selon l'une des revendications 9 à 16, comprenant en outre: 10 - Une étape de génération du signal à l'émetteur selon un code espace- temps par blocs de type Alamouti ou autre code espace-temps par bloc de type orthogonal, - Une étape d'estimation des paramètres des canaux entre les antennes de l'émetteur et celles du récepteur, 15 - Une étape du calcul des paramètres du récepteur dans le domaine fréquentiel, qui prennent en compte la structure du code espace-temps par blocs utilisé, - Une étape d'égalisation linéaire à suppression d'interférence dans le domaine fréquentiel en utilisant les filtres obtenus lors du calcul des 20 paramètres des égaliseurs dans le domaine fréquentiel, et qui combine les signaux reçus en entrée selon la méthode associé au code espace-temps par blocs utilisé par l'émetteur. 18 - Procédé pour supprimer des interférences au sein des signaux issus 25 d'un émetteur reçus sur un récepteur, dans un réseau de type MIMO où l'émetteur est doté de plusieurs antennes de transmissions et d'un entrelaceur qui répartit les bits codés sur les antennes disponibles, comprenant une étape d'exécution du procédé selon l'une des revendications 5 à 16 comprenant en outre : 67 - Une étape de traitement des signaux reçus en provenance des antennes de l'émetteur par un récepteur à suppression d'interférence (502), qui en utilisant ses paramètres et une estimation préalable des signaux transmis it par les antennes d'émission t = 1, ..., T1, génère 5 une autre estimation des signaux émis Xt (207-211, 205), - Une étape du calcul des paramètres du récepteur (503) qui prend en compte la présence de plusieurs antennes d'émission et les estimations des paramètres des canaux entre les antennes de l'émetteur et du récepteur, 10 - Une étape de démodulation souple (105Tu) où l'estimation du signal de chaque antenne d'émission est traitée indépendamment. 19 - Procédé pour supprimer des interférences au sein des signaux issus de plusieurs émetteurs reçus sur un récepteur, dans un réseau de type MU- 15 MIMO où les émetteurs sont dotés d'une ou plusieurs antennes de transmissions et d'un entrelaceur qui répartit les bits codés sur les antennes disponibles, comprenant une étape d'exécution du procédé selon l'une des revendications 5 à 16 comprenant: - Une étape de traitement des signaux reçus en provenance des 20 émetteurs par un récepteur à suppression d'interférence (502), qui en utilisant ses paramètres et une estimation préalable des signaux transmis ït,u par les antennes d'émission t = 1, ..., Tu des utilisateurs= 1, ..., U, génère une autre estimation des signaux émis x(204-21 1), - Une étape du calcul des paramètres du récepteur qui prend en 25 compte la présence de plusieurs utilisateurs, éventuellement avec plusieurs antennes d'émission, et d'éventuelles estimations des paramètres des canaux entre les antennes des émetteurs et du récepteur, 68 - Une étape de démodulation souple (105T) où l'estimation du signal de chaque antenne d'émission de chaque utilisateur est traitée indépendamment, - Une étape d'ordonnancement (311) des utilisateurs à décoder selon 5 un critère basé sur la puissance des signaux reçus, ou selon une métrique dérivée de la puissance des signaux reçus et de la fiabilité de leurs constellations, ou encore selon une métrique basée sur les écarts de puissances et fiabilité entre couples d'utilisateurs. 10 20 - Récepteur comportant au moins une antenne pour recevoir un signal émis par au moins un émetteur et un processeur configuré pour exécuter le procédé d'amélioration de calcul des symboles et leurs fiabilités selon l'une des revendications 1 à 17, lesdits symboles étant utilisés pour la modulation d'un signal. 15 21 - Récepteur selon la revendication 20 caractérisé en ce qu'il comporte un démodulateur souple comprenant au moins les éléments suivants : - Un module (300) de calcul des distributions a priori, - Un premier module (301) d'estimations souples, ayant une sortie 20 reliée à un module de calcul de distributions a posteriori (305), - Un module de calcul des vraisemblances gaussiennes (304) ayant une sortie reliée au module de calcul de distributions a posteriori (305), dont la sortie est reliée à un deuxième module d'estimations souples (301), 25 - Le deuxième module d'estimations souples (301) transmet ses résultats directement à un module de division gaussienne (306) et à un dispositif effectuant une moyenne dont la sortie est reliée au module de division gaussienne, 69 - Le module de division gaussienne est relié à un lissage adaptatif (307) fournissant les valeurs des symboles estimées ainsi que leurs variances. 5