TITRE : Procédé, dispositif et programme d'ordinateur d'estimation d'une vitesse d'un véhicule à roue

L’invention concerne un procédé, un dispositif et un programme d'ordinateur pour l'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol.

Le domaine d'application de l'invention est notamment les véhicules automobiles terrestres roulants, les aéronefs à train d'atterrissage ayant des roues roulant sur le sol au décollage ou à l'atterrissage.

La vitesse instantanée d'un véhicule est une donnée très importante dans le domaine de la navigation inertielle. Dans l'état de la technique, le signal de vitesse est généralement mesuré par un capteur supplémentaire de vitesse (odomètre) placé sur la roue du véhicule. L'installation d'un tel capteur présente un coût et n'est pas toujours envisageable dans certaines applications, particulièrement en aéronautique et en navigation terrestre militaire.

L'invention cherche à se dispenser d'un capteur direct de vitesse tel qu'un odomètre ou d'un capteur de position.

Un objectif de l'invention est de développer un procédé, un dispositif et un programme d'ordinateur pour l'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol, qui soient capables d'estimer avec précision la vitesse instantanée d'un véhicule à roue à partir d'au moins un capteur inertiel ou vibratoire (pouvant être par exemple un accéléromètre ou gyroscope ou microphone), en se dispensant d'un capteur direct de vitesse tel qu'un odomètre ou d'un capteur de position.

Dans l'état de la technique, des tentatives ont été faites pour surpasser le problème d'absence d'un capteur direct de vitesse en exploitant les informations présentes dans les signaux inertiels, comme par exemple dans le document « The Future of Automotive Localization Algorithms: Available, reliable, and scalable localization: Anywhere and anytime » de Rickard Karlsson et Fredrik Gustafsson, 2017. Selon ce document, l'estimation de la vitesse instantanée à partir d'un signal inertiel est calculée en sélectionnant à chaque instant la fréquence la plus énergétique d'un spectrogramme. Ce type d'approche peut être efficacement applicable dans certaines applications simples où le bruit est faible, mais pas dans le domaine d'application mentionné ci-dessus, où le rapport bruit/signal est souvent très fort. Le procédé décrit par ce document de Rickard Karlsson et Fredrik Gustafsson n'est pas robuste au bruit environnemental tel que la combustion d'un moteur et nécessite un spectre dans lequel les pics sont facilement identifiables, ce qui n'est pas possible en présence d'un moteur à combustion. De plus, le spectre du signal inertiel ou vibratoire est riche en harmoniques, et ces derniers subissent des fuites d'énergie considérables ce qui rend l'identification des pics par des transformées de Fourrier délicate. La grande majorité des systèmes de navigation inertielle de l'état de la technique ne sont pas utilisables dans un bruit environnemental prépondérant.

Dans certaines applications, notamment les applications aéronautiques et les véhicules terrestres militaires, où le nombre de capteurs est très restreint, l'installation d'un capteur de position est très difficile, voire impossible en aéronautique. La présence d'un signal de vitesse instantanée est de grande importance pour un système de navigation afin d'estimer la position d'un porteur de manière très précise. Pour cela, il y a un besoin de mesurer cette vitesse instantanée sans avoir besoin d'ajouter une instrumentation matérielle sur les roues du véhicule.

L'invention vise à obtenir un procédé, un dispositif et un programme d'ordinateur pour l'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol, qui pallient les inconvénients mentionnés ci-dessus et atteignent les objectifs mentionnés ci-dessus.

A cet effet, un premier objet de l'invention est un procédé d'estimation d'une vitesse de déplacement d'un véhicule ayant au moins une roue de roulement sur le sol, caractérisé en ce que on mesure un signal inertiel ou vibratoire à partir d'au moins un capteur inertiel ou vibratoire du véhicule, on calcule par au moins un calculateur une pluralité de spectres du signal sur une pluralité de fenêtres temporelles successives définies par rapport à des instants successifs, la pluralité des spectres formant au moins un spectrogramme en fonction des instants successifs, on filtre le spectrogramme par au moins un filtre d'atténuation de bruit et/ou d'extraction de caractéristiques déterministes, pour obtenir au moins un spectrogramme filtré du signal sur les instants successifs, on estime par le calculateur une trajectoire de fréquence, représentative de la vitesse de la roue du véhicule dans le spectrogramme filtré sur les instants successifs, par le fait que : on estime par le calculateur une probabilité d'observation de la trajectoire de fréquence à partir d'une amplitude calculée du spectrogramme filtré pour chaque instant, on estime par le calculateur une loi a posteriori d'observation de la trajectoire de fréquence comme étant proportionnelle au produit d'une probabilité de transition prescrite de la trajectoire de fréquence et de la probabilité d'observation de la trajectoire de fréquence, on estime par le calculateur la trajectoire de fréquence pour les instants successifs à partir de la loi a posteriori d'observation de la trajectoire de fréquence, on estime par le calculateur la vitesse de déplacement de la roue du véhicule pour les instants successifs à partir de la trajectoire de fréquence, ayant été estimée.

L'invention propose ainsi une méthode d'estimation du signal inertiel ou vibratoire, qui est robuste au bruit et qui utilise l'information de fréquence présente dans le signal (inertiel ou vibratoire) mesuré, et une méthodologie basée sur une approche probabiliste spécifiquement dédiée. L'invention, à la différence de l'état de la technique, permet de s'affranchir du bruit environnemental et d'extraire un signal de vitesse du véhicule sans utiliser d'odomètre et permet ainsi l'exploitation de la donnée dans des cas réels. Elle s'exprime en fonction de l'amplitude en chaque point du plan temps-fréquence et d'une loi de transition permettant de favoriser la continuité, la direction et la régularité de la trajectoire tout en exploitant la cinématique du système comme les harmoniques. Un signal inertiel est toujours présent dans le cas où un système de navigation du véhicule est présent. Dans le cas du capteur de vibrations, la fréquence des vibrations émises par les roues du véhicule est proportionnelle à la vitesse du véhicule (il suffit de multiplier par le périmètre de la roue).

Suivant un mode de réalisation de l'invention, l'on détermine par le calculateur un support de valeurs de fréquence, on estime par le calculateur la trajectoire de fréquence comme étant pour chaque instant dans le support la valeur de fréquence, qui maximise la loi a posteriori d'observation de la trajectoire de fréquence jusqu'à cet instant.

Suivant un autre mode de réalisation de l'invention, l'on détermine par le calculateur un support de valeurs de fréquence, on estime pour chaque instant par le calculateur la fréquence comme étant un barycentre de la loi a posteriori d'observation de la fréquence, prise en les valeurs de fréquence du support et affectée de ces valeurs de fréquence.

