FR 3099672 A1 20210205 FR 1908697 A 20190731 Description Titre de l'invention : Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire [0001] L'invention concerne un procédé de détermination, mis en oeuvre par ordinateur, d'une puissance maximale d'émission d'un satellite d'une constellation de satellites non-géostationnaires pour respecter un niveau de puissance (l'interférence réglementaire vis-à-vis de satellites géostationnaires [0002] Un angle topocentrique est un angle formé par les directions respectives entre un point de la Terre et un satellite non-géostationnaire, et le point sur Terre et un satellite géostationnaire, autrement dit un angle formé à partir d'un point sur Terre entre un satellite non-géostationnaire et un satellite géostationnaire. [0003] L'angle a est l'angle topocentrique minimal pour un satellite non-géostationnaire et un point sur Terre donné. Autrement dit, l'angle a est obtenu pour le satellite géostationnaire qui se trouve à la position sur l'arc géostationnaire qui minimise l'angle topocentrique. [0004] La puissance maximale d'émission est déterminée selon la valeur de l'angle a pour une station terrestre géostationnaire en un point de la Terre et pour un satellite non-géostationnaire de la constellation émettant vers cette station, de manière à respecter un niveau de puissance d'interférence réglementaire vis-à-vis des satellites géostationnaires. [0005] Afin de protéger les émissions et réceptions de signaux des satellites géosta- tionnaires, il est nécessaire de vérifier que les satellites d'une constellation non-géostationnaire n'interférent pas avec les satellites géostationnaires et leurs stations terrestres. Pour cela, la réglementation internationale impose, dans certaines bandes de fréquence, des contraintes de niveaux de brouillage émis par les satellites des constellations non-géostationnaires pouvant se traduire par des contraintes sur la puissance d'émission de ces constellations non-géostationnaires. [0006] Par ailleurs, plus un brouillage est important, moins il doit être fréquent, la régle- mentation internationale exprime donc ces contraintes sous forme de limites statistiques de puissance de brouillage. Or, limiter le niveau de brouillage revient à limiter la puissance d'émission des satellites de la constellation non-géostationnaire et donc revient à contraindre la capacité de la constellation. [0007] Afin de limiter la contrainte sur la capacité de la constellation, il est donc important de déterminer le plus précisément possible les contraintes opérationnelles de puissance pour minimiser l'impact sur la capacité de la constellation tout en assurant la protection des systèmes géostationnaires. 2 [0008] Actuellement, ce problème est formulé de la manière suivante : étant donnée une puissance rayonnée par un satellite non-géostationnaire appartenant à une constellation non-géostationnaire, quel doit être l'espacement minimal angulaire entre les satellites de la constellation et les satellites géostationnaires pour respecter le niveau de brouillage requis. [0009] [fig.1] décrit ce problème. Une station terrestre GSO _SOL reçoit des signaux d'un satellite géostationnaire GSO. Un satellite non-géostationnaire NGSO, se déplaçant sur une orbite ORB NGSO d'une constellation de satellites non-géostationnaires, peut interférer avec ces signaux et les brouiller. On cherche donc à connaître les dimensions d'une zone ZNA, zone non-adressable, dans laquelle le satellite NSGO peut émettre à une puissance réduite vers la station sol GS0 SOL. Actuellement, la définition de cette zone ZNA est assez conservatrice, et la puissance d'émission, vers la station sol GSO SOL, du satellite NGSO dans la zone ZNA est fortement limitée. [0010] L'invention vise à remédier aux inconvénients précités de l'art antérieur, plus parti- culièrement elle vise à calculer les puissances maximales d'émission de la constellation de satellites non-géostationnaires afin de minimiser l'impact sur sa capacité et d'assurer la protection des systèmes géostationnaires. [0011] Un objet de l'invention est donc un procédé, mis en oeuvre par ordinateur, de déter- mination d'une puissance maximale d'émission d'un satellite non-géostationnaire en direction d'une station terrestre, le procédé comprenant les étapes de : - déterminer la valeur minimale d'un angle topocentrique, formé entre d'une part la direction entre le satellite non-géostationnaire et la station terrestre et d'autre part entre la station terrestre et un point de l'arc géostationnaire correspondant à l'orbite d'une constellation de satellites géostationnaires ; - comparer, en valeur absolue, la valeur minimale de l'angle topocentrique à au moins deux valeurs seuils de sorte que : - si ladite valeur minimale est inférieure, en valeur absolue, au premier seuil, définir la puissance maximale d'émission du satellite non-géostationnaire à une première valeur, - si ladite valeur minimale est comprise, en valeur absolue, entre le premier seuil et le second seuil, définir la puissance maximale d'émission du satellite non-géostationnaire à une deuxième valeur, supérieure à la première valeur, - si ladite valeur minimale est supérieure, en valeur absolue au second seuil, définir la puissance maximale d'émission du satellite non-géostationnaire à une troisième valeur, supérieure à la deuxième valeur ; la première, la deuxième et la troisième valeur de puissance maximale d'émission, le premier et le second seuil étant déterminés de manière à minimiser l'écart entre une distribution des niveaux de puissance reçus par la station terrestre et cumulés sur un in- 3 tervalle de temps et une distribution de référence, la distribution des niveaux de puissance reçus et cumulés étant inférieure à la distribution de référence. [0012] Selon des modes de réalisation : - Les seuils et les valeurs de puissance sont déterminés par les étapes suivantes : 1) Pour une valeur du premier seuil angulaire et pour une valeur du second seuil angulaire : simulations des niveaux de puissance reçus par la station terrestre et cumulés sur un intervalle de temps, comparaison d'une distribution des niveaux de puissance cumulés à la distribution de référence et détermination de la valeur du troisième seuil de puissance de manière à minimiser l'écart entre la distribution de référence et la distribution des niveaux de puissance cumulés ; 2) Pour la même valeur de premier seuil angulaire de l'étape 1) et pour la valeur de troisième seuil de puissance déterminée à l'étape 1) : simulations, pour plusieurs valeurs de second seuil angulaire, des niveaux de puissance reçus par la station et cumulés sur un intervalle de temps, comparaison d'unc distribution des niveaux de puissance cumulés à la distribution de référence pour chaque valeur de second seuil angulaire et détermination de la valeur du second seuil angulaire parmi les valeurs simulées de manière à minimiser l'écart entre la distribution de référence et la distribution des niveaux de puissance cumulés, et de manière à ce que la différence entre le second seuil angulaire et le premier seuil angulaire soit minimale. - Les étapes 1) et 2) sont répétées successivement, la valeur du premier seuil angulaire étant différente entre chaque répétition, de manière à obtenir un triplet de valeurs du premier seuil, du second seuil et de la troisième valeur de puissance à chaque répétition. - Le procédé comprend l'étape suivante exécutée après la répétition des étapes 1) et 2) : 3) Détermination du triplet de valeurs du premier seuil, du second seuil et de la troisième valeur de puissance parmi les triplets obtenus par les répétitions des étapes 1) et 2) de manière à respecter un critère de qualité de service minimale rendue par un système non-géostationnaire auquel appartient le satellite non-géostationnaire - L'étape 3) comprend, pour chaque triplet de valeurs et pour chaque instant d'un intervalle de temps, la détermination des satellites non-géostationnaires ayant une valeur minimale d'angle topocentrique supérieure au premier seuil du triplet, et la simulation des puissances émises par ces satellites déterminés vers la station terrestre et des puissances reçues par la station terrestre afin de déterminer un débit reçu par la station terrestre à chaque instant et pour chaque triplet de manière à sélectionner l'un des triplets de valeurs. - Le triplet de valeurs déterminé est celui dont la somme des débits à chaque instant est la plus élevée. - Le triplet de valeurs déterminé est celui dont le débit est maximal pour une proportion prédéterminée de l'intervalle de temps. - Le procédé comprend les étapes suivantes effectuées après l'étape 3) : 4) Pour chaque instant de l'intervalle de temps, sélectionner des satellites non-géostationnaires ayant une valeur d'angle topocentrique minimale supérieure au premier seuil angulaire ; 5) simulation des niveaux de puissance émis par les satellites sélectionnés à l'étape 4) et reçus par la station terrestre cumulés sur l'intervalle de temps ; 6) identification des instants de l'intervalle de temps et ajustement de la troisième valeur maximale de puissance de manière à minimiser l'écart entre une distribution des niveaux de puissance reçus par la station terrestre et cumulés sur l'intervalle de temps et la distribution de référence, la distribution de référence étant supérieure aux niveaux de puissance cumulés. - Le procédé comprend également une étape de raffinement consistant à déterminer des seuils angulaires et des seuils de puissance supplémentaires de manière à réduire l'écart entre les niveaux de puissance cumulés de l'étape 3) et le niveau de référence, les niveaux de puissance cumulés étant inférieurs au niveau de référence. - Le procédé est mis en oeuvre pour une pluralité d'hypothèses de positions de stations terrestres. - Le procédé comprend également une sélection d'au moins un satellite d'une constellation de satellites non-géostationnaires à chaque instant, les niveaux de puissance reçus par la station résultant de signaux émis par les satellites sélectionnés ; et dans ce cas, le satellite sélectionné peut être celui présentant la puissance reçue par la station terrestre la plus élevée, ou celui situé au nadir de la station terrestre. - Le procédé comprend une étape de correction de la première, de la seconde et de la troisième valeurs maximales de puissance d'émission d'un satellite non-géostationnaire de la constellation, situé en un point diffèrent du nadir de la station terrestre, de manière à ce que la puissance reçue par la station terrestre de ce satellite de la constellation soit égale à la puissance reçue par la station terrestre du satellite sélectionné situé au nadir. [0013] Un autre objet de l'invention est un programme d'ordinateur comportant des ins- tructions pour l'exécution du procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'invention, lorsque le programme est exécuté par un processeur. [0014] Encore un autre objet de l'invention est un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'invention ; lorsque le programme est exécuté par un processeur. [0015] D'autres caractéristiques, détails et avantages de l'invention ressortiront à la lecture de la description faite en référence aux figures annexées données à titre d'exemple et qui représentent, respectivement : [0016] - [fig.1], déjà décrite, un schéma représentant les contraintes actuelles des constellations de satellites non-géostationnaires selon l'art antérieur ; et [0017] - [fig.2], un schéma du principe du procédé selon l'invention ; [0018] - [fig.3], une figure représentant un profil de puissance d'émission obtenu par l' invention ; [0019] - [fig.4], un diagramme d'antenne d'un satellite non-géostationnaire permettant de déterminer le troisième seuil de puissance du profil déterminé par le procédé selon l' invention ; [0020] - [fig.5], une représentation des zones d'émission du point de vue d'un satellite non- géostationnaire, les zones étant déterminées par le procédé selon l'invention ; [0021] - [fig.6], un schéma des étapes du procédé selon un premier mode de réalisation de l' invention ; [0022] - [fig.7a], un schéma détaillant la réalisation des étapes 1) et 2) du procédé selon le même mode de réalisation ; [0023] - [fig.7h], un schéma détaillant la réalisation de l'étape 4) du procédé ; [0024] - [fig.8], un schéma détaillant la réalisation de la troisième étape du procédé selon un mode de réalisation de l'invention ; et [0025] - [fig.9], une figure représentant un profil de niveau de puissance reçu par un point sur Terre en fonction des latitude et longitude de ce point sur Terre obtenu par l'invention. [0026] [fig.2] présente un schéma d'application du procédé selon l'invention. Le procédé a pour objectif de définir des zones (ZO, ZR, ZB) dans lesquelles la puissance d'émission vers la station terrestre GSO_SOL d'un satellite NSGO1 ou NGSO2 est limitée, les satellites NGS01 et NGSO2 appartenant à une constellation de satellites non-géostationnaires ayant une orbite définie ORB_NGSO et émettant vers un point sur Terre GSO_SOL, et également de définir les valeurs de puissance d'émission maximale au sein de ces zones. [0027] La définition de ces zones permet de minimiser les interférences entre des satellites de la constellation NOSO et les systèmes géostationnaires comprenant des satellites géostationnaires 050 de l'arc géostationnaire ARC_GSO et des stations terrestres GSO_SOL, et plus particulièrement de minimiser les interférences entre les satellites de la constellation NGS01 et NGSO2 émettant vers la station terrestre GSO_SOL et les satellites géostationnaires GS0 en liaison avec cette même station terrestre GSO_SOL. 6 [0028] Pour cela, on suppose que la constellation de satellites non-géostationnaires émet vers un point sur Terre sur lequel se trouve une station terrestre géostationnaire GS0 SOL, et on définit l'angle topocentrique minimal a pour cc point sur Terre et les satellites non-géostationnaires de la constellation. L'angle topocentrique est l'angle formé depuis un point sur Terre et entre un satellite non-géostationnaire et un satellite géostationnaire. L'angle a est l'angle topocentrique minimal parmi tous les angles topocentriques défini par rapport à tous les satellites géostationnaires de l'arc géostationnaire ARC G5O. [0029] Par exemple, sur la figure 2, pour le satellite non-géostationnaire NG501 émettant vers la station terrestre GS0 SOL, l'angle topocentrique minimal a est aNGS01, tandis que pour le satellite non-géostationnaire NGSO2 émettant vers la même station terrestre 050 SOL, topocentrique minimal a est aNG502. [0030] Les zones ZR, ZO et ZB définies par le procédé sont des zones angulaires définies en fonction de valeurs seuils angulaires ar et ao. Les valeurs seuils angulaires ar et ao représentent des angles topoccntriques minimaux a. Les seuils angulaires ar (premier seuil angulaire) et ao (second seuil angulaire) permettent de définir la première zone ZR, pour laquelle l'angle topocentrique minimal a, formé par le satellite non-géostationnaire et une station terrestre vis-à-vis de l'arc géostationnaire est inférieur à ar en valeur absolue ; la seconde zone ZO, pour laquelle ce même angle a est compris entre ar et ao en valeur absolue et la troisième zone ZB pour laquelle ce même angle a est supérieur à ao en valeur absolue. [0031] Dans chacune des zones ZR, ZO et ZB, la puissance d'émission des satellites non- géostationnaires, plus précisément la puissance isotrope rayonnée équivalente (PIRE) dans la direction de la station terrestre GSO_SOL, est limitée à une valeur particulière de puissance. Dans la première zone ZR, la puissance d'émission maximale est limitée à PR, dans la seconde zone ZO, la puissance d'émission maximale est limitée à PO supérieure à PR et dans la troisième zone ZB, la puissance d'émission maximale est limitée à Pmax supérieure à PO. [0032] Ainsi, sur la figure 2, les seuils angulaires ar et ao autour de deux satellites non- géostationnaires NGS01 et NGSO2 sont représentés. On compare les valeurs d'angle topocentrique minimales aNGS01 et aNGSO2 à ces seuils angulaires pour déterminer dans quelle zone se trouve le satellite non-géostationnaire et donc à quelle puissance maximale, il peut émettre en direction de la station terrestre GSO_SOL. [0033] Pour le satellite NGSOI, son angle topocentrique minimal aNGS01 est inférieur au premier seuil ar, il se trouve donc en première zone ZR et sa puissance d'émission est limitée à PR. [0034] Pour le satellite NGSO2, son angle topocentrique minimal aNGSO2 est supérieur au second seuil ao, il se trouve donc en troisième zone ZB et sa puissance d'émission est 7 limitée à Pmax. [0035] Un exemple de représentation graphique de la puissance maximale en fonction de l'angle topoccntrique minimal a est donné en figure 3. La figure 3 représente la puissance d'émission maximale d'un satellite non-géostationnaire émettant des signaux de brouillage en direction d'une station terrestre géostationnaire. La puissance maximale est fonction de l'angle a défini pour la station terrestre géostationnaire et le satellite non-géostationnaire de la constellation considéré. Les différents seuils angulaires as, ao associés aux seuils de puissance PR, PO et Pmax sont représentés sur la figure 3. Les seuils angulaires ar et ao permettent de définir les zones ZR, ZO et ZB dans lesquelles la puissance d'émission du satellite non-géostationnaire est limitée par les seuils de puissance PR, PO et Pmax. [0036] Pmax est définie selon les contraintes vis-à-vis du système géostationnaire grâce au procédé selon l'invention. La différence entre Pmax et PO, représentée en figure 3 par X dB, représente la diminution maximale de puissance PIRE que peut supporter le système non-géostationnaire pour garantir un service à un utilisateur sur Terre. La différence entre Pmax et PR, représentée en figure 3 par Y dB, représente la diminution de puissance PIRE du satellite non-géostationnaire liée à la réjection de son antenne. [0037] [fig.4] représente un diagramme d'antenne d'un satellite non-géostationnaire avec en abscisses un angle 0 en degrés et en ordonnées une puissance normalisée rayonnée en décibels dB. Celui-ci permet de déterminer la différence Y dB, donc la première valeur seuil PR. Sur le diagramme d'antenne, on relève la différence entre le lobe principal et la valeur plancher, représentée sur la figure 4 par Y dB. La valeur PR sera alors égale à la différence entre Pmax et Y, Y étant exprimé en décibels. [0038] La définition de ces seuils de puissance maximale PR, PO et Pmax en fonction de seuils angulaires ar et ao permet ainsi de minimiser les interférences d'un satellite non-géostationnaire émettant vers un point sur Terre sur la liaison entre une station terrestre placée en ce point sur Terre et un satellite géostationnaire. [0039] Une fois que les seuils de puissance PR, PO, Pmax et les seuils angulaires ar et ao sont déterminés, les zones déterminées, pour un point donné sur Terre, sont les mêmes pour chaque satellite non-géostationnaire de la constellation de satellites non-géostationnaires émettant vers ce point de la Terre. Pour un point donné sur Terre, on utilisera la même définition de zones (ar et ao) vis-à-vis de ce point pour tous les satellites de la constellation émettant vers ce point. Comme ces zones sont définies pour un seul point sur Terre, il faudra donc déterminer des nouvelles zones pour un autre point sur Terre. [0040] Selon le type de système isoFLUX ou isoPIRE, il faudra adapter les valeurs de puissance Pmax, PO et PR pour chaque satellite de la constellation non- 8 géostationnaire. [0041] Un système isoPIRE est un système dans lequel tous les satellites non- géostationnaires de la constellation, qui émettent en direction de la même station terrestre, sont configurés pour avoir la même puissance d'émission vers cette station terrestre. [0042] Un système isoFLUX est un système dans lequel tous les satellites non- géostationnaire de la constellation, qui émettent en direction de la même station terrestre, sont configurés pour avoir une puissance d'émission telle que la puissance reçue par cette station terrestre est la même pour chacun de ces satellites. [0043] Ainsi, pour un système isoPIRE et pour un point donné sur Terre, on utilisera également les mêmes puissances maximales (Pmax, PO et PR) pour tous les satellites de la constellation émettant vers ce point sur Terre ; tandis que pour un système isoFLUX et pour un point sur Terre, les puissances maximales déterminées par le procédé (Pmax, PO et PR) correspondent aux puissances maximales du satellite de la constellation situé au nadir de ce point sur Terre. Un autre satellite de la constellation, différent du satellite au nadir du point sur Terre, a des valeurs de puissance maximales (Pmax, PO et PR) ajustées de + x dB, avec x correspondant à la différence de perte de propagation entre ce satellite et le satellite au nadir, et la perte de propagation correspondant à l'atténuation de puissance entre la puissance reçue par la station terrestre venant du satellite au nadir et un autre satellite de la constellation. L'atténuation de puissance est due à la différence de trajet entre le satellite au nadir de la station et cet autre satellite de la constellation. [0044] [fig.5] représente les zones ZR, ZO et ZB du point de vue d'un satellite non- géostationnaire par rapport à plusieurs points sur Terre vers lesquels il émet. Pour un satellite non-géostationnaire NGSO fixé, on peut définir ces mêmes zones ZR, ZO et ZB par rapport à ce satellite émettant vers la Terre. [0045] Si le satellite NGSO émet en direction du point NI, l'angle topocentrique minimal aNI est inférieur au premier seuil, donc vers ce point NI, le satellite NGSO est en première zone ZR. [0046] Si le satellite NGSO émet en direction de la station terrestre GSO_SOL, l'angle topo- centrique minimal est quasiment nul (car le satellite NGSO, la station terrestre et un point de l'arc géostationnaire sont alignés), donc il est toujours inférieur au premier seuil ar, le satellite NGSO est toujours en première zone ZR. [0047] Si le satellite NGSO émet en direction du point N2, l'angle topocentrique minimal aN2 est supérieur au premier seuil ar mais est inférieur au second seuil ao, le satellite NGSO est donc en deuxième ZO. [0048] Enfin, si le satellite NGSO émet en direction du point N3, l'angle topocentrique minimal aN3 est supérieur au second seuil ao, le satellite NGSO est donc en troisième 9 zone ZB. [0049] La Terre visible depuis le satellite NGS0 peut être représentée comme une ellipse découpée en bandes où chacune des bandes représente une zone ZR, ZO ou ZB dans laquelle peut émettre le satellite NGS0 avec une valeur de puissance maximale d'émission définie. [0050] [fig.6] présente un schéma des étapes du procédé selon un premier mode de réa- lisation de l'invention. [0051] Dans une première étape 1) 301, on va déterminer la valeur du troisième seuil de puissance Pmax. Pour cela, on fixe arbitrairement une valeur du premier seuil angulaire ar et une valeur du second seuil angulaire no. Avec ces valeurs fixées, on réalise des simulations du niveau de puissance reçu par la station terrestre pour plusieurs instants successifs d'un intervalle de temps venant d'un satellite non-géostationnaire de la constellation. Ensuite, on détermine la distribution des niveaux de puissance reçus par la station et cumulés dans l'intervalle de temps, ce qui donne une courbe représentant la distribution des niveaux de puissance cumulés et reçus par la station avec en abscisse les niveaux de puissance reçus par la station au cours du temps et en ordonnées un pourcentage de temps, le pourcentage 100% représentant l'intervalle de temps complet. On compare alors cette distribution à une distribution de référence, pour déterminer la valeur du troisième seuil de puissance Pmax de manière à minimiser l'écart entre la distribution de référence et la distribution des niveaux de puissance cumulés, la distribution des niveaux de puissance cumulés restant inférieure à la distribution de référence. Un exemple de cette comparaison avec la distribution de référence est détaillé en référence à la figure 7a. [0052] La distribution de référence est définie par une norme comme un seuil acceptable pour ne pas interférer sur la liaison entre la station terrestre et le satellite géostationnaire, l'objectif de l'invention étant de s'approcher au plus près de ce seuil acceptable. [0053] Dans une seconde étape 2) du procédé 302, on va déterminer la valeur du second seuil angulaire no. Pour cela, on fixe la valeur du premier seuil angulaire ar qui sera égal à la même valeur que celle fixée dans l'étape 1) et on fixe la valeur du troisième seuil de puissance Pmax à la valeur déterminée à l'étape 1). Avec ces valeurs, on réalise des simulations, pour plusieurs valeurs de second seuil angulaire no, du niveau de puissance reçu par la station terrestre aux mêmes instants qu'à l'étape 1). On détermine ensuite la distribution des niveaux de puissance reçus par le point sur Terre et cumulés dans l'intervalle de temps pour chaque valeur de second seuil angulaire. On compare ensuite ces distributions de niveaux de puissance cumulés à la distribution de référence, pour déterminer quelle est la valeur de second seuil angulaire no qui permet de minimiser la différence no - ar et qui permet toujours de minimiser l'écart entre la distribution de référence et les distributions des niveaux de puissance cumulés. Lc second seuil angulaire no sera alors égal à cette valeur. [0054] Ces deux étapes sont réalisées pour une seule valeur de premier seuil angulaire ar. Afin de déterminer la valeur du premier seuil angulaire ar, on répète les étapes 1) et 2) en modifiant la valeur de ar à chaque répétition. En répétant ces étapes 1) et 2) pour différentes valeurs de ar, on obtient autant de triplets (ar, ao, Pmax) qu'il y a eu de répétitions des étapes 1) et 2). [0055] Puis, dans une troisième étape 3) du procédé 303, on va chercher à déterminer la valeur du premier seuil angulaire ar, qui sera égale à l'une des valeurs utilisées pour répéter les étapes 1) et 2), autrement dit, on va chercher à déterminer quel triplet (ar, ao, Pmax), parmi les triplets obtenus par les répétitions des étapes 1) et 2), est le plus intéressant pour la constellation. [0056] Pour cela, on peut étudier, par exemple, la qualité de service de la constellation de sa- tellites non-géostationnaires, en particulier le nombre de bits par seconde par hertz que l'on peut envoyer vers des stations terrestres, et la valeur de ar sera déterminée de manière à respecter un critère de qualité de service. Cette étape est détaillée en référence à la figure 8. [0057] [fig.8] représente un schéma détaillant la réalisation de l'étape 3) du procédé selon un mode de réalisation de l'invention. [0058] Dans cet exemple, on détermine le triplet (ar, no, Pmax) et donc la valeur du premier seuil angulaire ar en calculant par simulation la capacité de la constellation de satellites non-géostationnaires vis-à-vis de la station terrestre pour plusieurs instants. [0059] Dans une première étape 3031, on sélectionne des satellites non-géostationnaires pouvant émettre vers la station terrestre à chaque instant d'un intervalle de temps pour chaque triplet de valeurs (ar, ao, Pmax). Les satellites sélectionnés sont ceux capables de rendre un service en direction de la station terrestre, donc les satellites non-géostationnaires présentant un angle alpha a supérieur à la valeur ar du triplet. [0060] A l'issue de la première étape, il y a donc un ensemble de satellites non- géostationnaires sélectionnés pour chaque triplet (ar, ao, Pmax). [0061] Puis dans une seconde étape 3032, on détermine par simulations les niveaux de puissances reçus par la station et les niveaux de puissance émis par les satellites non-géostationnaires sélectionnés pour chaque instant et pour chaque triplet (ar, ao, Pmax), sachant que les satellites sélectionnés auront chacun une puissance d'émission maximale égale à PO ou Pmax selon leur valeur d'angle topocentrique minimal (inférieure ou supérieure à ao). [0062] Puis, dans l'étape suivante 3033, on calcule la capacité de la constellation pour chaque instant et pour chaque triplet (ar, ao, Pmax), autrement dit, on calcule l'efficacité spectrale reçue par la station terrestre pour chaque instant et pour chaque 11 triplet (ar, no, Pmax), l'efficacité spectrale étant exprimée en bits/s/Hz. A la place de l'efficacité spectrale, il est également possible de calculer le débit reçu par la station terrestre qui s'exprime en bits/s. [0063] Avec ccs valeurs de capacité pour chaque instant, on peut déterminer le triplet (nr, no, Pmax) et donc la valeur de ar de plusieurs manières. Les étapes 3034 et 3035 détaillent deux de ces manières. [0064] A la suite de l'étape 3033, on peut additionner toutes les valeurs de capacité obtenues aux différents instants pour chaque triplet (ar, no, Pmax) 3034 et le triplet (ar, no, Pmax) donnant la valeur de ar, qui sera retenue par le procédé, est celui donnant la somme des capacités la plus élevée. [0065] A la suite de l'étape 3033, on peut, à la place de l'étape 3034, calculer les capacités cumulées sur l'intervalle de temps pour chaque triplet (ar, no, Pmax) 3035. Le triplet (ar, no, Pmax) donnant la valeur de ar, qui sera retenue par le procédé, est celui donnant la capacité la plus élevée pour un certain pourcentage de temps, par exemple pour 95% du temps. [0066] Une étape supplémentaire 4) (304) de raffinement des limites de puissance maximale en fonction de l'angle n peut également être présente. Elle consiste à définir des valeurs seuils supplémentaires de puissance et des seuils angulaires supplémentaires de manière à réduire l'écart entre la distribution de référence et la distribution des niveaux de puissance reçus et cumulés par la station terrestre. [0067] L'objectif de cette étape de raffinement est d'ajouter des seuils de puissance et des seuils angulaires entre les seuils déjà déterminés de manière à minimiser les interférences du satellite non-géostationnaire vers un point sur Terre tout en ayant une puissance d'émission forte. En effet, quand on s'approche d'un seuil angulaire, le seuil de puissance associé peut être légèrement dépassé, puisqu'avec un angle égal à ce seuil angulaire, on passe à une limite de puissance égale à un autre seuil plus élevé. On cherche donc dans cette étape à relâcher les contraintes de niveau de puissance au sein d'une zone grâce à la détermination de ces seuils supplémentaires. [0068] Ainsi, dans cette étape, on cherche à déterminer au moins un seuil de puissance et un seuil angulaire intermédiaire avant le premier seuil angulaire et entre le premier et second seuil angulaire. Comme pour les étapes précédentes, pour faire cela, on fait varier une valeur de seuil angulaire nr_i associé à un seuil de puissance PR_1, étant inférieur à ar et PR_i étant supérieur à PR et inférieur au seuil de puissance PO, de manière à s'approcher le plus possible de la norme. On peut également déterminer une valeur de seuil angulaire no_i associé à un seuil de puissance P0_1 où ao_l est supérieur à ar et inférieur à no et P0_1 est supérieur à PO et inférieur à Pmax. [0069] Comme précédemment, cela est réalisé en effectuant des simulations des niveaux de puissance reçus par le point sur Terre à plusieurs instants dans l'intervalle de temps en 12 prenant les valeurs de Pmax, ar et ao déterminées dans les étapes précédentes et en faisant varier les valeurs des seuils supplémentaires. Puis on détermine les niveaux de puissance reçus et cumulés dans l'intervalle de temps pour chacune des valeurs de seuils supplémentaires simulés. On les compare au niveau de référence et on choisit les valeurs de seuils supplémentaires qui permettent dc minimiser la différence entre le niveau de référence et les niveaux cumulés tout en la gardant positive. [0070] Les étapes 1) à 4) sont réalisées pour une station terrestre géostationnaire en un point sur Terre. [0071] Selon un autre mode de réalisation, ces étapes sont réalisées pour une pluralité dc stations terrestres géostationnaires représentant N points sur Terre, N étant un entier. Les étapes 1) à 4) sont donc répétées N fois. A chaque fois que le point sur Terre varie, les zones ZR, ZO et ZB ainsi que leurs seuils angulaires ar et ao et leurs valeurs limite de puissance PR, PO et Pmax sont à nouveau définies par le procédé. En effet, pour un même satellite géostationnaire considéré, l'angle topocentrique minimal a varie quand le point sur Terre varie, comme le représente la figure 5. Pour N points sur Terre considérés, le procédé est donc répété N fois. [0072] Ainsi, le procédé de l'invention peut également être mis en oeuvre de manière à minimiser les interférences entre des satellites d'une constellation de satellites non-géostationnaires et un système géostationnaire comprenant l'ensemble des satellites géostationnaires de l'arc géostationnaire et toutes les stations terrestres géostationnaires. [0073] Selon un autre mode de réalisation, les simulations effectuées dans les étapes 1) et 2) ne déterminent non pas les niveaux de puissance reçus par la station terrestre émis par un seul satellite non-géostationnaire, mais les niveaux de puissance reçus par la station terrestre et émis par plusieurs satellites non-géostationnaires Cela permet de prendre en compte les cas où plusieurs satellites non-géostationnaires émettent en direction d'un même point sur Terre. [0074] Les niveaux de puissance limites PR, PO et Pmax prennent alors en compte toutes les émissions de signaux de la constellation non-géostationnaire en direction de ce point sur Terre et sont les mêmes pour chaque satellite de la constellation émettant vers ce point. [0075] Selon un autre mode de réalisation, l'étape 1) comprend également la sélection du ou des satellites de la constellation non-géostationnaire, utilisés pour la mise en oeuvre du procédé, notamment pour les simulations des niveaux de puissance reçus par la station terrestre, pour chaque instant de l'intervalle de temps. Les niveaux de puissance simulés et reçus par le point sur Terre sont alors issus uniquement de ces satellites sélectionnés. Cela permet de sélectionner les satellites non-géostationnaires les plus pertinents pour effectuer ces simulations et ainsi limiter les temps de calcul. 13 [0076] Les satellites sélectionnés seront par exemple, ceux présentant le niveau de puissance reçu par le point sur Terre le plus élevé. Ils pourront être, selon un autre exemple, les satellites ayant un angle topocentrique minimal inférieur au premier seuil angulaire ar fixé pour l'étape 1) du procédé. [0077] Il est important de conserver le même critère de sélection entre chaque répétition des étapes 1) et 2), c'est-à-dire entre chaque variation de l'angle ar. [0078] Après avoir déterminé les valeurs des seuils angulaires ar et no ainsi que les puissances maximales d'émission associées aux zones pour un point sur Terre vers lequel émet un satellite de la constellation, on peut déterminer la puissance maximale d'émission de ce satellite. Pour cela, on détermine la valeur minimale de l'angle topocentrique formé par cc point sur Terre, le satellite de la constellation émettant vers ce point et un point de l'arc géostationnaire. Puis on compare cette valeur minimale aux seuils angulaires : si la valeur minimale de l'angle topocentrique est inférieure à ar, alors le satellite aura une puissance d'émission limitée à PR, si la valeur minimale de l'angle topocentrique est comprise entre ar et no, alors le satellite aura une puissance d'émission limitée à PO, et si la valeur minimale de l'angle topocentrique est supérieure à Pmac, alors le satellite aura une puissance d'émission limitée à Pmax. [0079] Si plusieurs satellites de la constellation non-géostationnaire émettent vers le même point sur Terre, on applique cela à chaque satellite. Ainsi chaque satellite aura sa propre limite de puissance d'émission selon sa valeur minimale d'angle topocentrique. [0080] Selon un autre mode de réalisation de l'invention, le procédé comprend également des étapes supplémentaires effectuées après l'étape 3) qui permettent d'ajuster les valeurs de puissance maximale d'émission à une situation opérationnelle. [0081] Ces étapes supplémentaires consistent à sélectionner, dans une étape 4), pour chaque instant de l'intervalle de temps, des satellites non-géostationnaires ayant une valeur d'angle topocentrique minimale a supérieure au premier seuil angulaire ar défini à l'étape 3) du procédé vis-à-vis de la station terrestre considérée et de l'arc géostationnaire. Les satellites sélectionnés ont donc une puissance d'émission égale à PO si leur angle topocentrique minimal a est inférieur au second seuil ao ou une puissance d'émission égale à Pmax si leur angle topocentrique minimal a est supérieur au second seuil ao. [0082] Dans une étape suivante 5), on simule les niveaux de puissance émis par ces satellites sélectionnés et reçus par la station terrestre cumulés sur l'intervalle de temps. [0083] Puis dans une étape 6), on identifie des instants de l'intervalle de temps et on ajuste la troisième valeur maximale de puissance Pmax de ces instants de manière à minimiser l'écart entre une distribution des niveaux de puissance reçus par la station terrestre et cumulés sur l'intervalle de temps et la distribution de référence REF, la distribution de référence étant supérieure aux niveaux de puissance cumulés. 14 [0084] [fig.7a] présente un exemple de comparaison de la distribution des niveaux de puissance reçus par la station terrestre et cumulés issus des simulations des étapes 1) et 2) avec la distribution de référence (REF). Les niveaux de puissance reçus par la station sont représentés en abscisses et le pourcentage de temps est représenté en ordonnées. La distribution de référence REF donne les valeurs limite autorisées de puissance reçue par la station depuis le satellite non-géostationnaire en fonction du temps. Par exemple, le satellite ne peut jamais émettre avec un niveau de puissance reçu au sol au-delà de -160 dBW/m2/40kHz, tandis que 50% du temps, il doit présenter un niveau de puissance inférieur à -175 dBW/m2/40kHz. [0085] La distribution des niveaux de puissance cumulés obtenue à l'étape 1) est représentée par la courbe 1). En comparant cette distribution 1) à la référence REF, on peut déterminer la valeur du troisième seuil de puissance Pmax. L'objectif est de translater la courbe issue des simulations vers la courbe de référence tout en restant inférieur à cette référence. La translation effectuée permet de déterminer la valeur de Pmax. [0086] Après avoir déterminé la valeur seuil Pmax, on obtient la distribution 2) de niveaux de puissance cumulés. On remarque que cette distribution est toujours inférieure à la référence REF et qu'il est peut-être possible de s'en approcher un peu plus. Le but de l'étape 2) est donc de faire varier les valeurs du second seuil angulaire ao pour se rapprocher de la distribution de référence REF. La courbe 2) et 2bis) représentent deux distributions de niveaux de puissance reçus par la station terrestre et cumulés sur l'intervalle de temps pour deux valeurs de ao différentes. On remarque que la distribution 2bis) présente des niveaux de puissance supérieurs à la référence REF, donc la valeur de seuil angulaire donnant cette distribution 2bis) est éliminée [0087] [fig.7b] présente un exemple de comparaison de la distribution des niveaux de puissance reçus par la station et cumulés dans le temps avec la distribution de référence REF pour l'étape 4) de raffinement. Les niveaux de puissance reçus par la station terrestre sont représentés en abscisses et le pourcentage de temps est représenté en ordonnées. Comme pour les étapes précédentes, on cherche à se rapprocher de la distribution de référence REF. Lors de l'étape 4), on introduit des valeurs de seuils angulaires et de puissance supplémentaires, qui permettent de relâcher les contraintes sur les niveaux de puissance en bordure des zones ZR et ZO, et ainsi de passer de la distribution 3) à la distribution 4) pour minimiser l'écart entre la distribution de référence REF et la distribution des niveaux de puissance cumulés et reçus par la station. [0088] Si le procédé est réalisé pour une pluralité de stations terrestres, on peut déterminer des seuils de puissance Pmax, PO, PR et des seuils angulaires ar, ao selon la figure 3 pour chacune de ces stations terrestres, et ensuite représenter les puissances maximales d'émission du satellite non-géostationnaire en fonction de n'importe quel point sur Terre. [0089] [fig.9] présente une figure représentant les puissances maximales d'émission d'un satellite non-géostationnaire en fonction de n'importe quel point sur Terre, défini par sa latitude et longitude. Cette figure est obtenue en agrégeant des puissances maximales d'émission en fonction de l'angle a, du même type que celui de la figure 3, déterminées pour plusieurs stations terrestres. La position du point sur Terre (ou de la station terrestre) se trouve en abscisses (longitude et latitude) et la puissance maximale d'émission du satellite non-géostationnaire est donnée en ordonnées. Ainsi, on peut voir sur cette figure que pour un point se trouvant proche de l'Equateur, la constellation ne pourra émettre au maximum qu'avec un niveau de -160 dB, alors qu'en s'approchant des pôles, elle pourra émettre au maximum avec un niveau de puissance de -135 dB. [0090] L'invention peut être mise en oeuvre en tant que programme d'ordinateur comportant des instructions pour son exécution. Le programme d'ordinateur peut être enregistré sur un support d'enregistrement lisible par un processeur. Le support peut être électronique, magnétique, optique ou électromagnétique. [0091] En particulier, l'invention peut être implémentée par un dispositif comprenant un processeur et une mémoire. Le processeur peut être un processeur générique, un processeur spécifique, un circuit intégré propre à une application (connu aussi sous le nom anglais d'ASTC pour « Application-Specific Tntegrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour Field-Programmable Gate Array »). [0092] Le dispositif peut utiliser un ou plusieurs circuits électroniques dédiés ou un circuit à usage général. La technique de l'invention peut se réaliser sur une machine de calcul reprogrammable (un processeur ou un micro-contrôleur par exemple) exécutant un programme comprenant une séquence d'instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre module matériel). [0093] Selon un mode de réalisation, le dispositif comprend au moins un support de stockage lisible par ordinateur (RAM, ROM, EEPROM, mémoire flash ou une autre technologie de mémoire, CD-ROM, DVD ou un autre support à disque optique, cassette magnétique, bande magnétique, disque de stockage non transitoire lisible par ordinateur) codé avec un programme d'ordinateur (c'est-à-dire plusieurs instructions exécutables) qui, lorsqu'il est exécuté sur un processeur ou plusieurs processeurs, effectue les fonctions des modes de réalisation de l'invention décrits précédemment. [0094] A titre d'exemple d'architecture matérielle adaptée à mettre en oeuvre l'invention, un dispositif selon l'invention peut comporter un bus de communication auquel sont reliées une unité centrale de traitement ou microprocesseur (CPU, acronyme de Central Processing Unit » en anglais), une mémoire morte (ROM, acronyme (le 16 « Read Only Memory » en anglais) pouvant comporter les programmes nécessaires à la mise en oeuvre de l'invention ; une mémoire vive ou mémoire cache (RAM, acronyme de « Random Acccss Memory » en anglais) comportant des registres adaptés à enregistrer des variables et paramètres créées et modifiés au cours de l'exécution des programmes précités ; et une interface de communication ou E/S (VO, acronyme de « Input/Output » en anglais) adaptée à transmettre et à recevoir des données. [0095] La référence à un programme d'ordinateur qui, lorsqu'il est exécuté, effectue l'une quelconque des fonctions décrites précédemment, ne se limite pas à un programme d'application s'exécutant sur un ordinateur hôte unique. Au contraire, les termes programme d'ordinateur et logiciel sont utilisés ici dans un sens général pour faire référence à tout type de code informatique (par exemple, un logiciel d'application, un micro logiciel, un microcode, ou toute autre forme d'instruction d'ordinateur) qui peut être utilisé pour programmer un ou plusieurs processeurs pour mettre en oeuvre des aspects des techniques décrites ici. Les moyens ou ressources informatiques peuvent notamment être distribués (« Cloud computing »), éventuellement avec selon des technologies de pair-à-pair. Le code logiciel peut être exécuté sur n'importe quel processeur approprié (par exemple, un microprocesseur) ou coeur de processeur ou un ensemble de processeurs, qu'ils soient prévus dans un dispositif de calcul unique ou répartis entre plusieurs dispositifs de calcul (par exemple tels qu'éventuellement accessibles dans l'environnement du dispositif). Le code exécutable de chaque programme permettant au dispositif programmable de mettre en oeuvre les processus selon l'invention, peut être stocké, par exemple, dans le disque dur ou en mémoire morte. De manière générale, le ou les programmes pourront être chargés dans un des moyens de stockage du dispositif avant d'être exécutés. L'unité centrale peut commander et diriger l'exécution des instructions ou portions de code logiciel du ou des programmes selon l'invention, instructions qui sont stockées dans le disque dur ou dans la mémoire morte ou bien dans les autres éléments de stockage précités. 17 [Revendication I] [Revendication 2] REVENDICATIONSProcédé, mis en oeuvre par ordinateur, de détermination d'une puissance maximale d'émission (Pmax, PR, PO) d'un satellite non-géostationnaire (NSGO , NGSO2) en direction d'une station terrestre (GS0 SOL), le procédé comprenant les étapes de : - déterminer la valeur minimale d'un angle topocentrique (aNGS0 , aNGSO2), formé entre d'une part la direction entre le satellite non-géostationnaire et la station terrestre et d'autre part entre la station terrestre et un point de l'arc géostationnaire (ARC GSO) correspondant à l'orbite d'une constellation de satellites géostationnaires (GSO) ; - comparer, en valeur absolue, la valeur minimale de l'angle topocentrique à au moins deux valeurs seuils (ar, ao) de sorte que : - si ladite valeur minimale est inférieure, en valeur absolue, au premier seuil (ar), définir la puissance maximale d'émission du satellite non-géostationnaire à une première valeur (PR), - si ladite valeur minimale est comprise, en valeur absolue, entre le premier seuil et le second seuil (ao), définir la puissance maximale d'émission du satellite non-géostationnaire à une deuxième valeur (PO), supérieure à la première valeur, - si ladite valeur minimale est supérieure, en valeur absolue au second seuil, définir la puissance maximale d'émission du satellite non-géostationnaire à une troisième valeur (Pmax), supérieure à la deuxième valeur ; la première, la deuxième et la troisième valeur de puissance maximale d'émission (PR, PO, Pmax), le premier et le second seuil (ar, ao) étant déterminés de manière à minimiser l'écart entre une distribution des niveaux de puissance reçus par la station terrestre (GSO_SOL) et cumulés sur un intervalle de temps et une distribution de référence (REF), la distribution des niveaux de puissance reçus et cumulés étant inférieure à la distribution de référence (REF). Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 1 dans lequel les seuils (ar, ao) et les valeurs de puissance (PR, PO, Pmax) sont déterminées par les étapes suivantes 1) Pour une valeur du premier seuil angulaire (ar) et pour une valeur du second seuil angulaire (ao) : simulations des niveaux de puissance reçus par la station terrestre (GSO_SOL) et cumulés sur un intervalle de 18 [Revendication 3] [Revendication 4] [Revendication 5] temps, comparaison d'une distribution des niveaux de puissance cumulés à la distribution de référence (REF) et détermination de la valeur du troisième seuil de puissance (Pmax) de manière à minimiser l'écart entre la distribution de référence (REF) et la distribution des niveaux de puissance cumulés ; 2) Pour la même valeur de premier seuil angulaire (ar) de l'étape 1) et pour la valeur de troisième seuil de puissance (Pmax) déterminée à l'étape 1) : simulations, pour plusieurs valeurs de second seuil angulaire (no), des niveaux de puissance reçus par la station (GS0 SOL) et cumulés sur un intervalle de temps, comparaison d'une distribution des niveaux de puissance cumulés à la distribution de référence (REF) pour chaque valeur de second seuil angulaire et détermination de la valeur du second seuil angulaire (no) parmi les valeurs simulées de manière à minimiser l'écart entre la distribution de référence (REF) et la distribution des niveaux de puissance cumulés, et de manière à cc que la différence entre le second seuil angulaire (no) et le premier seuil angulaire (ar) soit minimale. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 2 dans lequel les étapes 1) et 2) sont répétées successivement, la valeur du premier seuil angulaire (ar) étant différente entre chaque répétition, de manière à obtenir un triplet de valeurs du premier seuil, du second seuil et de la troisième valeur de puissance (ar, ao, Pmax) à chaque répétition. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 3 comprenant l'étape suivante exécutée après la répétition des étapes 1) et 2) : 3) Détermination du triplet de valeurs du premier seuil, du second seuil et de la troisième valeur de puissance (ar, no, Pmax) parmi les triplets obtenus par les répétitions des étapes 1) et 2) de manière à respecter un critère de qualité de service minimale rendue par un système non-géostationnaire auquel appartient le satellite non-géostationnaire. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 4 dans lequel l'étape 3) comprend, pour chaque triplet de valeurs (ar, no, Pmax) et pour chaque instant d'un intervalle de temps, la détermination des satellites non-géostationnaires ayant une valeur minimale d'angle topocentrique (n, aNGS01, aNGS02) supérieure au premier seuil (ar) du triplet (3031), et la simulation des puissances émises par ces satellites dé- 19 [Revendication 6] [Revendication 7] [Revendication 8] [Revendication 9] [Revendication 10] terminés vers la station terrestre et des puissances reçues par la station terrestre (3032) afin de déterminer un débit reçu par la station terrestre à chaque instant (3033) et pour chaque triplet de manière à sélectionner (3034, 3035) l'un des triplets de valeurs (ar, ao, Pmax). Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 5 dans lequel le triplet de valeurs (ar, ao, Pmax) déterminé est celui dont la somme des débits à chaque instant est la plus élevée (3034). Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 5 dans lequel le triplet de valeurs (ar, ao, Pmax) déterminé est celui dont le débit est maximal pour une proportion prédéterminée de l'intervalle de temps (3035). Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une des revendications 4 à 7 comprenant les étapes suivantes effectuées après l'étape 3) : 4) Pour chaque instant de l'intervalle de temps, sélectionner des satellites non-géostationnaires ayant une valeur d'angle topocentrique minimale (a) supérieure au premier seuil angulaire (ar) ; 5) simulation des niveaux de puissance émis par les satellites sélectionnés à l'étape 4) et reçus par la station terrestre cumulés sur l'intervalle de temps ; 6) identification des instants de l'intervalle de temps et ajustement de la troisième valeur maximale de puissance (Pmax) de manière à minimiser l'écart entre une distribution des niveaux de puissance reçus par la station terrestre et cumulés sur l'intervalle de temps et la distribution de référence (REF), la distribution de référence étant supérieure aux niveaux de puissance cumulés. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une des revendications 4 à 8 comprenant également une étape de raffinement (304) consistant à déterminer des seuils angulaires (cxr_i, ao_1) et des seuils de puissance (PR_1, P0_1) supplémentaires de manière à réduire l'écart entre les niveaux de puissance cumulés de l'étape 3) et le niveau de référence (REF), les niveaux de puissance cumulés étant inférieurs au niveau de référence. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une des revendications 1 à 9 dans 20 [Revendication 11] [Revendication 12] [Revendication 13] [Revendication 14] [Revendication 15] [Revendication 16] lequel le procédé est mis en oeuvre pour une pluralité d'hypothèses de positions de stations terrestres. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une des revendications 1 à 10 comprenant également une sélection d'au moins un satellite d'une constellation de satellites non-géostationnaires à chaque instant, les niveaux de puissance reçus par la station résultant de signaux émis par les satellites sélectionnés. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 11 dans Icqucl le satellite sélectionné est celui présentant la puissance reçue par la station terrestre la plus élevée. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 11 dans Icqucl le satellite sélectionné est celui situé au nadir de la station terrestre. Procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon la revendication 13 comprenant une étape de correction de la première (PR), de la seconde (PO) et de la troisième (Pmax) valeurs maximales de puissance d'émission d'un satellite non-géostationnaire de la constellation, situé en un point différent du nadir de la station terrestre, de manière à ce que la puissance reçue par la station terrestre de ce satellite de la constellation soit égale à la puissance reçue par la station terrestre du satellite sélectionné situé au nadir. Programme d'ordinateur comportant des instructions pour l'exécution du procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une quelconque des revendications précédentes, lorsque le programme est exécuté par un processeur. Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de détermination d'une puissance maximale d'émission d'un satellite non-géostationnaire selon l'une quelconque des revendications 1 à 14 ; lorsque le programme est exécuté par un processeur.