FR 3080244 A1 20191018 FR 1800298 A 20180411 Procédé et système d'étalonnage de décalage temporel entre voies d'un système de réception multivoies de signaux radioélectriques La présente invention concerne un procédé d'étalonnage de décalage temporel entre voies d'un système de réception multivoies de signaux radioélectriques, et un système d'étalonnage de décalage temporel associé. Elle concerne également un procédé de synchronisation entre voies d'un système de réception multivoies de signaux radioélectriques, et un système de synchronisation associé. L'invention se situe dans le domaine général des systèmes de réception multivoies de signaux radioélectriques, utilisés dans de nombreux domaines tels que la goniométrie d'amplitude, l'interférométrie, les systèmes radar à synthèse d'ouverture multistatique, des mesures de fréquence et l'imagerie médicale. De manière connue, dans un système de réception multivoies, chaque voie de réception comprend une chaîne de réception délivrant un signal électrique en entrée d'un convertisseur analogique numérique ayant une fréquence d'échantillonnage associée, et un module d'interfaçage branché en sortie dudit convertisseur analogique numérique et adapté à fournir les échantillons de signal numérisé en entrée d'un module de traitement numérique du signal par paquets d'échantillons, c'est-à-dire avec un minimum de parallélisation. Cela tient au fait qu'un module de traitement numérique du signal comporte un ou plusieurs composants numériques programmables de type FPGA (Field Programmable Gate Array) dont la fréquence de travail est bien inférieure à la fréquence d'échantillonnage de la conversion analogique numérique. En particulier pour des applications qui exploitent le déphasage entre voies de réception, par exemple pour l'interférométrie, il est nécessaire que les échantillons issus de chaque voie portent sur la même tranche temporelle des signaux radioélectriques captés en entrée des voies de réception, ou, en d'autres termes, que les paquets d'échantillons soient synchronisés entre voies de réception. De plus, il est nécessaire que la synchronisation entre voies de réception soit stable et reproductible. Il existe diverses causes possibles de désynchronisation entre voies de réception. Dans un premier cas de figure, toutes les voies de réception comportent un convertisseur analogique numérique fonctionnant à une fréquence d'échantillonnage Fe. Il s'agit du cas des chaînes de réception respectant la condition de Shannon, effectuant un filtrage de largeur de bande inférieure à Fe/2 avant la conversion numérique-analogique. La fréquence Fe est très élevée et bien supérieure à la fréquence maximum de fonctionnement des composants numériques programmables. On utilise des modules d'interfaçage, par exemple des démultiplexeurs pouvant comporter un ou plusieurs étages 2 de démultiplexage, qui fournissent des paquets de P échantillons en parallèle, P étant un nombre entier non nul, à la fréquence Fe/P. Tout éventuel décalage temporel des démultiplexeurs, par exemple lors du démarrage, induit un décalage possible entre échantillons dans les paquets d'échantillons et les sorties des diverses voies de réception 5 ne sont plus synchrones. Dans ce cas de figure, on connaît dans l'état de la technique des méthodes de synchronisation basées sur l'utilisation d'architectures particulières permettant une distribution de signaux de synchronisation, ayant pour objectif d'effectuer une synchronisation des démultiplexeurs. Des architectures de distribution de signaux de 10 synchronisation en série ou en étoile ont été proposées. De telles solutions sont complexes, et la synchronisation obtenue dépend de la synchronisation d'horloges. La fréquence d'échantillonnage Fe étant très élevée, la période correspondante est très faible, et les temps de propagation et leur dérive ne sont pas négligeables devant cette période. Les systèmes de synchronisation ainsi obtenus peuvent donc devenir instables. 15 Dans un autre cas de figure, les chaînes de réception ne respectent pas la condition de Shannon, effectuant un filtrage de très large bande instantanée, et les convertisseurs analogique-numérique utilisent des fréquences d'échantillonnage Fe, différentes pour résoudre la problématique de repliement de spectre due au non respect de la condition de Shannon. La synchronisation entre voies de réception à base de 20 signaux de synchronisation n'est pas possible dans ce cas. L'invention a pour objectif de remédier aux inconvénients de l'art antérieur, en proposant un procédé et système moins complexe, permettant d'obtenir une synchronisation temporelle robuste et reproductible, et fonctionnant également lorsque les fréquences d'échantillonnage des convertisseurs analogiques numériques du système de 25 réception multivoies sont différentes. A cet effet, l'invention propose un procédé d'étalonnage de décalage temporel entre voies d'un système de réception multivoies de signaux radioélectriques, le système comportant un ensemble de voies de réception de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception délivrant un signal électrique en entrée 30 d'un convertisseur analogique numérique ayant une fréquence d'échantillonnage associée, et un module d'interfaçage branché en sortie dudit convertisseur analogique numérique et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module de traitement numérique de signal comportant un ou plusieurs composants numériques programmables. Ce procédé comporte suite à une mise sous tension du 35 système de réception : 3 - une injection sur toutes les voies de réception d'un même signal aléatoire représentant un bruit obtenu à partir d'une source de bruit, - pour au moins une paire de voies de réception comportant une première voie et une deuxième voie, une estimation d'un décalage temporel entre échantillons de signal 5 numérisés fournis par lesdites première et deuxième voies, en utilisant un calcul d'une fonction d'intercorrélation entre un premier ensemble d'échantillons obtenus par la première voie de réception et un deuxième ensemble d'échantillons obtenus par la deuxième voie de réception, -une mémorisation du décalage temporel estimé en association avec ladite paire 10 de voies. Avantageusement, le procédé de l'invention permet d'estimer un décalage temporel entre voies de réception sans fourniture d'un signal de synchronisation à la fréquence d'échantillonnage entre voies de réception lors d'une mise sous tension du système de réception multivoies. Ainsi, l'architecture du système est simplifiée. 15 Le procédé d'étalonnage de décalage temporel entre voies d'un système de réception multivoies de signaux radioélectriques selon l'invention peut également présenter une ou plusieurs des caractéristiques ci-dessous, prises indépendamment ou en combinaison. Le signal aléatoire représentant un bruit injecté est un bruit analogique thermique 20 injecté en entrée de chaque chaîne de réception. Le signal aléatoire représentant un bruit injecté est un bruit analogique obtenu à partir d'une séquence numérique pseudo-aléatoire et injecté en entrée du convertisseur analogique-numérique de chaque voie de réception. L'estimation dudit décalage temporel comporte une détection d'un maximum de 25 ladite fonction d'intercorrélation, ledit décalage temporel étant un décalage temporel numérique fonction d'un indice temporel correspondant audit maximum. Le procédé est mis en oeuvre dans un système de réception multivoies dans lequel chaque convertisseur analogique numérique a une même fréquence d'échantillonnage, et le calcul d'une fonction d'intercorrélation comporte des étapes de : 30 - acquisition d'un premier ensemble d'échantillons sur ladite première voie de réception et d'un deuxième ensemble de même nombre d'échantillons sur ladite deuxième voie de réception ; - application d'une transformation de Fourier discrète sur ledit premier ensemble d'échantillons, permettant d'obtenir un premier ensemble de valeurs complexes 35 d'échantillons transformés, et application d'une transformation de Fourier discrète sur ledit 4 deuxième ensemble d'échantillons, permettant d'obtenir un deuxième ensemble de valeurs complexes d'échantillons transformés, - calcul du produit complexe entre le premier ensemble de valeurs complexes et le complexe conjugué du deuxième ensemble de valeurs complexes, 5 - obtention de ladite fonction d'intercorrélation par application d'une transformation de Fourier discrète inverse audit produit complexe. Le procédé est mis en oeuvre dans un système de réception multivoies dans lequel chaque convertisseur analogique numérique a une fréquence d'échantillonnage propre associée, et le calcul d'une fonction d'intercorrélation comporte des étapes de : 10 - acquisition d'un premier ensemble d'un premier nombre M1 d'échantillons sur ladite première voie de réception pendant une durée prédéterminée, et acquisition d'un deuxième ensemble d'un deuxième nombre M2 d'échantillons, différent du premier nombre, sur ladite deuxième voie de réception, pendant ladite durée prédéterminée, - application d'une transformation de Fourier discrète sur ledit premier ensemble 15 d'échantillons, permettant d'obtenir un premier ensemble de M1 valeurs complexes d'échantillons transformés, et application d'une transformation de Fourier discrète sur ledit deuxième ensemble d'échantillons, permettant d'obtenir un deuxième ensemble de M2 valeurs complexes d'échantillons transformés, - obtention de deux ensembles de valeurs complexes de même cardinal par 20 suppression des valeurs complexes associées aux fréquences les plus élevées dans l'ensemble de plus grand cardinal, - calcul du produit complexe entre le premier ensemble de valeurs complexes et le complexe conjugué du deuxième ensemble de valeurs complexes de même cardinal, - obtention de ladite fonction d'intercorrélation par application d'une transformation 25 de Fourier discrète inverse audit produit. Le procédé est appliqué dans un système de réception multivoies comportant un nombre N de voies de réception supérieur ou égal à 2, comportant une sélection d'une voie de réception de référence, et l'estimation d'un décalage temporel pour chaque paire de voies composée de ladite voie de réception de référence et d'une deuxième voie parmi 30 les N-1 autres voies de réception. Selon un autre aspect, l'invention concerne un système d'étalonnage de décalage temporel entre voies d'un système de réception multivoies tel que brièvement décrit ci-dessus. Ce système d'étalonnage comporte : - une source de bruit adaptée à injecter un même signal aléatoire représentant un 35 bruit sur toutes les voies de réception, 5 - un module de calcul adapté à estimer, pour au moins une paire de voies de réception comportant une première voie et une deuxième voie, un décalage temporel entre échantillons de signal numérisés fournis par lesdites première et deuxième voies, en utilisant un calcul d'une fonction d'intercorrélation entre un premier ensemble 5 d'échantillons obtenus par la première voie de réception et un deuxième ensemble d'échantillons obtenus par la deuxième voie de réception, - une mémoire adaptée à mémoriser le décalage temporel estimé en association avec ladite paire de voies. Selon un autre aspect, l'invention concerne un procédé de synchronisation entre 10 voies de réception d'un système de réception multivoies de signaux radioélectriques comportant un ensemble de voies de réception de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception délivrant un signal électrique en entrée d'un convertisseur analogique numérique ayant une fréquence d'échantillonnage associée, et un module d'interfaçage branché en sortie dudit convertisseur analogique- 15 numérique et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module de traitement numérique de signal comportant un ou plusieurs composants numériques programmables. Ce procédé de synchronisation comporte, pour au moins une paire de voies de réception : 20 - une phase d'étalonnage de décalage temporel mettant en oeuvre un procédé d'étalonnage tel que brièvement décrit ci-dessus, et - une phase d'application d'une compensation du décalage temporel estimé en sortie du module d'interfaçage d'au moins une des voies de ladite paire de voies par rapport à une voie de référence, de manière à fournir des séries d'échantillons 25 numériques synchronisés en entrée du module de traitement numérique du signal, chaque série d'échantillons numériques synchronisés comportant un nombre d'échantillons correspondant à une même durée et chaque premier échantillon de chaque série correspondant temporellement à un même instant. Selon un autre aspect, l'invention concerne un système de synchronisation entre 30 voies de réception d'un système de réception multivoies de signaux radioélectriques, le système de réception comportant un ensemble de voies de réception de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception délivrant un signal électrique en entrée d'un convertisseur analogique-numérique ayant une fréquence d'échantillonnage associée, et un module d'interfaçage branché en sortie dudit 35 convertisseur analogique numérique et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module de traitement numérique de signal comportant un 6 ou plusieurs composants numériques programmables. Ce système de synchronisation comporte un système d'étalonnage de décalage temporel tel que brièvement décrit ci-dessus, mis en oeuvre dans une phase d'étalonnage de décalage temporel pour au moins une paire de voies de réception, et chaque voie de réception comporte un module de 5 compensation de décalage temporel, branché en sortie d'un module d'interfaçage correspondant, le module de compensation de décalage temporel d'une voie étant adapté à compenser le décalage temporel estimé par rapport à une voie de référence, de manière à fournir des séries d'échantillons numériques synchronisés en entrée du module de traitement numérique du signal, chaque série d'échantillons numériques synchronisés 10 comportant un nombre d'échantillons correspondant à une même durée et chaque premier échantillon de chaque série correspondant temporellement à un même instant. D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux 15 figures annexées, parmi lesquelles : - la figure 1 illustre schématiquement un système de synchronisation entre voies de réception d'un système de réception multivoies selon un premier mode de réalisation ; - la figure 2 illustre schématiquement un système de synchronisation entre voies de réception d'un système de réception multivoies selon un deuxième mode de 20 réalisation; - la figure 3 est un schéma synoptique d'un procédé d'étalonnage de décalage temporel entre voies selon mode de réalisation adapté au cas où les convertisseurs analogiques numériques sont cadencés à une même fréquence Fe ; - la figure 4 est un schéma synoptique des principales étapes d'un procédé de 25 synchronisation de voies selon mode de réalisation adapté au cas où les convertisseurs analogiques numériques sont cadencés à différentes fréquences d'échantillonnage Fe, ; - la figure 5 illustre schématiquement deux ensembles d'échantillons temporels acquis sur une période temporelle donnée ; - la figure 6 est un schéma synoptique des principales étapes d'un mode de 30 réalisation d'un procédé de synchronisation entre voies de réception ; - la figure 7 illustre en détail un mode de réalisation d'un dispositif de synchronisation pour une voie de réception ; - la figure 8 représente schématiquement des paquets d'échantillons issus d'un dispositif d'interfaçage d'une voie de réception ; 7 - les figures 9 et 10 représentent schématiquement un exemple de paquets d'échantillons issus de deux voies de réception, avant et après application de la synchronisation ; - la figure 11 illustre schématiquement des premiers et deuxièmes paquets 5 d'échantillons utilisés pour obtenir une série d'échantillons selon un mode de réalisation ; - les figures 12 et 13 illustrent, dans un exemple, des paquets d'échantillons issus des premier et deuxième ensembles de registres du dispositif de synchronisation, et des séries synchronisées obtenues, et - la figure 14 est un synoptique des principales étapes d'un procédé de 10 synchronisation selon un mode de réalisation. La figure 1 illustre schématiquement un système de synchronisation 1 de voies de réception d'un système de réception multivoies 2 de signaux radioélectriques. Le système de réception multivoies 2 de signaux radioélectriques comprend N 15 voies de réception V1 à VN, l'une quelconque des voies de réception étant référencée ci- après V' avec i compris entre 1 et N. Chacune des voies de réception V1,..., VN comprend une antenne 41,-, 4N, recevant un signal radioélectrique et délivrant un signal électrique à une chaîne de réception (CR) 61 à 6N. 20 Chaque chaîne de réception délivre un signal électrique en entrée d'un convertisseur analogique numérique ou CAN 81 à 8N ayant une fréquence d'échantillonnage associée. Chaque CAN 8, fournit des échantillons numériques en sortie, à une cadence dictée par la fréquence d'échantillonnage Fg, égale à une même fréquence Fe dans un 25 mode de réalisation. Ces échantillons sont fournis à un module d'interfaçage 10, correspondant. Ce module d'interfaçage a pour fonction d'assurer la compatibilité avec un module de traitement numérique de signal 12, associé à la voie de réception et faisant partie d'un module de traitement numérique de signal 12, comportant un ou plusieurs composants numériques programmables de type FPGA, lorsque la fréquence Fe est bien 30 supérieure à la fréquence maximale de fonctionnement du ou des composants numériques programmables. Dans un mode de réalisation, chaque module d'interfaçage 101,..., 10N est un démultiplexeur à un ou plusieurs étages de démultiplexage, fournissant des paquets de P échantillons en parallèle, P étant un nombre entier non nul, à la fréquence Fe/P. 35 Dans un mode de réalisation alternatif, les modules d'interfaçage 101 à 10N sont constitués de liens série rapides. 8 Deux types de systèmes de réception multivoies sont envisagés. Un premier type de système de réception est un système respectant la condition de Shannon, les chaînes de réception réalisant un filtrage de largeur de bande inférieure à Fe/2, chacun des CAN 81 à 8N fonctionnant avec une même fréquence 5 d'échantillonnage Fe. Dans les systèmes de réception de premier type, toutes les voies de réception V1 à VN fournissent chacune un même nombre d'échantillons numériques pour une même durée d'acquisition de signal. Un deuxième type de système de réception est un système ne respectant pas la condition de Shannon, les chaînes de réception réalisant un filtrage de largeur de bande 10 supérieure à Fe,/2, où Fe, est la fréquence d'échantillonnage du CAN 8,. La réalisation d'une synchronisation pour un système de réception de deuxième type est bien plus complexe que dans le cas d'un système de réception du premier type. Dans les systèmes de réception de deuxième type, toutes les voies de réception V1 à VN fournissent des nombres d'échantillons numériques différents pour une même durée d'acquisition de 15 signal. Le système de synchronisation 1 de la figure 1 est adapté pour les systèmes de réception de premier type. Le système de synchronisation 1 enrichit le système de réception multivoies 2 avec un système d'étalonnage de décalage temporel 3 pour réaliser l'estimation de 20 décalages temporels entre les voies de réception et des modules de compensation de décalage temporel pour la synchronisation des voies de réception. Le système d'étalonnage de décalages temporels 3 du système de synchronisation 1 comporte en plus une source de signal aléatoire représentant un bruit 14, branchée en entrée de modules de commutation 161 à 16N. 25 Dans le mode de réalisation de la figure 1, chaque module de commutation 161 à 16N est branché en entrée d'une chaîne de réception correspondante 61 à 6N, et reçoit sur une première entrée, la sortie de l'antenne 41 à 4N de la voie de réception V1 à VN correspondante, et sur une deuxième entrée, le bruit issu de la source signal aléatoire représentant un bruit 14. 30 Chaque module de commutation 161 à 16N est adapté à commuter entre une première position reliée à la première entrée et une deuxième position reliée à la deuxième entrée, pour transmettre soit le signal électrique délivré par l'antenne correspondante, soit le bruit issu de la source 14, en entrée du convertisseur CAN 81 à 8N correspondant. 35 Dans une phase d'étalonnage de décalage temporel, mise en oeuvre lors de la mise sous tension du système, et avant une phase de traitement des signaux 9 radioélectriques reçus, tous les modules de commutation 161 à 16N sont en deuxième position. Dans une phase de fonctionnement opérationnel, tous les modules de commutation 161 à 16N sont en première position. 5 La source de signal aléatoire représentant un bruit 14 est par exemple une source analogique de bruit thermique. La sortie de la source 14 est branchée en entrée de chaque module de commutation 161 à 16N, de préférence en respectant une longueur de câble sensiblement identique pour chaque liaison, les câbles étant de même type, de manière à présenter le 10 même bruit sans décalage temporel en entrée de chaque CAN 81 à 8N. Ainsi un même bruit, de même durée et sans décalage temporel est fourni en entrée de chaque CAN. Le système 1 comprend en outre un processeur 180 et une mémoire 20, qui sont, dans un mode de réalisation, intégrés dans un module de calcul 22 qui comprend également le ou les modules de traitement numérique de signal 12. 15 Le processeur 180 est adapté à mettre en oeuvre des instructions de code logiciel lorsqu'il est mis sous tension, permettant de réaliser des étapes du procédé d'étalonnage de décalage temporel selon l'invention. En phase d'étalonnage, le processeur 180 met en oeuvre un module 190 d'étalonnage de décalage temporel entre échantillons de signal numérisé fournis par une 20 première voie de réception V, et une deuxième voie de réception V2 qui sont distinctes, adapté à calculer un décalage temporel entre voies en utilisant un calcul d'une fonction d'intercorrélation, comme expliqué en détail ci-après. Le module 190 d'étalonnage de décalage temporel entre voies est par exemple mis en oeuvre par des instructions de code logiciel. En variante, le module 190 est réalisé 25 sous forme d'un composant logique programmable tel qu'un FPGA, ou encore sous forme de circuit intégré dédié tel qu'un ASIC (pour Application Specific Integrated Circuit). Les ensembles d'échantillons de signal numérisés 24,...24N sont mémorisés dans la mémoire 20. Chaque ensemble d'échantillons 241...24N comprend un nombre d'échantillons successifs correspondant à une même tranche de temps To 30 Le calcul de décalage temporel est de préférence effectué pour autant de paires de voies que nécessaire pour réaliser un étalonnage de décalage temporel ri,r pour chacune des voies de réception V, par rapport à une voie de référence Vr, comme expliqué plus en détail ci-après. Le processeur 180 est adapté à commander des modules de compensation de 35 décalage temporel 261 à 26N, pour réaliser le décalage temporel en phase opérationnelle, 10 de la valeur de décalage temporel T -1,r à - T N,r préalablement calculé et mémorisé pour chacune des voies de réception V1 à VN. Avantageusement, le système de l'invention permet de réaliser une synchronisation sans fourniture d'un signal de synchronisation entre voies de réception 5 lors d'une mise sous tension du système. La figure 2 illustre schématiquement un deuxième mode de réalisation d'un système de synchronisation 11 de système de réception multivoies. Les références communes avec les éléments de la figure 1 sont conservées. Le système de réception multivoies est analogue à celui décrit en référence à la figure 1. 10 La principale différence avec le premier mode de réalisation est que le système d'étalonnage de décalage temporel 30 comporte une source de signal aléatoire représentatif d'un bruit 29 branchée de manière à fournir un bruit analogique en entrée de chaque CAN 81 à 8N. Dans ce mode de réalisation les modules de commutation 161 à 16N sont insérés 15 entre la sortie de la chaîne de réception 61 à 6N et l'entrée du CAN 81 à 8N correspondant pour chaque voie de réception V1 à VN. En phase d'étalonnage, les modules de commutation 16, sont en deuxième position reliant la deuxième entrée de chaque module, branchée à la source de signal aléatoire représentatif d'un bruit 29, à l'entrée du CAN 8, correspondant. 20 Dans ce mode de réalisation, la source de bruit 29 est par exemple réalisée par un générateur pseudo-aléatoire utilisant par exemple des registres à décalage et un rebouclage par « ou exclusif », ou bien avec un bloc mémoire contenant une séquence numérique pseudo-aléatoire. Dans les deux cas, le signal numérique est transformé en signal analogique par un dispositif d'interfaçage numérique-analogique. Le signal 25 analogique issu de telles sources 29 a des caractéristiques de bruit blanc, à fonction d'autocorrélation de type Dirac. De préférence, la fréquence d'horloge de la séquence pseudo-aléatoire est égale à la fréquence d'échantillonnage Fe lorsque le système de réception 2 est un système de premier type. 30 Lorsque le système de réception 2 est un système de deuxième type, la fréquence d'horloge de la séquence pseudo-aléatoire est par exemple choisie égale à la plus petite fréquence d'échantillonnage Fg. Comme dans le mode de réalisation de la figure 1, le branchement de la source de bruit 29 en entrée des CAN 8, est effectué de manière à s'assurer de l'absence de 35 décalage temporel entre les bruits injectés en entrées des CAN 8, des différentes voies de réception. 11 Le système d'étalonnage de décalage temporel 30 et le système de synchronisation 11 sont adaptés à fonctionner à la fois avec un système de réception multivoies de premier type et avec un système de réception multivoies de deuxième type. La figure 3 est un schéma synoptique des principales étapes d'un procédé 5 d'étalonnage de décalage temporel entre voies de réception mis en oeuvre avec un système de réception multivoies de premier type. Comme expliqué ci-dessus, dans un système de réception multivoies de premier type, tous les convertisseurs analogique-numérique 8, utilisent une même fréquence d'échantillonnage Fe. 10 Les étapes de calcul de ce procédé sont mises en oeuvre par le module de calcul 22, par exemple en partie par un module de traitement numérique de signal de type FPGA, et en partie par un module d'instructions de code logiciel implémentées dans un processeur 180. Le procédé d'étalonnage est mis en oeuvre dans une phase d'étalonnage de 15 décalage temporel lors de la mise sous tension du système, avant son utilisation opérationnelle. Lors d'une première étape le procédé comprend la commande 40 de l'injection d'un bruit analogique, soit en entrée de la chaîne de réception correspondante selon le mode de réalisation de la figure 1, soit en entrée du convertisseur analogique-numérique 20 correspondant selon le mode de réalisation de la figure 2. Lors d'une première étape d'acquisition 42 est effectuée l'acquisition synchrone d'un ou plusieurs premiers ensembles d'échantillons numériques en sortie des modules d'interfaçage 101 à 10N. Le nombre d'échantillons est par exemple M=4096. A l'étape 44, par application d'une transformation de Fourier discrète (TFD), on 25 obtient, pour chaque ensemble d'échantillons acquis, un ensemble de M valeurs complexes d'échantillons transformés. De préférence le nombre M d'échantillons est de la forme M=2q, et une TFD optimisée est appliquée par un composant numérique programmable FPGA 12. De préférence, les TFD appliquées en mode opérationnel sont utilisées. 30 De préférence, les étapes d'acquisition 42 et d'application d'une TFD 44 sont appliquées simultanément pour toutes les voies de réception. Le procédé comprend ensuite une étape 46 de sélection d'une première voie de réception V1 et d'une deuxième voie de réception V2, formant une paire de voies de réception (V1, V2) à traiter. Lorsque N=2, le système comporte une seule paire de voies 35 qui forme la paire de voies à traiter. 12 Optionnellement, une voie de réception donnée est choisie comme voie de référence, qui est réutilisée dans chaque paire de voies à traiter. Pour la première voie de réception V1 on a obtenu et mémorisé un premier ensemble TFD, de M valeurs complexes d'échantillons transformés, notées TFD1,1 à 5 TFIkm. Pour la deuxième voie de réception V2 on a obtenu et mémorisé un deuxième ensemble TFD2 de M valeurs complexes d'échantillons transformés, notées TFD2,1 à TFD2,m. Le procédé comprend ensuite un calcul 48 du produit complexe entre le premier 10 ensemble TFD, de valeurs complexes {TFDia, ...,TFD1,m} et le complexe conjugué TFD2* du deuxième ensemble de valeurs complexes, puis l'application 50 d'une transformation de Fourier discrète inverse sur le résultat du calcul de l'étape 48. On entend ici par calcul du produit entre le premier ensemble de valeurs complexes TFD, et le complexe conjugué TFD2* du deuxième ensemble de valeurs 15 complexes le produit effectué entre valeurs de même rang : TFDtk x TDF",/, pour chaque indice k compris entre 1 et M. Par application de la TFD inverse, on obtient la fonction d'intercorrélation sur M points de la paire de voies de réception traitée. Le procédé comprend ensuite une détermination 52 du maximum de la fonction d'intercorrélation sur M points, en d'autres 20 termes, la détermination de l'indice temporel k correspondant à la valeur maximale de la fonction d'intercorrélation. Le décalage temporel entre les voies V1 et V2 est obtenu par multiplication de la valeur de l'indice temporel correspondant par la période d'échantillonnage 1/Fe ou To/M. L'indice k correspondant au décalage temporel obtenu pour la paire de voies 25 traitée est mémorisé à l'étape 54. En alternative, le décalage temporel obtenu est mémorisé à l'étape 54. Le procédé retourne à l'étape 46 pour sélection d'une nouvelle paire de voies à traiter, et les étapes 48 à 54 sont répétées jusqu'à ce qu'un décalage temporel soit déterminé pour chacune des voies par rapport à au moins une autre voie de réception, le 30 procédé d'étalonnage de décalage prend alors fin. La figure 4 est un schéma synoptique d'étapes d'un procédé d'étalonnage de décalage temporel entre voies de réception mis en oeuvre avec un système de réception multivoies de deuxième type, dans lequel les CAN, ont des fréquences d'échantillonnage Fe, différentes. 13 Le procédé comprend une première étape d'injection de bruit analogique 60, analogue à l'étape 40 précédemment décrite en référence à la figure 3. Optionnellement, une sélection de paire de voies de réception à traiter est également effectuée. Lorsque N=2, le système comporte une seule paire de voies qui 5 forme la paire de voies à traiter. Des étapes 62, 64 d'acquisition d'ensembles d'échantillons sur une première voie V1 et une deuxième voie V2 formant une paire de voies, sur une tranche temporelle d'acquisition To. On acquiert M1 échantillons sur la première voie et M2 échantillons sur la 10 deuxième voie, avec M1 différent de M2, et Mi=Fei.To et M2=Fe2.To. Ainsi, avantageusement, même si le nombre de points acquis respectivement sur la première voie de réception diffère du nombre de points acquis sur la deuxième voie de réception, la résolution spectrale de chaque TFD est égale à 1T1-0. La figure 5 illustre schématiquement un exemple de points correspondant à M1 et 15 M2 échantillons temporels acquis respectivement sur une première voie V1 avec une fréquence d'échantillonnage Fei et une deuxième voie V2 avec une fréquence d'échantillonnage Fe2, avec un nombre M2>Midans cet exemple. De retour à la figure 4, suite à l'acquisition de M1 échantillons sur la première voie V1 et de M2 échantillons sur la deuxième voie V2, on applique une transformation de 20 Fourier discrète sur chacune des voies (étapes 66, 68). La TFD appliquée à l'étape 66 aux M1 échantillons de la première voie permet d'obtenir un premier ensemble de M1 valeurs complexes d'échantillons transformés, noté TFD, : {TFD1,1,...,TFD1 M1}. La TFD appliquée à l'étape 68 aux M2 échantillons de la deuxième voie permet 25 d'obtenir un deuxième ensemble de M2 valeurs complexes d'échantillons transformés, noté TFD2: {TFD1,1, En variante, les TFD des étapes 66 et 68 ne retiennent que la moitié (M1/2 et M2/2) des valeurs complexes des TFD, étant donnée la propriété de symétrie hermitienne du spectre d'un signal réel. 30 Il est à noter que les transformations de Fourier discrètes appliquées aux étapes 66, 68 sont de tailles différentes, et sont de préférence appliquées par calcul logiciel par le processeur 180. Une étape suivante 70 consiste à déterminer le nombre le plus grand d'échantillons fournis, Mmax=Max(Mi, M2), qui est M2 dans l'exemple de la figure 4, et à 35 supprimer dans l'ensemble de valeurs complexes de TFD2, les M2-M1 valeurs correspondant aux fréquences les plus grandes du domaine réel. 14 On obtient un deuxième ensemble de valeurs complexes tronqué, TFD'2, de même cardinal Mr',n=M1 que le premier ensemble de valeurs complexes. Bien entendu, le procédé s'applique de manière analogue si le nombre d'échantillons M1 fournis sur la première voie est supérieur au nombre d'échantillons M2 5 fournis sur la deuxième voie, avec M --max=M1 et Mmin=M2. A l'issue de l'étape 70 on dispose de deux ensembles de valeurs complexes de même cardinal Mmin, égal au nombre le plus petit entre M1 et M2. On applique ensuite à l'étape 72 un calcul du produit complexe entre le premier ensemble de valeurs complexes TFDi et le complexe conjugué TFD'2* du deuxième 10 ensemble de valeurs complexes tronqué, puis une transformation de Fourier discrète inverse (étape 74) sur le résultat du calcul 72. De manière analogue au mode de réalisation de la figure 3, le produit est effectué entre valeurs de même rang : TFD1,k x TDFL1 pour chaque indice k compris entre 1 et Minn- 15 Finalement, une fonction d'intercorrélation 76, sur Mmir,=min(Mi, M2) points, est obtenue. Le procédé se poursuit alors pour rechercher le maximum de la fonction d'intercorrélation, comme décrit ci-dessus en référence à l'étape 52 de la figure 3. La figure 6 est un schéma synoptique des principales étapes d'un procédé de 20 synchronisation entre voies de réception d'un système de réception multivoies, dans un mode de réalisation. Le procédé comprend une étape 80 d'acquisition, en sortie des modules d'interfaçage 10' d'ensembles d'échantillons de signaux numériques fournis par les voies de réception V1 à VN, lorsque les modules de commutation 161 à 16N sont en première 25 position. L'étape 80 est suivie d'une lecture 82 des décalages temporels mémorisés à appliquer à chacune des voies de réception. Ensuite, une compensation de décalage temporel est appliquée à l'étape 84, de manière à fournir en entrée du module de traitement numérique des séries d'échantillons 30 numériques synchronisés temporellement sur une durée dite de synchronisation donnée. Par exemple, la compensation de décalage temporel est effectuée par la mise en oeuvre de dispositifs de synchronisation tels que décrits ci-dessous, adaptés pour fonctionner à la fois avec un système de réception de premier type et un système de réception de deuxième type. 35 En particulier, on décrit ci-dessous en référence aux figures 7 à 15 un mode de réalisation d'un module de compensation de décalage temporel 26i, appelé également 15 dispositif de synchronisation, adapté notamment pour un système de réception de deuxième type, dans lequel les fréquences d'échantillonnage Fei sont différentes. Chaque fréquence d'échantillonnage Fei est bien supérieure à la fréquence maximale de fonctionnement des composants numériques programmables (FPGA) situés 5 en aval des CAN, en particulier dans le module de traitement numérique de signal 12. Le module d'interfaçage 10, a pour fonction de transformer le flux d'échantillons série au rythme Fei en flux série de paquets de P échantillons en parallèle à la première fréquence F* P étant un nombre entier non nul, cette première fréquence étant plus lente et compatible avec la deuxième fréquence ou fréquence de travail des composants 10 numériques programmables F T . On désigne par Si(t) le signal numérique délivré par le module d'interfaçage 10,, issu d'un signal analogique si(t) correspondant. Le module d'interfaçage 10, est, dans ce mode de réalisation, un démultiplexeur à un ou plusieurs étages de démultiplexage, fournissant des paquets de P échantillons en parallèle. Les démultiplexeurs 10, des différentes voies de réception ne sont pas 15 synchronisés, et ne sont pas synchronisables car les fréquences Fei sont différentes, et par conséquent les valeurs sont également différentes. Le dispositif de synchronisation 26;, branché entre le module d'interfaçage 10, et le module de traitement numérique 12' a pour fonction, à partir d'un flux d'entrée série de paquets de P échantillons successifs passés en parallèle, les échantillons étant obtenus à 20 une fréquence d'échantillonnage Fei et ayant subi un possible décalage de temps lié à la voie 1, de former des séries de Ni échantillons successifs, toujours formatées en paquets de P échantillons successifs en parallèle, dont le premier échantillon desdites séries correspond, quel que soit i, au même temps absolu, c'est-à-dire au temps courant dans le flux d'entrée corrigé du décalage temporel à appliquer. 25 Chaque dispositif de synchronisation 26, reçoit en entrée une valeur de décalage temporel ATiy entre la voie V, et la voie de référence Vr, calculé par le procédé décrit ci-dessus. Le système de synchronisation comprend un module de comptage 100, qui est, dans le mode de réalisation préféré un module de comptage commun à tous les 30 dispositifs de synchronisation 26,, qui distribue une même valeur de compteur à un instant donné à chaque dispositif de synchronisation. En variante, le module de comptage 100 est adapté à effectuer une synchronisation des compteurs des dispositifs de synchronisation 26,, par exemple par une remise à zéro initiale, de façon à ce que chacun des compteurs délivre une même 35 valeur de compteur à chaque dispositif de synchronisation 26;. 16 Chaque module de traitement numérique 12, peut, par exemple, réaliser une transformée de Fourier discrète (TFD) sur les Ni échantillons successifs d'une série, correspondant à une durée de signal AT qui est la même pour toutes les voies de réception. Ainsi, quel que soit i : Ni 5 -Fe, = AT (EQ1) En particulier pour des applications interférométriques, il est nécessaire que les séries de Ni échantillons successifs pour des valeurs dei différentes, correspondent toutes au même intervalle de temps absolu pour que le déphasage entre fin de voies, entre deux voies données, soit représentatif de celui en sortie des antennes desdites 10 deux voies sur un même signal incident. La figure 7 illustre en détail un dispositif de synchronisation 26, faisant partie de la voie de réception V. Le dispositif de synchronisation 26, reçoit en entrée des paquets de P échantillons 15 à la première fréquence Fp , et fournit en sorties des séries de Ni échantillons, par paquets de P échantillons. La figure 8 illustre un premier problème à résoudre du fait que le nombre Ni n'est pas, généralement, un multiple de P. La figure 8 illustre la parallélisation du flux d'échantillons en P lignes parallèles, 20 réalisée par le démultiplexeur 10,. Chaque échantillon est symbolisé par une case 140. Ces échantillons sont illustrés dans un diagramme à deux axes temporels TI et T2. Les paquets 150 de P échantillons sont référencés par un indice qi. La durée correspondant à un paquet est de PTei = Les flèches recourbées 160 indiquent la continuité des Fei échantillons dans le temps. 25 Les séries 18 à fournir en sortie sont des séries de Ni échantillons, référencées par un indice k,, montrées en en traits épais. Cependant, le nombre Ni n'étant pas un multiple de P, chaque série 18 ne commencent pas systématiquement par un paquet de P échantillons lui appartenant intégralement, comme illustré à la figure 8 par l'échantillon d'indice 0 de chaque paquet.. La position du premier échantillon d'une série 18 dans un 30 paquet de P échantillons 150 évolue dans le temps, et dépend du reste Ri de la division euclidienne de Ni par P. Le nombre Ni s'écrit sous la forme : Ni-MxP+Ri (EQ 2) 17 avec 0 0 (EQ7) AM(q - 1) = -1 si Ati j = 0 (EQ8) Après initialisation, les formules récurrentes suivantes sont appliquées : AX (q) = Mod (Ent (cn x (P - Ri) - Mod(Ati, P) , P) (EQ9) IAM(q) = 0 si AX(q) AX(q - 1) (EQ10) 15 IAM(q) = 1 si AX(q) REVENDICATIONS1.- Procédé d'étalonnage de décalage temporel entre voies (V1,..., VN) d'un système de réception multivoies (2) de signaux radioélectriques, le système (2) comportant un ensemble de voies de réception VN) de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception (61,.., 6N) délivrant un signal électrique en entrée d'un convertisseur analogique numérique (ah._ 8N) ayant une fréquence d'échantillonnage associée, et un module d'interfaçage 10N) branché en sortie dudit convertisseur analogique numérique (81,..., 8N) et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module (121,...,12N) de traitement numérique de signal comportant un ou plusieurs composants numériques programmables, le procédé étant caractérisé en ce qu'il comporte suite à une mise sous tension du système de réception : - une injection (40, 60) sur toutes les voies de réception d'un même signal aléatoire représentant un bruit obtenu à partir d'une source de bruit (14, 29), - pour au moins une paire de voies de réception comportant une première voie et une deuxième voie, une estimation (52, 76) d'un décalage temporel entre échantillons de signal numérisés fournis par lesdites première et deuxième voies, en utilisant un calcul d'une fonction d'intercorrélation entre un premier ensemble d'échantillons obtenus par la première voie de réception et un deuxième ensemble d'échantillons obtenus par la deuxième voie de réception, - une mémorisation (54) du décalage temporel estimé en association avec ladite paire de voies. 2.- Procédé selon la revendication 1, dans lequel ledit signal aléatoire représentant un bruit injecté est un bruit analogique thermique injecté en entrée de chaque chaîne de réception (61,.., 6N). 3.- Procédé selon la revendication 1, dans lequel ledit signal aléatoire représentant un bruit injecté est un bruit analogique obtenu à partir d'une séquence numérique pseudoaléatoire et injecté en entrée du convertisseur analogique-numérique (81,.., 8N) de chaque voie de réception 22 4.- Procédé selon l'une quelconque des revendications 1 à 3, dans lequel l'estimation dudit décalage temporel comporte une détection d'un maximum de ladite fonction d'intercorrélation, ledit décalage temporel étant un décalage temporel numérique fonction d'un indice temporel correspondant audit maximum. 5 5.- Procédé selon l'une quelconque des revendications 1 à 4, mis en oeuvre dans un système de réception multivoies (2) dans lequel chaque convertisseur analogique numérique (81, ..., 8N) a une même fréquence d'échantillonnage (Fe), dans lequel le calcul d'une fonction d'intercorrélation comporte des étapes de : 10 - acquisition (42) d'un premier ensemble d'échantillons sur ladite première voie de réception et d'un deuxième ensemble de même nombre d'échantillons sur ladite deuxième voie de réception ; - application (44) d'une transformation de Fourier discrète sur ledit premier ensemble d'échantillons, permettant d'obtenir un premier ensemble de valeurs complexes 15 d'échantillons transformés, et application (44) d'une transformation de Fourier discrète sur ledit deuxième ensemble d'échantillons, permettant d'obtenir un deuxième ensemble de valeurs complexes d'échantillons transformés, - calcul (48) du produit complexe entre le premier ensemble de valeurs complexes et le complexe conjugué du deuxième ensemble de valeurs complexes, 20 - obtention (50) de ladite fonction d'intercorrélation par application d'une transformation de Fourier discrète inverse audit produit complexe. 6.- Procédé selon l'une quelconque des revendications 1 à 4, mis en oeuvre dans un système de réception multivoies dans lequel chaque convertisseur analogique 25 numérique (81,..., 8N) a une fréquence d'échantillonnage (Fei,...,FeN) propre associée, dans lequel le calcul d'une fonction d'intercorrélation comporte des étapes de : - acquisition (62) d'un premier ensemble d'un premier nombre M1 d'échantillons sur ladite première voie de réception pendant une durée (T0) prédéterminée, et acquisition (64) d'un deuxième ensemble d'un deuxième nombre M2 d'échantillons, différent du 30 premier nombre, sur ladite deuxième voie de réception, pendant ladite durée prédéterminée, - application (66) d'une transformation de Fourier discrète sur ledit premier ensemble d'échantillons, permettant d'obtenir un premier ensemble de M1 valeurs complexes d'échantillons transformés, et application (68) d'une transformation de Fourier discrète sur ledit deuxième ensemble d'échantillons, permettant d'obtenir un deuxième ensemble de M2 valeurs complexes d'échantillons transformés, 23 - obtention (70) de deux ensembles de valeurs complexes de même cardinal par suppression des valeurs complexes associées aux fréquences les plus élevées dans l'ensemble de plus grand cardinal, - calcul (72) du produit complexe entre le premier ensemble de valeurs complexes 5 et le complexe conjugué du deuxième ensemble de valeurs complexes de même cardinal, - obtention (74) de ladite fonction d'intercorrélation par application d'une transformation de Fourier discrète inverse audit produit. 7.- Procédé selon l'une quelconque des revendications 1 à 6, appliqué dans un 10 système de réception multivoies comportant un nombre N de voies de réception supérieur ou égal à 2, comportant une sélection d'une voie de réception de référence, et l'estimation d'un décalage temporel pour chaque paire de voies composée de ladite voie de réception de référence et d'une deuxième voie parmi les N-1 autres voies de réception. 15 8.- Procédé de synchronisation entre voies d'un système de réception multivoies de signaux radioélectriques, le système comportant un ensemble de voies de réception de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception délivrant un signal électrique en entrée d'un convertisseur analogique numérique ayant une fréquence d'échantillonnage associée, et un module d'interfaçage 20 branché en sortie dudit convertisseur analogique-numérique et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module de traitement numérique de signal comportant un ou plusieurs composants numériques programmables, comportant, pour au moins une paire de voies de réception, une phase 25 d'étalonnage de décalage temporel mettant en oeuvre un procédé conforme à l'une des revendications 1 à 7, et une phase d'application d'une compensation du décalage temporel estimé en sortie du module d'interfaçage d'au moins une des voies de ladite paire de voies par rapport à une voie de référence, de manière à fournir des séries d'échantillons 30 numériques synchronisés en entrée du module de traitement numérique du signal, chaque série d'échantillons numériques synchronisés comportant un nombre d'échantillons correspondant à une même durée et chaque premier échantillon de chaque série correspondant temporellement à un même instant. 35 9.- Système d'étalonnage de décalage temporel entre voies VN) d'un système de réception multivoies (2) de signaux radioélectriques comportant un ensemble 24 de voies de réception (V1,..., VN) de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception (61,.., 6N) délivrant un signal électrique en entrée d'un convertisseur analogique-numérique (81,..., 8N) ayant une fréquence d'échantillonnage associée, et un module d'interfaçage (101,..., 10N) branché en sortie dudit convertisseur 5 analogique-numérique 8N) et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module (12) de traitement numérique de signal comportant un ou plusieurs composants numériques programmables, le système d'étalonnage étant caractérisé en ce qu'il comporte : - une source de bruit (14, 29) adaptée à injecter un même signal aléatoire 10 représentant un bruit sur toutes les voies de réception (V1,..., VN), - un module de calcul (22) adapté à estimer, pour au moins une paire de voies de réception comportant une première voie et une deuxième voie, un décalage temporel entre échantillons de signal numérisés fournis par lesdites première et deuxième voies, en utilisant un calcul d'une fonction d'intercorrélation entre un premier ensemble 15 d'échantillons obtenus par la première voie de réception et un deuxième ensemble d'échantillons obtenus par la deuxième voie de réception, - une mémoire (20) adaptée à mémoriser le décalage temporel estimé en association avec ladite paire de voies. 20 10.- Système de synchronisation entre voies d'un système de réception multivoies de signaux radioélectriques, le système de réception comportant un ensemble de voies de réception de signaux radioélectriques, chaque voie de réception comportant une chaîne de réception délivrant un signal électrique en entrée d'un convertisseur analogique-numérique ayant une fréquence d'échantillonnage associée, et un module 25 d'interfaçage branché en sortie dudit convertisseur analogique numérique et adapté à fournir des ensembles d'échantillons de signal numérisé en entrée d'un module de traitement numérique de signal comportant un ou plusieurs composants numériques programmables, le système de synchronisation comportant un système d'étalonnage de décalage 30 temporel selon la revendication 9, mis en oeuvre dans une phase d'étalonnage de décalage temporel pour au moins une paire de voies de réception, et chaque voie de réception comportant un module de compensation (261, .- - , 26N) de décalage temporel, branché en sortie d'un module d'interfaçage (101,..., 10N) correspondant, le module de compensation de décalage temporel d'une voie étant adapté 35 à compenser le décalage temporel estimé par rapport à une voie de référence, de manière à fournir des séries d'échantillons numériques synchronisés en entrée du module 25 (12) de traitement numérique du signal, chaque série d'échantillons numériques synchronisés comportant un nombre d'échantillons correspondant à une même durée et chaque premier échantillon de chaque série correspondant temporellement à un même instant. 5