Suivant un mode de réalisation de l'invention, l'on détermine par le calculateur un support Ω de valeurs de fréquence, l'on estime par le calculateur la probabilité d'observation sachant une trajectoire comme dépendant de l'amplitude calculée du spectrogramme filtré et normalisé pour chaque instant. Suivant un mode de réalisation de l'invention, l'on estime par le calculateur la probabilité p _Zt ǀ _ξt d'observation selon l'équation suivante : où |S(f _t ,t)| est l'amplitude calculée du spectrogramme filtré pour chaque instant t, f _min (t) est une trajectoire minimale de fréquence ayant été mesurée, f _max (t) est une trajectoire maximale de fréquence ayant été mesurée, f _min (t) et f _max (t) SC telles que f _min (t) - f _max (t),

1 _{[ fmin(t), fmax(t)]} (ft) est une fonction prenant la valeur 1 pour f _t entre f _min (t) f _max (t) et prenant la valeur 0 pour f _t en-dessous de f _min (t) ou au-dessus de f _max (t)

Suivant un mode de réalisation de l'invention, l'on estime par le calculateur la probabilité p _Zt ǀ _ξt d'observation selon l'équation suivante :

OÙ y est uunn eennttiieerr oouu rrééeell positif, vvééririffiiaanntt où f _t désigne les fréquences de chaque spectrogramme filtré S(f _t , t) à l'instant t, et où F _s est une fréquence de seuil, ayant été prédéterminée en fonction du spectrogramme filtré S(f _t , t).

Suivant un mode de réalisation de l'invention, l'on détermine par le calculateur un support de valeurs de fréquence, l'on estime par le calculateur la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence selon les équations suivantes : pour ξ _t = f _t ∈ Ω avec dépendant de l'amplitude calculée du spectrogramme filtré pour chaque instant t, où Z _t est un vecteur binaire d'observation ayant une taille N _s prescrite pour l'instant t ∈ [1, ...,N _s ] tel que

Suivant un mode de réalisation de l'invention, où |S(f _t ,t)| est l'amplitude calculée du spectrogramme S(f _t , t) filtré pour chaque instant t, f _min (t) est une trajectoire minimale de fréquence ayant été mesurée, f _max (t) est une trajectoire maximale de fréquence ayant été mesurée, f _min (t) et sont telles que

1 _{[ fmin(t), fmax(t)]} est une fonction prenant la valeur 1 pour f _t entre et f _max (t) et prenant la valeur 0 pour f _t en-dessous de ou au-dessus de f _max (t).

Suivant un mode de réalisation de l'invention, l'on détermine par le calculateur un support de valeurs de fréquence, on prescrit par le calculateur la probabilité p _ξ (f ₁ , ... , f _NS ) de transition de la trajectoire ξ _t de fréquence pour l'instant N _s en partant de la probabilité initiale donnée p _ξ1 de l'instant 1, pour t ∈ {1, ...,N _s } et pour chaque fréquence f _t ∈ Ω selon l'équation suivante : où N _s est un nombre entier prescrit, supérieur ou égal à deux, d'instants prescrits successifs t, est la probabilité de transition de la fréquence ξ _t de la trajectoire de fréquence de l'instant présent t à la fréquence ξ _{t +1} de la trajectoire de fréquence de l'instant suivant t+1.

Suivant un mode de réalisation de l'invention, la probabilité de transition de la fréquence ξ _t de la trajectoire de fréquence de l'instant présent t à la fréquence ξ _{t +1} de la trajectoire de fréquence de l'instant suivant t+1 suit une loi conditionnelle gaussienne centrée sur la fréquence f _t du spectrogramme filtré à l'instant t.

Suivant un mode de réalisation de l'invention, la loi conditionnelle gaussienne centrée sur la fréquence f _t a un écart-type σ _t = Δt.y _t où Δt est une résolution temporelle calculée du spectrogramme de l'instant t, Yt est une vitesse moyenne calculée avec laquelle la trajectoire ξ _t de fréquence évolue entre l'instant t et l'instant suivant t+1.

Suivant un mode de réalisation de l'invention, l'on estime par le calculateur la vitesse de déplacement de la roue du véhicule en multipliant la trajectoire de fréquence par 2nR, où R est le rayon de la roue et est prédéterminé ou estimé.

Suivant un mode de réalisation de l'invention, le filtre d'atténuation de bruit et/ou d'extraction de caractéristiques déterministes comprend un filtre médian glissant et/ou un filtre de séparation de caractéristiques déterministes par rapport à des caractéristiques aléatoires. Un deuxième objet de l'invention est un dispositif d'estimation d'une vitesse de déplacement d'un véhicule ayant au moins une roue de roulement sur le sol, caractérisé en ce que le dispositif comporte : un capteur inertiel ou vibratoire du véhicule pour mesurer un signal inertiel ou vibratoire, un calculateur configuré pour calculer une pluralité de spectres du signal sur une pluralité de fenêtres temporelles successives définies par rapport à des instants successifs, la pluralité des spectres formant au moins un spectrogramme en fonction des instants successifs, un filtre d'atténuation de bruit et/ou d'extraction de caractéristiques déterministes dans le spectrogramme, pour obtenir au moins un spectrogramme filtré du signal sur les instants successifs, le calculateur étant configuré pour estimer une trajectoire de fréquence, représentative de la vitesse de la roue du véhicule dans le spectrogramme filtré sur les instants successifs, par: estimation, par le calculateur, d'une probabilité d'observation de la trajectoire de fréquence à partir d'une amplitude calculée du spectrogramme filtré pour chaque instant, estimation, par le calculateur, d'une loi a posteriori d'observation de la trajectoire de fréquence comme étant proportionnelle au produit d'une probabilité de transition prescrite de la trajectoire de fréquence et de la probabilité d'observation de la trajectoire de fréquence, estimation, par le calculateur, de la trajectoire de fréquence pour les instants successifs à partir de la loi a posteriori d'observation de la trajectoire de fréquence, le calculateur étant configuré pour estimer la vitesse de déplacement de la roue du véhicule pour les instants successifs à partir de la trajectoire de fréquence.

Un troisième objet de l’invention est un programme d’ordinateur, comportant des instructions de code pour la mise en œuvre du procédé d’estimation de la vitesse de déplacement du véhicule tel que décrit ci-dessus, lorsqu’il est exécuté dans un calculateur.

L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif en référence aux figures ci-dessous des dessins annexés.

[Fig. 1] représente un exemple d’organigramme du déroulement du procédé d’estimation d’une vitesse de déplacement d’un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l’invention. [Fig. 2] représente schématiquement un exemple d'étape de calcul de spectrogrammes d'un signal du procédé, du dispositif et du programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 3] représente un exemple de spectrogramme d'un signal temporel dans le plan temps-fréquence du procédé, du dispositif et du programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 4] représente schématiquement un exemple de dispositif d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 5] représente un exemple de signal brut pouvant être mesuré par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 6] représente un spectrogramme du signal de la figure 5, pouvant être calculé par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 7] représente un exemple de spectrogramme filtré du signal de la figure 5, pouvant être calculé par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 8] représente un exemple de probabilité d'observation de la trajectoire de fréquence, pouvant être calculée à partir du signal de la figure 5 par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 9] représente un exemple de probabilité de transition de la trajectoire de fréquence, pouvant être calculée à partir du signal de la figure 5 par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

