FR 2695211 A1 19940304 FR 9210290 A 19920826 DISPOSITIF ET PROCEDE D'ANALYSE DE SIGNAUX ILS La présente invention concerne un procédé et un dispositif d'analyse de signaux de guidage d'atterrissage connus sous l'appellation de signaux ILS "Instruments Landing System" Le dispositif est particulièrement destiné au contrôle au sol de toutes les caractéristiques des installations d'atterrissage aux instruments. Dans la demande de brevet française FR-A-1487 441, la demanderesse a décrit un procédé de mesure d'un signal ILS mettant en oeuvre une technique analogique évaluant l'écart de l'avion par rapport à l'axe d'atterrissage, à partir de la détermination de la différence de modulation existant entre les deux composantes à 90 et 150 Hz constituant le signal composite ILS. La demande de brevet française FR-A-2 596 547 montre un dispositif pour élaborer des informations de radionavigation reposant, dans un mode de réalisation, sur une architecture numérique mettant en oeuvre un ensemble à microprocesseur au niveau duquel le signal ILS est traité selon un algorithme faisant appel à un calcul de transformée de Fourier rapide récursive. Ces deux demandes, qui décrivent des dispositifs destinés à servir de récepteurs de navigation embarqués, engendrent des contraintes distinctes de celle d'un dispositif de mesure et d'analyse; notamment ils obligent à une intégration des défauts dus entre autres aux obstacles pouvant apparaître sur le trajet d'atterrissage, afin d'éviter tout louvoiement lors de l'approche Ce "lissage" des informations ne permet que des cadences de renouvellement des données faibles qui, par contre, autorisent le recours à des traitements tels que la transformée de Fourier rapide récursive, laquelle exige des temps d'observation longs. Au contraire, dans un dispositif d'analyse tel que celui de l'invention, il est essentiel de pouvoir appréhender toutes les imperfections du système à contrôler, ce qui nécessite des temps de réponse particulièrement rapides. Un but de l'invention est de réaliser un dispositif de mesure présentant des performances élevées, c'est-à-dire une grande précision et une rapidité de traitement non moins grande. Un autre but de l'invention est de réaliser un dispositif simple et compact mettant en oeuvre un nombre de composants restreint et autorisant, de ce fait, une maintenance rapide et réduite. Ces buts sont atteints par un dispositif d'analyse de signaux de guidage d'atterissage aux instruments (ILS) comportant un convertisseur analogique/numérique recevant le signal à analyser et délivrant une succession de valeurs aptes à être traitées par un ensemble numérique, caractérisé en ce que ledit 2 2695211 signal à analyser constitue en outre une entrée d'une boucle à verrouillage de phase dont la sortie constitue une entrée d'un synthétiseur de fréquence qui délivre un signal d'échantillonnage pour ledit convertisseur analogique/numérique. La boucle à verrouillage de phase est synchronisée sur la fréquence du signal à analyser extraite d'une sub-harmonique des composantes à 90 Hz ou 150 Hz constituant ce signal. Cette architecture spécifique permet de synchroniser le processus de calcul, mis en oeuvre au niveau de l'ensemble de traitement numérique, sur le signal de façon à appréhender constamment l'amplitude maximale des composantes à 90 et 150 Hz formant du signal ILS, et cela malgré la dérive en fréquence existant au niveau de ce signal. L'ensemble numérique comporte une unité de traitement, une mémoire vive, une mémoire non volatile et un module d'entrées/sorties et le traitement des valeurs issues du convertisseur analogique/numérique est réalisé en temps réel entre deux instants d'échantillonnage. Le traitement instantané des valeurs numériques acquises permet de limiter la capacité mémoire du dispositif et ainsi de lui assurer une compacité maximale. La mémoire non volatile comporte une première table dans laquelle sont inscrites des valeurs numériques correspondant à la valeur d'une fonction fenêtre de Hamming pour chaque instant d'échantillonnage, ainsi qu'au moins deux tables dans chacune desquelles sont inscrites des valeurs numériques représentatives, sur une période un quart, d'une fonction cosinus, de fréquence égale à celle d'une composante caractéristique du signal à analyser. Le recours à ces tables de valeurs évite un calcul, pour chaque échantillon acquis, des valeurs des fonctions mises en oeuvre dans le processus de calcul. De préférence, le dispositif d'analyse selon l'invention comporte un filtre numérique utilisant une fonction de pondération dont l'équation est: W(t) = 2 O la + ( 1-c) cos O l avec =T, et qui réalise une extraction d'un signal d'identification du signal ILS à analyser o T représente la durée d'observation du signal à analyser et k' et a sont des coefficients numériques. Le choix de cette fenêtre (et notamment avec k' = 4 et a = 0,54) permet l'obtention d'un filtre ayant une réponse en fréquence très sélective et quasi-linéaire dans la bande passante. Comme précédemment, cette fenêtre de pondération est avantageusement 3 2695211 emmagasinée sous la forme d'une table de valeurs numériques. Le dispositif selon l'invention est mis en oeuvre conformément à un procédé comportant les étapes suivantes: a) acquisition d'un échantillon S(n Te) du signal ILS à analyser par la conversion analogique/numérique de ce signal, Te représentant la période d'échantillonnage et N étant un nombre entier, b) multiplication de la valeur de l'échantillon prélevé par une première valeur numérique correspondant à la valeur d'une fonction fenêtre à l'instant d'acquisition: S(n Te) = Fen(n Te) x S(n Te) c) détermination, aux fréquences k des composantes caractéristiques du signal ILS, des parties réelles et imaginaires du spectre du signal S(t) échantillonné et pondéré: R(k) = R(k) + S(n Te) x cos 2 nk n Te et I(k) = I(k) + S(n Te) x sin 2,nk n Te d) répétition des étapes a) à c) sur une durée d'observation T du signal ILS, e) calcul des taux de modulation des composantes caractéristiques du signal ILS. M(k) = 4,r 4 R 2 (k) + (k), lHF étant la valeur moyenne de S(n Te), i) calcul des différence (DDM) et somme (SDM) des modulations (M( 90) et M( 150)) des composantes à 90 et 150 Hz du signal ILS: SDM = M( 90) + M( 150) et DDM = M( 90) M( 150). Dans un mode particulier de réalisation, il peut n'être prélevé qu'un échantillon du signal à analyser toutes les deux ou trois valeurs des tables de la fonction cosinus, ce qui permet ainsi de déterminer également l'amplitude des composantes harmoniques de ce signal. D'autres caractéristiques et avantages de la présente invention ressortiront mieux à la lecture de la description qui va suivre faite, à titre indicatif et non limitatif, en regard des dessins annexés, sur lesquels: la figure 1 montre la structure d'un dispositif d'analyse de signaux ILS selon l'invention, les figures 2 à 4 sont une approche graphique permettant d'observer les transformations réalisées sur le signal ILS dans le domaine temporel et dans le domaine fréquentiel ainsi que les relations entre ces deux domaines, la figure 5 est un organigramme décrivant le traitement numérique des signaux ILS dans le dispositif selon l'invention, la figure 6 montre une forme du signal emmagasiné dans les tables du dispositif et intervenant dans le traitement numérique décrit à la figure 4, la figure 7 est un organigramme décrivant le traitement du signal d'identification dans le dispositif selon l'invention, la figure 8 représente, dans le domaine temporel, une fonction "fenêtre de pondération" mise en oeuvre dans le traitement du signal d'identification, et la figure 9 représente la réponse en fréquence du filtre digital correspondant à la fenêtre de pondération de la figure 8. On sait que, dans un système de guidage d'atterrissage aux instruments, on dispose de deux émetteurs; l'un, à la fréquence de 110 M Hz, connu sous le nom de signal de radioalignement de piste ou "localizer" et émettant un faisceau de 2,5 ' de largeur de part et d'autre de l'axe de la piste d'atterrissage, et l'autre, à la fréquence de 332 M Hz, connu sous l'appellation de signal de radioalignement de descente ou "glide- path" et émettant un faisceau de 0,7 ' de largeur autour d'un axe de descente typique de 3 '. En pratique, ces faisceaux sont constitués chacun par la zone commune de deux faisceaux principaux modulés à basse fréquence, l'un à 90 Hz l'autre à 150 Hz. Le signal ILS résultant se présente donc en théorie sous la forme d'un signal haute fréquence ( 110 ou 332 M Hz) modulé en amplitude par deux tensions sinusoïdales 90 et 150 Hz. Dans le cas particulier du signal ILS d'alignement de piste, un signal d'identification à 1020 Hz vient en outre se superposer aux composantes à 90 et Hz. La figure 1 est un exemple préférentiel de réalisation d'un dispositif de traitement de tels signaux ILS, conforme à l'invention. Ce signal ILS constitue un signal d'entrée d'une part d'un filtre passebas antirepliement 1 et d'autre part d'une boucle à verrouillage de phase (PLL)2 dont la sortie est reliée à un synthétiseur de fréquence 3 La sortie du filtre passe-bas 1 et du synthétiseur 3 constituent des entrées d'un ensemble de traitement numérique 4. De façon connue, une boucle à verrouillage de phase est constituée d'un comparateur de phase dont la sortie commande un oscillateur contrôlé en tension (OCI) qui réagit à son tour sur l'entrée du comparateur Dans le cadre de la présente invention, la synchronisation de ce circuit est réalisée, non pas sur une fréquence fondamentale du signal analysé, comme il est fréquement pratiqué, mais sur une sub-harmonique de ce signal En l'espèce, il s'agit de la sub-harmonique 3 du signal Hz ou de la sub-harmonique 5 du signal 150 Hz, c'est-à-dire de la fréquence Hz. Egalement de façon connue, un synthétiseur de fréquence se présente comme une boucle à verrouillage de phase au sein de laquelle un diviseur programmable est inséré dans le circuit de contre-réaction. Dans l'exemple de l'invention, le diviseur programmable est un diviseur par 540 permettant d'obtenir une fréquence de 16,2 k Hz en sortie du synthétiseur 3 à partir de la fréquence stable de 30 Hz fournie par la boucle à verrouillage de phase 2. L'ensemble de traitement 4 comporte un convertisseur analogique/numérique qui reçoit le signal ILS issu du filtre antirepliement 1 et délivre sur 10 bits un signal échantillonné S*(t) au rythme d'un signal d'échantillonnage délivré par un module d'entrées/sorties 11 Le résultat de cette conversion est emmagasiné dans une mémoire vive 12 qui est reliée à une unité de traitement 13, laquelle est également reliée à une mémoire non volatile 14. Le module d'entrées/sorties 11 qui reçoit le signal de synchronisation de 16,2 k Hz généré par le synthétiseur de fréquence 3 est relié à l'unité de traitement 13 et délivre également des signaux d'exploitation du signal ILS analysé. Les transferts d'information entre l'unité de traitement et ses circuits associés, mémoires 12,14 ou module d'entrées/sorties 11 sont avantageusement réalisés sur 16 bits. La mémoire non volatile contient un programme de gestion de l'ensemble de traitement ainsi que plusieurs tables de valeurs dont l'utilité apparaîtra ci-après, au regard des figures 2 à 9 qui explicitent le fonctionnement de l'unité de traitement selon l'invention. Les figures 2 et 3 montrent les différentes opérations effectuées, par l'ensemble de traitement, sur le signal tant dans le domaine temporel que dans le domaine fréquentiel. La figure 2 a montre le signal ILS en entrée du dispositif selon l'invention. Dans un souci didactique, le signal analysé est le signal d'alignement de descente qui ne comporte que les composantes à 90 Hz et 150 Hz à l'exclusion du signal d'identification à 1020 Hz Mais il est évident qu'une analyse comparable peut être menée avec le système d'alignement de piste dont le signal démodulé ferait alors également ressortir ce signal à 1020 Hz En outre, le signal sera tout d'abord supposé parfait, sans harmoniques. La forme du signal démodulé est donnée par l'équation suivante: S(t) = HF ( 1 + M 90 sin 2 x 90 t + M 150 sin 2 iî 150 t) avec HF la composante continue image de la porteuse haute fréquence à 332 M Hz à 6 2695211 laquelle est superposée la somme des ondes 90 Hz et 150 Hz, M 90 le taux de modulation de l'onde 90 Hz, M 150 le taux de modulation de l'onde 150 Hz. Ces taux de modulation servent à reconstituer la différence des modulations (DDM) et la somme des modulations (SDM), paramètres qui, pour le premier, caractérise la position de l'avion par rapport à l'axe de descente (ou de piste) et, pour le second, sert à la validation de cette précédente mesure. L'échantillonnage du signal S(t) à une fréquence d'échantillonnage Fe = l/Te est obtenu par la multiplication de ce signal par un peigne de Dirac A(t) de période Te (voir figure 2 b) soit S*(t) = S(t) x A(t) (figure 2 c), S* étant une suite de valeurs de S(t) prises aux instants n Te. Dans le domaine fréquentiel, cet échantillonage correspond à une convolution entre le signal S(F) (figure 3 a) et le signal A(F) (figure 3 b), laquelle donne un spectre S*(F) qui est la répétition de S(F) autour des fréquences multiples de Fe (voir figure 3 c). On peut alors observer que pour échantillonner, sans perte d'information, ce signal S(t) dont le spectre est borné à une fréquence maximale Fm; il faut impérativement que la fréquence d'échantillonnage Fe soit supérieure à 2 Fm. Mais il faut être sûr que le spectre de S(t) soit borné à cette valeur Fm, or ce signal est soumis à des variations rapides et des parasites multiples dépendant de la configuration de la piste au moment de l'atterrissage C'est pourquoi, pour éviter ces inconvénients, l'échantillonnage a été précédé d'un filtre passe-bas qui limite obligatoirement le spectre à Fm, évitant ainsi toutes perturbations non désirées; c'est le filtre antirepliement 1 Le spectre des fréquences en sortie de ce filtre est identique à celui de S(t) et c'est ce signal, et non S(t), qui fait en réalité l'objet de l'échantillonnage précité. Les figures 2 d et 3 d représentent respectivement dans le domaine temporel et dans le domaine fréquentiel une fonction de pondération connue sous le nom de fenètre de Hamming. L'échantillonnage ne peut être réalisé sur l'ensemble du signal S(t), c'est-à- dire pour un nombre infini d'échantillons correspondant à un temps d'observation du signal infini Le nombre d'échantillons doit être fini Cette limitation réduit la durée d'observation à une valeur T, ce qui équivaut à multiplier la fonction échantillonnée S*(t) par une fonction fenêtre de largeur T Afin de réduire les effets dus à cette troncature, il est intéressant de recourir à une fenêtre spécifique telle que la fenêtre de Hamming qui présente, au niveau de son spectre fréquentiel, des lobes latéraux secondaires de très faibles amplitudes, le lobe principal ayant une largeur de 1/T. Cette fonction fenêtre est définie par l'équation suivante: ham(t) = a + ( 1-a) cos ( 2 z T t) De préférence, on choisira pour a la valeur 0,54, sans toutefois être limité à cette seule valeur (la valeur a= 0,5 qui correspond à une fenêtre de Hann pourrait, par exemple, être également envisagée). Les figures 2 e et 3 e montrent l'allure du signal S(t) après échantillonnage et troncature et celle du spectre correspondant. Celui-ci est donné par la convolution du spectre du signal échantillonné S*(F) et du spectre de la fenêtre ham (F) soit: S*(F) * ham (F) Dans le domaine fréquentiel, cette convolution revient à répéter le spectre ham (F) autour des différentes fréquences du spectre S*(F) On peut alors observer, que pour éviter toute perte d'information et obtenir une sélection parfaite des fréquences 90 et 150 Hz, il est nécessaire que la largueur fréquentielle du spectre de la fenêtre de pondération ne soit pas supérieure à 60 Hz Dans le cas de la présente invention, cela impliquerait i 1 une observation théorique du signal pendant 16,6 ms (r T = 30 Hz soit T = 60 seconde), soit une demi-période du signal ILS qui présente une périodicité de 33,3 ms ( 30 Hz). En pratique, il en va toutefois différemment En effet, du fait de défauts présents au niveau de la modulation et de la génération des signaux 90 et 150 Hz, le signal ILS comporte également les composantes harmoniques de ces signaux et notamment, comme le montre la figure 4, les harmoniques 2 du 150 Hz ( 75 Hz) et 2 du 90 Hz ( 180 Hz) Donc, pour obtenir en pratique une sélection parfaite des signaux utiles à 90 et 150 Hz, il sera nécessaire que la largeur fréquentielle de la fenêtre d'observation ne soit pas supérieure à 15 Hz, ce qui correspondant à une largeur temporelle de 66,6 ms ( 2 T = 7,5 Hz soit T = 1 seconde). La figure 5 montre, sous forme d'organigramme, les différentes opérations exécutées par l'ensemble de traitement 4 pour déterminer les paramètres différence de modulation DDM et somme des modulations SDM caractéristiques du signal ILS analysé. Ces paramètres sont fonction des taux de modulation M 90 et M 150 qui peuvent être déterminés simplement à partir de la transformée de Fourier du signal s(t). HF.M 90 En effet le module de S(F) pour 90 Hz est égal à 2 et, pour 150 Hz, il HF.M 150 est égal à 2, la composante continue ayant pour module HF. HF.M 90 HF M 150 Il vient donc IS( 90)1 = 2 et IS( 150)= 2 2 2 2 2 soit encore M 90 = WS( 90 " HF R 2 ( 90)+I 2 ( 90) et M 150 = O S( 150 = = 2 R 2 ( 150)+I 2 ( 150) R( 90), I( 90),R( 150), 1 ( 150) étant les parties réelles et imaginaires des modules de S(F) aux fréquences 90 et 150 Hz. Il convient donc de déterminer ces parties réelles R (F) et imaginaires I (F) pour retrouver ensuite simplement les taux de modulation Or R(F) et I(F) sont donnés par les formules suivantes: R(F) =f S(t) cos 2 z Ft dt I(F) = S(t) sin 2 z Ft dt Lesquelles peuvent également s'écrire, en tenant compte de l'échantillonnage et de la troncature, sur la période T, amenée par la fenètre de pondération: R(k) = YN S(n Te) cos 2 zkn Te I(k) = EN S(n Te) sin 2 zkn Te N étant le nombre d'échantillons prélevés sur la période T, Te étant la période d'échantillonnage et N représentant un nombre entier. Ces calculs seront effectués à partir de chaque échantillon du signal ILS et de tables stockées dans la mémoire non volatile 14 et comportant chacune une suite de valeurs définissant une courbe cosinus ou sinus à 90 ou 150 Hz. La figure 6 montre un exemple d'une telle table Les valeurs sont emmagasinées dans la table sous forme d'entiers de 16 bits et varient donc entre 215 et -215 Avantageusement, chaque table, 90 ou 150 Hz, comportera une période un quart de sinusoïde En effet, en commençant au début de la table, on obtient un cosinus sur N points, alors qu'en commençant après le premier quart de période, on obtient, également sur N points, un sinus. La fréquence d'échantillonnage du signal S(t) étant fixé à 16,2 k Hz, soit une acquisition d'un échantillon toutes les 61,728 is, le nombre d'échantillons prélevé pendant la durée d'observation du signal sera égal à 1080, cette durée ayant été fixée 9 2695211 préalablement à deux périodes du signal ILS soit 66,6 ms ( 15 Hz). Une période du signal 90 Hz pourra alors être décrite exactement sur 180 points, ( 16200/90), la table comportant alors au total 225 points ( 1 période + 1/4 de période) De même, une période du signal 150 Hz pourra être décrite exactement sur 108 points ( 16200/150), la table comportant alors un maximum de 135 points. Il peut être noté qu'il est possible, en n'extrayant de ces tables qu'un point sur deux ou trois, de recueillir des informations sur le taux d'harmoniques 2 ou 3 existant au niveau du signal, lorsque celui-ci ne se présente pas sous la forme d'un signal à spectre idéal tel que représenté à la figure 3 a. Il est fait à nouveau référence à la figure 5 qui décrit le programme de calcul des taux de modulation M 90 et M 150 Il s'agit d'une boucle synchronisée sur le signal d'échantillonnage à 16,2 k Hz. Après une première étape 20, consistant en l'initialisation des différents paramètres nécessaires aux calculs, il est procédé dans une seconde étape 21 à l'attente du signal de synchronisation commandant l'échantillonnage du signal S(t). A réception de cet ordre, il est procédé dans une étape 22 à l'acquisition d'un premier échantillon S(n Te) avec n= 1, Te étant la période d'échantillonnage. Dans une étape suivante 23, la valeur de S (n Te) obtenue est multipliée par la valeur de la fenêtre de pondération à l'instant correspondant Ham (n Te), cette dernière valeur étant extraite d'une première table 140 comportant des valeurs de la fonction Ham (t) prédéterminées pour chaque instant d'échantillonnage Il est procédé ensuite dans une étape 24 à une première détermination du niveau continu HF alors égal à l'amplitude A du signal relatif au premier échantillon acquis. L'étape 25 réalise le calcul des parties réelle et imaginaire de la composante à 90 Hz du signal ILS. La partie réelle R 90 est donnée en multipliant la valeur échantillonnée et pondérée déterminée au calcul précédent par cos(n Te) et la partie imaginaire I 90 en multipliant cette valeur par sin(n Te), les valeurs de cos(n Te) et sin(n Te) étant extraites d'une table unique 90 Hz 142. L'étape 26 réalise un calcul identique pour la composante à 150 Hz du signal ILS, les valeurs cos(n Te) et sin (n Te) étant extraites d'une table unique 150 Hz 143. Le nombre d'échantillons N est ensuite incrémenté à l'étape 27 et les opérations précédentes des étapes 21 à 26 peuvent être recommencées pour un nouvel échantillon et ainsi de suite jusqu'à achèvement à l'étape 28 du temps dobservation qui est limité avantageusement à 66,6 ms soit exactement 1080 points du signal S(t). Les taux de modulations sont alors calculés à l'étape 29 par les formules suivantes déjà décrites: 2 2 n M 90 = i(S( 90))"= \R 2 ( 90)+I 2 ( 90) et M 150 = = 2 n(S( 150)= R R 2 ( 150)+ 12 ( 150) Le rapport A/n donnant la valeur moyenne HF du signal S(t). Enfin, à l'étape 30, préalablement à l'étape 31 de fin de traitement, les différence et somme des modulations DDM et SDM sont déterminées par différence ou somme des taux de modulation calculées à l'étape précédente 29 Ces données seront avantageusement délivrées sur douze bits. Il est fondamental de noter que le signal échantillonné n'est en aucune façon mémorisé, les calculs se faisant pendant la période séparant deux acquisitions du signal Le procédé est réalisé en temps réel et de ce fait ne nécessite pas de capacité mémoire importante Ainsi, la mémoire vive contenant, notamment le résultat de conversion, pourra n'être que de 256 octets. Ces calculs devant être réalisés entre deux instants d'échantillonnage soit 61,728 Ms, l'unité de traitement autorisera de préférence des calculs sur des entiers de 16 et 32 bits ainsi qu'en virgule flottante Un fonctionnement à 16 M Hz sera avantageusement envisagé sans que cette fréquence de calcul soit limitative. La figure 7 décrit, à l'aide d'un organigramme, l'analyse du signal d'identification à 1020 Hz. Ce signal est une onde audio modulée selon le code morse et superposée au signal ILS de radioalignement de piste. Dans ce cas, le signal S(t) présentera donc une composante spectrale légèrement distincte de celle de la figure 3 b avec, en supplément, notamment une raie à la fréquence 1020 Hz d'amplitude HF M 1020 Il est donc possible, en reprenant les principes définis pour le calcul des taux de modulation 90 et 150 Hz, d'évaluer également le taux de modulation 1020 Hz Toutefois, le signal d'identification pouvant varier dans une plage de + 50 Hz par rapport à sa valeur nominale de 1020 Hz, la fenêtre de Hamming définie précédemment ne permet plus de recueillir valablement l'amplitude de cette raie de fréquence. Il est nécessaire de définir une nouvelle fenêtre de pondération qui permettra un filtrage de la raie de fréquence 1020 Hz avec un bande passante d'au moins Hz et si possible une ondulation quasi inexistante. La figure 8 montre la fenêtre de pondération développée pour le filtrage de la il 2695211 composante 1020 Hz. Cette fenêtre présente une forme en (sin n)/n amorti résultant de la combinaison d'une fonction en (sin n)/n avec une fonction de Hamming L'équation de cette fonction fenêtre est donnée par: W 1020 (t) = 2 O la + ( 1-a) cos O l avec y T o T représente la durée d'observation et a étant avantageusement choisi égal à 0,54 et k' est égal à 4. L'application de cette fenêtre spécifique W 1020 (t) au signal échantillonné S*(t) définit un filtre digital dont la réponse en fréquence est représentée figure 9, et sur laquelle on peut constater l'absence d'ondulation dans la bande passante ainsi que la raideur des flancs de ce filtre, qui montre bien sa très grande sélectivité. Il est fait de nouveau référence à la figure 7 qui explicite le calcul du taux de modulation M 1020 Comme pour la détermination des taux M 90 et M 150, le calcul est réalisé dans une boucle synchronisée à 16,2 k Hz. Après une première étape 40 d'initialisation des différents paramètres de calcul, il est procédé dans une seconde étape 41 à l'attente du signal de synchronisation commandant l'échantillonnage et donc la conversion analogique/numérique du signal ILS Dès réception de cet ordre de conversion, il est procédé dans une étape 42 à l'acquisition du premier échantillon, puis dans une étape 43 à la multiplication de la valeur de cet échantillon par la valeur de la fenêtre de pondération précitée W 1020 pour cet instant d'échantillonnage, cette dernière valeur étant extraite d'une seconde table 141 comportant l'ensemble des valeurs définissant la fonction W 1020 (t) Il est procédé ensuite dans une étape 44 à une première détermination du niveau continu HF L'étape suivante 45 réalise le calcul des parties réelle et imaginaire de la composante à 1020 Hz par le produit de la valeur échantillonnée et pondérée par des valeurs extraites d'une table unique 1020 Hz 144 délivrant des valeurs de cosinus et sinus pour chaque instant d'échantillonnage Une nouvelle évaluation est ensuite possible après passage, dans l'étape 46, à un échantillon suivant, ces calculs étant effectués sur toute la durée d'observation, c'est-à-dire sur 1080 points (étape 47). Le taux de modulation peut alors être calculé, à l'étape 48, à partir des valeurs de R( 1020) et I( 1020) et de la composante continue HF, l'étape 49 terminant le traitement de la composante à 1020 Hz Le taux de modulation à 1020 Hz est comme les valeurs DDM et SDM avantageusement délivré sur 12 bits. Il peut être noté que le contenu de la table sin/cos 1020 (il en est de même pour les tables 90 et 150 Hz) est accessible très facilement en utilisant deux pointeurs l'un se déplaçant cycliquement sur les N premierspoints de la table, définissant ainsi la fonction cosinus, et l'autre se déplaçant cycliquement sur N points à partir du premier quart de période de cette table, pour définir la fonction sinus. Contrairement aux tables 90 et 150 Hz pour lesquelles chaque période peut être définie avec un nombre entier multiple de la période d'échantillonnage, une période de la table 1020 Hz n'est pas égale à un nombre entier de fois cette période d'échantillonnage ( 16200/1020 = 15, 882) On choisira donc l'entier le plus proche soit 16 pour définir une période du signal d'identification, ce qui ramène en pratique la fréquence de ce signal à 1012,5 Hz, valeur encore comprise dans la plage de variation admissible ( 1020 + 50 Hz). Il en résulte que la table 1020 Hz comporte 16 + 4 soit 20 points au total. Il est clair que la structure ainsi décrite permet des calculs d'une grande précision et cela d'autant plus que ces calculs sont synchronisés sur le signal par l'ensemble de synchronisation à boucles à verrouillage de phase On est ainsi sûr de se positionner constamment sur le maximum des lobes principaux et donc, corrélativement, d'éviter toute erreur d'amplitude préjudiciable à la détermination des taux de modulation et donc à la définition des axes de descente et de piste. En outre, la précision de la détermination des paramètres DDM et SDM est encore amplifiée par le fait que la mesure est réalisée à partir des taux de modulation réels définis comme le rapport des amplitudes des signaux modulants par rapport à l'amplitude de la porteuse, et non en considérant cette dernière amplitude comme constante et en évaluant alors seulement l'amplitude des signaux modulants, comme il était courant dans les procédés de l'art antérieur, notamment dans les brevets décrits en tête de la description. REVENDICATIONS 1 Dispositif d'analyse de signaux de guidage d'atterrissage aux instruments (ILS) comportant un convertisseur analogique/numérique recevant le signal à analyser et délivrant une succession de valeurs aptes à être traitées par un ensemble numérique, caractérisé en ce que ledit signal à analyser constitue en outre une entrée d'une boucle à verrouillage de phase ( 2) dont la sortie constitue une entrée d'un synthétiseur de fréquence ( 3) qui délivre un signal d'échantillonnage pour ledit convertisseur analogique/numérique ( 10). 2 Dispositif d'analyse de signaux ILS selon la revendication 1, caractérisé en ce que ladite boucle à verrouillage de phase est synchronisée sur la fréquence du signal à analyser extraite d'une subharmonique des composantes à Hz ou 150 Hz constituant ce signal. 3 Dispositif d'analyse de signaux ILS, selon la revendication 1, caractérisé en ce que ledit ensemble numérique comporte une unité de traitement ( 13), une mémoire vive ( 12), une mémoire non volatile ( 14) et un module d'entrées/sorties ( 11) et en ce que le traitement des valeurs issues du convertisseur analogique/numérique ( 10) est réalisé en temps réel entre deux instants d'échantillonnage. 4 Dispositif d'analyse de signaux ILS, selon la revendication 3, caractérisé en ce que ladite mémoire non volatile comporte une première table ( 140) dans laquelle sont inscrites des valeurs numériques correspondant à la valeur d'une fonction fenêtre de Hamming pour chaque instant d'échantillonnage. Dispositif d'analyse de signaux ILS selon la revendication 3, caractérisé en ce que ladite mémoire non volatile comporte au moins deux tables ( 142,143) dans chacune desquelles sont inscrites des valeurs numériques représentatives, sur une période un quart, d'une fonction cosinus de fréquence égale à celle d'une composante caractéristique du signal à analyser. 6 Dispositif d'analyse de signaux de guidage d'atterrissage aux instruments (ILS) comportant un convertisseur analogique/numérique recevant le signal à analyser et délivrant une succession de valeurs aptes à être traitées par un ensemble numérique, caractérisé en ce qu'il comporte un filtre numérique utilisant une fonction de pondération dont l'équation est donnée par la relation suivante: sin k'0 2 xrt W(t) = 2 o la + ( 1-a) cos 6 l avec O = T, et qui réalise une extraction d'un signal d'identification du signal ILS à analyser, o T représente la durée d'observation du signal à analyser et k' et a sont des coefficients numériques. 7 Dispositif d'analyse de signaux ILS selon la revendication 6, caractérisé en ce que le coefficient a est égal à 0,54 et le coefficient k' égal à 4. 8 Dispositif d'analyse de signaux ILS selon la revendication 6 ou la revendication 7, caractérisé en ce que ladite fonction de transfert est emmagasinée sous forme numérique dans une table ( 141) contenue dans une mémoire non volatile ( 14) d'un ensemble de traitement ( 4) réalisant un traitement numérique du signal ILS. 9 Procédé d'analyse de signaux de guidage d'atterrissage aux instruments (ILS), caractérisé en ce qu'il comporte les étapes suivantes: a) acquisition d'un échantillon S(n Te) du signal ILS à analyser par la conversion analogique/numérique de ce signal, Te représentant la période d'échantillonnage et N étant un nombre entier, b) multiplication de la valeur de l'échantillon prélevé par une première valeur numérique correspondant à la valeur d'une fonction fenêtre à l'instant d'acquisition: S(n Te) = Fen(n Te) x S(n Te) c) détermination, aux fréquences k des composantes caractéristiques du signal ILS, des parties réelles et imaginaires du spectre du signal S(t) échantillonné et pondéré: R(k) = R(k) + S(n Te) x cos 2 tk n Te et I(k) = I(k) + S(n Te) x sin 2 xk n Te d) répétition des étapes a) à c) sur une durée d'observation T du signal ILS, e) calcul des taux de modulation des composantes caractéristiques du signal ILS. M(k) = V- RZ(k) + I 2 (k), HF étant la valeur moyenne de S(n Te). Procédé d'analyse de signaux ILS selon la revendication 9, caractérisé en ce qu'il comporte en outre l'étape suivante: f) calcul des différence (DDM) et somme (SDM) des modulations (M ( 90) et M( 150)) des composantes à 90 et 150 Hz du signal ILS: SDM = M( 90) + M( 150) et DDM = M( 90) M( 150). 11 Procédé d'analyse de signaux ILS selon la revendication 9, caractérisé en ce que ladite fonction fenêtre est une fonction de Hamming. 12 Procédé d'analyse de signaux ILS selon la revendication 9, caractérisé en ce que ladite fonction fenêtre est une fonction répondant à l'équation suivante: W(t) = 2 so O (a + ( 1-a) cos 0) avec O = T T étant la durée d'observation, k' = 4 et act = 0,54. 13 Procédé d'analyse de signaux ILS selon la revendication 9, caractérisé en ce que, pour procéder à la détermination des harmoniques de rang r des composantes caractéristiques du signal ILS, le calcul des parties réelle et imaginaire du spectre du signal S*(t) pondéré réalisé à l'étape c) est modifié comme suit: R(k) = R(k) + S(n Te) x cos 2 zkm Te I(k) = I(k) + S(n Te) x sin 2 krn Te 14 Procédé d'analyse de signaux ILS selon l'une quelconque des revendications 9 à 12, caractérisé en ce que lesdites fonctions cos 2 nkn Te et sin 2 x kn Te sont emmagasinées sous la forme de tables de valeurs, N variant de 1 à N, N étant le nombre total d'échantillons prélevés pendant la durée d'observation T.