FR 3096451 A1 20201127 FR 1905397 A 20190522 Description Titre de l'invention : FILTRAGE PARTICULAIRE ET CENTRALE DE NAVIGATION A CORRELATION DE MESURE Domaine technique [0001] La présente invention concerne un procédé de filtrage particulaire, ainsi qu'une unité de calcul et une centrale de navigation à corrélation de mesure qui mettent en oeuvre un tel procédé. Plus particulièrement, le procédé de filtrage particulaire est du type régularisé à boîtes. Technique antérieure [0002] La fonction de navigation d'un aéronef inclut l'estimation de sa position instantanée, sa vitesse instantanée et son orientation instantanée, appelée attitude, dans le référentiel de navigation, aussi appelé trièdre géographique local. L'ensemble des valeurs instantanées des coordonnées de position, c'est-à-dire la latitude, la longitude et l'altitude de l'aéronef, de vitesse, comprenant une coordonnée de vitesse en direction du Nord, une coordonnée de vitesse en direction de l'Est et une coordonnée de vitesse de descente, et des angles d'attitude de l'aéronef, comprenant un angle de roulis, un angle de tangage et un angle de lacet, constitue l'état instantané du système formé par l'aéronef. L'aéronef peut être un avion, un drone ou tout porteur aérien autopropulsé, sans limitation. [0003] L'accélération et la vitesse angulaire de l'aéronef sont mesurées répétitivement et à haute cadence chacune selon trois axes, par exemple à une fréquence de répétition d'environ 1000 Hz (hertz), à l'aide d'accéléromètres et de gyromètres d'une centrale inertielle qui est à bord de l'aéronef. La centrale de navigation délivre alors des estimations des coordonnées de position, de vitesses et des angles d'attitude de l'aéronef, par intégration des résultats des mesures accélérométriques et gyrométriques au cours du temps. Toutefois, les mesures d'accélération et de vitesse angulaire sont affectées d'une erreur, qui est principalement composée d'un biais, d'un facteur d'échelle et d'un processus aléatoire, et l'accumulation des erreurs de mesure résulte en une erreur qui affecte l'estimation de l'état instantané de l'aéronef Cette erreur sur l'état instantané qui est estimé croît en fonction du temps, et est couramment appelée dérive inertielle. Il est alors nécessaire d'associer la centrale inertielle à au moins un capteur additionnel, afin de corriger ou réduire la dérive inertielle. [0004] Des méthodes usuelles pour corriger ou réduire la dérive inertielle consistent à utiliser des signaux de géolocalisation, tels que des signaux de navigation par satellite, par exemple de type GNSS, pour «global navigation satellite system» en anglais, ou d'utiliser des signaux qui sont produits par des balises situées au sol, par exemple des 2 signaux de radionavigation ou de type GBAS, pour «ground-based augmentation system» ou système de renforcement au sol, ou encore de recevoir par radio une localisation de l'aéronef qui a été effectuée à l'aide d'un radar. Mais de telles méthodes sont sensibles à des brouillages ou à des leurres, à la disponibilité d'une couverture de la zone où se trouve l'aéronef en signaux de localisation ou en communication, etc. Il est alors souhaitable dans certaines circonstances de disposer à bord de l'aéronef d'une méthode qui soit autonome pour corriger et/ou réduire la dérive inertielle. Pour cela, il est courant d'associer à la centrale inertielle une sonde télémétrique, qui mesure une distance entre l'aéronef et le sol. Cette sonde télémétrique peut être un radio-altimètre, un télémètre laser, etc. sans limitation. Elle mesure la distance entre l'aéronef et le sol selon une direction qui peut être fixe ou non par rapport à l'aéronef. Lorsque cette direction de mesure de distance peut être variable, son orientation par rapport à l'aéronef est connue. La centrale inertielle est alors aussi associée à une unité de calcul qui corrèle des résultats de mesures successives qui sont effectuées par la sonde télé-métrique avec l'état instantané de l'aéronef tel qu'estimé par la centrale inertielle. Plus précisément, une caractérisation de la zone qui est survolée par l'aéronef est stockée à bord de l'aéronef, par exemple sous forme d'une carte de relief qui associe une valeur de hauteur de relief à chaque couple de valeurs de latitude et de longitude. Un tel enregistrement de carte de relief peut être sous la forme d'une table, avec la latitude et la longitude qui constituent les entrées de la table, et les valeurs de hauteur du relief qui constituent les réponses de lecture dans la table. Alternativement, la caractérisation de la zone survolée peut être stockée sous forme d'une fonction analytique qui permet de calculer les valeurs de hauteur de relief en fonction des valeurs de latitude ou de longitude, ou sous toute autre forme appropriée. Alors, à chaque nouvelle estimation de la position instantanée de l'aéronef qui est produite par la centrale inertielle, une valeur de la distance qui devrait exister entre l'aéronef et le sol est obtenue en interrogeant la caractérisation de la zone survolée telle que stockée à bord de l'aéronef, conformément à la position estimée de l'aéronef. Eventuellement, la valeur de distance aéronef-sol peut résulter d'un calcul qui combine la position estimée de l'aéronef avec le(s) résultat(s) d'une ou plusieurs interrogation(s) de la caractérisation de la zone survolée qui est stockée à bord de l'aéronef, notamment lorsque la direction de mesure de la sonde télémétrique est oblique par rapport à l'axe d'altitude. Un tel calcul est connu de l'Homme du métier, si bien qu'il n'est pas utile de le répéter ici. La valeur de distance aéronef-sol qui est ainsi estimée est ensuite comparée avec le résultat de mesure qui est délivré par la sonde télémétrique. Une telle méthode est couramment appelée navigation avec corrélation de terrain. Plusieurs variantes en ont été proposées, mais certaines d'entre elles présentent une grande sensibilité à l'existence de non-linéarités des profils de terrain. Autrement dit, elles présentent l'inconvénient d'un 3 manque de robustesse dans leur efficacité de convergence vers l'état vrai de l'aéronef, en fonction des profils de terrain possibles. [0005] Les méthodes de navigation à corrélation de terrain qui sont basées sur un filtre par- ticulairc régularisé à boîtes, ou BRPF pour «Box Regularized Particic Filter» en anglais, permettent de corriger la dérive inertielle avec une robustesse plus grande, compatible avec l'existence de non-linéarités et d'ambiguïtés de terrain. L'article de Merlinge, N., Dahia, K.. Piet-Lahanier, H., Bruscy, J., & Horri, N. qui est intitulé «A Box Regularized Particic Filter for state estimation with severely ambiguous and nonlincar nicasurements», Automatica (2019), Vol. 104, pp. 102-110, décrit une telle méthode. Chacune de ces méthodes consiste encore à calculer itérativement un état instantané de l'aéronef à partir d'un dernier état déterminé antérieurement, en utilisant les mesures qui sont effectuées par la centrale inertielle, et en corrélant le résultat obtenu avec la mesure de distance qui est effectuée par la sonde télémétrique, la mesure de distance étant reconstituée avec la caractérisation du relief terrestre qui est embarquée à bord de l'aéronef Mais, un filtre particulairc régularisé à boîtes procède en manipulant des intervalles d'état, de dimension neuf lorsque chaque état de l'aéronef est constitué par trois coordonnées de position, trois coordonnées de vitesse et trois angles d'attitude comme décrit plus haut. De plus, une valeur de poids à signification probabiliste est associée à chaque intervalle d'état : le poids de chaque intervalle d'état correspond à la probabilité pour que l'état vrai de l'aéronef se situe dans cet intervalle d'état. Toutefois, ces méthodes de navigation avec corrélation de terrain qui sont basées sur des filtres particulaires n'ont pas encore été mises en oeuvre pour des aéronefs réels, à cause des ressources de calcul très importantes qui sont nécessaires. En effet, il est requis pour de nombreuses applications aéronautiques que le procédé de navigation à corrélation de terrain qui est utilisé puisse être mis en oeuvre par une unité de calcul de type circuit à réseau de portes programmables, ou FGPA pour «Field-Programmable Gate Array» en anglais. Or les circuits de ce type ont des capacités qui sont encore trop limitées. [0006] De façon connue, un procédé de filtrage particulaire régularisé à boîtes comprend au moins une exécution d'une séquence qui inclut les étapes suivantes : /1/ une étape de prédiction, comprenant de prédire des intervalles d'état ultérieurs, chaque intervalle d'état ultérieur étant obtenu en appliquant au moins une règle de propagation à un d'une pluralité d'intervalles d'état antérieurs, chaque intervalle d'état étant possiblement multidimensionnel ; /2/ une étape de mesure de l'état vrai ; /3/ une étape de contraction de l'un au moins des intervalles d'état ultérieurs, en fonction d'au moins un résultat de mesure de l'état vrai qui a été obtenu à l'étape /2/ ; /4/ une étape de mise-à-jour de poids, comprenant d'affecter un poids à chaque in- tervallc d'état ultérieur en fonction d'une taille de cet intervalle d'état ultérieur tel que résultant de l'étape /3/, d'une taille du même intervalle d'état ultérieur tel que résultant de l'étape /1/ avant l'étape /3/, et d'un poids de l'intervalle d'état antérieur dont a résulté cet intervalle d'état ultérieur lors de l'étape /1/ ; et /5/ une étape de redistribution des intervalles d'état ultérieurs, comprenant de remplacer au moins un des intervalles d'état ultérieurs i par ni sous-intervalles provenant d'une division de l'intervalle d'état ultérieur i, i étant un indice entier de numérotation des intervalles d'état ultérieurs, et ni étant un nombre des sous-intervalles qui remplacent cet intervalle d'état ultérieur i, chaque sous-intervalle étant destiné à constituer un intervalle d'état antérieur pour une itération ultérieure de la séquence des étapes /1/ à /5/. [0007] L'étape /5/ est couramment exécutée en utilisant une méthode de tirage multinomial. Une telle méthode consiste à tirer des points aléatoirement et répétitivement à l'intérieur d'un segment unidimensionnel, pour déterminer le nombre ni des sous-intervalles qui sont destinés à remplacer l'intervalle d'état ultérieur i, en fonction des poids, notés wi, de tous les intervalles d'état ultérieurs tels que mis-à-jour à l'étape /4/. Le segment de tirage est constitué par une juxtaposition de segments élémentaires, dont les longueurs individuelles correspondent une-à-un aux poids des intervalles d'état ultérieurs. Le nombre ni des sous-intervalles qui remplaceront l'intervalle d'état ultérieur i est alors proportionnel au nombre de points tirés aléatoirement qui appartiennent au segment élémentaire dont la longueur est égale au poids de l'intervalle d'état ultérieur i. Si la valeur de ni qui est ainsi déterminée pour l'un des intervalles d'état ultérieurs est nulle, cet intervalle d'état ultérieur est supprimé, sinon l'intervalle d'état ultérieur i est divisé en ni sous-intervalles, la valeur obtenue pour ni étant statistiquement d'autant plus grande que le poids de l'intervalle d'état ultérieur i tel que résultant de l'étape /4/ est élevé. Mais une telle méthode de tirage multinornial ne peut pas être effectuée de façon parallélisée entre tous les intervalles d'état ultérieurs, de sorte qu'elle génère un goulet d'étranglement dans l'exécution de la séquence des étapes /1/ à /5/. En outre, elle ne permet pas de garantir une durée d'exécution pour l'étape /5/ qui soit constante. En effet, cette méthode de tirage multinomiale requiert de comparer chacun d'une série de nombres tirés aléatoirement, correspondant aux abscisses des points introduits précédemment, aux sommes progressives des poids Si \vi, dans l'ordre des valeurs croissantes de i jusqu'à identifier la première de ces sommes progressives qui est supérieure au nombre aléatoire. Par principe, le nombre de tests de comparaison de chaque nombre aléatoire à une des sommes progressives des poids n'est pas prévisible. Or, notamment pour une correction de dérive de la centrale inertielle qui est effectuée par corrélation de terrain, il est essentiel que la durée d'exécution de l'étape /5/ soit maîtrisée, pour garantir la périodicité d'actualisation de l'estimation de l'état de l'aéronef en temps réel. Problème technique 100081 A partir de cette situation, un but de la présente invention est d'optimiser l'utilisation d'un filtre particulaire régularisé à boîtes, notamment pour que ce filtre puisse être mis en oeuvre à bord d'un aéronef. Au-delà, l'invention a pour but de contribuer à ce qu'un procédé de filtrage particulaire régularisé à boîtes puisse être mis en oeuvre par un circuit de type PRIA. [0009] Un autre but dc l'invention est de permettre l'utilisation d'un filtre particulaire ré- gularisé à boîtes, en assurant que la durée dc calcul qui est nécessaire pour chaque itération d'estimation de l'état instantané du système soit inférieure à une durée maximale contrôlée. En particulier, l'invention a pour but de permettre de paralléliser au moins une partie de l'étape redistribution des intervalles d'état ultérieurs, qui est particulièrement exigeante en ressources de calcul. Résumé de l'invention [0010] Pour atteindre l'un au moins de ces buts ou un autre, un premier aspect de l'invention propose un procédé de filtrage particulaire à boîtes pour prédire un état d'un système, qui comprend les étapes /1/ à /5/ mentionnées ci-dessus, mais dans lequel l'étape /5/ est modifiée. Selon l'invention, à la place de la méthode de tirage multinomial, l'étape /5/ comprend d'initialiser le nombre de sous-intervalles ni, et comprend la séquence des sous-étapes suivantes /5-1/ à /5-4/ qui est effectuée pour chacun des intervalles d'état ultérieurs utilisé alternativement en tant qu'intervalle de référence : /5-1/ adopter initialement l'intervalle de référence comme intervalle d'état retenu ; /5-2/ obtenir un nombre qui est compris entre 0 et 1, appelé facteur aléatoire, par une méthode de tirage considérée comme aléatoire, et sélectionner un des intervalles d'état ultérieurs, différent ou non de l'intervalle de référence, en tant qu'intervalle de comparaison par une méthode de sélection considérée comme aléatoire ; et /5-3/ si le résultat d'une multiplication du facteur aléatoire par le poids de l'intervalle de référence, tel que mis-à-jour à l'étape /4/, est inférieur au poids de l'intervalle de comparaison, aussi tel que mis-à-jour à l'étape /4/, alors remplacer l'intervalle d'état retenu antérieurement par l'intervalle de comparaison ; les sous-étapes /5-2/ et /5-3/ étant répétées plusieurs fois pour chaque intervalle d'état ultérieur utilisé en tant qu'intervalle de référence, puis /5-4/ augmenter d'une unité le nombre des sous-intervalles destinés à remplacer l'intervalle d'état qui est retenu àl'issue de la répétition des sous-étapes /5-2/ et /5-3/. Chaque intervalle d'état ultérieur est ensuite remplacé lors d'une sous-étape /5-5/ par des sous-intervalles qui proviennent d'une division dc cet intervalle d'état ultérieur, conformément au nombre de sous-intervalles qui résulte pour cc même intervalle d'état 6 ultérieur de toutes les exécutions de la sous-étape /5-4/. [0011] Un tel mode d'exécution de l'étape /5/ peut facilement être parallélisé, car il ne procède que par des tests de comparaison qui utilisent les poids d'intervalles d'état pris par couples. En outre, cc mode d'exécution de l'étape /5/ garantit une durée d'exécution qui est déterministe, notamment parce que le nombre de tests de comparaison à effectuer ne dépend pas de la valeur de chaque facteur aléatoire qui est tiré à la sous-étape /5-2/. [0012] Dans le cadre de l'invention, on entend par méthode de tirage/sélection considérée comme aléatoire une méthode supposée fournir des résultats de tirage/sélection qui sont conformes ou sensiblement conformes à une loi de probabilité uniforme. De telles méthodes sont connues, y compris des méthodes couramment dites pseudo-aléatoires. [0013] Dans diverses mises en oeuvre de l'invention, les préférences et perfectionnements suivants peuvent avantageusement être utilisés, indépendamment les uns des autres ou en combinaison de plusieurs d'entre eux : - les sous-étapes /5-2/ et /5-3/ peuvent être répétées un nombre de fois qui est compris entre 1 et 100 pour chaque intervalle d'état ultérieur utilisé en tant qu'intervalle de référence ; - à l'étape /5/, le nombre ni des sous-intervalles qui proviennent d'une division de l'intervalle d'état ultérieur i peut être initialisé systématiquement à zéro avant d'effectuer la séquence des sous-étapes /5-1/ à /5-4/ pour cet intervalle d'état ultérieur i utilisé en tant qu'intervalle de référence. Alternativement, une initialisation à la valeur zéro du nombre ni peut être appliquée si le poids de l'intervalle d'état ultérieur i tel que résultant de l'étape /4/ est nul, sillon une initialisation du nombre ni à la valeur un peut être appliquée. Une telle initialisation différenciée des nombres ni assure que chaque intervalle d'état i soit conservé lors de l'étape de redistribution /5/ si son poids tel que résultant de l'étape /4/ est non-nul ; - l'étape /5/ peut n'être exécutée que si un critère de représentativité n'est pas satisfait par les poids des intervalles d'état ultérieurs tels que mis-à-jour à l'étape /4/. En particulier, un critère de représentativité d'un type dit N-efficace ou d'un type en-tropique peut être utilisé ; - à l'étape /4/, le poids mis-à-jour de chaque intervalle d'état ultérieur peut être égal au poids de l'intervalle d'état antérieur dont a résulté cet intervalle d'état ultérieur lors de l'étape /1/, multiplié par la taille de l'intervalle d'état ultérieur tel que résultant de l'étape /3/, et divisé par la taille de l'intervalle d'état ultérieur tel que résultant de l'étape /1/ avant d'appliquer l'étape /3/ à cet intervalle d'état ultérieur ; - l'étape /5/ peut comprendre en outre la sous-étape supplémentaire suivante, exécutée après que chaque intervalle d'état ultérieur a été remplacé par des sous-intervalles obtenus par division de cet intervalle d'état ultérieur : 7 /5-6/ une sous-étape de lissage par noyau, lors de laquelle l'un au moins des sous-intervalles et un poids qui est associé à ce sous-intervalle sont modifiés selon un noyau probabiliste pour former un intervalle d'état antérieur qui est destiné à être utilisé pour l'itération ultérieure de la séquence des étapes!]! à /5/ ; - à chaque exécution de la sous-étape /5-2/, le facteur aléatoire peut être obtenu par une méthode de type LFSR, pour «lincar feedback shift register» en anglais ou registre à décalage avec rétroaction linéaire. Une telle méthode LFSR peut aussi être utilisée pour sélectionner aléatoirement chaque intervalle de comparaison. Des méthodes LFSR peuvent notamment être utilisées pour générer un nombre pseudo-aléatoire selon un noyau uniforme, notamment une méthode LFSR de Galois ; - les sous étapes /5-1/ à /5-5/, et la sous-étape /5-6/ le cas échéant, peuvent être effectuées de façon parallélisée pour tous les intervalles d'état ultérieurs utilisés chacun en tant qu'intervalle de référence ; et - le système peut être un aéronef, l'état du système comprenant alors des coordonnées de position, de vitesse et d'attitude de l'aéronef en cours de vol, et chaque mesure de l'état vrai du système délivrant une valeur d'une distance qui existe entre l'aéronef et un point du sol. [0014] Un deuxième aspect de l'invention propose une unité de calcul, cette unité de calcul comprenant au moins une première entrée qui est adaptée pour recevoir des résultats de mesures répétées d'accélération et de vitesse angulaire d'un système, et une seconde entrée qui est adaptée pour recevoir des résultats de mesures répétées de l'état vrai du système, supplémentaires pai- rapport aux mesures d'accélération et de vitesse angulaire. L'unité de calcul est alors agencée pour exécuter un procédé de filtrage particulaire régularisé à boîtes qui est conforme au premier aspect de l'invention. De cette façon, l'unité de calcul produit en sortie une série d'intervalles d'état avec des poids respectifs, le poids qui est associé à chacun des intervalles d'état correspondant à une valeur de probabilité pour que l'état vrai du système soit dans cet intervalle d'état. [0015] Une telle unité de calcul peut être du type circuit à réseau de portes programmables, ou FPGA, du type circuit à réseau de portes fixes, connu sous l'acronyme DSP pour «Digital Signal Processor» en anglais, ou du type processeur d'unité centrale, connu sous les acronymes CPU pour «Computer Processing Unit» ou RISC pour «Reduced Instruction set Computer». [0016] L'unité de calcul peut comprendre plusieurs sous-unités de calcul qui sont destinées à être dédiées une-à-une aux intervalles d'état ultérieurs, et adaptées pour exécuter de façon parallélisée les sous-étapes /5-1/ à /5-5/, et le cas échéant aussi la sous-étape / 5-6/, pour tous les intervalles d'état ultérieurs. [0017] Un troisième aspect de l'invention propose une centrale de navigation à corrélation de mesure, qui est adaptée pour être embarquée à bord d'un véhicule et qui comprend : 8 - une centrale inertielle, adaptée pour mesurer itérativement des accélérations et des vitesses angulaires du véhicule, et pour déduire en utilisant des résultats de mesures des accélérations et des vitesses angulaires, des intervalles d'état ultérieurs respectivement à partir de plusieurs intervalles d'état antérieurs, chaque état du véhicule comprenant des coordonnées de position, de vitesse et d'attitude de ce véhicule ; - un système de mesure, adapté pour mesurer itérativement au moins une caractéristique d'un état vrai du véhicule ; et - une unité de calcul qui est conforme au deuxième aspect de l'invention, et qui est adaptée pour réduire au moins une dérive de position, de vitesse et/ou d'attitude de la centrale inertielle, en utilisant les résultats de mesure de la caractéristique de l'état vrai du véhicule qui sont délivrés par le système de mesure. [0018] Enfin, un quatrième aspect de l'invention propose un véhicule qui comprend une centrale de navigation à corrélation de mesure conforme au troisième aspect de l'invention. Un tel véhicule peut être un aéronef, notamment un avion, un drone volant ou tout porteur aérien autopropulsé, ou un véhicule capable de se déplacer sur un sol, notamment un dronc mobile sur sol, un navire, un sous-marin, un engin spatial, notamment une sonde, un satellite, etc. sans limitation. Scion le cas, le système de mesure peut être une sonde télémétrique destinée à mesurer une distance entre le véhicule et un point de sol, un système de repérage d'amères ou d'étoiles, un sonar pour mesurer une hauteur d'eau sous le navire ou le sous-marin, etc. Brève description des figures [0019] Les caractéristiques et avantages de la présente invention apparaîtront plus clairement dans la description détaillée ci-après d'exemples de mise en oeuvre non-limitatifs, en référence aux figures annexées parmi lesquelles : [0020] [fig.1[ montre un aéronef qui est équipé d'une centrale inertielle recalée par cor- rélation de terrain conforme à l'invention ; et [0021] [fig.21 est un diagramme qui montre l'enchaînement d'étapes d'un procédé conforme à l'invention, avec une exécution parallélisée des étapes. Description détaillée de l'invention 100221 Pour raison de clarté, les dimensions des éléments qui sont représentés symbo- liquement dans 1fig.11 ne correspondent ni à des dimensions réelles, ni à des rapports de dimensions réels. En outre, l'invention est décrite à titre d'exemple non-limitatif dans un cas d'application à un aéronef, mais il est entendu qu'elle peut être appliquée à tout véhicule qui est muni d'une centrale de navigation à corrélation de mesure, que ce véhicule soit terrestre, aérien, maritime, spatial, etc. 100231 Conformément àlfig.11, un aéronef 20 est équipé d'une centrale de navigation à cor- rélation de terrain, désignée par la référence 10. La centrale de navigation 10 comprend 9 une centrale inertielle 1, une sonde télémétrique 2 et une unité de calcul 3. De façon connue, la centrale inertielle 1 réalise répétitivement des mesures de trois coordonnées d'accélération et de trois coordonnées de vitesse angulaire de l'aéronef 20, à l'aide d'accéléromètres et de gyromètres non-représentés. Par ailleurs, la sonde télémétrique 2 réalise répétitivement des mesures de la distance H qui existe entre le sol 100 et l'aéronef 20. Cette distance H est mesurée selon une direction qui peut être fixe par rapport à l'aéronef 20, par exemple perpendiculairement à un plan de référence de l'aéronef Eventuellement, la direction de mesure télémétrique peut être variable par rapport à l'aéronef 20, mais dans un tel cas cette direction de mesure est contrôlée et prise en compte d'une façon adéquate, connue par ailleurs. La distance H varie donc en fonction du relief terrestre qui est survolé par l'aéronef 20, ainsi qu'en fonction de son altitude et de son attitude. Bien que la description soit limitée dans la suite au seul cas de valeurs de la distance H comme résultats de mesures supplémentaires par rapport à celles de la centrale inertielle 1, correspondant au cas d'une navigation avec corrélation de terrain, il est entendu que d'autres mesures supplémentaires peuvent être utilisées alternativement à celles de la distance H, ou en complément de celles-ci. Enfin, l'invention est compatible avec des modèles de centrale inertielle et de sonde télé-métrique tels que disponibles commercialement. En particulier, la centrale inertielle peut être d'un type à MEMS, pour «Micro-Electro-Mechanical System» en anglais, à quadrason, à gyrolascr, etc., et la sonde télémétrique peut être un radio-altimètre, un télémètre laser, etc. [0024] Chaque état possible pour l'aéronef 20 peut être composé de trois valeurs de co- ordonnées spatiales qui identifient une position pour l'aéronef, par exemple dans le référentiel terrestre, trois valeurs de vitesses chacune selon une des coordonnées spatiales, et trois valeurs angulaires pour identifier une orientation de l'aéronef, par exemple trois valeurs d'angles d'Euler, pour un total de neuf coordonnées d'état. Dans ces conditions, un intervalle d'état pour l'aéronef 20 est formé par une combinaison de neuf intervalles unidimensionnels qui sont relatifs chacun séparément à l'une des coordonnées d'état. Un tel intervalle d'état est appelé boîte dans le jargon de l'Homme du métier. [0025] Le procédé est initialisé par la fourniture de plusieurs intervalles d'état initiaux, chacun associé à un poids qui indique une valeur de probabilité pour que l'état vrai de l'aéronef 20 se trouve initialement dans cet intervalle d'état initial. Ainsi, N intervalles d'état initiaux sont fournis, N étant un nombre entier, par exemple égal à 16 ou 32, de préférence inférieur ou égal à 128. Chaque intervalle d'état initial est associé individuellement à une valeur de poids, qui peut être égale à 1/N. [0026] Le procédé est ensuite constitué par des itérations successives d'une séquence d'étapes, chaque nouvelle exécution de la séquence d'étapes produisant une actua- lisation des intervalles d'état, avec des valeurs actualisées de poids qui sont associées une-à-un aux intervalles d'état actualisés. En outre, chaque nouvelle exécution de la séquence d'étapes est effectuée à partir des intervalles d'état et leurs valeurs de poids associées tels que fournis par l'exécution juste antérieure de la séquence d'étapes. [0027] Chaque séquence d'étapes comprend une étape de prédiction, notée /1/ dans [fig.2], une étape de mesure notée /2/, une étape de contraction notée /3/, une étape /4/ de mise-à-jour du poids de chaque intervalle d'état, et une étape /5/ de redistribution des intervalles d'état. Les étapes /1/ et /3/ à /5/ sont exécutées pour chaque intervalle d'état, de préférence de façon parallélisée pour tous les intervalles d'état. Ainsi, au sein de l'unité de calcul 3, une première sous-unité de calcul, référencée 3-1 dans [Fig. 1] peut être dédiée à exécuter les étapes /1/ et /3/ à /5/ pour un premier des intervalles d'état, noté boîte _1 dans [Fig. 2], une deuxième sous-unité de calcul 3-2 peut être dédiée à exécuter les étapes /1/ et /3/ à /5/ pour un deuxième des intervalles d'état, noté boîte 2, une Nième sous-unité de calcul 3-N peut être dédiée à exécuter les étapes /1/ et / 3/ à /5/ pour le Nième des intervalles d'état, noté boîte N. Etant donné que l'étape /5/ produit une redistribution des intervalles d'état tels que résultant de l'étape /3/, elle est de préférence exécutée de façon à conserver un nombre constant d'intervalles d'état. Alors l'unité de calcul 3 peut être conçue et dimensionnée pour comprendre les N sous-unités de calcul 3-1, 3-2,..., 3-N. [0028] Pour chaque intervalle d'état, l'étape /1/, dite de prédiction, consiste à recueillir les résultats de dernières mesures d'accélération et de vitesse angulaire, tels que délivrés par la centrale inertielle 1. Eventuellement, les résultats de plusieurs mesures qui ont été effectuées par la centrale inertielle 1 depuis l'exécution précédente de la séquence des étapes /1/ à /5/ peuvent être recueillis. Pour chacun des intervalles d'état, l'étape /1/ comprend aussi de calculer une évolution de cet intervalle d'état pendant la période de temps qui s'est écoulée entre les deux exécutions de l'étape /1/, relatives à l'itération précédente et la nouvelle itération en cours de la séquence des étapes /1/ à /5/. Le principe de calcul d'une telle évolution de chaque intervalle d'état, à partir des résultats des mesures d'accélération et de vitesse angulaire, est supposé connu de l'Homme du métier. On pourra se référer à ce sujet à l'article de Merlinge, N., Dahia, K., PietLahanier, H., Brusey, J., & Horri, N. qui est intitulé «A Box Regularized Particle Filter for state estimation with severely ambiguous and non-linear measurements», Automatica (2019), Vol. 104, pp. 102-110. Cette étape /1/ aboutit à une translation, accompagnée le plus souvent d'un changement de longueur, de chaque intervalle unidimensionnel qui est relatif à l'une des coordonnées d'état de l'aéronef 20. Dans la partie générale de la présente description, chaque intervalle d'état tel qu'existant au moment de débuter l'exécution de l'étape /1/ a été appelé intervalle d'état antérieur, et cet intervalle d'état tel que modifié par l'étape /1/ a été appelé intervalle d'état ultérieur. 11 [0029] L'étape /2/, dite de mesure, consiste à recueillir le résultat de la dernière mesure de distance H tel que délivré par la sonde télémétrique 2. Eventuellement, les résultats de plusieurs dernières mesures qui ont été effectuées par la sonde télémétrique 2 depuis l'exécution précédente de l'étape /2/, lors de l'itération précédente de la séquence des étapes /1/ à /5/, peuvent être recueillis. [0030] L'étape /3/, dite de contraction, consiste à réduire la taille de chaque intervalle d'état en fonction d'incompatibilités qui existeraient entre des parties de cet intervalle d'état et le ou les dernier(s) résultat(s) de mesure de distance H recueillis à l'étape /2/. Typiquement, les six coordonnées de position et d'attitude de l'aéronef 20, telles que pouvant varier à l'intérieur de chaque intervalle d'état tel que re-déterminé à l'étape / 1/, sont combinées à des valeurs de hauteur de relief lues dans une carte de terrain qui est mémorisée, afin d'obtenir une estimation de la distance H associée à chaque état. Un tel calcul peut être effectué de la façon connue qui a été rappelée au début de cette description. La distance H ainsi estimée est comparée au résultat de la mesure de l'étape /2/. De façon générale, une équation implicite peut être mise en oeuvre pour convertir chaque état de l'aéronef 20 en une valeur de distance H, en utilisant la carte de terrain mémorisée. Mais une telle équation, appelée équation d'observation par l'Homme du métier, peut être difficile à inverser localement, si bien que modéliser sa fonction inverse par une partie de fonction analytique et une partie de fonction tabulée peut être avantageuse. L'Homme du métier pourra encore se reporter à l'article de Merlinge et al., Automatica 104(2019), cité plus haut, au sujet de l'étape /3/ qui n'est pas l'objet de l'amélioration du procédé que procure l'invention. [0031] L'étape /4/, dite de mise-à-jour, consiste à mettre-à-jour le poids de chaque intervalle d'état tel que résultant de l'étape /3/. Par exemple, une nouvelle valeur du poids de l'intervalle d'état boîte_l peut être égale à la valeur du poids de l'intervalle d'état boîte_l telle qu'existant avant d'effectuer cette mise-à-jour, multipliée par la taille de l'intervalle d'état boîte_l tel que résultant de l'étape de contraction /3/, et divisée par la taille de l'intervalle d'état boîte_l tel que résultant de l'étape de prédiction /1/ avant d'appliquer l'étape de contraction /3/. Des calculs analogues de mise-à-jour de poids peuvent être effectués séparément pour l'intervalle d'état boîte_2,..., boîte_N. D'autres formules de mise-à-jour des poids peuvent être adoptées alternativement. Eventuellement, les valeurs de poids telles que résultant de l'une de ces formules peuvent en plus être con-igées, par exemple en les multipliant par un facteur commun non nul pour assurer que leur somme soit égale à l'unité. [0032] Les étapes de prédiction /1/, de contraction /3/ et de mise-à-jour des poids /4/ qui viennent d'être décrites sont avantageusement effectuées de façon parallélisée pour les intervalles d'état boîte_1, boîte_2,..., boîte_N par les sous-unités de calcul 3-1, 3-2,...,3-N. Un gain de temps de calcul en résulte. 12 [0033] [fig.2]L'étape /5/ a pour but de redistribuer les intervalles d'état tels que résultant de l'étape /3/, pour obtenir une meilleure représentativité statistique des états qui sont possibles pour l'aéronef 20. L'exécution de cette étape /5/ peut être soumise au résultat d'un test optionnel désigné par CR dans . Ce test consiste à déterminer si un critère de représentativité est satisfait par les intervalles d'état avec leurs valeurs de poids respectives. Il porte sur tous les poids wi tels qu'actualisés à l'étape /4/, où wi est le poids qui est associé à l'intervalle d'état boîte_i, i étant un indice entier variant dc 1 à N. Un premier critère de représentativité qui peut être utilisé est celui connu sous l'appellation de critère N-efficace. Ce critère est satisfait si ( y N 1M-1 42.. ) Oera - N' °`'` 1. d'ajustement du critère entropique, compris entre 0 et 1. Mais d'autres critères de représentativité qui sont aussi connus de l'Homme du métier peuvent encore être utilisés alternativement. L'étape /5/ de redistribution des intervalles d'état est alors appliquée si le critère de représentativité n'est pas satisfait. [0034] L'étape /5/ comprend d'abord de déterminer un nombre de division ni qui est attribué à l'intervalle d'état boîte_i, avec l'indice i qui numérote encore les intervalles d'état de 1 à N. Elle débute alors par une initialisation à zéro de tous les nombres de division ni. Cette sous-étape d'initialisation est notée /5-0/ dans ifig. 21. Les sous-étapes /5-1/ à / 5-4/ à exécuter ensuite à l'intérieur de l'étape /5/ aboutissent à la détermination des nombres de division ni. Puis l'étape /5/ comprend les sous-étapes /5-5/ et /5-6/ pour déterminer ni sous-intervalles d'état qui remplaceront l'intervalle d'état boîte_i pour une nouvelle exécution complète de la séquence des étapes /1/ à /5/. [0035] La méthode de détermination des nombres de division ni est décrite maintenant pour l'intervalle d'état générique boîte_i. Elle est basée sur les poids de tous les intervalles d'état boîte_i, boîte_2,..., boîte_N, tels que mis-à-jour à l'étape /4/. indique cette méthode sous forme de pseudocode pour 1=1,1=2 et La valeur du poids wi est comparée à celles des poids respectifs de NB intervalles d'état qui sont sélectionnés de façon aléatoire ou pseudo-aléatoire, en étant multipliée à chaque fois par un facteur aléatoire ou pseudo-aléatoire noté u et appartenant à l'intervalle [0 ; 1]. L'intervalle d'état i a été appelé intervalle de référence dans la partie générale de la présente description, et chacun des NB intervalles d'état sélectionnés de façon aléatoire ou pseudoaléatoire a été appelé intervalle de comparaison, avec son poids qui est noté wj de façon générique, j étant un indice de numérotation des NB intervalles de comparaison. 13 Les intervalles de comparaison sont sélectionnés séparément et indépendamment pour chaque intervalle de référence boîte i. Alors chaque fois que pour un des intervalles de comparaison boîte_j l'inégalité suivante est vérifiée, U - Wi REVENDICATIONSProcédé de filtrage particulairc régularisé à boîtes, pour prédire un état d'un système, comprenant au moins une exécution d'une séquence qui inclut les étapes suivantes : /1/ une étape de prédiction, comprenant de prédire des intervalles d'état ultérieurs, chaque intervalle d'état ultérieur étant obtenu en appliquant au moins une règle de propagation à un d'une pluralité d'intervalles d'état antérieurs ; /2/ une étape de mesure d'un état vrai du système ; /3/ une étape de contraction de l'un au moins des intervalles d'état ultérieurs, en fonction d'au moins un résultat de mesure de l'état vrai qui a été obtenu à l'étape /2/ ; /4/ une étape de mise-à-jour de poids, comprenant d'affecter un poids à chaque intervalle d'état ultérieur en fonction d'une taille dudit intervalle d'état ultérieur tel que résultant de l'étape /3/, d'une taille dudit intervalle d'état ultérieur tel que résultant de l'étape /1/ avant l'étape /3/, et d'un poids de l'intervalle d'état antérieur dont a résulté ledit intervalle d'état ultérieur lors de l'étape /1/ ; et /5/ une étape de redistribution des intervalles d'état ultérieurs, comprenant de remplacer au moins un des intervalles d'état ultérieurs i par ni sous-intervalles provenant d'une division de l'intervalle d'état ultérieur i, i étant un indice entier de numérotation des intervalles d'état ultérieurs et ni étant un nombre des sous-intervalles qui remplacent ledit intervalle d'état ultérieur i, chaque sous-intervalle étant destiné à constituer un intervalle d'état antérieur pour une itération ultérieure de la séquence des étapes /1/ à /5/, le procédé étant caractérisé en ce que l'étape /5/ comprend d'initialiser le nombre de sous-intervalles ni, et comprend la séquence des sous-étapes suivantes /5-1/ à /5-4/ qui est effectuée pour chacun des intervalles d'état ultérieurs utilisé alternativement en tant qu'intervalle de référence : /5-1/ adopter initialement l'intervalle de référence comme intervalle d'état retenu ; /5-2/ obtenir un nombre qui est compris entre 0 et 1, appelé facteur aléatoire, par une méthode de tirage considérée comme aléatoire, et sélectionner un des intervalles d'état ultérieurs, différent ou non de l'intervalle de référence, en tant qu'intervalle de comparaison par une 17 [Revendication 2] [Revendication 3] [Revendication 4] [Revendication 5] méthode de sélection considérée comme aléatoire ; /5-3/ si un résultat d'une multiplication du facteur aléatoire par le poids de l'intervalle de référence, tel que mis-à-jour à l'étape /4, est inférieur au poids dc l'intervalle dc comparaison, aussi tel que mis-à-jour à l'étape /4/, alors remplacer l'intervalle d'état retenu antérieurement par l'intervalle dc comparaison ; les sous-étapes /5-2/ et /5-3/ étant répétées plusieurs fois pour chaque intervalle d'état ultérieur utilisé en tant qu'intervalle de référence, puis /5-4/ augmenter d'une unité le nombre des sous-intervalles destinés à remplacer l'intervalle d'état qui est retenu à l'issue de la répétition des sous-étapes /5-2/ et /5-3/ ; et chaque intervalle d'état ultérieur étant ensuite remplacé lors d'une sous-étape /5-5/ par des sous-intervalles provenant d'une division dudit intervalle d'état ultérieur, conformément au nombre de sous-intervalles qui résulte pour ledit intervalle d'état ultérieur de toutes les exécutions de la sous-étape /5-4/. Procédé selon la revendication 1, suivant lequel les sous-étapes /5-2/ et / 5-3/ sont répétées un nombre de fois qui est compris entre 1 et 100 pour chaque intervalle d'état ultérieur utilisé en tant qu'intervalle de référence. Procédé selon la revendication 1 ou 2, suivant lequel l'étape /5/ n'est exécutée que si un critère de représentativité n'est pas satisfait par les poids des intervalles d'état ultérieurs tels que mis-à-jour à l'étape /4/. Procédé selon l'une quelconque des revendications précédentes, suivant lequel, à l'étape /4/, le poids mis-à-jour de chaque intervalle d'état ultérieur est égal au poids de l'intervalle d'état antérieur dont a résulté ledit intervalle d'état ultérieur lors de l'étape /1/, multiplié par la taille de l'intervalle d'état ultérieur tel que résultant de l'étape /3/, et divisé par la taille de l'intervalle d'état ultérieur tel que résultant de l'étape /1/ avant d'appliquer l'étape /3/ audit intervalle d'état ultérieur. Procédé selon l'une quelconque des revendications précédentes, suivant lequel l'étape /5/ comprend en outre la sous-étape supplémentaire suivante, exécutée après que chaque intervalle d'état ultérieur a été remplacé par des sous-intervalles obtenus par division dudit intervalle d'état ultérieur : /5-6/ une sous-étape de lissage par noyau, lors de laquelle un au moins des sous-intervalles et un poids associé audit sous-intervalle sont modifiés selon un noyau probabiliste pour former un intervalle d'état 18 [Revendication 6] [Revendication 7] [Revendication 8] [Revendication 9] [Revendication 10] [Revendication 11] [Revendication 12] antérieur qui est destiné à être utilisé pour l'itération ultérieure de la séquence des étapes /1/ à /5/. Procédé selon l'une quelconque des revendications précédentes, suivant lequel, à chaque exécution de la sous-étape /5-2/, le facteur aléatoire est obtenu par une méthode de type à registre à décalage avec rétroaction linéaire, et/ou l'intervalle de comparaison est sélectionné une méthode de type à registre à décalage avec rétroaction linéaire. Procédé selon l'une quelconque des revendications précédentes, suivant lequel les sous étapes /5-1/ à /5-5/ sont effectuées de façon parallélisée pour tous les intervalles d'état ultérieurs utilisés chacun en tant qu'intervalle de référence. Procédé selon l'une quelconque des revendications précédentes, suivant lequel le système est un aéronef (20), l'état du système comprend des coordonnées de position, de vitesse et d'attitude de r aéronef en cours de vol, et suivant lequel chaque mesure de l'état vrai du système délivre une valeur d'une distance qui existe entre l'aéronef et un point de sol. Unité de calcul (3), comprenant au moins une première entrée adaptée pour recevoir des résultats de mesures répétées d'accélération et de vitesse angulaire d'un système, une seconde entrée adaptée pour recevoir des résultats de mesures répétées d'un état vrai du système, supplémentaires par rapport aux mesures d'accélération et de vitesse angulaire, et l'unité de calcul étant agencée pour exécuter un procédé de filtrage particulaire régularisé à boîtes qui est conforme à l'une quelconque des revendications précédentes, de façon à produire en sortie une série d'intervalles d'état avec des poids respectifs, le poids qui est associé à chacun des intervalles d'état correspondant à une valeur de probabilité pour que l'état vrai du système soit dans ledit intervalle d'état. Unité de calcul (3) selon la revendication 9, de type circuit à réseau de portes programmables, circuit à réseau de portes fixes, ou processeur d'unité centrale. Unité de calcul (3) selon la revendication 9 ou 10, comprenant plusieurs sous-unités de calcul (3-1, 3-2,..., 3-N) qui sont destinées à être dédiées une-à-une aux intervalles d'état ultérieurs, et adaptées pour exécuter de façon parallélisée au moins les sous-étapes /5-1/ à /5-5/, pour tous les intervalles d'état ultérieurs. Centrale de navigation à corrélation de mesure (10), adaptée pour être embarquée à bord d'un véhicule, comprenant : 19 [Revendication 13] - une centrale inertielle (1), adaptée pour mesurer itérativement des accélérations et des vitesses angulaires du véhicule, et pour déduire en utilisant des résultats de mesures des accélérations et des vitesses angulaires, des intervalles d'état ultérieurs respectivement à partir de plusieurs intervalles d'état antérieurs, chaque état du véhicule comprenant des coordonnées de position, de vitesse et d'attitude dudit véhicule ; - un système de mesure (2), adapté pour mesurer itérativement au moins une caractéristique d'un état vrai du véhicule ; et - une unité de calcul (3) conforme à l'une quelconque des revendication 9 à 11, et adaptée pour réduire au moins une dérive de la centrale inertielle (1), parmi une dérive de position, une dérive de vitesse et une dérive d'attitude, en utilisant des résultats de mesure de la caractéristique de l'état vrai du véhicule qui sont délivrés par le système de mesure (2). Véhicule, comprenant une centrale de navigation à corrélation de mesure (10) qui est conforme à la revendication 12.