FR 2667160 A1 19920327 FR 9011861 A 19900926 PROCEDE ET DISPOSITIF DE MESURE DE L'INTEGRITE D'UNE EMISSION L'invention se situe dans le domaine des procédés et dispositifs de mesure destinés à contrôler l'intégrité d'une émission en particulier l'intégrité de l'émission d'un signal émis par un radiophare. On sait qu'un radiophare est un dispositif radio émettant sur une fréquence porteuse en général V H F ou U H F. un ensemble de signaux Ces signaux sont émis de telle sorte qu'un mobile disposant d'un récepteur approprié calé sur la fréquence de la porteuse peut, à l'aide de ce signal, déterminer la position angulaire de l'axe issu de l'émetteur sur lequel il se trouve Les radiophares peuvent être omnidirectionnels ou directionnels Parmi les radiophares omnidirectionnels on trouve le radiophare connu internationalement sous le nom de V O R. (V.H F omnirange) qui permet à un avion équipé de façon ad hoc de déterminer l'angle de l'axe issu du radiophare sur lequel il se trouve par rapport à la direction du nord magnétique au lieu du phare Parmi les radiophares directionnels on trouve par exemple les dispositifs connus internationalement sous le nom d'I L S (Instrument landing system dispositif d'atterrissage aux instruments) Ces radiophares permettent de matérialiser dans l'espace un axe de descente vers une piste d'atterrissage Les signaux émis par les radiophares sont en général émis par un réseau d'antennes, chaque antenne du réseau étant alimentée par un signal H F modulé d'une façon spécifique. L'exactitude des informations de position que l'on peut tirer des radiophares est directement fonction de la capacité de ces équipements à émettre un signal restant à l'intérieur de normes déterminées Dans le cas de radiophares tels que le V O R. ou l'I L S ces normes sont définies par l'Organisation Internationale de l'Aviation Civile (O A C I) Ainsi par exemple le chapitre 3 1 du document l'O A C I, est consacré aux 1 'I L S. "& 3 1 Le système I L S coml ci-après: vertical, piste. dit "Annexe 10 ", édité par spécifications techniques de prend les éléments fondamentaux a) radiophare d'alignement de piste (dit Localizer) b) radiophare d'alignement de descente (dit Glide-path) c) radiobornes. Le radiophare d'alignement de piste définit un plan lieu des points de DDM nulle, et contenant l'axe de Le radiophare d'alignement de descente définit une surface, lieu des points de DDM nulle, dont l'intersection avec le plan du Localizer définit la trajectoire de guidage. Les radiobornes donnent des indications de position par rapport au seuil de piste. Le sigle DDM désigne la différence entre les profondeurs de modulation de la porteuse par les deux tensions BF Hz et 150 Hz La DDM constitue le signal de guidage exploité par l'aéronef. En dehors des performances proprement dites, les principales exigences portent sur les deux points suivants: continuité de service Qualité liée à la rareté des interruptions du signal rayonné. intégrité de l'I L S Qualité liée à la confiance que l'on peut avoir dans l'exactitude des renseignements fournis par l'intallation " Suivant les propres termes du document OACI, il est possible de dire que l'objectif essentiel d'une installation ILS est de réaliser et de maintenir: l l le plus haut degré possible d'intégrité, de "& 3 1 2 6. fiabilité, de stabilité de fonctionnement du système dans les conditions d'environnement les plus défavorables susceptibles d'être rencontrées Le Supplément C à la 1 ère partie contient des éléments indicatifs permettant de réaliser cet objectif l l 2 8 2 13 En général, la conception de l'équipement de contrôle repose sur le principe d'un contrôle continu des signaux électromagnétiques rayonnés à des points spécifiques à l'intérieur du volume de couverture, ce contrôle devant permettre de s'assurer qu'ils sont conformes aux normes énoncées en 3 1 3 11 et 3 1 5 7 ". Pour se conformer à ces directives le dispositif de contrôle de l'intégrité du signal du "glide path" comporte d'une part un capteur situé au voisinage du mât d'émission à une position représentative du signal émis sur l'axe de descente et d'autre part un dispositif constitué par des coupleurs mesurant les courants d'antennes et les recombinant de façon à recons- tituer le signal théoriquement reçu dans deux autres directions. En ce qui concerne le VOR le dispositif comporte au moins un capteur situé à proximité du réseau d'antennes. Les résultats des mesures effectuées à l'aide de ces dispositifs doivent rester à l'intérieur de fourchettes prédéterminées L'inconvénient principal de ce mode de contrôle de l'intégrité est signalé par l'OACI toujours dans l'annexe 10 à propos de l'ILS. "Bien que ce contrôle indique dans une certaine mesure que le signal électromagnétique partout ailleurs dans le volume de couverture se situe également dans les limites de tolérance, cette constatation est largement induite". Il est donc indispensable de procéder périodiquement à des vérifications minutieuses en vol et au sol afin d'assurer l'intégrité du signal électromagnétique dans tout le volume de couverture. En effet le champ au voisinage du récepteur de bord est perturbé par toutes les réflexions parasites en provenance de l'environnement proche ou lointain. Ceci signifie que la valeur du champ en un point donné dépend de la valeur des champs rayonnés dans toutes les autres directions: toute modification des diagrammes d'antennes est susceptible de modifier le signal de guidage au point d'utilisation. Il s'agit là du principal obstacle à la garantie de l'intégrité du signal dans l'espace, et de la principale insuffisance des dispositifs de contrôle actuels Ceux-ci ne mesurent en réalité que des échantillons représentatifs du champ émis dans deux ou trois directions prédéterminées Ces dispositifs, s'ils sont capables de prendre en compte globalement la plupart des paramètres, ne permettent ni d'évaluer avec précision leurs valeurs respectives, ni de déterminer le champ émis dans toutes les directions de l'espace. Le but de l'invention est de remédier à cet inconvénient en fournissant un procédé et un dispositif à l'aide desquels il est possible de s'assurer en permanence que le signal reçu en un point de l'espace à l'intérieur du volume de couverture est conforme à la norme définie par l'OACI dans le cas des radiophares ou par le constructeur d'un émetteur dans le cas le plus général Comme il n'est pas envisageable de disposer des capteurs dans tout l'espace le procédé selon l'invention repose essentiellement sur le calcul Pour être à même d'effectuer ces calculs, il faut tout d'abord disposer d'un modèle mathématique fiable représentant les paramètres caractéristiques du signal reçu en un point en fonction des coordonnées de ce point Il faut ensuite mesurer toutes les grandeurs nécessaires à la réalisation du calcul défini par le modèle Or dans l'état actuel des connaissances les modèles que l'on envisage font appel à des grandeurs que l'on ne sait pas mesurer de façon permanente sans interrompre l'émission De telles grandeurs sont par exemple les gains des antennes, variables en fonctions des conditions météorologiques, essentiellement la présence de givre sur les antennes ou les radômes, et dans certains cas les coefficients de réflexion du sol Si l'on savait mesurer en permanence la valeur de ces grandeurs, on serait mieux à même de calculer en fonction de la valeur de ces paramètres mesurés à tout instant et d'autres grandeurs que l'on sait déjà mesurer, la valeur du signal en tout point intéressant du volume de couverture L'emploi du procédé et du dispositif selon l'invention permet de s'assurer de la pertinence du modèle retenu, et donc de diminuer sensiblement les coûts d'exploitation tout en améliorant la qualité du service rendu L'invention permet aussi dans certains cas l'aide du diagnostic de la cause d'un disfonctionnement. A toutes ces fins l'invention a pour objet un procédé de contrôle permanent de l'intégrité d'une émission susceptible d'être affectée par des paramètres non directement accessibles à la mesure et dont on sait au moyen d'un modèle mathématique déterminer, l'influence sur le module et la phase du champ rayonné dans l'espace, l'émission étant réalisée à partir d'un réseau d'antennes, chaque antenne ou groupe d'antennes du réseau étant alimentée par une ou plusieurs fréquences porteuses H F, chacune de ces fréquences étant modulées d'une façon connue par un ou plusieurs signaux basse fréquence sinusoïdaux purs caractérisé en ce que: 1) on enregistre en au moins un point et pendant un temps T, la valeur de l'amplitude de la combinaison du signal reçu en ce point et d'un signal dit de référence dont la fréquence fondamentale est égale à l'une des fréquences H F du signal à analyser; 2) on extrait l'amplitude de composantes B F du signal ainsi obtenu et on les mémorise; 3) on recommence N fois (n désignant un entier) les étapes 1 et 2 en enregistrant le signal au même point et en déphasant à chaque fois le signal de référence d'une quantité a k 21 k étant un entier prenant successivement toutes les valeurs entre 1 et n 4) on traite les N amplitudes ainsi mémorisées comme les n valeurs d'une fonction F (ak) dont on extrait les valeurs de certains coefficients de sa transformée de Fourier. ) on recommence l'ensemble des opérations 1 à 4 au moins autant de fois qu'il y a de points de mesure; 6) on détermine les valeurs des paramètres inconnus par utilisation du modèle mathématique. Le point d'enregistrement pourra être soit à l'intérieur du radôme d'une antenne et dans ce cas le signal en ce point sera représentatif du signal émis par cette antenne, soit à l'extérieur et dans ce cas le signal en ce point sera représentatif du signal résultant de l'ensemble du réseau Dans ce dernier cas le point sera situé de préférence à proximité du réseau. Le temps T sera de préférence un nombre entier de périodes du signal B F à analyser. Les valeurs des coefficients de la transformée de Fourier des fonctions de k peuvent être en particulier la composante continue ou la phase et l'amplitude de la composante à la périodicité 21 E. Les principes théoriques ayant conduit à l'utilisation du procédé ci-dessus décrit seront explicités ci-après à l'aide d'un exemple fictif choisi simple à dessein. Considérons deux antennes A et B séparées d'une distance 1 L'antenne A est alimentée par un signal de la forme a cos (et + q 1) cos Q 1 t et l'antenne B par un signal de la forme b cos (tt + q 2) cos Q 22 t Dans cet exemple a et b représentent les modules des signaux, X la pulsation de la porteuse, Q 1 et 2 les pulsations de modulation, q 1 et e 2 la phase de la porteuse sur chacune des antennes A et B à l'instant t = O. On admet qu'un capteur placé en un point quelconque de l'espace recevra un signal en provenance de ces deux antennes qui est une fonction connue d'éléments connus tels la distance du capteur à l'antenne et le coefficient d'absorbtion de l'atmos- phère La distance au capteur est considérée comme un invariant et l'absorbtion comme affaiblissant de la même façon les signaux émis par chacune des antennes On admet par contre que les coefficients de transmission des antennes T A et TB seront variables de façon indépendante en fonction des conditions météorologiques (présence d'eau de neige ou de givre sur les radomes d'antenne). Le modèle que nous avons choisit conduit donc à admettre qu'en un point de l'espace il sera reçu un signal de la forme TA a cos (wt + ça) cos 21 t + TB b cos (wt + eb) cos Q 2 t Nous sommes donc bien dans les conditions requises pour appliquer le procédé Nous avons un signal dont le module et la phase en un point sont susceptibles d'être affectés par deux paramètres TA et TB que l'on ne peut mesurer en permanence au cours de l'émission On sait modéliser l'influence de ce paramètre sur les composantes du signal reçu. La première étape du procédé va être d'enregistrer en un point situé à proximité du réseau d'antennes un signal qui est la somme du signal transmis en ce point par le réseau d'antennes et d'un signal dit de référence à la même fréquence que la porteuse Lorsque l'on parle d'un point situé à proximité du réseau d'antennes, il s'agit d'un point qui est suffisamment proche pour que le signal capté en ce point puisse être ramené dans les installations concourant à l'émission (salle de contrôle par exemple) au moyen de guides ou de câbles coaxiaux tout en gardant une valeur suffisante pour les besoins de la mesure Dans l'exemple fictif ici décrit nous admettons que le signal de référence est de la forme: E = E O cos (wt + oek) Dans cette expression W est la fréquence de la porteuse commune, (o la phase pour t = O et k un déphasage qui prendra successivement toutes les valeurs comprises entre et 21 par n incrément de L'addition du signal de référence et du signal n de réseau suppose que le signal de référence remplit certaines conditions: être suffisant pour introduire une variation significative du niveau moyen; n'être ni trop faible, ni trop fort pour que le signal combiné ne soit jamais en état de surmodulation. Ces conditions, sans être trop contraignantes, demandent cependant une séquence d'étalonnage destinée à déterminer la valeur optimale à mettre sur un atténuateur programmable en fonction du signal à analyser et éventuellement de la valeur de la phase. Dans ces conditions le signal somme s'écrit E + E E =a' cos (wt+) cos t o A B a 1 + b' cos (wt + eb) cos Q 2 t + E cos (wt + ço + ik) avec a' = If A I 1, b' = ITB b 1 Pa = 1 + arg TA b = 92 + arg TB On sait que si z 1, z 2, z 3 désigne trois nombres complexes la somme z 1 z 2 z 3 aura pour module Il z 1 + z 2 + z 3 tel que: i 1 z 1 + Izlc 2 S + z 212 = 12 1 2 + 1211 2 + 1 Z 311 2 + l 1 + 2 +Z 12 12 1 1 2 1 1 il 1 lz 2 lc cos z 1 z 2 + 2 1 lzll l i 1 cos z 1 z 3 + A 2 lz 211 lz 31 cos z 2 z 3 Après détection l'expression du module du signal I s'écrit: 11 2 = a'2 CO 52 1 t +b'2 cos 2 Q t + E 2 l 2 O + 2 a'b' cos Ol t cos O 2 t cos (ça qb) + 2 a'E cos Ql t cos (qo + k a) + 2 b'E O cos O 2 t cos (qo + k b) Le signal somme IIl ci-dessus est le signal de combinaison décrit à l'étape 1 du procédé Pour réaliser l'étape 2 il faut encore filtrer le signal autour de fréquences BF qui apparaissent dans l'expression mathématique du signal de combinaison Le choix de filtrage à retenir est limité par la forme de l'expression mathématique du signal obtenu après détection Dans l'exemple simple que nous avons choisi les termes intéresssants pour l'emploi du procédé sont les termes qui sont des fonctions de Yk Ceci est également vrai dans le cas général. La forme mathématique de I N? comporte en outre des termes en 2 2 cos 2 Olt cos 2 Q 2 t qui peuvent s'analyser comme des signaux variant selon des pulsations 201 ou 2 02, des termes en cos Ol t cos Q 2 t qui peuvent s'analyser comme une somme de signaux variant seront des pulsations (Ol + O 2) et (O 1 O 2) et enfin des termes en cos l t et cos O 2 t Il est donc possible d'effectuer des filtrages centrés autour des pulsations 0, 2 Q 1 P 2 ' M 1 + 2 ' 2 1 2), 1 i 2 Nous voyons donc sur cet exemple simple que les choix possibles sont des combinaisons linéaires simples de O 1 et O 2. Nous aurons l'occasion de le vérifier à nouveau sur d'autres exemples Parmi les choix possibles on choisira de préférence ceux qui permettent d'accéder aux paramètres inconnus et ce, de façon simple Dans l'exemple fictif que nous avons choisi si par exemple nous filtrons sur 2 nous obtenons en sortie de filtre un signal SA de la forme: SA a a'E cos O 1 t cos (q, + k A) Si nous enregistrons ce signal pendant un nombre entier de périodes nous serons capable d'extraire de cet enregistrement l'amplitude de SA qui vaut dans ce cas: |SA| îe a'E cos (e + k - A Dans ce cas le temps T dont il est question pour la deuxième étape du procédé est égal à un nombre entier de périodes du signal à la pulsation Q 1 Dans l'exemple choisi pour obtenir TB il faudra également filtrer sur Q 2 pendant un temps T correspondant à un nombre entier de périodes du signal à la pulsation 02 * Dans le cas le plus général o on filtre sur une combinaison linéaire des différentes pulsations de modulation, le temps T sera un multiple entier d'un temps T égal au plus petit commun multiple PPCM des différentes périodes réelles intervenant dans la combinaison retenue Nous disons réelles pour prendre en compte le fait que la période correspondant par exemple au signal à la pulsation 201 est la moitié de celui correspondant à la période du signal à la pulsation O 1 De même s'il apparaît dans la combinaison un terme à la pulsation O 1 + O 2 ou O 1 22 la période réelle à prendre en compte pour ce terme est le PPCM des périodes correspondant à O 2 et O 2 Si l'on filtre en parallèle il sur plusieurs combinaisons le temps TM sera au moins égal au plus long des temps déterminés comme indiqué ci-dessus, mais il faudra en sortie de chacun des filtres disposer des séquenceurs interrompant l'enregistrement au bout d'un temps T multiple entier du temps Tm correspondant à la période de la combinaison filtrée Si l'on filtre en série chacun des temps T successifs sera égal à un multiple entier de la période du signal filtré et le séquencement des enregistrements devra tenir compte de ces différents temps Filtrer sur un nombre entier de périodes chacun des signaux est nécessaire pour avoir la meilleure précision possible sur la valeur de l'amplitude du signal filtré Ce point sera expliqué plus en détail plus loin à propos d'un exemple particulier de réalisation Nous avons vu plus haut que le signal en sortie du filtre centré sur la fréquence correspondant à la pulsation Q 1 était: SA a a'E cos cy + 'PA) Cos 2 t. Le module de ce signal est égal à |SA| a a'E cos (ç o + fa). Cette extraction termine la deuxième étape du procédé ci-avant décrit. Si maintenant nous recommencons la mesure N fois en incrémentant à chaque fois i d'une valeur 2 réalisant ainsi la troisième étape du procédé nous sommes en mesure de repré- senter SA en fonction de Yk Ces variations sont représentées A 211 figure la lorsque ak varie entre O et 211 par pas de n N étant dans l'exemple choisi égal à 16. On notera que dans l'exemple simple choisi, la différence entre le maximum et le minimum de cette fonction est égale à 2 a'E et que la phase de cette fonction est égale à PO ea' ço ça L'extraction de ces deux grandeurs constitue la quatrième étape du procédé ci-avant décrit. Ces deux grandeurs sont la phase et l'amplitude du coefficient pour la périodicité 211 de la transformée de Fourier de la fonction de Yk. Si maintenant nous faisons un filtrage sur une fréquence correspondant à la pulsation 22 la fonction obtenue de la même manière lorsque ik varie entre O et 2 H est représentée figure lb La différence entre le maximum et le minimum de cette fonction est égale à 2 b'E et la phase de cette fonction est égale à V o b La connaissance de chacune des grandeurs ISA| et |SB| en fonction de ak permet donc d'accéder en permanence au déphasage et à l'affaiblissement apportés par chacune des antennes A et B que nous avons appelé les gains Il suffit de reprendre notre modèle mathématique et d'écrire Sa a' = Ea = 1 ITA al o Sb b' = = IIT Bbi o q O -a = > l aarg T (Po e b = eo e 2 -arg T b pour réaliser la dernière étape du procédé On note que toutes les phases et tensions sont mesurées par rapport à l'amplitude et à la phase du signal de référence Les valeurs de a, b, el' P 2 seront obtenues en analysant par le même processus les signaux issus de capteurs au niveau de chacune des antennes A et B. Le procédé est dit permanent en ce sens que d'une part il n'est pas nécessaire d'interrompre l'émission pour effectuer la mesure et que d'autre part la mesure peut être renouvelée en permanence par répétition du procédé. L'invention est également relative à un dispositif pour effectuer la mesure de paramètres d'une émission non directement accessibles à la mesure dont on sait au moyen d'un modèle mathématique déterminer l'influence sur le module et la phase de l'émission en un point, l'émission étant réalisée à partir d'un réseau d'antennes, chaque antenne ou groupe d'antennes du réseau étant alimentée au moyen de un ou plusieurs émetteurs par une fréquence porteuse HF commune, cette fréquence étant pour chacune des antennes modulée par un ou plusieurs signaux sinusoïdaux B F. purs caractérisé en ce que le dispositif comprend des coupleurs recevant les courants d'alimentation des antennes; des capteurs placés à proximité du réseau d'antenne un sommateur comportant deux entrées et une sortie et recevant sur l'une de ses entrées par l'intermédiaire d'un commutateur programmable HF permettant le raccordement successif de chacun des capteurs et des coupleurs, l'un des signaux en provenance d'un capteur et sur l'autre entrée par l'intermédiaire d'un phaseur programmable un signal dit de référence capté par l'intermédiaire d'un coupleur en sortie de l'un des émetteurs concourant à l'émission, la sortie du sommateur étant raccordée à un détecteur; un calculateur recevant par l'intermédiaire d'une carte d'acquisition le signal en sortie du détecteur et envoyant par l'intermédiaire d'une carte interface des signaux de commande en direction du phaseur et en direction du commutateur haute fréquence. Ces deux cartes ne sont pas toujours nécessaires Leur présence est fonction du calculateur employé et des niveaux nécessaires pour les entrées et sorties La présence des coupleurs peut ne pas être nécessaire si le modèle mathématique n'exige pas la connaissance des courants d'antenne Le dispositif ci-dessus décrit qui est apte à réaliser le procédé selon l'invention ne comporte pas de filtres, car le filtrage est réalisé de façon numérique par le calculateur Il serait évidemment possible de placer en sortie du détecteur une série de filtres, chacun des filtres étant centré sur une fréquence qui est une combinaison linéaire de fréquences de modulation, un commutateur programmable dirigeant le signal en sortie de filtre vers la carte d'acquisition. Dans certains cas une solution mixte avec filtrage numérique par exemple pour les basses fréquences faibles et filtrage analogique pour les basses fréquences élevées pourra être un bon compromis si la fréquence d'acquisition des échantillons permise par la carte d'acquisition n'est pas assez élevée. De même il a été dit dans la description du procédé qu'il convenait que le signal de référence et le signal en provenance du capteur soient à des niveaux comparables Si ce n'est pas le cas on peut placer un atténuateur programmable entre la sortie du commutateur de capteurs et l'entrée du sommateur, le choix de l'atténuation est déterminé par étalonnage préalable. L'atténuateur est commandé par le calculateur par l'intermédiaire de la carte interface en fonction du capteur commuté. Le procédé et le dispositif selon l'invention seront maintenant illustrés en référence aux dessins annexés et à deux exemples de réalisation Le premier exemple est relatif à la mesure de paramètres susceptible d'affecter le signal émis par un glide d'I L S Le second qui sera décrit plus succintement est relatif à la mesure de paramètres susceptibles d'affecter l'émission d'un V O R Doppler. Dans les dessins annexés la figure 2 représente la forme générale du dispositif selon l'invention; la figure 3 représente une variante du dispositif général selon l'invention (avec commutateur et filtre analogique) la figure 4 représente une autre variante du dispositif général selon l'invention (avec atténuateur et interrupteur à diode); la figure 5 représente une autre variante du dispositif selon l'invention (avec plusieurs coupleurs d'introduction du signal de référence; la figure 6 représente une autre variante du dispositif selon l'invention (avec un mélangeur en lieu et place d'un sommateur et d'un détecteur). Sur la figure 2 représentant le dispositif selon l'invention on voit: des coupleurs 1 recevant des courants d'alimentation des antennes en sortie des émetteurs des capteurs 2 placés à proximité du réseau d'antenne un sommateur 3 comportant deux entrées 4, 5 et une sortie 6 et recevant sur l'une de ses entrées 5, par l'intermédiaire d'un commutateur programmable H F 7 permettant le raccordement successif de chacun des capteurs 2 et des coupleurs 1, le signal en provenance du capteur ou coupleur sélectionné, et sur l'autre entrée 4, par l'intermédiaire d'un phaseur programmable 8 un signal dit de référence, capté par l'intermédiaire d'un coupleur 9 en sortie de l'un des émetteurs concourant à l'émission, la sortie 6 du sommateur étant raccordée à un détecteur 10; un calculateur 11 recevant par l'intermédiaire d'une carte d'acquisition 12 le signal en sortie 16 du détecteur 10 et envoyant par l'intermédiaire d'une carte interface 13 des signaux de commande en direction du phaseur et en direction du commu- tateur haute fréquence 7 Les bus de transmission 17 des bits de commandes issus du calculateur sont représentés en grisé De préférence le phaseur programmable 8 est un phaseur à 2 P positions. La figure 3 représente une variante du dispositif selon l'invention dans laquelle le filtrage en sortie du détecteur 10 est réalisé par une série de filtres analogiques 14, chaque filtre étant centré sur une fréquence qui dans le cas le plus général est une combinaison linéaire de différentes fréquences de modulation de la porteuse H F Dans ce cas la sortie 16 du détecteur est raccordée à la série de filtre 14 et un commutateur programmable 15, commandé par le calculateur 11 par l'intermé- diaire de la carte interface d'adaptation de commande 13, sélectionne la sortie de filtre enregistrée Cette fonction commutation peut être incluse dans la carte d'acquisition. La figure 4 représente une autre variante de réalisation dans laquelle la sortie du commutateur 7 de signal d'entrée alimente un atténuateur programmable 18 dont la sortie est reliée à l'entrée 5 du sommateur 3 Un interrupteur à diode 19 commandé par le calculateur branché entre la sortie du phaseur 8 et l'entrée 4 du sommateur 3 permet d'interrompre l'arrivée du signal de référence Le commutateur comporte une position permettant de sélectionner une charge 20. Le dispositif de la figure 5 diffère des dispositifs précédents en ce que le signal de référence peut être choisi par le calculateur 11 Il comporte à cet effet plusieurs coupleurs 91 à 93 (dans le cas de la figure 5) La sortie de chacun de ces coupleurs peut être sélectionnée par l'intermédiaire d'un commutateur HF programmable 22 Un dispositif de synchronisation 21 permet l'envoi de signaux de synchronisation vers la carte d'acquisition 12 La sortie du détecteur 10 alimente la carte d'acquisition 12 d'une part directement et d'autre part par l'intermédiaire de filtres 14-1 14-2 dont la sortie est commandée par un commutateur programmable 15. Le dispositif de la figure 6 représente le dispositif selon l'invention dans lequel le sommateur a été remplacé par un mélangeur 23 recevant le signal de référence et le signal en provenance du commutateur 7 La sortie du mélangeur alimente un circuit d'adaptation 24 La sortie du circuit 24 est raccordée à la carte d'acquisition 12. Les montages représentés sur les figures 1 à 6 ne sont pas limitatifs de l'invention en particulier il est précisé que le nombre de coupleurs 9 peut être aussi grand que le nombre d'antennes du réseau Un seul capteur 2, sans coupleur 1, peut permettre de réaliser le dispositif selon l'invention. Le dispositif de la figure 4 a été utilisé pour contrôler l'intégrité de l'émission de l'émetteur glide path d'un ILS en utilisant le procédé selon l'invention Cet exemple de réalisation sera maintenant décrit de façon détaillée. On rappelle que le signal émis par un glide est une composition de trois signaux. Un signal dit CSB de la forme: CSB = Pl l 1 + m (sin Q 1 t + sin 22 t)l cos (wt + 1) Un signal dit SBO de la forme SBO = B (sin Q 1 t sin m 2 t) cos (wt + g 2) Un signal dit "'clearance" de la forme: C = C ( 1 + M 1 sin O 1 t + M 2 sin 02 t) cos l(w + A)t + 3) Dans ces expressions w, q 1 ' q 2 ' q 3 ' sont respectivement la pulsation de la porteuse et les déphasages de la porteuse à l'instant t = O. Q 1 et O 2 sont deux pulsations de modulation, Q 1 = 21 x 90 et 02 = 2 x 150. P Bl C sont les amplitudes des signaux. Le signal de guidage sur l'axe de descente est appelé DDM et sa valeur en un point a pour valeur: B 1 DDM = P cos (q 2 1) L'expression indique que le signal DDM est le produit scalaire des vecteurs représentant les champs CSB et SBO en ce point. Le signal de clearance est sur une porteuse séparée de la porteuse DDM d'un écart A égal à 18 K Hz Dans ces conditions on peut considérer en première approximation que ce signal est séparé par filtrage des signaux CSB et SBO. Le modèle mathématique retenu est constitué de la façon suivante: les deux vecteurs CSB et SBO sont des fonctions complexes des courants d'alimentation des antennes du réseau d'émission, ainsi que des éléments extérieurs pouvant influer sur les champs électromagnétiques au niveau du récepteur de bord. En particulier, le rayonnement réfléchi par le sol au voisinage immédiat des antennes joue un rôle essentiel On peut le caractériser par un coefficient de réflexion complexe p multipliant le champ rayonné par une antenne fictive, "image" de l'antenne réelle par le plan du sol. Si U représente le trajet direct entre l'antenne et le récepteur et v le trajet en provenance de l'antenne "image", le facteur global affectant le rayonnement d'une antenne s'obtient en multipliant son courant d'alimentation par un terme E = U + p V Si Bi représente le courant d'alimentation SBO de la ième antenne, Pl le courant CSB et Ti un coefficient représenta- tif de son état éventuel: présence de givre, d'eau etc on obtient pour les champs totaux: CSB = E Ti Pl Ei SBO = E Ti Bi Ei 1, représentant les champs rayonnés au voisinage immédiat du réseau d'antennes, c'est-à-dire sans tenir compte d'éventuelles réflexions de l'émission sur des reliefs éloignés. La somme porte sur toutes les antennes. Ces expressions permettent de mieux évaluer l'influence des différents facteurs sur le signal dans l'espace. Ces facteurs peuvent être classés en deux groupes principaux: facteurs internes * courants d'alimentation CSB, SBO, et Clearance facteurs externes: * position géométrique des antennes réelles et images * coefficient de réflexion du sol (caractéristique de la géométrie et de la nature du sol sur une distance d'environ 500 mètres devant le réseau) * coefficient de transmission des antennes. Dans le cas de l'émetteur glide de l'ILS le réseau est constitué de trois antennes placées sur un même mât vertical Les courants d'antenne sont captés par trois coupleurs 11 12 13 placés au niveau des sources antennes les paramètres variables que l'on doit mesurer pour faire fonctionner le modèle mathéma- tique sont donc les coefficients de transmission des antennes Ti qui sont donc au nombre de trois et le coefficient de réflexion du sol p Il a été admis que ce coefficient était le même pour chacune des trois antennes Il y a donc au total quatre paramètres à mesurer Pour la réalisation particulière à 1 'ILS il a été disposé quatre capteurs externes 21 à 24 Ces quatre capteurs ont été placés verticalement les uns au-dessus des autres au niveau d'un mât déjà existant sur lequel est monté le capteur faisant partie de l'installation de "monitoring" de l'ILS telle qu'elle est définie par l'OACI le capteur existant fait partie du groupe des quatre capteurs 2 à 24. Ce mât est situé à 85 m des antennes Le signal est ramené au niveau de l'installation de mesure par quatre câbles coaxiaux. Le signal de référence qui a été choisi et qui est introduit par l'intermédiaire du coupleur 9 au niveau du phaseur programmable 8 est le signal CSB Il a pour expression: Pref = Po cos lwt + a(k) + gol ( 1 + m (sin 1 t + sin 2 t)l Le signal reçu au niveau de la sortie 16 du détecteur 10 est donc la somme de ce signal et du signal en provenance de l'un des quatre capteurs. Le résultat de la détection va s'écrire comme étant la racine carrée de la somme des carrés des modules plus les doubles produits des modules et des phases respectives de chaque terme (Formule générale D = V 1 2 + V 2 + 2 V 1 V 2 cos) +V 2 12 On obtient: DP= lPo( 1 + msin 1 t + msin Q 2 t)l lP( 1 + msin Qt + msin O 2 t)l 2 + lB 1 (sin 21 t sin 22 t)l 2 + 2 P 1 Po(l+ msin 2 lt+msin Q 2 t) ( 1 +msin Q 1 t+msin Q 2 t) cos l 1-"(k)-Wol + 2 Po Bl(l+msin 2 Qlt+msin Q 2 t) (sinr Ql t-sinl 2 t) coslp 2-a(k)-ep Ol + 2 P 1 Bl(l+msin Qlt+msin Q 2 t) (sinl 1 t-sin Q 2 t) cos (e 1-P 2) En développant cette expression et en ne conservant que les termes en Q 1 t et O 2 t, les autres termes étant éliminés par le filtrage sur Q 1 puis sur Q 2 on obtient: 2 r 2 2 D 2 = lP 2 m + Pl 2 m + 4 m P 1 P cos l-1 (k) l + 2 P Bl cos lq 2-"(k) gol + 2 P 1 B 1 cos ( 1 2)l sin 1 t z l 2 + lPo 2 2 m + P 12 2 m + 4 m Pl Po cos le 1 a(k) eol - 2 Po B 1 cos lq 2a(k) aol 2 P 1 B 1 cos ( 1-2) l sin 2 t+ Le filtrage va nous donner les modules des deux vecteurs V 90 et V 150 correspondant à Qit et 02 t. Les expressions de la somme et de la différence sont donc: 2 2 V 90 + V 150 = 4 m P 2 + 4 m Pl + 8 m P 1 Po cos (e 1 ak go) = Sk V 90 V 150 = 4 Po B 1 cos ( 2ak go) + 4 P B cos (O 1 = D i O 1 111 D 2 k Avant de décrire la façon dont est effectué le filtrage numérique il convient d'expliciter comment est faite l'acquisition des données. L'acquisition consiste à prendre un certain nombre, dans le cas présent, 540 échantillons à 0,185 ms d'intervalle du signal issu du détecteur Les opérations de conversion analo- gique/digital sont prises en charge par la carte d'acquisition 12 pilotée directement par un programme du calculateur 11. L'échantillonnage couvre 3 périodes complètes du signal ILS composite soit 100 ms (la fréquence du signal 90 + 150 Hz est de Hz). Cette opération est renouvelée pour les 16 positions du phaseur. On obtient ainsi une mesure complète concernant l'un des capteurs à analyser. Les conditions de la mesure (valeur de l'atténuation, présence ou non du signal de référence) sont fixées par logiciel. Le filtrage consiste à extraire à partir de ces 540 échantillons le module du signal correspondant à la pulsation Q 1 = 2 ir x 90 et celui correspondant à la pulsation 02 = 2 x x 150. Considérons un tableau de 540 nombres complexes de module 1 et d'argument 6 = 0,185 x 90 x 360 x P dans lequel 0,185 est la période d'échantillonnage 90 est la fréquence du signal que l'on cherche à filtrer; P est le rang de l'élément dans le tableau. 6 exprime donc en degrés le pas d'échantillonnage par rapport au 90 Hz multiplié par le rang de l'échantillon. Chaque élément peut s'exprimer par Il faut également tenir compte du fait que le signal réellement analysé est limité dans le temps De façon classique on utilise à cette fin une fenêtre dite "fenêtre de pondération". La fenêtre utilisée ici est une fenêtre de Blackman. Le calcul réellement effectué prend en compte la multiplication du signal par la fenêtre de pondération. La moyenne sur une période du produit fenêtre x signal x e représente le coefficient de la composante à 90 Hz recherché. Le coefficient de la composante à 150 Hz est obtenu de façon analogue. Connaissant chacun de ces modules on peut calculer le module de la somme soit S et le module de la différence soit D. On recommence N fois l'opération en augmentant à chaque fois de 2 le retard ak introduit par le phaseur 8 Dans le cas de la réalisation N est égal à 16. Lorsque "k varie suivant les 16 valeurs du phaseur on obtient une variation périodique du module de S et D Si l'on associe à chacune des 16 valeurs de modules trouvées respectivement pour S et D l'argument oc k correspondant à chacune des 16 valeurs du phaseur, on définit deux tableaux de 16 nombres complexes représentant chacun une fonction périodique de oe k. Le terme de rang 2 de la transformée de Fourier rapide (FFT) de chacune des deux fonctions ainsi définies sera conformément à l'expression mathématique de S et D définie plus haut un nombre complexe Dans le cas de S, il aura pour module 8 m P 1 PO et dans le cas de D 4 P B 1 Les arguments seront égaux respectivement à (p 1 PO) et (V 2 eo) Si l'on suppose PO et m constant il est possible d'obtenir à un coefficient multiplicateur près P et B w O L'opération ci- dessus décrite sera répétée pour chacun des quatre capteurs 21 à 24 On dispose ainsi en quatre endroits différents de la valeur des vecteurs Pl et B 1. Pour atteindre les valeurs des coefficients p T 1 T 2 et T 3 il nous faut maintenant disposer de la valeur de ces mêmes grandeurs au niveau des antennes elles-mêmes. Pour cela l'opération d'acquisition sera recommencée au niveau de chacun des coupleurs d'antennes pour obtenir P et B Si maintenant nous connaissons les facteurs de transmission (déphasage et affaiblissements) introduits par les trajets U et V nous serons à même d'écrire les équations Pl Tî Cl Pl MI, + p Vil) 1 Ti C 1 Bl (Uil +P Vi) pour le capteur 21 et les mêmes équations pour chacun des trois autres capteurs. Dans chacune de ces équations le coefficient Ti représente le coefficient de transmission de l'antenne i Les facteurs C sont des coefficients propre à chacun des capteurs soit C 1 C 2 C 3 C 4 P Si nous faisons le rapport B nous obtenons l'équation 1 i i i il + il B Ti Bi (u + P V Nous disposons alors d'un système de quatre équations comportant les quatre inconnues p Tl T 2 T 3. Dans ces équations E et V ne dépendent que de la distance de chacune des antennes réelles ou image à chacun des capteurs. Les positions des antennes et capteurs sont entrées dans une mémoire du calculateur en coordonnées cartésiennes x, y, z dans un tableau nommé DONNEES GLD. Le point O est défini par l'intersection du plan de sol; du plan vertical perpendiculaire à l'axe de piste et passant par le centre de la première antenne; du plan vertical parallèle à l'axe de piste et passant par le centre de la première antenne. Les positions des antennes images sont introduites par le calcul en tenant compte de la pente du sol par rapport à l'horizontale qui doit être mesurée Les antennes images sont considérées comme symétriques des antennes réelles par rapport au sol réel. On calcule ensuite la distance entre chacune des antennes réelles ou image et chacun des capteurs soit au total pour 3 antennes réelles A 1 A 2 A 3, 3 antennes images A 4 A 5 A 6 et 4 capteurs 21 22 2 24 24 distances Y i désigne l'indice de l'antenne et j l'indice du capteur Si l'on désigne par ij l'angle entre l'axe de l'antenne i et l'axe du capteur j l'expression du champ reçu au niveau du capteur j en provenance de l'antenne i est: cos e E =E 1 ij Ei Yij L'expression du déphasage antenne capteur introduite par le temps de parcours de la distance Yij est: 2 TN ij = X ij Le mode de résolution des quatre équations sera maintenant abordé succintement. Le facteur de transmission de l'antenne est écrit sous la forme: Ti ='i ( 1 + (i) nominales Tn est le facteur de transmission mesuré en conditions i de fonctionnement, (i en est la variation recherchée. Avec les notations suivantes: Bl et Pl les courants SBO et CSB dans la i ème antenne B = Z B Tn Ei j aii et Pn = PT Ei J j L i D le rapport SBO/CSB mesuré sur le j ème capteur, la i détermination des (i s'effectue à partir du système matriciel de J équations à I inconnues: Yj = Wji (i avec yj= Dm pn Bn 1i et W Ji ={Bi Pi} T' Eij la matrice W i est en général dégénérée (déterminant nul), et la résolution implique un processus d'approximation avec un certain choix de contraintes Dans le cas actuel, compte tenu des conditions du problème, on choisit la minimisation de la fonctionnelle: 1 Y II + oc II 1 I I o le paramètre de régularisation a minimise le résidu Il W Y Il suivant la norme quadratique (métrique de L 2). Nous connaissons dans ces équations 3; c'est le second membre de l'équation ( 1) écrite plus haut Il nous faut connaître Bn e et T Le calcul d'obtention de ces valeurs est décrit ci-après. Ce calcul est fait à partir d'une phase d'étalonnage, chaque antenne étant successivement alimentée par un même signal de référence (CSB seul). Considérons le coefficient Tij qui est le résultat de l'émission envoyée vers le capteur j par l'antenne i et l'antenne image i+ 3 de l'antenne i Seule la connaissance en relatif de ces T est nécessaire, ils sont donc directement des complexes dont les modules sont les amplitudes des porteuses recueillies sur les capteurs et les arguments les phases de ces mêmes porteuses. Une acquisition préalable sur les coupleurs d'antennes va nous permettre de mesurer un facteur de couplage pour chacun de ces coupleurs, afin d'obtenir une détermination relative exacte en phase et en amplitude des courants d'alimentation des antennes. On veut dire par là, qu'il existe un rapport complexe de couplage entre l'antenne et son coupleur Connaissant le courant d'alimentation de l'antenne et le courant mesuré sur le coupleur on en déduit le coefficient de couplage Ce coefficient est ensuite considéré comme une constante. L'ensemble des coefficients complexes tli constituent une matrice Soit F et F 2 deux nombres complexes, on peut écrire que: t =F 1 U + F V ji 1 ij 2 (i+ 3) j Afin d'éliminer les incertitudes sur les phases relatives entre les capteurs nous allons traiter 4 groupes (dans le cas de quatre capteurs) de trois équations. Le groupe d'équations correspondant au capteur j s'écrit t U V F j il 4 j h i t U V F 2 j 2 j 5 j 2 j 3 j U V 3 j6 j La résolution de ce système d'équations par régression multilinéaire selon les moindres carrés, nous permet de calculer F les coefficients F et F et donc la valeur du coefficient 2 1 2 F 'j qui est égal au coefficient de réflexion pi du sol On considère que la moyenne des p 1 ainsi calculés constitue le coefficient de réflexion nominal moyen p du sol. A partir des tij mesurés nous pouvons calculer les signaux théoriques reçus sur les capteurs. Pour le signal CSB nous pouvons écrire CS Bj = p 1 Pl et pour le signal SBO SBOJ = Z Bl I Nous calculons ensuite le rapport CSB I Dn j j SB Oi Pour atteindre les grandeurs qui nous intéressent nous allons exploiter par le calcul les écarts entre les signaux théoriques ainsi calculés et les signaux mesurés sur les capteurs soit Dm. Nous avons donc bien dans le cas de l'ILS appliqué le procédé selon l'invention Nous avons créé un modèle mathématique qui, en fonction des courants d'alimentation des trois antennes du réseau de l'émetteur glide, nous permet de calculer le champs en un point dont on connaît les coordonnées Ce champs dépend de quatre paramètres que nous ne savons mesurer directement. Dans ce cas, avant de commencer la phase de mesure proprement dite il convient d'effectuer les phases d'étalonnage que nous venons de décrire. Cette phase d'étalonnage est effectuée dans les conditions de calibration et permet de définir des conditions nominales de fonctionnement Les mesures effectuées au cours de cette phase préalable utilisent le procédé selon l'invention. L'étape 5 du procédé qui consiste à exécuter les opérations l à 4 un certain nombre de fois est ensuite exécuté. Les étapes 2 à 5 du procédé selon l'invention sont exécutées une fois pour chacun des quatres capteurs situés à proximité du réseau d'antennes et une fois pour chacun des coupleurs d'antennes Le signal de référence introduit au niveau du coupleur 9 est le signal CSB. L'étape 6 est ensuite exécutée le calcul permettant de trouver les paramètres inconnus du modèle mathématique. Dans ce qui précède nous avons négligé l'influence du signal de clearance sur les niveaux des signaux CSB et SBO. Afin de connaître les caractéristiques du signal de clearance on peut recommencer l'ensemble des étapes 1 à 4 en prenant comme signal de référence le signal de clearance et en mesurant les courants d'antennes à l'aide des coupleurs 11 et 13 (l'antenne 2 ne reçoit pas de signal de clearance). Ces résultats joints à ceux obtenus par l'emploi du procédé peuvent être utilisés pour calculer la valeur du signal ILS en des points de coordonnées connues. Il est également possible de déterminer le lieu des points o le SBO a une valeur déterminée et en particulier la valeur O. La comparaison de ce lieu avec l'axe de descente permet de contrôler les indications données par le réseau dans tout l'espace de couverture et non plus en quelques points seulement comme c'était le cas auparavant. Le procédé général qui vient d'être décrit comporte des variantes de réalisation Ainsi par exemple si l'on dispose de six capteurs extérieurs, il sera possible en repétant le processus d'acquisition défini par les étapes 1 à 4 pour les deux capteurs supplémentaires de définir des coefficients de réflexion du sol différents pour chacune des antennes d'émission. Si par contre on ne dispose que d'un capteur extérieur, qui est le capteur existant on pourra définir un coefficient de réflexion du sol et admettre que les gains d'antenne sont à leur valeur nominale Cette dernière configuration qui est la configuration la moins onéreuse pour le contrôle de l'intégrité de l'émission du glide d'un ILS en utilisant le procédé selon l'invention permet déjà une bonne amélioration par rapport au dispositif existant. L'emploi du procédé et du dispositif pour le contrôle d'un radiophare VOR va maintenant être décrit. Le signal émis par une station de radionavigation VOR comprend: une première composante contenant l'information de l'azimut de référence (Nord Magnétique de la station) une deuxième composante contenant l'information indicatrice de l'azimut e du récepteur par rapport à l'azimut de référence Le signal total peut s'écrire sous la forme d'une somme de trois vecteurs complexes modulant une fréquence porteuse de pulsation W 0: S = M + A e lest + A e-j Qst o p m Suivant le mode de réalisation, les vecteurs représentatifs s'écrivent (dans l'hypothèse idéale) dans le cas du VOR "conventionnel" M = l 1 + M cos (M t + N sin (QO t))l ejo o i A M e JE p 2 F A M e-je a 2 s o dans le cas VOR Doppler M = ( 1 + M cos Q t)ej o o o A = M ejm sin (ot a) p 2 A M -jn sin (Mt a) Qs i 1 o Les valeurs numériques sont les suivantes: = 2 x 30 Hz 2 = 2 X 9960 Hz *M = 0, 3 X = 16 On constate facilement qu'au renversement près des signaux à phase variable et de référence, le signal total est identique pour les deux modes de réalisation. Les conditions idéales de fonctionnements sont obtenues si: a les vecteurs Ap et An sont complexes conjugués b leur somme est toujours en phase avec la porteuse W c la variation de phase en fonction de l'azimut est aussi proche que possible des valeurs idéales. La réalisation dans le cas du VOR "conventionnel" comprend pour l'essentiel: un émetteur à W modulé M alimentant une antenne o o omnidirectionnelle en amplitude et en phase, un émetteur à (w O + QO) alimentant une antenne omnidirectionnelle en amplitude et dont la phase varie en azimut comme a, un émetteur à (w Q) alimentant une deuxième antenne omnidirectionnelle en amplitude et dont la phase varie en azimut comme -e ces deux antennes sont en général constituées par une combinaison en quadrature de phase de deux antennes ayant des diagrammes de forme cos ( 0) et sin (e), soit: cos ( 0) + j sin ( 0) et cos (e) j sin ( 0) Dans le cas du VOR Doppler la réalisation comprend pour l'essentiel: un émetteur à w> modulé par Mo alimentant une antenne omnidirectionnelle en amplitude et en phase, un émetteur à (w + Q 1) alimentant cycliquement une série d'antennes omnidirectionnelles en amplitude et en phase et disposées sur un cercle de périmètre -q x (A est la longueur d'onde du signal à c). un émetteur à (wo Ol) alimentant cycliquement avec un décallage d'un demi-tour la même série d'antennes. On voit que les performances du VOR dépendent d'un nombre très important de paramètres: électroniques: stabilité en amplitude et en phase des émetteurs antennes: stabilité des diagrammes et des champs rayonnés. Toutes variations des paramètres modifie les valeurs des vecteurs Mo, Ap et Am, la mesure précise de ceux-ci permet de s'assurer en permanence des caractéristiques rayonnées. Les principes de la mesure est le suivant. On considère un signal de la forme M = Rei 1 loo + na) t + &kl a o 2 peut prendre les valeurs 0, DO _ ni, et ak k valeurs équi réparties entre O et 2 ir Il s'agit en pratique d'un signal prélevé à la sortie de chaque émetteur et passant à travers un phaseur possédant k positions de phase. Ce signal est ajouté au signal reçu par un certain nombre de capteurs placés au voisinage des antennes, la détection du signal composite est donc le module du vecteur global, soit: 1 e + Ma + Apejs -Js t avec S prenant les valeurs O ou 1 suivant qu'il s'agit d'un VOR conventionnel ou Doppler. En ne retenant que les quantités ne contenant que a et a k, on obtient les termes suivants: l O 0 RlMo e( at) at + txki)l+ RlAp ejl(Qa s)t l+Ae-(t +k+ A elJ(a t+ ll + RlAm* ejl( a + s)t + k l +m e-Jl(a + s)t + kll Un premier calcul de Fourier sur les N valeurs de la composante continue de la tension détectée va fournir une détermination des vecteurs: M pour 2 a= O; A pour Qa = S et A pour Q = s oa p a S m a s Un second calcul des composantes de Fourier sur les n valeurs de la tension détectée à S pour Qa = 0, c'est-à-dire du module du vecteur: * * Ap (M + Rexp( Ji k)) + Am ( M O + exp (-j k)) va donner le module et la phase de la combinaison de A et A p m Les modalités de mesure et d'exécution de ces calculs sont explicités ci-après dans le cas du VOR Doppler. Dans le cas du VOR Doppler le dispositif employé est celui de la figure 5. La carte d'acquisition 12 comporte une entrée M de synchronisation extérieure, trois coupleurs 1 12 13 permettent de capter le courant d'émission de chacun des émetteurs du VOR et un commutateur HF 22 permet de sélectionner l'un des coupleurs. On dispose d'un seul capteur 2 Le fonctionnement est le suivant le coupleur 91 étant sélectionné par le commutateur 22 commandé par le calculateur 11, le signal à W (correspondant à 2 a = o) en provenance de l'émetteur alimentant l'antenne omnidirectionnelle est introduit sur l'entrée 4 du sommateur 3 L'entrée 5 de ce sommateur est alimenté par le signal en provenance du capteur 2. Ce capteur est placé sur ce qu'on appelle le contrepoid du VOR à l'endroit o se trouve le capteur de "monitoring" actuel le signal reçu par ce capteur est composé du signal émis par l'antenne omnidirectionnelle et du signal émis tour à tour par deux antennes diamétralement opposées émettant un signal non modulé à la pulsation W + Q 1 pour l'une et W O Q 1 pour l'autre. Il y a au total 48 antennes Au cours d'un tour de S chaque 1 30 antenne est alimentée deux fois à ème de seconde d'inter- valle, une fois par un signal à la pulsation W + Q 1 et une fois par un signal à la pulsation W O Q Un organe de synchro- nisation interne au VOR permet la commutation régulière des alimentations des paires d'antenne Le dispositif 21 est un organe de capture d'une partie de ce signal de synchronisation. Ce signal est utilisé par la carte d'acquisition pour déterminer l'instant d'acquisition d'un échantillon qui corres- pondra à l'instant o le signal émis par l'antenne qui émet est à son maximum On obtient ainsi 48 échantillons par tour Chacun des échantillons reçoit une affectation, il y a donc autant d'affectations que d'antennes équiréparties On recommence cette opération N fois en ajoutant à chaque fois au signal de référence un déphasage de - On a obtenu ainsi pour chaque antenne une affectation qui contient N valeurs correspondant à chacune des valeurs de Y On recommence l'opération de la même façon en sélectionnant le coupleur 9 2 qui fournit un signal qui est le signal non modulé à la fréquence W + 2 puis en sélectionnant le coupleur 9 3 qui fournit un signal non modulé à la fréquence &> Q 2 o 1. On obtient donc au total trois affectations par antennes chacune des affectations contenant N valeurs. Nous allons pouvoir effectuer les premiers et seconds calculs définis plus haut. Le calcul de Fourier sur la composante continue de la tension détectée sur les N valeurs de chacune des trois affectations correspondant à la première antenne respectivement pour a = 0, = 2 etn = 2 donneront les valeurs de M, a a i a i A et A. p m Le calcul de Fourier sur la composante à O 1 pour Qa = O sur les N valeurs tirées de la première affectation relative à cette antenne fournira pour cette antenne le module et la phase de la combinaison de A et A p m La même série de deux calculs sur chacune des trois affectations des 47 autres antennes fournira les mêmes éléments pour chacune de ces antennes. Les comparaisons des valeurs obtenues à des gabarits de tolérance permettront de vérifier l'intégrité du signal transmis par la station VOR L'examen de l'ensemble des valeurs de Ap, A et Mo permettra dans certains cas de localiser un disfonction- nement au niveau de l'un des émetteurs, de une ou plusieurs antennes ou au niveau du réflecteur constitué par le contrepoids. REVEND ICATI ONS i Procédé de contrôle permanent de l'intégrité d'une émission susceptible d'être affectée par des paramètres non directement accessibles à la mesure et dont on sait au moyen d'un modèle mathématique déterminer, l'influence sur le module et la phase du champ rayonné dans l'espace, l'émission étant réalisée à partir d'un réseau d'antennes, chaque antenne ou groupe d'antennes du réseau étant alimentée par une ou plusieurs fréquences porteuses H F, chacune de ces fréquences étant modulées d'une façon connue par un ou plusieurs signaux sinusoïdaux basse fréquence purs caractérisé en ce que: 1) on enregistre en au moins un point et pendant un temps T, la valeur de l'amplitude de la combinaison du signal reçu en ce point et d'un signal dit de référence dont la fréquence fondamentale est égale à l'une des fréquences H F du signal à analyser; 2) on extrait l'amplitude de composantes B F du signal ainsi obtenu et on les mémorise; 3) on recommence N fois (n désignant un entier) les étapes i et 2 en enregistrant le signal au même point et en déphasant à chaque fois le signal de référence d'une quantité a k = a k étant un entier prenant successivement toutes les valeurs entre 1 et n 4) on traite les N amplitudes ainsi mémorisées comme les n valeurs d'une fonction F (Cok) dont on extrait les valeurs de certains coefficients de sa transformée de Fourier. ) on recommence l'ensemble des opérations 1 à 4 au moins autant de fois qu'il y a de points de mesure; 6) on détermine les valeurs des paramètres inconnus par utilisation du modèle mathématique. 2 Procédé selon la revendication 1, caractérisé en ce que au moins un des points d'enregistrement se situe à l'in- térieur du radôme d'une antenne, le signal enregistré étant représentatif du courant d'alimentation de cette antenne. 3 Procédé selon la revendication 1, caractérisé en ce que au moins un des points d'enregistrement se situe à proximité du réseau d'antenne. 4 Procédé selon la revendication 1, caractérisé en ce que le temps T d'enregistrement est égal à un nombre entier de périodes du signal à analyser. Procédé selon la revendication 1, caractérisé en ce que les coefficients de la transformée de Fourier de la fonction (yk) comportent au moins l'un des coefficients suivants valeur de la composante continue amplitude et phase de la fréquence fondamentale. 6 Procédé selon la revendication 1, caractérisé en ce que les étapes 1 à 4 sont effectuées avec des signaux de référence différents. 7 Procédé selon la revendication 1 et destiné à contrôler l'intégrité du signal émis par l'émetteur "Glide path" d'un système ILS lequel comporte un réseau constitué de trois antennes, A 1 A 2 A 3 émettant un signal dit CSB, un signal dit SBO et un signal dit clearance, caractérisé en ce que les étapes 1 à 4 sont effectuées en plaçant le point d'enregistrement à l'intérieur de l'une des antennes puis successivement à l'intérieur de chacune des deux autres. 8 Procédé selon la revendication 7, caractérisé en ce que les étapes 1 à 4 sont effectuées en outre en plaçant au moins une fois le point d'enregistrement à proximité du réseau d'antennes. 9 Procédé selon la revendication 8, caractérisé en ce que les étapes 2 à 5 sont effectuées en plaçant le point d'enregistrement à proximité du réseau d'antennes successivement en quatre endroits différents. 10 Procédé selon la revendication 8, caractérisé en ce que les étapes 1 à 4 sont effectuées en plaçant le point d'enregistrement à proximité du réseau d'antennes en six endroits différents. 11 Procédé selon la revendication 8, caractérisé en ce que le signal de référence est le signal dit CSB prélevé au niveau d'un émetteur alimentant le réseau d'antennes. 12 Procédé selon la revendication 7, caractérisé en ce que le temps T d'enregistrement est de 100 millisecondes. 13 Procédé selon la revendication 7, caractérisé en ce que le filtrage défini à l'étape deux est effectué sur la fréquence de 90 Hz. 14 Procédé selon la revendication 13, caractérisé en ce que le filtrage est effectué en outre sur la fréquence de Hz. Procédé selon la revendication 7, caractérisé en ce que le nombre N défini à l'étape 4 est égal à 16. 16 Procédé selon les revendications 7, 8, 11 à 15, caractérisé en ce qu'on crée deux fonctions S(ak) et D(ak) dont les 16 valeurs sont respectivement la somme et la différence des valeurs obtenues pour chacune des valeurs de Ok en filtrant à 90 Hz et 150 Hz et que les coefficients extraits de la trans- formée de Fourier de chacune de ces deux fonctions sont l'amplitude et la phase de la composante à la périodicité 211. 17 Procédé selon les revendications 7-9 et 11 à 15, caractérisé en ce qu'on crée deux fonctions S(k) et D(ak) dont les 16 valeurs sont respectivement la somme et la différence des valeurs obtenues pour chacune des valeurs de k en filtrant à Hz et 150 Hz et que les coefficients extraits de la trans- formée de Fourier de chacune de ces deux fonctions sont l'amplitude et la phase de la composante à la périodicité 211. 18 Procédé selon la revendication 17, caractérisé en ce que les étapes 1 à 4 sont effectuées en outre en prenant le signal de clearance pour signal de référence le point d'enregis- trement étant placé successivement au niveau des antennes A et A 3. 19 Procédé selon la revendication 1 et destiné à contrôler l'intégrité du signal émis par une station VOR Doppler lequel comporte un réseau d'antennes constitué d'une antenne centrale et de 48 antennes équiréparties sur un cercle ayant pour centre l'antenne centrale, l'antenne centrale étant alimentée par un émetteur produisant une fréquence porteuse de pulsation pulsation W O Q 1, caractérisé en ce que le point d'enregis- trement est situé à proximité du réseau d'antenne. Procédé selon la revendication 19, caractérisé en ce que le temps T d'enregistrement est de n/30 ème de seconde. 21 Procédé selon la revendication 19, caractérisé en ce que les valeurs enregistrées au cours d'une des périodes T sont réparties en 48 affectations, les changements d'affectation étant synchronisés avec les changements d'alimentation des paires d'antennes équiréparties. 22 Procédé selon la revendication 19, caractérisé en ce que le nombre N défini à l'étape 3 est de 16. 23 Procédé selon la revendication 19, caractérisé en ce que les étapes 1 à 4 sont effectuées au moins une fois avec chacun des signaux produits par les trois émetteurs comme signal de référence. 24 Procédé selon la revendication 23, caractérisé en ce que le filtrage défini à l'étape 2 est effectué sur la composante continue et la composante Q 1 Procédé selon la revendication 24, caractérisé en ce que les coefficients extraits de la transformée de Fourier de i comportent au moins la composante continue la composante à la fondamentale 21 26 Procédé selon les revendications 19, 20, 21, 23, 24 et 25, caractérisé en ce qu'on crée pour chacune des 48 antennes de la périphérie 3 affectations chacune des affectations contenant les N valeurs enregistrées avec chacun des signaux de référence et qu'on extrait de l'affectation réalisée avec le signal de référence à la pulsation W O la composante continue et la composante à la fréquence fondamentale de la transformée de Fourier de la fonction de ak et qu'on extrait de chacune des deux autres affectations la composante continue de la transformée de Fourier de la fonction de Yk. 27 Dispositif pour effectuer la mesure de paramètres d'une émission non directement accessibles à la mesure dont on sait au moyen d'un modèle mathématique déterminer l'influence sur le module et la phase de l'émission en un point, l'émission étant réalisée à partir d'un réseau d'antennes, chaque antenne ou groupe d'antennes du réseau étant alimentée au moyen de un ou plusieurs émetteurs par une fréquence porteuse HF commune, cette fréquence étant pour chacune des antennes modulée par un ou plusieurs signaux sinusoïdaux BF purs caractérisé en ce que le dispositif comprend: au moins un capteur ( 2) placés à proximité du réseau d'antenne; un sommateur ( 3) comportant deux entrées ( 4, 5) et une sortie ( 6) et recevant sur l'une de ses entrées ( 5) par l'intermédiaire d'un commutateur programmable HF ( 7) permettant le raccordement successif de chacun des capteurs ( 2),l'un des signaux en provenance d'un capteur et sur l'autre entrée ( 4) par l'intermédiaire d'un phaseur programmable ( 8) un signal dit de référence capté par l'intermédiaire d'un coupleur ( 9) en sortie de l'un des émetteurs concourant à l'émission, la sortie du sommateur étant raccordée à un détecteur ( 10); un calculateur ( 11) recevant par l'intermédiaire d'une carte d'acquisition ( 12) le signal en sortie ( 16) du détecteur ( 10) et envoyant par l'intermédiaire d'une carte interface ( 13) des signaux de commande en direction du phaseur et en direction du commutateur haute fréquence au moyen d'un bus ( 17). 28 Dispositif selon la revendication 27, caractérié en ce que la sortie ( 16) du détecteur ( 10) est appliqué à l'entrée d'une pluralité de filtres analogiques ( 14) la sortie de chacun des filtres ( 14) alimentant un commutateur programmable ( 15) dont la sortie alimente la carte d'acquisition ( 12). 29 Dispositif selon la revendication 27, caractérisé en ce que un atténuateur programmable ( 18) est branché entre la sortie du commutateur ( 7) et l'entrée du sommateur ( 3). Dispositif selon la revendication 29, caractérisé en ce que un interrupteur à diode ( 19) est inséré entre la sortie du phaseur ( 8) et l'entrée ( 4) du sommateur ( 3). 31 Dispositif selon la revendication 30, caractérisé en ce que l'entrée du commutateur ( 7) est raccordée en outre à trois coupleurs ( 11, 12, 13) et une charge ( 20) 32 Dispositif selon la revendication 31, caractérisé en ce qu'il comporte un commutateur programmable ( 22) dont la sortie est raccordée à l'entrée du phaseur ( 8) et dont les entrées sont raccordées à au moins deux coupleurs ( 91 92). 1 33 Dispositif selon les revendications 28 et 32, caractérisé en ce que la carte d'acquisition ( 12) est alimentée par un organe ( 21) de synchronisation. 34 Dispositif selon l'une des revendications 27 à 33, caractérisé en ce que le sommateur ( 3) et le détecteur ( 10) sont remplacés par un mélangeur ( 23) recevant sur l'une de ses entrée ( 5) le signal en sortie de l'atténuateur ( 18) et sur l'autre entrée ( 4) le signal en sortie du phaseur ( 8), et un circuit d'adaptation ( 24) recevant en entrée le signal en provenance du mélangeur ( 23) et alimentant en sortie la carte d'acquisition ( 12).