FR 3014269 A1 20150605 FR 1361885 A 20131129 Procédé de coordination d'émetteurs radio basé sur un codage de niveau de puissance émise et émetteur correspondant Domaine de l'invention La présente invention se rapporte au domaine des télécommunications. Au sein de ce domaine, l'invention se rapporte aux techniques de coordination entre émetteurs d'une même technologie ou de technologies différentes qui utilisent des ressources radio communes sur une même zone géographique, c'est-à-dire sur une zone où les communications de l'un perturbent les communications d'un autre, on dit qu'il y a un phénomène d'interférence. L'invention se place dans le contexte d'un système de télécommunication à K émetteurs et K récepteurs qui communiquent respectivement avec les K émetteurs. Un émetteur peut tout aussi bien être une station de base qu'un point d'accès d'un réseau d'accès. Le phénomène d'interférence peut être limité par des mécanismes de coordination pilotés par une entité de communication centrale à condition qu'il y ait une infrastructure (backhaul selon la terminologie anglo-saxonne) pour véhiculer les signaux de contrôle et de signalisation entre les émetteurs, les récepteurs et l'entité centrale. Ces mécanismes sont toutefois consommateur en bande passante. En outre, une telle infrastructure est absente lorsqu'un émetteur exploite une bande de fréquence sans licence dite bande libre (white space) et donc il ne peut y avoir de coordination possible. En ce qui concerne en particulier les transmissions basées sur la technologie sans fil de type WiFi qui s'appuie sur des standards de la série IEEE 802.11, leur qualité est menacée par l'explosion des objets communicants liés à l'Internet des objets (IoT) et l'explosion des réseaux de distribution avec remontée automatique d'information (smart grids). Les réseaux cellulaires ne sont pas exempts de risque car la tendance des déploiements à venir consiste à déployer des cellules dites femto (small cells) entre des stations de base pour augmenter la capacité tout en limitant l'infrastructure (backhaul) ce qui est contraire au besoin des mécanismes de coordination pour lutter contre les interférences. Art antérieur Dans le contexte de l'invention, la transmission des communications entre émetteurs et récepteurs est considérée faite par un canal dit équivalent. Par canal équivalent, il faut comprendre le canal de propagation lui-même ainsi que les opérations de filtrages linéaires tant à l'émission qu'à la réception qui mènent à la structure mathématique du rapport signal à bruit plus interférence retenue. Dans la suite, le terme générique canal est utilisé en lieu et place du terme canal équivalent. On désigne par gain de canal (équivalent), le carré du module du coefficient dominant de la réponse impulsionnelle du canal, le carré du module de la réponse en fréquence du canal pour la fréquence m ou bien encore le rapport entre la puissance du signal reçu par le récepteur sur toute la bande m et la puissance émise, rapport qui serait mesuré en l'absence d'interférence des autres émetteurs. L'indice m dit de canal ou de fréquence peut prendre une valeur parmi M valeurs. Par M on désigne le nombre de voies de retour disponibles pour le système de communication considéré, chaque voie de retour étant caractérisée par un rapport signal à bruit plus interférence SINR qui lui est propre. De [1], il est connu un algorithme d'allocation de puissance émise par un émetteur lors d'une communication avec un récepteur via une liaison point à point. Cet algorithme est connu sous l'acronyme IWFA (Iterative WaterFilling Algorithm). L'allocation de puissance entre les M canaux orthogonaux de la liaison point à point est déterminée en exploitant uniquement la connaissance du rapport signal à interférence plus bruit ou SINR (Signal to Interference plus Noise Ratio) connu pour chacun de ces canaux fréquentiels, Les gains/coefficients de canaux sont supposés soit constants soit constant par morceau. Dans le dernier cas, les gains de canaux sont donc stationnaires sur une durée appelée temps Tc de cohérence. L'IWFA opére sur une durée T inférieure ou égale à ce temps de cohérence, T selon un mode privilégié. L'algorithme considère un système de télécommunication à K émetteurs et K récepteurs, 2. L'algorithme IWFA découpe le temps en trames successives et chaque trame comprend K sous-trames. La durée d'une trame correspond à la durée d'une itération de l'algorithme IWFA. Un émetteur ne peut mettre à jour la puissance allouée qu'une seule fois par trame. La trame correspond donc à la durée minimale nécessaire pour que les K émetteurs mettent à jour leur puissance allouée. L'algorithme considère un ordre arbitraire entre les émetteurs pour la mise à jour de leur puissance. Sans perte de généralité, l'émetteur 1 est celui qui met à jour en premier et ainsi de suite. L'algorithme IWFA suppose implicitement un mécanisme de synchronisation au niveau trame (boucle externe) et au niveau intervalle de temps correspondant à une sous-trame (boucle interne). La figure 1 illustre un système de télécommunication avec deux émetteurs et deux récepteurs, K=2. Le nombre de canaux par liaison est M=2. Le coefficient de canal appelé aussi gain de canal entre l'émetteur i et le récepteur j pour le canal m est noté il est supérieur ou égal à zéro. La puissance allouée par l'émetteur i pour le canal m est notée Pi . Le SINR évalué par le récepteur j pour le canalm m est noté y.J . Le côté émission est noté Em et le côté réception est noté Re. L'expression du SINR est celle correspondant au signal reçu par un récepteur i du type sans annulation d'interférence servi par l'émetteur i et interféré par les émetteurs j avec j # i : (1) j#i L'algorithme comprend une phase d'initialisation qui correspond à l'indice de sous trame t égal à zéro. Pendant cette phase, chaque émetteur i E {1,2,...,10 alloue la puissance P (0) au canal me {1,2,...,M}. Pour chaque émetteur, les M puissances allouées forment un vecteur d'allocation. Les vecteurs d'allocation peuvent être choisis arbitrairement à condition que chacun respecte la contrainte de puissance totale Pimax disponible à l'émetteur i considéré : M (0) max (2) m=1 L'algorithme comprend ensuite une phase d'itérations sur l'indice k de trames (dite boucle externe), de k=0 à un critère d'arrêt. A chaque itération k, l'algorithme comprend une phase d'itérations sur l'indice i d'émetteurs (dite boucle interne), de i=1 à K. A chaque itération i, l'algorithme calcule l'indice de sous trame t = kK +i et met à jour le vecteur de puissance de l'émetteur i : Pi,m(t) = où est (3) un paramètre à régler de manière à ce que la somme des puissances mP m du 1 13,,,'(t -1) y, ,,,(t -1) vecteur égale la puissance max Pi. et où [x]+ = max (0, x) Le critère d'arrêt de la phase d'itérations sur l'indice k de trames est par exemple déterminé tel que : ViE {1,2,...,K},VmE {1,2,...,M} Pi',(t)- Pi',(t -1)> e (4) où 0 est un paramètre à fixer. L'algorithme se déroule de la façon suivante. A chaque itération i, un seul émetteur met à jour le vecteur de puissance d'émission allouée, puis les émetteurs émettent simultanément et réceptionnent ensuite les SINR respectifs remontés par les récepteurs. L'algorithme IWFA ne peut être utilisé lorsque l'ensemble des configurations de l'émetteur est fini.. L'IWFA ne peut être utilisé dans le cas où l'émetteur a un nombre fini de possibilités d'allocation/contrôle de puissance c'est-à-dire le cas où le niveau de puissance total ou par bande est discret. En effet, tout comme la solution dite de water-filling, l'IFWA repose sur une opération de différentiation (dérivée) par rapport une quantité continue, cette opération ne peut être réalisée par rapport à une quantité discrète. De plus, l'IWFA est un cas particulier de dynamique séquentielle de meilleure réponse connue pour ne pas converger lorsque l'espace d'optimisation est discret. Une telle contrainte n'est pas compatible des réseaux de télécommunications dits 3G et de ses versions ultérieures normalisés au 3GPP et connus sous les acronymes UMTS, LTE, LTE advanced,...,imposent que les niveaux de puissance soient discrets. L'algorithme IWFA ne converge que pour des fonctions d'utilité très particulières. En effet, la convergence est garantie pour le scénario des récepteurs co-localisés et des fonctions d'utilité de type débits de transmission de Shannon dont l'expression correspond à l'équation (5). Il peut cependant être utilisé pour n'importe quelle fonction d'utilité mais sans aucune garantie sur les performances obtenues. En effet, la formule (3) de mise à jour permet d'assurer la convergence pour les fonctions d'utilité suivantes : ui = Ilog2(1+ (5) F7,1 uniquement si le SINR correspond au cas d'un canal à accès multiple multi-bande ou multi-voie i.e., au cas d'un système de communication à K émetteurs et un seul récepteur. Pour le cas à K récepteurs non co-localisés, la convergence n'est assurée pour les fonctions d'utilités selon la relation (5) que lorsque les niveaux d'interférence sont faibles.. Exposé de l'invention L'invention propose une technique permettant d'améliorer la coordination des émetteurs d'un système de télécommunication en déterminant l'allocation de puissance émise par les émetteurs qui minimise les interférences au niveau des récepteurs. Ainsi, l'invention a pour objet un procédé de coordination d'un émetteur avec K-1 émetteurs, les K émetteurs étant en communication respectivement avec K récepteurs via un lien de communication décrit pour une bande m de fréquence entre l'émetteur j et le récepteur i par des gains de canaux de transmission directs lorsque j = i et des gains de canaux de transmission interférents lorsque j# i, i E {1, ..., K}, j E {1, ..., K}, m E {1, ..., M}, K>1, chaque émetteur disposant d'estimations du gain de canal direct et des gains de canaux interférant son récepteur. Le procédé comprend : une phase dite d'échange pendant laquelle les K émetteurs émettent simultanément respectivement K séquences de T2 sous trames de respectivement T2 niveaux de puissance par bande m codant les gains de canaux estimés et pendant laquelle les K émetteurs réceptionnent après chaque sous trame des estimations de SINR basées sur des observations d'un signal reçu remontées par les K récepteurs à leur émetteur respectif, T2 > 1. Il existe des techniques connues permettant à un émetteur i de connaître l'estimation du gain de canal direct gii,m et des gains de canaux interférant son récepteur i, j # i. Lors de la phase d'échange les émetteurs échangent leurs gains de canaux estimés en les encodant selon un codage de niveaux de puissance. Le coefficient gii,,' giLmpeut modéliser aussi bien le coefficient d'affaiblissement lié à la propagation (path loss) que le coefficient d'affaiblissement lié aux évanouissements lents (shadowing) ou rapides (fast fading). Un émetteur émet donc successivement T2 sous trames dont les niveaux de puissance associés codent les estimations des gains de canaux dont l'émetteur a la connaissance. Le récepteur remonte une estimation d'un rapport signal à interférence plus bruit (SINR) à son émetteur à la fin d'une sous trame. L'estimation est effectuée à partir d'observations du signal reçu pendant que les K émetteurs émettent simultanément une sous trame. L'estimation du SINR remontée peut tout aussi bien être exprimée sous la forme d'une estimation d'une puissance utile et d'une estimation d'une puissance d'interférence comprenant le bruit, que sous la forme d'une estimation du rapport de ces puissances ou sous toute autre forme permettant d'en déduire le SINR. L'estimation est généralement bruitée notamment compte tenu que l'information remontée par la voie de retour est quantifiée. L'estimation remontée peut être le résultat d'un moyennage. Par exemple dans le cas d'une bande m qui regroupe plusieurs sous porteuses bande étroite le moyennage peut consister à effectuer l'estimation sur l'ensemble ou un sous ensemble des sous porteuses. Un émetteur dispose ainsi à la fin de la succession de T2 sous trames de T2 estimations d'un rapport signal à interférence plus bruit évalué par son récepteur. Après la dernière sous-trame, l'émetteur i dispose d'un système de T2 équations avec T2x(K-1) inconnues qui sont les puissances émises par les K-1 émetteurs et reçues par son récepteur. L'estimation des niveaux de puissances émises par les autres émetteurs peut consister pour l'émetteur i à choisir la séquence la plus vraisemblable au sens du Maximum de vraisemblance. Connaissant les puissances émises par les autres émetteurs, l'émetteur peut déterminer les gains de canaux correspondant par décodage des niveaux de puissance des puissances émises. Finalement à la fin de la phase d'échange, chaque émetteur connait tous les gains de canaux. Ainsi chaque émetteur peut déterminer sa puissance d'émission selon un critère de minimisation d'interférence globale pour l'ensemble des K émetteurs. Ainsi, le procédé selon l'invention converge vers une allocation de puissance dont l'efficacité est globale c'est-à-dire évaluée au niveau du système contrairement à l'algorithme IWFA qui converge vers une allocation de puissance à efficacité locale c'est-à-dire évaluée au niveau d'un émetteur. Contrairement à l'algorithme IWFA, le procédé selon l'invention ne repose pas sur une méthode itérative ce qui permet d'éviter le problème de convergence de l'algorithme IWFA. En outre, l'invention est compatible de toutes sortes de fonction d'utilités dont le choix dépend uniquement des politiques de contrôle de puissance et des gains de canaux du système de communication considéré. Selon un mode de réalisation, l'estimation des gains de canaux est effectuée pendant une phase d'acquisition de gains de canaux de transmission sur une bande m de fréquence entre l'émetteur j et le récepteur i pendant laquelle les K émetteurs émettent simultanément respectivement des séquences de T1 sous trames de respectivement T1 niveaux de puissance déterminés et pendant laquelle les K émetteurs réceptionnent après une sous trame des estimations d'un SINR remontées par les K récepteurs à leur émetteur respectif. Selon ce mode de réalisation, le procédé comprend une phase d'acquisition particulière. Selon cette phase d'acquisition, les émetteurs participant à la coordination émettent une succession de T1 niveaux de puissance qui sont l'équivalent d'une séquence d'apprentissage. Ces séquences d'apprentissage sont connues des émetteurs, c'est-à-dire qu'un émetteur connait les niveaux de puissance émis successivement pendant les T1 sous trames par chacun des K-1 émetteurs. Un récepteur estime Ti fois le rapport signal à interférence plus bruit (SINR) reçu. Le récepteur remonte à son émetteur à la fin d'une sous trame, la valeur estimée du SINR. Un émetteur dispose ainsi à la fin de la séquence d'apprentissage de T1 estimations du rapport signal à interférence plus bruit estimé par son r , récepteur. Le SINR peut être exprimé sous la forme générale suivante : y = f ({P, mi ,tg K ) t=i m j_i ou 21m = f ({ P m 1,_1 , E { gai mi ) ou sous la forme : 1+fi,mP ,,m gii,mPi,m . L'expression du SINR j#i dépendant du canal direct entre l'émetteur i et son récepteur i, des canaux interférents le récepteur i et des puissances d'émission des K émetteurs, l'émetteur i dispose d'un système de Ti équations avec K inconnues qui sont les gains de canaux pour j E {1, ..., K}. L'émetteur peut donc obtenir une estimation des K gains de canaux en résolvant le système d'équations. Si la matrice des puissances émises est inversible, il est aisé de résoudre le système d'équations en inversant cette matrice. Selon un mode de réalisation, le nombre Ti de sous trames est égal à K. Ce choix judicieux permet de minimiser le nombre de sous trames émises et le nombre de calculs. Selon un mode de réalisation, le nombre Ti de sous trames est égal à K-X avec X le nombre de gains de canaux connus d'un émetteur. Selon un mode de réalisation, les rapports signal à interférence plus bruit remontés par les récepteurs sont normalisés. Ce mode a pour avantage de simplifier le procédé en évitant de considérer le bruit observé par les récepteurs lors de l'estimation des gains de canaux et lors de la phase d'échange. Selon ce mode, le rapport signal à interférence plus bruit peut s'exprimer sous la forme suivante : mP 7_m n avec g if' = . Ce mode réduit avantageusement le coût de signalisation puisqu'il ± ni:Pm (Fi j#i permet d'éviter de remonter les variances de bruit 6,2 aux émetteurs. Selon un mode de réalisation, la détermination des niveaux de puissance de la séquence de Ti sous trames d'un émetteur est effectuée sous contrainte que la matrice S soit pseudo-inversible au sens de Moore Penrose avec une probabilité proche de un, cette matrice S i étant : -2,,,,'(1)Pt',(1) -y-,,,'(1)P2,',(1) P'(1) -y-,,,'(1)PK,',(1) (2) -y,,,'(2)P2,,,, (2) ... Pm (2) ... (2) PK (2) Si ril = -7i,',(T1)P2',(T1) 1 (T1) -yi,',(T1)PK,',(T1)) avec y,,,,(t) l'estimation du rapport signal à interférence plus bruit remonté par le récepteur i à l'émetteur i, Pm(t) les niveaux déterminés de puissance de la séquence de Ti sous trames, tc 11,2,...,T11. Ce mode de détermination des niveaux de puissance est avantageux en ce que la matrice S est pseudo inversible et donc la résolution du système d'équations à K inconnues est aisé. Les K gains de canaux peuvent être estimés selon la relation suivante : g- Ki,m 1 - = (s,s m Selon un mode de réalisation, la détermination des niveaux de puissance de la séquence de T1 sous trames d'un émetteur est effectuée selon un tirage aléatoire jusqu'à ce que la matrice ni', soit de rang plein : ( Pm(1) P2, (1) ... Pm(1) ... PK, (1) avec P (t) les niveaux ni,m = Pi (2) P2, (2) ... Pm(2) PK (2) _Pi,,,,(T1) P2, (T1) ... P,,m(T1) ... PK m (T1) j déterminés de puissance de la séquence de Ti sous trames, t c 11, 2, ..., T11 . Ce mode de réalisation permet de déterminer simplement les niveaux de puissance des séquences de Ti sous trames de chaque émetteur. Selon un mode de réalisation, le codage en niveaux de puissance est un codage source-canal conjoint de type spatio-temporel dont la dimension spatiale correspond au nombre M de canaux m et la dimension temporelle correspond au nombre T2 de sous-trames. Ce mode de réalisation permet de bénéficier d'un codage optimal. Selon un mode de réalisation, le codage en niveaux de puissance est un codage source-canal séparé. Ce mode de réalisation permet de bénéficier d'un codage moins complexe qu'un codage source-canal conjoint. Selon un mode de réalisation, l'estimation des rapports signal à interférence plus bruit est effectuée avant filtrage, ce filtrage étant associé à une détection du signal reçu. Une alternative à ce mode de réalisation consiste à ce que le récepteur estime le SINR après le filtrage associé à la détection du signal reçu. L'estimation avant filtrage a pour avantage d'être indépendante des caractéristiques du filtre qui peuvent dépendre des puissances des interféreurs. L'invention a en outre pour objet un émetteur destiné à mettre en oeuvre un procédé de coordination avec K-1 émetteurs, les K émetteurs étant en communication respectivement avec K récepteurs via un lien de communication décrit pour une bande m de fréquence entre l'émetteur j et le récepteur i par des gains de canaux de transmission directs lorsque j = i et des gains de canaux de transmission interférents lorsque j # i, l'émetteur disposant d'estimations du gain de canal direct et des gains de canaux interférant son récepteur. L'émetteur comprend : - une unité de traitement pour coder les gains de canaux estimés en une séquence de T2 sous trames de respectivement T2 niveaux de puissance par bande m, - un module d'émission de la séquence, - un module de réception de valeurs estimées d'un SINR remontées par le récepteur après chaque sous-trame. L'invention a en outre pour objet un système de télécommunication comprenant K émetteurs selon l'objet précédent. Selon une implémentation préférée, les étapes du procédé de coordination selon l'invention sont déterminées par les instructions d'un programme incorporé dans un circuit électronique telle une puce elle-même pouvant être disposée dans un dispositif électronique tel une entité de communication (station de base, point d'accès, etc). Le procédé de coordination selon l'invention peut tout aussi bien être mis en oeuvre lorsque ce programme est chargé dans un organe de calcul tel un processeur ou équivalent dont le fonctionnement est alors commandé par l'exécution du programme En conséquence, l'invention s'applique également à un programme d'ordinateur, notamment un programme d'ordinateur sur ou dans un support d'informations, adapté à mettre en oeuvre l'invention. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable pour implémenter un procédé selon l'invention. Le support d'informations peut être n'importe quel dispositif capable de stocker le programme Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM, une clé USB ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par un disque dur. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. D'autre part, le programme peut être traduit en une forme transmissible telle qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Liste des figures D'autres caractéristiques et avantages de l'invention apparaîtront lors de la description qui suit d'exemples particuliers faite en regard de figures annexées données à titre d'exemples non limitatifs. La figure 1 est une représentation schématique du lien entre deux émetteurs et deux récepteurs, ce lien étant décrit par des gains g Jon de canaux de transmission directs lorsque j = i et des gains g Jon de canaux de transmission interférents lorsque j # i entre l'émetteur j et le récepteur i, i E {1,2}, j E {1,2}, La figure 2 est un schéma des trois phases d'un procédé selon l'invention selon un mode de réalisation particulier. La figure 3 est un schéma de la structure simplifiée d'un émetteur selon l'invention. La figure 4 représente des résultats de simulation comparant la méthode IWFA de l'art antérieur et le procédé selon l'invention. Description d'un mode de réalisation de l'invention Les gains/coefficients des canaux sont supposés soit constants soit constants par morceaux. Dans le dernier cas, les gains de canaux sont donc stationnaires sur une durée Tc appelée durée de cohérence. Le procédé selon l'invention opére sur une durée T inférieure ou égale à cette durée de cohérence, T Tc, selon un mode privilégié, de même que l'algorithme IWFA de l'art antérieur. Dans une prise en compte d'évanouissements rapides, les canaux correspondant à une sous bande m sont supposés stationnaires sur la durée T avec T Tc ce qui permet d'omettre l'indice temporel (ou du symbole OFDM dans le cas d'une émission multi porteuses) pour les coefficients du canal. On note C le nombre de porteuses consécutives définissant la Bande de Cohérence. Soient les signaux reçus ou observations disponibles correspondant à une sous-trame de durée To T T (To symboles OFDM dans le cas d'une émission multi porteuses ) en) (k) E CNR par le récepteur i équipé de NR antennes pour la sous bande m comprenant N porteuses bandes étroites, n étant l'indice de numérotation des porteuses à l'intérieur de la sous bande m, et k l'indice du symbole OFDM (les observations/échantillons utiles pour estimer le SINR sont liées à l'ensemble des indice n, k disponibles en temps et en fréquence pour un récepteur i et une sous bande m), ils peuvent s'exprimer sous la forme : en) (k)- hçi m) .s,!nm)(k)+Ihi(ni),,si(nm) (k)+w,!nm)(k) V n =0,- ,N -1V k =0,- ,T,' -1 (6) jai (n ou hji)m E (CNR est le canal à entrée unique et sortie multiple de l'émetteur j vers le récepteur i pour la sous-bande m, où w(n) (k) est un échantillon de bruit blanc gaussien additif de variance 62 etoù est la puissance émise associée au signal émis s,(,n)(k) qui est supposée constante pour toutes les porteuses de la sous bande m, ei = 1, . , K . Sans perte de généralité, il est supposé que N est divisible par C, c.a.d., Q=N/C. Il vient que la sous bande m peut être divisée en Q blocs où les coefficients du canal sont approximativement constant (Q-Block Fading en anglais). Soit n = Cq +c , la division entière de n par C ayant pour dividende q et reste c avec c variant de 0 à C-1 par définition, l'expression (6) peut s'écrire sous la forme : eq±`)(k)- h,'qm, s;`,7q7+0(k)+Ilii(7),,si(c:7±`)(k)+w,(r)(k) q =0,- ,Q- 1 et c = 0,...0 -1 (7) J#1 1ci où h(q) =-Ih(cq+c) avec q 0,- --,Q -1 i'm Cc=0 Le SINR associé à la sous bande m pour le récepteur i est estimé sur la base des observations disponibles y(ine) (k) sur une durée d'observation To d'une sous-trame (To Le SINR remonté par le récepteur peut être pris en compte avant ou après filtrage. Après filtrage, le SINR prend en compte le filtre de réception associé au détecteur du signal reçu. Dans ce cas, il est considéré un filtre linéaire t(igni) E CNR appliqué en réception sur les échantillons n = Cq + c avec c variant de 0 à C-1. La sous bande m peut être définie par les Q SINRs théoriques, avec 1« l'opérateur transposé conjugué : v(q) _ / i ,m CF (fiqr)n fi(qr)n j#i La sous bande m peut être définie par l'unique SINR théorique équivalent : (f(q) h(q) i,m ii,m 2 Pi,m (Éq, ti,m ji,m 2 Plrn Q Q g-0 rtm= Ti 2 1 1 r; 1( f(qt,± IP/ m g-0 j#, Q g-0 2 fiimJ hjim 2 (f,q, ii,m J ii,m Q-1 (fi(ig,r1 h(qi) rn q=0 Q 1 62 (fi (q)ft(tq) m m g-0 Q-1 fiimJ ji,m 0-i2I(e) ) (q) ,M ,M q=0 gii,mPi,m 1+ jimP]m j#i Q-1 2 ( fi(ig,r)n (q 2 1+ IP. q=0 j,m Q-1 2 (8) Q-1 et Vj q- Où g m e 0-,21(tqrn) f/grn) q=0 q=0 d ( f.(.q) f.`.q) I _0 II,M II,M (4 E R+ sont définis comme les «,m = Q-1 (eq, )t iim ii,m 2 canaux équivalents (appartenant à l'ensemble des réels Ife ) entre l'émetteur j et le récepteur i pour la bande m. Dans le cas où AIR =1 et Q=1, il vient en posant fit) = hi(,°,7," (filtre adapté) : 2 te 2 gll,,, = 2 E Ik+ et Vj#i i'm2 E CF, CF, Le SINR équivalent estimé au niveau du récepteur i pour la sous bande m peut être estimé dans un mode privilégié comme suit. Les échantillons in) (k) après filtrage s'expriment par la relation : 12 en) ,(k)= t;qm)th(ilq)s;nn),(k)+If,;qm)th(;),,,s(inm) (k)+ f,;qm)tw;nn),(k) avec n=Cq+c j=i Une estimation du SINR équivalent basée sur la connaissance des échantillons 4nm) (k) (après filtrage) et des symboles transmis s,(cmq+c)(k) (symboles pilotes) sur la durée d'observation T0, ainsi que du canal estimé i^(cqm+c) (calculé à partir des symboles pilotes) s'exprime sous la forme : 9> ,m = 1 , 1(z i(c:±`) (k) - fii,m hll,m si,m (k)) (zi( ,cmq±`)(k)- fi(i7m)tili(Çqm+c) si( Cmq+c)(k)) (q)t (Cq+c) (Cq+c) ToCQ (f(q) le) ii,m ,m 1 isi(q)t f(q) z(Cq+c)tios(Cq+c)t k) ii,m ii,m i,m 1 i,m ToCQ 1'; P Q q=0 2 (9) 2 Q-1 t l 2 0- Q1( figt)// figt)// q=0 j#1 q=0 (f.(q) h(q) ii,m ji,m Dans le cas d'une estimation avant filtrage, l'estimation du SINR équivalent peut être réalisée à partir des échantillons avant filtrage y(in) (k) sur la durée To de la manière suivante en supposant les C symboles transmis scq fiq ±`)(k) et le canal estimé +c) connus : ,m 7,,m 1 1:1(q)t ',(Cq+c) s(Cq+c)t k) ,m Y i,m i,m ToCQ 1 V ((Cq+c) _ci(Cq+c),(n) (,,(Cq+c)/LA ci(Cq+c),(n) LA) (10) Y m 1'1 nnii,m / Y i,m "ii,m ToCQ Qc2-1 ht )t r 1 vc2-1(h(q) 1 h(q) m f - -ji,m j#i Q q=0 10 Soit le SINR théorique associé: 1 c2-1 h(q) ,m q-0 ri,. = 1 Q-1 62 +ri , rn -101() ) W.?) Q ',. gii,mPi,. 1+ 1g ,Pj 1 Q-1 1 C2-1 )t ii,m ji,m Q(4_0 ,m ,m Ît où gr , E le et Vj # i g ji, q-° E Ik+ sont définis comme les 2 sgi o_i2 canaux équivalents (appartenant à l'ensemble des réels Ife ) entre l'émetteur j et le récepteur i pour la bande m. Dans le cas où AIR =1 et Q=1, il vient : e) 2 1/(0) ii'm2 E le et Vj # i g 2 " E 6i 6i Une coordination sur une durée très longue, supérieure à la durée Tc de cohérence du canal, permet de prendre en compte des évanouissements lents. Les coefficients g ji,n, peuvent alors être moyennés (ou filtrés dans le temps) pour accéder à une échelle de fluctuation plus lente du type « shadowing ». Le SINR équivalent s'écrit formellement : Ym (11) j#i où E tg jim1 est la moyenne stochastique (ou temporelle par ergodicité) du coefficient gjim sur la distribution des évanouissements rapides. Finalement, l'expression du SINR équivalent peut se mettre sous la forme d'une fonction f des canaux équivalents ( g ou E tg jim1) de l'émetteur j vers le récepteur i pour la sous-bande m et des puissances d'émissions des K émetteurs 1 m : yi = f ( ,{g ) ou yim = f(fri )1( Etg J'A (12) expression dans laquelle g est supposé être un réel non-négatif (ou quantifié avec un nombre suffisant de bits) et dans laquelle la puissance 1 m émise sur le canal m par le point d'accès i appartient à un espace discret. L'invention est décrite en détail en considérant à titre d'illustration un système de communication sans fil dit système WiFi comprenant K point d'accès et K stations qui communiquent respectivement avec les points d'accès, K>1. Selon le standard IEEE 802.11n, chaque point d'accès émet sur un canal qui occupe une bande B autour d'une fréquence porteuse, par exemple B=20MHz occupée par 52 sous porteuses OFDM. Selon le standard IEEE802.11n, on distingue dans certaines configurations 13 fréquences porteuses WiFi définissant 13 canaux décalés de 5MHz. Le nombre de canaux pouvant être sélectionnés est donc restreint. Dans cette configuration, il y a donc une forte probabilité de sélection du même canal WiFi ou de canaux voisins pour des points d'accès différents et donc d'interférence entre canaux utilisés par des points d'accès géographiquement voisins. Deux points d'accès sont voisins si une même station parmi les K stations détecte les deux points d'accès. Pour exemple, en prenant comme hypothèse que le canal WiFi est déjà sélectionné par les points d'accès, les M canaux pris en compte par le procédé sont alors ceux correspondant aux sous porteuses OFDM. Le procédé selon l'invention se propose donc de déterminer une allocation de puissance pour les M=52 sous porteuses pour les K points d'accès du système de télécommunication qui émettent sur le même canal WiFi à destination respectivement de K stations WiFi. Pour chaque canal m parmi les M canaux associés aux M=52 sous porteuses, la qualité du lien de communication entre le point d'accès i E 11, 2, ..., K} et la station j E 11, 2, ..., K} est représentée par le coefficient 0. Dans ce cas, N=1 (porteuse) dans l'équation (6). Selon un autre exemple, le procédé a pour objectif de participer à la sélection du canal WiFi par les points d'accès. Les M canaux pris en compte par le procédé sont alors ceux correspondant aux fréquences porteuses WiFi qui définissent 13 canaux décalés de 5MHz. Dans ce cas, N=52 porteuses dans l'équation (6). Le procédé selon l'invention se propose donc de déterminer une allocation de puissance pour les M=13 canaux pour les K points d'accès du système de télécommunication qui émettent sur un des 13 canaux WiFi à destination respectivement de K stations WiFi. Pour chaque canal m parmi les M=13 canaux, la qualité du lien de communication entre le point d'accès i E 11, 2, ..., K} et la station j e 11, 2, ..., KI est représentée par le coefficient 0 . Le SINR est supposé avoir une structure du type correspondant à celle d'un récepteur sans annulation d'interférence, son expression générale (12) peut être donnée par une des relations (8)-(11). Toutefois, le procédé peut tout aussi bien être utilisé avec des récepteurs à annulation d'interférence. Dans ce cas, seule l'expression du SINR est modifiée pour prendre en compte les interférences annulées ; c'est uniquement le nombre de signaux interféreurs qui change dans le dénominateur du SINR. Le système peut être schématisé par la figure 1 lorsque K=2 et M=2. Les deux émetteurs sont ceux des deux points d'accès AP1 et AP2. Les deux récepteurs sont ceux des deux stations STA1 et STA2. Les deux canaux ou bandes considérés sont désignés par Chl et Ch2. Pour un K quelconque, chaque émetteur APi dispose d'une estimation du gain de canal direct et des gains de canaux interférant son récepteur STAi, c'est-à-dire des gji,m pour tout j E {1, ..., K} et pour tout me 11, 2, ..., MI . L'estimée du gain de canal est notée giu..' La connaissance des gains de canaux peut être acquise selon plusieurs modes de réalisation. Selon un premier mode de réalisation du procédé, l'acquisition par l'émetteur i de la connaissance des gains de canaux pour tout j E fl, ...,KI et pour tout me 11, 2, ..., M} est effectuée selon des méthodes connues d'estimation des gains de canaux. Selon un deuxième mode de réalisation du procédé, l'émetteur i est supposé ne pas connaître un certain nombre L de canaux incidents à son récepteur i, ce nombre L étant inférieur ou égal à K. L'acquisition est alors effectuée pendant une phase Phl dont le déroulement temporel est illustré sur la figure 2. Pendant cette phase Phl, les K émetteurs émettent simultanément respectivement K séquences. Une séquence est une succession de Tl sous trames (de durée To) de respectivement Ti vecteurs de niveaux déterminés de puissance, la dimension d'un vecteur étant égale au nombre M de bandes. Ces séquences sont comparables à des séquences d'apprentissage, elles sont connues des émetteurs. Après l'émission d'une sous trame, chaque émetteur i réceptionne par bande m le SINR normalisé ou pas remonté par son récepteur i. Cette phase Phl dite d'acquisition permet ainsi à l'émetteur i d'obtenir Ti estimations du SINR évalué par son récepteur pour une bande m donnée. T1 étant strictement supérieur à 1 un émetteur i réceptionne au moins deux estimations d'un SINR remonté par son récepteur pour une bande m donnée. L'indice t de sous trame varie de 1 à Ti, t E {1, 2, ..., Tl } . Les puissances déterminées successives d'une séquence émise sont notées : Pim(t) t E {1, 2, ..., T1 } . Les SINR estimés et remontés à chaque sous-trame t par la station i sont notés : y, m (t) . En considérant que le SINR est normalisé, l'émetteur i résout alors le système de Ti équations par bande m : P,m(1)gii,m ,m(1)yi,m(1)g (1) ;#i (14) p,m(T1)gj,m - r3j,,'(T1)7,,,'(T1)gii,m = (T1) J#i à L inconnues qui sont les gains de canaux gii,m pour j E {1, ..., L} pour obtenir une estimation de ces L gains de canauxeit,m. Le choix du nombre Ti tel que T1K assure que le système puisse être résolu selon l'équation (17). Un choix de Ti quelconque est possible pour la résolution de (14) par exemple si une estimation plus grossière est permise ou si des a priori sont disponibles. Par exemple, en TDD le coefficient g, m peut s'avérer disponible à l'émetteur i, ce qui élimine une inconnue et dans ce cas il est possible de choisir T1=K-1. Plus généralement, si l'émetteur i connait X coefficients alors il est possible de choisir T1=K-X. En outre, si la matrice des puissances émises est inversible, il est aisé de résoudre le système d'équations en inversant cette matrice. Cette phase Phl dite d'acquisition permet d'obtenir une estimation d'au maximum K gains de canaux de transmission sur une bande m de fréquence entre le point d'accès j et la station i lorsque T1>K. Selon un mode de réalisation, la détermination des niveaux de puissance de d'apprentissage d'un émetteur est effectuée sous contrainte que la matrice Sim suivante inversible au sens de Moore Penrose, pour un canal m, avec une probabilité proche de un : -2,,,,'(1)/1,,'(1) -2,,,,'(1)P2,,'(1) Pm(1)-2,,,,'(1)P',,'(1) (2)/1 m (2)-Yi,m (2)P2,', (2) ... P m (2) ...-Yi,m (2)PK,,,, (2) Sim = la séquence soit pseudo- (15) (T1) -yi'(T1)P2'(T1) ... (T1) ...-Yi,m(T1)PKm(T1)) Selon un mode de réalisation, les niveaux de puissance Pi' (t), t e 11, 2, ..., T11 déterminés des séquences d'apprentissage des émetteurs i forment une matrice s-2, ,' de rang plein : ( p'(1) P2'(1) p'(1) PK'(1) i ,m = p' (2) P2 ,m (2) ... Pi m (2) ... PK m (2) _Pin(T1) P2 (T1) ... (T1) ... PK m (T 1) j Cette détermination peut consister à tirer aléatoirement chacune des séquences de puissances jusqu'à ce que la matrice n m soit de rang plein. Les gains de canaux peuvent alors être estimés selon la relation suivante : (16) g li,m ( g Ki ,m j y, (1) \1 - (17) =i , mSi ,m ) S X M (T1) j où Si, désigne la matrice transposée de Sim . Lorsque le niveau de bruit sur le SINR tend vers zéro ou lorsque T1 devient grand alors les gains estimés g tendent vers leurs vraies valeurs g ji', . Le procédé comprend une phase Ph2 dite d'échange qui comprend une étape de codage, une étape d'estimation et une étape de décodage. Pendant la phase Ph2 d'échange, les K émetteurs émettent simultanément respectivement des séquences de T2 sous trames. Une séquence est une succession de T2 sous trames (de durée To) de respectivement T2 vecteurs de puissance codant les gains estimés g . Les vecteurs successifs de puissances codées d'une séquence émise sont notés : P, (t) P i(t) = (Pi i(t), ..., Pm (t)) . L'indice t de sous trame varie donc de 1 à T2, t e {1, 2,..., T2} , T2 1. Pour l'émetteur i , le codage des niveaux de puissances consiste, dans sa formulation générale, à coder la connaissance de cet émetteur en termes d'état de canal en une suite de vecteurs de niveaux de puissances : Si pMT2 si 1-> (Lji (T1+1), . , P (T1+ T2)) où si E Si représente la connaissance que l'émetteur a sur le super-vecteur de l'état global du canal [In' 212,'- - -' IKK1 avec g =(gij et P représente l'alphabet auquel appartiennent les niveaux de puissance. Le codage O est de manière générale un codage source-canal conjoint et spatio-temporel. La dimension spatiale est donnée par M le nombre de canaux alors que la dimension temporelle est donnée par le nombre T2 de sous-trames pendant lesquelles les connaissances individuelles sont échangées entre les émetteurs. Le codage O peut tout aussi bien être un codage source-canal séparé, par sous bande m, i.e. ° Oism par sous bande m. Dans ce cas, la fonction Osm de codage source est choisie comme une opération de quantification vectorielle opérant sur le vecteur (gi,''..., uni) E R qui fait correspondre par l'emetteur i la connaissance des canaux incidents au recepteur i à un indice de cellule de quantification vectorielle k (L est le nombre de cellules) : 17 (18) R+K gli,m,---, g Ki,m) {1, . , L1 H> k VE - (19) Le choix des cellules de quantification est fait par utilisation de l'algorithme de Lloyd-Max [2]. est l'opération de codage canal qui fait correpondre à un indice de cellule de quantification vectorielle une suite de puissances : {1,---,L} pT2 (20) k H> (P,', (T1+1) , , m (T1 + T2)) Soit p la cardinalité de P, le nombre de sequences possibles de puissance est donc pT2 =z . Le codage canal consiste à selectionner un sous ensemble C, de cardinalité 2k , de ces séquences ce qui k définit un codage de rendement r=-. Le nombre de cellules de quantification L donne la distorsion n dûe à la quantification. Selon un mode de réalisation, la distorsion est réduite en choisissant un nombre maximal de cellules, c'est-à-dire L = pT2 2n . Le codage canal O c consiste simplement à faire correspondre à im ,m - Yi chaque indice de cellule une séquence de puissance. On a ainsi un codage de canal de rendement 1. Le représentant de chaque cellule de quantification est choisi de la manière suivante. Les cellules sont indicées de 1 à p T 2 , p étant la cardinalité de P. Pour chaque cellule, on peut calculer la probabilité que le vecteur de gain (git ,..., g Ki, )G R appartienne à celle-ci. Il est alors possible d'ordonner ces cellules selon leur probabilité d'utilisation. Le procédé associe à la cellule de probabilité la plus élevée la séquence de niveaux de puissance (Pim(T1+1),..., (T1+ T2)) dont la somme des composantes v t=T1+1Tl+T 2 (t) est la plus élevée, c'est-à-dire la séquence des niveaux maximum de puissances sur la s bande. Le procédé considère alors la deuxième cellule, soit de probabilité maximale, soit de probabilité juste inférieure, et lui associe la séquence de niveaux de puissance dont la somme est immédiatement inférieure et ainsi de suite. Ainsi l'ensemble des suites de niveaux de puissances IT1+T2 (T 1 + T 2)) se retrouve ordonné de manière à ce que la somme Lt=T1+1P m (t) soit décroissante. Après l'émission d'une sous trame, chaque émetteur i réceptionne par bande m le SINR normalisé ou pas remonté par son récepteur i. Les SINR estimés et remontés à chaque sous-trame sont notés : yi (t) . Un émetteur dispose ainsi à la fin de la succession de T2 sous trames de T2 estimations du SINR déterminées par son récepteur par bande m. Selon un mode de réalisation, les canaux sont supposés quasi statiques, ils ne varient pas sur la durée de la transmission. L'émetteur i possède des estimées des canaux dont son SINR dépend et les niveaux de puissances sont discrets. L'étape d'estimation se déroule alors de la manière suivante. Connaissant Pi (t) l'émetteur i calcule la quantité : (t)) = ,ril P1,111 (t) (21) 1+1g mPi,m(t) où x. (t) = (t) - - - P. (t) P (t) - - - P (01 est le vecteur des puissances associées aux m 1,m K ,m autres émetteurs. On définit le vecteur xi sur l'ensemble de la trame T2 : xi, =[Xim (T1+1),---,xim(T1+T2)1 On définit le vecteur ÿi m des SINR estimés : jri, = Efi, (T1+1), - - - , j>i,m(T1+T 2)1 L'estimation des niveaux de puissance consiste à choisir la séquence la plus vraisemblable (au sens du Maximum de Vraisemblance) c'est-à-dire la séquence qui maximise la distribution de probabilité P(Srijnxi,m)Lx,n, E Cl où [x] est égal à 1 si x est vrai et égal à 0 sinon. Le crochet de Iverson [.] impose que chaque séquence de puissance appartienne au dictionnaire C des séquences possibles (dictionnaire supposé le même pour chaque émetteur). A une séquence choisie correspond un indice de cellule, c'est-à-dire au décodage. Soit le vecteur : yim= (x (T1+1)),- - -, (x,, (T1+ T2))1 Si l'on suppose que : n avec n qui suit une distribution Gaussienne dont la covariance est un multiple de l'identité. Alors, l'estimation selon le maximum de vraisemblance s'écrit : T1+T2 2 2 (22) X. = arg min 1lYi, - yi,11 = arg min (t)- (xi (t))1 i,mI,ECK I t=T1+1 Si le nombre K d'interféreurs effectifs est faible (typiquement un à quelques unités), l'émetteur i peut en général approximer son SINR et tester de manière exhaustive les pl?T2 configurations possibles. Ainsi pour p= 8 niveaux de puissances, K =1 interféreur dominant et T2 = 3 sous-trames, il y a 512 configurations à tester. Selon un mode de réalisation, le procédé effectue une approximation du SINR en ne tenant compte que d'un nombre réduit d'interféreurs. Cette approximation découle d'une réduction de modèle. Cette réduction consiste à choisir une valeur pour K , par exemple K =1 . L'expression du SINR estimé par le récepteur i devient : g ii,mPi,m ,m 1+ g .imPi,m- Seul l'ensemble Um est considéré comme influant pour l'estimation du SINR. Cet ensemble comprend le voisinage de la liaison ii sur le canal m. S'il n'y a qu'un seul intérfereur dominant, l'ensemble Um des voisins ne comprend qu'un seul interféreur (Um est un singleton). Si les ensembles 1V m sont de cardinalités différentes alors selon un mode de réalisation, le procédé considère uniquement l'ensemble V = max(i ) I V, m I, où la notation 1.1 désigne la cardinalité d'un ensemble, pour l'expression des SINR estimés par les différents récepteurs. En reprenant l'exemple illustré par la figure 1, le nombre d'interféreur est V = 1 , il n'y a pas de besoin d'effectuer une réduction de modèle. Le nombre de canaux M = 2. Chaque émetteur est supposé avoir le choix entre seize niveaux de puissances d'émission (i.e., I P 1= p= 16) et chaque (23) coefficient/gain de canal est supposé avoir été quantifié par une procédure quelconque et peut prendre quatre valeurs notées g(I) ... g(4) tij ,,m. Pendant la phase d'acquisition, en fixant T1=2, les deux émetteurs AP1, AP2 transmettent en même temps pour chacun des deux canaux Chl, Ch2 une première sous-trame puis une seconde sous- trame qui forment une séquence d'apprentissage. Pour chaque sous-trame, les émetteurs ont connaissance des niveaux de puissances utilisés par l'autre émetteur (phase d'apprentissage). A la fin d'une sous trame, chacun des points d'accès réceptionne le SINR remonté par sa station. Pour les deux points d'accès AP1 et AP2 les SINR remontés sont respectivement les suivants : et g11,1PL1 ru= l+ g , P, 71 2 = - 21,1 -,1 gl 1 2P1,2 1+ g21,2P2,2 l p = g22,' 2 2 1 1 p g2 1 2 p2 72 = g221 1 g1 2P2 p2 A la fin de la séquence d'apprentissage en ce qui concerne le point d'accès AP1, pour le canal Chl (m=1), le système (14) linéaire d'équations à résoudre devient : 1/1,1(1)g-11,i P2,1(1)211,1(1)g21,1 = 211,1(1) Pi,i (2)g11,1 P2,1(2) 211,1(2)g- 21,1 = 211,1(2) où Pia (1) est la puissance allouée par l'émetteur AP1 sur le canal Chl pour la sous-trame t=1, P21(1) est la puissance allouée par l'émetteur AP2 sur le canal Chl pour la sous-trame t=1, etc. Les inconnues du système sontg111 et g211 , les autres quantités sont connues de l'émetteur AP1. Pour le canal Ch2 (m=2), le système (14) linéaire d'équations à résoudre devient : {P1,2 (1)g11,2 P2,2 (1)Î11,2 ( 1)g 21,2 = 211,2 ( 1) 2D P1,2 (2) g 11,2 P2,2(2) 2/1,2(2) g 21,2 = il1,2 I(2) OÙ /12(1) est la puissance allouée par l'émetteur AP1 sur le canal Ch2 pour la sous-trame t=1, P2,2 ( 1) est la puissance allouée par l'émetteur AP2 sur le canal Ch2 pour la sous-trame t=1, etc. Les inconnues du système sont 811,2 et g21,2 , les autres quantités sont connues de l'émetteur AP1. 25 Pour le point d'accès AP1 la relation (17) pour estimer les gains de canaux devient : g 11,m g 21,m = (s1,ms1,m) x avec S1,m = -r1,m(1)p2,,,(1)` -ri,,,(2)P2,,, (2) Et pour le point d'accès AP2 le système d'équations (14) devient : {P2,m (1)g 22,m - Pl,m(l)r2,m(1)g12,m = 72,in (1) Et la relation (17) pour estimer les gains de canaux devient : g 22,m (p2,,, -r2,,,(1)pi,,,(1)` = r2,.(1) avec S2m = vg12,m r2,m(2)) (2) -r2,,,(2)1tin (2)) Pendant la phase d'échange, le paramètre T2 est fixé à un, les gains estimés pendant la phase d'acquisition sont échangés pendant une sous trame. Le codage en niveau de puissance se déroule de la façon suivante. Pour tout (ij,m), les valeurs possibles pour gis m sont ordonnées de la plus faible à la plus forte : g(1) g(4) . De même les niveaux de puissances sont notés de manière à ce que la suite P(1), ---, P(16) soit une suite croissante. Pour l'émetteur AP1, le codage 0 est choisi comme une application bijective (g (£ 1) g(£2 )) 1-> P(£) pour me 11, 21, 11, ..., 41 , £2 E 11,..., 41, 11m 21,m 1,m ee {1,...,16}. L'émetteur AP2 fait de même pour les deux canaux Chl et Ch2. L'observation par l'émetteur APi du niveau de puissance de l'émettteur APj se fait au moyen de la remontée du SINR : P= Pi mg ii,m (24) L'APi détermine alors les gains de canaux de l'APj en appliquant : (g11,m, g21,m = j,m) (25) où 0;1' est l'application inverse de . Ainsi, l'API détermine Qg et l'AP2 221, gl22 g 222 détermine g111 , g21,1, g11,2 g21,2 Selon un mode de réalisation du procédé, le procédé comprend en outre une phase Ph3 dite d'allocation de puissance pendant laquelle les K émetteurs déterminent leur puissance à émettre qui satisfasse un critère de maximisation des SINR des différentes stations : K M (El arg max m(ri,m) (Pi ,, POEPMK 1=1 m=1 (26) où fi m est une fonction qui définit le critère de performance de l'émetteur i, par exemple : fi,m(x) = log(1+ x) . (27) La structure simplifiée d'un émetteur mettant en oeuvre un procédé de coordination selon l'invention est décrit en relation avec la figure 4. L'émetteur est une partie d'un point d'accès WiFi, d'une station de base d'un réseau cellulaire, d'une femto cell, etc. L'émetteur fait partie d'un système P2,m (2)g 22,m Pt!), (2) r2,', (2)g 12,m = (2) de télécommunication tel que celui représenté à la figure 2. L'émetteur est en communication avec un récepteur du système. Les liens de communication sont décrits pour une bande m de fréquence par des gains g Jon de canaux de transmission directs lorsque j = i et des gains g de canaux de transmission interférents lorsque j # i entre l'émetteur j et le récepteur i. L'émetteur dispose d'une estimation du gain de canal direct et des gains de canaux interférant son récepteur. Un tel émetteur AP1, AP2 comprend un module de mémorisation 100 comprenant une mémoire tampon RAM, une unité de traitement 101, équipée par exemple d'un microprocesseur 1.1.13 et pilotée par le programme d'ordinateur 102, mettant en oeuvre le procédé de coordination selon l'invention. A l'initialisation, les instructions de code du programme d'ordinateur 102 sont par exemple chargées dans la mémoire RAM avant d'être exécutées par le processeur de l'unité de traitement 101. Le microprocesseur de l'unité de traitement 101 met en oeuvre les étapes du procédé de coordination décrit précédemment, selon les instructions du programme d'ordinateur 102, pour émettre une trame de T2 sous trames de respectivement T2 niveaux de puissance (P,,,,(1), , Fm (T 2)) . L'unité de traitement 101 exploitent en entrée les estimations des gains de canaux et les rapports de signal à interférence plus bruit remontés par le récepteur. Pour cela, l'émetteur comprend : - un module d'émission d'une séquence de T2 sous trames de respectivement T2 niveaux de puissance par bande m, - un module de réception d'estimations d'un SINR en provenance de son récepteur. Lors de l'exécution par le processeur de l'unité de traitement 101 des instructions de code du programme d'ordinateur 102 : - l'unité de traitement 101 récupère les gains gji,m de canaux estimés, par exemple depuis le module de mémorisation 100, - l'unité de traitement 101 code ces estimations en une séquence de T2 niveaux de puissance (P, m(1), , (T2)) par bande m, - le module d'émission émet la séquence, - le module de réception réceptionne T2 estimations d'un SINR en provenance du récepteur. Selon un exemple de mise en oeuvre, le procédé selon l'invention est utilisé dans un contexte de sélection de canal par des points d'accès AP1, AP2 et des stations SAT1, STA2 compatibles d'une norme IEEE 802.11xx dite WiFi. Quelle que soit la version considérée de la norme WiFi, le récepteur STA1, STA2 peut effectuer une estimation du SINR. Le récepteur i remonte à son émetteur respectif via une voie de retour une estimation de ce SINR. La voie de retour primitive et toujours disponible est la voie ACK/NACK qui peut être utilisée pour coder le SINR sur plusieurs bits, un ACK/NACK représentant un bit. Les versions avancées du standard IEEE 802.11 telles que IEEE 802.11n et suivantes spécifient une voie de retour plus évoluée (de type quantized CSI, 3 de ces modes sont spécifiés dans la version IEEE 802.11n) qui peut être utilisée pour remonter le SINR périodiquement. A l'issue des phases d'acquisition Phl et d'échange Ph2, chaque émetteur peut tester de manière exhaustive l'ensemble des meilleures sélections de canaux. Par exemple, s'il y a M=13 canaux possibles et K=5 émetteurs actifs, cela revient à tester 135 = 371293 configurations possibles pour tout le système, ce qui est très facile à faire avec la technologie actuelle. Chaque émetteur peut par exemple choisir le super-vecteur de puissance : X [P11, PlM P2M PK1, PKM K M qui maximise le critère global de performance : G = lu. avec u; =Ilog2(1 + Y, m) sous les 1=1 m=1 contraintes C1, C2 (les SINR étant définis comme précédemment). La contrainte C1 est telle que chaque composante du vecteur ne peut prendre que deux valeurs (0 ou Pmax). La contrainte C2 est telle que le vecteur d'allocation de puissance d'un émetteur donné ne peut avoir qu'une seule composante non nulle. Selon un mode de réalisation plus particulier, le procédé prend en compte en outre la contrainte C3. La contrainte C3 est telle qu'une bande ne peut être allouée qu'à un seul émetteur. Alors le gain optimal G.; devient sous les contraintes C1, C2, C3 : K M Gmax = max G = maxIIlog (1+ gil,f(i)Pmax) avec s.t. Cl nC2nC3 signifiant sous les x s.t. CinC2nC 3 i=1 m contraintes conjointes C1, C2, C3. La maximisation se fait sur toutes les injections f de {1,...,K} vers M. 20 {1,...,M} soit sur AK = M= 154440 injections f. Choisir une telle injection revient en effet à (M - K)! choisir un vecteur de puissance dont les composantes vérifient les trois contraintes considérées. A noter que selon ce mode particulier, seule la connaissance de l'ensemble des tg" 1K Vm est nécessaire à ,m 1=1 chaque émetteur. Des résultats de simulation sont illustrés par la figure 4. Les simulations ont été faites avec les 25 paramètres suivants de configuration : - le système comprend deux émetteurs et deux récepteurs, - largeur d'un canal Et, =10 kHz, - SNR = SINR (en l'absence d'interférence) = 40 dB, - coefficients de canal : gll' =1 and gum, = 0.5 pour i # j , - critère de performance selon (5) plus moyennage sur les coefficients de canal tirés selon une loi de Rayleigh. Selon les résultats présentés sous forme de deux courbes, une pour l'algorithme IWFA, une autre pour le procédé selon l'invention, on constate que le procédé selon l'invention permet de se rapprocher de la limite ultime prise comme référence et obtenue avec un système centralisé. Le gain en débit de transmission est d'environ 60% par rapport à la solution IWFA. Références : [1] : G. Scutari, D. P. Palomar, and S. Barbarossa, "The MIMO iterative waterfilling algorithm," IEEE Trans. Signal Process., Vol. 57, No. 5, pp. 1917-1935, May 2009. [2] : Robert M. Gray, David L. Neuhoff, "Quantization", IEEE Trans. on Inf. Theory, 1998 REVENDICATIONS1. Procédé (1) de coordination d'un émetteur avec K-1 émetteurs, les K émetteurs étant en communication respectivement avec K récepteurs via un lien de communication décrit pour une bande m de fréquence entre l'émetteur j et le récepteur i par des gains (gel,,,) de canaux de transmission directs lorsque j = i et des gains (g ji,m) de canaux de transmission interférents lorsque j# i, i E {1, , K}, j E {1, , K}, K>1, m E {1, , M}, chaque émetteur disposant d'estimations du gain de canal direct et des gains de canaux interférant son récepteur, caractérisé en ce qu'il comprend : une phase dite d'échange (Ph2) pendant laquelle les K émetteurs émettent simultanément respectivement K séquences de T2 sous trames de respectivement T2 niveaux de puissance par bande m codant les gains de canaux estimés (gnon) et pendant laquelle les K émetteurs réceptionnent après chaque sous trame des estimations de SINR basées sur des observations d'un signal reçu remontées par les K récepteurs à leur émetteur respectif, T2 > 1. 2. Procédé (1) de coordination de K émetteurs selon la revendication 1 dans lequel l'estimation des gains de canaux est effectuée pendant une phase d'acquisition (Phl) de gains de canaux de transmission sur une bande m de fréquence entre l'émetteur j et le récepteur i pendant laquelle les K émetteurs émettent simultanément respectivement des séquences de T1 sous trames de respectivement T1 niveaux de puissance déterminés et pendant laquelle les K émetteurs réceptionnent après une sous trame des estimations d'un SINR remontées par les K récepteurs à leur émetteur respectif. 3. Procédé (1) de coordination de K émetteurs selon la revendication 2 dans lequel le nombre T1 de sous trames est égal à K. 4. Procédé (1) de coordination de K émetteurs selon la revendication 2 dans lequel le nombre T1 de sous trames est égal à K-X avec X le nombre de gains de canaux connus d'un émetteur. 5. Procédé (1) de coordination de K émetteurs selon l'une des revendications 2 à 4 dans lequel la détermination des niveaux de puissance de la séquence de T1 sous trames d'un émetteur est effectuée sous contrainte que la matrice S i m soit pseudo-inversible au sens de Moore Penrose avec une probabilité proche de un, cette matrice S i étant :mSi ril = -2,,,,'(1)Pt',(1) -y-,,,'(1)P2,',(1) Pm(1) -y-,,,'(1)P',',(1) (2) -y,,,'(2)P2,,,, (2) ... P m (2) ... -Yi,m (2) PK (2) -yi,',(T1)P2,',(T1) p,m(T1) -yi,',(T1)PK,',(T1)) avec y,,,,(t) l'estimation du rapport signal à interférence plus bruit remonté par le récepteur i à l'émetteur i, Pm(t) les niveaux déterminés de puissance de la séquence de T1 sous trames, tc 11,2,...,T11. 6. Procédé (1) de coordination de K émetteurs selon l'une des revendications 2 à 4 dans lequel la détermination des niveaux de puissance de la séquence de T1 sous trames d'un émetteur est effectuée selon un tirage aléatoire jusqu'à ce que la matrice ni', soit de rang plein : 1 pi (1) P2 (1) ... Pm(1) ... PK (1) ni,m = Pi,m(2) P2 ,m (2) ... /3,,m(2) ... PK m (2) avec P m (t) les niveaux Pi,m(T1) P2,m (T1) ... P,,m(T1) ... PK (T1) j déterminés de puissance de la séquence de T1 sous trames, t c 11, 2, ..., T11 . 7. Procédé (1) de coordination de K émetteurs selon l'une des revendications précédentes dans lequel le codage en niveaux de puissance est un codage source-canal conjoint de type spatio-temporel dont la dimension spatiale correspond au nombre M de canaux m et la dimension temporelle correspond au nombre T2 de sous-trames. 8. Procédé (1) de coordination de K émetteurs selon l'une des revendications 1 à 6 dans lequel le codage en niveaux de puissance est un codage source-canal séparé. 9. Procédé (1) de coordination de K émetteurs selon l'une des revendications 1 à 8 dans lequel l'estimation des rapports signal à interférence plus bruit est effectuée avant filtrage, ce filtrage étant associé à une détection du signal reçu. 10. Emetteur (AP1, AP2), destiné à mettre en oeuvre un procédé de coordination avec K-1 émetteurs, les K émetteurs étant en communication respectivement avec K récepteurs via un lien de communication décrit pour une bande m de fréquence entre l'émetteur j et le récepteur i par des gains (gnon) de canaux de transmission directs lorsque j = i et des gains (gnon) de canaux detransmission interférents lorsque j# i, i E {1, ... , K}, j E {1, ... , K}, m E {1, ... , M}, K>1, l'émetteur disposant d'estimations du gain de canal direct et des gains de canaux interférant son récepteur, caractérisé en ce qu'il comprend : - une unité de traitement (101) pour coder les gains de canaux estimés en une séquence de T2 sous trames de respectivement T2 niveaux de puissance par bande m, - un module d'émission de la séquence, - un module de réception de valeurs estimées d'un SINR remontées par le récepteur après chaque sous-trame. 11. Système de télécommunication comprenant K émetteurs selon la revendication précédente. 12. Programme d'ordinateur sur un support d'information, ledit programme comportant des instructions de programme adaptées à la mise en oeuvre d'un procédé de coordination de K émetteurs selon l'une quelconque des revendications 1 à 9, lorsque ledit programme est chargé et exécuté dans les K émetteurs destinés à mettre en oeuvre le procédé de coordination. 13. Support d'informations comportant des instructions de programme adaptées à la mise en oeuvre d'un procédé de coordination selon l'une quelconque des revendications 1 à 9, lorsque ledit programme est chargé et exécuté dans K émetteurs destinés à mettre en oeuvre le procédé de coordination.