Brevet FR2963194A1 — publié le 20120127

MÉTHODE DISTRIBUÉE D'ALLOCATION DE PUISSANCE SOUS CONTRAINTE DE QUALITÉ DE SERVICE DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne de manière générale le domaine des systèmes des télécommunications sans fil et plus particulièrement une méthode d'allocation de 10 puissance dans un tel système. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Un des problèmes affectant une communication au 15 sein d'un système de télécommunication cellulaire est celui de l'interférence générée par les autres communications de la cellule ou des cellules voisines. On distingue classiquement l'interférence inter- cellulaire due aux communications des cellules voisines 20 et l'interférence intra-cellulaire due aux communications de la même cellule où se trouve le terminal. De nombreuses techniques ont été proposées et mises en oeuvre pour réduire l'interférence intra- 25 cellulaire. La plupart de ces techniques reposent sur une allocation de ressources de transmission orthogonales, par exemple d'intervalles de temps de transmission (TDMA), de fréquences de transmission (FDMA), d'intervalles de sous-porteuses d'un multiplex 30 OFDM (OFDMA), de codes de transmission (CDMA), de faisceaux de transmission (SDMA), voire d'une combinaison de telles ressources, de manière à séparer les différentes communications d'une même cellule. Les ressources de transmission étant rares, celles-ci sont généralement réutilisées, au moins partiellement, d'une cellule à l'autre. Un module de gestion des ressources de transmission ou RRM (Radio Resource Management) est alors en charge d'allouer de manière statique ou dynamique les ressources de transmission aux différentes cellules. Il est notamment connu de réutiliser de manière statique des fréquences de transmission selon un schéma bi-dimensionnel (Frequency Reuse Pattern). Du fait de la réutilisation des ressources de transmission, une première communication entre un premier terminal et une première station de base d' une cellule peut être interférée par une seconde communication, utilisant la même ressource de transmission, entre un second terminal et une seconde station de base d'une cellule voisine. La situation est encore davantage critique lorsque les cellules sont adjacentes et les terminaux sont en bordure de cellule. Dans ce cas, les terminaux doivent émettre à pleine puissance et le niveau d'interférence est alors élevé. Pour une communication donnée, dite ici première communication, on appelle communément interférence intra-bande l'interférence induite par une seconde communication utilisant la même ressource de transmission que la première. Par opposition, on appelle interférence inter-bande celle induite par une seconde communication utilisant une ressource de transmission distincte (par exemple une fréquence de transmission voisine ou un autre intervalle de transmission) de celle utilisée par la première. La Fig. 1 représente un système de télécommunication cellulaire très simplifié, comprenant deux cellules 151 et 152. La première cellule 151 contient un premier couple de terminaux constitué d'un premier terminal émetteur 110 et d'un premier terminal récepteur 120. De manière similaire, la seconde cellule 152 comprend un second couple de terminaux constitué d'un second terminal émetteur 130 et d'un second terminal récepteur 140 Par terminal, on entend ici un terminal mobile ou une station de base, voire un terminal relais dans le cas d'un canal relayé. En particulier, on comprendra que l'on considère ici aussi bien les communications sur voie montante (uplink) que sur voie descendante (downlink). On suppose en outre que la première communication entre les terminaux 110 et 120, utilise la ou les même(s) ressource(s) de transmission que la seconde transmission entre les terminaux 130 et 140 de sorte que les deux communications interfèrent entre elles. Le traitement et la réduction de l'interférence 25 inter-cellulaire ont fait l'objet de très nombreuses recherches. La méthode de traitement la plus simple est de considérer l'interférence comme un simple bruit thermique. Cette méthode de traitement n'est toutefois 30 acceptable que si le niveau d'interférence est faible. Il est à noter que la plupart des algorithmes d'allocation de puissance (power allocation) sont basés sur cette hypothèse. D'autres méthodes de traitement permettent de réduire l'interférence en décodant le signal d'information de la (ou les) communication(s) interférente(s). Cela suppose que le terminal récepteur considéré connaisse les codes ayant servi à les coder. Parmi ces méthodes, on connaît notamment les schémas de réduction d'interférence de type parallèle ou PIC (Parallel Interference Canceller) ou série (Successive Interference Canceller), bien connus de l'homme du métier. Une autre approche classique pour réduire le niveau d'interférence est de mettre en oeuvre une méthode adaptative d'allocation de puissance (adaptive power control). Une telle méthode permet de contrôler les niveaux de puissance des différents terminaux émetteurs de manière à garantir une qualité de service prédéterminée aux différents utilisateurs. Cette qualité de service peut être mesurée selon le cas en termes de débit, de temps de latence, de taux d' erreur paquet, de couverture spatiale, etc. On appelle conventionnellement métrique de qualité de service, le ou les paramètre(s) utilisé(s) pour la mesurer. En règle générale, une communication d'un utilisateur requiert une qualité de service minimale qui est prise en compte ou négociée lors de la procédure d' admission de l'utilisateur dans la cellule. Cette qualité de service minimale s'exprime sous la forme d'une contrainte sur la métrique de qualité de service latence inférieure à un seuil, débit supérieur à un minimum garanti, etc.. L'allocation de puissance est alors effectuée de manière à respecter la contrainte sur la métrique de qualité de service. Enfin, l'allocation de puissance peut être prise en charge de manière centralisée (Centralized Power Allocation) par un noeud spécifique du réseau, NC (Network Controller), ou bien de manière distribuée (Distributed Power Allocation) au sein de chaque terminal. Le principal inconvénient de l'approche centralisée est de requérir la transmission par les terminaux d'une importante quantité d'information au noeud de contrôle, notamment celle relative aux coefficients de canal entre terminaux émetteur et terminaux récepteur ou CSI (Channel State Indicator). Une méthode distribuée d'allocation de puissance est connue de l'article de Jong-Shi Pang et al. intitulé « Distributed power allocation with rates constraints in Gaussian parallel interference channels publié dans IEEE Trans. on Information Theory, Vol. 54, N° 8, Août 2008, pages 3471-3489. Les méthodes d'allocation de puissance, distribuée ou non, font généralement l'hypothèse que l'interférence est assimilable à du bruit thermique. Toutefois cette hypothèse est bien souvent pessimiste de sorte que les puissances d'émission allouées peuvent être sensiblement supérieures à celles qui seraient nécessaires pour obtenir la qualité de service souhaitée. Cette allocation conduit par conséquent à une consommation d'énergie inutile et, pour les communications montantes, à une réduction de l'autonomie des terminaux. Un premier but de l'invention est par conséquent de proposer une méthode distribuée d'allocation de puissance qui permette de réduire les puissances d'émission des terminaux pour une contrainte de qualité de service donnée. Un second but de l'invention est de proposer une méthode distribuée d'allocation de puissance qui tienne compte des types d'interférence entre communications. EXPOSÉ DE L'INVENTION La présente invention est définie par une méthode d'allocation de puissance pour système de télécommunication sans fil comprenant au moins un premier couple de terminaux constitué d'un premier terminal émetteur et d'un premier terminal récepteur et d'un second couple de terminaux constitué d'un second terminal émetteur et d'un second terminal récepteur, une première communication entre les terminaux du premier couple utilisant les mêmes ressources de transmission qu'une seconde communication entre les terminaux du second couple, de sorte que les deux communications interfèrent entre elles, la première communication devant respecter un premier débit de consigne et la seconde communication devant respecter un second débit de consigne. Selon cette méthode on initialise les puissances d'émission des premier et second terminaux émetteur à des valeurs prédéterminées ; on effectue une pluralité d'itérations 30 jusqu'à ce qu' un critère d' arrêt soit satisfait, chaque itération comprenant une première pluralité d'étapes effectuées par le premier couple de terminaux suivie d'une seconde pluralité d'étapes effectuées par le second couple de terminaux, la première pluralité d'étapes comprenant : (al) une estimation par le premier terminal récepteur du rapport signal sur bruit ainsi que du rapport interférence sur bruit, le signal étant celui de la première communication et l'interférence étant induite par la seconde communication ; (bl) un calcul de la puissance minimale d'émission du premier terminal émetteur autorisant un premier débit de consigne, à partir des rapports estimés à l'étape précédente ; (cl) une actualisation de la puissance 15 d'émission du premier terminal émetteur à la valeur ainsi calculée ; la seconde pluralité d'étapes comprenant : (a2) une estimation par le second terminal récepteur du rapport signal sur bruit ainsi que du 20 rapport interférence sur bruit, le signal étant celui de la seconde communication et l'interférence étant induite par la première communication ; (b2) un calcul de la puissance minimale d'émission du second terminal émetteur autorisant le 25 second débit de consigne, à partir des rapports estimés à l'étape précédente ; (c2) une actualisation de la puissance d'émission du second terminal émetteur à la valeur ainsi calculée. 30 La puissance minimale d'émission du premier terminal récepteur peut être obtenue à l'étape (bl) en effectuant une première partition d'un diagramme d'interférence de la première communication en différentes zones, chaque zone correspondant à un type de traitement distinct de l'interférence, en calculant, à partir des premier et second débits de consigne, une fonction donnant le rapport signal sur bruit minimal permettant d'effectuer le traitement de l'interférence dans chacune desdites zones, et en déterminant l'image par cette fonction du rapport interférence sur bruit mesurée à l'étape (a1). La puissance minimale d'émission, Pin), du premier terminal émetteur peut être calculée à la neme itération 15 par : P(n) =21 (Ô-1(n» No g l 1 où 2l est ladite première fonction, No est la puissance 20 de bruit thermique et g11 est le coefficient de canal de la première communication, estimés par le premier terminal récepteur et e) est le rapport interférence sur bruit mesurée à l'étape (al) de la ne,ne itération. Le calcul de la puissance minimale d'émission en) 25 est effectué par le premier terminal récepteur. Alternativement, ce calcul peut être effectué par le premier terminal émetteur. De manière similaire, la puissance minimale d'émission du second terminal récepteur est obtenue à l'étape (b2) en effectuant une seconde partition d'un diagramme d'interférence de la seconde communication en différentes zones, chaque zone correspondant à un type de traitement distinct de l'interférence, en calculant, à partir des premier et second débits de consigne, une seconde fonction donnant le rapport signal sur bruit minimal permettant d'effectuer le traitement de l'interférence dans chacune des zones de la seconde partition, et en déterminant l'image par cette seconde fonction du rapport interférence sur bruit mesurée à l'étape (a2). La puissance minimale d'émission, el), du second terminal émetteur peut être calculée à la Hème itération par : PZ(n)_22(62 J No 2 Ig22 où 22 est ladite seconde fonction, No est la puissance de bruit thermique et g22 est le coefficient de canal de la seconde communication, estimés par le second terminal récepteur et 6 est le rapport interférence 25 sur bruit mesurée à l'étape (a2) de la en' itération. Le calcul de la puissance minimale d' émission peut être effectué par le second terminal récepteur. Alternativement, ce calcul peut être effectué par le second terminal émetteur. L'invention concerne également une méthode d'allocation de puissance pour système de télécommunication sans fil comprenant une pluralité de couples de terminaux, chaque couple étant constitué d'un terminal émetteur et d'un terminal récepteur, dans laquelle on alloue initialement la même puissance d' émission (P') aux terminaux émetteurs desdits couples, et pour une communication, dite première communication entre les terminaux émetteur et récepteur d'un couple, dit premier couple, on détermine la communication, dite seconde communication, entre les terminaux émetteur et récepteur du couple, dit second couple, qui induit le plus fort niveau d'interférence sur la première communication, lesdits premier et second couples formant une paire, et l'on effectue ensuite l'allocation de puissance aux premier et second terminaux au moyen de la méthode d'allocation telle qu'exposée ci-dessus. On peut alors allouer une puissance d'émission réduite (P") aux terminaux émetteurs autres que ceux appartenant à ladite paire pendant le temps où les première et seconde communications sont actives et l'on intègre alors dans l'estimation du bruit thermique l'interférence induite par les communications extérieures à ladite paire. Alternativement, on peut allouer une puissance d'émission nulle aux terminaux émetteurs autres que ceux appartenant à ladite paire pendant le temps où les première et seconde communications sont actives. Dans les deux cas, le premier couple de terminaux peut être choisi tour à tour parmi les couples de terminaux de ladite pluralité. Alternativement, le premier couple de terminaux peut être choisi parmi les couples de terminaux de ladite pluralité selon un degré d'urgence, un ordre de priorité ou un facteur de qualité de canal de transmission. BRÈVE DESCRIPTION DES DESSINS D' autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention, fait en référence aux figures jointes parmi lesquelles : La Fig. 1 représente schématiquement une situation d'interférence inter-cellulaire intra-bande 20 dans un système de communication cellulaire ; La Fig. 2 représente schématiquement un modèle du système de communication cellulaire de la Fig. 1 ; La Fig. 3 représente schématiquement un premier diagramme d'interférence pour une première 25 communication ; La Fig. 4 représente schématiquement un second diagramme d'interférence pour cette même communication ; La Fig. 5 représente schématiquement un diagramme 30 d'interférence pour une première et une seconde communications ; La Fig. 6 représente schématiquement un ordinogramme de la méthode distribuée d'allocation de puissance selon un premier mode de réalisation de l'invention ; La Fig. 7 représente la suite des points de fonctionnement des premier et second terminaux récepteur dans le diagramme de la Fig. 5 lorsque la méthode de la Fig. 6 est mise en oeuvre ; La 10 échanges récepteur oeuvre ; La échanges 15 récepteur Fig. 8A représente schématiquement les de messages entre terminaux émetteur et lorsque la méthode la Fig. 6 est mise en Fig. 8B représente schématiquement les de messages entre terminaux émetteur et lorsqu'une variante de la méthode de la Fig. 6 est mise en oeuvre. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Nous considérerons ici un système de communication 20 sans fil comprenant une pluralité de couples de terminaux émetteur et de terminaux récepteur et supposerons qu'une communication entre les terminaux d'un couple peut interférer avec une communication entre les terminaux d'un autre couple, par exemple par 25 ce que ces deux communications utilisent des ressources de transmission communes (interférence intra-bande). Le système de communication sans fil peut être par exemple un système de communication cellulaire ou un système de communication ad hoc. 30 Dans un but de simplification de la présentation, nous nous placerons tout d'abord dans le cas d'un système ne comportant que deux couples de terminaux, comme illustré en Fig. 1. La Fig. 2 représente un modèle du système de 5 télécommunication de la Fig. 1. Si l'on note x,, x2 les signaux respectivement transmis par les terminaux émetteur 210 et 230 et yl,y2 les signaux respectivement reçus par les terminaux récepteurs 220 et 240, on a : 10 Yi=gilxi+g,2x2+z1 (1) Y2 = gnxl + g22x2 + z2 où z'z2 sont des échantillons de bruit blanc gaussien, g1l,g21 sont les coefficients de canal entre le terminal 15 émetteur 210 et les terminaux récepteurs 220 et 240, respectivement et g22,g,2 sont les coefficients de canal entre le terminal émetteur 230 et les terminaux récepteurs 240 et 220, respectivement. On comprendra que l'interférence entre 20 communications est représentée en (1) par les termes croisés. Pour la première communication, les puissances respectives du signal d'information et de l'interférence due à la seconde communication sont 25 respectivement 2P et 2 g121 P2 . g> > On notera que ce modèle est similaire au modèle classique de canal multi-accès à deux utilisateurs dit MAC (Multiple Access Channel) tel que décrit par 30 exemple dans l'article de Shih-Chun Chang et al. intitulé « Coding for T-User Multiple-Access Channels, » IEEE Transactions on Information Theory, vol. IT-25, no. 6, pp. 684-691, Nov. 1979. Il en diffère cependant fondamentalement au sens où dans le modèle MAC, les terminaux récepteurs reçoivent une information utile et non de l'interférence sur les canaux « croisés ». Le rapport signal à bruit (SNR) au niveau du 10 premier récepteur peut s'exprimer sous la forme : YI= g1 ~ 2 P No (2) De manière similaire, le rapport interférence à 15 bruit au niveau du premier terminal récepteur n'est autre que : g12 2 P2 No (3) 20 De la même façon, les rapports signal à bruit et interférence à bruit au niveau du second terminal récepteur peuvent s'écrire respectivement comme : 2 P2 g22 N O Yz - 25 (4) 82 = 12P, T7(5) On remarquera que les relations suivantes sont et vérifiées = f27282 =1i Y~ (6) (6°) 2 g21 ! _ g12 OÙ f - glll2 et Îz g22 On notera dans la suite p, et p2 les débits respectifs sur la première et la seconde communications 10 et l'on introduit les variables Cl =2n-1, C2=2'2-1 et c12 =2'l"2 -1. Pour un couple de débits (pl,p2) donnés, on peut distinguer plusieurs régimes d'interférence, chaque 15 régime donnant lieu à un traitement distinct. Plus précisément, si l'on considère la première communication, on peut classer l'interférence en trois régimes possibles Schématiquement, dans un premier régime, la 20 puissance de l'interférence due à la seconde communication est plus faible que la puissance du signal d'information reçue au niveau du premier terminal récepteur. Plus précisément, si l'on raisonne en termes de capacité au sens de Shannon, le rapport 25 signal sur bruit plus interférence sur le canal direct entre le terminal émetteur 210 et le terminal récepteur 220 permet de passer le débit p, alors que le rapport signal sur bruit sur le canal « croisé » entre le terminal émetteur 230 et le terminal récepteur 220 ne 30 permet pas de passer le débit p2 (le canal « croisé » 2 2 est en situation de coupure au sens de la théorie de l'information), autrement dit ( p _ p2 >log2(1+INR12)=1og2(1+8~) (8) où SINR' et 1NR12 sont respectivement le rapport signal à bruit plus interférence et le rapport interférence sur bruit au niveau du premier terminal récepteur. On déduit de (7) et (8) les contraintes sur 81 et relatives au premier régime : Cl (1 +81) (9) 8, où. SINR,2 est le rapport « signal sur bruit au niveau du terminal récepteur 220 dans lequel on entend le signal comme le signal d'information de la seconde communication. Dans une seconde étape, une fois la contribution de la seconde communication soustraite du signal reçu, on se retrouve dans le cas d'un signal simplement bruité par un bruit thermique, autrement dit log2 (l+SNRII)= log2 (1+A) (12) où SNRn est le rapport signal sur bruit après élimination de l'interférence due à la seconde 25 communication. Les expressions (11) et (12) se traduisent par les contraintes suivantes sur A et SI : Y, 10 conjointement le signal d'information de la première communication et le signal d' information de la seconde communication au niveau du premier terminal récepteur. Le décodage conjoint des deux signaux d'information pourra, par exemple, être effectué au moyen d'un schéma 15 PIC ou d'une méthode de décodage à maximum de vraisemblance de type MMSE-GDFE (Minimum Mean Square Error - Generalized Decision Feedback Equalizer), de manière connue en soi. Ce régime d'interférence est intermédiaire entre 20 le premier et le troisième dans la mesure où le débit p2 ne vérifie plus (8) et (11), autrement dit : 18 (13) i ô 10g21-1- 1+71/ Cz et 10 y1>C12-81 (18) La Fig. 3 montre un diagramme d'interférence dans lequel on a représenté en abscisses le rapport 15 interférence sur bruit 81 et en ordonnées le rapport puissance sur bruit y1. Ce diagramme est obtenu pour des valeurs données de débit pl et p2, et par conséquent pour des valeurs données de C1, C2, C12 . 20 Le rapport y1 varie de 0 à e = 2 p'nax /No et le rapport 8 varie de 0 a 81'nax 2 Pmax/N0 où P1 PI' et P2' = SI2 2 sont respectivement les puissances d'émission maximales des terminaux 210 et 230. Les droites Al et A2 définies par les équations 25 81 =C2 (cf . (10)) et yl = 1 -1 (cf . (13)) délimitent les z trois régimes d'interférence. Les droites D, , D2, D3 , respectivement définies par les équations y1=C1(1+81) (cf. (9)) ; 71=C12-81 (cf. (18)) ; y1=C1 (cf. (14)) ; donnent la frontière inférieure de puissance, notée A, pour chacun de ces régimes. La zone 310 correspondant au premier régime d'interférence est délimitée par les droites D, et Al ainsi que l'axe des ordonnées, celle correspondant au second régime d'interférence, 320, est délimitée par les droites A,, A2 et D2, et enfin celle, 330, correspondant au troisième régime d'interférence est délimitée par A, et D3. Au-dessous de la frontière inférieure A se trouve une quatrième zone 340 dans laquelle il n'est pas possible de traiter l'interférence pour la qualité de service demandée, ici pour les débits pl et p2. On pourra alors utiliser une autre ressource de transmission, par exemple un autre intervalle de temps de transmission pour éliminer l'interférence entre les deux communications. Si le débit de la première ou de la seconde communication varie, les paramètres des équations des droites D, , D2, D3 et A,, A2 varient également et, par conséquent, les zones correspondant aux différents régimes d'interférence sont modifiées. Pour des débits pl et p2 donnés, on peut déterminer à partir d'une estimation, au niveau du récepteur, de la puissance du signal d'information, de la puissance de l'interférence et de celle du bruit thermique, dans quel régime d'interférence on se trouve, et effectuer le traitement afférent à cette zone. En outre, comme indiqué sur la figure, pour des débits pl et p, donnés, on peut déterminer pour chaque puissance d'interférence, 8;, la plus faible valeur de y, désignée ici par y', permettant de traiter cette interférence. Autrement dit, pour une puissance d'émission du terminal 230 donnée, on peut allouer au terminal émetteur 210 la puissance minimale PI = N2YI ° On remarquera que si l'on se trouve dans la zone 320 ou 330, cette valeur de puissance allouée est bien inférieure à celle, désignée ici par y", qu'il aurait été nécessaire d'allouer si l'on avait assimilé l'interférence à du simple bruit thermique. La Fig. 4 montre un diagramme d'interférence dans lequel on a représenté cette fois-ci en abscisses le rapport puissance sur bruit et en ordonnées le rapport puissance sur bruit A. Le rapport y, varie de 0 à mmax 2 m"/No et le rapport y2 varie de 0 à yrn2 ax = g22 2 p2i.nax /No . Les inégalités (9), (13) et (18) peuvent être réécrites en fonction de rapports y, et y2 (19) 1'2 c2 72 -1 Cl2 -f2Y2 De manière similaire à la Fig. 3, les droites et ,à2 définies par les équations y2=C2/f2 (cf. (10)) et f2 A= y2-1 (cf. (20)), délimitent les trois régimes C2 (20) (21) d'interférence et les droites D1, D2, D3 , respectivement définies par les équations y,=C,Q+f2y2) (cf. (19)) ; y1=C,2-f2y2 (cf. (21)) ; y,=C, (cf. (14)) ; donnent la frontière inférieure de puissance, A. Les zones correspondant aux différents régimes d'interférence sont désignées par 410, 420 et 430. La Fig. 5 représente un diagramme d'interférence identique à celui de la Fig. 4 dans lequel on a indiqué les différents régimes d'interférence pour les deux communications. Le rapport signal à bruit y, est porté en abscisses et le rapport signal à bruit 72 en ordonnées. On obtient neuf zones distinctes, obtenues par intersection des trois zones relatives à la première communication et des trois zones relatives à la seconde communication. Ces neuf zones sont notées S2pq , 1,193, 1 .} (23) y; >(y,+1)cj/fj} (24) où ji est l'indice de la communication interférente. Cette partition s'exprime de manière équivalente 20 au moyen des coordonnées (y,,(5,) (22' ) (23' ) (2 4' ) 25 Dans tous les cas, la partition du diagramme d'interférence pour l'ensemble des deux communications est obtenue au moyen des intersections SZPq = CCDP (1 GJq (25) Chaque zone est associée à un régime d'interférence (indexé par p) de la première communication et à un régime d'interférence (indexé par q) de la seconde communication. La Fig. 6 représente une méthode distribuée d'allocation de puissance selon un premier mode de réalisation de l'invention. Ce mode de réalisation suppose que deux couples de terminaux T =(E1,R1), T2=(E2,R2) aient été préalablement sélectionnés parce que leurs communications interfèrent entre elles. La méthode d'allocation de puissance procède de manière itérative, chaque itération étant composée d'une première demi-itération effectuée par un des deux couples de terminaux et d'une seconde demi-itération effectuée par l'autre couple. Plus précisément, dans une première étape, 610, les terminaux émetteur E, et E2 diffusent les valeurs de débit de consigne p1,p2 ou, de manière équivalente, les valeurs des variables C1 et C2 aux autres terminaux. Les terminaux récepteur R1 et R2 procèdent à une estimation du bruit thermique No et à une estimation de canal. Le terminal R1 transmet au terminal E1 les valeurs No et g,,. De même, le terminal R2 transmet au terminal E2 les valeurs No et g22. Dans un seconde étape, 615, les puissances 5 d'émission des terminaux émetteur E1 et E2 sont initialisées à des valeurs données, Pr et Pr, par exemple aux valeurs maximales de puissance d'émission p" et Pr L' indice n donnant le nombre d' itérations est initialisé à 1. On se place dans la suite à la lierne 10 itération. A l'étape 620, le terminal récepteur R1 estime le rapport signal sur bruit, (n) et le rapport interférence sur bruit e), à l'aide de séquences de symboles pilote transmises par exemple tour à tour par 15 le terminal et le terminal E2. Le terminal R1 classifie l'interférence, c'est-à-dire détermine dans quelle zone wl, se trouve le point de fonctionnement p (ri(n),gl(n) A l'étape 630, le terminal R1 notifie au terminal 20 E1 la valeur du rapport interférence sur bruit el) et le régime d'interférence correspondant, p(''). Cette notification est effectuée à l'aide d'un message spécifique. Selon une variante, seule la valeur du rapport interférence sur bruit 61('') est notifiée au 25 terminal E1, ce dernier se chargeant alors de la classification et de la détermination de l'indice p (n). A l'étape 640, le terminal E1 actualise la puissance d'émission p(n) en fonction de e) et p(") Plus précisément, la puissance e est obtenue au moyen de la fonction donnant la borne inférieure de 5 puissance, définie par intervalles : (1(n) ) pour s(1n) ew1' ( (e» C12 -81(0 pour (5I(n) e wl (51)=C1 pour 81(n) e w3 au moyen de : p(n) = 21 (8 n) 1 No Selon une variante, le calcul de l'étape 640 est effectué par le terminal récepteur R1 et le résultat de (27) est ensuite transmis au terminal El. Dans ce cas, 20 le terminal récepteur R1 n'a pas à transmettre la puissance du bruit thermique No et le coefficient de canal g11 à l'étape 610. Quelle que soit l'option retenue, le terminal E1 émet à la puissance d'émission actualisée Pl('). 25 On comprendra que les étapes 620 à 640 correspondent à la première demi-itération. 10 15 (26) (27) A l'étape 650, le terminal récepteur R, estime le rapport signal sur bruit, y2 , ainsi que le rapport interférence sur bruit P)1 au moyen de séquences de 2 symboles pilote, transmises par exemple tour à tour par 5 le terminal E2 et le terminal El. Le terminal R2 en déduit dans quelle zone co se trouve le point de q(n) fonctionnement (y; el . A l'étape 660, le terminal R2 notifie au terminal E2 la valeur du rapport interférence sur bruit 8' et 10 le régime d'interférence correspondant, n A l'étape 670, le terminal E2 actualise la puissance d' émission e') en fonction de e) et e). Plus précisément, la puissance .F) est obtenue au moyen de la fonction 22 donnant la borne inférieure de 15 puissance, définie par intervalles : (8n» = c2 (1+ pour 82n) E col 20 (e» = c 12 - 82 pour (5n) E GJZ 25 22 (e» = c'2 pour 8(n) E c03 au moyen de : :P2(n) ) N° gn 2 (28) (29) Selon la variante précitée, le calcul de l'étape 670 est effectué par le terminal récepteur R2 et le résultat de {29) est ensuite transmis au terminal E2. Dans ce cas, le terminal récepteur R2 n'a pas à transmettre la puissance du bruit thermique No et le coefficient de canal g22 à l'étape 610. Quelle que soit l'option retenue, le terminal E2 émet à la puissance d' émission actualisée P. On comprendra que les étapes 650 à 670 10 correspondent à la seconde demi-itération. Au terme des deux demi-itérations, on vérifie en 680 si un critère d'arrêt est satisfait. Ce critère peut être par exemple libellé sous la forme d'un nombre maximum d'itérations K. au-delà duquel l'algorithme 15 s'arrête. Alternativement le critère d'arrêt pourra être exprimé sous la forme d'une mesure de convergence. Par exemple, si l'on note e, k=1,2, la puissance d'émission allouée au terminal émetteur Ek à l'itération n, ce critère pourra être exprimé par : p(n+I) p(n) où el) est le vecteur des puissances allouées aux terminaux émetteur à l'itération n , à savoir 25 P(n) = (en) p(n) )7' et e est un seuil prédéterminé. Le critère d'arrêt pourra être encore une combinaison des deux critères précédents, au sens où les itérations 20 5 15 20 25 seront arrêtées dès lors que l'un ou l'autre des critères sera satisfait. Si le critère d'arrêt est satisfait, l'allocation de puissance se termine en 690. A défaut, on incrémente le compteur d'itérations en 685 et on retourne à l'étape 620. La Fig. 7 illustre sur le diagramme d'interférence, les points de fonctionnement des terminaux récepteur R1 et R2 à la suite de chaque demi-itération. Les points sont notés ici (n°,n°°) où n° et n" sont respectivement les nombres de demi-itérations effectuées par 1, et 1' 2. On remarque que ces points sont obtenus par des projections alternées sur les courbes n1 et A2. On peut montrer qu'en règle générale les fonctions et 22 sont k-lipschitzienne avec k La Fig. 8A représente de manière schématique les échanges de messages entre les terminaux émetteur et récepteur des couples T1 = (E1,R1) et T2 = (E2, R2) lorsque la méthode d'allocation de puissance de la Fig. 6 est mise en oeuvre. Lors de l'étape préliminaire 610, terminaux E, et E2 diffusent les valeurs de débit P~ et P2 . 5 Alternativement, ceux-ci diffusent les valeurs des variables correspondantes C1 et C2 . Les terminaux récepteur RI et R2 envoient respectivement aux terminaux E, et E2 les valeurs No , g11 et No , g22. A chaque itération, le terminal R, transmet au 10 terminal E, le niveau d' interférence e et l'indice du régime correspondant, p("). De même le terminal R2 transmet au terminal E2 le niveau d'interférence 5n) et l'indice du régime correspondant, q(n) . 15 La Fig. 8B représente de manière schématique les échanges de messages entre les terminaux émetteur et récepteur des couples T =(E,,Rl) et T2=(E2,R2) lorsque la méthode d'allocation de puissance est mise en oeuvre selon la variante précitée. 20 Selon cette variante, les calculs de Pin) et 4(n) , respectivement aux étapes 640 et 670, sont effectués par les terminaux R, et R2 en lieu et place des terminaux E, et E2. Dans ce cas, les valeurs No , g' et No , g22 ne sont 25 pas transmises lors de l'étape préliminaire 610. A chaque itération, le terminal R1 transmet au terminal E, la valeur de puissance actualisée Pin) à l'étape 640, et le terminal R2 transmet au terminal E, la valeur de puissance actualisée à l'étape 670. La variante de la Fig. 8B sera préférée à celle de la Fig. 8A lorsque les terminaux récepteur sont des stations de base, c'est-à-dire des communications montantes. En revanche, pour des communications descendantes, on pourra opter pour la variante de la Fig. 8A pour soulager la charge de calcul des terminaux mobiles. La méthode d'allocation de puissance selon le premier mode de réalisation de l'invention se généralise à un nombre quelconque K de couples de terminaux et donc à un nombre K de communications correspondantes. Le diagramme d'interférence de la Fig. 4 ou 5 est alors construit dans un espace à K2 dimensions et l'on a PK combinaisons possibles de P régimes d'interférence pour l'ensemble des K communications du système. Les consignes de débit des K communications forment un vecteur de débit p = (PI -1 PK )r Selon un second mode de réalisation, non représenté, on regroupe les couples de terminaux émetteur et récepteur, 1; =(Ek,Rj k=t",K par paires selon un critère prédéterminé. Par exemple, ce critère peut être un ordre de priorité, un degré d' urgence, un facteur de qualité de canal, etc. ou encore une combinaison de tels critères. De préférence, dans une phase préalable, chaque terminal récepteur Rk détermine quelle est la communication k' induisant le plus fort niveau d'interférence. Pour ce faire, on active par exemple tour à tour les différents terminaux émetteur, autres que Ek, en leur allouant la même puissance d'émission P", et l'on compare les niveaux de réception obtenus. Les couples Tk et Tk, sont alors appariés. On notera que cette phase préalable nécessite au plus C. mesures, étant donné que celles-ci peuvent être effectuées en activant tour à tour les terminaux émetteur Ek et en mesurant simultanément les niveaux reçus par les différents terminaux récepteur autres que Rk. Pour une paire (Tk,Tk,) ainsi constituée, on pourra considérer que les autres paires ne sont que faiblement interférentes. L'interférence « résiduelle » due aux autres communications sera traitée comme du bruit thermique. Ce bruit thermique sera pris en compte dans la valeur de No pour le calcul des puissances Pk et Pk 20 à allouer aux terminaux Ek et Ek, . La méthode d'allocation de puissance décrite en relation avec la Fig. 6 peut alors s'appliquer aux paires ainsi constituées. Plus précisément, les étapes 610 à 670 sont effectuées pour chacune de ces paires. 25 On comprendra ainsi qu'une itération de l'algorithme est constituée encore de deux demi-itérations, les premiers couples des différentes paires effectuant le traitement de la première demi-itération et les seconds 10 15 20 25 couples de ces mêmes paires effectuant le traitement de la seconde demi-itération. Le critère d'arrêt est commun à l'ensemble des paires. Il peut être exprimé comme précédemment sous la 5 forme d'un nombre maximal d'itérations et/ou d'une mesure de convergence. Dans ce dernier cas, on utilisera l'expression (30) au vecteur en),(e),,ef des puissances allouées respectivement aux K terminaux d' émission à l'itération n . Selon une variante non représentée, le critère d'arrêt peut être spécifique à chaque paire et être vérifié, de manière indépendante, pour chacune d'entre elles. Ceci permet notamment de ne pas retarder la convergence de la méthode d'allocation de puissance pour les paires bien conditionnées. Avantageusement, on vérifiera également si les puissances allouées remplissent en outre un critère d'homogénéité. Il faut en effet éviter que certains terminaux transmettent à pleine puissance alors que d'autres sont réduits à être quasi-silencieux. Le critère d'homogénéité pourra notamment être libellé en termes de puissances d'émission minimale et maximale ou bien encore en termes de variance des puissances allouées. Si les puissances allouées ne sont pas homogènes, celles qui ne remplissent pas le critère d'homogénéité sont identifiées et les couples de terminaux correspondants sont rejetés. Ces couples pourront par exemple se voir allouer des ressources de transmission différentes. 30 Selon un troisième mode de réalisation, également non représenté, les couples de terminaux Tk=(Ek,Rk), k=1,...,K sont d' abord appariés comme dans le second mode de réalisation. Toutefois, à la différence du second mode de réalisation, on procède à une allocation de puissance paire par paire, chaque paire étant traitée de manière indépendante. Selon une première variante de ce mode de réalisation, on considère que pour toute paire (Tk,Tk,), l'interférence induite par les autres communications v k,k' est négligeable. Afin de s' assurer que cette hypothèse est bien vérifiée, on pourra allouer une puissance nulle, aux autres terminaux Ev, v k,k' pendant que les communications k et k' seront actives. Le procédé de sélection de la paire active (TI,Tj.) peut être répété périodiquement, notamment pour tenir compte de l'évolution des coefficients de canal. Tl peut être également répété pour respecter un principe d'équité, les couples constitutifs de la paire sélectionnée en dernier lieu étant alors exclus du processus de sélection pendant une période donnée. Cette période d'exclusion peut être fonction de la qualité de service requise pour la communication en question, une qualité de service d'autant plus élevée conduisant à une période d'exclusion d'autant plus courte. Selon une seconde variante de ce mode de réalisation, on considère que pour toute paire (Tk,Tk,), l'interférence induite par les communications extérieures à celles de la paire en question peut être assimilée à du bruit thermique. Plus précisément, la puissance de bruit No, transmise le cas échéant aux terminaux émetteur, en 610, et entrant dans le calcul des rapports signal sur bruit et interférence sur bruit, en 620, est augmentée de la puissance de l'interférence induite par les communications vk,k'. Afin de s'assurer que cette hypothèse est bien vérifiée, on pourra allouer une puissance minimale Pm', aux autres terminaux Ev , v k,k' pendant que les communications k et k' seront actives. La sélection de la paire active pourra être effectuée comme dans la première variante. La présente invention a été décrite comme une méthode distribuée d'allocation de puissance respectant une contrainte de débit donnée. L'homme du métier comprendra cependant qu'elle peut s'appliquer à une allocation de puissance respectant une contrainte portant sur un paramètre quelconque de qualité de service, par exemple un temps de latence maximal, un rayon de couverture minimal, un taux d'erreur maximal, etc. sans pour autant sortir du cadre de la présente invention. REVENDICATIONS1. Méthode d'allocation de puissance pour système de télécommunication sans fil comprenant au moins un premier couple de terminaux constitué d'un premier terminal émetteur et d'un premier terminal récepteur et d'un second couple de terminaux constitué d'un second terminal émetteur et d'un second terminal récepteur, une première communication entre les terminaux du premier couple utilisant les mêmes ressources de transmission qu'une seconde communication entre les terminaux du second couple, de sorte que les deux communications interfèrent entre elles, la première communication devant respecter un premier débit de consigne et la seconde communication devant respecter un second débit de consigne, caractérisée en ce que : on initialise (615) les puissances d'émission des premier et second terminaux émetteur à 20 des valeurs prédéterminées ; on effectue une pluralité d'itérations jusqu'à ce qu'un critère d'arrêt soit satisfait (680), chaque itération comprenant une première pluralité d'étapes (620, 630, 640) effectuées par le premier 25 couple de terminaux suivie d'une seconde pluralité d'étapes (650, 660, 670) effectuées par le second couple de terminaux, la première pluralité d'étapes comprenant : (al) une estimation (620) par le premier 30 terminal récepteur du rapport signal sur bruit ainsi que du rapport interférence sur bruit, le signal étantcelui de la première communication et l'interférence étant induite par la seconde communication ; (bl) un calcul (630, 640) de la puissance minimale d'émission du premier terminal émetteur 5 autorisant un premier débit de consigne, à partir des rapports estimés à l'étape précédente ; (cl) une actualisation (640) de la puissance d'émission du premier terminal émetteur à la valeur ainsi calculée ; 10 la seconde pluralité d'étapes comprenant (a2) une estimation (650) par le second terminal récepteur du rapport signal sur bruit ainsi que du rapport interférence sur bruit, le signal étant celui de la seconde communication et l'interférence 15 étant induite par la première communication ; (b2) un calcul (650, 660) de la puissance minimale d'émission du second terminal émetteur autorisant le second débit de consigne, à partir des rapports estimés à l'étape précédente ; 20 (c2) une actualisation (670) de la puissance d'émission du second terminal émetteur à la valeur ainsi calculée. 2. Méthode d'allocation de puissance selon 25 la revendication 1, caractérisée en ce que la puissance minimale d'émission du premier terminal récepteur est obtenue à l'étape (bl) en effectuant une première partition d'un diagramme d'interférence de la première communication en différentes zones, chaque zone 30 correspondant à un type de traitement distinct de l'interférence, en calculant, à partir des premier etsecond débits de consigne, une fonction (ÂI) donnant le rapport signal sur bruit minimal permettant d'effectuer le traitement de l'interférence dans chacune desdites zones, et en déterminant l'image par cette fonction du rapport interférence sur bruit mesurée à l'étape (al). 3. Méthode d'allocation de puissance selon la revendication 2, caractérisée en ce que la puissance minimale d'émission, Pin), du premier terminal émetteur est calculée à la Hème itération par pl(n) ( (5°1(n> No 2 gl I où 2 est ladite première fonction, No est la puissance de bruit thermique et g11 est le coefficient de canal de la première communication, estimés par le premier terminal récepteur et e) est le rapport interférence sur bruit mesurée à l'étape (al) de la nome itération. 4. Méthode d'allocation de puissance selon la revendication 3, caractérisée en ce que le calcul de la puissance minimale d'émission en) est effectué par le premier terminal récepteur. 5. Méthode d'allocation de puissance selon la revendication 3, caractérisée en ce que le calcul de la puissance minimale d'émission J- est effectué par le premier terminal émetteur. 6. Méthode d'allocation de puissance selon l'une des revendications 1 à 5, caractérisée en ce que la puissance minimale d'émission du second terminal récepteur est obtenue à l'étape (b2) en effectuant une seconde partition d'un diagramme d'interférence de la seconde communication en différentes zones, chaque zone correspondant à un type de traitement distinct de l'interférence, en calculant, à partir des premier et second débits de consigne, une seconde fonction (.4) donnant le rapport signal sur bruit minimal permettant d'effectuer le traitement de l'interférence dans chacune des zones de la seconde partition, et en déterminant l'image par cette seconde fonction du rapport interférence sur bruit mesurée à l'étape (a2). 7. Méthode d'allocation de puissance selon la revendication 6, caractérisée en ce que la puissance minimale d'émission, du second terminal émetteur est calculée à la neme itération par : P 2 ( n ) _ 22 (8 ) No l ! g22 où .4 est ladite seconde fonction, No est la puissance de bruit thermique et g22 est le coefficient de canal de la seconde communication, estimés par le second terminal récepteur et ô est le rapport interférence sur bruit mesurée à l'étape (a2) de la nème itération. 8. Méthode d'allocation de puissance selon la revendication 7, caractérisée en ce que le calcul de la puissance minimale d'émission PZ~n~ est effectué par le second terminal récepteur. 9. Méthode d'allocation de puissance selon la revendication 7, caractérisée en ce que le calcul de la puissance minimale d'émission e est effectué 10 par le second terminal émetteur. 10. Méthode d'allocation de puissance pour système de télécommunication sans fil comprenant une pluralité de couples de terminaux, chaque couple étant 15 constitué d'un terminal émetteur et d'un terminal récepteur, caractérisée en ce qu'on alloue initialement la même puissance d'émission (PCOm) aux terminaux émetteurs desdits couples, et que pour une communication, dite première communication entre les 20 terminaux émetteur et récepteur d'un couple, dit premier couple, on détermine la communication, dite seconde communication, entre les terminaux émetteur et récepteur du couple, dit second couple, qui induit le plus fort niveau d'interférence sur la première 25 communication, lesdits premier et second couples formant une paire, et que l'on effectue ensuite l'allocation de puissance aux premier et second terminaux au moyen de la méthode d'allocation selon l'une des revendications précédentes. 30 11. Méthode d'allocation de puissance selon la revendication 10, caractérisée en ce que l'on alloue une puissance d'émission réduite (1") aux terminaux émetteurs autres que ceux appartenant à ladite paire pendant le temps où les première et seconde communications sont actives et que l'on intègre dans l'estimation du bruit thermique l'interférence induite par les communications extérieures à ladite paire. 12. Méthode d'allocation de puissance selon la revendication 10, caractérisée en ce que l'on alloue une puissance d'émission nulle aux terminaux émetteurs autres que ceux appartenant à ladite paire pendant le temps où les première et seconde communications sont actives. 13. Méthode d'allocation de puissance selon l'une des revendications 10 à 12, caractérisée en ce que le premier couple de terminaux est choisi tour à tour parmi les couples de terminaux de ladite pluralité. 14. Méthode d'allocation de puissance selon l'une des revendications 10 à 12, caractérisée en ce que le premier couple de terminaux est choisi parmi les couples de terminaux de ladite pluralité selon un degré d'urgence, un ordre de priorité ou un facteur de qualité de canal de transmission.