FR 3099673 A1 20210205 FR 1908699 A 20190731 Description Titre de l'invention : Procédé de détermination des contraintes d'un système non-géostationnaire vis-à-vis d'un autre système non- géostationnaire [0001] L'invention concerne un procédé de détermination de contraintes opérationnelles à respecter pour une première constellation de satellites non-géostationnaires émettant en direction d'une station terrestre en un point sur Terre vis-à-vis d'une seconde constellation de satellites non-géostationnaires en liaison avec la même station terrestre. Ce procédé permet également de déterminer les contraintes opérationnelles à respecter pour une première station terrestre en liaison avec une première constellation de satellites non-géostationnaires émettant en direction d'un satellite appartenant à une seconde constellation de satellites non-géostationnaires et en liaison avec une seconde station terrestre en un point sur Terre. [0002] La règlementation internationale impose aux différents systèmes non- géostationnaires de se coordonner afin d'éviter de se brouiller les uns les autres. Cette coordination se traduit généralement par la définition et la mise en place de niveaux d'interférences acceptables et par des contraintes opérationnelles qui peuvent induire une réduction de la capacité de certains systèmes. [0003] Le besoin de coordination est une véritable contrainte car il impose de protéger et/ou de se protéger de tous les types de stations des autres constellations sur les zones de service de ces autres constellations, si les zones sont connues, voire même de protéger ou de se protéger pour tous les points sur Terre, sans connaître la direction dans laquelle pointe le terminal victime d'une interférence. [0004] [fig.1] représente un schéma des solutions existantes pour coordonner des systèmes non-géostationnaires. Un satellite NGSO1 d'une première constellation non-géostationnaire CONS1 est en liaison avec une station terrestre S1 en un point sur Terre. Un angle topocentrique 0 est défini entre le satellite NGSO1, la station terrestre S1 et un satellite NGSO2 d'une seconde constellation non-géostationnaire CONS2. Le satellite NGSO2 ne pourra émettre en direction de la station terrestre SI que si l'angle topocentrique 0 est supérieur à un angle prédéterminé, par exemple égal à 10°. [0005] Cette condition sur l'angle topocentrique impose de fortes contraintes sur la constellation interférente CONS2, qui sont supérieures aux besoins du système subissant les interférences de la constellation CONS2, ce système comprenant la première constellation CONS1 et la station terrestre SI. Ces contraintes sont donc assez strictes et non flexibles. [0006] De plus, cela implique un nombre important de contraintes opérationnelles à prendre 2 en compte par le module opérationnel de management des ressources radio. [0007] L'invention vise à remédier aux inconvénients précités de l'art antérieur, plus parti- culièrement elle vise à proposer un procédé de détermination des contraintes opérationnelles pour une première constellation de satellites non-géostationnaires émettant en direction d'une station terrestre en liaison avec une seconde constellation de satellites non-géostationnaires. Elle vise également à proposer un procédé de détermination des contraintes opérationnelles pour une première station terrestre émettant en direction d'un satellite non-géostationnaire d'une constellation de satellites non-géostationnaires en liaison avec une seconde station terrestre. Le procédé selon l'invention permet notamment d'adapter le seuil angulaire topocentrique au-dessous duquel l'émission des satellites de la première constellation ou de la première station terrestre est interdite, et ainsi d'obtenir des contraintes s'adaptant à une situation opérationnelle précise. [0008] Un objet de l'invention est donc un procédé, mis en oeuvre par ordinateur, de déter- mination de contraintes opérationnelles à respecter pour une première constellation de satellites non-géostationnaires émettant en direction d'une station terrestre en un point sur Terre vis-à-vis d'une seconde constellation de satellites non-géostationnaires en liaison avec la station terrestre, les contraintes opérationnelles comprenant au moins une puissance d'émission maximale des satellites de la première constellation, le procédé comprenant les étapes de : - détermination de triplets de valeurs limites de deux angles et d'un rapport interférence à bruit, les deux angles définissant une position d'un satellite de la première constellation par rapport à un axe formé par la station terrestre et un satellite de la seconde constellation et le rapport interférence à bruit étant le rapport entre des interférences émises par la première constellation sur une liaison entre la station terrestre et le satellite de la seconde constellation et du bruit de la liaison, la détermination des triplets étant réalisée de manière à ce qu'une distribution de rapports signal à bruit et interférence cumulés sur un intervalle de temps soit supérieure à une distribution de référence, les rapports signal à bruit et interférence étant les rapports entre un signal utile de la liaison et le bruit et les interférences ; - détermination d'au moins la puissance maximale d'émission d'au moins un satellite de la première constellation à partir des triplets de valeurs. [0009] Selon des modes de réalisation : - Les triplets de valeurs limites sont déterminés par les étapes suivantes 1) sélectionner, pour chaque instant d'un intervalle de temps, des satellites de la première constellation et un satellite de la seconde constellation et déterminer, pour chaque instant de l'intervalle de temps, un triplet de valeurs d'angles définissant une position des satellites sélectionnés de la première constellation vis-à-vis d'un axe 3 formé par la station terrestre et le satellite sélectionné de la seconde constellation et de rapport signal à bruit et interférence, le rapport signal à bruit et interférence étant le rapport signal à bruit et interférence des satellites sélectionnés de la première constellation vis-à-vis de la liaison entre la station terrestre et le satellite sélectionné de la seconde constellation ; 2) déterminer les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence du triplet de ces instants de manière à ce qu'une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps soit supérieure à une distribution de référence ; 3) déterminer des triplets de valeurs aux instants déterminés et une équation de surface paramétrée par les triplets, les triplets de valeurs (0, 4), TIN) étant les angles (0, 4)) définissant la position d'un satellite sélectionné de la première constellation et I/N étant un rapport interférence à bruit sur la liaison déterminé par les rapports signal à bruit et interférence ajustés, les points de cette équation de surface représentant les triplets de valeurs limites (0, (1), I/N). - La sélection des satellites dans l'étape 1) est faite de manière à minimiser un rapport signal à bruit et interférence sur une liaison entre la station terrestre et un satellite de la seconde constellation. - Le procédé comprend les étapes suivantes effectuées après l'étape 3) : 4) Pour chaque instant de l'intervalle de temps et pour un satellite de la seconde constellation en liaison avec la station terrestre, sélectionner des satellites de la première constellation ayant des valeurs d'angle et un rapport interférence à bruit sur la liaison entre la station terrestre et le satellite de la seconde constellation, tels que pour ces valeurs d'angle, le rapport interférence à bruit est inférieur ou égal au rapport interférence à bruit obtenu par l'équation de surface pour ces mêmes valeurs d'angle ; 5) déterminer les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence des satellites de la première constellation sélectionnés de ces instants de manière à minimiser l'écart entre une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps des satellites sélectionnés et une distribution de référence, la distribution des rapports signal à bruit et interférence cumulés étant supérieure à la distribution de référence ; 6) déterminer une deuxième équation de surface paramétrée par les valeurs d'angle et de rapport interférence à bruit ajustés à partir des rapports signal à bruit et interférence ajustés des satellites sélectionnés aux instants déterminés à l'étape 5), les points de cette deuxième équation de surface représentant les triplets de valeurs limites (0, (L.), I/ N). - Le procédé comprend les étapes suivantes effectuées après l'étape 3) : 3') déterminer les valeurs d'angle définissant la position d'un satellite de la première constellation émettant en direction de la station terrestre en liaison avec un satellite de la seconde constellation et de rapport interférence à bruit du satellite de la première constellation sur la liaison entre la station terrestre et le satellite de la seconde constellation ; 3") comparer le rapport interférence à bruit au rapport interférence à bruit déterminé par l'équation de surface pour les valeurs d'angle déterminées en 3'), de sorte que : - si le rapport interférence à bruit est inférieur ou égal à celui déterminé par l'équation de surface, conserver ou augmenter la puissance d'émission du satellite de la première constellation de manière à ce que son rapport interférence à bruit reste inférieur à celui de l'équation de surface ; - si le rapport interférence à bruit est supérieur à celui obtenu avec l'équation de surface, diminuer la puissance d'émission du satellite de la première constellation de manière à ce que son rapport interférence à bruit soit supérieur ou égal à celui de l'équation de surface. - Le procédé est mis en oeuvre pour une pluralité d'hypothèses de positions de stations terrestres. - Les satellites de la première constellation et le satellite de la seconde constellation sélectionnés à l'étape 1) sont ceux minimisant la valeur d'un angle définissant la position des satellites de la première constellation par rapport à l'axe formé entre la station terrestre et un satellite de la seconde constellation. [0010] Un autre objet de l'invention est un procédé, mis en oeuvre par ordinateur, de déter- mination de contraintes opérationnelles à respecter pour une première station terrestre en un point sur Terre, émettant en direction d'un satellite non-géostationnaire d'une constellation de satellites non-géostationnaires en liaison avec une seconde station terrestre vis-à-vis de la liaison entre le satellite et la seconde station terrestre, les contraintes opérationnelles comprenant au moins une puissance d'émission maximale de la première station terrestre, le procédé comprenant les étapes de : - détermination de triplets de valeurs limites (0, (1), I/N) de deux angles et d'un rapport interférence à bruit, les deux angles définissant une position de la première station terrestre par rapport à un axe formé par la seconde station terrestre et le satellite non-géostationnaire et le rapport interférence à bruit étant le rapport entre des interférences émises par la première station terrestre sur la liaison entre la seconde station terrestre et le satellite non-géostationnaire et du bruit de la liaison, la détermination des triplets étant réalisée de manière à ce qu'une distribution de rapports signal à bruit et interférence cumulés sur un intervalle de temps soit supérieure à une distribution de référence, les rapports signal à bruit et interférence étant les rapports entre un signal utile de la liaison et le bruit et les interférences ; - détermination d'au moins la puissance maximale d'émission de la première station terrestre à partir des triplets de valeurs. [0011] Selon des modes de réalisation : - Les triplets de valeurs limites (0, 4), T/N) sont déterminés par les étapes suivantes : 1) sélectionner, pour chaque instant d'un intervalle de temps, des premières stations terrestres et un satellite non-géostationnaire en liaison avec la seconde station terrestre, et déterminer, pour chaque instant de l'intervalle de temps, un triplet de valeurs d'angles définissant une position des premières stations terrestres sélectionnées par rapport à un axe formé par la seconde station terrestre et le satellite non-géostationnaire et de rapport signal à bruit et interférence, le rapport signal à bruit et interférence étant le rapport signal à bruit et interférence des premières stations sélectionnées vis-à-vis de la liaison entre la seconde station terrestre et le satellite non-géostationnaire ; 2) déterminer les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence du triplet de ces instants de manière à ce qu'une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps soit supérieure à une distribution de référence ; 3) déterminer des triplets de valeurs (0, 4), UN) aux instants déterminés et une équation de surface paramétrée par les triplets (0, 4), TIN), les triplets de valeurs étant les angles définissant la position des premières stations terrestres sélectionnées par rapport à l'axe formé par la seconde station terrestre et le satellite non-géostationnaire et I/N étant un rapport interférence à bruit sur la liaison entre la seconde station terrestre et le satellite non-géostationnaire déterminé par les rapports signal à bruit et interférence ajustés, les points de cette équation de surface représentant les triplets de valeurs limites (0, 4), I/N). - La sélection des stations terrestres et du satellite non-géostationnaire à l'étape 1) est réalisée de manière à minimiser un rapport signal à bruit et interférence sur une liaison entre la seconde station terrestre et le satellite non-géostationnaire sélectionné. - Le procédé comprend les étapes suivantes effectuées après l'étape 3) : 4) Pour chaque instant de l'intervalle de temps et pour la seconde station terrestre, sélectionner des premières stations terrestres en des points sur Terre et un satellite non-géostationnaire en liaison avec la seconde station terrestre, tels que les première stations terrestres ont des valeurs d'angle et un rapport interférence à bruit sur la liaison entre la seconde station terrestre et le satellite non-géostationnaire, et que pour ces valeurs d'angle, le rapport interférence à bruit est inférieur ou égal au rapport interférence à bruit obtenu par l'équation de surface pour ces mêmes valeurs d'angle; 5) déterminer les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence des premières stations terrestres sélectionnées de ces instants, de manière à minimiser l'écart entre une distribution des rapports signal à 6 bruit et interférence cumulés sur l'intervalle de temps des premières stations sélectionnées et une distribution de référence, la distribution des rapports signal à bruit et interférence cumulés étant supérieure à la distribution de référence ; 6) déterminer une deuxième équation de surface paramétrée par les valeurs d'angle et de rapport interférence à bruit ajustés à partir des rapports signal à bruit et interférence ajustés des premières stations terrestres sélectionnées aux instants déterminés à l'étape 5), les points de cette deuxième équation de surface représentant les triplets de valeurs limites (0, 4), TIN). - Le procédé comprend les étapes suivantes effectuées après l'étape 3) : 3') déterminer les valeurs d'angle définissant la position d'une première station terrestre émettant en direction d'un satellite non-géostationnaire en liaison avec la seconde station terrestre et le rapport interférence à bruit de la première station terrestre sur la liaison entre la seconde station terrestre et le satellite non-géostationnaire ; 3") comparer le rapport interférence à bruit au rapport interférence à bruit déterminé par l'équation de surface pour les valeurs d'angle déterminées en 3'), de sorte que : - si le rapport interférence à bruit est inférieur ou égal à celui déterminé par l'équation de surface, conserver ou augmenter la puissance d'émission de la première station terrestre de manière à ce que son rapport interférence à bruit reste inférieur à celui de l'équation de surface ; - si le rapport interférence à bruit est supérieur à celui obtenu avec l'équation de surface, diminuer la puissance d'émission de la première station terrestre de manière à ce que son rapport interférence à bruit soit supérieur ou égal à celui de l'équation de surface. - Le procédé est mis en oeuvre pour une pluralité d'hypothèses de positions de la seconde station terrestre. - Les premières stations terrestres et le satellite non-géostationnaire sélectionnés à l'étape 1) sont ceux minimisant une valeur d'angle définissant la position des premières stations terrestres par rapport à un axe formé par la seconde station terrestre et le satellite non-géostationnaire sélectionné. [0012] 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 de triplets de valeurs limites (0, (1), I/N) selon l'invention, lorsque le programme est exécuté par un processeur. [0013] 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 de triplets de valeurs limites (0, 4, TIN) selon l'invention ; lorsque le programme est exécuté par un processeur. [0014] 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 7 qui représentent, respectivement : [0015] [fig.1], déjà décrite, un schéma représentant les contraintes actuelles des constellations de satellites non-géostationnaires vis-à-vis d'un autre système non-géostationnaire selon l'art antérieur ; [0016] [fig.2a] et [0017] [fig.2h], deux schémas du principe du procédé selon l'invention pour deux cas d'application ; [0018] [fig.3], un schéma des étapes du procédé scion un premier mode de réalisation ; [0019] [fig.4], un schéma des étapes du procédé scion un second mode de réalisation ; [0020] [fig.5] et [0021] [fig.6], deux figures représentant le principe de l'étape 2) du procédé selon l' invention ; et [0022] [fig.7], une figure représentant une équation surfacique obtenue par le procédé selon l' invention. [0023] [fig.2a] présente un schéma du principe du procédé selon l'invention dans un premier cas d'application. Cc premier cas d'application représente une liaison descendante. Une constellation CONS V, comprenant plusieurs satellites situés sur une orbite non-géostationnaire et comprenant des moyens aptes à communiquer avec une ou plusieurs stations terrestres, comprend un satellite non-géostationnaire NGSO_V en liaison avec une station terrestre SV située en un point sur Terre. Le satellite NGSO_V envoie et/ou reçoit donc des signaux C dits utiles à la station terrestre SV. [0024] Une autre constellation CONS_I comprend un satellite non-géostationnaire NGSO_I émettant des signaux I vers la station terrestre SV. Ces signaux peuvent brouiller la liaison entre la station terrestre SV et le satellite victime NGSO_V et constituer une source d'interférence pour cette liaison. La constellation CONS_V et le satellite NOSO_V subissant ces interférences sont notés constellation victime et satellite victime dans la suite de la description. [0025] On définit l'angle topocentrique 0 comme l'angle formé entre le satellite NGSO_V et le satellite NGSO_I à partir de la station terrestre SV, et l'angle d'élévation ci) comme l'angle formé entre le plan TAN tangent au sol à la station terrestre SV et l'axe formé par la station terrestre SV et le satellite NGSO_I. Ces deux angles el et 4) permettent de définir une position du satellite NGSO_I par rapport au satellite NGSO_V et à la station SV. [0026] L'objectif du procédé est de définir la valeur maximale de puissance du signal d'interférence I que peut émettre le satellite NGSO_I vers la station SV pour une valeur d'angle topocentrique 0 et une valeur d'angle d'élévation 4) fixées, en générant un niveau d'interférences acceptable vis-à-vis de la constellation CONS_V. Autrement dit, l'objectif du procédé est de déterminer une puissance maximale pour chaque 8 satellite NGSO Ide la constellation CONS_I en tenant compte des interférences générées éventuelles vis-à-vis de satellites NGSO_V d'une autre constellation CONS_V situés dans son voisinage. Pour cela, on détermine des triplets de valeurs limites (0, cp, R) d'un angle topocentrique 0, d'un angle d'élévation (1) et d'un rapport signal à bruit et interférence R de manière à cc qu'une distribution de rapports signal à bruit et interférence R cumulés sur un intervalle de temps soit supérieure à une distribution de référence REF. [0027] Le rapport signal à bruit et interférence R est le rapport entre un signal utile C de la liaison entre la station terrestre SV et le satellite NGSO V de la constellation victime CONS_V et le bruit N et les interférences I émises par la constellation interférente CONS_I sur cette liaison. Autrement dit, le rapport signal à bruit et interférence R est égal à C / [0028] A partir du rapport signal à bruit et interférence R, il est possible de déduire la puissance maximale d'émission d'un satellite NGS0 1 de la constellation interférente CONS_I à partir des triplets de valeurs déterminés et à l'aide de relations théoriques connues du domaine. Plus particulièrement, on peut déterminer pour une valeur donnée d'angle topocentrique 0 et une valeur donnée d'angle d'élévation 4), la puissance maximale d'émission à partir du rapport signal à bruit et interférence du triplet de valeurs comprenant les valeurs d'angles données. [0029] [fig.2b] présente un schéma du principe du procédé selon l'invention dans un second cas d'application, représentant une liaison montante. Une constellation CONS_V, comprenant plusieurs satellites situés sur une orbite non-géostationnaire et comprenant des moyens aptes à communiquer avec une ou plusieurs stations terrestres, comprend un satellite non-géostationnaire NGSO_V en liaison avec une station terrestre SV située en un point sur Terre. Le satellite NGSO_V envoie et/ou reçoit donc des signaux C dits utiles à la station terrestre SV. [0030] Une autre station terrestre SI, dite station interférente, est en liaison avec une seconde constellation CONS_I de satellites non-géostationnaires, notamment avec un satellite non-géostationnaire NGSO_I. Cette autre station terrestre SI émet également des signaux I vers le satellite NGSO_V en liaison avec la station terrestre SV. Ces signaux I peuvent brouiller la liaison entre la station terrestre SV et le satellite NGSO_V, dit victime, et constituer une source d'interférence pour cette liaison. [0031] Afin d'identifier la position de la station terrestre interférente SI, on définit un angle 0 formé entre la station SV et la station interférence SI à partir du satellite non-géostationnaire victime NGSO_V, et un angle d'élévation (L) comme l'angle formé entre le plan TAN tangent à l'orbite de la constellation victime CONS_V et au satellite NGSO_V et l'axe formé par la station terrestre interférente SI et le satellite NGSO_V. Ces deux angles 0 et ci) permettent de définir la position de la station SI par rapport au 9 satellite NGSO V et à la station SV. [0032] L'objectif du procédé, pour ce second cas d'application, est de définir la valeur maximale de puissance du signal d'interférence i que peut émettre la station SI vers le satellite NGSO_V pour une valeur d'angle 0 et une valeur d'angle 4) fixées, en générant un niveau d'interférences acceptable vis-à-vis de la constellation CONS_V. Autrement dit, l'objectif du procédé est de déterminer une puissance maximale pour chaque station terrestre SI en tenant compte des interférences générées éventuelles vis-à-vis de satellites NGS0 V d'une constellation CONS V en liaison avec une station terrestre SV située dans son voisinage. Pour cela, on détermine des triplets de valeurs limites (0, cp, R) des angles 0, (I) et d'un rapport signal à bruit et interférence R de manière à ce qu'une distribution de rapports signal à bruit et interférence R cumulés sur un intervalle de temps soit supérieure à une distribution de référence REF. [0033] Le rapport signal à bruit et interférence R est le rapport entre un signal utile C de la liaison entre la station terrestre SV et le satellite NGSO V de la constellation victime CONS_V et le bruit N et les interférences I émises par la station interférente SI sur cette liaison. Autrement dit, le rapport signal à bruit et interférence R est égal à C / (N-FI). [0034] A partir du rapport signal à bruit et interférence R, il est possible de déduire la puissance maximale d'émission de la station SI à partir des triplets de valeurs déterminés et à l'aide de relations théoriques connues du domaine. Plus particulièrement, on peut déterminer pour des valeurs données d'angles 0 et 4), la puissance maximale d'émission à partir du rapport signal à bruit et interférence du triplet de valeurs comprenant les valeurs d'angles données. [0035] [fig.3] présente un schéma des étapes du procédé selon un premier mode de réa- lisation. On se place dans le premier cas d'application (figure 2a) pour la description des figures 3 à 7. [0036] Dans une première étape 201, on sélectionne, pour chaque instant d'un intervalle de temps et pour une station terrestre SV en un point sur Terre, des satellites de la constellation interférence CONS_I et un satellite NGSO_V de la constellation victime CONS_V selon un critère de sélection. Dans l'exemple décrit, on sélectionne les satellites de la constellation interférente CONS_I et un satellite NGSO_V qui minimisent le rapport signal à bruit et interférence R sur la liaison entre la station terrestre SV et le satellite de la constellation victime CONS_V. Cependant, on peut également sélectionner le satellite NGSO_V de la constellation victime présentant la plus haute élévation dans cette constellation vis-à-vis de la station terrestre SV. On peut également sélectionner un satellite interférent CONS_I ayant la plus longue durée de visibilité vis-à-vis de la station terrestre SV et/ou d'un satellite particulier de la constellation victime CONS_V. [0037] On détermine ensuite, pour ces mêmes instants, un triplet de valeurs (0, cp, R) d'angle topocentrique 0, d'angle d'élévation 4) et de rapport signal à bruit et interférence R pour ces satellites sélectionnés par rapport à la station terrestre SV. L'angle topocentrique 0 est l'angle minimal formé par le satellite NGSO V de la constellation victime CONS V sélectionné, la station terrestre SV et les satellites sélectionnés NGSO 1 de la constellation interférente CONS I. L'angle d'élévation (1) est l'angle minimal formé entre le plan tangent TAN au sol à la station terrestre SV et l'axe formé entre la station terrestre SV et les satellites sélectionnés NGSO I de la constellation interférente CONS I. Le rapport signal à bruit et interférence R est le rapport signal à bruit et interférence des satellites interférents NGSO I de la constellation interférente CONS I sur la liaison entre la station terrestre SV et le satellite victime NGSO V de la constellation victime CONS V. Ce rapport R est déterminé par simulation. [0038] Puis dans une seconde étape 202, on détermine les instants de l'intervalle de temps pour lesquels on peut ajuster la valeur de rapport signal à bruit et interférence R du triplet de valeurs R) de ces instants de manière à ce qu' une distribution des rapports signal à bruit et interférence R cumulés sur l'intervalle de temps soit supérieure à la distribution de référence REF. [0039] [fig.5] et [Fig. 6] illustrent l'étape 202 du procédé. [0040] La figure 5 représente le pourcentage de temps %t de l'intervalle de temps en fonction du rapport signal à bruit et interférence R = C / La courbe REF re- présente la distribution de référence et la courbe Ni représente les rapports signal à bruit et interférence cumulés des satellites des constellations victime et interférents sélectionnés vis-à-vis de la station terrestre, tandis que la courbe N2 représente ces mêmes rapports signal à bruit et interférence après l'identification des instants de temps dont la valeur du rapport signal à bruit et interférence est ajustée pour que la courbe NI soit supérieure à la distribution de référence REF. [0041] La figure 6 représente le rapport signal à bruit et interférence R en fonction des différents instants t de l'intervalle de temps. Ce sont ces rapports signal à bruit et interférence cumulés sur l'intervalle de temps (donc cumulés sur tous les instants t) qui permettent d'obtenir la courbe NI de la figure 5. Les instants de temps identifiés INST sont représentés sur la figure 6 et leur valeur de rapport signal à bruit et interférence R est modifiée de manière obtenir la courbe N2 de la figure 5. [0042] Pour déterminer ces instants t de l'intervalle de temps à partir des figures 5 et 6, on suppose que la courbe Ni de la figure 5 suit une loi de probabilité suivante : [0043] PN1 ( X = pn (I) [0044] où X représente une variable aléatoire du rapport signal à bruit et interférence R, x' une valeur particulière d'un rapport signal à bruit et interférence, pi, la probabilité de temps correspondant à la valeur de rapport signal à bruite et interférence xi, et n est un 11 entier inférieur ou égal au nombre total d'instants de temps de l'intervalle de temps. [0045] On suppose également que la distribution de référence REF suit la loi de probabilité suivante : [0046] Parr' ( X x) p ( 2 ) [0047] avec p une probabilité de temps correspondant à la valeur de rapport signal à bruit et interférence x et X une variable aléatoire représentant un rapport signal à bruit et interférence. [0048] Si la courbe NI devient inférieure à la distribution de référence REF pour toutes les probabilités p, alors cela signifie qu'il existe un entier i et ni tels que : [0049] PNI(X -..2; xi) = pi et PREF(X x',) = pi avec xi X,') = pin avec pi ) et un rapport interférence à bruit (I/N) sur la liaison entre la station terrestre SV et le satellite non-géostationnaire de la constellation victime, et que pour ces valeurs d'angle, le rapport interférence à bruit est inférieur ou égal au rapport interférence à bruit obtenu par l'équation de surface pour ces mêmes valeurs d'angle. [0085] Puis dans l'étape 205, on détermine les instants de l'intervalle de temps et on ajuste la valeur de rapport signal à bruit et interférence (R) des premières stations terrestres 16 sélectionnées de ces instants, de manière à minimiser l'écart entre une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps des premières stations sélectionnées et une distribution de référence (REF), la distribution des rapports signal à bruit et interférence cumulés étant supérieure à la distribution de référence. [0086] Enfin, dans l'étape 206, on détermine une deuxième équation de surface paramétrée par les valeurs d'angle et de rapport interférence à bruit ajustés à partir des rapports signal à bruit et interférence ajustés des premières stations terrestres sélectionnées aux instants déterminés à l'étape 205. Les points de cette deuxième équation de surface représentent les triplets de valeurs limites (0, 4), TIN). [0087] Pour le cas de la figure 3, on détermine dans l'étape 203', les valeurs d'angle (0, 4)) définissant la position d'une première station terrestre émettant en direction d'un satellite non-géostationnaire en liaison avec la station terrestre SV et le rapport interférence à bruit (I/N) de la première station terrestre sur la liaison entre la station terrestre et le satellite non-géostationnaire. [0088] Puis on compare dans l'étape 203" le rapport interférence à bruit au rapport in- terférence à bruit déterminé par l'équation de surface pour les valeurs d'angle déterminées dans l'étape 203', de sorte que : - si le rapport interférence à bruit est inférieur ou égal à celui déterminé par l'équation de surface, conserver ou augmenter la puissance d'émission de la première station terrestre de manière à ce que son rapport interférence à bruit reste inférieur à celui de l'équation de surface ; - si le rapport interférence à bruit est supérieur à celui obtenu avec l'équation de surface, diminuer la puissance d'émission de la première station terrestre de manière à ce que son rapport interférence à bruit soit supérieur ou égal à celui de l'équation de surface. [0089] Comme pour le premier cas d'application, on peut mettre en oeuvre le procédé pour une pluralité d'hypothèses de positions de la station terrestre SV. [0090] De même, les premières stations terrestres et le satellite non-géostationnaire sé- lectionnés à l'étape 201 peuvent être ceux minimisant une valeur d'angle, par exemple 0, définissant la position des premières stations terrestres par rapport à un axe formé par la seconde station terrestre et le satellite non-géostationnaire sélectionné. [0091] Le procédé selon l'invention a été décrit dans le cadre de deux constellations non- géostationnaires. Cependant, le procédé n'est pas limité à cela, il peut également être utilisé dans un contexte plus général faisant intervenir une constellation non-géostationnaire et une constellation géostationnaire. Cela peut par exemple être dans le cas où une constellation géostationnaire (satellite de la constellation ou station terrestre en liaison avec la constellation géostationnaire) émet des signaux de brouillage sur une 17 liaison entre une station terrestre sur Terre et une constellation non-géostationnaire. Cela peut également être dans le cas où une constellation non-géostationnaire (satellite de la constellation ou station terrestre en liaison avec la constellation non-géostationnaire) émet des signaux de brouillage sur une liaison avec une station terrestre sur Terre et une constellation géostationnaire. [0092] 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. [0093] 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-Spccific Tntegrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). [0094] 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). [0095] 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. [0096] 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 de - 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 Access 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 (I/O, acronyme de - Input/Output » en anglais) adaptée à transmettre et à recevoir des données. 18 [0097] 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. 19 [Revendication I] [Revendication 2] REVENDICATIONSProcédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles à respecter pour une première constellation de satellites non-géostationnaires (CONS_I) émettant en direction d'une station terrestre (SV) en un point sur Terre vis-à-vis d'une seconde constellation de satellites non-géostationnaires (CONS_V) en liaison avec la station terrestre, les contraintes opérationnelles comprenant au moins une puissance d'émission maximale des satellites de la première constellation, le procédé comprenant les étapes de : - détermination de triplets de valeurs limites (0, 4), I/N) de deux angles (B, 4)) et d'un rapport interférence à bruit, les deux angles (B, 4)) définissant une position d'un satellite (NGSO _I) de la première constellation (CONS_I) par rapport à un axe formé par la station terrestre (SV) et un satellite (NGSO_V) de la seconde constellation (CONS_V) et le rapport interférence à bruit (UN) étant le rapport entre des interférences (I) émises par la première constellation sur une liaison entre la station terrestre (SV) et le satellite (NGSO_V) de la seconde constellation (CONS_V) et du bruit (N) de la liaison, la détermination des triplets étant réalisée de manière à ce qu'une distribution de rapports signal à bruit et interférence (R) cumulés sur un intervalle de temps soit supérieure à une distribution de référence (REF), les rapports signal à bruit et interférence (R) étant les rapports entre un signal utile (C) de la liaison et le bruit (N) et les interférences (I) ; - détermination d'au moins la puissance maximale d'émission d'au moins un satellite de la première constellation à partir des triplets de valeurs. Procédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles selon la revendication 1 dans lequel les triplets de valeurs limites (0, (I), UN) sont déterminés par les étapes suivantes : 1) sélectionner (201), pour chaque instant d'un intervalle de temps, des satellites de la première constellation et un satellite de la seconde constellation et déterminer, pour chaque instant de l'intervalle de temps, un triplet (B, 4), R) de valeurs d'angles (0, (k) définissant une position des satellites sélectionnés de la première constellation (CONS_I) vis-à-vis d'un axe formé par la station terrestre (SV) et le satellite sélectionné de la seconde constellation (CONS_V) et de rapport signal à bruit et interférence (R), le rapport signal à bruit et interférence étant le 20 [Revendication 3] [Revendication 4] rapport signal à bruit et interférence des satellites sélectionnés de la première constellation vis-à-vis de la liaison entre la station terrestre et le satellite sélectionné de la seconde constellation ; 2) déterminer (202) les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence du triplet (0, (1), R) de ces instants de manière à ce qu'une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps soit supérieure à une distribution de référence (REF) ; 3) déterminer (203) des triplets de valeurs (0, 4), 1/N) aux instants déterminés et une équation de surface paramétrée par les triplets (0, (1), I/ N), les triplets de valeurs (0, (1), I/N) étant les angles (0, t) définissant la position d'un satellite sélectionné de la première constellation (CONS I) et I/N étant un rapport interférence à bruit sur la liaison déterminé par les rapports signal à bruit et interférence (R) ajustés, les points de cette équation de surface représentant les triplets de valeurs limites (0, (1), TIN). Procédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles selon la revendication 2 dans lequel la sélection des satellites dans l'étape 1) est faite de manière à minimiser un rapport signal à bruit et interférence (R) sur une liaison entre la station terrestre et un satellite de la seconde constellation. Procédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles selon l'une des revendications 2 ou 3 comprenant les étapes suivantes effectuées après l'étape 3) : 4) Pour chaque instant de l'intervalle de temps et pour un satellite de la seconde constellation en liaison avec la station terrestre, sélectionner (204) des satellites de la première constellation ayant des valeurs d'angle (0, ci)) et un rapport interférence à bruit (T/N) sur la liaison entre la station terrestre (SV) et le satellite de la seconde constellation (NGSO_V), tels que pour ces valeurs d'angle, le rapport interférence à bruit est inférieur ou égal au rapport interférence à bruit obtenu par l'équation de surface pour ces mêmes valeurs d'angle; 5) déterminer (205) les instants de l'intervalle de temps et ajuster la valeur de rapport signal à bruit et interférence (R) des satellites de la première constellation sélectionnés de ces instants de manière à minimiser l'écart entre une distribution des rapports signal à bruit et interférence cumulés sur l'intervalle de temps des satellites sélectionnés et une distribution de référence (REF), la distribution des rapports signal à 21 [Revendication 5] [Revendication 6] [Revendication 7] bruit et interférence cumulés étant supérieure à la distribution de référence ; 6) déterminer (206) une deuxième équation de surface paramétrée par les valeurs d'angle et de rapport interférence à bruit ajustés à partir des rapports signal à bruit et interférence ajustés des satellites sélectionnés aux instants déterminés à l'étape 5), les points de cette deuxième équation de surface représentant les triplets de valeurs limites (0, (1), 1/ N). Procédé, mis en oeuvre par ordinateur, de détermination dc contraintes opérationnelles selon l'une des revendications 2 ou 3 dans lequel le procédé comprend les étapes suivantes effectuées après l'étape 3) : 3') déterminer (203') les valeurs d'angle (0, (p) définissant la position d'un satellite de la première constellation émettant en direction de la station terrestre en liaison avec un satellite de la seconde constellation et de rapport interférence à bruit (UN) du satellite de la première constellation sur la liaison entre la station terrestre et le satellite de la seconde constellation ; 3") comparer (203") le rapport interférence à bruit au rapport interférence à bruit déterminé par l'équation dc surface pour les valeurs d'angle déterminées en 3'), de sorte que : - si le rapport interférence à bruit est inférieur ou égal à celui déterminé par l'équation de surface, conserver ou augmenter la puissance d'émission du satellite de la première constellation de manière à ce que son rapport interférence à bruit reste inférieur à celui de l'équation de surface ; - si le rapport interférence à bruit est supérieur à celui obtenu avec l'équation de surface, diminuer la puissance d'émission du satellite de la première constellation de manière à ce que son rapport interférence à bruit soit supérieur ou égal à celui de l'équation de surface. Procédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles selon l'une des revendications précédentes dans lequel le procédé est mis en oeuvre pour une pluralité d'hypothèses de positions de stations terrestres. Procédé, mis en oeuvre par ordinateur, de détermination de contraintes opérationnelles selon l'une des revendications 2 à 6 dans lequel les satellites de la première constellation et le satellite de la seconde constellation sélectionnés à l'étape 1) sont ceux minimisant la valeur d'un angle (0) définissant la position des satellites de la première [Revendication 8] constellation par rapport à l'axe formé entre la station terrestre et un satellite de la seconde constellation. [Revendication 9] Programme d'ordinateur comportant des instructions pour l'exécution du procédé de détermination de triplets de valeurs limites (0, (1), I/N) selon l'une quelconque des revendications 1 à 7, 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 de triplets de valeurs limites (0, ch, TIN) selon l'une quelconque des revendications 1 à 7, lorsque le programme est exécuté par un processeur.