[Fig. 10] représente un exemple de trajectoire de fréquence, pouvant être calculée à partir du signal de la figure 5 par le procédé, le dispositif et le programme d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention. [Fig. 11] représente schématiquement un exemple de dispositif d'estimation d'une vitesse de déplacement d'un véhicule à roue de roulement sur le sol suivant des modes de réalisation de l'invention.

On décrit ci-dessous plus en détail en référence aux figures 1 à 10 des modes de réalisation et des exemples de réalisation du procédé, du dispositif 1 et du programme d'ordinateur pour l'estimation d'une vitesse v instantanée de déplacement d'un véhicule VE à roue RO de roulement sur le sol.

Dans le domaine d'application de l'invention, mentionné ci-dessus, dans le cas où des signaux vibratoires issus de systèmes mécaniques complexes tels que les véhicules terrestres ou les systèmes aéronautiques sont utilisés, ces signaux sont fortement affectés par du bruit (les bruits ambiants, de combustion et aérodynamique, etc), d'autant plus pour les sources rotatives interférantes et asynchrones. Ainsi, l’identification de la bonne trajectoire représentant la vitesse v de déplacement du véhicule VE dans le plan temps- fréquence est particulièrement difficile dans des conditions de faible rapport signal sur bruit : quelques pics liés à la fréquence instantanée peuvent être discriminés et de faux pics dus au bruit peuvent en prendre la place localement. Cette première difficulté peut être résolue en ajoutant une étape E4 de prétraitement de la représentation temps-fréquence par filtrage du bruit pour améliorer l'émergence de la signature des roues par rapport au bruit de fond. La méthode E4 de dé-bruitage proposée permet non seulement d'atténuer le bruit de fond mais aussi d'atténuer les composantes quasi-stationnaires dans le plan temps- fréquence dans le cas où la composante recherchée est fortement non-stationnaire (montée de régime, descente de régime) et inversement d'atténuer les composantes non- stationnaires dans le cas où la composante recherchée est quasi-stationnaire.

Ensuite, la nature multi-composante du signal et la possibilité d'existence de composantes asynchrones dans le moteur comme dans un moteur d'aéronef, fait que dans la même bande de fréquence, plusieurs harmoniques peuvent coexister, se croiser, être très proches ou même se confondre. Ainsi, en suivant simplement les maximums à chaque instant, on peut dévier de la bonne trajectoire et donc fournir une estimation erronée de la vitesse v. Pour pallier ce problème, l'invention prévoit dans la procédure d'estimation, en plus de l'amplitude du spectrogramme, des contraintes supplémentaires sur la trajectoire optimale se traduisant par l'estimation d'une probabilité de transition. La prise en compte de plusieurs harmoniques et de la direction de la trajectoire d'intérêt permet par exemple de donner plus de poids à des trajectoires plus proches de celle de la fréquence instantanée recherchée ce qui est avantageux pour les moments de couplage entres différentes composantes. L'ajout de ces contraintes de régularité par l'estimation d'une probabilité de transition permet également d'éviter les sauts et les discontinuités en présence de pics isolés liés au bruit résiduel après l'étape E4 de prétraitement.

L'objectif est la recherche, parmi toutes les trajectoires possibles dans les spectrogrammes au cours du temps, d'une trajectoire optimale ξ _t ou £ de fréquence, vérifiant les contraintes considérées et représentative de la vitesse v de la roue RO du véhicule VE. L'optimalité d'une trajectoire par rapport à une autre est mesurée par une fonction de coût ou d'une façon équivalente par une densité de probabilité. Pour ce faire, les informations du spectrogramme et les contraintes sont exprimées sous forme de probabilités. Plus particulièrement, la trajectoire ξ _t ou ξ ₁ ,..,ξ _Ns de fréquence peut être modélisée par une chaîne de Markov avec une densité de transition pξ _t+1 | ξ _t et une probabilité initiale donnée p _ξ1 . La densité de transition est définie en fonction des contraintes ajoutées (continuité, harmoniques, direction, etc.). A chaque pas de temps t, une observation du spectrogramme est fournie avec une probabilité p _Zt ǀ _ξt qui dépend de l'amplitude du spectrogramme. La trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE peut être choisie comme étant celle qui maximise la loi a posteriori Pξ ₁ ,..,ξ _Ns |z ₁ ,...,z _Ns (ou en calculant sa moyenne). On peut également calculer les incertitudes sur cette estimation en s'intéressant aux moments d'ordres 2. Le nombre d'harmoniques à prendre en compte peut être entré par l'utilisateur sur l'interface physique INT1 d'entrée de données. Aux figures, le procédé d'estimation d'une vitesse v de déplacement d'un véhicule ayant une (ou plusieurs) roue RO de roulement sur le sol comporte un première étape E1 de mesure d'un signal temporel X, qui est inertiel ou vibratoire à partir d'un (ou plusieurs) capteur CAP inertiel ou vibratoire du véhicule. Ce capteur CAP inertiel ou vibratoire peut être situé en dehors du moteur du véhicule. Le capteur CAP inertiel ou vibratoire peut être présent sur la roue RO de roulement au sol du véhicule VE. Du fait d'un environnement extrême, un capteur est souvent sujet à des défauts ou à des erreurs de mesure. Le (ou les) capteur CAP peut être ou comprendre un (ou plusieurs) accéléromètre et/ou un (ou plusieurs) gyroscope et/ou un (ou plusieurs) capteur acoustique et/ou un (ou plusieurs) microphone. L'invention peut utiliser des capteurs inertiels qui ne sont pas situés dans le moteur. De plus, dans le cas où le véhicule VE comporte une centrale inertielle composée de 6 capteurs inertiels (3 accéléromètres et 3 gyroscopes), l'invention peut utiliser la redondance de capteurs CAP formés par plusieurs de ou tous ces 6 capteurs inertiels pour améliorer la précision de l'information de vitesse v (le risque d'erreur de mesure simultanée sur 6 capteurs étant plus faible). Un exemple de signal X brut en fonction du temps en abscisse est illustré à la figure 5. Au cours d'une deuxième étape E2 postérieure à la première étape E1 , on calcule par un (ou plusieurs) calculateur CAL plusieurs spectres SP _t du signal temporel X sur respectivement plusieurs fenêtres temporelles successives (ou glissantes) W _t , qui sont définies par rapport à des instants t successifs, ainsi que représenté par exemple à la figure 2. Chaque fenêtre temporelle W _t peut être par exemple centrée sur l'instant t et a une largeur temporelle finie prédéterminée par rapport à l'instant t (par exemple les fenêtres temporelles successives W _t peuvent avoir des largeurs temporelles identiques). Le calculateur CAL calcule plusieurs signaux temporels partiels successifs Y (ou segments temporels Y), égaux respectivement au signal X multiplié par les fenêtres temporelles successives W _t pour les instants t successifs. Par exemple, les instants t successifs appartiennent à l'ensemble [1, ...,N _s ], où N _s est un nombre prédéfini, plus grand que 1. Les instants t peuvent être multiples d'un même pas d'échantillonnage à une constante prescrite près et être exprimés dans ce qui suit par un nombre entier de pas d'échantillonnage à une constante prescrite près. Le pas d'échantillonnage et/ou la fréquence d'échantillonnage peut être entré par l'utilisateur sur l'interface physique INT1 d'entrée de données ou être préenregistré dans le calculateur CAL. Le pas d'échantillonnage et/ou la fréquence d'échantillonnage définit le nombre d'échantillons acquis par seconde et les instants t.

