FR 3147470 A1 20241004 FR 2302810 A 20230328 L’invention concerne une méthode (100) pour contrôler les performances, en termes de puissance de bruit rayonnée, d'une charge utile d'un satellite de télécommunications basée sur une antenne active multifaisceaux. La méthode comporte notamment la détermination (105) d'une densité spectrale de puissance de bruit d'intermodulation rayonnée au sol par l'ensemble des sources de l'antenne à partir : du plan d'allocation des porteuses et des puissances RF associées permettant de remplir la mission de télécommunications du satellite, des coefficients de pondération du réseau de formation de faisceaux permettant de former les faisceaux radio répondant à la mission, les caractéristiques de linéarité des amplificateurs associés aux différentes sources de l'antenne, et les diagrammes de rayonnement des différentes sources. La méthode permet de déterminer (109) un rapport de puissance rayonnée entre signal utile et bruit d'intermodulation ("C/IM rayonné") pour au moins une porteuse et pour au moins une zone géographique d'intérêt. Figure pour l’abrégé : Fig. 4 Validation de la charge utile d’un satellite de télécommunications basée sur un calcul de rapport de puissance rayonnée entre signal utile et bruit d'intermodulation Domaine de l’invention La présente invention appartient au domaine du contrôle des performances de la charge utile d’un satellite de télécommunications. Plus particulièrement, l’invention concerne une méthode de validation de la configuration de la charge utile d’un satellite de télécommunications basée sur une estimation d’un rapport de puissance rayonnée entre signal utile et bruit d'intermodulation (« C/IM rayonné »). Etat de la technique Certains satellites de télécommunication sont affectés à la transmission d'informations vers des points de réception fixes, on parle alors de « service fixe par satellite » (FSS pour « Fixed-Satellite Service » dans la littérature anglo-saxonne) ou de diffusion DTH (acronyme anglais pour "Direct-To-Home", c'est à dire « directement vers le particulier »). C’est le cas par exemple pour la diffusion de programmes de télévision ou pour la téléphonie internationale. D’autres satellites de télécommunications sont affectés à des services de connexion à Internet à haut débit correspondant à des missions de type HTS (acronyme anglais pour « High-Throughput Satellite », « satellite haut débit » en français). Pour une même quantité de spectre radio alloué, un satellite répondant à une mission de type HTS fournit nettement plus de débit qu'un satellite répondant à une mission de type FSS. Cette augmentation significative du débit pour les missions de type HTS est obtenue par une réutilisation de canaux fréquentiels dans des faisceaux radio étroitement focalisés et agencés pour couvrir une région géographique d'intérêt. L’augmentation du débit s’accompagne également d’une augmentation de la flexibilité de couverture (variations des faisceaux au cours du temps). Les satellites répondant à une mission de type HTS embarquent généralement une antenne active multifaisceaux. Une telle antenne comporte un réseau de formation de faisceaux (BFN, pour « Beam-Forming Network » en anglais) relié à un réseau de sources. Des coefficients de pondération en phase et en amplitude sont associés respectivement aux différentes sources et aux différents canaux fréquentiels afin de former les faisceaux radio souhaités. Chaque source est associée à un amplificateur adapté pour amplifier le signal devant être émis par la source. La charge utile d’un satellite de télécommunications peut être contrôlée et mise à jour par un centre de contrôle de mission (MCC pour « mission control center » en anglais), notamment pour répondre à des changements dans une mission en cours, voire pour réponde à une nouvelle mission significativement différente de la mission en cours. Pour contrôler les performances de la charge utile d’un satellite de télécommunications, il convient de contrôler le bilan de liaison radio entre le satellite et une station de réception située dans une zone géographique d’intérêt à la surface de la Terre. La non-linéarité des amplificateurs de l’antenne active multifaisceaux génère un bruit d’intermodulation. Le rapport de puissance rayonnée entre le signal utile et le bruit d'intermodulation (« rapport C/IM ») est un élément important contribuant au calcul du bilan de liaison radio d’une communication entre un satellite et une station terrestre. Toutefois, il n’existe pas actuellement de solution satisfaisante pour calculer rapidement un rapport C/IM pour une antenne active multifaisceaux comportant un grand nombre de sources (plusieurs dizaines voire plusieurs centaines de sources) et configurée pour former un grand nombre de faisceaux (plusieurs centaines voire plusieurs milliers de faisceaux). Le document « Calculation of third order intermodulation spectrum », Y.H. LAU et al., décrit une méthode de calcul rapide d’une densité spectrale de bruit d’intermodulation de manière conduite (et non pas de manière rayonnée). La présente invention a pour objectif de remédier à tout ou partie des inconvénients de l’art antérieur, notamment ceux exposés ci-avant. A cet effet, et selon un premier aspect, il est proposé par la présente invention, une méthode mise en œuvre par un ordinateur pour valider une configuration d’une charge utile d’un satellite de télécommunications. Le satellite est configuré pour émettre un signal avec une antenne active multifaisceaux comportant un réseau de formation de faisceaux. Le réseau de formation de faisceaux est relié à une pluralité de sources et adapté à former simultanément plusieurs faisceaux d'intérêt desservant des zones géographiques à la surface de la Terre. Le réseau de formation de faisceaux est configuré pour recevoir en entrée, pour chaque faisceau d’intérêt à former, un signal d'entrée comportant un ou plusieurs canaux fréquentiels avec des niveaux de puissance respectifs, et pour transmettre à chaque source une combinaison linéaire des signaux d’entrée pondérés par des coefficients de pondération, chaque source étant reliée à un amplificateur. La méthode comporte pour chaque amplificateur : une détermination d'un « masque normalisé » de densité spectrale de puissance de bruit d'intermodulation à partir des canaux fréquentiels composant les signaux d'entrée et des niveaux de puissance respectifs desdits canaux fréquentiels, une détermination d'un point d'opération correspondant à une puissance de sortie de l'amplificateur, le point d'opération étant déterminé en fonction des niveaux de puissance des canaux fréquentiels des signaux d'entrée, et en fonction des coefficients de pondération associés audit amplificateur, une détermination d'une puissance totale de bruit d'intermodulation en fonction du point d'opération et de caractéristiques de linéarité de l'amplificateur, une détermination d'un « masque absolu » de densité spectrale de puissance de bruit d'intermodulation par une pondération du masque normalisé avec la puissance totale de bruit d'intermodulation. La méthode comporte pour chaque source : une détermination d'un « masque rayonné » de densité spectrale de puissance de bruit d'intermodulation à partir du masque absolu de densité spectrale de puissance de bruit d'intermodulation de l'amplificateur auquel la source est reliée, et à partir d'un diagramme de rayonnement de la source. La méthode comporte en outre : une détermination d'un « masque rayonné total » de densité spectrale de puissance de bruit d'intermodulation par une sommation en puissance des masques rayonnés des différentes sources. Pour au moins un canal fréquentiel d'intérêt, et pour au moins une zone géographique d'intérêt : une détermination d'une puissance rayonnée de bruit d'intermodulation par une intégration, sur une sous-bande de fréquences correspondant audit canal fréquentiel, du masque rayonné total de densité spectrale de puissance de bruit d'intermodulation pour la zone géographique d'intérêt, une détermination, pour le canal fréquentiel d'intérêt, d'une puissance isotrope équivalente rayonnée par l'antenne active multifaisceaux, une détermination d'un rapport « C/IM rayonné » entre la puissance isotrope équivalente rayonnée et la puissance rayonnée de bruit d'intermodulation, une vérification d’un critère prédéterminé en fonction du rapport « C/IM rayonné » ainsi déterminé. La méthode selon l’invention fournit ainsi une méthode rapide et efficace pour une évaluation analytique de la performance du rapport de puissance rayonnée entre signal utile et bruit d'intermodulation. Une configuration de la charge utile d’un satellite de télécommunications peut ainsi être validée comme satisfaisant à un cahier des charges. Plusieurs configurations peuvent avantageusement être testées et comparées pour sélectionner la configuration répondant au mieux à un besoin. La solution proposée permet avantageusement de calculer une performance rayonnée au sol par porteuse en prenant en compte : le plan de porteuses de la mission de télécommunication considérée (les canaux fréquentiels composant les signaux d’entrée du réseau de formation de faisceaux), les poids utilisés par le réseau de formation de faisceaux de l’antenne active (les coefficients de pondération associés aux différents canaux fréquentiels et aux différentes sources), les caractéristiques de linéarité des amplificateurs associés respectivement aux différentes sources de l’antenne active, et les diagrammes de rayonnement des différentes sources de l’antenne active. La méthode selon l’invention est applicable à tout type d’antenne active multifaisceaux. Il peut s’agir d’une antenne d’émission ou d’une antenne de réception. L’antenne peut comporter un ou plusieurs réseaux de formation de faisceaux de type analogique ou de type numérique. Il peut également s’agir d’une antenne avec ou sans réflecteur. L’antenne active peut présenter plusieurs réflecteurs et/ou plusieurs réseaux de sources. Les amplificateurs associés aux différentes sources peuvent être de type SSPA (acronyme anglais pour « Solid State Power Amplifier », « amplificateur à état solide » en français) ou à base de tubes. Avantageusement, le temps d’exécution de la méthode reste très court et peu impacté par le nombre de canaux fréquentiels, le nombre de sources et le nombre de faisceaux formés. La méthode selon l’invention peut être appliquée à des outils de conception, de paramétrage ou de validation de la configuration de la charge utile d’un satellite de télécommunications, ou à des outils de modélisation des performances de la charge utile d’un satellite de télécommunications. La méthode peut être mise en œuvre pour des satellites évoluant sur différents types d’orbite : orbite géostationnaire (« Geostationary Orbit » ou GEO en anglais), orbite moyenne de type MEO (acronyme anglais pour « Medium Earth Orbit »), orbite basse de type LEO (acronyme anglais pour « Low Earth Orbit »), orbite de Molnia, etc. Dans des modes particuliers de mise en œuvre, l’invention peut comporter en outre l’une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. Dans des modes particuliers de mise en œuvre, le masque normalisé de densité spectrale de puissance de bruit d'intermodulation est commun à tous les amplificateurs. Dans des modes particuliers de mise en œuvre, le masque normalisé de densité spectrale de puissance de bruit d'intermodulation est déterminé pour un amplificateur en fonction des coefficients de pondération associés audit amplificateur. Dans des modes particuliers de mise en œuvre, la détermination du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de produits d'intermodulation d'ordre trois générés par les canaux fréquentiels des signaux d'entrée. Dans des modes particuliers de mise en œuvre, la détermination du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte en outre un calcul de produits d'intermodulation d'ordre impair strictement supérieur à trois générés par les canaux fréquentiels des signaux d'entrée. Dans des modes particuliers de mise en œuvre, la détermination du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte l’application d’une valeur de correction visant à compenser des produits d'intermodulation d'ordre strictement supérieur à trois non pris en compte. Dans des modes particuliers de mise en œuvre, la détermination du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de normalisation en fonction du résultat d’une intégration, sur une bande de fréquences, de la densité spectrale de puissance de bruit d'intermodulation. Dans des modes particuliers de mise en œuvre, la détermination du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de normalisation en fonction du résultat d’une intégration, sur une bande de fréquences, de la densité spectrale de puissance de bruit d'intermodulation pour des valeurs de puissance comprises dans une plage de valeurs prédéterminée. Dans des modes particuliers de mise en œuvre, la méthode comporte en outre l’application d’une valeur de correction au rapport « C/IM rayonné ». La valeur de correction est déterminée en fonction d'un facteur de réutilisation des canaux fréquentiels pour former les faisceaux d'intérêt. Dans des modes particuliers de mise en œuvre, la bande de fréquences est l’une des bandes W, Q, V, Ka, Ku, L, S, C. Selon un deuxième aspect, la présente invention concerne un produit programme d’ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes d’une méthode de validation d’une configuration d’une charge utile d’un satellite de télécommunications selon l’un quelconque des modes de mise en œuvre précédents. Selon un troisième aspect, la présente invention concerne un support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes d’une méthode de validation d’une configuration d’une charge utile d’un satellite de télécommunications selon l’un quelconque des modes de mise en œuvre précédents. Présentation des figures L’invention sera mieux comprise à la lecture de la description suivante, donnée à titre d’exemple nullement limitatif, et faite en se référant aux figures suivantes : une représentation schématique d’un système de télécommunications par satellite, une représentation schématique d’une antenne active multifaisceaux d’un satellite de télécommunications, une représentation schématique d’un ordinateur configuré pour mettre en œuvre une méthode selon l’invention pour valider la configuration de la charge utile d’un satellite de télécommunications, une représentation schématique des principales étapes d’un premier mode de mise en œuvre de la méthode selon l’invention, une représentation schématique des principales étapes d’un deuxième mode de mise en œuvre de la méthode selon l’invention. Dans ces figures, des références identiques d’une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté, les éléments représentés ne sont pas nécessairement à une même échelle, sauf mention contraire. Description détaillée d’un mode de réalisation de l’invention La représente schématiquement un système de télécommunications par satellite. Tel qu’illustré sur la , le système comporte au moins un satellite 30 de télécommunications en orbite autour de la Terre 50. Dans l’exemple considéré, la mission du satellite 30 de télécommunications correspond à un service de connexion à Internet à haut-débit (service de type HTS). On se place de manière non limitative dans le cas où le satellite 30 est en orbite géostationnaire (orbite de type GEO). Rien n’exclut cependant, suivant d’autres exemples, de considérer une orbite terrestre non géostationnaire, telle qu’une orbite de basse altitude (orbite de type LEO) ou une orbite de moyenne altitude (orbite de type MEO). De façon conventionnelle, le satellite 30 de télécommunications peut permettre des échanges de données entre deux stations passerelles 40 via des liens 41 de communication radio établis respectivement entre chaque station passerelle 40 et le satellite 30 de télécommunications. Selon un autre exemple, le satellite 30 de télécommunications peut permettre des échanges de données entre un utilisateur 42 et une station passerelle 40 via d’une part un lien 41 de communication radio établi entre le satellite 30 et l’utilisateur 42 et d’autre part un lien 43 de communication radio établi entre le satellite 30 et la station passerelle 40. Selon encore un autre exemple, le satellite 30 de télécommunications peut permettre des échanges de données entre deux utilisateurs 42 via des liens 43 de communication radio établis respectivement entre chaque utilisateur 42 et le satellite 30 de télécommunications. Chaque lien 41, 43 de communication radio peut supporter une direction descendante (« downlink » dans la littérature anglo-saxonne) pour l’émission de données par le satellite 30 et/ou une direction montante (« uplink » dans la littérature anglo-saxonne) pour la réception de données par le satellite 30. Les émissions et/ou les réceptions sur les liens 41, 43 de communication radio s’effectuent par exemple dans une ou plusieurs bandes de fréquences parmi les bandes W, Q, V, Ka, Ku, L, S, C, etc. Dans l’exemple considéré et illustré à sur la , le satellite 30 de télécommunications est contrôlé par un centre 17 de contrôle de mission localisé au sol, à la surface de la Terre 50. Le centre 17 de contrôle de mission permet notamment de mettre à jour la charge utile du satellite 30 de télécommunications pour répondre à d’éventuels changements dans la mission en cours, voire pour répondre à une nouvelle mission significativement différente de la mission en cours. De façon conventionnelle, la charge utile du satellite 30 de télécommunications peut être configurée à distance par l'envoi de commandes à travers des signaux radio émis sur un lien 16 de contrôle établi entre le satellite 30 de télécommunications et le centre 17 de contrôle de mission. Tel qu’illustré sur la , le centre 17 de contrôle de mission comporte une antenne 15 pour communiquer avec le satellite 30 de télécommunication sur le lien 16 de contrôle. Le centre 17 de contrôle de mission comporte également un ordinateur 10 configuré pour mettre en œuvre une méthode de validation de la charge utile du satellite 30 de télécommunications. Il est en effet intéressant de pouvoir contrôler les performances de la charge utile du satellite 30 de télécommunications en termes de bilan de liaison radio entre le satellite 30 et une station de réception (correspondant par exemple à une station passerelle 40 ou à un utilisateur 42) située dans une zone géographique d’intérêt à la surface de la Terre 50. Plus particulièrement, il est intéressant de pouvoir évaluer, pour une porteuse donnée et pour une zone géographique d’intérêt donnée, un rapport de puissance rayonnée entre signal utile et bruit d'intermodulation (« rapport C/IM »). Le rapport C/IM est en effet un élément important dans le bilan de liaison radio d’une communication établie entre un satellite et une station de réception terrestre. Pour répondre aux besoins de sa mission HTS, notamment en termes de flexibilité de débit et de couverture géographique, le satellite 30 de télécommunications embarque une antenne 20 active multifaisceaux. La représente schématiquement une telle antenne 20 active multifaisceaux. L’antenne 20 comporte notamment un réseau 25 de formation de faisceaux relié à un réseau 23 de sources 24. La méthode selon l’invention est particulièrement bien adaptée au cas où le nombre de sources 24 est important, par exemple au moins égal à cinquante (M ≥ 50), voire au moins égal à cent (M ≥ 100). Rien n’empêcherait toutefois d’appliquer l’invention à des cas où le nombre de sources serait moins important. Dans l’exemple considéré, le réseau 23 de sources comporte deux-cent sources 24 (M = 200). Dans l’exemple considéré et illustré sur la , l’antenne 20 comporte également un réflecteur 22. Le réseau 23 de sources 24 est par exemple positionné au niveau du point focal du réflecteur 22 (antenne de type AFR, acronyme anglais de « Array Fed Reflector », en français « réflecteur alimenté par un réseau »). Selon un autre exemple, le réseau 23 de sources 24 est décalé par rapport à un point focal du réflecteur 22 (antenne de type DAFR, acronyme anglais de « Defocused Array Fed Reflector », en français « réflecteur alimenté par un réseau défocalisé »). Rien n’exclut toutefois d’utiliser une antenne 20 sans réflecteur (antenne de type DRA, acronyme anglais de « Direct Radiating Array », en français « réseau à rayonnement direct »), ou une antenne comportant plusieurs réflecteurs. Rien n’empêcherait non plus de considérer une antenne 20 avec plusieurs réflecteurs. Le réseau 25 de formation de faisceaux peut être « adaptatif », c’est-à-dire que les faisceaux formés par le réseau 25 de formation de faisceaux peuvent varier au cours du temps. Le réseau 25 de formation de faisceaux peut être analogique ou numérique. Un réseau de formation de faisceaux analogique effectue des traitements sur des signaux analogiques (électriques ou micro-ondes), par exemple avec des circuits électroniques discrets tels que des atténuateurs, des déphaseurs ou des circuits à retards temporels. Un réseau de formation de faisceaux numérique effectue quant à lui des traitements sur des signaux numériques. Le réseau 25 de formation de faisceaux peut être de tout type connu de l’homme du métier, y compris de type déphaseur (« phase shifter » en anglais) ou de type à retard réel (« true time delay » en anglais). Dans la suite de la description, on se place dans le cas où l’antenne 20 est une antenne en émission utilisant un réseau 25 de formation de faisceaux adaptatif analogique de type déphaseur. Tel qu’illustré sur la , le réseau 25 de formation de faisceaux reçoit en entrée une pluralité de signaux d’entrée comportant chacun un ou plusieurs canaux fréquentiels avec des niveaux de puissance respectifs. Dans l’exemple considéré et illustré sur la , on considère M signaux d’entrée. Chaque signal d’entrée correspond à un faisceau d’intérêt à former. Il peut y avoir plusieurs centaines, voire plusieurs milliers de faisceaux à former simultanément. La méthode selon l’invention est particulièrement bien adaptée au cas où le nombre de faisceaux à générer est important, par exemple au moins égal à cinq-cents (L ≥ 500), ou au moins égal à mille (L ≥ 1000). Rien n’empêcherait toutefois d’appliquer l’invention à des cas où le nombre de faisceaux serait moins important. On se place à titre d’exemple nullement limitatif dans le cas où il y a deux mille faisceaux à former (L = 2000). Ces signaux d’entrée correspondent par exemple à un plan de porteuses défini par un opérateur du satellite 30 de télécommunications. Chaque porteuse correspond à un canal fréquentiel occupant une sous-bande Δf d’une bande de fréquences ΔF. La bande de fréquences ΔF correspond par exemple à la bande Ka qui s’étend entre 27,5 et 31 GHz pour les communications en émission. La sous-bande Δf occupée par un canal fréquentiel présente par exemple une largeur de bande de 10 MHz. Chaque canal fréquentiel est associé à un niveau de puissance prédéterminé. Un faisceau peut être formé par une ou plusieurs porteuses (un ou plusieurs canaux fréquentiels par signal d’entrée). Il est à noter qu’un même canal fréquentiel peut être utilisé pour différents faisceaux (un même canal fréquentiel peut être présent dans différents signaux d’entrée). Les niveaux de puissance (puissance RF en conduit) des canaux fréquentiels dans les signaux d’entrée peuvent être déterminés à partir de spécifications de PIRE (acronyme de « Puissance Isotrope Rayonnée Equivalente ») fournies par l’opérateur. Ces spécifications fournissent le niveau minimal de PIRE que l’antenne doit fournir pour chaque canal fréquentiel. La PIRE est une combinaison de la directivité antenne et de la puissance conduite RF. La directivité antenne est connue via les diagrammes de rayonnement des sources 24 et les coefficients de pondération utilisés par le réseau 25 de formation de faisceaux (voir ci-après). On peut donc en déduire la puissance conduite RF nécessaire pour tenir les spécifications. De façon conventionnelle, le réseau 25 de formation de faisceaux est configuré pour appliquer des coefficients de pondération en phase et en amplitude aux différentes sources 24 et aux différents canaux fréquentiels afin de former les faisceaux radio souhaités pour desservir des zones géographiques d’intérêt à la surface de la Terre 50. Autrement dit, le réseau 25 de formation de faisceaux est configuré pour transmettre à chaque source 24 une combinaison linéaire des signaux d’entrée pondérés par les coefficients de pondération en phase et en amplitude afin de former les faisceaux radio souhaités. Les coefficients de pondération permettent de moduler la position et la forme de chacun des faisceaux à former. Dans ce but, le réseau 25 de formation de faisceaux peut comporter, de façon connue, des circuits duplicateurs, des déphaseurs et des atténuateurs permettant d’appliquer un coefficient de pondération en amplitude et en phase à chaque canal fréquentiel pour chaque source 24 du réseau 23. Il convient de noter qu’à un instant donné, toutes les sources 24 du réseau 23 ne sont pas nécessairement utilisées pour former les faisceaux souhaités. Si tel est le cas, il suffit de prendre en compte les sources 24 qui sont effectivement utilisées. Chaque source 24 est reliée à un amplificateur. Il peut s’agir par exemple d’un amplificateur de type SSPA ou un amplificateur à base de tubes. Dans l’exemple considéré, un amplificateur est associé à une et une seule source 24 (bijection entre les amplificateurs et les sources 24). Rien n’empêcherait toutefois de considérer le cas où un même amplificateur serait associé à plusieurs sources 24. Un amplificateur associé à une source 24 a pour rôle d’amplifier la puissance du signal transmis à la source 24 par le réseau 25 de formation de faisceau avant l’émission radio dudit signal par la source 24. Un amplificateur reçoit donc en entrée la combinaison linéaire des signaux d’entrée pondérés destinée à la source 24 auquel (ou aux sources 24 auxquelles) l’amplificateur est associé. Un amplificateur présente généralement des caractéristiques non linéaires lorsque le niveau de puissance attendue en sortie de l’amplificateur est élevé. Plus l’amplificateur doit fournir de la puissance en sortie, plus il est non-linéaire, et plus il va générer des produits d'intermodulation. Les produits d’intermodulation issus de la non-linéarité des amplificateurs impactent le bilan de liaison d’une communication établie entre le satellite 30 de télécommunications et une station terrestre de réception. Cet impact peut être évalué en calculant un rapport entre signal utile et bruit d’intermodulation (« rapport C/IM ») au niveau du signal rayonné au sol dans une zone géographique d’intérêt (zone géographique desservie par un faisceau radio émis par l’antenne 20 active multifaisceaux). Ce rapport C/IM rayonné peut notamment être calculé pour une porteuse donnée (c’est-à-dire pour un canal fréquentiel particulier). Les figures 4 et 5 représentent les principales étapes de deux modes de mise en œuvre d’une méthode 100 visant à contrôler la configuration de la charge utile d’un satellite 30 de télécommunications. La méthode 100 selon l’invention permet de calculer un rapport C/IM rayonné pour au moins une porteuse et pour au moins une zone géographique d’intérêt. La méthode 100 permet ensuite de vérifier un critère prédéterminé en fonction du rapport C/IM rayonné calculé. Dans l’exemple considéré, cette méthode est mise en œuvre par un ordinateur 10 du centre 17 de contrôle de mission. Autrement dit, dans l’exemple considéré, cette méthode est utilisée pendant une phase opérationnelle du satellite 30 de télécommunications. Il convient toutefois de noter que cette méthode pourrait également être utilisée pendant des phases amont, par exemple pendant une phase de conception ou de dimensionnement de la charge utile du satellite 30 de télécommunications. La représente schématiquement un ordinateur 10 configuré pour mettre en œuvre la méthode 100 selon l’invention pour valider la configuration de la charge utile d’un satellite 30 de télécommunications. Tel qu’illustré sur la , l’ordinateur 10 comporte au moins un processeur 11 et au moins une mémoire 12. La mémoire 12 mémorise des instructions de code d’un programme 13 informatique qui, lorsqu’elles sont exécutées par le processeur 11, configure le processeur 11 pour mettre en œuvre la méthode 100 selon l’invention. L’ordinateur 10 peut également comporter une interface utilisateur 14 (aussi connue sous le terme IHM, acronyme de « Interface Homme-Machine ») pour recevoir d’éventuels paramètres d’entrée. Alternativement ou en complément, l’ordinateur 10 peut également être connecté à une autre machine configurée pour fournir tout ou partie des paramètres d’entrée. Tel qu’illustré sur les figures 3 et 4, la méthode 100 prend en entrée différents paramètres. Un premier ensemble de paramètres correspond aux niveaux de puissance des canaux fréquentiels composant les signaux d’entrée. Cela correspond à un plan de fréquences alloué pour la mission affectée au satellite 30 de télécommunications. Ce plan de fréquences est par exemple fourni par un opérateur du système de télécommunications décrit en référence à la . Le plan de fréquence indique le niveau de puissance en conduit (« RF conducted power » en anglais) de chaque canal fréquentiel contenu dans les signaux d’entrée. Un deuxième ensemble de paramètres correspond aux coefficients de pondération associés aux différents canaux fréquentiels et aux différentes sources 24. Ces coefficients de pondération comprennent notamment des poids en amplitude appliqués par le réseau 25 de formation de faisceaux aux différents canaux fréquentiels pour chaque source 24. Ces poids correspondent à une configuration particulière du réseau 25 de formation de faisceaux de l’antenne 20 pour former des faisceaux radio permettant de répondre à la mission du satellite 30 de télécommunications à un instant donné. Dans la présente demande, l’expression « les coefficients de pondération associés à un amplificateur » doit être comprise comme « les coefficients de pondération associés à une source 24 (ou aux sources 24 auxquelles) l’amplificateur est associé ». Un troisième ensemble de paramètres correspond à des caractéristiques de linéarité des amplificateurs associés aux différentes sources 24. Les caractéristiques de linéarité d’un amplificateur sont par exemple définies par une courbe décrivant un rapport de puissance de bruit (NPR pour « Noise Power Ratio » en anglais) en fonction d’une puissance en sortie de l’amplificateur (point d’opération, ou OBO pour « output back-off » en anglais). Cette courbe NPR = f(OBO) est par exemple fournie par le constructeur de l’amplificateur. L’OBO correspond à une puissance de sortie relativement à point de saturation de l’amplificateur. Un quatrième ensemble de paramètres correspond aux diagrammes de rayonnement des différentes sources 24 de l’antenne 20. Chaque diagramme de rayonnement est une représentation de la distribution dans l’espace tridimensionnel d’une grandeur représentative de la puissance rayonnée par une source 24. Tel qu’illustré sur les figures 3 et 4, la méthode 100 comporte une étape de détermination 101 d'un « masque normalisé » de densité spectrale de puissance (DSP) de bruit d'intermodulation (IM) à partir des niveaux de puissance des canaux fréquentiels composant les signaux d’entrée. La densité spectrale de puissance du bruit engendré par les produits d’intermodulation des canaux fréquentiels contenu dans le plan de fréquences peut être calculée de façon analytique tel que décrit dans le document « Calculation of third order intermodulation spectrum », Y.H. LAU et al., ELECTRONICIS LETTERS, October 11th, 1990, Vol. 26, No. 21. Cette méthode permet de calculer la densité spectrale de puissance des produits d’intermodulation en convoluant le spectre en puissance du signal présent en entrée d’un amplificateur. Il est envisageable de considérer uniquement les produits d'intermodulation d'ordre trois (ce sont ceux qui dominent le spectre de puissance du bruit d’intermodulation). Les produits d'intermodulation d'ordre trois comportent les produits triples (du type , formés à partir de trois canaux fréquentiels distincts centrés respectivement sur les fréquences , et ), et les produits doubles (du type , formés à partir de deux canaux fréquentiels distincts centrés respectivement sur les fréquences , et ). La puissance d’un produit triple est 6 dB plus forte que celle d’un produit double. La densité spectrale de puissance du bruit engendré par les produits d’intermodulation d’ordre trois peut alors s’écrire sous la forme : Dans cette expression, l’opérateur correspond à une opération de convolution. N correspond au nombre de canaux fréquentiels différents contenus dans les signaux d’entrée. Le terme correspond au spectre de puissance d’un canal fréquentiel d’indice . On considère par exemple qu’il prend la forme d’un créneau centré à la fréquence et de largeur de bande Δf. L’amplitude du créneau correspond au niveau de puissance attribué au canal fréquentiel (différents niveaux de puissance sont attribués aux différents canaux fréquentiels, selon le plan de fréquences considéré). On a alors ) et . Le terme correspond au spectre de puissance du signal d’entrée : . Les convolutions impliquant le spectre de puissance du signal d’entrée ( ) peuvent être calculées en utilisant des routines FFT (acronyme anglais de « Fast Fourier Transform », en français « transformation de Fourier rapide »). Les convolutions impliquant des canaux fréquentiels de façon individuelle ( ) peuvent être obtenues de façon numérique ou analytique. Le « masque normalisé » de densité spectrale de bruit d’intermodulation peut alors être déterminé en normalisant le résultat par rapport à la puissance de bruit d'intermodulation totale sur la bande de fréquences ΔF. La puissance de bruit d'intermodulation totale est obtenue en intégrant la densité spectrale de puissance de bruit d'intermodulation sur la bande de fréquences ΔF. Selon un autre exemple, le résultat peut être normalisé par rapport à la puissance de bruit d'intermodulation agrégée sur la bande de fréquences ΔF pour des valeurs de puissance comprises dans une plage de valeurs prédéterminée (par exemple pour les valeurs de puissance comprises entre un seuil inférieur, par exemple -10 dB, et une valeur de puissance maximale). D’autres méthodes de normalisation de la densité spectrale de bruit d’intermodulation pourraient être envisagées (par exemple une normalisation par rapport à une valeur de puissance maximale plutôt que par rapport à une valeur de puissance agrégée). Le choix d’une méthode de normalisation n’est qu’une variante de l’invention. Il convient de noter qu’il est également envisageable de prendre en compte des produits d'intermodulation d'ordre impair supérieur à trois (par exemple les produits d’intermodulation d’ordre cinq ou sept). Ces produits d'intermodulation peuvent être calculés par une méthode analytique similaire à celle décrite ci-avant pour les produits d'intermodulation d’ordre trois. La prise en compte des produits d'intermodulation d'ordre impair supérieur à trois permet d’obtenir une meilleure précision pour la détermination 101 du masque normalisé de densité spectrale de bruit d’intermodulation. Toutefois, cela implique un plus grand nombre de calculs et par conséquent une augmentation du temps nécessaire pour exécuter la méthode 100 selon l’invention. Au lieu de calculer précisément les produits d'intermodulation d'ordre impair supérieur à trois, il est envisageable d’appliquer à la densité spectrale de bruit d’intermodulation une valeur de correction visant à compenser les produits d'intermodulation non pris en compte. Cette valeur de correction peut par exemple être estimée par simulation en comparant pour plusieurs exemples la densité spectrale de bruit d’intermodulation calculée en prenant en compte les produits d’intermodulation d’ordre trois et cinq (voire les produits d’intermodulations d’ordre trois, cinq et sept) avec la densité spectrale de bruit d’intermodulation calculée en prenant en compte uniquement les produits d’intermodulation d’ordre trois. La valeur de correction peut alors être déterminée en fonction d’une différence moyenne observée entre les densités spectrales ainsi calculées. Dans le mode de mise en œuvre décrit en référence à la , le masque normalisé de densité spectrale de puissance de bruit d'intermodulation est le même pour tous les amplificateurs. Autrement dit, pour calculer le masque normalisé de densité spectrale de puissance de bruit d'intermodulation d’un amplificateur, on considère l’ensemble de tous les canaux fréquentiels présents dans les signaux d’entrée, et on ne prend pas en compte les coefficients de pondération. Cela ne devrait théoriquement pas être le cas car les canaux fréquentiels en entrée de l’amplificateur et les niveaux de puissance desdits canaux fréquentiels ne sont pas les mêmes d’un amplificateur à l’autre (car les coefficients de pondération associés aux différents canaux fréquentiels ne sont pas les mêmes d’un amplificateur à l’autre). L’expérience montre toutefois que cette approximation est acceptable car l’erreur introduite dans le calcul du rapport C/IM avec cette hypothèse reste relativement faible. Cette approximation est d’autant plus valable que le nombre de sources 24 est grand (la puissance est moyennée sur l’ensemble des sources). Il est cependant également possible de ne pas faire cette approximation et de prendre en compte, pour chaque amplificateur, les coefficients de pondération associés aux différents canaux fréquentiels dans le calcul du masque normalisé de densité spectrale de puissance de bruit d’intermodulation. C’est ce qui est proposé par le mode de réalisation décrit en référence à la . Dans ce cas, un masque normalisé différent est calculé pour chaque amplificateur. Ce mode de réalisation permet d’augmenter la précision du rapport C/IM calculé, mais il nécessite un plus grand nombre de calculs et donc un temps d’exécution plus important. Les étapes suivantes de la méthode 100 selon l’invention (étapes 102 à 110) sont identiques pour les deux modes de réalisation décrits respectivement aux figures 4 et 5. La méthode 100 comporte une étape de détermination 102, pour chaque amplificateur, d'un point d'opération de l'amplificateur. Il est déterminé à partir des niveaux de puissance respectifs des canaux fréquentiels des signaux d'entrée, et à partir des coefficients de pondération associés respectivement audit amplificateur (c’est-à-dire aux coefficients de pondération associés à la source 24 à laquelle est associé ledit amplificateur). Le point d’opération correspond à la puissance agrégée en sortie de l’amplificateur sur l’ensemble de la bande de fréquences ΔF considérée. En sortie de l’étape 102, on connaît donc le point d’opération (c’est-à-dire la puissance en sortie) de chaque amplificateur. La méthode 100 comporte ensuite une étape de détermination 103, pour chaque amplificateur, d'une puissance totale de bruit d'intermodulation en fonction du point d'opération et en fonction des caractéristiques de linéarité de l'amplificateur. Comme indiqué précédemment, les caractéristiques de linéarité de l’amplificateur sont par exemple définies par une courbe NPR = f(OBO) fournie par le constructeur de l’amplificateur. Il est toutefois également envisageable d’obtenir les caractéristiques de linéarité de l’amplificateur par des méthodes de simulation. Cette simulation peut notamment être mise en œuvre en modélisant l’amplificateur avec un modèle de Shimbo. Les caractéristiques de linéarité de l’amplificateur peuvent alors être définies par une courbe de type C/IM = f(OBO). En sortie de l’étape 103, on connaît le niveau de puissance avec lequel chaque amplificateur contribue au bruit d’intermodulation, mais on ne connaît pas la répartition spectrale des produits d’intermodulation (les produits d’intermodulation ne sont pas répartis uniformément en termes de puissance). La répartition spectrale des produits d’intermodulation est en revanche indiquée par le masque normalisé déterminé à l’étape 101. La méthode 100 comporte une étape de détermination 104, pour chaque amplificateur, d'un « masque absolu » de densité spectrale de puissance de bruit d'intermodulation en pondérant le masque normalisé (déterminé à l’étape 101) avec la puissance totale de bruit d'intermodulation (déterminée à l’étape 103). Un « masque absolu » de densité spectrale de puissance de bruit d’intermodulation est ainsi déterminé pour chaque amplificateur. On obtient ainsi en sortie de l’étape 104, pour chaque amplificateur, un masque de densité spectrale de puissance de bruit d’intermodulation qui est représentatif à la fois de la répartition spectrale et des niveaux de puissance des produits d’intermodulation propres à l’amplificateur. La méthode 100 comporte ensuite une étape de détermination 105, pour chaque source 24, d'un « masque rayonné » de densité spectrale de puissance de bruit d'intermodulation en pondérant le masque absolu (déterminé à l’étape 104) de l'amplificateur auquel la source est reliée par le diagramme de rayonnement de la source 24. L’étape 105 permet ainsi d’ajouter la dimension spatiale, via le diagramme de rayonnement. Le « masque rayonné » permet ainsi de déterminer en un point donné au niveau de la surface terrestre, pour chaque source 24, une densité spectrale de puissance de bruit d'intermodulation rayonnée en ce point. La méthode 100 comporte ensuite une étape de détermination 106 d'un « masque rayonné total » de densité spectrale de puissance de bruit d'intermodulation par une sommation en puissance des masques rayonnés de densité spectrale de puissance de bruit d'intermodulation des différentes sources 24. Le masque rayonné total permet de déterminer en un point donné au niveau de la surface terrestre, pour l’ensemble des sources 24, une densité spectrale de puissance de bruit d'intermodulation rayonnée en ce point. L’utilisation de la somme en puissance des masques rayonnés repose sur une simplification qui revient à considérer que les bruits d’intermodulation entre les différentes sources sont décorrélés, ce qui n’est théoriquement pas le cas. Toutefois, l’expérience montre que cette simplification est acceptable car l’erreur introduite dans le calcul du rapport C/IM avec cette hypothèse reste relativement faible en comparaison avec des résultats fournis par modèles théoriquement exacts. Cette simplification permet de réduire les temps de calcul pour fournir une modélisation des performances en termes de C/IM rayonné en un temps limité. Il est à noter que plus le nombre de produits d’intermodulation est élevé, et plus les bruits d’intermodulation entre les différentes sources sont décorrélés. La méthode 100 selon l’invention est donc particulièrement bien adaptée à des antennes permettant de former un grand nombre de faisceaux simultanément. De façon optionnelle, une valeur de correction peut être appliquée au rapport C/IM calculé pour prendre en compte l’erreur introduite par l’approximation liée à la sommation en puissance des masques rayonnés. Cette valeur de correction peut notamment être déterminée en fonction d'un facteur de réutilisation des canaux fréquentiels pour former les différents faisceaux radio souhaités. Il peut par exemple s’agir d’un taux de réutilisation moyen (correspondant à la moyenne des taux de réutilisation des différents canaux fréquentiels). Il a en effet été observé que plus le taux de réutilisation moyen est important, et plus il convient d’appliquer une valeur de correction élevée au rapport C/IM calculé pour compenser l’approximation liée à la sommation en puissance des masques rayonnés. Selon un autre exemple, le facteur de réutilisation peut correspondre à une différence entre le taux de réutilisation maximal et le taux de réutilisation minimal parmi les différents canaux fréquentiels. Il a en effet été observé que plus cette différence est grande (c’est-à-dire plus le taux de réutilisation est hétérogène pour les canaux fréquentiels) et plus il convient d’appliquer une valeur de correction élevée au rapport C/IM calculé pour compenser l’approximation liée à la sommation en puissance des masques rayonnés. On considère par la suite un canal fréquentiel particulier et une zone géographique d’intérêt particulière (zone desservie par un faisceau), et on cherche à déterminer le rapport C/IM rayonné dans cette zone géographique d’intérêt pour le canal fréquentiel considéré. Par exemple, on détermine le rapport C/IM rayonné en un point correspondant au centre de la zone géographique. Alternativement, il est envisageable de déterminer un rapport C/IM rayonné en plusieurs points de la zone géographique d’intérêt (par exemple selon un quadrillage de la zone géographique d’intérêt). Dans ce but, la méthode 100 comporte une étape de détermination 107 d'une puissance rayonnée de bruit d'intermodulation pour la zone géographique d’intérêt et pour le canal fréquentiel considérés. Cette étape de détermination 107 est mise en œuvre par une intégration, sur la sous-bande de fréquences Δf correspondant audit canal fréquentiel, du masque rayonné total de densité spectrale de puissance de bruit d'intermodulation pour la zone géographique d'intérêt considérée. La méthode 100 comporte ensuite une étape de détermination 108, pour le canal fréquentiel considéré, d'une puissance isotrope équivalente rayonnée par l'antenne 20 active multifaisceaux (PIRE de la porteuse correspondant au canal fréquentiel). Cette PIRE peut être calculée à partir de la puissance RF en conduit de la porteuse et à partir des diagrammes de rayonnement et des coefficients de pondération. La méthode 100 comporte ensuite une étape de détermination 109, pour la zone géographique d’intérêt et pour le canal fréquentiel considérés, d'un rapport « C/IM rayonné » entre la puissance isotrope équivalente rayonnée (déterminée à l’étape 108) et la puissance rayonnée de bruit d'intermodulation (déterminée à l’étape 107). Enfin, la méthode 100 comporte une étape de vérification 110 d’un critère prédéterminé en fonction du rapport « C/IM rayonné » ainsi déterminé. Cette étape permet de valider la configuration de la charge utile du satellite 30 de télécommunications en termes de performances du C/IM rayonné. Il est par exemple envisageable de lever une alerte (par exemple via l’interface utilisateur 14) si le C/IM rayonné est inférieur à un seuil prédéfini. De telles dispositions permettent de vérifier s’il est nécessaire de prendre des mesures pour modifier la configuration de la charge utile du satellite 30 de télécommunications pour qu’il puisse accomplir sa mission de façon performante. Selon un autre exemple, il est possible de définir une fonction de répartition (CDF, pour « Cumulative Distribution Function » en anglais) d’un rapport C/IM rayonné sur l’ensemble de la bande de fréquences ΔF considérée et pour l’ensemble de la région de couverture considérée. Par exemple, pour un point de la région de couverture, il est possible de déterminer par quelle porteuse il est illuminé, et de calculer le rapport C/IM rayonné pour la porteuse en question au point considéré. Si le point considéré est illuminé par plusieurs porteuses, il est possible de calculer une valeur moyenne, ou une valeur correspondant à un pire cas, de rapport C/IM rayonné au point considéré. Une telle fonction de réparation permet par exemple de déterminer quel pourcentage de la région de couverture présente un rapport C/IM rayonné supérieur ou égal à un seuil prédéfini. La description ci-avant illustre clairement que, par ses différentes caractéristiques et leurs avantages, la présente invention atteint les objectifs fixés. En particulier, la solution proposée permet de calculer un rapport C/IM rayonné dans un temps limité afin notamment de permettre son utilisation lors de phases opérationnelles, et ce tout en assurant un niveau de précision suffisant. La complexité de l’évaluation du rapport C/IM est liée au nombre élevé de porteuses et de faisceaux nécessaires pour assurer la mission du satellite 30 de télécommunications. La méthode proposée est donc particulièrement bien adaptée pour évaluer le rapport C/IM dans les cas de charge utile complexe (les autres méthodes de calculs disponibles dans l’art antérieur étant alors trop longues à mettre en œuvre). La méthode proposée utilise en effet une approche globale et analytique qui permet de s’affranchir de nombreux calculs unitaires et ainsi conserver un temps d’exécution quasi négligeable. Méthode (100) mise en œuvre par un ordinateur (10) pour valider une configuration d’une charge utile d’un satellite (30) de télécommunications, ledit satellite (30) étant configuré pour émettre un signal avec une antenne (20) active multifaisceaux comportant un réseau (25) de formation de faisceaux, le réseau (25) de formation de faisceaux étant relié à une pluralité de sources (24) et adapté à former simultanément plusieurs faisceaux d'intérêt desservant des zones géographiques à la surface de la Terre, le réseau (25) de formation de faisceaux étant configuré pour recevoir en entrée, pour chaque faisceau d’intérêt à former, un signal d'entrée comportant un ou plusieurs canaux fréquentiels avec des niveaux de puissance respectifs, et pour transmettre à chaque source (24) une combinaison linéaire des signaux d’entrée pondérés par des coefficients de pondération, chaque source (24) étant reliée à un amplificateur, la méthode (100) comporte pour chaque amplificateur : une détermination (101) d'un « masque normalisé » de densité spectrale de puissance de bruit d'intermodulation à partir des canaux fréquentiels composant les signaux d'entrée et des niveaux de puissance respectifs desdits canaux fréquentiels, une détermination (102) d'un point d'opération correspondant à une puissance de sortie de l'amplificateur, le point d'opération étant déterminé en fonction des niveaux de puissance des canaux fréquentiels des signaux d'entrée, et en fonction des coefficients de pondération associés audit amplificateur, une détermination (103) d'une puissance totale de bruit d'intermodulation en fonction du point d'opération et de caractéristiques de linéarité de l'amplificateur, une détermination (104) d'un « masque absolu » de densité spectrale de puissance de bruit d'intermodulation par une pondération du masque normalisé avec la puissance totale de bruit d'intermodulation, la méthode (100) comporte pour chaque source (24) : une détermination (105) d'un « masque rayonné » de densité spectrale de puissance de bruit d'intermodulation à partir du masque absolu de densité spectrale de puissance de bruit d'intermodulation de l'amplificateur auquel la source (24) est reliée, et à partir d'un diagramme de rayonnement de la source, la méthode (100) comporte en outre : une détermination (106) d'un « masque rayonné total » de densité spectrale de puissance de bruit d'intermodulation par une sommation en puissance des masques rayonnés des différentes sources (24), pour au moins un canal fréquentiel d'intérêt, et pour au moins une zone géographique d'intérêt : une détermination (107) d'une puissance rayonnée de bruit d'intermodulation par une intégration, sur une sous-bande de fréquences correspondant audit canal fréquentiel, du masque rayonné total de densité spectrale de puissance de bruit d'intermodulation pour la zone géographique d'intérêt, une détermination (108), pour le canal fréquentiel d'intérêt, d'une puissance isotrope équivalente rayonnée par l'antenne (20) active multifaisceaux, une détermination (109) d'un rapport « C/IM rayonné » entre la puissance isotrope équivalente rayonnée et la puissance rayonnée de bruit d'intermodulation, une vérification (110) d’un critère prédéterminé en fonction du rapport « C/IM rayonné » ainsi déterminé. Méthode (100) selon la revendication 1 dans laquelle le masque normalisé de densité spectrale de puissance de bruit d'intermodulation est commun à tous les amplificateurs. Méthode (100) selon la revendication 1 dans laquelle le masque normalisé de densité spectrale de puissance de bruit d'intermodulation est déterminé pour un amplificateur en fonction des coefficients de pondération associés audit amplificateur. Méthode (100) selon l’une quelconque des revendications 1 à 3 dans laquelle la détermination (101) du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de produits d'intermodulation d'ordre trois générés par les canaux fréquentiels des signaux d'entrée. Méthode (100) selon la revendication 4 dans laquelle la détermination (101) du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte en outre un calcul de produits d'intermodulation d'ordre impair strictement supérieur à trois générés par les canaux fréquentiels des signaux d'entrée. Méthode (100) selon la revendication 4 dans laquelle la détermination (101) du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte une application d’une valeur de correction visant à compenser des produits d'intermodulation d'ordre strictement supérieur à trois non pris en compte. Méthode (100) selon l’une quelconque des revendications 1 à 6 dans laquelle la détermination (101) du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de normalisation en fonction du résultat d’une intégration, sur une bande de fréquences, de la densité spectrale de puissance de bruit d'intermodulation. Méthode (100) selon l’une quelconque des revendications 1 à 6 dans laquelle la détermination (101) du masque normalisé de densité spectrale de puissance de bruit d'intermodulation comporte un calcul de normalisation en fonction du résultat d’une intégration, sur une bande de fréquences, de la densité spectrale de puissance de bruit d'intermodulation pour des valeurs de puissance comprises dans une plage de valeurs prédéterminée. Méthode (100) selon l’une quelconque des revendications 1 à 8 comportant une application d’une valeur de correction au rapport « C/IM rayonné », la valeur de correction étant déterminée en fonction d'un facteur de réutilisation des canaux fréquentiels pour former les faisceaux d'intérêt. Méthode (100) selon l’une quelconque des revendications 8 ou 9 dans laquelle la bande de fréquences est l’une des bandes W, Q, V, Ka, Ku, L, S, C. Produit programme (13) d’ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur (10), conduisent celui-ci à mettre en œuvre les étapes de la méthode (100) selon l’une quelconque des revendications 1 à 10. Support d'enregistrement (12) lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes de la méthode (100) selon l’une quelconque des revendications 1 à 10.