FR 2937423 A1 20100423 FR 0857181 A 20081022 Dispositif de détermination d'une trajectoire formée de positions successives sensiblement coplanaires d'un accéléromètre triaxial lié de manière solidaire à un élément mobile L'invention porte sur un dispositif de détermination d'une trajectoire formée de positions successives sensiblement coplanaires d'un accéléromètre triaxial lié à un élément mobile. De manière générale, l'invention est liée aux domaines de la capture de mouvement, et de la géolocalisation d'individus en déplacement. La capture de mouvement concerne des applications grand public telles des applications de loisir (consoles de jeu interactives, suivi du geste sportif, réalité virtuelle ou réalité augmentée), des applications d'aide à la navigation de piétons (les plus utilisés à l'heure actuelle étant les systèmes de navigation par satellites tel le GPS), des applications 15 d'aide à la mobilité de personnes vulnérables ou temporairement fragilisées par leur environnement (personnes handicapées ou plongées dans l'obscurité), et des applications de fitness (podomètre, calcul de la dépense énergétique ou de la distance parcourue). La capture de mouvement concerne également des applications médicales (suivi des 20 personnes âgées et/ou dépendantes, analyse de la marche pour la rééducation posturale ou aide au diagnostic), des applications de sécurité ou de secours (localisation de pompiers à l'intérieur d'un bâtiment en feu, suivi opérationnel de militaires, ou surveillance de prisonniers), mais aussi des applications à visée commerciale (statistiques sur les trajectoires 25 suivies par les consommateurs dans les centres commerciaux ou les supermarchés, définition d'archétypes d'usage, ou proposition de services commerciaux topo-dépendants). Il est notamment connu de reconstituer le déplacement d'un objet ou d'une personne munie d'un objet comprenant un émetteur de 30 signaux reconnus par un système d'aide à la navigation par satellites, tel le système GPS, en extérieur, et par un système de radiolocalisation à base de transmissions Ultra Large Bande (ULB) ou WiFi en intérieur. Dans de nombreuses zones géographiques, externes ou internes, la navigation à base de système d'aide à la navigation par 35 satellites ou de radiolocalisation s'avère néanmoins très délicate, en raison de l'obstruction des signaux radio nécessaires à la mesure de métriques topo-dépendantes, par exemple en cas d'indisponibilité des signaux émis par un ou plusieurs satellites du système d'aide à la navigation, nécessaires à la mesure de pseudo-distances. En pratique, ces situations de blocage ou d'obstruction dégradent fortement la précision de l'information de localisation, et rendent même parfois le service d'aide à la navigation indisponible, comme illustré dans le document "An evaluation of UWB localization under non line-of-sight (NLOS) propagation" (A. Maali, A. Ouldali, H. Mimoun, et G. Baudoin), ~o publié dans Wireless Pervasive Computing, ISWPC, 3rd International Symposium on, pages 379û382, en mai 2008. Généralement, les réalisations citées précédemment présentent un coût, une consommation, un encombrement, et/ou une infrastructure inadaptés, voire rédhibitoires pour des applications grand public. 15 D'autres systèmes, ont pour but d'élaborer un trajet à partir d'un point de départ donné, en utilisant des capteurs d'attitude (accéléromètre(s), magnétomètre(s), gyromètre) délivrant le plus souvent des mesures inertielles. Ces mesures permettent d'assurer la navigation par des moyens plus ou moins précis et compliqués, tels que décrits, par 20 exemple dans les documents "Inertial head-tracker sensor fusion by a complementary separate-bias kalman filter" (E. Foxlin) publié en mars 1996 dans Virtual Reality Annual International Symposium, Proceedings of the IEEE 1996, pages 185û194; "Detection of spatio-temporal gait parameters by using wearable motion sensors" (K. Kogure L. Seon-Woo, 25 et K. Mase), publié en 2005 dans Engineering in Medicine and Biology Society. IEEE-EMBS 2005. 27th Annual International Conference of the, pages 6836û6839, 2005; "Pedestrian tracking with shoe-mounted inertial sensors" (E. Foxlin), publié en novembre 2005 dans Computer Graphics and Applications, IEEE, 25 :38û46, Nov-Dec 2005; ou "Integration of foot- 30 mounted inertial sensors into a Bayesian location estimation framework" (P. Robertson B. Krach), publié en mars 2008 dans Positioning, Navigation and Communication, 2008. WPNC 2008. 5th Workshop on, pages 55û61. De tels systèmes sont de coût et de complexité élevés, et manquent souvent de précision. Des consoles de jeu vidéo, telle la Wii, utilisent des capteurs optiques et/ou à ultrasons pour déterminer la trajectoire d'un élément de commande de jeu. Ces systèmes sont couteux et limités. Il existe également, comme présentés dans les documents "Assessment of walking features from foot inertial sensing" (S. Scapellato F. Cavallo A.M. Sabatini, et C. Martelloni), publié en mars 2005 dans Biomedical Engineering, IEEE Transactions on, pages 486û 494, ou "Multisensor approach to walking distance estimation with foot inertial sensing" (D. Alvarez A.M Lopez J. Rodriguez-Uria J.C. Alvarez, et R.C. Gonzalez), publié en août 2007 dans Engineering in Medicine and Biology Society, 2007, EMBS 2007. 29th Annual International Conference of the IEEE, pages 5719û5722, des systèmes dans lesquels on suppose deux accéléromètres et un gyromètre disposés dans le plan sagittal d'un utilisateur, qui permettent de déterminer l'attitude de ce capteur dans le plan. L'accélération est intégrée dans l'axe du pied. Ces systèmes ont un coût réduit, mais présentent d'autres inconvénients. D'une part, ils supposent que les capteurs sont parfaitement placés dans le plan sagittal, ce qui est quasiment impossible à réaliser concrètement par l'utilisateur, et provoque une erreur d'estimation liée au mauvais positionnement du ou des capteurs. D'autre part, ils supposent que la marche a lieu dans le plan sagittal, ce qui peut ne pas être le cas, lors d'une marche sur le côté, par exemple. Le document, "An innovative shoe-mounted pedestrian navigation system" (K. Fyfe, Gérard Lachapelle, R. Stirling, et J. Collin), publié en avril 2003 dans Proc. European Navigation Conf. (GNSS), CD-ROM, Austrian Inst. of Navigation, divulgue un système muni de trois accéléromètres et d'un quatrième accéléromètre afin de calculer l'angle subit par le capteur. Cet accéléromètre supplémentaire est placé sur le même élément mobile, mais à distance de ce dernier. Les deux capteurs voient donc la même rotation et le même déplacement. Une telle structure permet d'estimer la vitesse de rotation de l'élément mobile, en rotation autour d'une direction sensiblement constante, en se passant d'un gyromètre. Ce système a un coût réduit du au remplacement d'un gyromètre mono-axe par un couple d'accéléromètres par rapport aux systèmes avec des gyromètres proposés par Alvarez et al et Scapellato et al. Cependant, un tel système présente les mêmes inconvénients que l'exemple cité ci-dessus, car seule la rotation autour de l'axe orthogonal au plan sagittal est prise en compte. En outre, l'estimation de la rotation par un couple d'accéléromètres, tout comme celle réalisée à l'aide d'un gyromètre, présente des dérives importantes. Un cap est également déterminé à l'aide d'un magnétomètre, par exemple lorsque le pied est posé, ce qui donne la direction du pied et non pas la direction du déplacement recherchée. L'invention a pour but de déterminer une trajectoire à coût réduit, et de résoudre les problèmes précédemment cités. Selon un aspect de l'invention, il est proposé un dispositif de détermination d'une trajectoire formée de positions successives sensiblement coplanaires d'un accéléromètre triaxial lié de manière solidaire à un élément mobile d'un corps humain, animal ou artificiel, entre un premier instant d'immobilité et un deuxième instant d'immobilité de l'accéléromètre triaxial, consécutif audit premier instant d'immobilité, l'accéléromètre triaxial étant en rotation par rapport à une direction sensiblement constante entre lesdits premier et deuxième instants d'immobilité. Ledit dispositif comprend, en outre, un capteur additionnel de mesure d'un vecteur de champ vectoriel sensiblement constant entre lesdits premier et deuxième instants d'immobilité, dans un repère global fixe lié au repère terrestre, ledit capteur additionnel étant lié de manière solidaire audit élément mobile et d'axes de mesures sensiblement parallèles aux axes de mesures de l'accéléromètre, et des moyens de commande. Lesdits moyens de commande comprennent : des premiers moyens de détermination desdits premier et deuxième instants d'immobilité de l'accéléromètre triaxial ; des deuxièmes moyens de détermination du plan de ladite trajectoire de l'accéléromètre triaxial, à partir des vecteurs délivrés par l'accéléromètre triaxial, ou à partir des vecteurs délivrés par le capteur additionnel, à différents instants successifs entre le premier et le deuxième instants d'immobilité, dans un repère mobile lié à l'accéléromètre triaxial et au capteur additionnel ; des premiers moyens de calcul, auxdits instants successifs, des projections orthogonales du vecteur délivré par l'accéléromètre triaxial et du vecteur délivré par le capteur additionnel dans ledit repère mobile, sur ledit plan de ladite trajectoire de l'accéléromètre triaxial ; des troisièmes moyens de détermination, auxdits instants successifs, de la rotation de passage dudit repère mobile audit repère global fixe, à partir desdites projections orthogonales des vecteurs délivrés par le capteur additionnel ; des deuxièmes moyens de calcul, auxdits instants successifs, de la projection orthogonale dans ledit repère global fixe, du vecteur délivré par l'accéléromètre triaxial fournie par lesdits premiers moyens de calcul ; des troisièmes moyens de calcul pour soustraire au vecteur délivré par lesdits deuxièmes moyens de calcul auxdits instants successifs, le vecteur moyen sur lesdits instants successifs, de manière à obtenir les accélérations centrées dans ledit plan de ladite trajectoire de l'accéléromètre triaxial, dépourvues de l'influence de la gravité terrestre et de dérives dudit dispositif, dans ledit repère global fixe ; et des quatrièmes moyens de calcul de la double intégration temporelle desdites accélérations centrées délivrées par lesdits troisièmes moyens de calcul, de manière à obtenir la trajectoire de l'accéléromètre triaxial dans ledit plan, exprimée dans ledit repère global fixe, entre lesdits premier et deuxième instants d'immobilité. Un tel dispositif permet, à coût réduit de déterminer une trajectoire planaire d'un accéléromètre triaxial lié à un élément mobile, de manière précise et à coût réduit. L'élément mobile peut être une partie du corps humain comme un pied ou une main, ou une partie d'un corps animal, ainsi qu'une partie d'un corps artificiel tel un robot ou une souris informatique. Dans la suite de la description, l'invention sera décrite pour un pied d'un corps humain, mais peut s'appliquer par analogie à tout élément mobile d'un corps humain, animal ou artificiel. Aussi, dans un mode de réalisation, ledit élément mobile est un pied, l`intervalle de temps séparant lesdits premier et deuxième instants d'immobilité de l'accéléromètre triaxial correspondant à un pas dudit pied. Selon un mode de réalisation, lesdits moyens de commande comprennent des premiers moyens d'élaboration d'un vecteur de direction, adapté pour relier les positions respectives de l'accéléromètre triaxial audit premier instant d'immobilité et audit deuxième instant d'immobilité, exprimées dans ledit repère global fixe. La direction de ce vecteur dans le repère global fixe donne la direction du déplacement, et sa longueur la distance parcourue. Dans un mode de réalisation, ledit plan de la trajectoire étant vertical, lesdits moyens de commande comprennent, en outre, des quatrièmes moyens de calcul d'une projection orthogonale, dudit vecteur de champ vectoriel mesuré par le capteur additionnel, sur un plan horizontal orthogonal audit plan de la trajectoire, et des deuxièmes moyens d'élaboration d'un cap de déplacement par rapport audit champ vectoriel mesuré par le capteur additionnel. Si le capteur additionnel est un magnétomètre, on obtient la direction du déplacement par rapport au Nord. On peut, par exemple, tracer la direction du déplacement sur une carte ou donner les coordonnées de ce vecteur, indiquant ainsi dans quelle direction l'utilisateur s'est déplacé. Dans un mode de réalisation, ledit plan de la trajectoire étant vertical, lesdits deuxièmes moyens de détermination du plan de ladite trajectoire de l'accéléromètre triaxial sont adaptés pour effectuer une décomposition en valeurs singulières de la matrice à trois colonnes formée par les vecteurs délivrés par l'accéléromètre triaxial auxdits instants successif, dans ledit repère mobile, et extraire les deux vecteurs correspondant aux deux plus grandes valeurs singulières. Le dispositif est en translation dans un plan et effectue une rotation autour d'une direction sensiblement constante. C'est un moyen de calcul peu coûteux en nombre d'opérations et robuste. D'autre part on s'affranchit de l'orientation du capteur sur l'élément mobile. On ne fait pas d'hypothèse sur la direction de la trajectoire qui peut être quelconque. Dans un mode de réalisation, ledit plan de la trajectoire étant vertical, lesdits deuxièmes moyens de détermination du plan de ladite trajectoire de l'accéléromètre triaxial sont adaptés pour effectuer une décomposition en valeurs propres et vecteurs propres, dans ledit repère mobile, la matrice d'autocorrélation de la matrice à trois colonnes formée par les vecteurs délivrés par l'accéléromètre triaxial auxdits instants successifs, et extraire les deux vecteurs propres correspondant aux deux plus grandes valeurs propres. De même, le dispositif est en translation dans un plan et effectue une rotation autour d'une direction sensiblement constante. C'est un moyen de calcul peu coûteux en nombre d'opérations et robuste. D'autre part on s'affranchit de l'orientation du capteur sur l'élément mobile. On ne fait pas d'hypothèse sur la direction de la trajectoire qui peut être quelconque. Selon un mode de réalisation, lesdits deuxièmes moyens de calcul comprennent des cinquièmes moyens de calcul de passage dans un repère intermédiaire, lié au repère terrestre, défini par un premier vecteur délivré par lesdits premiers moyens de calcul à partir d'un vecteur 15 fourni par le capteur additionnel, et par un deuxième vecteur, orthogonal audit premier vecteur et appartenant audit plan de la trajectoire. On minimise ainsi le nombre de calculs à effectuer puisque le premier vecteur est disponible et on utilise seulement quelques multiplications et additions pour passer au repère global. 20 Selon un mode de réalisation, dans lequel lesdits troisièmes moyens de détermination de la rotation de passage dudit repère mobile audit repère global fixe comprennent des moyens d'estimation d'un angle, ladite rotation de passage étant définie par une projection orthogonale sur ledit champ vectoriel sensiblement constant, suivie d'une rotation dudit 25 angle estimé dans ledit plan de trajectoire de l'accéléromètre triaxial. Cette estimation permet d'orienter la trajectoire dans un repère global connu, tel le repère NED ou "North-East-Down" en langue anglo-saxonne. Par exemple lesdits moyens d'estimation sont adaptés pour 30 estimer ledit angle à partir de la pente de la surface de contact de l'élément mobile, et du vecteur transmis par lesdits cinquièmes moyens de calcul à partir du vecteur transmis par l'accéléromètre triaxial. La pente peut être connue par une table de correspondance prédéterminée fournissant la pente de la surface de contact de l'élément 35 mobile à partir de la position de l'accéléromètre. Cette méthode est robuste dès lors que la pente est précisément connue (dans les bâtiments par exemple la pente est souvent nulle, un autre exemple est une marche sur une pente inclinée de pente constante). Grâce à cette méthode, dans la trajectoire estimée, le pied est au sol à chaque instant d'immobilité Par exemple, lesdits moyens d'estimation sont adaptés pour estimer ledit angle à partir d'une décomposition en valeurs singulières de la matrice à trois colonnes fournie par lesdits cinquièmes moyens de calcul à partir de la matrice formée par les vecteurs délivrés par l'accéléromètre triaxial auxdits instants successifs. En ce cas, aucune connaissance de la pente n'est nécessaire, on se base uniquement sur la statistique du signal. Par exemple, lesdits moyens d'estimation sont adaptés pour estimer ledit angle comme l'angle, audit premier ou deuxième instant 15 d'immobilité, entre le vecteur délivré par lesdits premiers moyens de calcul à partir du vecteur fourni par le capteur additionnel, et le vecteur délivré par lesdits premiers moyens de calcul à partir du vecteur fourni par l'accéléromètre triaxial. Dans ce cas, on n'a pas besoin de cartographie, ni d'hypothèse 20 sur la trajectoire de type trajectoire plus longue que haute. Dans un mode de réalisation, lesdits moyens de commande comprennent des troisièmes moyens d'élaboration de la distance parcourue entre lesdits premier et deuxième instants d'immobilité, et/ou lors d'un ensemble d'intervalles de temps successifs séparant deux 25 instants successifs d'immobilité de l'accéléromètre triaxial, à partir des données fournies par lesdits quatrièmes moyens de calcul. On calcule facilement la distance parcourue entre les deux instants successifs d'immobilité de l'accéléromètre triaxial, et donc, dans le cas du pied, la distance parcourue par le pied entre deux appuis 30 successifs. Selon un mode de réalisation, lesdits moyens de commande comprennent des quatrièmes moyens d'élaboration d'une trajectoire durant un ensemble d'intervalles de temps successifs séparant deux instants successifs d'immobilité de l'accéléromètre triaxial, à partir des 35 données fournies par lesdits quatrièmes moyens de calcul. L'invention permet d'élaborer, avec un nombre de calculs réduits, la trajectoire d'un utilisateur durant un déplacement. Dans un mode de réalisation, ledit capteur additionnel est un magnétomètre triaxial de mesure du champ magnétique terrestre, 5 sensiblement constant entre lesdits premier et deuxième instants d'immobilité. Ainsi le coût est réduit, par exemple en comparaison de l'utilisation de gyromètres. Selon un mode de réalisation, ledit accéléromètre triaxial est 10 disposé de manière à se positionner sensiblement au niveau de la surface de contact de l'élément mobile, à un instant d'immobilité dudit accéléromètre triaxial. En ce cas, l'accélération est bien nulle lors d'un appui du pied, ce qui améliore la précision du dispositif. 15 Dans un mode de réalisation, ledit capteur additionnel est disposé de manière à se positionner à distance de la surface de contact dudit élément mobile à un instant d'immobilité dudit accéléromètre triaxial. On évite les perturbations de mesure du champ magnétique terrestre au niveau du sol, ce qui améliore la précision du dispositif. 20 Selon un mode de réalisation, lesdits moyens de commande sont liés à l'élément mobile ou disposés à distance de l'élément mobile, et adaptés pour fonctionner en temps réel ou en différé. On peut réaliser un dispositif dans lequel lesdits moyens de commande se déplacent avec l'élément mobile ou non, et qui traite les 25 informations en temps réel ou non. Selon un autre aspect de l'invention, il est proposé un système de jeu vidéo comprenant un indicateur visuel ou auditif, tel un écran ou un haut-parleur, pour indiquer au moins un pas à effectuer par le joueur, comprenant au moins un dispositif tel que décrit précédemment. 30 Il est ainsi possible de contrôler la bonne exécution de la consigne et, par exemple, établir un score. Selon un autre aspect de l'invention, il est proposé un système d'analyse de saut, comprenant au moins un dispositif tel que décrit précédemment. On peut ainsi estimer la hauteur, la longueur, la durée, la qualité d'un saut. Selon un autre aspect de l'invention, il est proposé un système de géolocalisation en extérieur ou en intérieur, comprenant au moins un 5 dispositif tel que décrit précédemment. Ce système est un système bas coût qui peut être complété par une initialisation manuelle (l'utilisateur indique où il se trouve sur une carte) ou à l'aide d'un autre système de géolocalisation (GPS, ULB ...). Dans ce dernier cas, le système proposé viens améliorer les 10 performances de l'autre système ou se substituer à lui lorsqu'il n'est plus disponible (exemple le GPS dans les bâtiments). Selon un autre aspect de l'invention, il est proposé un système de calcul de dépense énergétique d'un utilisateur, comprenant au moins un dispositif tel que décrit précédemment. 15 Ce dispositif apporte une information importante sur le déplacement du pied. On peut ainsi intégrer ces informations sur la durée du pas afin d'obtenir une évaluation fine de la dépense énergétique qui tient compte de la hauteur du pied à chaque instant. L'invention sera mieux comprise à l'étude de quelques modes 20 de réalisation décrits à titre d'exemples nullement limitatifs et illustrés par les dessins annexés sur lesquels : - la figure 1 illustre schématiquement un mode de réalisation d'un dispositif, selon un aspect de l'invention ; et - la figure 2 illustre schématiquement un mode de réalisation 25 des moyens de commande d'un dispositif de la figure 1, selon un aspect de l'invention. Dans l'ensemble de figures, les éléments ayants les mêmes références sont similaires. Tel qu'illustré sur la figure 1, un module de commande CMD 30 reçoit les mesures transmises par un accéléromètre triaxial 3A et un capteur additionnel de mesure d'un vecteur de champ vectoriel sensiblement constant localement, en l'occurrence un magnétomètre triaxial 3M de mesure du champ magnétique terrestre. L'accéléromètre triaxial 3A est lié à un élément mobile EM, et le magnétomètre 3M est 35 également lié au même élément mobile EM, directement ou indirectement, et ses axes de mesures sont sensiblement parallèles aux axes de mesures de l'accéléromètre. L'élément mobile EM peut être, par exemple, de manière non limitative, un pied d'un corps humain. Un premier module de détermination MD1 permettent, à pa des mesures transmises par l'accéléromètre triaxial 3A et le magnétomètre triaxial 3M, de déterminer un premier instant d'immobilité to et un deuxième instant d'immobilité tä de l'accéléromètre triaxial, consécutif audit premier instant d'immobilité to. Lorsque l'élément mobile EM est un pied, un instant d'immobilité correspond à un appui du pied sur la surface de contact. La technique utilisée par le premier module de détermination MD1 peut, par exemple être celle décrite dans le document "Robust step detection method for pedestrian navigation systems" (D.H. Hwang, J.H. Jang, et J.W. Kim) publié dans Electronics Letters, 43(14), en juillet 2007, dans lequel les auteurs cherchent à détecter une forme de base, correspondant à un signal typique pour un pas, dans le signal de la norme des accéléromètres. Lorsque cette forme est détectée, le début correspond à l'instant to et la fin, à l'instant tn. Les instants successifs to, t1, ...tn sont situés dans le temps, de manière chronologique entre le premier instant d'immobilité to et le deuxième instant d'immobilité tn. Les différents instants successifs to, ti, ...tn correspondent à une numérisation et à un échantillonnage des signaux issus des capteurs suivie d'une calibration pour remonter aux grandeurs physiques mesurées. Un deuxième module de détermination MD2 permet de 25 déterminer le plan de ladite trajectoire de l'accéléromètre triaxial 3A, à partir des vecteurs délivrés par l'accéléromètre triaxial 3A, aux différents instants successifs to, t1, ...,tn, dans un repère mobile LF lié à l'accéléromètre triaxial 3A et au magnétomètre triaxial 3M. En variante, comme illustré sur la figure 2, le deuxième module 3o de détermination MD2 peut déterminer le plan P1,1-i2] de ladite trajectoire à partir des vecteurs délivrés par magnétomètre triaxial 3M. La méthode utilisée par cette variante est analogue, en remplaçant les mesures de l'accéléromètre triaxial 3A par les mesures du magnétomètre triaxial 3M. Le plan [û,,i,] correspond au plan contenant au maximum les 35 mesures, considérées comme des points ou vecteurs dans un espace à trois dimensions en considérant les mesure transmis pour chaque axe de l'accéléromètre triaxial 3A, ou respectivement du magnétomètre triaxial, et l'axe porté par le vecteur est l'axe sur lequel la projection des mesures varie le moins. Le vecteur est le vecteur normal au plan de la trajectoire mais aussi l'axe de rotation de l'accéléromètre triaxial 3A. Pour déterminer ces trois vecteurs iil,Ü2,û3, les deuxièmes moyens de détermination forment une matrice ALE à n+1 lignes, n+1 étant le nombre d'échantillons pris entre to et compris, et à 3 colonnes correspondant aux trois mesures des trois axes x, y, z de l'accéléromètre LF (t0) A LE (t0) Azu�to)- triaxial 3A. On a : A LF(t,, A Li'(t ALE t, kt,, En, variante, on peut ne prendre qu'un sous-ensemble des échantillons aux instants successifs to, Plusieurs modes de réalisation sont possibles pour déterminer 15 les trois vecteurs ûä ûä , , à partir de la matrice A Lf�= ALI�to �. Un mode de réalisation consiste à effectuer la décomposition en valeurs singulières SVD de la matrice ALE: A"=USV. V est une matrice à trois lignes et trois colonnes dont chaque colonne correspond à un des trois vecteurs La matrice S contient les trois valeurs 20 singulières si, s2, s3. Ces trois valeurs singulières si, s2, s3, doivent êtres classées par ordre décroissant. Le vecteur Ü3 est la colonne de la matrice V correspondant à la plus petite valeur singulière. Les deux autres colonnes donnent les coordonnées des vecteurs formant le plan [à,,-12]. Un autre mode de réalisation consiste à décomposer la matrice 25 de corrélation R=ALLT A"' en trois vecteurs propres et valeurs propres (R=uTSu) et à les classer par valeurs propres décroissantes s1,s2,s3 pour obtenir les trois vecteurs T correspond à la matrice transposée de A'' . Les deux vecteurs propres correspondant aux deux plus grandes valeurs propres correspondent au plan le vecteur 30 propre ïii est celui correspondant à la plus petite valeur propre. Le module de commande CMD comprend également un premier module de calcul MC1 pour calculer, aux instants successifs to, tn, des projections orthogonales du vecteur délivré par l'accéléromètre triaxial 3A et du vecteur délivré par le magnétomètre triaxial 3M dans ledit repère mobile LF, sur le plan de la trajectoire [û,,û2] de l'accéléromètre triaxial 3A, on (- obtient les matrices des projections: AL-F (tl~ ...tn) = [A(t0...tr,) ul , A(tll...tn) u2 ] B ' (t0...t,,) = [B(tl B(t( ..tr,) u2 Le module de commande CMD comprend, en outre, un troisième module de détermination MD3, aux instants successifs to, de la rotation de passage du repère mobile LF à un repère global fixe GF, lié au repère terrestre, à partir des projections orthogonales des ~o vecteurs délivrés par le magnétomètre triaxial 3M. Le troisième module de détermination MD3 détermine avantageusement un plan intermédiaire IF formé par les vecteurs bp (t), et bpi (t) , bp (t) étant le vecteur unitaire du champ magnétique terrestre BLF(t) et bpi(t) étant obtenu par la rotation de 90' du vecteur hp(t) dans le plan de la trajectoire ou de la marche [ûäâ72]. 15 Ainsi, le troisième module de détermination élabore les vecteurs du repère intermédiaire IF: LF hp(t) = Bp (t) Un deuxième module de calcul MC2, du module de commande CMD, effectue le calcul, aux instants successifs to, t1, ..., ta,, de la 20 projection orthogonale, dans le repère global fixe GF, du vecteur délivré par l'accéléromètre triaxial 3A fournie par le premier module de calcul MC1. Le deuxième module de calcul MC2 comprend un cinquième module de calcul MC5 de calcul de passage dans le repère intermédiaire 25 IF, lié au repère terrestre : En d'autres termes, on projette orthogonalement A`� (t0...tä) dans le plan de la trajectoire [ûl,û2] sur les vecteurs hp(t), et h, (t) du repère intermédiaire IF pour obtenir les mesures de l'accéléromètre projetées sur le plan de la trajectoire [ûl,û2 ] et exprimées dans le repère 30 intermédiaire IF : A" (t) = [A i" (t) h(t) (t) . h (t)~ = [A A Jt. J Le troisième module de détermination MD3 de la rotation de passage du repère mobile LF au repère global fixe GF comprend un et bp1(t) = ~LI90bp(t) _ .bp(t) 1 -1 0 30 module d'estimation MEST d'un angle a, ladite rotation de passage étant définie par une projection orthogonale sur le champ magnétique terrestre sensiblement constant entre les deux instants d'immobilité to et t,,, suivie d'une rotation de l'angle estimé a dans le plan de trajectoire de l'accéléromètre triaxial 3A. En d'autres termes, le module d'estimation MEST permet de trouver la rotation a (caractérisé par une matrice de rotation Mu) à effectuer pour passer du repère intermédiaire IF {b,,(t),hr (t)} au repère global GF dans le plan de la marche porté par les vecteurs Vh,i , le troisième vecteur du repère global GF étant i7 Dans le cas d'une navigation d'un piéton, est l'horizontale dans le plan de la marche et v la verticale. Pour déterminer a, plusieurs modes de réalisation du module d'estimation MEST sont possibles. Une premier mode de réalisation peut consister, connaissant la pente 8 entre les deux points d'appui du pied aux instants to et tä par rapport à l'horizontale, à estimer a par l'équation suivante : (r)drdt -tan 1" AbF (z)(1 zdt + tan 8 f_ Abi (r)d rdt Cette pente 8 peut être estimée grâce à l'inclinaison des capteurs sur le pied lorsqu'il est en appui, à to ou t,,, en considérant que , la pente est constante pendant la durée du pas, ou à l'aide d'une cartographie mémorisée délivrant en sortie la pente de la surface de contact de l'endroit où la personne (l'élément mobile EM) est située, Dans un autre mode réalisation, en admettant que la plus grande variance de l'accélération est appliquée suivant le déplacement du pied, la décomposition en valeurs singulières ou SVD de la matrice d'accélération AIF détermine la matrice Ma , le vecteur singulier correspondant à la plus grande valeur singulière correspond au vecteur h . La décomposition en valeurs singulières s'écrit : A =USV où V est la matrice Ma, Cette hypothèse n'est valable que dans le cas de grands pas pour lequel le déplacement horizontal est plus grand que le déplacement vertical sur un sol plat, de pente nulle. Si la pente n'est pas nulle, on peut utiliser cette méthode mais ïl ne sera pas l'horizontale mais le vecteur reliant la position du capteur au premier instant d'immobilité to à sa position au deuxième instant d'immobilité tn. Un autre mode de réalisation peut également être envisagé, en 10 partant du fait qu'aux instants ta et tn le pied est posé, la verticale est suivant l'accélération mesurée à cet instant par l'accéléromètre triaxial 3A et a est l'angle entre ;Ip''(to) et hp'(to), ou entre Âpi'(tn) et fipi'(tn) Le deuxième module de calcul MC2 comprend également un multiplicateur pour projeter les mesures Ap dans le repère global : -À"F (t) A "''�t) A"f�t) = [Aa'F (t) Aa",�t) Mu(,) = An t) .1)(0 Au(t) . b (t 15 cos(a) sin(a) et Mu = û sin(a) cos(a) Le module de commande CMD comprend, en outre, un troisième module de calcul MC3 pour soustraire au vecteur délivré par 20 lesdits deuxièmes moyens de calcul MC2 auxdits instants successifs to, tn, le vecteur moyen sur lesdits instants successifs, de manière à obtenir les accélérations centrées dans le plan [iil,1 21 de la trajectoire de l'accéléromètre triaxial 3A, dépourvues de l'influence de la gravité terrestre et de dérives du dispositif, dans le repère global fixe GF. Le troisième 25 module de calcul réalise les calculs suivants de retraits des moyennes horizontales et verticales : A"' (t) - JAfcr' (r)dr A F (t) - h tn ù tO r=to 1,, JA1'-- (r )dr ùfo A"1-(t) n Le module de commande CMD comprend également un quatrième module de calcul MC4 de la double intégration temporelle des 30 accélérations centrées délivrées par le troisième module de calcul MC3, de manière à obtenir la trajectoire de l'accéléromètre triaxial 3A c'est à dire du pied de l'utilisateur, dans le plan de la trajectoire du pied durant un pas, exprimée dans le repère global fixe GF, entre le premier appui et le deuxième appui, ou instants d'immobilité, to et tn. Le quatrième module de calcul MC4 permet d'intégrer deux fois, par rapport au temps, l'accélération pour calculer la trajectoire du pied, ou plus précisément de l'accéléromètre lié au pied, dans le plan de la trajectoire ou de la marche. L'intégration numérique peut être réalisée par mise en oeuvre des relations suivantes, utilisant l'intégration par la méthode des rectangle 17(t)=v(tù1)+AG� /Fe (t)=[x,, x,]=z(t-1)+ii(t)/Fe+AGNI dans laquelle ; (to) = [0 0] et .i(to) = [0 0 Dans lesquelles Fe est la fréquence d'échantillonnage exprimée en Hz, et les composantes de la matrice ÀCG sont exprimées en m/s2. v et x sont respectivement les vecteurs vitesse et position selon les deux axes du repère global {h,v}. 15 En variante, le quatrième module de calcul MC4 peut mettre en oeuvre la double intégration par des modes de calculs numériques plus complexes, par exemples par utilisation de Splines. Le module de commande CMD peut comprendre un premier module d'élaboration MEL1 d'un vecteur de direction, adapté pour relier 20 les positions respectives de l'accéléromètre triaxial 3A au premier instant d'immobilité to et au deuxième instant d'immobilité t,,, exprimées dans ledit repère global fixe GF. Le module de commande CMD peut également comprendre un deuxième module d'élaboration MEL2 d'un cap de déplacement par 25 rapport audit champ vectoriel mesuré par le capteur additionnel, en l'occurrence par rapport au champ magnétique terrestre mesuré par le magnétomètre triaxial 3M. Le deuxième module d'élaboration MEL2 peut élaborer le cap de la manière qui suit. Le cap du déplacement est défini comme l'angle entre le pâle 30 nord magnétique terrestre et la direction du déplacement. Il ne faut pas le confondre avec le cap de la personne, qui est son orientation par rapport au Nord. C'est souvent ce dernier qui est estimé dans la littérature. Le cap de la personne est souvent supposé être le même que celui du déplacement. Or, cela est le cas uniquement dans le cas d'une marche vers l'avant mais pas lorsque l'on marche en arrière, sur le côté, ou un peu en biais. Aussi, il est proposé pour estimer le cap du déplacement et plus généralement pour estimer la trajectoire dans le repère NED (North,East,Down) un deuxième module d'estimation MES2 mettant en oeuvre ce qui suit. Dans ce qui précède, la trajectoire est estimée dans un plan correspondant au plan de la trajectoire du pied. Le troisième vecteur ii, qui permet de former un repère à trois dimensions {Ï1,17,ii, correspondant respectivement à l'horizontale dans le plan de la marche, à la verticale dans le plan de la marche et à un vecteur normal au plan de la marche, qui est également supposé horizontal. Pour calculer le cap, on effectue le calcul : = Bi, + B,, + Bä3 dans lequel : [cosa 0 0] est la projection des mesures du magnétomètre triaxial 3M BGF sur h={1,0,0} [0 ûsin a 0] est la projection des mesures du magnétomètre triaxial 3M BGF sur v={0,1,0} ; et GF Gr = Bu- (t).ui" est la projection des mesures du magnétomètre triaxial 3M BGF sur LI,' ={0,0,1} et 20 113 est le vecteur orthogonal au plan de la marche, exprimé dans le repère des capteurs et calculé précédemment. La composante B,,3 doit être sensiblement invariante en fonction du temps, ce qui revient à faire l'hypothèse d'une marche dans un 25 plan et plus généralement d'une rotation autour d'un axe sensiblement invariant et que le champ magnétique terrestre B n'est pas perturbé. Comme ce n'est généralement pas le cas en pratique, le deuxième module d'élaboration, peut prendre la moyenne temporelle sur le pas (entre les instants d'immobilité tO et tn) comme estimation de cette 30 composante. Pour déterminer le cap seuls BhetB,,3 sont utiles. Il reste des indéterminations sur le signe des vecteurs 'i et ïï. En effet, tels qu'ils ont été définis ci-avant, ïl peut être dans le sens inverse du sens du déplacement et peut être dirigé vers le haut ou vers le bas. Ba- Bh = Bu' = Le deuxième module d'estimation MES2 construit un repère {P1,P2,P3} tel que : : est un vecteur horizontal dans le plan de la marche (parallèle à ïl ) et dans le sens du déplacement des capteurs ; p2 est un vecteur vertical (parallèle à ï;.) dirigé vers le bas 7), est un vecteur perpendiculaire au plan de la marche (parallèle à formant un repère direct }p,, p,} Les vecteurs 1-51,152,153 sont déterminés comme suit. est un vecteur horizontal qui correspond au sens du déplacement de l'utilisateur du dispositif. Il est obtenu grâce à au déplacement horizontal xh(t) qui correspond à la première composante du vecteur de déplacement -(1) donné par le quatrième module de calcul MC4. Si xh(t) est positif, alors = h ù' Il 0 0} ,sinon 15, = ùh = 1ù1 0 01, le déplacement devant être positif dans le sens du déplacement. est vertical et orienté vers le bas. Si la projection de B(t) sur v" est positive alors 1-52 -1)ù 10 1 01, sinon pz = -v = {o -1 o}, le champ magnétique terrestre B mesuré devant toujours être orienté vers le bas. est le troisième vecteur formant un repère direct {A,Ï52,15,} avec les deux autres. L'angle q) entre la direction du déplacement pl et le nord magnétique dans le repère {Nord, Est} peut être obtenu en inversant les formules suivantes : BGF• 3I cosO BGI 2 BGl .P j BG.� (BGF + BGF .Î)3 )2 Le module de commande CMD peut aussi comprendre un troisième module d'élaboration MEL3 de la distance parcourue entre les premier et deuxième instants d'immobilité to et tn, i.e. un pas, et/ou lors d'un ensemble d'intervalles de temps successifs séparant deux instants successifs d'immobilité to et tn de l'accéléromètre triaxial, i.e. une succession de pas, à partir des données fournies par le quatrième module de calcul. Par exemple la distance horizontale parcourue durant un pas est donnée par : sin q, = d1, (to , t,,) = x1, (ti, ) Pour une série de pas, la distance parcourue est donnée par la somme des distances calculées pour chaque pas. Le module de commande CMD peut également comprendre un quatrième module d'élaboration MEL4 d'une trajectoire durant un ensemble d'intervalles de temps successifs séparant deux instants successifs d'immobilité tO et tn de l'accéléromètre triaxial 3A, en l'espèce durant un ensemble de pas successifs, à partir des données fournies par le quatrième module de calcul. Par exemple, au moyen des mêmes formules que celles utilisées par le deuxième module d'élaboration MEL2, ou directement à partir des résultats du deuxième module d'élaboration MEL2 si le dispositif les comprend. Le déplacement dans le plan horizontal correspondant est donné par xh(t). Le déplacement selon l'axe dirigé vers le Nord est donné par la relation suivante : (BGF ..Ï-51 .Î)3 Le déplacement selon le bas est donné par la relation suivante : (t) �x,, (t) , le signe étant donné le signe de 32 qui est orienté vers le bas. On obtient ainsi la trajectoire estimée dans le plan horizontal (xN(t), xE(t)) pour une série de pas. On peut également estimer le déplacement dans tout repère global GF pour lequel on connait la matrice de rotation permettant de passer du repère NED au repère voulu. Par exemple, lorsqu'une personne est devant sa console de jeu. On peut faire une acquisition des mesures lorsque la personne est immobile et face à son écran. Le cap de la personne qui est calculé sert de référence pour l'orientation de son personnage virtuel. En variante, la prise en compte des deux instants d'immobilité successifs, peuvent être remplacés par la prise en compte de deux 19 (B"F.p BGF .p-3 Le déplacement selon l'axe dirigé vers l'Est est donné par la relation suivante : ù x,,(t el" x N (t) = x,, (t) BGF .P3 instants successifs pour lesquels les vecteurs vitesses dans le plan de la trajectoire sont connus, ainsi que de l'angle ci à ces deux instants, ces données pouvant être mesurées par des capteurs additionnels ou estimées à partir de mesures fournies par des capteurs additionnels. Dans un tel cas, le quatrième module de calcul doit initialiser la vitesse avec la vitesse du premier instant et ajouter une constante dépendant des vitesses des deux instants à l'accélération centrée. La présente invention permet, à coût réduit de déterminer une 10 trajectoire planaire d'un accéléromètre triaxial lié à un élément mobile, notamment un pied, de manière précise et à coût réduit. L'élément mobile peut être une autre partie du corps humain, comme une main, une partie d'un corps animal, ou une partie d'un corps artificiel tel un robot ou une souris informatique. 15 REVENDICATIONS1. Dispositif de détermination d'une trajectoire formée de positions successives sensiblement coplanaires d'un accéléromètre triaxial (3A) lié de manière solidaire à un élément mobile (EM) d'un corps humain, animal ou artificiel, entre un premier instant d'immobilité (t0) et un deuxième instant d'immobilité (tn) de l'accéléromètre triaxial (3A), consécutif audit premier instant d'immobilité (to), l'accéléromètre triaxial (3A) étant en rotation par rapport à une direction sensiblement constante entre lesdits premier et deuxième instants d'immobilité (ta,tn), et ledit dispositif comprenant, en outre, un capteur additionnel de mesure d'un vecteur de champ vectoriel sensiblement constant entre lesdits premier et deuxième instants d'immobilité Mn), dans un repère global fixe (GF) lié au repère terrestre, ledit capteur additionnel étant lié de manière solidaire audit élément mobile (EM) et d'axes de mesures sensiblement parallèles aux axes de mesures de l'accéléromètre (3A), et des moyens de commande (CMD), caractérisé en ce que lesdits moyens de commande (CMD) comprennent : des premiers moyens de détermination (MD1) desdits premier et deuxième instants d'immobilité (to,tn) de l'accéléromètre triaxial (3A) ; des deuxièmes moyens de détermination (MD2) du plan de ladite trajectoire de l'accéléromètre triaxial (3A), à partir des vecteurs délivrés par l'accéléromètre triaxial (3A), ou à partir des vecteurs délivrés par le capteur additionnel, à différents instants successifs (ta,t1,...,tn) entre le premier et le deuxième instants d'immobilité, dans un repère mobile (LF) lié à l'accéléromètre triaxial (3A) et au capteur additionnel ; des premiers moyens de calcul (MC1), auxdits instants successifs (to,ti,...,tn), des projections orthogonales du vecteur délivré par l'accéléromètre triaxial (3A) et du vecteur délivré par le capteur additionnel dans ledit repère mobile (LF), sur ledit plan de ladite trajectoire de l'accéléromètre triaxial (3A) ; des troisièmes moyens de détermination (MD3), auxdits instants successifs (to,t1,...,tn), de la rotation de passage dudit repère mobile(LF) audit repère global fixe (GF), à partir desdites projections orthogonales des vecteurs délivrés par le capteur additionnel ; des deuxièmes moyens de calcul (MC2), auxdits instants successifs (tO,tl,...,tn), de la projection orthogonale dans ledit repère global fixe (GF), du vecteur délivré par l'accéléromètre triaxial (3A), fournie par lesdits premiers moyens de calcul (MC1) ; des troisièmes moyens de calcul (MC3) pour soustraire au vecteur délivré par lesdits deuxièmes moyens de calcul (MC2) auxdits instants successifs (ta,ti,...,tn), le vecteur moyen sur lesdits instants successifs (to,ti,...,tn), de manière à obtenir les accélérations centrées dans ledit plan de ladite trajectoire de l'accéléromètre triaxial (3A), dépourvues de l'influence de la gravité terrestre et de dérives dudit dispositif, dans ledit repère global fixe (GF) ; et des quatrièmes moyens de calcul (MC4) de la double intégration temporelle desdites accélérations centrées délivrées par lesdits troisièmes moyens de calcul (MC3), de manière à obtenir la trajectoire de l'accéléromètre triaxial (3A) dans ledit plan, exprimée dans ledit repère global fixe (GF), entre lesdits premier et deuxième instants d'immobilité (to,tn). 2. Dispositif selon la revendication 1, dans lequel ledit élément mobile (EM) est un pied, l'intervalle de temps séparant lesdits premier et deuxième instants d'immobilité (tO,tn) de l'accéléromètre triaxial (3A) correspondant à un pas dudit pied. 3. Dispositif selon la revendication 1 ou 2, dans lequel lesdits moyens de commande (CMD) comprennent des premiers moyens d'élaboration (MEL1) d'un vecteur de direction, adapté pour relier les positions respectives de l'accéléromètre triaxial (3A) audit premier instant d'immobilité (4)) et audit deuxième instant d'immobilité exprimées dans ledit repère global fixe (GF). 4. Dispositif selon la revendication 3, dans lequel ledit plan de la trajectoire étant vertical, lesdits moyens de commande (CMD) comprennent, en outre, des quatrièmes moyens de calcul (MC4) d'uneprojection orthogonale, dudit vecteur de champ vectoriel mesuré par le capteur additionnel, sur un plan horizontal orthogonal audit plan de la trajectoire, et des deuxièmes moyens d'élaboration (MEL2) d'un cap de déplacement par rapport audit champ vectoriel mesuré par le capteur additionnel. 5. Dispositif selon l'une des revendications 1 à 4, dans lequel ledit plan de la trajectoire étant vertical, lesdits deuxièmes moyens de détermination (MD2) du plan de ladite trajectoire de l'accéléromètre triaxial (3A) sont adaptés pour effectuer une décomposition en valeurs singulières de la matrice à trois colonnes (A`' (t0...t,,)) formée par les vecteurs délivrés par l'accéléromètre triaxial (3A) auxdits instants successif (to,t1,...,tn), dans ledit repère mobile, et extraire les deux vecteurs correspondant aux deux plus grandes valeurs singulières. 6. Dispositif selon l'une des revendications 1 à 4, dans lequel ledit plan de la trajectoire étant vertical, lesdits deuxièmes moyens de détermination (MD2) du plan de ladite trajectoire de l'accéléromètre triaxial (3A) sont adaptés pour effectuer une décomposition en valeurs propres et vecteurs propres, dans ledit repère mobile (LF), la matrice d'autocorrélation de la matrice à trois colonnes formée par les vecteurs délivrés par l'accéléromètre triaxial (3A) auxdits instants successifs (to,t1,...,tn), et extraire les deux vecteurs propres correspondant aux deux plus grandes valeurs propres. 7. Dispositif selon l'une des revendications 1 à 6, dans lequel lesdits deuxièmes moyens de calcul (MC2) comprennent des cinquièmes moyens de calcul (MC5) de passage dans un repère intermédiaire (IF), lié au repère terrestre, défini par un premier vecteur délivré par lesdits premiers moyens de calcul (MC1) à partir d'un vecteur fourni par le capteur additionnel, et par un deuxième vecteur, orthogonal audit premier vecteur et appartenant audit plan de la trajectoire. 8. Dispositif selon l'une des revendications 1 à 7, dans lequel, lesdits troisièmes moyens de détermination (MD3) de la rotation de passage dudit repère mobile (LF) audit repère global fixe (GF) comprennent des moyens d'estimation (MEST) d'un angle (a), ladite rotation de passage étant définie par une projection orthogonale sur ledit champ vectoriel sensiblement constant, suivie d'une rotation dudit angle estimé (a) dans ledit plan de trajectoire de l'accéléromètre triaxial (3A). 9. Dispositif selon la revendication 8, dans lequel lesdits moyens d'estimation (MEST) sont adaptés pour estimer ledit angle (a) à partir de la pente (0) de la surface de contact de l'élément mobile (EM), et du vecteur transmis par lesdits cinquièmes moyens de calcul (MC5) à partir du vecteur transmis par l'accéléromètre triaxial (3A). 10. Dispositif selon la revendication 8, dans lequel lesdits moyens d'estimation (MEST) sont adaptés pour estimer ledit angle (c) à partir d'une décomposition en valeurs singulières de la matrice à trois colonnes fournie par lesdits cinquièmes moyens de calcul (MC5) à partir de la matrice (AiF) formée par les vecteurs délivrés par l'accéléromètre triaxial (3A) auxdits instants successifs (tO,tl,...,tn). 11. Dispositif selon la revendication 8, dans lequel lesdits moyens d'estimation (MEST) sont adaptés pour estimer ledit angle (a) comme l'angle, audit premier ou deuxième instant d'immobilité (tO,tn), entre le vecteur délivré par lesdits premiers moyens de calcul (MC1) à partir du vecteur fourni par le capteur additionnel, et le vecteur délivré par lesdits premiers moyens de calcul (MC1) à partir du vecteur fourni par l'accéléromètre triaxial (3A). 12. Dispositif selon l'une des revendications précédentes, dans lequel lesdits moyens de commande (CMD) comprennent des troisièmes moyens d'élaboration (MEL3) de la distance parcourue entre lesdits premier et deuxième instants d'immobilité (to,tn), et/ou lors d'un ensemble d'intervalles de temps successifs séparant deux instantssuccessifs d'immobilité (to,tn) de l'accéléromètre triaxial (3A), à partir des données fournies par lesdits quatrièmes moyens de calcul (MC4). 13. Dispositif selon l'une des revendications précédentes, dans lequel lesdits moyens de commande (CMD) comprennent des quatrièmes moyens d'élaboration (MEL4) d'une trajectoire durant un ensemble d'intervalles de temps successifs séparant deux instants successifs d'immobilité (to,tn) de l'accéléromètre triaxial (3A), à partir des données fournies par lesdits quatrièmes moyens de calcul (MC4). 14. Dispositif selon l'une des revendications précédentes, dans lequel ledit capteur additionnel est un magnétomètre triaxial (3M) de mesure du champ magnétique terrestre, sensiblement constant entre lesdits premier et deuxième instants d'immobilité (to,tn). 15. Dispositif selon l'une des revendications précédentes, dans lequel ledit accéléromètre triaxial (3A), est disposé de manière à se positionner sensiblement au niveau de la surface de contact de l'élément mobile (EM) à un instant d'immobilité dudit accéléromètre triaxial (3A). 16. Dispositif selon l'une des revendications précédentes, dans lequel ledit capteur additionnel est disposé de manière à se positionner à distance de la surface de contact dudit élément mobile (EM), à un instant d'immobilité dudit accéléromètre triaxial (3A). 1 7. Dispositif selon l'une des revendications précédentes, dans lequel lesdits moyens de commande (CMD) sont liés à l'élément mobile (EM) ou disposés à distance de l'élément mobile (EM), et adaptés pour fonctionner en temps réel ou en différé. 18. Système de jeu vidéo comprenant un indicateur visuel ou auditif pour indiquer au moins un pas à effectuer par le joueur, caractérisé en ce qu'il comprend au moins un dispositif selon l'une des revendications 2 à 16.19. Système d'analyse de saut caractérisé en ce qu'il comprend au moins un dispositif selon l'une des revendications 2 à 16. 20. Système de géolocalisation en extérieur ou en intérieur, caractérisé en ce qu'il comprend au moins un dispositif selon l'une des revendications 2 à 16. 21. Système de calcul de dépense énergétique d'un utilisateur, caractérisé en ce qu'il comprend au moins un dispositif selon l'une des revendications 2 à 16.