FR 3057723 A1 20180420 FR 1659993 A 20161014 1 PRECODEURS CONSTRUITS A PARTIR DE PARAMETRES D'ESTIMATION DE CANAUX DE TRANSMISSION La présente invention concerne le domaine des télécommunications, notamment les systèmes de communication comportant une pluralité de dispositifs émetteurs et une pluralité de dispositifs récepteurs, tels que les systèmes entrées multiples sorties multiples, ou MIMO pour « Multiple- Input Multiple-Output » en anglais. Plus particulièrement, l'invention propose un procédé de traitement d'un signal à émettre par un dispositif émetteur vers un dispositif récepteur, incluant une construction efficace de matrices de précodage. De manière non limitative, l'invention a une application privilégiée dans le cas d'un réseau de télécommunications sans fil multicellulaire, en voie descendante, dans lequel les dispositifs émetteurs sont par exemple des stations de base, chaque station de base étant associée à une cellule du réseau et pouvant servir une pluralité d'utilisateurs munis de dispositifs récepteurs tels que des terminaux. Un tel réseau est notamment soumis à des interférences multi-utilisateurs, qui représentent un obstacle majeur dans les réseaux de télécommunications sans fil. Il est connu de modéliser les liaisons entre les stations de base et les terminaux, ainsi que la topologie d'interférences en résultant, par un canal à interférences, ou IBC pour « Interférence Broadcast Channel » en anglais. Un canal IBC modélise un système (ou de façon équivalente, un réseau) de communication dans lequel plusieurs dispositifs émetteurs essaient de communiquer avec plusieurs dispositifs récepteurs mais chaque dispositif récepteur n'est lié qu'à un unique dispositif émetteur. Des techniques d'élimination de telles interférences existent, telles que les méthodes d'alignement d'interférence, ou IA pour « Interference Alignment » en anglais, qui proposent de confiner les interférences vues par un dispositif récepteur particulier dans un sous-espace vectoriel de dimension réduite en appliquant au signal utile, au niveau du dispositif émetteur servant ce dispositif récepteur, une matrice de précodage. Ces techniques nécessitent généralement une connaissance parfaite de l'état de canal, ou CSI pour « Channel State Information » en anglais, au niveau de l'émetteur. En pratique, l'état de canal de l'émetteur n'est pas connu parfaitement, et on ne dispose que de certains paramètres d'estimation de cet état de canal, tels que l'estimé de la matrice de canal (que 3057723 2 l'on appelle « estimé de canal ») ou la matrice de covariance de l'erreur d'estimation. Des méthodes ont été développées dans le cas où l'état de canal n'est pas connu parfaitement, qui reposent sur l'utilisation de paramètres d'estimation. Toutefois, ces méthodes ne fournissent généralement pas des résultats complètement satisfaisants, notamment car elles n'exploitent pas 5 suffisamment les informations contenues dans la matrice de covariance de l'erreur d'estimation. Par ailleurs, certaines de ces méthodes ne peuvent pas être utilisées dans le cas où l'estimé de canal est nul et où on dispose seulement des matrices de covariance de l'erreur d'estimation. Il y a ainsi un besoin pour des méthodes utilisant de manière optimale les informations issues des différents paramètres d'estimation de canaux de transmission. 10 La présente invention vient améliorer la situation. A cet effet, la présente invention propose une méthode de traitement d'un signal pour un système de communications pour lequel la connaissance de l'état de canal est imparfaite, qui repose sur l'utilisation conjointe d'un estimé du canal de communication et de la matrice de covariance de l'erreur d'estimation du canal. 15 La présente invention vise alors un procédé de traitement de signal pour un système de communication comprenant des stations de base, chaque station de base étant apte à servir un ensemble respectif de terminaux d'utilisateurs, le procédé comportant : - la réception d'une pluralité de paramètres, relatifs chacun à une estimation du canal de transmission d'une station de base vers un terminal ; 20 - le calcul itératif d'une pluralité de matrices de covariance de transmission, chaque matrice de covariance de transmission étant relative à un terminal, le calcul comprenant au moins une itération des étapes suivantes : pour chaque station de base parmi les stations de base : a) pour chaque terminal parmi l'ensemble respectif de terminaux, calcul d'une première matrice 25 en fonction d'au moins un paramètre parmi la pluralité de paramètres et d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission ; b) pour chaque terminal parmi l'ensemble respectif de terminaux, mise à jour de la matrice de covariance de transmission associée audit terminal, à partir d'au moins un paramètre parmi la pluralité de paramètres, d'au moins une matrice de covariance de transmission parmi la pluralité de 30 matrices de covariance de transmission et de la première matrice calculée pour ledit terminal ; 3057723 3 - prise en compte des matrices de covariance de transmission mises à jour pour émettre des signaux par les stations de base vers les terminaux. Les paramètres d'estimation reçus peuvent être, pour chaque station de base bi et pour chaque terminal i : 5 - un estimé Hi,bi du canal de transmission de ladite station de base bi vers ledit terminal i ; et - une matrice de covariance Cp j,bi d'une erreur d'estimation du canal de transmission de ladite station de base bi vers ledit terminal i. Un tel procédé permet de calculer des matrices de précodage dans le cas d'un système de télécommunication dans lequel les canaux ne sont pas connus parfaitement, en exploitant au 10 maximum les informations contenues dans les paramètres d'estimation de ces canaux. En particulier, la méthode fait apparaître des équations dans lesquelles figurent un certain nombre de termes dépendant de Cp j,bi, ce qui fait que l'information contenue dans les matrices de covariance des erreurs d'estimation est avantageusement utilisée. En outre, l'étape a) peut comprendre, pour chaque terminal parmi les ensembles de terminaux, 15 le calcul d'au moins une matrice de covariance relative à des interférences entre ledit terminal et des terminaux, tous distincts dudit terminal, parmi les ensembles de terminaux. Le calcul de la première matrice pour un terminal parmi l'ensemble respectif de terminaux peut alors être fonction d'au moins une matrice de covariance relative à des interférences entre ledit terminal et des terminaux, tous distincts dudit terminal. 20 Ainsi, les différentes interférences possibles sont avantageusement prises en compte dans le calcul des matrices de précodage. Avantageusement, les matrices de covariance relatives à des interférences entre un terminal donné et une pluralité de terminaux sont des fonctions d'au moins un paramètre d'estimation et d'au moins une matrice de covariance de transmission. 25 Dans un mode de réalisation, deux matrices de covariance relatives à des interférences sont calculées pour chaque terminal k parmi les K terminaux. Ainsi, le procédé comporte : - le calcul d'une première matrice de covariance iirk relative à des interférences entre ledit terminal k et une pluralité de terminaux et - le calcul d'une deuxième matrice de covariance relative à des interférences entre ledit 3057723 4 terminal k et une pluralité de terminaux. Selon un mode de réalisation, pour chaque terminal k, la première matrice de covariance iirk relative à des interférences est fonction de Eri,j,k tr {QiCp j,bi} et la deuxième matrice de covariance Ri k relative à des interférences est fonction de Ek tr {QiCp j,bi}, où Qi est la matrice 5 de covariance de transmission associée à un terminal j et trtM1 est une trace d'une matrice M. Plus particulièrement, la première matrice de covariance iirk relative à des interférences peut être fonction de EiKi,j,k + tr {QiCp j,bi} Kif et la deuxième matrice de covariance relative à des interférences peut être fonction de EjK,k [Hi,bi QiHt;Ibi + tr {QiCp K3], où Qi est la matrice de covariance de transmission associée à un terminal j, et K1 et K3 sont des matrices qui 10 ne dépendent pas d'une station de base ou d'un terminal. En particulier, les matrices K1 et K3 peuvent être choisies comme étant égales à un scalaire (qui peut être différent pour K1 et K3) multiplié par une matrice identité. Selon ce mode de réalisation, pour une station de base bk donnée et pour un terminal k parmi l'ensemble respectif de terminaux associé à bk, la première matrice peut être exprimée comme 15 étant une fonction de kk et kk, avec : KG ITJH FA61v Ir "i,bk-`rl("1,bk ; ar,e i,k i,bk ik 1,bk f -t ; où K2 et K4 sont des matrices qui ne dépendent pas d'une station de base ou d'un terminal, MEI est une matrice adjointe d'une matrice M et M-1 est une matrice inverse d'une matrice M. Là encore, les matrices K2 et K4 peuvent être choisies comme étant égales à un scalaire (qui peut être différent pour K2 et K4) multiplié par une matrice identité. 20 Selon un mode de réalisation particulier, le procédé comprend en outre la réception d'une pluralité de coefficients de pondération, chaque coefficient de pondération ui étant associé à un terminal i parmi les ensembles de terminaux. La première matrice peut alors s'écrire : K f Àk = [K,k (IMbk Qkfil,k) 1 OMbk Qke,k) 1] i#k Les coefficients de pondération peuvent tous être égaux, si aucun terminal ne doit être privilégié et 3057723 5 par rapport à un autre. Alternativement, il est possible d'attribuer des poids plus grands à des terminaux auxquels on souhaite attribuer un débit plus important. En variante au mode de réalisation précédent, le procédé peut en outre comprendre, pour chaque terminal i parmi les ensembles de terminaux, une réception d'une matrice de covariance Cr 5 relative à une réception de signal par ledit terminal i. Dans ce cas, la première matrice de covariance iirk relative à des interférences peut s'exprimer comme une fonction de EiK,i,j,k [Hi,bi QiHt;Ibi + tr {QiCp,i,bi Crl, et la deuxième matrice de covariance Ri k relative à des interférences peut s'exprimer comme une fonction de : bk [Hi,b,Q;Kbi + tr Le procédé peut en outre comprendre la réception d'une valeur G2 relative à une variance 10 associée à un bruit de signal. Dans ce cas la première matrice de covariance iirk relative à des interférences peut s'écrire iirk = + tr {QiCp,i,bi cr2 I et la deuxième matrice de covariance rti,k relative à des interférences peut s'écrire Ri,k = Ebk tr{Q~Cp+G21, où I est une matrice identité. Pour une station de base bk et pour chaque terminal k parmi l'ensemble respectif de terminaux, 15 la première matrice peut alors s'écrire comme une fonction de kk et kk, avec : k,k = H bkRrkHi,bk trtrtdiCr ,ijCp,i,bk et kk = Htibkrt2,Hi,bk + Selon un mode de réalisation particulier, le procédé comprend en outre la réception d'une pluralité de coefficients de pondération, chaque coefficient de pondération ui étant associé à un terminal i parmi les ensembles de terminaux. La première matrice peut alors s'écrire : K f Àk = [K,k (IMbk Qkfil,k) 1 OMbk Qke,k) 1] i#k Les coefficients de pondération peuvent tous être égaux, si aucun terminal ne doit être privilégié 20 par rapport à un autre. Alternativement, il est possible d'attribuer des poids plus grands à des terminaux auxquels on souhaite attribuer un débit plus important. Un dispositif destiné au traitement d'un signal pour un système de communication comprenant des stations de base, chaque station de base étant apte à servir un ensemble respectif de terminaux d'utilisateurs peut être avantageux, en lui-même, dès lors qu'il permet de calculer des matrices de 25 covariance de transmission qui peuvent être utilisées pour traiter et émettre des signaux par les 3057723 6 stations de base vers les terminaux. Ainsi, la présente invention vise également un dispositif comportant : - une interface d'entrée pour recevoir une pluralité de paramètres relatifs chacun à une estimation d'un canal de transmission d'une station de base parmi les stations de base vers un 5 terminal parmi les ensembles de terminaux ; - un circuit pour un calcul itératif d'une pluralité de matrices de covariance de transmission, chaque matrice de covariance de transmission étant relative à un terminal parmi les ensembles de terminaux, le calcul comprenant au moins une itération comportant : 10 pour chaque station de base parmi les stations de base : a) pour chaque terminal parmi l'ensemble respectif de terminaux, calcul d'une première matrice en fonction d'au moins un paramètre parmi la pluralité de paramètres et d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission ; 15 b) pour chaque terminal parmi l'ensemble respectif de terminaux, mise à jour de la matrice de covariance de transmission associée audit terminal, à partir d'au moins un paramètre parmi la pluralité de paramètres, d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission et de la première matrice calculée pour ledit terminal ; 20 - une interface de sortie pour fourniture des matrices de covariance de transmission mises à jour aux stations de base. Un programme informatique, mettant en oeuvre tout ou partie du procédé décrit ci-avant, installé sur un équipement préexistant, est en lui-même avantageux. Ainsi, la présente invention vise également un programme informatique comportant des 25 instructions pour la mise en oeuvre du procédé précédemment décrit, lorsque ce programme est exécuté par un processeur. Ce programme peut utiliser n'importe quel langage de programmation (par exemple, un langage objet ou autre), et être sous la forme d'un code source interprétable, d'un code partiellement compilé ou d'un code totalement compilé. 3057723 7 La figure 2 décrite en détails ci-après, peut former l'organigramme de l'algorithme général d'un tel programme informatique. D'autres caractéristiques et avantages de l'invention apparaîtront encore à la lecture de la 5 description qui va suivre. Celle-ci est purement illustrative et doit être lue en regard des dessins annexés sur lesquels : - la figure 1 illustre un système de télécommunications dans lequel peut être utilisée l'invention, selon un mode de réalisation ; - la figure 2 illustre un ordinogramme d'une réalisation possible dans un mode de 10 réalisation de l'invention ; - la figure 3 illustre une mise à jour des matrices de covariance de transmission, selon un mode de réalisation de l'invention ; - la figure 4 illustre un dispositif de traitement de signal selon un mode de réalisation de l'invention ; 15 - les figures 5a et 5b illustrent des résultats de simulations, pour différents algorithmes existants et pour la présente invention réalisée selon un mode de réalisation particulier. La figure 4 représente un exemple de dispositif de traitement d'un signal à émettre depuis dispositif émetteur vers un dispositif récepteur, dans un mode de réalisation de l'invention. Dans ce mode de réalisation, le dispositif comporte un ordinateur 400, comprenant une mémoire 20 405 pour stocker des instructions permettant la mise en oeuvre du procédé, les données de mesures reçues, et des données temporaires pour réaliser les différentes étapes du procédé tel que décrit précédemment. L'ordinateur comporte en outre un circuit 404. Ce circuit peut être, par exemple : - un processeur apte à interpréter des instructions sous la forme de programme 25 informatique, ou - une carte électronique dont les étapes du procédé de l'invention sont décrites dans le silicium, ou encore - une puce électronique programmable comme une puce FPGA (pour « Field-Programmable Gate Array » en anglais). 3057723 8 Cet ordinateur comporte une interface d'entrée 403 pour la réception de données 401 relatives à une estimation de l'état de canal, et une interface de sortie 406 pour la fourniture de matrices de transmission à un dispositif émetteur 402. Ce dispositif émetteur peut être par exemple une station de base. Dans un mode de réalisation de l'invention, l'ordinateur 400 est inclus dans le dispositif 5 émetteur 402. Alternativement, l'ordinateur 400 est externe au dispositif émetteur 402. Par ailleurs, le schéma fonctionnel présenté sur la figure 2 est un exemple typique d'un programme dont certaines instructions peuvent être réalisées auprès du dispositif décrit. A ce titre, la figure 2 peut correspondre à l'organigramme de l'algorithme général d'un programme 10 informatique au sens de l'invention. La figure 1 illustre un système de télécommunications dans lequel peut être utilisée l'invention, selon un mode de réalisation. Le système IBC représenté à la figure 1 comporte plusieurs cellules 101, 10C, chaque cellule 15 étant composée d'une station de base 11, 1C et d'un ensemble d'équipements d'utilisateurs 111, 111\11, 11Nc, 11K, tels que des terminaux. Par la suite, on peut aussi désigner ces équipements d'utilisateurs par « terminaux utilisateurs » ou « utilisateurs ». Une station de base d'une cellule peut envoyer des signaux aux utilisateurs de la même cellule. Pendant qu'elle envoie ces signaux, il peut y avoir des « fuites », appelées interférences inter- 20 cellules, vers les utilisateurs des autres cellules. Par exemple, en référence à la figure 1, la station de base 11 de la cellule 101 peut envoyer un signal aux utilisateurs 111, 111\11 appartenant à la cellule 101 (flèches en traits pleins partant de 11), et les interférences inter-cellules avec les utilisateurs 11Nc, 11K d'une autre cellule 10C sont représentées par les flèches en traits discontinus partant de 1C. 25 On suppose que le système comporte C stations de base, servant un nombre total d'utilisateurs égal à K. On note bk la station de base servant un utilisateur donné k, où k E {1, ..., K}. Le signal, vecteur de dimension dk x 1 reçu par l'utilisateur k de la cellule bk s'écrit : Yk = iik,bkGkXkHk,bkGiXi GiXi Vk (1) i#k,bi=bk j#bk i:bi=i où, pour tout i E fi, ..., KI et pour tout j E fi, ..., KI : 3057723 9 Xi représente le signal de dimensions di x 1 transmis pour l'utilisateur i, supposé être un signal blanc, de variance égale à l'unité ; di est le nombre de flux ; Hi,i, matrice de dimensions Ni x Mi (où Mi est le nombre d'antennes de la station de base j), 5 représente le canal de communication de la station de base j à l'utilisateur i ; Gi est une matrice de dimension Mbi x di représentant le filtre spatial de précodage (appelée aussi « matrice de précodage » ou « précodeur »). Enfin, vk est un bruit de signal, supposé suivre une loi normale complexe centrée et de variance G2 : vk-C,W(0, 152 INk). 10 Dans l'équation (1), le terme Hkbk Gkxk représente la composante du signal reçu émis par la base bk correspondant à la cellule à laquelle appartient l'utilisateur k. Le terme Ej=k, bi =bk Hk,bk GiXi représente les interférences intra-cellules, et Ei#bk Ei:bi=j Hk,j GiXi représente les interférences inter-cellules. 15 Le but de l'invention est de construire des précodeurs qui maximisent la somme pondérée des débits, ou WSR pour « weighted sum rate » en anglais. La WSR s'écrit : WSR = WSR(Q) = ukln det(INk + Rkillk,bkQkKbk) (2) k=1 où : pour i E {1, , K}, Qi est la matrice de covariance de transmission associée à l'utilisateur i, reliée au précodeur Gi par l'équation Qi = GiGti 20 Q = QK} représente l'ensemble des matrices de covariances de transmissions Qi, pour i E fi, , KI ; les uk sont les poids des débits ; et Rk, Rk sont respectivement la matrice de covariance totale et la matrice de covariance d'interférence-plus-bruit côté récepteur, définies par 3057723 10 Rk = et Rk = Hkj,iQiF111--(,bi + G2INk i#k où MH représente la matrice adjointe (ou matrice transconjuguée) d'une matrice M. La fonction WSR doit être maximisée avec les contraintes de puissance : tr[Qkl 0 : 1 (ukzi(c1)(1,1) Pk0,1) = 1 E (1) (1'1) k (E(2)0 ' 0 + Àbk k où le scalaire z+ représente le maximum entre un scalaire z et 0. 5 A l'issue de l'étape 304, on réalise à nouveau le deuxième test d'arrêt. Si le deuxième critère d'arrêt n'est pas vérifié, les étapes 302 à 304 sont réitérées. Si le deuxième critère d'arrêt est vérifié, les matrices Qk sont mises à jour pour tout k E {1, ..., K} lors d'une étape 306, dans laquelle : - On calcule P = Pk ; 10 - Si tr(P) > PUBS, où j est la cellule courante, on pose A = A ; sinon on pose A = A ; - Enfin, pour chaque utilisateur k appartenant à la cellule courante j, on pose Qk = Gkr Pk(Gkr )11. L'algorithme itératif complet de calcul des matrices de covariance de transmission Qk selon un mode de réalisation particulier de l'invention est donné ci-dessous : Pour k = 1, K, initialiser Qk Répéter jusqu'à un premier critère d'arrêt : Pour j = 1, C Poser A = 0 et A = Pour k tel que bk = j Calculer ;(k selon les équations (4) à (8) Incrémenter k Répéter jusqu'à un deuxième critère d'arrêt : = + Pour k tel que bk = j 3057723 15 Calculer Bk selon les équations (9) et (10) Calculer Gk =vecteurs propres(ijk, Ak ÀImbk) Normaliser Gk pour obtenir Gk Calculer Eci-) = (G)Hi3kG, et Ec2) = (G)1-1-ÀkG, Calculer Pk selon l'équation (11) Incrémenter k Calculer P = Pk Si tr(P) > PUBS, poser À = À ; sinon poser À = À Pour k tel que bk = j, poser Qk = G;(Pk(GOH Incrémenter j Les figures 5a et 5b illustrent des résultats de simulations, pour différents algorithmes existants et pour la présente invention réalisée selon un mode de réalisation particulier. Plus précisément, les figures 5a et 5b représentent la somme pondérée des débits obtenue avec 5 trois méthodes de calcul des précodeurs en fonction du rapport signal-sur-bruit, ou SNR pour « signal-to-noise ratio » en anglais, pour un système entrées multiples / sortie unique, ou MISO pour « Multiple-Input Single Output » en anglais, composé de 2 cellules, avec un total de 8 utilisateurs dont chacun dispose d'une seule antenne seulement. Les matrices de covariance de l'erreur d'estimation Cp,i,bk (pour tout i E {1, ..., K} et bk E 10 [1, ..., C}) sont considérées comme des matrices identité pour la figure 5a et comme des matrices à rangs réduits (rang = 2) pour la figure 5b. Dans le cas d'un système MISO, les matrices de covariance côté récepteur disparaissent dans les différentes équations. Les courbes 501a et 501b correspondent aux résultats de la méthode proposée par un mode de réalisation de l'invention. Les courbes 503a et 503b correspondent aux courbes obtenues à partir 15 d'une méthode conçue pour le cas d'une connaissance parfaite du canal, dans laquelle les matrices Hi,bk ont été remplacées par les estimés de canal Hi,bk . Enfin, les courbes 502a et 502b ont été obtenues à partir de la méthode présentée dans l'article « Sum Rate Maximization in the Noisy MIMO Interfering Broadcast Channel with Partial CSIT via the Expected Weighted MSE» de F. Negro et al., publié dans IEEE International Symposium on Wireless Communication 20 Systems (ISWCS), Paris, France, Août 2012. Cette méthode maximise le WSR par une minimisation de l'erreur quadratique du signal qui est en générale linéaire en Hi,i et dont 3057723 16 l'espérance converge ainsi vers une quantité linéaire en Hi,i, et prend pas suffisamment en compte les matrices Cp,i,bk Comme le montrent les courbes des figures 5a et 5b, la présente méthode fournit de meilleurs résultats que les deux autres algorithmes, notamment dans le cas où les matrices de covariance sont 5 de rangs réduits. Cela est en partie dû à une meilleure exploitation des informations contenues dans les matrices de covariance de l'erreur d'estimation Cp,i,bk Bien entendu, la présente invention ne se limite pas aux formes de réalisation décrites ci-avant à titre d'exemples ; elle s'étend à d'autres variantes. Par exemple, la présente méthode peut être utilisée dans n'importe quel système multi-cellules 10 ou multi-utilisateurs. Elle pourrait également être mise en place dans un système où une station de base peut servir plusieurs utilisateurs en même temps. REVENDICATIONS1. Procédé de traitement de signal pour un système de communication comprenant des stations de base (11, 1C), chaque station de base étant apte à servir un ensemble respectif de terminaux d'utilisateurs (111, 11N1, 11Nc, 11K), le procédé comportant : - réception d'une pluralité de paramètres (401) relatifs chacun à une estimation d'un canal de transmission d'une station de base (11, 1C) parmi les stations de base vers un terminal (111, 11N1, 11Nc, 11K) parmi les ensembles de terminaux ; - calcul itératif d'une pluralité de matrices de covariance de transmission (205), chaque matrice de covariance de transmission étant relative à un terminal (111, 11N1, 11Nc, 11K) parmi les ensembles de terminaux, le calcul comprenant au moins une itération comportant : pour chaque station de base (11, 1C) parmi les stations de base : a) pour chaque terminal parmi l'ensemble respectif de terminaux, calcul (203) d'une première matrice en fonction d'au moins un paramètre parmi la pluralité de paramètres (401) et d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission ; b) pour chaque terminal parmi l'ensemble respectif de terminaux, mise à jour (204) de la matrice de covariance de transmission associée audit terminal, à partir d'au moins un paramètre parmi la pluralité de paramètres (401), d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission et de la première matrice calculée pour ledit terminal ; - prise en compte des matrices de covariance de transmission mises à jour pour émettre des signaux par les stations de base (11, 1C) vers les terminaux (111, 11N1, 11Nc, 11K). 2. Procédé selon la revendication 1, dans lequel : l'étape a) comprend en outre, pour chaque terminal (111, 11N1, 11Nc, 11K) parmi les ensembles de terminaux, calcul d'au moins une matrice de covariance relative à des interférences entre ledit terminal et des terminaux, tous distincts dudit terminal, parmi les ensembles de terminaux ; et 3057723 18 le calcul (203) de la première matrice pour un terminal parmi l'ensemble respectif de terminaux est en outre fonction d'au moins une matrice de covariance relative à des interférences entre ledit terminal et des terminaux, tous distincts dudit terminal. 3. Procédé selon la revendication 2, dans lequel les au moins une matrice de covariance relative 5 à des interférences entre ledit terminal et une pluralité de terminaux sont des fonctions d'au moins un paramètre parmi la pluralité de paramètres (401) et d'au moins une matrice de covariance de transmission parmi les matrices de covariance de transmission. 4. Procédé selon l'une des revendications précédentes, dans lequel la pluralité de paramètres reçus (401) comprend, pour chaque station de base bi parmi les stations de bases (11, 1C) et pour 10 chaque terminal i parmi les ensembles de terminaux (111, 111\11, 11Nc, 11K) : - un estimé Hi,bi du canal de transmission de ladite station de base bi vers ledit terminal i ; et - une matrice de covariance Cp,i,bi d'une erreur d'estimation du canal de transmission de ladite station de base bi vers ledit terminal i. 5. Procédé selon la revendication 4, dans lequel tous les ensembles de terminaux représentent un 15 nombre total de K terminaux, et comprenant en outre, pour chaque terminal k parmi les ensembles de terminaux (111, 111\11, 11Nc, 11K) : calcul d'une première matrice de covariance iirk relative à des interférences entre ledit terminal k et une pluralité de terminaux et calcul d'une deuxième matrice de covariance Ri k relative à des interférences entre ledit 20 terminal k et une pluralité de terminaux. 6. Procédé selon la revendication 5, dans lequel : ladite première matrice de covariance iirk relative à des interférences est fonction de tr {Qi j#i,j#k et ladite deuxième matrice de covariance Ri k relative à des interférences est fonction de tr j#k 3057723 19 où Qi est la matrice de covariance de transmission associée à un terminal j, et trtM1 est une trace d'une matrice M. 7. Procédé selon la revendication 6, dans lequel : ladite première matrice de covariance iirk relative à des interférences est fonction de : i,bjQjHbj tr Ki] j#i,j#k 5 ladite deuxième matrice de covariance rtik relative à des interférences est fonction de : [Fli,biiibj tr {Qi K3] ; j#k où Qi est la matrice de covariance de transmission associée à un terminal j, et K1 et K3 sont des matrices qui ne dépendent pas d'une station de base ou d'un terminal. 8. Procédé selon l'une des revendications 6 et 7, dans lequel, pour une station de base bk et pour chaque terminal k parmi l'ensemble respectif de terminaux, la première matrice est fonction de 10 Ak 1,c et Ak, P avec : i, KG ITJH _L frfejtv "i,bk-`fk--1,bk ; r, -te i,k i,bk i,k 1,bk ^-.1 ; où : K2 et K4 sont des matrices qui ne dépendent pas d'une station de base ou d'un terminal ; le est une matrice adjointe d'une matrice M, et M-1 est une matrice inverse d'une matrice M. 9. Procédé selon la revendication 8, comprenant en outre : 15 - réception d'une pluralité de coefficients de pondération, chaque coefficient de pondération ui étant associé à un terminal i parmi les ensembles de terminaux (111, 111\11, 111\1c, 11K) ; dans lequel ladite première matrice pour ladite station de base bk et ledit terminal k est : et 3057723 20 10. Procédé selon la revendication 5, comprenant, pour chaque terminal i parmi les ensembles de terminaux (111, 111\11, 11Nc, 11K), une réception d'une matrice de covariance Cr,i relative à une réception de signal par ledit terminal i, dans lequel : ladite première matrice de covariance iirk relative à des interférences est fonction de : K + tr{Qj; ji,jk 5 ladite deuxième matrice de covariance rtik relative à des interférences est fonction de : tr j#k où trtM1 est une trace d'une matrice M. 11. Procédé selon la revendication 10, comprenant en outre une réception d'une valeur G2 relative à une variance associée à un bruit de signal, dans lequel : ladite première matrice de covariance iirk relative à des interférences est : K [Hi bj Qi H bj + tr {Qi + G21 ; j#i,j#k 10 ladite deuxième matrice de covariance relative à des interférences est : tr{Qj+ +G21 j#k où I est une matrice identité. 12. Procédé selon l'une des revendications 10 et 11, dans lequel, pour une station de base bk et pour chaque terminal k parmi l'ensemble respectif de terminaux (111, 111\11, 11Nc, 11K), la première matrice est fonction de kk et kk, avec : - ii bkRi kHi,bk + trtridiCr,i1Cp j,bk ; et ÀD _RH k +trtiidiCrjCp,i,bk i,k i,bk i, 1,bk 15 où le est une matrice adjointe d'une matrice M et M-1 est une matrice inverse d'une matrice M. 3057723 21 13. Procédé selon la revendication 12, comprenant en outre : - réception d'une pluralité de coefficients de pondération, chaque coefficient de pondération ui étant associé à un terminal i parmi les ensembles de terminaux (111, 111\11, 11Nc, 11K) ; dans lequel ladite première matrice pour ladite station de base bk et ledit terminal k est : K f Àk = [K,k (IMbk Qkfil,k) 1 OMbk Qke,k) 1] i#k 5 14. Dispositif de traitement de signal pour un système de communication comprenant des stations de base (11, 1C), chaque station de base étant apte à servir un ensemble respectif de terminaux d'utilisateurs (111, 111\11, 11Nc, 11K), le dispositif comportant : - une interface d'entrée (403) pour recevoir une pluralité de paramètres (401) relatifs chacun à 10 une estimation d'un canal de transmission d'une station de base (11, 1C) parmi les stations de base vers un terminal 111, 111\11, 11Nc, 11K) parmi les ensembles de terminaux ; - un circuit pour un calcul itératif d'une pluralité de matrices de covariance de transmission (205), chaque matrice de covariance de transmission étant relative à un terminal parmi les ensembles de terminaux, 15 le calcul comprenant au moins une itération comportant : pour chaque station de base parmi les stations de base : a) pour chaque terminal parmi l'ensemble respectif de terminaux, calcul (203) d'une première matrice en fonction d'au moins un paramètre parmi la pluralité de paramètres et d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de 20 covariance de transmission ; b) pour chaque terminal parmi l'ensemble respectif de terminaux, mise à jour (204) de la matrice de covariance de transmission associée audit terminal, à partir d'au moins un paramètre parmi la pluralité de paramètres, d'au moins une matrice de covariance de transmission parmi la pluralité de matrices de covariance de transmission et de la 25 première matrice calculée pour ledit terminal ; - une interface de sortie (406) pour fourniture des matrices de covariance de transmission mises 3057723 22 à jour aux stations de base. 15. Produit programme informatique comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 13, lorsque ce programme est exécuté par un processeur.