Le (ou les) calculateur CAL peut être ou comprendre un ou plusieurs ordinateurs, un ou plusieurs serveurs, une ou plusieurs machines, un ou plusieurs processeurs, un ou plusieurs microprocesseurs. Le calculateur CAL peut comprendre une ou plusieurs interfaces physiques INT1 d'entrée de données, une ou plusieurs interfaces physiques INT2 de sortie de données. Cette ou ces interfaces physiques INT1 d'entrée de données peuvent être ou comprendre un ou plusieurs claviers d'ordinateur, un ou plusieurs ports physiques de communication de données, un ou plusieurs écrans tactiles, ou autres. Cette ou ces interfaces physiques INT2 de sortie de données peuvent être ou comprendre un ou plusieurs ports physiques de communication de données, un ou plusieurs écrans, ou autres. Le dispositif 1 d'estimation de la vitesse v de déplacement du véhicule ayant une ou plusieurs roues RO de roulement sur le sol comporte notamment le (ou les) capteur CAL, le (ou les) calculateur CAL et le (ou les) filtre FIL, ainsi qu'illustré à titre d'exemple à la figure 11. Le filtre FIL peut être compris dans le calculateur CAL ou être distinct du calculateur CAL. Le capteur CAL est configuré pour exécuter l'étape E1 mentionnée ci-dessus. Le filtre FIL est configuré pour exécuter l'étape E4 mentionnée ci-dessous. Le calculateur CAL est configuré pour exécuter les étapes E2, E3, E5, E6 (et éventuellement E4) et les sous-étapes E51 , E52, E53 mentionnées ci- dessous. Le calculateur CAL et/ou le filtre FIL et/ou le capteur CAP peut être programmé par un programme PR d'ordinateur comportant des instructions de code mettant en œuvre le procédé d'estimation de la vitesse v de déplacement du véhicule, à savoir pour effectuer l'étape E1 et/ou E2 et/ou E3 et/ou E4 et/ou E5 et/ou E6 et la sous-étape E51 et/ou E52 et/ou E53. Le calculateur CAL comporte une mémoire MEM non volatile (ou un média MEM d'enregistrement non volatile), dans laquelle le programme PR d'ordinateur est enregistré d'une manière permanente. La mémoire MEM non volatile ou le média MEM d'enregistrement non volatile est lisible par un ordinateur et comporte des instructions enregistrées dessus, qui lorsqu'elles sont exécutées par un processeur, amènent le processeur à exécuter l'étape E1 et/ou E2 et/ou E3 et/ouE4 et/ou E5 et/ou E6 et la sous-étape E51 et/ou E52 et/ou E53. Le signal X et les données mentionnées ci-dessus et ci-dessous, utilisées par le procédé, le dispositif 1 et le programme PR sont enregistrées dans la mémoire MEM.

Le dispositif 1 d'estimation de la vitesse v de déplacement du véhicule VE à roue RO peut être embarqué sur le véhicule VE. Le dispositif 1 d'estimation de la vitesse v de déplacement peut être un outil d'aide à la navigation inertielle et peut être intégré à un dispositif de navigation inertiel, qui est inclus dans une centrale inertielle présente sur le véhicule. La centrale inertielle peut comprendre 3 accéléromètres dans 3 directions non parallèles, 3 gyroscopes autour de 3 directions non parallèles, le dispositif de navigation inertiel et des fonctions de navigation.

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la deuxième étape E2 pour chaque signal temporel partiel Y le spectre SPt égal à une transformée de Fourier (pouvant être par exemple une transformée de Fourier rapide FFT), limitée à la fenêtre temporelle W _t respective, de ce signal temporel partiel Y (transformée de Fourier à courte durée, également appelée STFT en anglais). Le calculateur CAL calcule ainsi la transformée de Fourier sur les fenêtres glissantes W _t du signal temporel X, et donc la pluralité de spectres SPt pour respectivement la pluralité d'instants t successifs. Chaque spectre SPt donne une amplitude du spectre en fonction de la fréquence f.

Pour chaque capteur CAP fournissant un signal X de mesure, la série des spectres SPt ainsi obtenus forme pour les instants t successifs un spectrogramme SP. Le spectrogramme SP est calculé par le calculateur CAL. Chaque spectre SPt forme une colonne du spectrogramme SP, ainsi qu'illustré à titre d'exemple à la figure 3. Ainsi, si une multiplicité de capteurs CAP fournissent une multiplicité de signaux X, le calculateur CAL calcule de manière correspondante une multiplicité de spectrogrammes SP et une multiplicité de spectrogrammes S(f _t , t). Dans le cas de plusieurs spectrogrammes SP et S(f _t , t), ces spectrogrammes SP et S(f _t , t) peuvent subir des traitements différents pour le même capteur.

Une représentation en plan temps-fréquence, telle que par exemple celle des figures 3 et 6 à partir du signal de la figure 5, représente sous la forme d'une image : sur l'axe des temps en abscisses : les instants successifs t, sur l'axe des ordonnées : les fréquences f _t des spectres respectifs SPt pour chaque instant t, en niveaux de gris et en échelle logarithmique sur les axes verticaux des temps t: la valeur absolue de l'amplitude du spectre respectif SPt pour chaque instant t, la valeur absolue de cette amplitude étant rapportée à l'échelle ECH de niveaux de gris.

La figure 3 représente à titre d'exemple des spectres SPt dans le plan temps-fréquence d'un signal X inertiel obtenu par un capteur CAP inertiel formé par un accéléromètre latéral enregistré sur une voiture VE. Les exemples d'enregistrements des figures 5 à 10 ont été obtenus sur une voiture VE à partir d'un capteur CAP inertiel. La réponse du système mécanique dépend fortement du régime dans lequel il opère d'où l'intérêt de conditionner les résultats d'analyse du signal X à la vitesse v à chaque instant. Une représentation temps- fréquence permet de révéler les signatures liées aux sources mécaniques. Une telle représentation fournit une idée sur l'évolution du contenu fréquentiel au cours du temps. D'autres représentations temps-fréquence peuvent aussi être utilisées telles que la représentation de Wigner ou autres.

