FR 2828621 A1 20030214 FR 0110720 A 20010810 Procédé et dispositif de création de matrices de contraintes La présente invention concerne un procédé et dispositif de création de matrices de contraintes en vue de leur utilisation dans la planification des fréquences porteuses à attribuer aux stations de base prévues pour gérer les cellules d'un réseau de radiotéléphonie cellulaire en cours de définition. Un réseau radiocellulaire, par exemple le réseau GSM, est constitué d'une pluralité de stations de base radio terrestres qui sont interconnectées à travers le réseau téléphonique filaire et auxquelles les terminaux mobiles peuvent accéder lorsqu'ils se trouvent dans la cellule radio d'une station. Une liaison dite montante, c'est-à-dire du terminal vers le réseau, s'effectue dans une première bande de fréquence, tandis que la liaison descendante du réseau vers le terminal s'effectue dans une deuxième bande de fréquence. Chaque bande de fréquence est partagée en plusieurs canaux, par exemple 62, chacun de largeur déterminée, par exemple 200 kHz. Un espacement spectral déterminé appelé "écart duplex" sépare la bande de fréquence de la liaison montante de celle de la liaison descendante. Chaque station de base doit avoir une portée suffisante pour que la cellule considérée déborde sur les voisines, afin d'éviter tout risque de coupure de communication lorsqu'un terminal mobile change de cellule. Cela nécessite d'augmenter la puissance d'émission des stations au-delà du strict nécessaire. Un terminal radio mobile reçoit donc, à un niveau utile, les émissions de deux stations lorsqu'il est dans la zone de recouvrement de deux cellules. En pratique, il en reçoit beaucoup plus et, s'il est alors mis sous tension, il se rattache à la station la plus proche, au sens radio, c'est-à-dire lui fournissant le niveau radio le plus élevé parmi les six meilleures. Dans la zone de recouvrement, les émissions de la station non sélectionnée par le terminal constituent un signal parasite d'interférence, du même ordre de grandeur que le signal utile de la cellule de rattachement. En outre, les signaux des deux stations étant constitués de trames de même structure, le terminal ne pourrait pas effectuer un filtrage spécifique de spectre pour atténuer le signal parasite. Pour établir le réseau radio, l'opérateur dispose d'une plage de fréquences qui lui a été allouée et qui constitue une ressource chère. Dans cette plage, il répartit régulièrement des fréquences porteuses, ou canaux de transmission, le pas de fréquence séparant deux porteuses étant suffisant pour que l'interférence intercanaux, compte tenu de la sélectivité des récepteurs, reste inférieure à un seuil de bon fonctionnement déterminé. Le nombre des valeurs de fréquences porteuses disponibles est cependant bien inférieur au nombre de canaux nécessaires dans la totalité du réseau. Il faut donc utiliser plusieurs fois la même valeur de fréquence, mais il faut vérifier que, dans chaque cellule prévue, le niveau de parasites à la fréquence de sa ou ses porteuses, en provenance des stations des autres cellules du voisinage, reste en dessous du seuil de bon fonctionnement. En d'autres termes, en tout point d'une cellule, il faut garantir une marge de sécurité entre le niveau de signal utile reçu et le niveau de signal parasite à la fréquence considérée. Comme un réseau comporte quelques milliers de stations qui se perturbent mutuellement, un calcul de ces perturbations nécessite un grand nombre d'heures de calcul par un calculateur classique, de puissance moyenne. Le problème se complique en outre par le fait qu'un terminal, déjà rattaché à une station, tolère une inversion entre les niveaux relatifs des signaux reçus de sa station et d'une autre station, et ne se rattache à cette dernière que si l'écart des niveaux dépasse une valeur de seuil d'hystérésis. En outre, les niveaux reçus en provenance des diverses stations sont susceptibles d'évoluer aléatoirement dans le temps, en fonction des conditions de transmission radio et de la position du terminal mobile. De plus, les signaux parasites de trames des stations distantes ont des énergies moyennes qui évoluent avec le trafic et, en outre, les phases de leurs porteuses reçues peuvent mutuellement glisser et ainsi produire une modulation aléatoire de l'amplitude instantanée du signal parasite composite reçu, qui en est la somme vectorielle. Lorsque l'opérateur veut planifier le réseau en déterminant les positions des stations et leurs caractéristiques radio, le plan d'attribution des fréquences doit donc être le résultat de calculs d'optimisation visant à attribuer chaque valeur de porteuse un nombre maximal de fois dans le réseau, tout en garantissant au mieux la marge de sécurité sur le signal utile, compte tenu des aléas mentionnés ci-dessus. La marge de sécurité va entre autre être définie par l'utilisation de matrices de contraintes exprimant des contraintes à respecter entre les différentes cellules du réseau. En effet il était déjà connu d'utiliser des matrices de contraintes dites "dures" pour les calculs de planification, mais ces matrices étaient souvent déterminées dès le départ par des critères fixés par le programmeur pour les utiliser ensuite pour la planification. A titre d'exemple le programmeur décidait de sa matrice en affectant un coefficient, par exemple, la valeur zéro à l'intersection de la colonne Cj et de la ligne Ci de la matrice de contrainte, lorsque l'écart entre les deux canaux de deux cellules Ci, Cj était égal à zéro (c'est-à-dire lorsque les deux cellules avaient la même fréquence), la valeur "un" lorsque l'écart entre les fréquences respectives avait une première valeur déterminée et la valeur "deux" lorsque l'écart entre les fréquences respectives avait une deuxième valeur déterminée. Une telle définition était figée une fois pour toute et de plus de façon subjective. La présente invention vise à fournir un outil pour permettre cette optimisation par la création de matrices de contraintes qui soient le plus proche possible de la réalité, en vue de leur utilisation pour la planification d'un réseau de radiotéléphonie cellulaire. A cet effet, l'invention concerne un procédé de création de matrices de contraintes en vue de leur utilisation dans la planification d'un réseau de radiotéléphonie cellulaire, dans lequel, disposant d'un nombre déterminé de valeurs de fréquences porteuses prédéterminées à attribuer à une pluralité de stations radio de base, de positions et puissances d'émission déterminées dans une zone de service, prévues pour gérer une pluralité correspondante de cellules radio qui jouent le rôle de cellule serveuse ou brouilleuse, comportant les étapes suivantes : - D'après les positions de chaque maille d'une cellule serveuse le système de calcul élabore une prédiction représentative de la valeur C du champ reçu dans cette maille de la zone de service à l'une des fréquences prédéterminées de la cellule serveuse, - Dans chacune de ces mailles, le système de calcul calcule la valeur I du champ d'interférence créé par chacune des autres cellules que le système considère comme cellule brouilleuse vis-à-vis de la cellule supposée comme serveuse par le système de calcul, - Dans chacune de ces mailles, le système de calcul calcule la valeur du rapport C/l pour la cellule serveuse où I représente le champ d'interférence créé par chacune des autres cellules brouilleuses de la zone de service et dont la valeur de champ brouilleur I reçu dans la maille est supérieure à une valeur prédéterminée CHAMP~MINI BROUILLEUR mémorisée au départ par le système de calcul, - Par une opération d'intégration sur l'ensemble de la zone de service, le système de calcul détermine le taux de brouillage global tbij de la cellule serveuse par les cellules brouilleuses, ces valeurs peuvent alors être affichées ou mémorisées dans une matrice dite initiale de contraintes en représentant ou mémorisant des niveaux (Ki, Li, Mi) aux intersections des lignes de désignation des cellules serveuses et des colonnes de désignation des cellules brouilleuses. Selon une autre particularité du procédé une première matrice ou tableau est déterminée en supposant chaque cellule serveuse et chaque cellule brouilleuse en cocanal. Selon une autre particularité du procédé le système détermine ainsi à partir de cette matrice de contrainte dite "initiale" des probabilités de service à l'initialisation qui représentent la probabilité pour qu'une cellule donnée soit la cellule serveuse fournissant au mobile le champ le plus élevé. La probabilité de service est déduite de l'écart entre la valeur du champ prédite et la valeur CHAMP~MINI (modifiable) fournie par un fichier Frequences.txt ,cette valeur CHAMP~MINI étant la valeur de champ qui correspond à une probabilité de service de 50% : - le système de calcul sélectionne ainsi les trois cellules constituant les meilleurs serveurs par comparaison des valeurs de la matrice, - puis détermine des probabilités d'assignation par calcul et stocke ces probabilités dans une matrice de contrainte d'assignation. Selon une autre particularité du procédé dans chaque maille, le système de calcul applique à chacune de ces cellules serveuses potentielles un coefficient de pondération qui comporte trois composantes. - probabilité d'assignation - probabilité de service - pondération de sursol, ces trois coefficients sont combinés au niveau de chaque maille pour affecter d'un poids chacune des trois cellules meilleurs serveurs et rendre compte de la probabilité qu'a chacune d'entre elles d'être la cellule serveuse. Selon une autre particularité du procédé le système calcule, à l'aide des trois probabilités des trois meilleurs serveurs, et d'une relation qui les lient aux probabilité pij de transition d'une cellule cell~ivers une autre cellule cell~j en tenant compte de la marge d'hystérésis du "handover" Hom, une seule valeur Proba~init qui permet ainsi de constituer de proche en proche et déterminer une matrice de probabilité initiale. Selon une autre particularité du procédé la probabilité d'assignation rend compte de la répartition probable du temps alloué à chacune des trois cellules serveuses potentielles sur la durée d'une communication, cette probabilité est calculée sur la base de la différence de champ prédit entre les cellules serveuses potentielles et tient compte d'une valeur moyenne de marge de basculement "HO~MARGIN" égale à 5 dB en fait. Selon une autre particularité du procédé les valeurs de probabilité sont mémorisées dans le fichier "Proba~assig.txt". Selon une autre particularité du procédé le coefficient de sursol (ou coefficient morphologique) permet de rendre compte des flux naturels de "handover" qui suivent la morphologie du terrain. Les cellules ayant le même sursol que la cellule serveuse courante sont favorisées. Ce coefficient est fourni par un fichier Poids.txt qui se présente sous la forme d'une matrice faisant intervenir la classe de sursol et le type d'axe de circulation auxquels appartient la maille Selon une autre particularité du procédé le système de calcul effectue les mêmes calculs que dans l'une des revendications 1,3 à 4 pour les cellules utilisant des fréquences à canal adjacent au canal de fréquence de la cellule serveuse. Selon une autre particularité du procédé le système de calcul prend également en compte les constats de brouillage et les voisinages sont pris en compte en ajoutant leur influence (en brouillage cocanal et canal adjacent) aux résultats obtenus à partir de la couverture, cette influence est donnée sous la forme d'un taux fourni par le fichier Fréquences.txt enregistré en mémoire. Selon une autre particularité du procédé le système de calcul disposant de tables de statistiques de fluctuations de conditions de propagation radio affectant les niveaux de réception des émissions des stations, le système détermine des probabilités correspondantes pour qu'un signal reçu présente un niveau le plus élevé de tous sur une durée déterminée. Un autre but de l'invention est de proposer un dispositif de création de matrices de contraintes permettant de mieux tenir compte de la réalité. Ce but est atteint par le fait que le dispositif de création de matrices de contraintes en vue de leur utilisation dans un système permettant la planification d'un réseau de radiotéléphonie cellulaire, dans lequel, un système informatique dispose de moyens d'attribuer un nombre déterminé de valeurs de fréquences porteuses prédéterminées à attribuer à une pluralité de stations radio de base, de positions et puissances d'émission déterminées dans une zone de service et mémorisées dans une base de données, prévues pour gérer une pluralité correspondante de cellules radio qui jouent le rôle de cellule serveuse ou brouilleuse, comportant des moyens d'élaborer, d'après les positions de chaque maille (P) d'une cellule serveuse, une prédiction représentative de la valeur C du champ reçu dans cette maille de la zone de service à l'une des fréquences prédéterminées de la cellule serveuse, - Des moyens de calculer pour chacune de ces mailles, la valeur I du champ d'interférence créé par chacune des autres cellules que le système considère comme cellule brouilleuse vis-à-vis de la cellule supposée comme serveuse par le système de calcul et des moyens de mémorisation de ces valeurs I. - Des moyens de calculer pour chacune de ces mailles, la valeur d'un rapport C/l pour la cellule serveuse où 1 représente le champ d'interférence I créé par chacune des autres cellules brouilleuse de la zone de service et dont la valeur de champ brouilleur 1 reçu dans la maille est supérieure à une valeur prédéterminée CHAMP MINI~BROUILLEUR mémorisée au départ par le système de calcul. - Des moyens d'effectuer une opération d'intégration sur l'ensemble de la zone de service, le système de calcul détermine le taux de brouillage global tbij de la cellule serveuse par les cellules brouilleuses, ces valeurs peuvent alors être affichées ou mémorisées par des moyens dans une matrice dite initiale de contraintes en représentant ou mémorisant des niveaux (Ki, Li, Mi) aux intersections des lignes de désignation des cellules serveuses et des colonnes de désignation des cellules brouilleuses. Selon une autre particularité du dispositif les moyens sont utilisés pour calculer une première matrice ou tableau déterminée en supposant chaque cellule serveuse et chaque cellule brouilleuse en cocanal, puis une deuxième matrice en supposant chaque cellule brouilleuse en canal adjacent de la cellule serveuse. Selon une autre particularité du dispositif il comporte des moyens de déterminer ainsi à partir de cette matrice de contrainte dite "initiale" des probabilités de service à l'initialisation qui représentent la probabilité pour qu'une cellule donnée soit la cellule serveuse fournissant au mobile le champ le plus élevé, la probabilité de service étant déduite de l'écart entre la valeur du champ prédite et la valeur CHAMP MINI (modifiable) fournie par un fichier Frequences.txt ,cette valeur CHAMP~MINI étant la valeur de champ qui correspond à une probabilité de service de 50%, des moyens de sélection du système de calcul sélectionne ainsi les trois cellules constituant les meilleurs serveurs par comparaison des valeurs de la matrice, et des moyens de détermination des probabilités d'assignation par calcul et de stockage de ces probabilités dans une matrice de contrainte d'assignation. Selon une autre particularité du dispositif, dans chaque maille, le système de calcul comporte des moyens de mémorisation de trois coefficients constitués par un coefficient de: - probabilité d'assignation - probabilité de service - pondération de sursol, et des moyens de combiner ces coefficients au niveau de chaque maille pour affecter d'un poids chacune des trois cellules meilleurs serveurs et rendre compte de la probabilité qu'a chacune d'entre elles d'être la cellule serveuse. Selon une autre particularité du dispositif le système calcule, à l'aide des trois probabilités des trois meilleurs serveurs extraites des matrices calculées, et d'une relation qui lie ces probabilités aux probabilité pij de transition d'une cellule cell~i vers une autre cellule cell~jen tenant compte de la marge d'hystérésis du "handover" Hom, une seule valeur Probajnit qui permet ainsi de constituer de proche en proche et déterminer une matrice de probabilité initiale. Selon une autre particularité du dispositif, les valeurs de probabilité d'assignation sont mémorisées dans le fichier "Proba~assig.txt". Selon une autre particularité du dispositif, le coefficient de sursol est fourni par un fichier Poids.txt qui se présente sous la forme d'une matrice faisant intervenir la classe de sursol et le type d'axe de circulation auxquels appartient la maille. Selon une autre particularité du dispositif, le système de calcul comporte des moyens interactifs permettant la prise en compte des constats de brouillage et des voisinages en ajoutant leur influence (en brouillage cocanal et canal adjacent) aux résultats obtenus des calculs précédents, cette influence étant mémorisée sous la forme d'un taux fourni par un fichier Fréquences.txt enregistré en mémoire. Ainsi, l'environnement radio de la cellule est défini avec une meilleure précision, ce qui permet par la suite de mieux pouvoir planifier les fréquences tout en ne nécessitant qu'une puissance de calcul limitée. Ainsi, il est facile de sélectionner un nombre limité de stations, les plus perturbatrices, et on limite donc la puissance de calcul requise, tout en fournissant une bonne précision par le fait que la matrice de contrainte tient mieux compte de la réalité du terrain dans les prévisions de parasites. La valeur du signal composite a ainsi une valeur qui représente bien la diversité des conditions radio aux divers points de la cellule, et fournit donc un modèle fiable pour la comparaison. La modélisation évoquée ci-dessus est ainsi encore améliorée par prise en compte des perturbateurs à fréquences proches, compte tenu de la qualité de réjection de ceux-ci prévue pour les mobiles récepteurs. L'invention sera mieux comprise à l'aide de la description suivante d'un mode préféré mais non limitatif de mise en #uvre du procédé de l'invention, en référence au dessins annexés, dans lesquels : - la figure 1 représente schématiquement une partie d'un réseau de radiotéléphonie cellulaire en cours d'élaboration, - la figure 2 représente une matrice de coefficients de perturbation entre cellules, servant de base de données et obtenue à partir de la mise en oeuvre du procédé de l'invention, - la figure 3. 1 représente respectivement les influences de trois cellules entre elles pour l'élaboration des probabilités de perturbation entre cellules, servant de base de données à la détermination des matrices de contrainte, - la figure 3. 2 représente le diagramme à états pour le calcul des probabilités d'assignation entre les trois cellules, - la figure 3. 3 représente les probabilités d'assignation entre les 3 cellules serveuses, - La figure 4. 1 représente un graphique sur l'évolution des valeurs de probabilité d'assignation en fonction de l'écart de valeur de champ entre CO et C2 (quand C1=C2). - la figure 4. 2 représente l'évolution de la probabilité de service pour des valeurs d'écart données. - la figure 4. 3 représente la matrice faisant intervenir la classe de sursol et le type d'axe de circulation. La figure 1 représente une carte 31 d'une zone géographique d'un réseau de radiotéléphonie cellulaire en cours de définition, sur laquelle ont été reportées les positions prévues de quatre cellules radio, respectivement référencées 1,2, 3 et 4, entourant des stations de base respectives 11,12, 13 et 14. Comme représenté, les cellules 1 à 4 sont proches les unes des autres, des zones de bord de la cellule 1 étant même en recouvrement avec ceux des trois autres cellules 2 à 4. La carte 31 est une représentation de données d'une carte numérique 32 stockée dans une base de données d'un calculateur 30. La carte numérique 32 spécifie les reliefs naturels et artificiels, ainsi que leur nature, telle que forêt, immeubles et autres, ce qui permet de calculer une estimation de l'atténuation radio de la liaison en fonction des reliefs définis par la carte numérique 32. L'opérateur du réseau à implanter dispose d'une plage de fréquences et y a défini une pluralité, par exemple 62, de valeurs de fréquences de porteuses constituant des canaux de transmission, numériques dans cet exemple, tels que ceux du réseau GSM. Les valeurs des fréquences sont également réparties dans la plage considérée selon un pas de fréquence déterminé en fonction de la sélectivité en fréquence des stations de base 11 à 14 et des terminaux mobiles de radiotéléphonie prévus, comme celui référencé 21, dessiné dans la cellule 1. Comme l'illustrent les flèches 41,42, 43 et 44, le terminal 21 peut recevoir les émissions des quatre stations ou cellules 11 à 14, les réceptions respectives représentant un signal utile fourni par la cellule serveuse ou un parasite d'interférence délivré par une des autres cellules dites brouilleuses. La cellule serveuse en un point du réseau est la cellule qui assure la communication. Il s'agit donc de la cellule qui a été sélectionnée par le terminal mobile de radiotéléphonie en veille lors de l'établissement de la communication ou celle qui a été choisie comme cible lors d'une commutation de cellule serveuse dite "handover" Un bilan de liaison radio correspondante peut être estimé d'après les positions relatives des stations 11 à 14 par rapport à celle du terminal mobile 21, c'est-à-dire un calcul de l'atténuation de transmission. Cette atténuation dépend en particulier de la distance entre l'émetteur considéré et le terminal récepteur 21 et peut facilement être calculée pour une propagation en espace libre dans l'air, à partir d'une valeur d'atténuation linéique qui vaut environ 20 dB/km pour la plage de fréquences considérée et pour les 500 premiers mètres, et 30 dB/km au-delà. Le cas échéant, en cas d'obstacles radio, les renseignements sur les reliefs fournis par la base de données 32 permettent d'estimer le bilan de liaison sur une trajectoire de propagation déterminée. Lorsque l'endroit de réception est choisi sensiblement en vue directe des stations émettrices 12 à 14, un calcul d'atténuation de la propagation, pour déterminer le niveau de réception, peut s'effectuer très simplement, en considérant que l'atténuation est celle d'une propagation en espace libre et est représentée par le nombre prédéterminé de décibels par kilomètre (20 ou 30 dB/km). On peut aussi se procurer les valeurs de niveaux de réception par des mesures sur le terrain. Par vue directe, on entend le fait que le premier ellipsoïde de Fresnel n'est pas sensiblement obturé par un obstacle radio. On rappellera que le premier ellipsoïde de Fresnel se définit comme étant le volume engendré par rotation autour de l'axe longitudinal d'une ellipse comportant l'émetteur et le récepteur 21 à ses foyers respectifs, à une distance déterminée, l'excentricité étant définie par le fait que la somme des deux distances de tout point de la surface de l'ellipsoïde aux foyers n'excède pas la distance inter-foyers de plus d'une longueur d'onde de la porteuse. En d'autres termes, la majeure partie de l'énergie reçue à un foyer est transmise à l'intérieur du premier ellipsoïde de Fresnel. De ce fait une absence d'obstacle important dans cet ellipsoïde correspond à une transmission dans l'air, sans atténuation supplémentaire. La carte 31 ainsi renseignée permet donc au calculateur 30 de calculer en tout point un niveau prévisionnel de champ C radioélectrique de réception des signaux des stations 11 à 14, si l'on a préalablement fixé leur niveau d'émission. Ayant choisi a priori les positions des stations 11 à 14 et leurs puissances d'émission, il s'agit alors de déterminer le degré de confinement géographique acceptable pour une même valeur de fréquence porteuse. En effet, comme exposé plus haut, il est souhaitable de réutiliser une même valeur de fréquence porteuse un nombre maximal de fois, sans toutefois que les interférences entre cellules ne dégradent la qualité du service radio au-delà d'un seuil déterminé. Afin de déterminer si ce seuil est atteint dans la cellule 1 par exemple, il est effectué une simulation, en pratique au moyen du calculateur 30, pour estimer le niveau de champ I dû aux parasites sur une fréquence porteuse déterminée. A cet effet, pour planifier le réseau de radiotéléphonie cellulaire, dans lequel on dispose d'un nombre déterminé, donc limité, de valeurs de fréquences porteuses à attribuer à la pluralité de stations radio de base 11 à 14, le procédé de l'invention permet d'évaluer des matrices de contraintes par les étapes suivantes : Les premiers éléments constitutifs de la matrice de contraintes sont élaborés à partir de la prédiction calculée par le système de calcul pour chacune des cellules ayant un impact dans la zone de service. Cette prédiction est décrite, pour une cellule, par la valeur C du champ reçu dans chaque maille p (figure 4.1) de la zone de service. Dans chacune de ces mailles, le système de calcul calcule la valeur du rapport C/l pour la cellule serveuse , par exemple Cell-1, et chacune des autres cellules brouilleuses de la zone de service, par exemple Cell-2, Cell-3 fig. 4. 1, (dont la valeur de champ reçu dans la maille est supérieure à une valeur prédéterminée CHAMP MINI~BROUILLEUR mémorisée au départ par l'utilisateur dans un fichier Frequences. txt du système de calcul 30. La valeur I représente le champ brouileur correspondant à une cellule brouilleuse(par exemple Cell-j) de la cellule serveuse (par exemple Cell-i). De ces valeurs, on peut déduire l'existence d'un brouillage de la cellule serveuse par les autres cellules en cocanal. Ce rapport est également calculé en canal adjacent pour chaque maille P. Par une opération d'intégration sur l'ensemble de la zone de service, il est alors possible de déterminer le taux de brouillage global tbij de la cellule serveuse par les cellules brouilleuses. Ces valeurs peuvent alors être reportées dans la matrice de contraintes. La figure 2 représente une matrice carrée 50i de valeurs d'influences ou perturbations mutuelles entre toutes les cellules 1 à 4 du réseau, appelée matrice de contraintes, ceci pour l'une, fi, des valeurs de fréquences porteuses. La matrice 50i permet ainsi de présenter de façon synthétique, à un logiciel de planification du réseau exploité par le système de calcul 30 ou à un opérateur humain utilisant ce système de calcul, le taux de brouillage tbij en chaque intersection ligne-colonne de la matrice de la cellule serveuse Cj correspondant à la ligne j par la cellule brouilleuse C, correspondant à la colonne i pour l'utiliser dans le choix d'attribution des fréquences aux stations 11à 14. Cette évaluation des taux de brouillage tbij permet d'évaluer deux matrices, une première d'évaluation des brouillages cocanal et une deuxième d'évaluation des brouillages en canal adjacent. Un brouillage cocanal est une interférence due à la réutilisation de la même porteuse dans des cellules voisines. Un brouillage en canal adjacent est une interférence due à l'utilisation de porteuses adjacentes dans des cellules voisines. Le procédé de l'invention permet éventuellement de tenir compte d'un constat de brouillage. Un constat de brouillage en un point entre deux cellules déterminées est un brouillage identifié notamment par des mesures terrain . Le système de calcul prend en compte un constat de brouillage en ajoutant la valeur du champ Taux~Exception d'un fichier FREQUENCE.TXT mémorisé dans le système de calcul à l'élément correspondant dans la matrice de contraintes. Ainsi, les niveaux de perturbation de la cellule 1 par les cellules 2,3, 4 atteignent les valeurs respectives Ki, Li et Mi. De même, la cellule 2 est perturbée à des niveaux Xi et Yi par les cellules respectives 3 et 4, ces deux dernières se perturbant mutuellement à un niveau Zi. Des matrices du même type, comme la matrice 50j pour la fréquence fj de la figure 4, peuvent être établies pour les diverses valeurs de fréquences porteuses, ce qui permet de disposer d'une vue d'ensemble des contraintes d'interférences radio Kj, Lj, Mj, Xj, Yj, Zj. Chaque matrice 50i ou 50j est en fait un bloc de données utilisable par le calculateur 30, la représentation qui en est ici faite étant à but illustratif. Si, dans chaque maille P de la zone de service, la cellule serveuse Celli était systématiquement la cellule "best-server" c'est-à-dire ayant la meilleure probabilité, les calculs précédents seraient immédiats et simples. Toutefois, il ne rendraient pas bien compte de la réalité de fonctionnement du réseau. Les variations du champ reçu et les déplacements locaux du mobile durant une communication entraînent des fluctuations et par voie de conséquence le changement de cellule serveuse (handover). Pour une maille donnée, la cellule serveuse peut, selon les principes de handover de la norme GSM être la cellule serveuse courante ou l'une de ses six meilleures voisines désignées. Pour des raisons pratiques (temps de calcul et ressources informatiques limitées), pour chacune des mailles de la zone de service, le système de calcul ne prend en compte comme cellules serveuses potentielles que les cellules Cell~0, Cell~1, Cell~2 ayant une première valeur de champ (C0),une seconde valeur (C1) et une troisième valeur (C2). L'objectif est d'avoir un calcul de probabilités traduisant au mieux la réalité du système lorsque l'on considère les trois cellules qui peuvent constituer chacune un meilleur serveur (best server). Ce calcul est basé sur l'utilisation du diagramme à états pour les trois cellules représenté à la figure 3. 2, l'état initial correspondant à la mise sous tension du mobile. Pour effectuer ce calcul on utilise les notations suivantes: Proba~Ci(P) : est la probabilité pour que le pixel P de la zone de couverture soit assigné à la cellule cell~i. Proba~Ci~init(P) : est la probabilité (Proba~Ci(P) ) pour que le pixel P de la zone de couverture soit assigné à la cellule cellj,à l'allumage du mobile. HOm : la marge d'hystérésis du "handover" Considérons trois cellules Cell~0, Cell~1 et Cell~2 (figure 3.1). Le champ mesuré en un pixel P de la zone couverte par ces trois cellules est CO(P) issu de Cell~0, C1 (P) issu de Cell-1 et C2 (P) de Cell~2. Un mobile qui se trouve en P capte un champ qui évolue avec une statistique gaussienne de valeur moyenne CO, respectivement C1 et C2 qui varie en log normal. La probabilité pour que ce mobile soit assigné à la cellule i, (ie{0,1,2}), à sa mise sous tension (probabilité d'assignation "à l'init") est la probabilité pour que le champ ci (P) soitle plus grand. La statistique étant la même pour les trois champs en P, cette probabilité ne dépend que de la valeur moyenne Ci des champs reçus et de l'écart type du modèle de propagation qui est pris à 8 dB dans notre cas. Le système de calcul 30 manipule les valeurs moyennes. Ainsi, en chaque pixel P, nous avons pour les trois meilleures serveuses, trois valeurs moyennes CO, C1 et C2 ordonnées. Aussi pour calculer les probabilités d'assignation, le système fait intervenir uniquement l'écart entre les trois valeurs moyennes. Ce qui nous permet toujours de nous ramener à un des trois cas d'écart type suivant: # 0 ; # champ~plus~fort-champ~plus~petit; # champ~intermédiaire-champ~plus~petit). Par exemple, supposons que l'ordonnancement entre les valeurs moyennes soit le suivant : C2#C1#C0, le système calculera pour chaque cellule serveuse Ci la probabilité Proba~Ci~init(C2-C0,C1-C0). Ainsi # Proba~C0~init(C2-C0,C1-C0) traduit le nombre de fois où le champ issu de la gaussienne de valeur moyenne la plus petite des trois a été le plus grand, # Proba~C1~init(C2-C0,C1-C0) traduit le nombre de fois où le champ issu de la gaussienne de valeur moyenne intermédiaire a été le plus grand, # Proba~C2~init(C2-CO,C1-CO) traduit le nombre de fois le champ issu de la gaussienne de valeur moyenne la plus grande des trois a été le plus grand. Les résultats obtenus représentés sur la figure 3. 3 sous l'intitulé Proba~Ci~init montrent pour une probabilité donnée l'assignation d'une cellule selon les variation de C2-CO comprises entre 0 et 30 dB et pour une valeur C1CO maintenue fixe égale à 5 dB. Ainsi pour une probabilité inférieure à 0,2 la cellule serveuse sera celle qui génère le champ CO jusqu'à un écart C2-CO > 12dB. Pour une probabilité comprise entre 0,4 et 0,6 la cellule serveuse sera celle qui génère le champ C1 jusqu'à un écart de 7dB puis deviendra celle qui génère le champ C2. Le procédé de l'invention permet de calculer les probabilités d'assignation d'une cellule en tant que serveur avec prise en compte de la marge d'hystérésis du "handover" Hom. Les probabilités Proba~Ci~init calculées précédemment ne représentent que l'état " initial " pour le calcul des probabilités d'assignation. On utilise le diagramme à états de la figure 3. 2. Les pij sont des probabilités de transition d'une cellule cell~ivers une autre cellule cell~jen tenant compte de la marge d'hystérésis du "handover" HOm. Les probabilités en cours de fonctionnement à un instant n dépendent de la probabilité de l'instant n-1 et de la probabilité de transition selon les équations ci-après: Proba~CO(n) = Proba~CO(n-1)*pOO + Proba~C1(n-1)*p10 + Proba~C2(n-1)*p20 Proba~C1(n) = Proba~C1(n-1)*p11 + Proba~CO(n-1)*p01 + Proba~C2(n-1)*p21 Proba~C2(n) = Proba~C2(n-1)*p22 + Proba~CO(n-1)*p02 + Proba~C1(n- 1)*p12 Avec p01 = proba(C1>CO+HOm && C1>C2) p02 = proba(C2>CO+HOm && C2>C1) pOO = proba(CO>C1-HOm && CO>C2-HOm) = 1-p01-p02 p10 = proba(CO>C1+HOm && CO>C2) p12 = proba(C2>C1+HOm && C2>CO) p11 = proba(C1>CO-HOm && C1>C2-HOm) = 1-p10-p12 p20 = proba(CO>C2+HOm && CO>C1) p21 = proba(C1>C2+HOm && C1>CO) p22 = proba(C2>CO-HOm && C2>C1-HOm) = 1-p21-p20 Les probabilités d'assignation sont obtenues lorsque Proba~Ci(n) Proba~Ci(n-1) c'est-à-dire lorsqu'elles convergent après plusieurs itérations vers une asymptote. Elles sont représentées sur la figure 3. 3 pour C1-CO = 5 dB. Nous pouvons remarquer que la somme des probabilités est effectivement égale à un tout le temps. D'autre part, pour une même valeur moyenne, nous avons la même probabilité. Enfin, nous notons par ailleurs comme prévu l'effet de l'hystérésis qui privilégie la serveuse courante (par rapport aux probabilités à " l'init "). Les probabilités d'assignation calculées précédemment dépendent de la valeur paramétrable de la marge d'hystérésis du basculement de cellule serveuse ("handover") Hom. Il est important de conserver cette possibilité de paramétrage sur le système de calcul La procédure pourrait être la suivante : 1- En entrée, le système de calcul dispose des probabilités à l'initialisation" l'init ". Il s'agit des probabilités qui ont été calculées au début par le système. Elle ne dépendent pas de la marge d'hystérésis du "handover" HOm mais uniquement de la valeur moyenne du champ. Avec les relations qui lient les trois probabilités des trois meilleurs serveurs, une seule valeur permet de déterminer les deux autres. Une matrice Proba~init[36][36] sera fournie, sachant que nous nous sommes limités à un écart maximum de 36 dB entre les C, C1 et C2, au-delà les probabilités deviennent proches de zéro. 2 - A partir de Proba~init[36][36] et de HOm, les probabilités d'assignation seront calculées au début du module CONTRAINTES et stockées. 3 - Ensuite, pour chaque pixel, on ira récupérer la probabilité correspondante en fonction des trois meilleures serveuses. Dans chaque maille, le système de calcul applique à chacune de ces cellules serveuses potentielles un coefficient de pondération qui comporte trois composantes. - probabilité d'assignation - probabilité de service - pondération de sursol La probabilité d'assignation rend compte de la répartition probable du temps alloué à chacune des trois cellules serveuses potentielles sur la durée d'une communication. Cette probabilité est calculée sur la base de la différence de champ prédit entre les cellules serveuses potentielles et tient compte d'une valeur moyenne de marge de basculement HO~MARGIN égale à 5 dB. Les valeurs de probabilité sont obtenues dans le fichier Proba~assig.txt . Le graphique de la figure 4. 1 donne l'évolution des valeurs de probabilité d'assignation en fonction de l'écart de valeur de champ entre CO et C2 (quand C1=C2). La probabilité de service fait prendre en compte la capacité d'une cellule à être en mesure d'assurer effectivement le service (fonction de la valeur du champ reçu de cette cellule dans la maille). La probabilité de service est déduite de l'écart entre la valeur du champ prédite et la valeur CHAMP MINI (modifiable) fournie par le fichier Frequences.txt . CHAMP MINI est la valeur de champ qui correspond à une probabilité de service de 50%. Ce coefficient a pour effet de renforcer le taux de brouillage en fonction du niveau de champ (plus fort est le champ, plus fort est le brouillage). Le graphique représenté en figure 4. 2 montre l'évolution de cette courbe pour des valeurs de l'écart C-CHAMP~MINI. On note que pour des valeurs d'écart inférieures à -12dB, la probabilité est forcée à 0. Enfin, le coefficient de sursol (ou coefficient morphologique) permet de rendre compte des flux naturels de "handover" qui suivent la morphologie du terrain. Les cellules ayant le même sursol que la cellule serveuse courante sont favorisées. Ce coefficient est fourni par un fichier Poids.txt qui se présente sous la forme d'une matrice faisant intervenir la classe de sursol et le type d'axe de circulation auxquels appartient la maille. Cette matrice est reproduite en figure 4.3. Ces trois coefficients sont combinés au niveau de chaque maille pour affecter d'un poids chacune des trois cellules meilleurs serveurs et rendre compte de la probabilité qu'a chacune d'entre elles d'être la cellule serveuse. Les constats de brouillage et les voisinages sont pris en compte en ajoutant leur influence (en brouillage cocanal et canal adjacent) aux résultats obtenus à partir de la couverture. Cette influence est donnée sous la forme d'un taux fourni par le fichier Fréquences.txt . Ces valeurs sont les suivantes : - TAUX~EXCEPTION pour les brouillages - TAUX~ADJ~VOIS~REEL pour les voisinages en canal adjacent et TAUX~COC~VOIS~REEL pour les voisinages en cocanal. Pour compléter cette première évaluation ou la rectifier en fonction de la réalité du réseau, le système de calcul 30 permet d'incorporer dans la matrice de contraintes les informations suivantes : - constats de brouillage : qui informent sur la réalité terrain et permettent de pallier les brouillages non prédits ou sous-estimés par le calcul de couverture, - voisinages : qui impliquent forcément des zones de recouvrement ; voisines sont celles que détermine le système de calcul, mais également les voisines déclarées dans le réseau déclarées par l'utilisateur, - zones neutres : qui permettent de ne pas prendre en compte inutilement des brouillages dans certaines zones (internes ou externes à la zone de service, comme par exemple les forêts). Les fonctions de le système de calcul 30 qui permettent de préparer, de calculer puis de vérifier la matrice de contraintes font l'objet de cette fiche. 1) Prise en compte des constats de brouillage 2) Prise en compte des voisines 3) Prise en compte des zones neutres 4) Calcul de la matrice des contraintes La prise en compte des constats de brouillage a pour objectif de savoir pondérer le calcul théorique par l'ajout de constats de brouillage mesurés sur le réseau (terrain, statistiques...). II s'agit ici de renseigner le système de calcul, avec les brouillages constatés sur le terrain avant même leur calcul théorique. Cela évite une réaffectation possible de canaux sur des cellules dont on connaît déjà l'incompatibilité. Le système de calcul 30 détermine les brouillages à partir des prédictions de couverture. Celles-ci peuvent être entachées d'erreurs liées à l'imprécision des bases de données terrain et à la méthode de calcul. Certains brouillages ne sont donc pas détectés par le système de calcul 30. Les mesures terrain et l'expérience de l'ingénieur radio permettent de corriger ces erreurs de prédiction. Les informations correspondantes, appelées constats de brouillage, doivent être fournies à l'outil, afin qu'il en tienne compte dans son calcul des contraintes de brouillage. Il existe deux types de brouillage : - le brouillage cocanal - le brouillage canal adjacent L'outil permet dans un même fichier (extension .CBR) de déclarer ces deux types de brouillage. Cette déclaration se fait au niveau de l'interface du système de calcul Ce fichier doit être mis à jour par l'ajout ou la suppression des constats de brouillage avant le calcul de la matrice des contraintes et des plans de fréquences successifs. La visualisation par matrice thématique doit être effectuée pour chacune des cellules de la zone de service. Elle permet de visualiser par exemple les contraintes entre une cellule brouilleuse et des cellules brouillées pour des interférences cocanales ou adjacentes. Cette visualisation est un bon moyen de vérifier les contraintes. Par exemple, deux cellules peuvent être en regard et ne pas avoir de contraintes calculées alors que l'expérience terrain ou les mesures démontrent le contraire. Dans ce cas, il faudra déclarer cette nouvelle contrainte en l'ajoutant dans la liste existante. Un constat de brouillage en un point entre deux cellules déterminées est un brouillage identifié notamment par des mesures terrain . Le système de calcul 30 prend en compte un constat de brouillage en ajoutant la valeur du champ Taux Exception du fichier FREQUENCE.TXT à l'élément correspondant dans la matrice de contraintes. Pour cela l'utilisateur du système de calcul utilise une interface logicielle lui permettant de Cliquer dans le menu Fréquences , la commande Planification puis la sous commande Constats de brouillage Une boite de dialogue Constats de Brouillage apparaît à l'écran faisant apparaître une liste "Fréquences / Planification / Constats de Brouillage". L'utilisateur à l'aide de moyens de pointage interactif clique sur le bouton liste . La boite de dialogue Constats apparaît à l'écran. Si un fichier de constat existe déjà, sélectionner le et cliquer sur le bouton Valider , sinon taper un nom de fichier dans Votre Sélection et appuyer sur la touche Entrée Dans ce dernier cas, le système de calcul 30 demande de confirmer la création de ce fichier par l'interrogation "Création ?". De retour dans la boîte de dialogue Constats de brouillage (soit après avoir Validé un fichier existant, soit après avoir Créé un nouveau fichier), cliquer sur le bouton Ajouter . Une boite de dialogue Ajouter Un Constat apparaît à l'écran. Dans le champ N de site (zone Utile Public ), taper le numéro de site dont une cellule est brouillée et valider par la touche Entrée. La liste des cellules du site s'affiche (ainsi que le nom du site). L'utilisateur sélectionne dans cette liste la cellule brouillée, puis dans le champ N de site (zone Brouilleur Public ),il tape le numéro de site dont une cellule est brouilleuse et valide par la touche Entrée La liste des cellules du site s'affiche (ainsi que le nom du site). L'utilisateur sélectionne dans cette liste la cellule brouilleuse puis la Nature du Brouillage , soit Cocanal , soit Adjacent . Après avoir saisi le Nom de Responsable et Validé , l'utilisateur de retour dans la boite de dialogue Constats de Brouillage , répète (en cliquant de nouveau sur le bouton Ajouter ) cette opération autant de fois que nécessaire. Le système de calcul 30 permet également la prise en compte des cellules voisines ou encore des voisinages. Deux cellules sont déclarées voisines afin de permettre des "handovers"; elles ont nécessairement un recouvrement de leur couverture. L'objectif de cette prise en compte est de savoir pondérer le calcul théorique par l'ajout de données de voisinage(fichier de voisines permettant une diminution des HO Losses ). Il s'agit ici de renseigner l'outil sur les voisines déclarées (au sens réseau et pas forcément au sens géographique). Cette déclaration de cellules voisines entre elles, permet d'éviter une affectation de fréquences cocanales ou adjacentes sur ces mêmes cellules, évitant ainsi les risques de brouillage sur les zones de recouvrement de la couverture de ces cellules. L'enjeu est évidemment de limiter autant que possible ce type d'interférences, cause de nombreux HO Losses (indicateur comptabilisant soit les handovers perdus, soit les appels coupés) pour améliorer la qualité des résultats du plan de fréquences. Pour la déclaration des voisines, le système de calcul 30 propose trois méthodes distinctes : 1. création d'un fichier personnel, 2. importation de données, 3. calcul par le menu Voisinage du système de calcul 30. Le système de calcul prend également en compte des zones neutres. L'objectif est de savoir identifier les zones sans abonné afin de ne pas pénaliser les zones de trafic le cas échéant. Afin de faciliter et d'améliorer encore le résultat du calcul des fréquences, il est utile de définir les zones neutres dont les cellules n'interviennent pas ou peu dans le trafic écoulé. Effectivement, il est préférable d'ignorer ces zones sans abonné, en acceptant les interférences éventuelles, afin de privilégier la zone à fort trafic et libérer ainsi des fréquences. Ces zones neutres peuvent par exemple permettre de ne pas polluer le calcul d'un plan de fréquences par la protection inutile d'une montagne surplombant une ville. Autres exemples de zones neutres : - des lacs (Etang de Berre) - la mer - des montagnes - des forêts - des zones sans trafic - les pays étrangers Un autre exemple est celui des frontières. Effectivement, en dehors du territoire couvert il est inutile de s'inquiéter d'éventuelles interférences. Les zones neutres ne sont donc pas forcément contenues dans la zone de service. La modélisation ci-dessus peut aussi prendre en compte des aléas de transmission radio. En pareil cas, disposant de tables de statistiques de fluctuations des conditions de propagation radio affectant les niveaux de réception des émissions des stations voisines 12 à 14, par exemple les conditions météo, le calculateur 30 détermine des probabilités correspondantes pour qu'une inversion de niveaux relatifs se produise entre deux des niveaux de signaux sélectionnés et modifie l'ordre de classement, et la composition est effectuée, dans le calculateur 30, par pondération des signaux reçus, en fonction des probabilités respectives pour qu'un signal reçu présente un niveau le plus élevé de tous sur une durée déterminée. Il doit être évident pour les personnes versées dans l'art que la présente invention permet des modes de réalisation sous de nombreuses autres formes spécifiques sans l'éloigner du domaine d'application de l'invention comme revendiqué. Par conséquent, les présents modes de réalisation doivent être considérés à titre d'illustration, mais peuvent être modifiés dans le domaine défini par la portée des revendications jointes, et l'invention ne doit pas être limitée aux détails donnés ci-dessus. Revendications 1. Procédé de création de matrices de contraintes en vue de leur utilisation dans la planification d'un réseau de radiotéléphonie cellulaire, dans lequel, disposant d'un nombre déterminé de valeurs de fréquences porteuses prédéterminées (fi, fj) à attribuer à une pluralité de stations radio de base (11 - 14), de positions et puissances d'émission déterminées dans une zone de service, prévues pour gérer une pluralité correspondante de cellules radio (1 - 4) qui jouent le rôle de cellule serveuse ou brouilleuse, comportant les étapes suivantes : - D'après les positions de chaque maille (P) d'une cellule serveuse le système de calcul élabore une prédiction représentative de la valeur C du champ reçu dans cette maille de la zone de service à l'une (fi) des fréquences prédéterminées de la cellule serveuse, - Dans chacune de ces mailles, le système de calcul calcule la valeur 1 du champ d'interférence créé par chacune des autres cellules que le système considère comme cellule brouilleuse vis-à-vis de la cellule supposée comme serveuse par le système de calcul, - Dans chacune de ces mailles, le système de calcul calcule la valeur du rapport C/l pour la cellule serveuse où I représente le champ d'interférence créé par chacune des autres cellules brouilleuse de la zone de service et dont la valeur de champ brouilleur 1 reçu dans la maille est supérieure à une valeur prédéterminée CHAMP MINI~BROUILLEUR mémorisée au départ par le système de calcul, - Par une opération d'intégration sur l'ensemble de la zone de service, le système de calcul détermine le taux de brouillage global tb,j de la cellule serveuse par les cellules brouilleuses, ces valeurs peuvent alors être affichées ou mémorisées dans une matrice dite initiale de contraintes. en représentant ou mémorisant des niveaux (Ki, Li, Mi) aux intersections des lignes de désignation des cellules serveuses et des colonnes de désignation des cellules brouilleuses. 2. Procédé selon la revendication 1, dans lequel une première matrice ou tableau est déterminée en supposant chaque cellule serveuse et chaque cellules brouilleuses en cocanal. 3. Procédé selon l'une des revendications 1 et 2, dans lequel le système détermine ainsi à partir de cette matrice de contrainte dite "initiale" des probabilités de service à l'initialisation qui représentent la probabilité pour qu'une cellule donnée soit la cellule serveuse fournissant au mobile le champ le plus élevé, La probabilité de service est déduite de l'écart entre la valeur du champ prédite et la valeur CHAMP~MINI (modifiable) fournie par un fichier Frequences.txt ,cette valeur CHAMP MINI étant la valeur de champ qui correspond à une probabilité de service de 50%, - le système de calcul sélectionne ainsi les trois cellules constituant les meilleurs serveurs par comparaison des valeurs de la matrice, - puis détermine des probabilités d'assignation par calcul et stocke ces probabilités dans une matrice de contrainte d'assignation. 4. Procédé selon l'une des revendications 1 à 3, dans lequel dans chaque maille, le système de calcul applique à chacune de ces cellules serveuses potentielles un coefficient de pondération qui comporte trois composantes. - probabilité d'assignation - probabilité de service - pondération de sursol, ces trois coefficients sont combinés au niveau de chaque maille pour affecter d'un poids chacune des trois cellules meilleurs serveurs et rendre compte de la probabilité qu'a chacune d'entre elles d'être la cellule serveuse. 5. Procédé selon l'une des revendications 1 à 4, dans lequel le système calcule, à l'aide des trois probabilités des trois meilleurs serveurs, et d'une relation qui les lient aux probabilité pij de transition d'une cellule cell~ivers une autre cellule cell~j en tenant compte de la marge d'hystérésis du "handover" Hom, une seule valeur Probajnit qui permet ainsi de constituer de proche en proche de déterminer une matrice de probabilité initiale. 6. Procédé selon l'une des revendications 1 à 5, dans lequel la probabilité d'assignation rend compte de la répartition probable du temps alloué à chacune des trois cellules serveuses potentielles sur la durée d'une communication, cette probabilité est calculée sur la base de la différence de champ prédit entre les cellules serveuses potentielles et tient compte d'une valeur moyenne de marge de basculement "HO~MARGIN" égale à 5 dB en fait. 7. Procédé selon l'une des revendications 1 à 6, dans lequel Les valeurs de probabilité sont mémorisées dans le fichier "Proba~assig.txt". 8. Procédé selon la revendications 1 à 7, dans lequel, le coefficient de sursol (ou coefficient morphologique) permet de rendre compte des flux naturels de "handover" qui suivent la morphologie du terrain. Les cellules ayant le même sursol que la cellule serveuse courante sont favorisées. Ce coefficient est fourni par un fichier Poids. txt qui se présente sous la forme d'une matrice faisant intervenir la classe de sursol et le type d'axe de circulation auxquels appartient la maille. 9. Procédé dans lequel, le système de calcul effectue les mêmes calcul que dans l'une des revendications 1,3 à 4 pour les cellules utilisant des fréquences à canal adjacent au canal de fréquence de la cellule serveuse. 10. Procédé selon l'une des revendications 1 à 9, dans lequel, le système de calcul prend également en compte les constats de brouillage et les voisinages sont pris en compte en ajoutant leur influence (en brouillage cocanal et canal adjacent) aux résultats obtenus à partir de la couverture, cette influence est donnée sous la forme d'un taux fourni par le fichier Fréquences.txt enregistré en mémoire. 11. Procédé selon l'une des revendications 1 à 10, dans lequel, le système de calcul disposant de tables de statistiques de fluctuations de conditions de propagation radio affectant les niveaux de réception des émissions des stations (12 - 14), le système détermine des probabilités correspondantes pour qu'un signal reçu présente un niveau le plus élevé de tous sur une durée déterminée. 12. Dispositif de création de matrices de contraintes en vue de leur utilisation dans un système permettant la planification d'un réseau de radiotéléphonie cellulaire, dans lequel, un système informatique dispose de moyens d'attribuer un nombre déterminé de valeurs de fréquences porteuses prédéterminées (fi, fj) à attribuer à une pluralité de stations radio de base (11 - 14), de positions et puissances d'émission déterminées dans une zone de service et mémorisées dans une base de donnée, prévues pour gérer une pluralité correspondante de cellules radio (1 - 4) qui jouent le rôle de cellule serveuse ou brouilleuse, comportant des moyens d'élaborer, d'après les positions de chaque maille (P) d'une cellule serveuse, une prédiction représentative de la valeur C du champ reçu dans cette maille de la zone de service à l'une (fi) des fréquences prédéterminées de la cellule serveuse, - Des moyens de calculer pour chacune de ces mailles, la valeur I du champ d'interférence créé par chacune des autres cellules que le système considère comme cellule brouilleuse vis-à-vis de la cellule supposée comme serveuse par le système de calcul et des moyens de mémorisation de ces valeurs I. - Des moyens de calculer pour chacune de ces mailles, la valeur d'un rapport C/l pour la cellule serveuse où 1 représente le champ d'interférence créé par chacune des autres cellules brouilleuse de la zone de service et dont la valeur de champ brouilleur I reçu dans la maille est supérieure à une valeur prédéterminée CHAMP MINI BROUILLEUR mémorisée au départ par le système de calcul, - Des moyens d'effectuer une opération d'intégration sur l'ensemble de la zone de service, le système de calcul détermine le taux de brouillage global tbij de la cellule serveuse par les cellules brouilleuses, Ces valeurs peuvent alors être affichées ou mémorisées par des moyens dans une matrice dite initiale de contraintes en représentant ou mémorisant des niveaux (Ki, Li, Mi) aux intersections des lignes de désignation des cellules serveuses et des colonnes de désignation des cellules brouilleuses 13. Dispositif selon la revendication 12, caractérisé en ce que les moyens sont utilisés pour calculer une première matrice ou tableau déterminée en supposant chaque cellule serveuse et chaque cellule brouilleuse en cocanal, puis une deuxième matrice en supposant chaque cellule brouilleuse en canal adjacent de la cellule serveuse. 14. Dispositif selon la revendication 13 caractérisé en ce qu'il comporte des moyens de déterminer ainsi à partir de cette matrice de contrainte dite "initiale" des probabilités de service à l'initialisation qui représentent la probabilité pour qu'une cellule donnée soit la cellule serveuse fournissant au mobile le champ le plus élevé, la probabilité de service étant déduite de l'écart entre la valeur du champ prédite et la valeur CHAMP~MINI (modifiable) fournie par un fichier Frequences.txt ,cette valeur CHAMP~MINI étant la valeur de champ qui correspond à une probabilité de service de 50%, des moyens de sélection du système de calcul sélectionne ainsi les trois cellules constituant les meilleurs serveurs par comparaison des valeurs de la matrice et des moyens de détermination des probabilités d'assignation par calcul et de stockage de ces probabilités dans une matrice de contrainte d'assignation 15. Dispositif selon l'une des revendications 12 à 14, dans lequel dans chaque maille, le système de calcul comporte des moyens de mémorisation de trois coefficients constitués par un coefficient de: - probabilité d'assignation - probabilité de service - pondération de sursol, et des moyens de combiner ces coefficients au niveau de chaque maille pour affecter d'un poids chacune des trois cellules meilleurs serveurs et rendre compte de la probabilité qu'a chacune d'entre elles d'être la cellule serveuse. 16. Dispositif selon l'une des revendications 12 à 14, dans lequel le système calcule, à l'aide des trois probabilités des trois meilleurs serveurs extraites des matrices calculées, et d'une relation qui lie ces probabilités aux probabilité pij de transition d'une cellule cell~ivers une autre cellule cellj en tenant compte de la marge d'hystérésis du "handover" Hom, une seule valeur Proba~init qui permet ainsi de constituer de proche en proche et déterminer une matrice de probabilité initiale. 17. Dispositif selon l'une des revendications 12 à 16, dans lequel les valeurs de probabilité d'assignation sont mémorisées dans le fichier "Proba~assig.txt". 18. Dispositif selon la revendications 12 à 17, dans lequel, le coefficient de sursol est fourni par un fichier Poids.txt qui se présente sous la forme d'une matrice faisant intervenir la classe de sursol et le type d'axe de circulation auxquels appartient la maille. 19. Dispositif selon l'une des revendications 12 à 17, dans lequel, le système de calcul comporte des moyens interactifs permettant la prise en compte des constats de brouillage et des voisinages en ajoutant leur influence (en brouillage cocanal et canal adjacent) aux résultats obtenus des calculs précédents, cette influence étant mémorisée sous la forme d'un taux fourni par un fichier Fréquences.txt enregistré en mémoire.