Dans les véhicules VE à roue RO, les roues RO peuvent émettre des vibrations qui se manifestent dans le plan temps-fréquence par une suite de pics (trajectoire ξ _t de fréquence). Les sources liées aux roues se manifestent par des trajectoires continues et lisses (trajectoire ξ _t de fréquence) sans sauts brusques de la vitesse v. La fréquence fondamentale (trajectoire ξ _t de fréquence) et ses harmoniques (autres composantes fréquentielles COMPF) sont les plus énergétiques dans le plan temps-fréquences. La fréquence de la vitesse v de la roue RO du véhicule VE se manifeste dans les spectres successifs par une trajectoire régulière variant selon l'état du système (accélération, décélération, arrêt etc.). Ainsi, une façon simple d'estimer la fréquence instantanée du signal X revient alors à chercher les positions des crêtes les plus énergétiques dans la bande de fréquence de la composante à suivre (trajectoire ξ _t de fréquence). La prise en compte des harmoniques multiples de la fondamentale des roues et l'ajout de contraintes de continuité et les informations sur la direction permettent d'atténuer les trajectoires continues liées à des composantes parasites telles que la combustion du moteur du véhicule.

Au cours d'une quatrième étape E4 postérieure à la deuxième étape E2, on filtre les spectres SPt par un (ou plusieurs) filtre FIL d'atténuation de bruit et/ou d'extraction de caractéristiques déterministes, pour obtenir un spectrogramme S(f _t , t) filtré du signal X sur les instants t successifs. Dans le spectrogramme S(f _t , t), la fréquence pour l'instant t est désignée par f _t . Ce prétraitement E4 permet de fournir une représentation temps-fréquence plus propre. Dans cette étape E4, on peut différencier les deux cas où la composante cible (vitesse v) est quasi-stationnaire ou fortement non-stationnaire. Dans le cas de composante cible (vitesse v) quasi-stationnaire, le filtre FIL peut comprendre un filtre de séparation (en anglais DRS : deterministic random separation) de caractéristiques déterministes (et donc prédictibles) par rapport à des caractéristiques aléatoires (et donc non prédictibles), y compris le bruit de fond et les composantes interférantes non-stationnaires. Dans le cas de composante cible (vitesse v) fortement non-stationnaire, le filtre FIL peut être configuré pour diminuer le bruit de fond et les fréquences invariantes en temps (par exemple bruit de ventilation, composantes électriques) et peut comprendre à cet effet un filtre médian glissant. Un exemple de spectrogramme S(f _t , t) filtré à partir du signal de la figure 5 est illustré à la figure 7. La variation maximum de fréquence (pour assurer la continuité du signal) du spectre filtré S(f _t , t) peut être déterminée par le calculateur CAL ou avoir été préenregistrée dans le calculateur CAL ou être prescrite par l'utilisateur sur l'interface physique INT1 d'entrée de données.

Au cours d'une cinquième étape E5 postérieure à la quatrième étape E4, on estime par le calculateur CAL une trajectoire ξ _t de fréquence, représentative de la vitesse v instantanée de la roue RO du véhicule VE dans le spectrogramme S(f _t , t) filtré sur les instants t successifs. Cette étape E5 comporte les sous-étapes E51, E52, E53 qui seront décrites ci- dessous. Dans ce qui suit, la trajectoire ξ _t de fréquence est également notée ξ.

En effet, la vitesse v de déplacement de la roue RO du véhicule VE se traduit par une certaine composante fréquentielle dans le plan temps-fréquence du spectrogramme S(f _t , t), cette composante fréquentielle pouvant avoir une certaine fréquence (égale dans ce qui suit à la trajectoire ξ _t de fréquence) variant dans le temps et une certaine amplitude variant dans le temps. Cette fréquence ξ _t variant au cours du temps t dans le spectrogramme S(f _t , t) forme donc une trajectoire ξ _t de fréquence de la vitesse v de la roue RO du véhicule VE en fonction du temps t dans le spectrogramme S(f _t , t) (et donc dans le plan temps-fréquence de la figure 3). Cette trajectoire ξ _t de fréquence est mélangée à des trajectoires d'autres fréquences d'autres composantes fréquentielles COMPF dans le spectrogramme S(f _t , t), ainsi qu'illustré à la figure 3. Suivant un mode de réalisation, pour la roue RO roulant sur le sol, la fréquence instantanée de rotation de la roue RO (nombre de tours par secondes) aux instants successifs t est la trajectoire ξ _t de fréquence en fonction du temps t. Suivant un autre mode de rélaisation, pour la roue RO roulant sur le sol, la vitesse tangentielle de la roue RO (nombre de tours par secondes multiplié par le rayon de la roue RO et multiplié par 2n) aux instants successifs t est la trajectoire ξ _t de fréquence en fonction du temps t. Soit N _s le nombre (supérieur ou égal à 2) de segments temporels ou d'instants successifs t dans le spectrogramme S(f _t , t) (axe des abscisses à la figure 2), c'est-à-dire que les instants t ∈ {1, ...,N _s }. Soit L _s le nombre (supérieur ou égal à 2) de bacs de fréquences dans le spectrogramme filtré S(f _t , t). On note f _t les fréquences de chaque spectrogramme S(f _t , t) à l'instant t. On note également par la vraie trajectoire ξ _t inconnue de fréquence à estimer. Le nombre N _s d'instants successifs t peut être déterminé par le calculateur CAL ou avoir été préenregistré dans le calculateur CAL ou être prescrit par l'utilisateur sur l'interface physique INT1 d'entrée de données.

Suivant un mode de réalisation de l'invention, on détermine au cours de la deuxième étape E2 par le calculateur CAL un support Q de valeurs {f ¹ , ...,f ^Ls } de fréquence, à savoir Ω = {f ¹ , ...,f ^Ls } Le support Q de valeurs {f ¹ , ...,f ^Ls } est le support des fréquences dus spectrogramme S(f _t , t) et peut comporter toutes les valeurs de fréquences comprises entre zéro et une fréquence maximum, par exemple des fréquences réparties régulièrement, pouvant être multiples d'un pas de fréquence prescrit à une constante prescrite près. Le nombre Ls de valeurs de fréquence du support Q est fini et peut être déterminé par le calculateur CAL ou avoir été préenregistré dans le calculateur CAL ou être prescrit par l'utilisateur sur l'interface physique INT1 d'entrée de données. Le pas de fréquence peut être calculé par le calculateur CAL en fonction de la fréquence d'échantillonnage et/ou en fonction de la durée de la fenêtre W _t .

On note également par la trajectoire ξ _t de fréquence à estimer. On a donc pour les instants successifs t allant de 1 à N _s , la trajectoire ξ _t de fréquence égale à un ensemble de fréquences ... ,ξ _Ns c'est-à-dire ξ _t = (ξ ₁ ,..,ξ _Ns ), qui est une variable aléatoire modélisant la trajectoire ξ _t à estimer. La trajectoire ξ _t de fréquence peut prendre les valeurs de fréquence ...,ξ _Ns dans le support Ω = {f ¹ , ...,f ^Ls } de valeurs de fréquence du spectrogramme S(f _t , t).

La chaîne des fréquences ...,ξ _Ns successives de la trajectoire ξ _t de fréquence au cours du temps t peut être décrite par une probabilité de transition pξ _t+1 |ξ _t et une probabilité initiale donnée pξ ₁ .

Suivant un mode de réalisation, au cours de la troisième étape E3 postérieure à la deuxième étape E2 et antérieure à la quatrième étape E4, on prescrit par le calculateur CAL une probabilité de transition pξ de la trajectoire ξ _t de fréquence selon une chaîne de Markov liant la probabilité de transition de la fréquence ξ _t de cette trajectoire à l'instant présent t à la fréquence ξ _{t +1} de cette trajectoire à l'instant suivant t+1, en partant de la probabilité initiale donnée pξ ₁ . Cette probabilité de transition pξ suit une loi a priori. Dans ce qui suit, la probabilité de transition pξ de la trajectoire ξ _t de fréquence est notée pξ (f ₁ , ...,f _Ns )- Un exemple de probabilité de transition p _ξ calculée à partir du signal de la figure 5 est illustré à la figure 9. On a spécifié à l'avance la probabilité de transition décrivant les évolutions possibles de la vitesse au cours du temps. Concrètement, une grande variation de fréquence a une probabilité faible et une petite variation a une probabilité haute. Cette probabilité de transition définit ce qu'on appelle une chaîne de Markov.

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de l'étape E3 la probabilité p _ξ (f ₁ ,..,f _Ns ) de transition de la trajectoire ξ _t de fréquence pour l'instant N _s en partant de la probabilité initiale donnée pξ ₁ de l'instant 1 , pour t ∈ {1, ..., N _s } et pour chaque fréquence f _t ∈ Ω (c'est-à-dire chaque f ₁ ,..,f _Ns ∈ Ω) selon l'équation suivante : pξ ₁ (f ₁ ) est la probabilité initiale donnée de l'instant t=1.

Dans cette équation, Pξ _t+1 |ξ _t =f _t (f _t+1 ) est la probabilité de transition de la fréquence ξ _t de la trajectoire de l'instant présent t à la fréquence ξ _{t +1} de la trajectoire de l'instant suivant t+1.

On a donc la loi jointe de la trajectoire ξ de fréquence, qui est donnée par :

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de l'étape E3 la probabilité de transition de la fréquence ξ _t de la trajectoire de l'instant présent t à la fréquence ξ _{t +1} de la trajectoire de l'instant suivant t+1 selon une loi conditionnelle gaussienne N centrée sur la fréquence f _t du spectrogramme filtré (S(f _t , t)) à l'instant t, c'est-à-dire est la variance de la loi conditionnelle gaussienne N et la fréquence est la moyenne de la loi conditionnelle gaussienne N. Cela permet de maintenir une certaine régularité de l'harmonique à suivre et d'ajouter des contraintes de continuités entre les instants consécutifs t et t+1.

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de l'étape E3 l'écart-type σ _t de la loi conditionnelle gaussienne N comme étant σ _t = Δt.Yt, où Δt est une résolution temporelle calculée du spectrogramme filtré (S(f _t , t)), Yt est une vitesse moyenne calculée avec laquelle la trajectoire ξ _t de fréquence évolue entre l'instant t et l'instant suivant t+1. Ainsi, en diminuant la variance , l'effet de lissage sera plus fort. Le calcul de la loi de transition Pξ _t+1 |ξ _t =f _t (f _t+1 ) par le calculateur CAL peut aussi tenir compte de la direction d'évolution (par exemple montée, descente, stabilité) de la trajectoire ξ _t de fréquence lorsque cette direction est fournie à chaque instant t. Par exemple, le calculateur CAL peut utiliser une loi de transition qui est tronquée sur ξ _{t +1} > ξ _t dans le cas d'une montée de la trajectoire ξ _t de fréquence et qui est tronquée sur ξ _{t +1} < ξ _t dans le cas d'une descente de la trajectoire ξ _t de fréquence. Dans le cas où la trajectoire ξ _t de fréquence devient constante (arrêt de descente et de montée), par exemple à une valeur nulle (arrêt du véhicule VE correspondant à une vitesse v nulle), il faut mettre y _t à une valeur assez faible.

A la première sous-étape E51 de la cinquième étape E5, on estime par le calculateur CAL une probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE à partir d'une amplitude calculée du spectrogramme S(f _t , t) filtré pour chaque instant t. Dans ce qui suit, la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t ) de fréquence de vitesse v de déplacement du véhicule VE est notée Un exemple de probabilité p _Zt ǀ _ξt d'observation calculée à partir du signal de la figure 5 est illustré à la figure 8.

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la sous-étape E51 la probabilité p _Zt ǀ _ξt d'observation sachant une trajectoire comme dépendant de l'amplitude calculée du spectrogramme S(f _t , t) filtré et normalisé pour chaque instant t.

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la sous-étape E51 la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence pour t ∈ {1, , N _s } et pour chaque fréquence f _t ∈ Ω (c'est-à-dire chaquef ₁ ,..,f _Ns ∈ Ω) selon l'équation suivante :

25

Le calculateur CAL calcule au cours de la sous-étape E51 le vecteur binaire Z ou Z _t d'observation ayant une taille N _s prescrite pour l'instant t ∈ {1, , N _s } (et pour chaque fréquence f _t ∈ Ω, c'est-à-dire chaquef ₁ ,..,f _Ns E n) et étant tel que chacune des coordonnées Z _t [f] du vecteur binaire Z _t d'observation (les coordonnées Z _t [f] du vecteur binaire Z _t d'observation étant donc (z _1, ...,z _Ns )) est égale à :

Dans ce qui précède, Z = (z _1, ...,z _Ns ) est l'observation.

La densité ou probabilitép _Zt ǀ _ξt d'observations s'écrit alors comme suit : oi

5 d'observation (1))

La probabilité d'observation p _Zt ǀ _ξt est définie à partir du l'amplitude |S(f _t ,t)| du spectrogramme pour f _t ∈ Ω. En effet, on peut avoir une idée sur la trajectoire able ft et alors confirmer si ξ _t = f _t est très probable d'être observé ou non selon le niveau de l'amplitude du spectrogramme au point par rapport aux autres points (f, t), f ∈ Ω.

L'idée principale de la méthode utilisée est de calculer une densité de probabilité à partir du spectrogramme décrivant la répartition des observations sachant la fréquence cible f _t à chaque instant t. En effet, pour ξ _t = f _t ∈ Ω donné, le coefficient de Z _t pour la fréquence f _t a une forte probabilité d'être égale à 1 à l'instant t, c'est-à-dire il est très probable que la vraie fréquence est égale à f _t .

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la sous-étape E51 la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence selon les équations suivantes : pour avec = 1) dépendant de l'amplitude calculée du spectrogramme (S(f _t , t)) filtré pour chaque instant (t),

Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la sous-étape E51 la probabilité p _Zt ǀ _ξt d'observation selon l'équation suivante : p _Zt ǀ _ξt = ft Z _t où iSCft.Ol est l'amplitude calculée du spectrogramme filtré pour chaque instant t, f _min (t) est une trajectoire i minimale de fréquence ayant été mesurée, f _max (t) est une trajectoire i maximale de fréquence ayant été mesurée, f _min (t) et f _max (t) sont telles que est une fonction prenant la valeur 1 pour f _t entre f prenant la valeur 0 pour f _t en-dessous de ou au-dessus de f _max (t). Suivant un mode de réalisation de l'invention, le calculateur CAL calcule au cours de la sous-étape E51 où |S(f _t ,t)| est l'amplitude calculée du spectrogramme filtré pour chaque instant t. Le calculateur CAL calcule au cours de la sous-étape E51 le coefficient diviseur c _t comme étant le maximum de z _t [f] pour la fréquence f située entre la trajectoire minimale de fréquence et la trajectoire maximale f _max (t) de fréquence, c'est-à-dire :

Le calculateur CAL mesure ou détermine au cours de la sous-étape E51 la trajectoire minimale de fréquence. Le calculateur CAL mesure ou détermine au cours de la sous- étape E51 la trajectoire maximale f _max (t) de fréquence. La trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE est comprise entre la trajectoire minimale de fréquence et la trajectoire maximale f _max (t) de fréquence. On a donc f _max (t)*

Le calculateur CAL calcule au cours de la sous-étape E51 la fonction comme étant une fonction à créneau unique entre f _min (t) c'est-à-dire la fonction prenant la valeur 1 pour f _t entre et f _max (t) et prenant la valeur 0 pour f _t en-dessous de ou au-dessus de

La trajectoire minimale de fréquence et la trajectoire maximale f _max (t) de fréquence permettent de ne tenir compte que des pics des composantes au voisinage de la trajectoire ξ _t de fréquence à suivre et définissent une bande fréquence dans laquelle se trouve la trajectoire ξ _t de fréquence.

Suivant un mode de réalisation de l'invention, la trajectoire minimale de fréquence est fixe dans le temps t, c'est-à-dire = f _min .

Suivant un mode de réalisation de l'invention, la trajectoire maximale f _max (t) fréquence est fixe dans le temps t, c'est-à-dire f _max (t) = f _max .

Suivant un autre mode de réalisation de l'invention, la trajectoire minimale de fréquence est variable dans le temps t.

Suivant un autre mode de réalisation de l'invention, la trajectoire maximale f _max (t) de fréquence est variable dans le temps t.

Dans le cas où la trajectoire ξ _t de fréquence cible s'étale sur une large bande de fréquence (par exemple rampe) où d'autres harmoniques parasites peuvent coexister, il est préférable de choisir une bande de fréquence variable dans le temps t. Cela revient à définir deux trajectoires f _min (t) et f _max (t) enveloppant la trajectoire ξ _t de fréquence. Par exemple, la bande de fréquence choisie ne contient qu'une seule trajectoire continue complète ξ _t sur toute la période du temps t considérée. Les contraintes de continuité qui sont introduites par la loi de transition pξ ont comme effet de donner plus de poids aux pics liés à une courbe régulière de la trajectoire ξ _t de fréquence. Suivant un mode de réalisation de l'invention, les trajectoires f _min (t) et f _max (t) peuvent être calculées automatiquement par le calculateur CAL. Suivant un mode de réalisation de l'invention, les trajectoires f _min (t) et f _max (t) peuvent être prescrites par un utilisateur au calculateur CAL par sélection aléatoire des points du spectrogramme S(f _t , t) enveloppant la suite des crêtes à suivre (maximums d'amplitude) du spectrogramme S(f _t , t), puis par interpolation de ces points. La définition de la bande de fréquence est optionnelle mais permet d'améliorer considérablement les performances de la méthode et de minimiser le coût de calcul surtout dans le cas où le signal observé contient un grand nombre d'harmoniques. On peut également identifier les intervalles d'arrêt et inclure cette information dans la définition de f _min et f _max et aussi dans la loi de transition par la suite. Si cette information n'est pas disponible, une façon simple d'estimer les instants d'arrêt à partir du spectrogramme S(f _t , t) est de comparer la somme des amplitudes du spectrogramme S(f _t , t) à chaque instant t ₀ sur toutes les fréquences f à un seuil donné ( Si , alors arrêt à l'instant t ₀ , où E > 0 est un seuil à fixer (par exemple la médiane de toutes les valeurs sur tout l'axe de temps)). La trajectoire minimale f _min (t) de fréquence et/ou la trajectoire maximale f _max (t) de fréquence peut être entrée par l'utilisateur sur l'interface physique INT1 d'entrée de données ou avoir été préenregistrée dans le calculateur CAL. Dans le cas où le dispositif 1 est intégré à un dispositif de navigation du véhicule VE, des données a priori possibles tels que par exemple f _min = 0 et f _max peuvent être entrées par l'utilisateur sur l'interface physique INT1 d'entrée de données ou avoir été préenregistrée dans le calculateur CAL, pour une navigation en temps réel.

La probabilité de l'observation Z _t à chaque instant t est calculée par le calculateur CAL sachant la fréquence cible ξ _t = f _t ∈ Ω.

Suivant un mode de réalisation de l'invention, dans le cas d'un signal monocomposant avec une Fl , si on dispose du spectrogramme S(f _t , t) du signal sans bruit alors c'est-à-dire le spectre Z _t ne peut être observé que pour la fréquence f. Dans le cas contraire, Pξ _t |ξ _t =ft est importante pour toutes les fréquences f _t où l'amplitude du spectrogramme S(f _t , t) est maximale à l'instant t. Ainsi, la fonction f _t → Pξ _t |ξ _t =ft(Zt) Peut être multimodale dans le cas où plusieurs pics (le plus souvent liés au bruit) sont présents à l'instant t dans f _min (t) f _max (t) - C'est-à-dire ce même spectre observé peut correspondre à plusieurs solutions de fréquences plausibles f _t . En considérant l'amplitude du spectrogramme S(f _t , t) pour plusieurs multiples de l'harmonique cible, nous pourrons favoriser un mode par rapport à un autre. Suivant un mode de réalisation de l'invention, Le calculateur CAL calcule au cours de la sous-étape E51 la probabilité p _Zt ǀ _ξt d'observation selon l'équation suivante : (la probabilité s'écrit dans ce cas comme le produit de plusieurs probabilités) : où J peut être entier ou réel positif quelconque, vérifiant et où F _s est une fréquence de seuil, ayant été prédéterminée en fonction du spectrogramme S(ft, t).

A la deuxième sous-étape E52 de la cinquième étape E5, on estime par le calculateur CAL une loi a posteriori p _ξ , _z d' observation de la trajectoire ξ _t de fréquence ou de la fréquence ξ _t . Cette loi a posteriori p _ξ , _z d'observation est proportionnelle au produit de la probabilité de transition pξ de la trajectoire ξ _t de fréquence et de la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence, c'est-à-dire Cette loi a posteriori p _ξ , _z d'observation de la trajectoire ξ _t de fréquence est également notée

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation égale au produit de la probabilité de transition pξ de la trajectoire ξ _t de fréquence, de la probabilité p _Zt ǀ _ξt d'observation de la trajectoire ξ _t de fréquence et d'un facteur multiplicatif C calculé par le calculateur CAL, c'est-à-dire Suivant un mode de réalisation de l'invention, le facteur multiplicatif C peut être pris égal à 1.

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule le facteur multiplicatif C égal à

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation égale à

Suivant un mode de réalisation de l'invention, à la deuxième sous-étape E52 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation

A la troisième sous-étape E53 de la cinquième étape E5, le calculateur CAL estime la trajectoire ξ _t de fréquence pour les instants t successifs à partir de la loi a posteriori p _ξ , _z d'observation de la trajectoire ξ _t de fréquence. Le calculateur CAL calcule comme trajectoire ξ _t de fréquence la trajectoire optimale parmi toutes les trajectoires possibles dans le ou les spectrogrammes S(f _t , t) à partir de la loi a posteriori p _ξ , _z d'observation de la trajectoire ξ _t de fréquence, sachant les observations Z données par la probabilité p _Zt ǀ _ξt d'observation. Un exemple de trajectoire ξ _t de fréquence calculée à partir du signal de la figure 5 est illustré à la figure 10. La trajectoire de fréquence ξ _t ou / et/ou une ou des valeurs de fréquences de la trajectoire de fréquence ξ _t ou / (pouvant être égale à la fréquence instantanée de la roue RO) peut être fournie par le calculateur CAL sur l'interface physique INT2 de sortie de données, par exemple en étant affichée sur un écran. Un intervalle de confiance du calcul, ayant été calculé par le calculateur, peut être fourni par le calculateur CAL sur l'interface physique INT2 de sortie de données, par exemple en étant affiché sur un écran. Dans le cas où le dispositif 1 est intégré à un dispositif de navigation du véhicule VE, la trajectoire de fréquence ξ _t ou / et/ou une ou des valeurs de fréquences de la trajectoire de fréquence ξ _t ou / peut être stockée dans la mémoire M ou dans une autre mémoire et être évaluée en temps réel dans un autre filtre pour une navigation en temps réel. Suivant un mode de réalisation de l'invention, à la troisième sous-étape E53 de la cinquième étape E5, le calculateur CAL estime pour chaque instant t la trajectoire ξ _t de fréquence comme étant dans le support Q celle des valeurs f ¹ , ...,f ^Ls de fréquence, qui maximise la loi a posteriori p _ξ , _z d'observation de la trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE à cet instant t. Dans ce cas, le calculateur CAL sélectionne pour chaque instant t la trajectoire ξ _t de fréquence dans le support Q des valeurs de fréquence. On a donc dans ce cas la trajectoire ξ _t de fréquence, qui est égale à Le calculateur comporte un estimateur du maximum a posteriori maximisant La résolution de ce problème d'optimisation peut être facilement effectuée par exemple en utilisant la programmation dynamique.

Suivant un mode de réalisation de l'invention, à la troisième sous-étape E53 de la cinquième étape E5, le calculateur CAL calcule la loi a posteriori p _ξ , _z d'observation par convolution (par exemple bidimensionnelle) de la probabilité p _Zt ǀ _ξt d'observation par des noyaux Gaussiens de la probabilité de transition pξ.

Suivant un mode de réalisation de l'invention, à la troisième sous-étape E53 de la cinquième étape E5, le calculateur CAL calcule récursivement la loi a posteriori selon l'équation suivante :

Où * est l'opérateur de convolution et N est un noyau gaussien représentant la densité Ces opérations ne sont pas coûteuses.

Il estime ensuite pour chaque instant t la trajectoire ξ _t de fréquence ou la fréquence ξ _t comme étant un barycentre de la loi a posteriori qui est prise en respectivement les valeurs f ¹ , ...,f ^Ls de fréquence du support Q et qui est affectée respectivement de ces valeurs f ¹ , ...,f ^Ls de fréquence. La fréquence ξ _t , ainsi obtenue, est dans ce cas une combinaison linéaire des valeurs f ¹ , ...,f ^Ls de fréquence du support Q, dont les poids sont formés par la loi a posteriori p _ξ , _z , prise en respectivement les valeurs f ¹ , ...,f ^Ls de fréquence du support Q. On a donc dans ce cas la trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE, qui est égale à

La trajectoire ξ _t de fréquence est ainsi une sorte de moyenne a posteriori, qui génère des résultats plus lisses. Au cours d'une sixième étape E6 postérieure à la cinquième étape E5, le calculateur CAL estime la vitesse v de déplacement de la roue RO du véhicule VE pour les instants t successifs à partir de la trajectoire ξ _t de fréquence. La vitesse v de déplacement de la roue RO du véhicule VE peut être fournie par le calculateur CAL sur l'interface physique INT2 de sortie de données, par exemple en étant affichée sur un écran. Dans le cas d'un dispositif de navigation, la vitesse v de déplacement de la roue RO du véhicule VE peut être stockée dans la mémoire MEM ou dans une autre mémoire et être évaluée en temps réel dans un autre filtre pour une navigation en temps réel.

Suivant un mode de réalisation de l'invention, à la sixième étape E6, le calculateur CAL estime la vitesse v de déplacement de la roue RO du véhicule VE pour l'instant t en multipliant la trajectoire ξ _t de fréquence de vitesse v de déplacement du véhicule VE par le périmètre 2nR de la roue RO, où R est le rayon de la roue RO et est prédéterminé ou estimé, par exemple au cours de la navigation par un procédé quelconque. La vitesse v est également notée et est dans ce cas une vitesse tangentielle. On a donc Le rayon R et/ou le diamètre de la roue RO et/ou le périmètre 2nR de la roue RO peut avoir été entré par l'utilisateur sur l'interface physique INT1 d'entrée de données ou avoir été préenregistré dans le calculateur CAL.

Suivant un autre mode de réalisation de l'invention, à la sixième étape E6, le calculateur CAL estime la vitesse v de déplacement de la roue RO du véhicule VE pour l'instant t comme étant la fréquence de rotation de la roue RO du véhicule VE (nombre de tours par secondes), égale à la trajectoire ξ _t de fréquence à l'instant t.

L'invention présente notamment les avantages suivants :

- Suppression d'un capteur, plus besoin d'un odomètre. Cela est très intéressant, surtout en aéronautique où les capteurs sont très chers non seulement pour leur résistance à un environnement extrême (forte pression et température) mais aussi pour une raison logistique où le moindre capteur change les propriétés modales de la structure.

- Une estimation précise de la vitesse instantanée v en utilisant l'information présente dans une multitude d'harmoniques linéairement liés, permettant d'estimer les variations fines de cette vitesse au lieu d'une estimation grossière.

- Le problème d'interception de fréquences ou de régimes (communément présente dans les moteurs d 'aéronef) est également considéré, ainsi que la continuité de cette variation.

- La méthode offre la possibilité d'améliorer de plus en plus la précision de l'estimation en permettant à l'utilisateur de donner des a priori sur la vitesse (par exemple si on est en montée de régime ou en régime quasi-stationnaire, selon la figure 4). Bien que ces informations ne soient pas obligatoires, elles peuvent améliorer significativement la qualité du signal obtenu.

Bien entendu, les modes de réalisation, caractéristiques, possibilités et exemples décrits ci-dessus peuvent être combinés l'un avec l'autre ou être sélectionnés indépendamment l'un de l'autre.